Skip to content

2D beam compression validation, plane stress (SOFA vs FreeFEM) - #37

Merged
th-skam merged 2 commits into
main-freefem-sofafrom
2d-compressionCP
Jul 17, 2026
Merged

2D beam compression validation, plane stress (SOFA vs FreeFEM)#37
th-skam merged 2 commits into
main-freefem-sofafrom
2d-compressionCP

Conversation

@Fimache

@Fimache Fimache commented Jul 17, 2026

Copy link
Copy Markdown

Adds uniaxial compression validation for the 2D cantilever beam
(fully clamped left edge, traction load on the right edge), plane
stress. Same setup as the plane strain counterpart (PR: 2d-compressionDP),

  • Elements: P1 triangles (beam2d_tri.msh, L=1, H=0.2)
  • BCs: left edge fully clamped (ux=uy=0); traction (-q, 0) applied
    on the right edge via consistent nodal forces
  • Plane stress: template="Vec2d" lambda = E*nu/(1-nu^2)` on the FreeFEM side)

@th-skam th-skam left a comment

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Your u_y seems to be asymmetric:

Image

Check that you don't have gravity activated.

@th-skam th-skam left a comment

Copy link
Copy Markdown
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Same as #37

@th-skam
th-skam merged commit 65d0f66 into main-freefem-sofa Jul 17, 2026
1 of 3 checks passed
@th-skam
th-skam deleted the 2d-compressionCP branch July 17, 2026 12:42
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants