L02 at125khz radialcrack=8le.txt

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

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

3、2=l1lcrack!缺陷在管道軸向的結(jié)束位置l=l1l3!管道長(zhǎng)度t1=deltnint(0.625e3delt)1delt!總的計(jì)算時(shí)間T1=7.620000000E04時(shí)間太短了!!a=2e6!激勵(lì)信號(hào)位移幅值ri=47e3!內(nèi)徑ro=51e3!外徑h=4e3!壁厚nv1=24!周向劃分的份數(shù)，也即一個(gè)圓周上節(jié)點(diǎn)的個(gè)數(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é)點(diǎn)坐標(biāo)系到柱坐標(biāo)系!csys0!施加軸對(duì)稱邊界條件nselslocy0nselalocy180dsymsymmy0allselallsavemodeleddb!!網(wǎng)格劃分結(jié)束!!!!!加載，求解!!!!!!!!SOLUANTYPE

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

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