[PYTHON] Stress analysis of torus under internal pressure using axisymmetric stress analysis program
This is an application example of an axisymmetric stress analysis program created in Python.
Theoretical solution
A torus (toroidal shell) is a donut-shaped solid in which the annulus is rotated around the central axis (here, the vertical axis), as shown on the left in the figure below.
Consider a torus that receives an equal internal pressure p, as shown on the right of the above figure. In this case, the circumferential axial force $ N_ {\ varphi} $ with radius a and its axial axial force $ N_ {\ theta} $ are given by Timoshenko (Theory of Plates and Shells) by the following equation. ..
From the above equation, it can be seen that $ N_ {\ varphi} $ changes depending on the location, but $ N_ {\ theta} $ takes a uniform value on the circle with radius a. Now, if you rewrite $ N_ {\ varphi} $ at the representative point,
It can be seen that the stress in the circumferential direction inside the annulus, that is, on the rotation axis side ($ r_0 = b --a $) increases. In addition, the circumferential stress at the top and bottom of the annulus with radius a is the same as the circumferential axial force of a straight annulus with radius a that receives normal uniform internal water pressure.
Analysis by axisymmetric FEM
The generated stress of the torus shown below was predicted by axisymmetric FEM analysis. The analysis conditions are as follows. Here, t indicates the plate thickness of the material that makes up the torus.
| E (MPa) | po | p (MPa) | a (mm) | b (mm) | t (mm) |
|---|---|---|---|---|---|
| 200,000 | 0.3 | 1.0 | 2,000 | 4,000 | 10 |
| Number of elements: 360 (a circle with radius a is divided by a central angle of 1 degree) | |||||
In the axisymmetric FEM, the load is input per rad with respect to the rotation axis. Note that water pressure is no exception, and it is necessary to create and input loads faithfully to the definition.
In the stress distribution map at the bottom right of the above figure, For the angle of the horizontal axis, refer to the displacement diagram. The outermost point (horizontal left end) of the torus is set to zero, the top is 90 degrees, the innermost point on the rotation axis side is 180 degrees, and the bottom is 270 degrees.
The table below shows a comparison between the analysis results and Timoshenko's solution.
| Location | FEM | Timoshenko | ||
|---|---|---|---|---|
| (φ) | σφ | σθ | σφ | σθ |
| 0 (r=b+a) | 166.59 | 99.54 | 166.66 | 100.00 |
| 90 (r=b) | 198.43 | 105.62 | 200.00 | 100.00 |
| 180 (r=b-a) | 300.54 | 99.55 | 300.00 | 100.00 |
| The unit of stress is MPa | ||||
program
The program is shown by a link to Gist.
- Axisymmetric FEM analysis program (py_fem_asnt.py)
- Data creation program for drawing (py_fig_data.py)
- Displacement mode / stress distribution map creation script by GMT (a_gmt_gra.txt)
- Torus model creation script (a_gmt_model.txt)
- Input data for torus analysis (inp_torus.txt)
- Torus drawing script by gnuplot (gpl_torus.txt)
that's all