项目介绍
实验中固定一块薄板,左端固定并施加均布压力。 这里考虑了小变形,并且自重忽略不计。 下面的参数集描述了该问题。 计算最大挠度 uz,max 。
材料 | 弹塑性 | 弹性模量 | E | 210000.000 | MPa |
泊松比 | ν | 0.000 | - | ||
剪切模量 | G | 105000.000 | MPa | ||
屈服强度 | fy | 40,000 | MPa | ||
几何尺寸 | 墙板 | 周长 | L | 1.000 | m |
宽度 | w | 0.050 | m | ||
厚度 | t | 0.005 | m | ||
荷载 | 均匀压力 | p | 2.750 | kPa |
解析解
讨论了荷载的数量。 首次屈服时的弯矩 Me和塑性铰发展时的极限弯矩 Mp计算如下:
板在压力 p 作用下进入弹塑性状态。弯曲应力按照下面的公式定义:
式中 ka 是曲率。 弹塑性区的长度用参数xp描述。 面的弯曲应力等于面上的塑性强度 fy在点 xp ,见下面的示意图。
弹塑性弯矩 Mep (内力)必须等于弯矩 M(外力)。 弹塑性区中的曲率 kp 。
RFEM 设置
- 在 RFEM 5.26 和 RFEM 6.01 中建模
- 单元尺寸 lFE =0.020 m
- 实体模型沿厚度方向进行网格细化(每个厚度6个单元)
- 考虑几何线性分析
- 增量数目为 5
- 忽略杆件的抗剪刚度
结果
模型 | 解析解 | RFEM 5 | RFEM 6 | ||
uz,max [mm] | uz,max [mm] | 比值 [-] | uz,max [mm] | 比值 [-] | |
一维各向同性塑性 | 166,234 | 166.214 | 1.000 | 166.018 | 0.999 |
二维/三维各向同性塑性,板 | 162.987 | 0.980 | 162.960 | 0.980 | |
二维/三维各向同性非线性弹性,板,von Mises | 165.730 | 0.997 | 165.700 | 0.997 | |
二维/三维,板,Tresca | 166.998 | 1.005 | 166.969 | 1.004 | |
二维/三维各向同性塑性,实体 | 160.601 | 0.966 | 162.429 | 0.977 | |
各向同性非线性弹性二维/三维,实体,von Mises | 163.003 | 0.981 | 165.593 | 0.996 | |
各向同性非线性弹性 二维/三维, Solid, Tresca | 168.725 | 1.015 | 169.691 | 1.021 | |
一维各向同性非线性弹性 | 166.214 | 1.000 | 166.018 | 0.999 |
注释: 解析得出的扭转常数和数值计算得出的扭转常数也是导致结果偏差的原因。