描述
一个薄壁球形容器承受内压 p。忽略自重时,目标是确定容器的 von Mises 应力 σMises 和径向挠度 uR。
| 材料 | 各向同性线弹性 | 弹性模量 | E | 210000.0 | MPa |
| 泊松比 | ν | 0.296 | - | ||
| 几何 | 球壳 | 半径 | R | 0.500 | m |
| 壳厚 | t | 5.000 | mm | ||
| 荷载 | 压力 | 内压 | p | 5.000 | MPa |
解析解
解析解基于薄壁容器理论。该理论假定壳体处于膜应力状态;因此,必须满足以下条件:
- 壳体厚度不能发生不连续变化。
- 分布荷载不能发生不连续变化。
- 曲率半径和圆心位置不能发生不连续变化。
- 包括反力在内的外力必须与壳表面相切。
应力状态由拉普拉斯方程描述:
其中 σ1、σ2 分别为子午方向和环向应力,R1、R1 为相应方向上的半径。对于球形容器,由于对称性(σ1 = σ2 = σ,R1 = R2 = R),可简化为:
von Mises 应力 σMises 可由主应力确定:
容器的径向挠度 uR 由胡克定律得出:
RFEM 设置
- 建模于 RFEM 5.39 和 RFEM 6.13
- 单元尺寸 lFE = 0.010 m
- 使用各向同性线弹性材料
- 使用完整模型和八分之一模型
结果
| 计算模型 | 理论 σMises [MPa] |
RFEM 6 σMises [MPa] |
比值 [-] |
RFEM 5 σMises [MPa] |
比值 [-] |
| 完整模型 | 250.000 | 249.984 | 1.000 | 249.987 | 1.000 |
| 八分之一模型 | 249.984 | 1.000 | 249.984 | 1.000 |
针对这两个计算模型,von Mises 应力均在测试点 A 处读取。由于网格拓扑,壳体表面存在轻微的应力偏差。
| 计算模型 | 理论 uR [mm] |
RFEM 6 uR [mm] |
比值 [-] |
RFEM 5 uR [mm] |
比值 [-] |
| 完整模型 | 0.419 | 0.419 | 1.000 | 0.419 | 1.000 |
| 八分之一模型 | 0.419 | 1.000 | 0.419 | 1.000 |