The high strain mechanical response of the wet Kelvin model for open-cell foams

Surface Evolver software was used to create the three-dimensional geometry of a Kelvin open-cell foam, to simulate that of polyurethane flexible foams. Finite Element Analysis (FEA) with 3D elements was used to model large compressive deformation in the [0 0 1] and [1 1 1] directions, using cyclic boundary conditions when necessary, treating the polyurethane as an elastic or elastic–plastic material. The predicted foam Young’s moduli in the [0 0 1] direction are double those of foams with uniform Plateau border cross-section edges, for the same foam density and material properties. For compression in the [1 1 1] direction, the normalized Young’s modulus increases from 0.9 to 1.1 with foam relative density, and the predicted stress–strain relationship can have a plateau, even for a linearly-elastic polymer. As the foam density increases, the predicted effects of material plasticity become larger. For foam of relative density 0.028, edge-to-edge contact is predicted to occur at a 66% strain for [1 1 1] direction compression. The foam is predicted to contract laterally when the [1 1 1] direction compressive strain exceeds 25%.