报错弹窗又弹出来了,negative eigenvalues,明明加载速率已经压到0.01mm/s了还是算不过去,上次调这个模型花了整整三天,那天下午办公室空调坏了闷得像蒸笼,旁边工位老张一直在啃辣条那个味道混着汗味飘过来,我盯着屏幕上红彤彤的abort提示心里烦躁得不行,鼠标点关闭都带着气,后来发现根本不是速率的问题而是接触定义里忘了勾选adjust to remove overclosure,初始穿透直接导致第一步就发散,这种低级错误在赶项目的时候特别容易犯,脑子一热就把所有参数都怀疑一遍唯独漏了最基础的几何检查。
Step模块里把Nlgeom打开,准静态压缩不开几何非线性等于白算,大变形下网格畸变严重的话还得开自适应网格重划分,ALE那个选项藏在Mesh模块下面不太起眼但关键时刻能救命,之前用Ansys做类似工况时默认就开了大变形所以没踩过这个坑,换到Abaqus反而要手动设置容易忘。
载荷施加别直接用Displacement/Rotation边界条件硬压,改用Smooth Step幅值曲线让位移从零平滑过渡到目标值,避免初始瞬间的冲击效应引发高频振荡,动能占比超过内能5%就算不上准静态了,History Output里一定要监控ALLKE和ALLIE的比值,这条曲线要是前几个增量步就飙升后面再怎么调都是白费功夫。
接触属性里的摩擦系数别随便填0.3,实测数据和经验公式差很多,钢对钢干摩擦0.15-0.2之间比较合理,润滑状态下可能低到0.05,罚函数摩擦比拉格朗日乘子法收敛性好但精度稍差,如果模型里有多个接触对优先把主面设成刚度大网格粗的那个面,从面选软材料细网格,反过来定义的话穿透量会大得离谱。
输出请求别贪多,Field Output里只留应力应变和接触状态就够了,历史输出按需添加,全开的话dmp文件膨胀到几十G后处理卡得要命,尤其是显式动力学求解器跑准静态时节点数乘以时间步数再乘以输出变量数量,硬盘读写直接成为瓶颈,有次提交作业后等了六个小时才发现结果文件写满了临时盘空间导致任务中断,从那以后养成了先估算文件大小再提交的习惯。
网格类型选C3D8R就行,减缩积分单元抗剪切自锁性能好适合压缩工况,完全积分单元C3D8在塑性大变形下容易出现体积锁定导致应力虚高,沙漏控制用Enhanced或者Relax Stiffness,默认的Stiffness有时候压制过度反而扭曲真实响应,Hourglass Energy占Internal Energy比例控制在1%以内才算合格,超了就得加密网格或者换单元公式。
材料本构别直接用弹性+理想塑性凑合,真实应力应变曲线至少拉伸到颈缩之后,压缩试验数据如果有最好单独定义Compression Only的塑性段,金属板材各向异性明显的话还得加Hill48或者Barlat屈服准则, isotropic hardening在循环加载下包辛格效应捕捉不到,kinematic hardening参数标定又麻烦,看项目精度要求取舍吧。

各位彦祖大佬们,我用abaqus进行准静态压缩仿真,显式动力学,时间为250,压缩位移为50,这速度为0.2mm/s,但是动内能之比从大概28s才小于5%,也就是大概压了5mm多动内能之比才小于5%,我的压头和样品是没有缝隙的。
我问AI,它说可以认为是准静态,但我感觉有点不严谨啊,请问这可以算是准静态压缩吗?
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删