ABAQUS运行到一半中断

持续记录自己在ABAQUS中踩过的坑

1.加载到一半中断：可以用平滑分析步，或者设多几个分析步，逐步加载

2.一开始就直接5U，可以检查一下网格，接触的地方尽量严丝合缝，每个网格对整齐

3.结果中，接触出现穿透，点了接触后可分离，还是黏在一起了，可以在面面接触里面设置容差

4.大荷载的话，一定要多设几个分析步，慢慢加上去。无论是对于收敛还是模拟正确结果，都有帮助

5.小滑移有助于收敛，但是对得出正确结果帮助不大

6.若组合荷载加不上去（有弯矩 水平 垂直等），可以试试单独加，看一下是什么情况

7.算盘坊小胖说，土体开挖不要单独建part，否则难以收敛。当时看不以为然，事到如今经历了才有体会，收敛全靠天意。建议采用分割之后用生死单元。生死单元开挖后，也要在内部施加边界条件或者荷载进行约束，不然不收敛。

8.适当放大网格+平滑分析步，可以避免一开始5U

9.用“创建surface”来确定部件内接触面的时候，正反方向很玄学，错了的话直接反过来

10.过了段时间，对第九点有更深的理解。如果ABAQUS提示是选择边（黄色，深红）与选择壳的内外（棕色，紫色）是不同的。

11.用关键字修改场变量来控制材料参数在不同分析步中的数值。见下图。

12.网格的数量影响精确度，改泥皮宽度时，泥皮对应网格的数量也在变化

13.加密网格时，若不收敛，适当放大底部的网格（桩土模型，轴对称，图中为桩土底部交界处，隐藏了桩）

14.自动化建模提取结果：如何寻找historyRegions？

from odbAccess import openOdb
from abaqusConstants import *
import os

# 指明odb文件所在位置并打开
file_name = 'Job-1.odb'
dir_path = r'E:\ABAQUSmodel\machine_learing'
file_path = os.path.join(dir_path, file_name)
# 获取odb对象
my_odb = openOdb(file_path)
# 获取分析步
step = my_odb.steps['Step-13']
# 获取名称
print(step.historyRegions) 

15.如何在编辑器上运行脚本？pip install -U abqpy==2021.*

16.我最近做了个随机变化桩长桩径的项目，写脚本的时候有一句话很关键，把坐标修改为序号，这样才能方便随机变化桩长桩径session.journalOptions.setValues(replayGeometry=INDEX,recoverGeometry=INDEX)

17.要尽量大的边界，避免边界效应。宽度至少与桩长一样。长100m的取30m宽和100m宽计算结果差别挺大

18.模型简单的话，开多核会降低速度，单核最快

19.出现下列问题：The area of 3 elements is zero, small, or negative. Check coordinates or node numbering, or modify the mesh seed. The elements have been identified in element set ErrElemAreaSmallNegZero.

2 nodes have negative radial coordinate values. Negative values are illegal for axisymetric elements. The nodes have been identified in node set ErrNodeRadCurCoord.

解决方法：检查接触

20.FORCE equilibrium accepted using the alternate tolerance.网上说要检查接触和网格，我是在桩土模型地应力平衡就5U了，是荷载中初始预应力场中横向系数太小导致

21.不收敛，一定要去结果看看云图，有无单元畸变过大？针对性进行修改

22.批量切割三维部件的线，在草图里切割：首先，打开pythonreader，确保你这次草图的model是s1（abaqus不是s就是s1，如果是s，退出草图再进一次就是s1），我需要的是51段

import sketch
import part

s1 = mdb.models['Model-1'].sketches['__profile__']
segment_length = 86.366/51.0  # 约1.693米

for i in range(51):
    x1 = i * segment_length
    x2 = (i + 1) * segment_length
    s1.Line(point1=(x1, 0.0), point2=(x2, 0.0))

23.出现沙漏现象：网格很乱，很尖锐。并且明显变形分层，变形只在局部。细化网格

24.建了个整风机模型，报错：Three factorizations in a row failed. Check the model. It is possible that the model contains the kinematic coupling definition set up in a way that a degree of freedom has neither mass nor stiffness.原因：叶片缺乏约束。我原本是将叶片的一个端点与RNA耦合。报错后，改成整个叶片与RNA耦合

25.整机模态分析，需要开多核计算

26.沙漏网格解决方法：二次单元

27.桩土模型：加了初始应力场了吗？初始应力场似乎与重力荷载不兼容，可以试试关掉重力荷载。文中用的是渐进的E与su，还是平均的？要彻底搞清楚