附錄  程序中的變量說明

輸入文件名：curve.dat，單位采用10³kg、m、s，去掉注釋行以后存為fort.16。

注釋行前四個字符為“C***”或“C”，終結行前四個字符為“STOP”。

以下標◎的欄目未做到前處理界面中，采用默認值的處理辦法。

以下標·的欄目為通過前處理界面的給定值計算得到，不必另行輸入。

以下標＃的欄目為前處理界面中要求輸入的內(nèi)容，同時附有“描述”項。

Ⅰ.計算控制參數(shù)

CARD   1        無注釋行

#Columns    1-80    Fortran Name      TITLE

描述：項目主題。

長度限制為80字符，作為對該計算問題的文字描述，并在輸出中作為標識。

如該行內(nèi)容為“RESTART”，則程序讀取fort.13中存取的斷點信息繼續(xù)上次未完成的計算。

CARD   2      注釋行：GENERALDATA

·Columns    l-5    Fortran Name      NNODE

模型中的總節(jié)點數(shù)，其中包括主節(jié)點和子節(jié)點，但不包括定義的梁單元方向節(jié)點。

在普通管道裂紋擴展問題中，取真實模型對稱部分的1/2進行計算。如徑向取單節(jié)點（默認）；周向?qū)δＰ妥?6等分，即11.25°；縱向沿管長250等分，則總節(jié)點的數(shù)目為：l×(16+l) ×(250+l)=4267。

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

·Columns    6-10  Fortran Name        NPRI

模型中的主節(jié)點數(shù)目，默認等于總節(jié)點數(shù)NNODE。

◎Columns    11-15  Fortran Name        NAXOR

梁單元的方向節(jié)點數(shù)，用于確定梁單元的局部y軸方向，管道問題中默認為零。

·Columns    16-20    Fortran Name        NELE

模型中的單元數(shù)目。如徑向取單節(jié)點（默認）；周向?qū)δＰ妥?6等分，即11.25°；縱向沿管長250等分，則模型單元數(shù)目為：1×16×250=4000。

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

◎Columns    21-25    Fortran Name        NUMMAT

單元組（即一起定義或生成的單元群體）的數(shù)目，默認為1。與CARD4第一項對應。

·Columns    26-30    Fortran Name        NUMDIS

給定位移邊界的節(jié)點數(shù)，由CARD7中參數(shù)統(tǒng)計確定。

#Columns    31-35    Fortran Name        MXSTEP

描述：計算時間步數(shù)。

計算中選取的時間步總數(shù)，由用戶指定。

◎Columns    36-40    Fortran Name      NDGREE

每節(jié)點的自由度數(shù)目，默認為6。

#Columns    41-50    Fortran Name     DELT

描述：時間步長（單位：S）。

計算選取的時間步長，由用戶指定。

·Columns    51-55    Fortran Name     NODMAX

模型中最大節(jié)點編號，管道問題中與本行第一項NNODE同。

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

描述：輸出歷史文件時間步數(shù)間隔。

和時間歷史有關的變形，應力等輸出文件的輸出步數(shù)間隔，由用戶指定。

CARD  3        注釋行：CONTROL DATA

◎Columns    1-5       Fortran Name        KONTRL（l)

有無子節(jié)點（局部坐標）。為O表示沒有；為l表示有。默認為0。

◎Columns        6-10      Fortran Name        KONTRL(2)

管道防護參數(shù)。為O表示沒有止裂措施；為l或2表示有剛性套箍（1型套箍位置為同心，2型套箍位置為內(nèi)切）；為3表示有柔性止裂環(huán)；為5表示埋置于泥土之中：為6表示各向異性材料（如纖維增強聚乙烯管道）。默認為0。

◎Colunlns     11-15      Fortran Name        KONTRL(3)

初始變形體軸線位置。為O表示以變形前的慣量主軸作為變形體軸線，為1表示軸線與整體坐標軸相一致。默認為1。

#Columns      16-20     Fortran Name          KONTRL(4)

描述：輸出圖線參數(shù)（選項形式）。

為2表示繪制全部節(jié)點，為l表示繪制沿軸線各節(jié)點，為O表示不繪制。

#Columns     21-25     Fortran Name          KONTRL(5)

描述：輸出繪圖數(shù)據(jù)到文件fort.15（選項形式）。

為1表示輸出，為0表示不輸出。

◎Columns     26-30     Fortran Name          KONTRL(6)

暫時無意義。默認為O。

◎Columns     31-35     Fortran Name          KONTRL(7)

默認節(jié)點類型。為1表示主節(jié)點，為0表示子節(jié)點。默認為1。

#Columns     36-40     Fortran Name          KONTRL(8)

描述：氣體壓力模式參數(shù)（選項形式）。

以此項表述裂紋后面氣體壓力衰減模型。為8表示采用線性衰減模式，衰減長度等于2倍（可調(diào)）管徑；為7表示采用線性衰減模式，衰減長度等于1.5倍（可調(diào)）管徑；為6表示采用階躍衰減模式，以裂尖為分界線；為5表示采用Fourier級數(shù)衰減模式，衰減長度隨裂紋擴展速度變化；為4表示采用Fourier級數(shù)衰減模式，衰減長度等于1.25倍直徑；為l,2，或3分別表示沿x,y,z向采用拋物線衰減模式；為O或小于0表示不采用氣體壓力模式。默認為7。為9表示選擇回流模式。

注意此項大于0時要添加CARD4*。

#Colunmns     41-45     Fortran Name          KONTRL(9)

描述：應力輸出點個數(shù)。

每個單元內(nèi)部應力輸出點的個數(shù)：為O則不輸出。

#Columns     46-50     Fortran Name          KONTRL(10)

描述：起算步輸出設置。

在時間步數(shù)為N=KONTRL(10）的地方輸出斷點信息到文件fort.13，便于下次從這里開始計算。為O表示不設置輸出點。

#Columns     51-55     Fortran Name          KONTRL(11)

描述：校核參數(shù)（選項形式）。

簡化輸出，以便數(shù)據(jù)校核。為-3表示僅輸出模型狀態(tài)及時間步分配情況待查：為-2表示略去所有過程輸出結果；為-1表示僅輸出基本模型數(shù)據(jù)待查；為0表示正常輸出。

◎Columns     56-60     Fortran Name          KONTRL(12)

暫時無意義。默認為O。

#Columns     61-65     Fortran Name          KONTRL(13)

描述：輸出變量存儲參數(shù)（選項形式）。

為1表示將輸出變量放入文件fort.14存儲，為O表示用內(nèi)存數(shù)組進行存儲。

◎Columns     66-70     Fortran Name          KONTRL(14)

暫時無意義。默認為0。

◎Columns     71-75     Fortran Name          KONTRL(15)

為1表示裂紋持續(xù)擴展，并用節(jié)點力釋放法計算裂紋驅(qū)動力；為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

描述：穩(wěn)態(tài)裂紋擴展速度（單位：m/s)(KONTRL(17）>=1時表示裂紋擴展的初速度）。

用于氣體壓力模式的計算，由用戶指定。

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

描述：速度衰減參數(shù)。

為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

描述：衰減速度參數(shù)（單位：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表示彈塑性平面曲梁單元（對應于某種止裂環(huán)結構）；為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）。

由用戶指定。

上一頁

下一頁