测试单元:
TMT,TMQ,MITC4;
软件版本:
Simdroid v3.1,Simdroid v4.0;
参考文献:
[1] Sze, Liu, Lo. Popular benchmark problems for geometric nonlinear analysis of shells.
说明:
本次测试共包含8个问题,几何和网格均完全采用Simdroid软件的对应模块建模和剖分。求解器采用Simdroid软件内置的结构求解器(YDFeapSTR)计算。
因为采用粗糙的网格尺寸会导致case-4和case-5结果有较大的偏差,且导致case-6~case-8无法收敛,本次测试采用了足够精细的网格尺寸(与参考文献[1]中比较接近)以获得准确解,并且所有四边形网格均为结构化网格。
case-1~case-7直接采用Simdroid软件的通用静力分析计算完成,case-8采用Simdroid软件的非线性屈曲分析计算完成,并考虑了两种厚度的壳,随着厚度的减半,可以看出壳的几何非线性响应也发生了显著的变化。所有分析采用自动时间增量,因自动时间增量会自适应的调整时间步长,绘制曲线图时采样点可能较少,但并不影响最终的收敛结果的精度。除case8以外,绘制曲线图时荷载全部除以Pmax进行了归一化,曲线纵轴表示荷载正比例系数,曲线横轴为各点对应的位移。
case-7和case-8在文献中同时采用了各向同性材料和复合材料进行测试,因Simdroid软件暂时没有复合材料的求解功能,本次测试只选择各向同性材料进行分析。测试中因riks分析的特性,case-8的两个情况分别限制了最大弧长增量为5000和1000,如不限制可能导致在曲线拐点处因搜索的弧长半径过大而使分析沿原路径返回。
此外,考虑到文献中的单位为无量纲单位,本次测试时人为的将所有几何、材料和荷载参数进行了无量纲处理,比较后处理云图时可以忽略单位,其中case-8的Simdroid结果应放大1000倍后与文献进行比较。
因为分析均为几何非线性分析,所有后处理云图的变形放大系数为1,未进行任何放大,图中所现变形即为真实变形。
测试结论:
每个case均分别采用TMT、TMQ和MITC4三种单元求解,所有单元均能通过所有8个测试,且所有单元的计算结果与参考文献[1]的结果一致,其中参考文献[1]的结果均由某商业软件软件的S4R单元计算得到,此处用作与Simdroid软件的参考对比。
1. Cantilever subjected to end shear force
文献[1]中的模型和结果如下:
Simdroid软件计算的结果如下:
Simdroid软件计算的TMT单元的变形云图为:
文献[1]中的模型和结果如下:
Simdroid软件计算的结果如下:
Simdroid软件计算的TMQ单元(加密网格后)的变形云图为:
文献[1]中的模型和结果如下:
Simdroid软件计算的结果如下:
Simdroid软件计算的TMT单元的变形云图为:
Simdroid软件计算的MITC4单元的变形云图为:
4. Hemispherical shell subjected to alternating radial forces
文献[1]中的模型和结果如下:
Simdroid软件计算的结果如下:
Simdroid软件计算的TMQ单元的变形云图为:
Simdroid软件计算的TMT单元的变形云图为:
5. Pullout of an open-ended cylindrical shell
文献[1]中的模型和结果如下:
Simdroid 软件计算的结果如下:
Simdroid 软件计算的TMT单元的变形云图为:
文献[1]中的模型和结果如下:
Simdroid软件计算的结果如下:
Simdroid软件计算的MITC4单元的变形云图为:
Simdroid软件计算的TMT单元的变形云图为:
文献[1]中的本例还包括复合材料的分层壳测试,Simdroid软件暂不支持复合材料,故只选择各向同性材料进行测试。文献[1]中的模型和结果如下:
Simdroid软件各向同性材料的计算的结果如下:
Simdroid软件计算的TMT单元的变形云图为:
文献[1]中的模型如下:
文献[1]中壳厚度为12.7的荷载位移曲线如下:
因为单位转换的问题,Simdroid软件计算时采用的模型尺寸为文献中的1/1000,弹性模量为文献[1]中的1e6倍,计算得到的厚度和位移值均为文献的1/1000,Simdroid软件采用0.0127厚度的壳计算结果如下:
文献[1]中壳厚度为6.35的荷载位移曲线如下:
Simdroid软件采用0.00635厚度的壳计算结果如下:
Simdroid软件计算的TMQ单元(厚度为0.0127)的变形云图为:
Simdroid软件计算的TMT单元(厚度为0.00635)的变形云图为: