SIGINI子程序是一类比较冷门的子程序,因此,公开能查询到的资料不多。它的作用主要是定义材料积分点上的初始预应力,初始预应力场可以定义为积分点坐标、单元号、积分点数等信息的函数。一般最常见的就是预应力场随着空间位置的变化。 这类子程序的接口如下: 其中必须要定义的变量为SIGMA数组,其维度及存储顺序取决于单元类型,NTENS即所用单元的应力张量中直接分量与间接分量之和,以C3D8R单元为例,NTENS=6,即SIGMA数组中有六个元素: SIGMA(1)为第一个应力分量,即S11 SIGMA(2)为第一个应力分量,即S22 SIGMA(3)为第一个应力分量,即S33 SIGMA(4)为第一个应力分量,即S12 SIGMA(5)为第一个应力分量,即S13 SIGMA(6)为第一个应力分量,即S23 另一个重要变量是COORDS,它代表的是积分点的初始位置,NCRDS一个积分点的坐标点个数;其他变量还有单元编号NOEL、单个单元中积分点的数量NPT、复合材料单元中铺层的数量LAYER、复合材料单元中单个铺层截面点的数量。 采用SIGINI子程序定义预应力场时,Abaqus/CAE无法直接操作,需要修改关键字,可以在Model→Edit Keywords中选择对应的模型对应的位置进行命令添加,也可以在inp文件中直接修改。 关键字中添加的命令如下,注意其修改位置: 然后就可以在子程序中去控制预应力的变化了,以下是示例代码: 另外,需要注意的是,施加的预应力场要和初始状态的外部载荷相互平衡,否则在初始预应力的作用下,结构会发生变形,关于内应力平衡的问题本文不做详述。 案例 — 本案例旨在说明如何在CAE中设置SIGINI子程序的调用以及如何编写SIGINI子程序,模型及源程序代码仅供参考,无特定物理意义,且模型中未添加外力平衡,仅供参考。 |
|