碰撞模擬中的控制卡片設置及意義

1、1指定輸出文件Control Cards for ASCI I output(BNDOUTJNTFORC TSU 祚0.100 RCFORCRWFORC| 0.1000 100 I TPRIMTTRHIST$tDATABASE_OPTION -ABSTAT AVSFLTBinary Options | ASCIIBINAR¥MOVIE S0TOUT(SPHOLTTDEFGEO DEFORC fELOUT GCEOUT (GLSTAT' I 0.101 -O-tTOMPGS NCFORC NODFOR (tNODOUT RBDOUTSPCFORC(SSSTAT SWFORCre

2、jectdefaultABSTAT】一一氣囊統(tǒng)計表。輸出體積、壓強、內能、氣體質量流入率、氣體質量流出率、質量、溫度、密度?！続VSFLT】 AVS數(shù)據(jù)【BNDOUT】一一邊界環(huán)境的力和能量。輸出三個方向的力。DEFGEO】變形的幾何體的文件DEFORC】一一離散單元。輸出三個方向的力?！綞LOUT 】 單元數(shù)據(jù)。(見 DATABASE HISTORY OPTION )梁單元平面應力塊平面應變軸向合力xx,yy,zz 應力xx,yy,zz 應力xx,yy,zz 應變S方向剪切合力xy,yz,zx 應力xy,yz,zx 應力xy,yz,zx 應變T方向剪切合力塑性應變有效應力下表面應變S方向合力

3、矩屈服函數(shù)上表面應變T方向合力矩扭力合力GCEOUT】一一幾何接觸實體。包含三個方向力和力矩。GLSTAT】一一總體數(shù)據(jù)。 【JNTFORC】一一運動副力文件 NATSUM】材料能量。GLSTATJNTFORCMATSUM動能x,y,z三方向的力動能內能x,y,z三方向的力矩內能總能量沙漏能比率x,y,z三方向的動量剛性墻能量px,y,z三方向的剛體速度彈簧和阻尼能量總動能沙漏能總內能阻尼能總沙漏能滑移面能量外功x,y,z二方向速度時間步單元ID號控制的時間步MOVIE】【MPGS】【NCFORC】一一接觸面節(jié)點力【NODFOR】一一節(jié)點力組【NODOUT】一一節(jié)點數(shù)據(jù)NCFORCNODOUT

4、NODFORx方向力位移x,y,z二方向力y方向力速度z方向力加速度轉動量角速度角加速度【RBDOUT】一一剛體數(shù)據(jù) 【RCFORC】一一接觸面合成力【RWFORC】一一剛性墻所受的力RBDOUTRCFORCRWFORC二方向合位移二方向合力法向力二方向合速度:二方向合力二方向合加速度【SBTOUT】一一安全帶輸出文件【SECFORC 】一一橫截面通過的力(見 DATABASE_CROSS_SECTION_OPTION ) 【SLEOUT】一一滑移面的能量SPCFORC】一一單點約束的反作用力【SPHOUT 】 SPH 數(shù)據(jù)(見 DATABASE_HISTORY_OPTION )【SSSTAT

5、】一一子系統(tǒng)數(shù)據(jù)【SWFORC】一一節(jié)點約束反力(焊點和鉚釘)SECFORCSLEOUTSPCFORCSWFORCx,y,z二方向力Slave能量x,y,z二方向力軸向力x,y,z二方向力矩Master能量x,y,z二方向力矩剪切力x,y,z二方向中心面積合力【TPRINT】一一結構對的熱量輸出【TRHIST】一一追蹤質點時間歷程信息2時間步控制"CONTROL_TI MESTERDTINITTZ3rACISDOTSUMIT DT2MSLCTM | ERODE MS1ST10 000 |0 900o|Oo|-7.000e-00 0o【DTNIT】一一初始時間步長。EQ.0.0 :由L

6、S-DYNA決定時間步長?！綯SSFAC】一一計算時間步長的縮放因子。默認為0.9。【ISDO】一一4節(jié)點殼單元時間步長計算的根據(jù)。EQ.0 :特征長度=面積/min最長邊，最長對角線EQ.1 :特征長度=面積/最長對角線EQ.2:時間步長取決于 條波速度（bar wave speed和MAX最短邊，面積/min （最長 邊，最長對角線）。該選項提供的時間步長相對很大，可能導致計算的不穩(wěn)定，尤其是在 應用三角形單元時。EQ.3 :時間步長取決于最大特征值。該選項適用于材料的聲音傳播速度漸變的結構。 用于計算最大特征值的計算開銷是很有意義的，但時間步長的增長通常考慮不用質量縮放的較短的計算周期。

7、【TSUMIT】一一指定殼單元最小時間步長。當某一單元的時間步長小于給定值時，該單元 的材料屬性（彈性模量而不是質量）將被調整，使其時間步長不低于給定值。該選項只適用 于以下材料：MAT_PLASTIC_KINEMA TIC， MAT_POWER_LAW_PLASTICITY ， MAT_STRAIN_RATE_DEPENDENT_PLASTICITY，MAT_PIECEWISE_LINEAR_PLASTICITY。不推薦所謂的剛度縮放選項。下面的 DT2MS選項適用于所有材料和所有單元類型， 并且是首選的。如果TSUMIT和DT2MS兩個選項都 被激活并且TSUMIT值為正，則TSUMIT的

8、值自動置為1E-18，使其功能被屏蔽。如果其 值為負并且其絕對值大于IDT2MS 則| TSUMIT |優(yōu)先應用到質量縮放中，如果其絕對值小于|DT2MS |,貝U TSUMIT的值自動置為 1E-18?！綝T2MS】 控制質量縮放的時間步長。DT2MS 值應用范圍大于0用于初始影響無關緊要的準靜態(tài)分析和時間歷程分析等于0默認小于0允許的最小時間步長為TSSFAC* IDT2MS | ,當且僅當時間步長小于判斷標準時，質量縮放才會進仃。該選項可用于質量增加影響不大的瞬態(tài)分析。警告：超單元和 ELMENT DIRECT MA TRIX INPUT 不進行質量縮放，所以DT2MS不影響他們的時間步

9、長。這種情況下計算會出錯終止，DT2MS應輸入一個較小的值。【LCTM】一一限制最大時間步長的曲線。【ERODE】到達 TSMIN （見下面卡片 CONTROL_TERMINATION ）時，實體單元和t- 殼單元的侵蝕標記。如果此項不設，計算會終止。EQ.0:無侵蝕EQ.1:有侵蝕【MS1ST】一一限制第一步的質量縮放并且根據(jù)之前的時間步確定質量矢量。EQ.0 :否EQ.1 :是3計算終止時間設置* CONTROLTERMI NAT I ONENDTIMENDCYCIDTMIN ENDENGENDMASS| 100.000O.ODO|o.oo|0.00【ENDTIM】 強制的計算結束時間。E

10、NDTIM指定的計算結【ENDCYC】一一計算循環(huán)次數(shù)。當達到指定循環(huán)次數(shù)而沒有到達 束時間時，同樣終止計算。循環(huán)次數(shù)等于時間步的數(shù)目。【DTMIN】 初始時間步長的縮放因子，用以決定最小時間步長（TSMIN），TSMIN=DTSTART*DTMIN 。式中 DTSTART是由LS-DYNA 決定的初始時間步長。當達到 TSMIN時，LS-DYNA計算終止并輸出一個重啟動文件?！綞NDENG】一一控制計算結束的能量比例變化?！綞NDMASS】一一控制計算結束的質量比例變化。當且僅當用質量縮放控制最小時間步長時，該選項才起作用。4.單元控制* CONTROL.SHELLWRPANGESORTIR

11、NX< ISTUPD THEORY BWCMITERPROJ|25.QOOJ 12? 1 1o【W(wǎng)RPANG】一一殼單元翹曲角度。當某個翹曲角度大于給定值時，會輸出警告信息。默 認值為20?！綞SORT】一一自動挑選退化的四邊形單元，并處理為CO三角形單元公式，以保證求解穩(wěn)EQ.0:不挑選（默認）EQ.1 :完全挑選并處理【IRNXX】 單元法線更新選項。該選項影響Hughes-Liu， Belytschko-Wong-Chiang ，和Belytschko-Tsay單元公式。當且僅當翹曲剛度選項被激活時， 即BCW=1時，以上單元公 式才受影響。對于Hughes-Liu殼單元類型1，6

12、和7，IRNXX必須設為-2以調用上表面或下 表面作為參考面。EQ.-2 :EQ.-1 :每個循環(huán)都重新計算法線方向。EQ.0 :自動設為-1。EQ.1 :重啟動時計算。EQ.n :每n個循環(huán)重新計算法線方向。（只適用于Hughes-Liu殼單元類型）【ISTUPD】一一單元厚度變化選項。該選項影響所有單元公式EQ.0 :不變化。EQ.1 :膜變形引起厚度改變。該選項對金屬板料成型或所有膜片拉伸作用很大的情況 都很重要?！綯HEORY】缺省的殼單元理論。默認值為2, Belytschko-Tsay單元公式?！綛WC】Belytschko-Tsay單元公式的翹曲剛度選項。EQ.1 :增加 Bel

13、ytschko-Wong-Chiang 公式的翹曲剛度。EQ.2: Belytschko-Tsay單元公式。不增加翹曲剛度。（默認）【MITER】一一平面應力塑性選項。（適用于3, 18, 19和24號材料）EQ.1: 3 次交叉迭代（iterative plasticity with 3 iterations ）（默認） EQ.2 :完全迭代。EQ.3:不迭代?？赡軐е洛e誤，慎用?！綪ROJ】5沙漏控制*C0N7R0L_H0URGLASSIFQQH百C.020【IHQ】一一沙漏控制類型：EQ.1 :標準LS-DYNA類型。（默認）EQ.2: Flanagan-Belyschko 積分類型。E

14、Q.3:有精確體積的 Flanagan-Belyschko積分類型。EQ.4 :類型2的剛度形式。EQ.5 :類型3的剛度形式。EQ.6:EQ.8:適用于單元類型為16的全積分殼單元。當IHQ=8時，激活翹曲剛度，以得到精 確解。該選項會增加 25%的計算開銷。在殼單元中，IHQ<4的是基于Belyschko-Tsay公式的粘性沙漏控制模式，IHQ=4 , 5, 6 為剛度控制模式。剛度控制模式在大變形問題中可能使響應變得過于剛硬，使用時要注意。在高速問題中推薦采用粘性模式，在低速問題中推薦采用剛度模式。對于大變形問題，推薦使用選項3或5。0.15可能引起計算不穩(wěn)定。缺省值為0.1，可適

15、用于【QH】一一沙漏系數(shù)。該值如果超過 除IHQ=6以外的所有選項。備注：對個別組件的沙漏控制，可通過先建立沙漏屬性集合器，再從 組件集合器中調用沙漏屬性的方法實現(xiàn)。6.能量耗散控制"CONTROL.ENERGYHGEN RWEN SLNTEN RYUEN2 | 2 2 | HGEN】一一沙漏能計算選項。該選項需要大量存儲空間，并增加10%的計算開銷。EQ.1 :不計算沙漏能。（默認）EQ.2 :計算沙漏能并包含在能量平衡中，計算結果寫入GLSTAT和MATSUM文件中?！綬WEN】一一剛性墻能量耗散選項。EQ.1 :不計算剛性墻能量耗散。EQ.2 :計算剛性墻能量耗散并包含在能量平

16、衡中，計算結果寫入GLSTAT文件中。（默認）SLNTEN】一一滑移面能量耗散控制選項。EQ.1 :不計算滑移面能量耗散。（當接觸激活時，該選項自動設為2）GLSTAT 和 SLEOUTEQ.2 :計算滑移面能量耗散并包含在能量平衡中，計算結果寫入 文件中?！綬YLEN】一一阻尼衰減能量耗散控制選項EQ.1 :不計算阻尼衰減能量耗散。（默認）EQ.2 :計算阻尼衰減能量耗散并包含在能量平衡中，計算結果寫入GLSTAT文件中。7.CPU時間控制尢8忖 ITOL.CHJCPUTIM | o.o do"【CPUTIM】一一用于電流相位分析或重啟動。EQ.0 :沒有CPU時間限制。8接觸控制

17、SLSFAC】滑移面懲罰因子。EQ.0:缺省值=0.1【RWPNAL】一一剛性墻懲罰因子，用來處理剛體與固定剛性墻的相互作用。EQ.0.0 :不考慮剛體與剛性墻的相互作用。GT.0.0 :剛體與剛性墻進行相互作用，推薦值為1。每個從節(jié)點將保存 7個變量。如果所有節(jié)點對于剛性墻都是從節(jié)點，則會大大增加對內存的要求。【ISLCHK】一一接觸面初始穿透檢查。EQ.0 :自動設為1EQ.1:不檢查。EQ.2 :完全檢查初始穿透?！維HLTHK】一一在面-面接觸和點-面接觸類型中考慮殼單元厚度的選項。選項1和2會激活新的接觸算法。厚度偏置通常包括在單面接觸、約束算法、自動面面接觸和自動點面接觸類型中。E

18、Q.0 :不考慮厚度偏置。EQ.1 :考慮厚度偏置但剛體除外。EQ.2 :考慮厚度偏置，包括剛體?！綪ENOPT】一一懲罰剛度值選項。EQ.0 :自動設為1EQ.1 :主動截面和從節(jié)點剛度的最小值。（默認）EQ.2 :用主動截面的剛度值。（過去的方法）EQ.3:用從節(jié)點的剛度值。EQ.4:用從節(jié)點的剛度值，面積或質量加權。EQ.5 :與4相同，但是厚度加權。通常不推薦使用。 選項4和5推薦在金屬成型計算中使用。【THKCHG】一一單面接觸中考慮殼單元厚度變化的選項。EQ.0 :不考慮。（默認）EQ.1 :考慮殼單元厚度變化?！綩RIEN】一一初始化過程中接觸面截面自動再定位選項。EQ.0 :自

19、動設為1。EQ.1 :僅自動（part）輸入時激活。接觸面由part定義。EQ.2 :手動（segment）和自動輸入（part）都激活。EQ.3 :不激活?！綞NMASS】一一對接觸過程中銷蝕掉的節(jié)點的質量的處理。該選項影響所有當周圍單元失效而自動移除相應節(jié)點的接觸類型。通常，銷蝕掉的節(jié)點的移除會使計算更穩(wěn)定，但是質量的減少會導致錯誤的結果。EQ.0:從計算中移除銷蝕的節(jié)點。（默認）EQ.1 :保留體單元銷蝕的節(jié)點并在接觸中繼續(xù)起作用。EQ.2 :保留體單元和殼單元銷蝕的節(jié)點并在接觸中繼續(xù)起作用?！綰SRSTR】【USRFRC】【NSBCS-【INTERM 【XPENE 【SSTHK】【EC

20、DT】【TIEDPRJ 】9.體積粘度控制。 體積粘度是為了解決沖擊波。* CONTROL. BUL K_VI SOOS I TY01Q2|I0O|1舫 D（/（IGF -1|【Q1】一一缺省的二次粘度系數(shù)（1.5）?！綫2】一一缺省的線性粘度系數(shù)（0.06）。【IBQ】一一體積粘度類型。EQ.-1 :標準。（對于單元類型為 2，10，16的殼單元）EQ.+1 :標準。（默認）15.精確性控制【OSU】第二目標應力（objective stress）更新的全局標志。一般來講，只有像輪胎那樣有大的轉動的顯式計算才會需要該選項。EQ.0 :關。（默認）EQ.1 :開?！綢NN】一一殼單元和體單元節(jié)

21、點的編號方式。EQ.1 :關。（默認為顯式計算）EQ.2 :開。僅對殼單元。（默認為隱式計算）EQ.3 :開。僅對體單元。EQ.4 :開。對殼單元和體單元?！綪IDOSU】一一需要目標應力更新的組件集的 ID號。如果該選項給定，則只有相應的組 件集會應用目標應力更新，OSU選項將設置被忽略。11.全部數(shù)據(jù)輸出* DATABASE_BIWARY_D3PL0TDTLCDT | 忖 O0EAMNPLTC PSET1D ISTATS T57ART LWG| ooo'op o-!rIOOPT0 I【DT】一一輸出的時間間隔。【CYCL】一一時間步內的輸出間隔（一個時間步長是一個循環(huán)）。一般不用。

22、【NOBEAM 】關于是 DATABESE-BINARY-D3POLT 或 DATABESE-BINARY-D3PART 的選擇標志。EQ.0 :被描述成 beam單元的離散的彈簧和減震器單元添加到D3POLT或D3PART的數(shù)據(jù)中。單元的球形坐標xyz還有合力也添加進去。EQ.1 :非不連續(xù)的彈簧和減震器單元添加進去。 舊的數(shù)據(jù)要轉換的 KEYWORD 中時這個選擇被選擇。在舊的數(shù)據(jù)中沒有必須為beam和彈簧單元創(chuàng)建單獨的 ID號，然而彈簧和減震器單元也用 beam表示，這樣就會出錯。EQ.2 :同0 樣。在beam中可以同時出現(xiàn)合力和軸力?！綨OBEAM】EQ.0 :離散彈簧單元和阻尼單元

23、寫入到D3PLOT或D3PART數(shù)據(jù)中，并顯示為梁單元。單元的整體坐標和合力都被寫入到數(shù)據(jù)中。EQ.1 :離散彈簧單元和阻尼單元不會被寫入到D3PLOT或D3PART數(shù)據(jù)中。EQ.2 :與等于0時類似?！綨PLTC】僅用于 D3POLT或D3PART中DT=ENDTIME/NPLTC。這個優(yōu)先于 DT?！綪SETID】 僅用于 D3PART的SET-ID號。ISTATS】一一為選擇的數(shù)據(jù)設定級別。僅用于 D3BEAM【TSTAPT】一一僅用于 D3BEAM，設定模擬開始的時間。默認為0。【IAVG】一一設定寫出平均數(shù)據(jù)的間隔，僅用于 D3BEAM。默認為100【IOOPT】一一僅用于D3PLO

24、T的選擇。EQ.1 :在這時刻每個 plot產(chǎn)生，載荷曲線的值也被加進來到當前的時刻，來決 定下一個plot的時間。這個為默認的。EQ.2:在這時刻每個plot產(chǎn)生，下一個plot的時間T被算出來，T=當前的時間+ 載荷曲線值在T時刻。EQ.3:載荷曲線里每個縱坐標都產(chǎn)生一個plot。曲線的準確值被忽略。12.定義關于動力釋放的卡片。對壓力初始值設定很重要。-CCNrROL.DYNAMl C,_RE LAX AT 1 ONNRCYCKDRTOLDRFCTRDRTERMTSSFDRIREUM_EDT1LIRDFLG10 I0.000 I0.000 10 000 IQ更客解禪桃（fed-1TSSF

25、DR擔嶄！本地段廢詞采劃薛*建這你劉麗緡査詢.【NRCYCK】一一檢驗收斂是的迭代數(shù)目。默認=250?！綝RTOL】收斂公差默認 0.001【DRFCTR】動力釋放因子。默認0.995。【DRTERM】一一終止時間選項為動力釋放。到達這個時間時終止或當收斂時終止?！綯SSFDR】一一在動力釋放中用于計算時間步的比例因子。如果是0,該比例因子為在CONTROL-TIMESTEP 定義的TSSFAC值。收斂后，比例因子被重新設定在TSSFAC中?！綢RELAL】一一基于 Papadrakakis運算法則的自動控制動力釋放選項。EQ.0 :不激活。EQ.1 :激活?！綞DTTL】一一自動控制動力釋放

26、時的收斂公差?！綢DRFLG】一一壓力初始值設定的動力釋放標記。EQ.-999 :即使在一個載荷曲線上指定，動力釋放也不激活，參加DEFINE-CURVE EQ.-1 :動力釋放激活并且時間歷程輸出被延長。EQ.0 :不激活。EQ.1 :激活動力釋放。EQ.2 :指定幾何值的初始化。13.二進制輸出重啟動文檔。定義每個循環(huán)輸出次數(shù)。*DATABASE,BINARY_RUNRSFDT| 10000,00014輸出數(shù)據(jù)控制。指定要輸入到D3PLOT、D3PART、D3THDT文件中的二進制數(shù)據(jù)。*DATABASE,EXTENT.B1 NARYNEIPHNEIPSMAXINT STRFLAG SIG

27、FLGEPSFLG RLTFLG ENGFLG1oi0|2I 1 I 2 I 1 ICMPFLGIEVERPBEAMIP DCXJMPSHGESTSSZ N37HDT0i 1 ro 111I 11NINTSLD【NEIPH】一一寫入二進制數(shù)據(jù)的實體單元額外積分點時間變量的數(shù)目【NEIPS】一一寫入二進制數(shù)據(jù)的殼單元和厚殼單元每個積分點處額外積分點時間變量的數(shù)目?！綧AXINT】一一寫入二進制數(shù)據(jù)的殼單元積分點數(shù)。如果不是默認值3，則得不到中面的結果?！維TRFLAG】一一設為1會輸出實體單元、殼單元、厚殼單元的應變張量，用于后處理繪圖。對于殼單元和厚殼單元，會輸出最外和最內兩個積分點處的張量，

28、對于實體單元，只輸出一個應變張量。【SIGFLG】一一殼單元數(shù)據(jù)是否包括應力張量。EQ.1 :包括（默認）EQ.2 :不包括【EPSFLG 殼單元數(shù)據(jù)是否包括有效塑性應變。EQ.1 :包括（默認）EQ.2 :不包括【RLTFLG 殼單元數(shù)據(jù)是否包括合成應力。EQ.1 :包括（默認）EQ.2 :不包括【ENGFLG 殼單元數(shù)據(jù)是否包括內能和厚度。EQ.1 :包括（默認）EQ.2 :不包括【CMPFLG 實體單元、殼單元和厚殼單元各項異性材料 應力應變輸出時的局部材料坐 標系。EQ.0 :全局坐標EQ.1 :局部坐標【IEVERP 限制數(shù)據(jù)在lOOOstate之內。EQ.0:每個圖形文件可以有不止 1個stateEQ.1:每個圖形文件只能有 1個state【BEAMIP 用于輸出的梁單元的積分點數(shù)。【DCOMP 數(shù)據(jù)壓縮以去除剛體數(shù)據(jù)。EQ.1 :關閉（默認）。沒有剛體數(shù)據(jù)壓縮。EQ.2 :開啟。激活剛體數(shù)據(jù)壓縮。EQ.3 :關閉。沒有剛體數(shù)據(jù)壓縮，但節(jié)點的速度和加速度被去除。EQ.4 :開啟。激活剛體數(shù)據(jù)壓縮，同時節(jié)點的速度和加速度被去除?！維HGE 輸出殼單元沙漏能密度。EQ.1 :關閉（默認）。不輸出沙漏能。EQ.2 :開啟。輸出沙漏能?！維TSSZ輸出殼單元時間步、質量和增加的質量。EQ.1 :關閉。（默認）EQ.2 :只輸出時間步