Note: This discussion is about an older version of the COMSOL Multiphysics® software. The information provided may be out of date.

Discussion Closed This discussion was created more than 6 months ago and has been closed. To start a new discussion with a link back to this one, click here.

Uniaxial compression test on representative foam: symmetries!

Please login with a confirmed email address before reporting spam

Hello everyone,
I came trough this problem and since I'm quite a new COMSOL user perhaps someone can help me.

I’m trying to reproduce an uniaxial compression test (structural mechanics) on a simplified cellular material (cell is represented by tetrakaidecahedron). Cells edges are modeled as beams.

In order to reproduce the periodic conditions on the unloaded faces I have to apply symmetries boundary conditions. At the first three orthogonal faces generic symmetries were applied (see attached file) . At the upper face a displacement is imposed (load condition).

At the other two (lateral) faces I would like to apply a symmetry-like condition (in plane translations and orthogonal rotation allowed) but with a moving reference plane (translated in the perpendicular to the sample direction). Practically I would reproduce the warping effect keeping the lateral faces planar. Simple symmetries may not reproduce this condition since lateral displacement would be locked.

Could anyone please help me with this problem?
Thank you!

Claudio D'Angelo


2 Replies Last Post 26 juil. 2012, 01:52 UTC−4
Nagi Elabbasi Facebook Reality Labs

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 25 juil. 2012, 09:55 UTC−4
Hi Claudio,

Instead of symmetry condition where u_Normal = 0, set u_Normal everywhere on the face equal to say u_1 where u_1 is the displacement of one node. You should use a Model Coupling operator for that. That way the whole face will move together in the normal direction. You can add rotational terms as well if you want that plane to also rotate.

Note that a general periodic constraint is more complex because the nodes on one face are not totally free to move within the plane. They are also constrained by the motion of the nodes on the opposing face in order not to violate periodic continuity.

Nagi Elabbasi
Veryst Engineering
Hi Claudio, Instead of symmetry condition where u_Normal = 0, set u_Normal everywhere on the face equal to say u_1 where u_1 is the displacement of one node. You should use a Model Coupling operator for that. That way the whole face will move together in the normal direction. You can add rotational terms as well if you want that plane to also rotate. Note that a general periodic constraint is more complex because the nodes on one face are not totally free to move within the plane. They are also constrained by the motion of the nodes on the opposing face in order not to violate periodic continuity. Nagi Elabbasi Veryst Engineering

Please login with a confirmed email address before reporting spam

Posted: 1 decade ago 26 juil. 2012, 01:52 UTC−4
Dear Nagi,
thank you for your kind reply.
By defining a coupling factor "average" on the boundary nodes and applying that value to the prescribed displacement it seems to work. I will also apply rotational constrain to that faces in order to reproduce continuity.
Thank you again.

Claudio D'Angelo
Dear Nagi, thank you for your kind reply. By defining a coupling factor "average" on the boundary nodes and applying that value to the prescribed displacement it seems to work. I will also apply rotational constrain to that faces in order to reproduce continuity. Thank you again. Claudio D'Angelo

Note that while COMSOL employees may participate in the discussion forum, COMSOL® software users who are on-subscription should submit their questions via the Support Center for a more comprehensive response from the Technical Support team.