L02 at125khz radialcrack=8le.txt

1、!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!003橫向裂紋radialcrack用軸對稱模型，只用建一半的模型!!!!!confignres15e4!!!用戶輸入開始!JPGPRF5001001!MACROTOSETPREFSFJPEGPLOTSFILNAMEL02at125khz8leradialcrack1!修改為nle.寬2mm長為nl_crackleprep7title004!圖片左下角顯示的圖標f=1

2、25e3!頻率nl_crack=25!裂紋軸向長度占幾個軸向單元長度，即寬nc_crack=2!裂紋周向長度占幾個周向單元，即長nr_crack=1!裂紋徑向長度占管道壁厚的比例，即高nsin=6!正弦激勵周期數(shù)delt=3e7!deltat時間步長!!le=2e3!單元長度軸向!!vmin=3238.06l1=450e3!缺陷在管道軸向的開始位置l3=550e3!從缺陷開始左端完整管道長度lcrack=nl_crackle!缺陷長度l

3、2=l1lcrack!缺陷在管道軸向的結(jié)束位置l=l1l3!管道長度t1=deltnint(0.625e3delt)1delt!總的計算時間T1=7.620000000E04時間太短了!!a=2e6!激勵信號位移幅值ri=47e3!內(nèi)徑ro=51e3!外徑h=4e3!壁厚nv1=24!周向劃分的份數(shù)，也即一個圓周上節(jié)點的個數(shù)。letheta=360nv1!管道周向劃分，單元角度數(shù)nr=1!徑向劃分份數(shù)!pc1=0!缺陷在管道周向的開始位

4、置!pc2=90!缺陷在管道周向的結(jié)束位置theta=nc_crackletheta!缺陷周向角度thick=nr_crackh!裂紋深度穿孔!et145mpex12.1e11!彈性模量mpnuxy10.29!泊松比mpdense17800!密度!!!!!!!!!!!!!!用戶輸入結(jié)束!!!!!!!!!!!!!!!!!lesizeallnrcsys1!裂紋處的上半部分管道徑向劃分lselsloczl1lselaloczl2lselrlo

5、cythetaletheta5thetaletheta5lesizeallnrvselall!選擇所有的體vatt111MSHAPE03Dallselallvsweepallcsys1NROTATall!轉(zhuǎn)換節(jié)點坐標系到柱坐標系!csys0!施加軸對稱邊界條件nselslocy0nselalocy180dsymsymmy0allselallsavemodeleddb!!網(wǎng)格劃分結(jié)束!!!!!加載，求解!!!!!!!!SOLUANTYPE

6、TRANSnlgeomoff!大變形選項eqslvpcg1e5!設(shè)置求解器為PCG收斂精度為1e5kbc1!設(shè)置載荷類型，為階躍載荷timet1!設(shè)置載荷終止時間deltimdelt!設(shè)置時間步長OUTPRbasic!Controlsthesolutionprintout.輸出所有的的基本求解結(jié)果OUTRESallall!Controlsthesolutiondatawrittentothedatabase.寫所有的結(jié)果到數(shù)據(jù)庫中nse

7、lsloczlle10lle10!約束另一端的所有自由度dallall0!DNODELabVALUEVALUE2NENDNINCLab2Lab3Lab4Lab5Lab6NSELSLOCZ0dimdisuytableNsub11!位移disuz(i)tst=delt!初始加載時間doi1Nsub1!分成了130段，則有131個數(shù)據(jù)t=tst(i1)delt!table的0行和0列要賦值這不同于數(shù)組array.disuy(i01)=t!di