Circular Inclusion Problem - GCMLab/GCMLab-FEM GitHub Wiki

Parameter | Material 1 | Material 2 |
---|---|---|
Young's Modulus | 1 | 10 |
Poisson's ratio | 0.25 | 0.3 |
The body is composed of two different materials.
Plane strain condition is applied.
A radial uniform displacement of 15 is applied to the outer boundary of the circular plate.
Symmetry conditions are applied so that only one quarter is modeled.
The problem can be solved analytically uder plane strain conditions leading the radial and hoop strains in polar coordinates [1]:
where
The FEM model is solved by a fine mesh. Then, FEM solution () is compared with analytical solution.
As you can see, there is a good agreement between FEM result and analytical one.
![]() |
![]() |
---|
Note: In current version of FEM code, we have replaced this fine mesh with the coarse one. In the below figure, the coarse mesh is depicted.

[1] Sukumar, Natarajan, David L. Chopp, Nicolas Moës, and Ted Belytschko. "Modeling holes and inclusions by level sets in the extended finite-element method." Computer methods in applied mechanics and engineering 190, no. 46-47 (2001): 6183-6200.