Skip to content

2D beam MMS: plane stress validation under distributed load (SOFA vs FreeFEM)#35

Merged
th-skam merged 1 commit into
main-freefem-sofafrom
2d-planestress
Jul 17, 2026
Merged

2D beam MMS: plane stress validation under distributed load (SOFA vs FreeFEM)#35
th-skam merged 1 commit into
main-freefem-sofafrom
2d-planestress

Conversation

@Fimache

@Fimache Fimache commented Jul 16, 2026

Copy link
Copy Markdown

Summary

Adds the plane stress counterpart to the existing plane strain distributed-load
validation for the 2D cantilever beam. No analytical solution is available for
this loading case, so validation is performed via cross-comparison between SOFA
and FreeFEM++ on the shared beam2d_tri.msh mesh.

  • Elements: P1 triangles (beam2d_tri.msh)
  • Load case: constant self-weight-type distributed load (0, -q) over the domain
  • Plane stress: obtained on the SOFA side via LinearSmallStrainFEMForceField
    with template="Vec2d" , and on the FreeFEM side via lambda = E*nu / (1 - nu^2) in the weak form.
  • BCs: fixed (Dirichlet) on the left edge (label 1); distributed load applied
    as consistent nodal forces (1/3 per node per triangle, integrated over element area).

Results

| Quantity | RMS (SOFA vs FreeFEM) |
| u_x | 1.006825e-12 |
| u_y | 5.512254e-12 |

@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.

Ok works
Image

@th-skam
th-skam merged commit 67749e5 into main-freefem-sofa Jul 17, 2026
1 of 3 checks passed
@th-skam
th-skam deleted the 2d-planestress branch July 17, 2026 08:28
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