先说说单位制,mm-tonne-s这套,钢密度7.85e-9,杨氏210000,重力9810,声速343000,空气密度1.225e-12,体积模量0.144,这套记住了后面省很多事,SI制和m-tonne-s也顺手记一下,SI的密度7850,杨氏2.1e11,重力9.81,m-tonne-s的密度7.85,杨氏2.1e8,力用kN,应力用kPa,换算脚本网上下一个,长度乘1000质量除1000密度乘1e-12,自己写也行,D:\abaqus\unit_table.xlsx存着。
刚开始别急着建复杂模型,拉伸一个方块,一端固定一端加100N,C3D8R网格,静力学跑通,看应力云图和位移,确认standard.exe在跑,odb里field output有数据,这个过程就是验证solver装没装,装的时候典型安装只勾CAE和Viewer,Solver在自定义里,results can not be found九成是这个问题。
静力学收敛是第一个坎,msg文件里zero pivot,数值奇异的先看约束够不够,刚体位移没锁死,再看接触,tie主从面距离太远容差调到0.5还不行就检查法向,接触定义general contact摩擦给0.1法向硬接触过盈容差0.01,增量步初始给0.01最小1e-9最大0.1,不收敛自动折半,打开几何非线性NLGEOM,大变形不开算出来是错的。
非线性材料分几类,弹塑性金属用plastic,屈服应力塑性应变表,超弹性橡胶用hyperelastic,neo-hookean或者mooney-rivlin,参数从单轴拉伸曲线拟合,混凝土用CDP,损伤因子表格单调加载下可以留空但软化段刚度不折减,反复荷载必须填,dt和dc从0到1,非弹性应变扣掉峰值弹性应变,刚度恢复因子wt=0,wc=1默认,热膨胀1.2e-5钢,混凝土1e-5,橡胶2e-4。

UP住这个会不会悬臂梁施加二次函数载荷。
动力学用frequency提取模态,number of eigenvalues给20,最小频率0最大20000,看dat里耦合频率,buckling做特征值屈曲,explicit做跌落碰撞,time period设0.01到0.02秒,重力9810方向-Z,初始速度-4429mm/s对应1米跌落,接触general contact,ALE自适应网格region选可能大变形的区域,frequency设10,mesh motion默认,distortion control和ALE互斥,sta文件时间增量掉到1e-8以下说明网格太密或者波速太高。
热分析mm-tonne-s下导热mW/(mm·K),比热mJ/(tonne·K),热膨胀1/K,热流mW/mm²,辐射常数5.67e-14mW/(mm²·K⁴),顺序耦合先算温度场再导入应力,声学AC3D4单元,声速343000,无限元ACIN3D4包外围,PML在step里other→PML创建指定声学域,耦合面tie容差0.1。
子程序用VUMAT写本构,UMAT写塑性,DLOAD写移动载荷,命令行abaqus cae noGUI=script.py传参,vf=0.35 radius_min=3 radius_max=8,random.seed固定值做参数化对照,工作目录D:\temp\abaqus_work里inp odb sta msg res全留着,重装软件不用重新建模,直接提交inp。
理论书翻翻塑性力学和有限元基础,应力更新算法return mapping,一致切线刚度,收敛判据L2范数1e-5,接触力学法向约束拉格朗日乘子,显式中心差分稳定步长跟最小单元尺寸和波速挂钩,声固耦合Helmholtz方程,这些看懂了再看帮助文档。
免责声明:本文系网络转载或改编,未找到原创作者,版权归原作者所有。如涉及版权,请联系删