Thermal Deformations of a 3-D Brick

Let us consider a block. A length of the block is L, a width is b, a height is h.
Sought quantities are the absolute elongations of the block along axes X,Y,Z because of a temperature change.
Let us use the following initial data: L = 0.3 m, b = 0.2 m, h = 0.1 m.
Material characteristics: the Young's modulus E = 2.1E+011 Pa, Poisson's ratio ν = 0.28, linear expansion coefficient α = 1.3E-005 K -1.

The temperature change ΔT is 100o .

The finite element model with applied loads and restraints

The analytical solutions are calculated by the formulas:
Δx = α L ΔT
Δy = α b ΔT
Δz = α h ΔT

Thus,
Δx = 3.90000000E-004 m
Δy = 2.60000000E-004 m
Δz = 1.30000000E-004 m .

After carrying out calculation with the help of AutoFEM, the following results are obtained:

Table 1.Parameters of the finite element mesh

Finite Element Type

Number of Nodes

Number of Finite Elements

quadratic tetrahedron

1367

6282

Table 2. Result "Displacement"

Numerical Solution
Displacement Δ*, m

Analytical Solution
Displacement Δ, m

Error δ =100%*|Δ* - Δ |/| Δ |

3.90000147E-004

3.90000000E-004

0.38E-004

2.60000059E-004

2.60000000E-004

0.23E-004

1.30000117E-004

1.30000000E-004

0.90E-004

 

Conclusions:

The relative error of the numerical solution compared to the analytical solution for displacements is equal to 0.0001% for quadratic finite elements.

*The results of numerical tests depend on the finite element mesh and may differ slightly from those given in the table.

 

Read more about AutoFEM Static Analysis

autofem.com

Return to contents