附录  程序中的变量说明

输入文件名：curve.dat，单位采用10³kg、m、s，去掉注释行以后存为fort.16。

注释行前四个字符为“C***”或“C”，终结行前四个字符为“STOP”。

以下标◎的栏目未做到前处理界面中，采用默认值的处理办法。

以下标·的栏目为通过前处理界面的给定值计算得到，不必另行输入。

以下标＃的栏目为前处理界面中要求输入的内容，同时附有“描述”项。

Ⅰ.计算控制参数

CARD   1        无注释行

#Columns    1-80    Fortran Name      TITLE

描述：项目主题。

长度限制为80字符，作为对该计算问题的文字描述，并在输出中作为标识。

如该行内容为“RESTART”，则程序读取fort.13中存取的断点信息继续上次未完成的计算。

CARD   2      注释行：GENERALDATA

·Columns    l-5    Fortran Name      NNODE

模型中的总节点数，其中包括主节点和子节点，但不包括定义的梁单元方向节点。

在普通管道裂纹扩展问题中，取真实模型对称部分的1/2进行计算。如径向取单节点（默认）；周向对模型做16等分，即11.25°；纵向沿管长250等分，则总节点的数目为：l×(16+l) ×(250+l)=4267。

一般用下式计算：NNODE=(NODGE1+l) ×(NODGE2+1）。

·Columns    6-10  Fortran Name        NPRI

模型中的主节点数目，默认等于总节点数NNODE。

◎Columns    11-15  Fortran Name        NAXOR

梁单元的方向节点数，用于确定梁单元的局部y轴方向，管道问题中默认为零。

·Columns    16-20    Fortran Name        NELE

模型中的单元数目。如径向取单节点（默认）；周向对模型做16等分，即11.25°；纵向沿管长250等分，则模型单元数目为：1×16×250=4000。

一般用下式计算：NNODE=NODGE1×NODGE2。

◎Columns    21-25    Fortran Name        NUMMAT

单元组（即一起定义或生成的单元群体）的数目，默认为1。与CARD4第一项对应。

·Columns    26-30    Fortran Name        NUMDIS

给定位移边界的节点数，由CARD7中参数统计确定。

#Columns    31-35    Fortran Name        MXSTEP

描述：计算时间步数。

计算中选取的时间步总数，由用户指定。

◎Columns    36-40    Fortran Name      NDGREE

每节点的自由度数目，默认为6。

#Columns    41-50    Fortran Name     DELT

描述：时间步长（单位：S）。

计算选取的时间步长，由用户指定。

·Columns    51-55    Fortran Name     NODMAX

模型中最大节点编号，管道问题中与本行第一项NNODE同。

#Columns    31-35    Fortran Name     NHISTORY（后添加）

描述：输出历史文件时间步数间隔。

和时间历史有关的变形，应力等输出文件的输出步数间隔，由用户指定。

CARD  3        注释行：CONTROL DATA

◎Columns    1-5       Fortran Name        KONTRL（l)

有无子节点（局部坐标）。为O表示没有；为l表示有。默认为0。

◎Columns        6-10      Fortran Name        KONTRL(2)

管道防护参数。为O表示没有止裂措施；为l或2表示有刚性套箍（1型套箍位置为同心，2型套箍位置为内切）；为3表示有柔性止裂环；为5表示埋置于泥土之中：为6表示各向异性材料（如纤维增强聚乙烯管道）。默认为0。

◎Colunlns     11-15      Fortran Name        KONTRL(3)

初始变形体轴线位置。为O表示以变形前的惯量主轴作为变形体轴线，为1表示轴线与整体坐标轴相一致。默认为1。

#Columns      16-20     Fortran Name          KONTRL(4)

描述：输出图线参数（选项形式）。

为2表示绘制全部节点，为l表示绘制沿轴线各节点，为O表示不绘制。

#Columns     21-25     Fortran Name          KONTRL(5)

描述：输出绘图数据到文件fort.15（选项形式）。

为1表示输出，为0表示不输出。

◎Columns     26-30     Fortran Name          KONTRL(6)

暂时无意义。默认为O。

◎Columns     31-35     Fortran Name          KONTRL(7)

默认节点类型。为1表示主节点，为0表示子节点。默认为1。

#Columns     36-40     Fortran Name          KONTRL(8)

描述：气体压力模式参数（选项形式）。

以此项表述裂纹后面气体压力衰减模型。为8表示采用线性衰减模式，衰减长度等于2倍（可调）管径；为7表示采用线性衰减模式，衰减长度等于1.5倍（可调）管径；为6表示采用阶跃衰减模式，以裂尖为分界线；为5表示采用Fourier级数衰减模式，衰减长度随裂纹扩展速度变化；为4表示采用Fourier级数衰减模式，衰减长度等于1.25倍直径；为l,2，或3分别表示沿x,y,z向采用抛物线衰减模式；为O或小于0表示不采用气体压力模式。默认为7。为9表示选择回流模式。

注意此项大于0时要添加CARD4*。

#Colunmns     41-45     Fortran Name          KONTRL(9)

描述：应力输出点个数。

每个单元内部应力输出点的个数：为O则不输出。

#Columns     46-50     Fortran Name          KONTRL(10)

描述：起算步输出设置。

在时间步数为N=KONTRL(10）的地方输出断点信息到文件fort.13，便于下次从这里开始计算。为O表示不设置输出点。

#Columns     51-55     Fortran Name          KONTRL(11)

描述：校核参数（选项形式）。

简化输出，以便数据校核。为-3表示仅输出模型状态及时间步分配情况待查：为-2表示略去所有过程输出结果；为-1表示仅输出基本模型数据待查；为0表示正常输出。

◎Columns     56-60     Fortran Name          KONTRL(12)

暂时无意义。默认为O。

#Columns     61-65     Fortran Name          KONTRL(13)

描述：输出变量存储参数（选项形式）。

为1表示将输出变量放入文件fort.14存储，为O表示用内存数组进行存储。

◎Columns     66-70     Fortran Name          KONTRL(14)

暂时无意义。默认为0。

◎Columns     71-75     Fortran Name          KONTRL(15)

为1表示裂纹持续扩展，并用节点力释放法计算裂纹驱动力；为O表示止裂。默认为1。

◎Columns     76一80     Fortran Name          KONTRL（16)

暂时无意义。默认为0。

Ⅱ.材料与单元属性

CARD 4*          注释行：PIPE DIAMETER

存在前提：(KONTRL（8)>0）注：不满足该条件时不将CARD4*写入输入文件，但以下输入项在前处理中照常。

Line l

#Columns          l-10          Fortran Name     PDIA

描述：管道直径（单位：m）。

一般建议此项取外径，由用户指定。

#Columns         N/A         Name         LENGTH

描述：管道全长（单位：m）。

由用户指定，在前处理中输入，不写入计算输入文件。

#Columns         11-20         Fortran Namne        VSS

描述：稳态裂纹扩展速度（单位：m/s)(KONTRL(17）>=1时表示裂纹扩展的初速度）。

用于气体压力模式的计算，由用户指定。

#Columns         21-25         Fortran Name      KONTRL(17)（后添加）

描述：速度衰减参数。

为2表示选择韧性减速模型；为1表示选择速度衰减模型：为0表示按恒定速度计算。

Line 2存在前提：KONTRL(17)＞=l（速度衰减模型或者韧性减速模型）

#Columns        l-10          Fortran Name        LR（后添加）

前提：KONTRL(17)=l

描述：速度衰减距离（单位：m）。

用于速度衰减模型的计算，由用户指定。

#Columns        11-20       Fortran Name        VRO（后添加）

前提：KONTRL(17)=l

描述：衰减速度参数（单位：m/s）。

用于速度衰减模型的计算，由用户指定。

#Columns        l-10        Fortran Name        TOU（后添加）

前提：KONTRL(17)=2

描述：管道的断裂韧性（单位：KJ/m）。

用于韧性减速模型的计算，由用户指定。

#Columns        11-20      Fortran Name        CRKBGT（后添加）

前提：KONTRL(17)=2

描述：起裂时间（单位：s）。

用于韧性减速模型的计算，由用户指定。

CARD  4            注释行：MATERIALANDELEMENT

Line l

◎Columns    l-5      Fortran Name      MTYP

单元组序号。默认为l，与NUMMAT对应。

◎Columns    6-10      Fortran Name        LTYP

单元类型。为1表示弹簧单元；为2表示弹性梁单元；为3表示弹塑性梁单元：为4表示弹性或弹塑性圆管单元；为5表示弹塑性平面曲梁单元（对应于某种止裂环结构）；为6表示三角形板单元；为7表示四边形板单元；为8表示六面体流体／固体单元。默认为7。

Line 2

#Columns        l-10       Fortran Name       E(l)

描述：材料密度（单位：103kg/m3）。

由用户指定。

#Columns        11-20      Fortran Name       E(2)

描述：杨氏模量（单位：KPa）。

由用户指定。

#Columns        21-30     Fortran Name        E(3)

描述：屈服应力（单位：KPa）。

由用户指定。对于弹性问题，应设为O。

#Columns        31-40     Fortran Name        E(4)

描述：泊松比。

由用户指定。

#Columns        41-50     Fortran Name       E(5）

描述：管壁厚度（单位：m）。

由用户指定。

上一页

下一页