引言：iSolver为一个完全自主的通用结构有限元软件，对标国际主流结构CAE商业软件Abaqus、Ansys、Nastran，支持结构分析的常用功能，精度、效率和Abaqus相当， iSolver自带友好的三维可视化前后处理界面，也可作为一个轻量化插件集成到Abaqus/FEMAP或者自主软件中。本帖以地基中应力波的传播为例，将iSolver和Abaqus计算结果进行对比，比对两种有限元软件的计算结果。
问题描述：
如图1所示，一地基表面中作用有一冲击荷载，荷载幅值曲线见图2所示。地基动弹性模量为720MPa,泊松比为0.33,密度为2.842t/m3,对应纵波波速为616.91m/s,剪切波速310.75m/s。
分析中利用对称性取一半进行分析，分析区域宽度为80m,高度为50m,荷载作用宽度为32.5m,大小在幅值曲线的基础上乘以100kPa,为清晰辨识波的传播，分析中不考虑阻尼。
图1 计算图示
图2 荷载幅值曲线
1 Abaqus操作及计算结果
（1）建立几何模型
考虑结构及荷载的对称性，建立一半的几何模型，长为80m，宽度为50m，如图3所示。为方便荷载施加，将左上顶点偏移(35.0,0,0)建立参考点，并切分几何体。
图3 几何模型
（2）材料及截面赋予
使用线弹性材料本构，地基弹性模量720MPa,泊松比为0.33,密度为2.842t/m3。创建solid,homogeneous截面并赋予地基。
（3）创建分析步
建立名为pulse的Dynamic,Implicit分析步，在EditStep对话框的Basic选项卡中将时间总长设为0.13，Incrementation选项卡中选择步长控制方去类型（Type）为Fixed，将增量步步长设为0.001，最大允许增量数设为400，在Other选项卡中将数值积分算法中的Alpha设为0,取消数值阻尼，接受其余默认选项后退出。
图4 分析步设置
（4）荷载及边界条件
图5 边界与荷载设置
在Load模块中，执行【BC】/【Create】命令，在initial(初始）分析步中约束模型左侧（对称轴）上的位移U1,其余边界均不约束，意味着在动力荷载下为自由边界。
从文件中读取数据，创建图2所示的幅值曲线。在pulse分析步中对所给区域施加荷载，在EditLoad对话框中将Magnitude设为100,在Amplitude下拉列表中选择刚才所创建的幅值曲线。
（5）网格划分
图6 网格划分
（6）作业提交
Abaqus创建并提交作业，开始分析，
图 7 Abaqus作业提交
（7）计算结果
图 8 Abaqus位移U1云图（放大10000倍，t=0.08s)
图9网格变形图（放大10000倍，t=0.13s)
2 iSolver操作及计算结果
（1）导入模型
打开iSolver，点击File->Importing->Model，
图 10 iSolver界面
选择Abaqus生成的inp并确定，
图 11 iSolver导入模型
此时iSolver会显示该模型，
图 12 iSolver显示模型
（2）作业提交
iSolver切换到Job模块，点击Job->Create…
图 13 iSolver创建作业
输入作业名称、选择已导入的模型，点击OK
图 14 iSolver设置作业
点击Job->Manager…，打开作业管理器，然后选择上述作业，然后点击Submit，iSolver将在后台进行分析，
图 15 iSolver作业提交
（3）计算结果
iSolver计算结果如下
图 16 iSolver位移U1云图（放大10000倍，t=0.08s)
图 17 iSolver位移U1云图（放大10000倍，t=0.13s)
模型界面上的波动荷载将在分析区域中产生纵波和剪切波，同时在荷载边缘处由于荷载突变产生面波也会向左、右两侧传递。这些波型在变形后的网格中清晰可见。根据所给数据，纵波约0.08s到达底部边界，此时的网格如图 8图 16所示。之后，由于波在分析区域边界上发生反射（包括纵波、剪切波和面波）,重新回到分析区域内部，波的相互叠加在图 9图 17中清晰可见。
iSolver各增量步云图动画如下，
(参考：费康, 彭劼《ABAQUS岩土工程实例详解》)
============================================================
软件免费下载地址
iSolver为免费软件，且无license限制，最新版免费下载地址如下：https://www.jishulink.com/content/post/337351
微信
QQ
微博