描述
一个悬臂梁在左端 (x = 0) 固定不动,右端承受一个横向力 F 和一个轴向力 Fa。抗拉强度为零,压缩行为保持弹性。问题由以下一组参数描述。考虑小变形,在此示例中忽略自重。确定最大挠度 uz,max。
| 材料 | 弹塑性 | 弹性模量 | E | 210000.000 | MPa |
| 泊松比 | ν | 0.000 | - | ||
| 剪切模量 | G | 105000.000 | MPa | ||
| 抗拉塑性强度 | ft | 0.000 | MPa | ||
| 几何形状 | 悬臂梁 | 长度 | L | 2.000 | m |
| 宽度 | w | 0.005 | m | ||
| 厚度 | t | 0.005 | m | ||
| 载荷 | 横向力 | F | 4.000 | N | |
| 轴向力 | Fa | 5000.000 | N | ||
解析解
横向力 F 与轴向力 Fa共同作用使悬臂梁处于弹塑性状态,参见下面的示意图。弹塑性区长度由参数 xp 描述。应力 σx 由弯曲应力 σb 和压缩应力 σc 组成,其根据以下公式定义:
其中 κ(x) 为曲率,参数 z0(x) 的定义满足 σx(x,z0)=0。当固定端顶部表面的弯曲应力达到压缩应力的值时,会发生第一次屈服。
得出横向力 F=2.083 N。因此,横向力 F=4.000 N 使悬臂梁处于弹塑性状态。在弹塑性区 (x < xp),必须满足弯矩和平衡的轴向力之间的平衡。
解这些方程,得到弹塑性区的曲率 κp 和参数 z0。
根据最后一个方程,在条件 z0(xp)=-t/2 下,可以求得弹塑性区长度 xp。
弹性区 (x > xp) 的曲率 κe 按照伯努利-欧拉公式描述:
最后,可以根据以下公式计算最大挠度 uz,max。
RFEM 设置
- 在 RFEM 5.26 和 RFEM 6.11 中建模
- 元素尺寸为 lFE= 0.020 m
- 考虑几何线性分析
- 增量数为 10
- 忽略构件的剪切刚度
结果
| 材料模型 | 解析解 | RFEM 6 | RFEM 5 | ||
| uz,max [m] | uz,max [m] | 比率 [-] | uz,max [m] | 比率 [-] | |
| 各向同性非线性弹性 1D | 1.232 | 1.231 | 0.999 | 1.231 | 0.999 |
| 各向同性砌体 2D | 1.237 | 1.004 | 1.237 | 1.004 | |
| 非线性弹性 2D/3D, Mohr - Coulomb | 1.237 | 1.004 | 1.237 | 1.004 | |
| 非线性弹性 2D/3D, Drucker - Prager | 1.237 | 1.004 | 1.237 | 1.004 | |
| 各向同性塑性 2D/3D, Mohr - Coulomb | 1.237 | 1.004 | 1.237 | 1.004 | |
| 各向同性塑性 2D/3D, Drucker - Prager | 1.237 | 1.004 | 1.236 | 1.003 | |