基于ABAQUS平台的截面特性计算PYTHON脚本

在构件计算过程中，不可避免需要计算截面特性，常见的特性值计算可采用CAD/MIDAS/ANSYS等软件计算。但有时我们需计算一些不常见的截面特性值，如截面不对称系数，这就带来一些困难，因为常见方法的计算结果中并未给出这些值。对于常见的形状规则的截面，我们可以根据公式进行手算积分计算，但对于形状较为复杂的截面，我们难以手算。此时，采用数值计算方法显得非常重要。常见的方法为：将截面离散为若干单元，将理论积分公式离散为各单元数值之和，如硕士阶段学过的条带法。

帖子内容是基于ABAQUS平台编写PYTHON脚本，以计算所需的截面特性值。

具体思路如下：

（1）将绘制的截面形状以IGS格式（也可以是其他能导入ABAQUS的格式）导入ABAQUS中；

（2）以导入的截面形状为草图，在PART中建立壳部件；

（3）装配并划分网格，以离散截面为若干三角形单元（划分三角形单元的目的是适应复杂截面的网格划分）；

（4）读取部件单元节点坐标；

（5）调用截面特性计算函数，以计算形心坐标；

（6）根据所计算的形心坐标移动部件，使得坐标原点位于截面几何形心；

（6）调用截面特性计算函数，以计算所需截面特性。

程序实现及各步骤解释如下图所示。

完整代码如下：