abaqus隐式结果导入显式+预应力弹片跌落+变形网格导出

先跑隐式，standard算弹片装配压紧的过程。弹片被压弯或者卡到某个位置，这个状态下应力分布不均匀，靠近根部的地方已经进塑性了，但整体还没到极限。分析步用static general，打开几何非线性NLGEOM，打开重启动请求，restart requests频率设1或者设一个具体的增量步号，后面import靠的就是这个。输出里面勾选stress、strain、displacement、reaction force这些基本量，如果需要接触力也把CFORCE加上。隐式算完之后job manager里会生成重启动文件，.res和.abq之类的文件全在工作目录里放着别动。

进入第二个模型，import要手动建。Load模块里predefined field → create，step选initial，category选other，类型选Initial state，然后选需要继承状态的部件，填job name、step number和frame number。job name就是隐式分析那个job的名字，不带扩展名。step填Last或者具体数字，frame一般填Last。这一步做完模型树里会出现一个黄色圆点符号表示初始场已经挂上去了，但此时网格看着没变，正常现象，abaqus是提交分析的时候才去读那个odb。文件路径这里有个容易翻车的地方：隐式分析的工作目录和显式分析必须一致，或者至少隐式的.odb、.res这些文件要拷到显式分析的目录下。之前有一次把odb留在上一级目录，导入的时候直接报找不到job文件。

Part名字必须对得上。用file → import → part从隐式odb里导入变形后的网格时，abaqus会自动在part名后面加“-1”，比如原来叫bimetal_spring，导入后变成bimetal_spring-1。这个后缀必须删掉，改成跟原始分析完全一样的名字，否则predefined field那里选不到，或者选上了也传不过去。导入类型选deformed shape，不要选undeformed。导入之后这个part就是带着变形形状的新part，不能再给它赋新的截面属性或者材料方向，所有和材料相关的东西已经跟着initial state一起锁死了。

ALE那块在step模块里设。显式分析步建好之后，other → ALE adaptive mesh domain → manager，先create一个domain，region选弹片主体，不要整个模型都选，只选可能发生大变形或者网格扭曲的区域。frequency设10或者20，remeshing sweeps per increment一般填1就行。ALE adaptive mesh controls里面，initial remeshing sweeps保持默认，mesh motion选默认的，不要随便改。弹片本身如果材料是超弹性或者泡沫，官方文档不推荐用ALE，但普通弹塑性金属没问题。另外注意一个细节：如果后面还要用distortion control来防止单元负体积，那ALE域里的单元不能同时开启distortion control，二者是互斥的。

显式跌落分析步，time period根据跌落高度定，比如从1米高度自由落体，到接触地面大概0.45秒左右，但显式里如果把时间设得太长计算量爆炸。一般设个0.01到0.02秒就够捕捉碰撞过程了，具体看弹片尺寸和接触时长。重力通过load模块的gravity施加，注意方向。初始速度用predefined field里的velocity，比如沿-Z方向给个-4429 mm/s，对应1米跌落。接触用general contact，property里给切向摩擦和法向硬接触，法向可以加一点过盈容差。弹片如果跟其他部件有预紧接触关系，接触定义要跟隐式里保持一致，不然导入之后接触状态会跳变。

弹片跌落的边界条件尽量少加，除了地面刚体固定，其他部件不要加多余的约束，让它自然飞。如果模型里有螺栓或者铆钉连接，用connector或者tie绑住就行，不需要单独建螺栓实体再导应力，那种做法是给电池包或者大型装配体用的，弹片这种薄片件没必要。输出设置里勾选energy，ALLKE和ALLIE对比一下，如果动能占比突然飙升超过内能的5%到10%，说明有数值不稳定，回去看是不是接触穿透或者单元畸变太严重了。

最后跑的时候如果报“Excessive distortion”那个经典错误，先别急着重画网格，检查ALE域有没有覆盖到畸变最严重的位置。弹片跌落到地上的时候根部圆角处最容易压溃，ALE域region要把那个圆角包含进去。如果ALE开了还是报错，把remeshing frequency从10降到5试试，让remesh更频繁一些。job提交之后看sta文件，时间增量步如果掉到1e-8以下，说明有单元尺寸太小或者材料波速太高，检查一下网格有没有过密的局部区域。变形网格导出这块，显式分析跑完如果想拿最终变形后的网格去重新建模，可以用file → import → part从显式odb里导deformed shape，同样记得删掉part名后缀。

工作目录里最终应该有的文件：隐式job的.odb/.res/.prt/.mdl，显式job的.inp/.odb/.sta/.msg/.res。如果中间要改弹片的几何尺寸重新跑隐式，那前面显式的所有initial state定义都要重新指向新的job name，frame号也要重新确认，因为增量步变了之后frame的编号也变了。参数写到这里，下次试一下frame填具体数字而不是Last看看有没有区别，之前一直用的Last，万一Last对应的是重启动文件里最后一个可用的而不是odb里最后

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删