This interactive simulation demonstrates 2D Finite Element Analysis (FEA) using Linear Triangle (T3 / Constant Strain Triangle) elements in plane stress. A rectangular mesh is built from nodes and triangular elements; the global stiffness matrix is assembled from each element’s B (strain–displacement) and D (constitutive) matrices. Fixed nodes (left edge by default) are enforced with the penalty method; the governing equation K d = F_ext + F_p is solved with the Conjugate Gradient method. Displacements are exaggerated for visualization; triangles are colored by Von Mises stress (blue = low, red = high).

Adjust Young’s Modulus E and Poisson’s ratio ν in the sidebar. Use Apply force (drag) to drag a node and set force; use Toggle fix node to fix or free a node. Use Mesh Density to change grid size; Interpolation for None/Linear/Shape Functions; Principal Stresses, Nodal Influence, Residuals for extra visuals. The Math Monitor shows K, d, F_ext, F_p; click K for assembly animation. Reset restores mesh and left-edge fix. The gray outline is the original mesh; the colored mesh is the deformed configuration (displacement scaled for visibility).

 

Math behind the simulation

1. Goal

Solve for nodal displacements d given applied forces and boundary conditions: K d = F_ext + F_p, where K is the global stiffness matrix, F_ext is external force, and F_p is plastic force.

2. Plane stress and D matrix

For plane stress, stress and strain are related by σ = D ε. With Young’s modulus E and Poisson’s ratio ν, D = E/(1−ν²) · [1, ν, 0; ν, 1, 0; 0, 0, (1−ν)/2].

3. Constant Strain Triangle (CST)

Each triangle has linear displacement; strain is constant in the element. The B matrix (strain–displacement) depends on the node coordinates and area: ε = B u_e where u_e is the 6×1 vector of nodal displacements (u,v for each of 3 nodes).

4. Element stiffness and assembly

Local stiffness k_e = B^T D B A t, where A is area and t is thickness (1 here). k_e is assembled into the global K by mapping local DOFs to global node indices. Click K in the equation box to animate this assembly.

5. Boundary conditions and solver

Fixed degrees of freedom are enforced with the penalty method: a large value is added to the diagonal of K for that DOF and the force is set to zero. The system is solved with the Conjugate Gradient method (iterative; scales with mesh size).

6. Von Mises stress and plasticity

For each element, strain ε = B u_e, stress σ = D ε. Von Mises (plane stress): σ_v = √(σ_x² + σ_y² − σ_xσ_y + 3τ_xy²). Colors map low (blue) to high (red) stress. With plasticity enabled, excess stress accumulates permanent strain and F_p is added to the right-hand side.

 

Usage

1. Apply force: Set mode to “Apply force (drag)”, then click and drag a node. The drag vector sets the force at that node; the mesh deforms and stress updates.
2. Fix / free nodes: Set mode to “Toggle fix node”, then click a node to fix or free it in X and Y. By default the left edge (x = 0) is fixed.
3. Material: Change Young’s Modulus E and Poisson’s ratio ν. Use Mesh Density (4–30) to change grid size (e.g. 10×5).
4. Plasticity: Enable and set Yield strength; drag hard to see permanent deformation. Clear Plasticity resets.
5. Interpolation: None (flat), Linear (Gouraud), or Shape Functions (with sampling markers).
6. Math Monitor: Click K in the equation to run assembly animation. Hover the K grid for K[i,j]; hover a node (with Nodal Influence on) to see its DOF rows.
7. Reset: Restores mesh and left-edge fix; keeps current Mesh Density; clears applied forces.

Parameters

• Young’s Modulus E: Stiffness (e.g. 1e6). Higher E ⇒ smaller displacements for the same force.
• Poisson’s ratio ν: Lateral strain / axial strain (0 < ν < 0.5). Larger ν ⇒ more “bulge” perpendicular to load.
• Mesh Density: Grid size (e.g. 10 → 10×5). Range 4–30; finer mesh ⇒ more DOFs and larger K matrix.
• Plasticity / Yield strength: When enabled, stress above yield accumulates permanent strain; F_p appears in the equation monitor.
• Interpolation: None (flat), Linear (Gouraud), or Shape Functions (sampling markers).

Visualization

• Gray outline: Undeformed mesh.
• Colored triangles: Deformed mesh; color = Von Mises stress (blue low → red high). Interpolation controls smooth vs. flat shading.
• Nodes: Blue = unloaded; gray = fixed; red = has applied force.
• Principal Stresses: Green/red crosses on elements (tension/compression).
• Nodal Influence: Hover a node to highlight connected elements and its global DOF rows.
• Residuals: Red arrows showing F − K d at each node.
• Math Monitor: Governing equation, dimensions, sparsity %, and K×d=F_ext+F_p canvases (K sparsity grid and d, F_ext, F_p vector ribbons).