基于VCCT的裂紋擴(kuò)展模擬

..基于VCCT的裂紋擴(kuò)展模擬虛擬裂紋閉合技術(shù)(VCCT)最初用于計(jì)算裂紋體的能量釋放率。因此被廣泛用于層合復(fù)合材料的界面裂紋擴(kuò)展模擬，并假定裂紋擴(kuò)展總是沿著預(yù)先定義的路徑，特別是在界面處?；赩CCT的裂紋擴(kuò)展模擬,當(dāng)前可用的線性單元如plane182和solid185?；赩CCT的裂紋擴(kuò)展模擬包含下述假定：·裂紋擴(kuò)展沿著預(yù)先定義的路徑·路徑通過(guò)界面單元來(lái)定義·分析為準(zhǔn)靜態(tài)分析，不考慮瞬態(tài)效應(yīng)·材料為線彈性材料，可以是各向同性，正交各向異性,各向性。裂紋可以位于一種材料或者兩種材料的界面。斷裂準(zhǔn)那么基于采用VCCT方法計(jì)算的能量釋放率?？刹捎枚喾N斷裂準(zhǔn)那么或自定義的準(zhǔn)那么。同一分析中可定義多條裂紋。VCCT裂紋擴(kuò)展模擬使用：·界面單元INTER202〔2D〕和INTER205〔3D〕·CINT命令計(jì)算能量釋放率·CGROW命令定義裂紋擴(kuò)展集，斷裂準(zhǔn)那么，裂紋擴(kuò)展路徑和求解控制參數(shù)。12.1.1VCCT裂紋擴(kuò)展模擬過(guò)程基于VCCT的裂紋擴(kuò)展模擬假定為準(zhǔn)靜態(tài)模擬。下面為進(jìn)展模擬的主要步驟：Step1：建立預(yù)先定義裂紋路徑的有限元模型Step2：進(jìn)展能量釋放率的計(jì)算Step3：進(jìn)展裂紋擴(kuò)展計(jì)算裂紋擴(kuò)展模擬為非線性構(gòu)造分析，這里詳述了一些特點(diǎn)，特別是裂紋擴(kuò)展的分析細(xì)節(jié)。12.1.1.1Step1：建立預(yù)先定義裂紋路徑的有限元模型標(biāo)準(zhǔn)的非線性求解過(guò)程需要建立有限元模型，有正確的求解控制設(shè)置，載荷和邊界條件。預(yù)先定義的裂紋路徑離散為界面單元，并分為一個(gè)單元組，如以下圖所示：圖12.1采用界面單元離散裂紋路徑界面單元可以通過(guò)CZMESH命令劃分或者能生成界面單元的第三方工具劃分。MPC約束單元選項(xiàng)(KEYOPT(2)=1)在裂紋擴(kuò)展前把潛在的裂紋面綁定在一起。當(dāng)滿(mǎn)足斷裂準(zhǔn)那么時(shí)，MPC約束隨后釋放，從而裂紋擴(kuò)展。在二維問(wèn)題中，裂紋尖端后的一個(gè)界面單元如果在一個(gè)指定的子步滿(mǎn)足斷裂準(zhǔn)那么那么可能開(kāi)。在三維問(wèn)題中，裂紋前沿后的所有界面單元如果滿(mǎn)足斷裂準(zhǔn)那么可能開(kāi)。裂紋尖端/前沿周?chē)膯卧叽缬绊懩芰酷尫怕实挠?jì)算精度。當(dāng)程序采用修正算法,可能不能產(chǎn)生準(zhǔn)確的結(jié)果。改為沿著裂紋擴(kuò)展路徑使用一樣單元尺寸的網(wǎng)格。12.1.1.2Step2：進(jìn)展能量釋放率計(jì)算基于VCCT的裂紋擴(kuò)展模擬，必須首先進(jìn)展能量釋放率的計(jì)算。計(jì)算能量釋放率，采用CINT,TYPE,VCCT命令，隨后使用CINT命令指定其它選項(xiàng)比方裂紋尖端節(jié)點(diǎn)組和裂紋面/邊的法向。VCCT計(jì)算采用下述假定：●當(dāng)裂紋增加一個(gè)小量時(shí)，釋放的應(yīng)變能等于裂紋閉合一樣的小量所需的能量●當(dāng)裂紋擴(kuò)展一個(gè)小量時(shí)，裂紋尖端〔前沿〕位置的裂尖場(chǎng)〔變形〕不變。當(dāng)裂紋擴(kuò)展接近邊界或者兩條裂紋彼此接近時(shí)假定不在適用。因此，使用VCCT計(jì)算要仔細(xì)檢查分析結(jié)果。12.1.1.3.Step3進(jìn)展裂紋擴(kuò)展計(jì)算裂紋擴(kuò)展計(jì)算在應(yīng)力計(jì)算之后，solution步進(jìn)展。為了進(jìn)展裂紋擴(kuò)展計(jì)算，必須先定義裂紋擴(kuò)展集合，然后指定裂紋路徑，斷裂準(zhǔn)那么，裂紋擴(kuò)展求解控制。求解命令CGROW定義裂紋擴(kuò)展計(jì)算所有必需的參數(shù)。進(jìn)展裂紋擴(kuò)展計(jì)算步驟如下：Step3a：初始裂紋擴(kuò)展集Step3b：指定裂紋路徑Step3c：指定裂紋計(jì)算的ID和斷裂準(zhǔn)那么Step3d：指定裂紋擴(kuò)展的求解控制12.1.1.3.1Step3a：初始裂紋擴(kuò)展集定義裂紋擴(kuò)展集，使用CGROW,NEW,n命令，其中n是裂紋擴(kuò)展集的編號(hào)12.1.1.3.1Step3b：指定裂紋路徑定義裂紋路徑，采用CGROW,PATH,cmname,其中cmname是界面單元組的名稱(chēng)。12.1.1.3.1Step3c：指定裂紋計(jì)算的ID和斷裂準(zhǔn)那么指定裂紋計(jì)算ID,通過(guò)CGROW,CID,n命令，其中n是采用VCCT計(jì)算能量釋放率的裂紋計(jì)算(CINT)的ID?！睠INT命令定義的參數(shù)和斷裂參數(shù)計(jì)算一致〕對(duì)于簡(jiǎn)單的斷裂準(zhǔn)那么，比方臨界能量釋放率，可以通過(guò)命令CGROW,FCOPTION,GC,指定，其中value為臨界能量釋放率。對(duì)于一些更復(fù)雜的斷裂準(zhǔn)那么，可以通過(guò)材料數(shù)據(jù)表定義斷裂準(zhǔn)那么。采用CGROW，F(xiàn)COPTION，MTAB，matid命令，其中matid是材料表的材料ID號(hào)。有多種斷裂準(zhǔn)那么可以使用，比方linear,bilinear,B-K,修正B-K，PowerLaw，和用戶(hù)自定義的準(zhǔn)那么。更多信息請(qǐng)參閱TB，CGCR命令和Fracturecriterial.對(duì)于每個(gè)裂紋擴(kuò)展集，只可以定義一條斷裂準(zhǔn)那么和一個(gè)單元組。可以采用不同的斷裂準(zhǔn)那么來(lái)定義多個(gè)裂紋擴(kuò)展集。多條裂紋可以同時(shí)擴(kuò)展或者彼此獨(dú)立。當(dāng)多條裂紋位于同一界面時(shí)，也可以合并為一條裂紋，如以下圖所示：圖12.2裂紋擴(kuò)展和合并也可以在各自的斷裂擴(kuò)展集中，對(duì)同一裂紋定義不同的斷裂準(zhǔn)那么。裂紋可以基于不同的準(zhǔn)那么〔根據(jù)哪條準(zhǔn)那么滿(mǎn)足〕擴(kuò)展，并且彼此獨(dú)立。這種方法對(duì)于比擬斷裂機(jī)理很有幫助。使用CGROW命令定義求解控制參數(shù)如下:指定求解控制參數(shù)使用CGROW命令:\o"12.1.3.

FractureCriteria"斷裂準(zhǔn)那么系數(shù)(fc)CGROW,FCRAT,value,Value為該系數(shù)初始時(shí)間步長(zhǎng)〔裂紋擴(kuò)展開(kāi)場(chǎng)時(shí)〕CGROW,DTIME,value,Value為初始時(shí)間步長(zhǎng)為防止過(guò)度預(yù)測(cè)承載能力,指定一個(gè)小的初始時(shí)間步長(zhǎng)隨后裂紋擴(kuò)展的最小時(shí)間步長(zhǎng)CGROW,DTMIN,value,value為最小時(shí)間步長(zhǎng)大小隨后裂紋擴(kuò)展的最大時(shí)間步長(zhǎng)CGROW,DTMAX,value,value為最大時(shí)間步長(zhǎng)大小裂紋擴(kuò)展前沿節(jié)點(diǎn)允許的最大裂紋擴(kuò)展量CGROW,STOP,CEMX,value,value為允許的最大的裂紋擴(kuò)展量裂紋擴(kuò)展模擬很耗時(shí),當(dāng)?shù)竭_(dá)感興趣的指定裂紋擴(kuò)展量時(shí),使用該命令來(lái)終止分析當(dāng)裂紋迅速擴(kuò)展時(shí)(比方,裂紋擴(kuò)展不穩(wěn)定),使用較小的DTMAX和DTMIN來(lái)容許載荷重新平衡。當(dāng)裂紋不再增長(zhǎng)時(shí)，指定的時(shí)間步長(zhǎng)控制被忽略，結(jié)果依賴(lài)于標(biāo)準(zhǔn)時(shí)間步長(zhǎng)控制。12.1.1.4.例如:裂紋擴(kuò)展集定義下面的輸入例如定義一個(gè)裂紋擴(kuò)展集：CGROW,NEW,1CGROW,CPATH,cpath1CGROW,FCOPTION,MTAB,5CGROW,DTIME,1.0e-4CGROW,DTMIN,1.0e-5CGROW,DTMAX,2.0e-3...裂紋擴(kuò)展在裂紋擴(kuò)展模擬中，一個(gè)關(guān)心的量是裂紋擴(kuò)展量。VCCT方法測(cè)量裂紋擴(kuò)展是基于已經(jīng)開(kāi)的界面單元的長(zhǎng)度，如下述方程和圖片所示：圖12.3二維和三維裂紋擴(kuò)展對(duì)于二維問(wèn)題，裂紋擴(kuò)展是當(dāng)前已經(jīng)開(kāi)的界面單元的長(zhǎng)度之和〔a〕。對(duì)于三維問(wèn)題，裂紋擴(kuò)展在每個(gè)裂紋前沿節(jié)點(diǎn)進(jìn)展測(cè)量，為沿著裂紋擴(kuò)展方向的界面單元邊長(zhǎng)的和〔b).裂紋擴(kuò)展量〔CEXT)是裂紋求解結(jié)果的一局部，和裂紋計(jì)算的ID號(hào)一致，可以和能量釋放率一樣通過(guò)POST1和POST26后處理命令〔比方PRCINT,PLCINT,和CISOL)進(jìn)展后處理。斷裂準(zhǔn)那么為建立裂紋擴(kuò)展,必須定義裂紋開(kāi)場(chǎng)和隨后裂紋擴(kuò)展的斷裂準(zhǔn)那么。對(duì)于線彈性斷裂力學(xué)，斷裂準(zhǔn)那么通常假定為三種斷裂模式的臨界能量釋放率的函數(shù)。表達(dá)為：對(duì)于一些模型可能需要其它的參數(shù)。當(dāng)斷裂準(zhǔn)那么滿(mǎn)足時(shí)，發(fā)生斷裂，表述為：其中QUOTE為斷裂準(zhǔn)那么比率。推薦值為0.95到1.05，默認(rèn)為1.0?？梢允褂靡韵聰嗔褱?zhǔn)那么：·臨界能量釋放率準(zhǔn)那么·線性斷裂準(zhǔn)那么·雙線性斷裂準(zhǔn)那么·B-K斷裂準(zhǔn)那么·修正B-K斷裂準(zhǔn)那么·指數(shù)斷裂準(zhǔn)那么·自定義斷裂準(zhǔn)那么用戶(hù)自定義選項(xiàng)需要提供子程序來(lái)定義你自己的斷裂準(zhǔn)那么。12.1.3.1.臨界能量釋放率準(zhǔn)那么臨界能量釋放率準(zhǔn)那么使用總的能量釋放率〔GT〕作為斷裂準(zhǔn)那么?？偟哪芰酷尫怕适侨N模式的能量釋放率的和，表述為：QUOTE其中QUOTE為臨界能量釋放率。對(duì)于I型斷裂模式，斷裂準(zhǔn)那么簡(jiǎn)化為：QUOTE能量釋放率準(zhǔn)那么是最簡(jiǎn)單的斷裂準(zhǔn)那么，適用于2D和3D的斷裂擴(kuò)展模擬。例如12.1臨界能量釋放率輸入gtcval=10.0CGROW,FCOPTION,GTC,gtcval12.1.3.2.線性斷裂準(zhǔn)那么線性選項(xiàng)假設(shè)斷裂準(zhǔn)那么是三種模式的能量釋放率的線性函數(shù)。表述為：其中QUOTE，QUOTE，QUOTE分別為I型，II型，III型斷裂模式的臨界能量釋放率。這三個(gè)值通過(guò)命令TBDATA輸入，如下：ConstantTBDATAInputmentsQUOTE

C1ModeI的臨界能量釋放率,

QUOTE>0

QUOTEC2ModeII的臨界能量釋放率,

QUOTE>0C3ModeIII的臨界能量釋放率,

QUOTE>0例如12.2線性準(zhǔn)那么輸入g1c=10.0g2c=20.0g3c=25.0TB,CGCR,1,,,LINEARTBDATA,1,g1c,g2c,g3c三個(gè)參數(shù)不能同時(shí)為零。如果其中一個(gè)設(shè)為零，相應(yīng)的項(xiàng)被忽略。當(dāng)三個(gè)臨界能量釋放率相等，那么線性斷裂準(zhǔn)那么簡(jiǎn)化為臨界能量釋放率準(zhǔn)那么。線性斷裂準(zhǔn)那么適用于當(dāng)三種斷裂模式的臨界能量釋放率明顯存在時(shí)的三維混合斷裂模式。12.1.3.3.雙線性斷裂準(zhǔn)那么雙線性斷裂選項(xiàng)假設(shè)斷裂準(zhǔn)那么是I型和II型斷裂模式的能量釋放率的線性函數(shù)，表達(dá)式為：其中QUOTE，QUOTE分別為I型，II型斷裂模式的臨界能量釋放率，ξ和ζ為材料常數(shù)。四個(gè)值都能通過(guò)TBDATA命令定義為溫度的函數(shù)。如下所示：Example12.3雙向性準(zhǔn)那么輸入g1c=10.0g2c=20.0x=2y=2TB,CGCR,1,,,BILINEARTBDATA,1,g1c,g2c,x,y雙線性斷裂準(zhǔn)那么適用于二維混合斷裂模式的模擬。12.1.3.4.B-K斷裂準(zhǔn)那么B-K選項(xiàng)表述為：其中QUOTE，QUOTE，分別為I型，II型斷裂模式的臨界能量釋放率,η為材料常數(shù)。三個(gè)值都能夠通過(guò)TBDATA命令定義為溫度的函數(shù)。如下所示：B-K準(zhǔn)那么本來(lái)用于復(fù)合材料界面斷裂，適用于三維混合斷裂模式的模擬。Example12.4B-K準(zhǔn)那么輸入g1c=10.0g2c=20.0h=2TB,CGCR,1,,,BKTBDATA,1,g1c,g2c,h12.1.3.5.修正B-K斷裂準(zhǔn)那么修正B-K選項(xiàng),表述為:其中QUOTE，QUOTE，QUOTE分別為I型，II型，III型斷裂模式的臨界能量釋放率,η為材料常數(shù)。四個(gè)值都能夠通過(guò)TBDATA命令定義為溫度的函數(shù)。如下所示：QUOTE=QUOTE時(shí)，修正B-K準(zhǔn)那么簡(jiǎn)化為B-K準(zhǔn)那么。修正B-K準(zhǔn)那么用于復(fù)合材料界面斷裂，考慮明顯的II型和III型臨界能量釋放率，適用于三維混合斷裂模式的模擬。Example12.5修正B-K準(zhǔn)那么輸入g1c=10.0g2c=20.0g3c=25.0h=2TB,CGCR,1,,,MBKTBDATA,1,g1c,g2c,g3c,h12.1.3.6.指數(shù)斷裂準(zhǔn)那么指數(shù)準(zhǔn)那么選項(xiàng)假設(shè)斷裂準(zhǔn)那么是三種模式的能量釋放率的指數(shù)函數(shù)。表述為：其中QUOTE，QUOTE，QUOTE分別為I型，II型，III型斷裂模式的臨界能量釋放率。QUOTE,QUOTE,QUOTE是指數(shù)，為常數(shù)。六個(gè)量都可以通過(guò)TBDATA命令定義為溫度的函數(shù)，如下所示:三個(gè)臨界能量釋放率不能同時(shí)為零。如果其中一個(gè)設(shè)為零，相應(yīng)的項(xiàng)被忽略。當(dāng)指數(shù)QUOTE,QUOTE,QUOTE設(shè)置為1時(shí)，指數(shù)準(zhǔn)那么簡(jiǎn)化為線性斷裂準(zhǔn)那么。指數(shù)準(zhǔn)那么適用于當(dāng)三種斷裂模式的臨界能量釋放率明顯存在時(shí)的三維混合斷裂模式。例如12.6指數(shù)準(zhǔn)那么輸入g1c=10.0g2c=20.0g3c=25.0n1=2n2=2n3=3TB,CGCR,1,,,POWERLAWTBDATA,1,g1c,g2c,g3c,n1,n2,n312.1.3.7.用戶(hù)自定義斷裂準(zhǔn)那么一個(gè)自定義的斷裂準(zhǔn)那么可以表述為：其中斷裂準(zhǔn)那么是三種斷裂模式的能量釋放率和材料常數(shù)的函數(shù)。所有的值通過(guò)TBDATA命令輸入。你必須提供一個(gè)子程序。下面為一個(gè)子程序定義線性斷裂準(zhǔn)那么的例子：*deck,user_cgfcritoptimizeSUBROUTINEuser_cgfcrit(cgi,cid,kct,&nprop,prop,fcscl,&var1,var2,var3,var4)c*****************************************************************cc***primaryfunction:cputefacturecriterionforcrackgrowthcuserfracturecriterionexamplec***notice:thisroutinecontainssasiconfidentialinformation***c#include"imp.inc"#include"ansysdef.inc"ccinputargumentsc===============ccgi(int,sc,in)CGROWsetidccid(int,sc,in)CINTIDtobeusedckct(int,sc,in)Currentcracktipnodecnprop(int,sc,in)numberofpropertiescprop(dp,ar(*),in)propertyarrayccOutputargumentsc===============cfcscl(dp,sc,ou)fracturecriterioncareturnvalueofoneorbiggercindicatesfractureccMisc.argumentsc===============cvar1(,,)notusedcvar2(,,)notusedcvar3(,,)notusedcvar4(,,)notusedcc*****************************************************************cc***subroutines/functionc***get_cgfpar:APItoaccessfracturedatac***wrinqr:ansysstandardiofunctionc***externalget_cgfparexternalwrinqrintegerwrinqrc***argumentcINTEGERcgi,cid,kct,npropdoubleprecisionfcscl,&var1,var2,var3,var4doubleprecisionprop(nprop)cc***localvariablecintegerdebugflag,iottintegernndoubleprecisiong1c,g2c,g3c,g1,g2,g3doubleprecisiongs(4),da(1)cc***localparametersDOUBLEPRECISIONZERO,ONEparameter(ZERO=0.0d0,&ONE=1.0d0)cc*****************************************************************c***initializationfcscl=ZEROc***retrieveenergy-releaseratesc***forcrackcidandcracktipnodekctc***gs(1:3)willbereturnedasG1,G2,G3c***getenergy-releaseratesnn=3gs(1:nn)=ZEROcallget_cgfpar('GS',cid,kct,0,nn,gs(1))c***getcrackextensionnn=1da(1)=ZEROcallget_cgfpar('DA',cid,kct,0,nn,da(1))c***energy-releaseratesg1=abs(gs(1))g2=abs(gs(2))g3=abs(gs(3))c***inputpropertyfromTBDATA,1,c1,c2,c3g1c=prop(1)g2c=prop(2)g3c=prop(3)c***linearfracturecriterionfcscl=ZEROif(g1c.gt.TINY)fcscl=fcscl+g1/g1cif(g2c.gt.TINY)fcscl=fcscl+g2/g2cif(g3c.gt.TINY)fcscl=fcscl+g3/g3cc***userdebugoutputdebugflag=1if(debugflag.gt.0)theniott=wrinqr(WR_OUTPUT)write(iott,1000)cgi,cid,kct,da(1),fcscl,gs(1:3)1000format(5x,'userfracturecriterion:'/&5x,'crackgrowthsetID=',i5/&5x,'crackID=',i5/&5x,'cracktipnode=',i5/&5x,'crackextension=',g11.5/&5x,'calculatedfractureparameter=',g11.5/&5x,'energy-releaseratesGs(1:3)=',3g12.5)endifreturnend裂紋擴(kuò)展模擬例題該例子采用雙懸臂梁,一端有線裂紋。在梁端裂紋上下施加相反方向的一樣位移，來(lái)開(kāi)裂紋，如下圖：圖12.4雙懸臂梁的裂紋擴(kuò)展以下圖顯示了有限元網(wǎng)格：采用PLANE182單元，翻開(kāi)enhancedstrain選項(xiàng)〔keyopt〔1〕=2〕來(lái)建立實(shí)體模型。INTER202單元來(lái)建立裂紋路徑。假設(shè)為平面應(yīng)變條件。在垂直方向采用6個(gè)單元，水平方向200個(gè)單元。以下圖顯示了預(yù)測(cè)的荷載-位移曲線。圖12.6雙懸臂梁荷載-位移曲線在裂紋開(kāi)場(chǎng)擴(kuò)展前，反力隨著位移增大并迅速到達(dá)峰值。然后在裂紋擴(kuò)展的初始階段，反力迅速下降，隨著裂紋增長(zhǎng)下降速度變慢。結(jié)果和文獻(xiàn)結(jié)果符合的很好。下面給出最大主應(yīng)力的云圖。圖12.7雙懸臂梁的云圖下面給出該雙懸臂梁的裂紋擴(kuò)展模擬的輸入文件：/BATCH/TITLE,,CRACKGROWTHSIMULATIONOFADCBPROBLEM-2DPLANESTRAIN/PREP7DIS1=0.9DIS2=12.0N1=1000N2=1000N3=10DL=100DH=3A0=30NEL=200NEH=6TOLER=0.1E-5ET,1,182!*2D4-NODESTRUCTURALSOLIDELEMENTKEYOPT,1,1,2!*ENHANCESTRAINFORMULATIONKEYOPT,1,3,2!*PLANESTRAINET,2,182KEYOPT,2,1,2KEYOPT,2,3,2ET,3,202!*2D4-NODECOHESIVEZONEELEMENT!KEYOPT,3,2,2!*ELEMENTFREEOPTIONKEYOPT,3,3,2!*PLANESTRAINMP,EX,1,1.353E5!*E11=135.3GPAMP,EY,1,9.0E3!*E22=9.0GPAMP,EZ,1,9.0E3!*E33=9.0GPAMP,GXY,1,5.2E3!*G12=5.2GPAMP,PRXY,1,0.24MP,PRXZ,1,0.24MP,PRYZ,1,0.46G1C=0.28!*CRITICALENERGY-RELEASERATEG2C=0.80G3C=0.80TB,CGCR,1,,3,LINEAR!*LINEARFRACTURECRITERIONTBDATA,1,G1C,G2C,G3C!FEMODELRECTNG,0,DL,DH/2!*DEFINEAREASRECTNG,0,DL,0,-DH/2LSEL,S,LINE,,2,8,2!*DEFINELINEDIVISIONLESIZE,ALL,DH/NEHLSEL,INVELESIZE,ALL,,,NELALLSEL,ALLTYPE,1!*MESHAREA2MAT,1LOCAL,11,0,0,0,0ESYS,11AMESH,2CSYS,0TYPE,2!*MESHAREA1ESYS,11AMESH,1CSYS,0NSEL,S,LOC,X,A0-TOLER,DLNUMMRG,NODESESLNTYPE,3MAT,5CZMESH,,,1,Y,0,!*GENERATEINTERFACEELEMENTSALLSEL,ALLNSEL,S,LOC,X,DL!*APPLYCONSTRAINTSD,ALL,ALLNSEL,ALL!ESEL,S,ENAME,,202!*SELECTINTERFACEELEMENTTOCM,CPATH,ELEM!*DEFINECRACKGROWTHPATHNSLENLISTNSEL,S,LOC,X,A0NSEL,R,LOC,Y,0NLISTESLNELISTCM,CRACK1,NODE!*DEFINECRACKTIPNODEPONENTNLISTALLSFINISH/SOLURESC,,NONEESEL,S,TYPE,,2NSLE,SNSEL,R,LOC,XNSEL,R,LOC,Y,DH/2!*APPLYDISPLACEMENTLOADINGONTOPD,ALL,UY,DIS1NSEL,ALLESEL,ALLESEL,S,TYPE,,1NSLE,SNSEL,R,LOC,XNSEL,R,LOC,Y,-DH/2!*APPLYDISPLACEMENTLOADINGONBOTTOMD,ALL,UY,