1、LS-DYNA常見問題匯總1.0資料來源：網(wǎng)絡(luò)和自己的總結(jié)yuminhust2005CopyrightoforiginalEnglishversionownedbyrelativeauthor.ChineseversionownedbyA/Kevin目錄TOC o 1-5 h z HYPERLINK l bookmark4 o Current Document Consistentsystemofunits單位制度2 HYPERLINK l bookmark6 o Current Document MassScaling質(zhì)量縮放4 HYPERLINK l bookmark8 o Current

2、 Document Longruntimes長分析時間9 HYPERLINK l bookmark10 o Current Document Quasi-static準(zhǔn)靜態(tài)11 HYPERLINK l bookmark12 o Current Document 5.Instability計算不穩(wěn)定14 HYPERLINK l bookmark14 o Current Document NegativeVolume負(fù)體積17 HYPERLINK l bookmark16 o Current Document Energybalance能量平衡20 HYPERLINK l bookmark18 o

3、 Current Document Hourglasscontrol沙漏控制27 HYPERLINK l bookmark20 o Current Document Damping阻尼32ASCIIoutputforMPPviabinout37 HYPERLINK l bookmark22 o Current Document ContactOverview接觸概述41 HYPERLINK l bookmark24 o Current Document ContactSoft1接觸Soft=145 HYPERLINK l bookmark26 o Current Document LS-DYN

4、A中夾層板(sandwich)的模擬47怎樣進行二次開發(fā)50Consistentsystemofunits單位制度相信做仿真分析的人第一個需要明確的就是一致單位系統(tǒng)(ConsistentUnits)。計算機只認(rèn)識0&1、只懂得玩數(shù)字，它才不管你用的數(shù)字的物理意義。而工程師自己負(fù)責(zé)單位制的統(tǒng)一，否則計算出來的結(jié)果沒有意義，不幸的是大多數(shù)老師在教有限元數(shù)值計算時似乎沒有提到這一點。見下面LS-DYNAFAQ中的定義：Definitionofaconsistentsystemofunits(requiredforLS-DYNA):1forceunit=1massunit*1accelerationu

5、nit1力單位=1質(zhì)量單位x1加速度單位1accelerationunit=1lengthunit/(1timeunit)人21加速度單位=1長度單位/1時間單位的平方Thefollowingtableprovidesexamplesofconsistentsystemsofunits.Aspointsofreference,themassdensityandYoungsModulusofsteelareprovidedineachsystemofunits.“GRAVITY”isgravitationalacceleration.MASSLENGTTIMFORCESTRESSENERGDEN

6、SITYOUNGsVelocity(56.3KMPH)GRAVITkgkgcmkgcmkgcmkggmgmgmgmtonlbfs2/inslugkgfs2/mkgmmcmcmmmmmmminftmmmm1.e-021.e+04Pamsusmsusms1.e+10KNdyne1.e+071.e-06lbflbfkgfmNJoule7.83E+032.07E+1115.659.8067.83E-032.07E+091.56E+039.81E+02GPady/cm2MbarPaMPaMPapsipsfKN-mmerg1.e7Ncm7.83E-037.83E-037.83E-067.83E+007.8

7、3E+007.83E-03N-mmN-mmlbf-inlbf-ft7.83E-037.83E-097.33E-0415.22.07E+032.07E-032.07E+022.07E+122.07E+002.07E+112.07E+052.07E+053.00E+074.32E+091.561.56E-0315.651.56E+031.56E-031.56E+0415.651.56E+046.16E+0251.339.81E-049.81E-109.81E-039.81E+029.81E-109.81E+039.81E-039.81E+0338632.17kgf/mm2kgf-mm8.02E-1

8、07.00E+021.56E+04(Japan)1000Pa7.83E-062.07E+089.81E+02100000Pgmcmms783E+002.07E+06aMassScaling質(zhì)量縮放質(zhì)量縮放指的是通過增加非物理的質(zhì)量到結(jié)構(gòu)上從而獲得大的顯式時間步的技術(shù)。在一個動態(tài)分析中，任何時候增加非物理的質(zhì)量來增大時間步將會影響計算結(jié)果（因為F=ma）。有時候這種影響不明顯，在這種情況下增加非物理的質(zhì)量是無可非議的。比如額外的質(zhì)量只增加到不是關(guān)鍵區(qū)域的很少的小單元上或者準(zhǔn)靜態(tài)的分析（速度很小，動能相對峰值內(nèi)能非常?。？偟膩碚f，是由分析者來判斷質(zhì)量縮放的影響。你可能有必要做另一個減小或消除了質(zhì)

9、量縮放的分析來估計質(zhì)量增加對結(jié)果的靈敏度。你可以通過人工有選擇的增加一個部件的材料密度來實現(xiàn)質(zhì)量縮放。這種手動質(zhì)量縮放的方法是獨立于通過設(shè)置*Control_timestep卡DT2MS項來實現(xiàn)的自動質(zhì)量縮放。當(dāng)DT2MS設(shè)置為一個負(fù)值時，質(zhì)量只是增加到時間步小于TSSFAC*IDT2MSI的單元上。通過增加這些單元的質(zhì)量，它們的時間達(dá)到TSSFAC*IDT2MSI。有無數(shù)種TSSFAC和DT2MS的組合可以得到同樣的乘積，因而有相同的時間步，但是對于每一種組合增加的質(zhì)量將是不一樣的。一般的趨勢是TSSFAC越小，增加的質(zhì)量越多。作為回報，當(dāng)TSSFAC減小時計算穩(wěn)定性增加（就像在沒有做質(zhì)量縮

10、放的求解中一樣）。如果TSSFAC缺省的值0.9會導(dǎo)致穩(wěn)定性問題，可以試試0.8或者0.7。如果你減小TSSFAC，你可以相應(yīng)增加IDT2MSI，這樣還是可以保證時間步乘積不變。為了確定什么時候和位置質(zhì)量自動增加了，可以輸出GLSTAT和MATSUM文件。這些文件允許你繪出完整的模型或者單獨部件所增加的質(zhì)量對時間的曲線。為了得到由殼單元組成的部件增加的質(zhì)量云圖，將*database_extent_binary卡的STSSZ項設(shè)置為3。這樣你可以用ls-prepost繪出每個單元的質(zhì)量增加量的云圖，具體方法是通過選擇FcompMisctimestepsize。在*control_timestep

11、中設(shè)置DT2MS正值和負(fù)值的不同之處如下：負(fù)值：初始時間步將不會小于TSSFAC*-DT2MS。質(zhì)量只是增加到時間步小于TSSFAC*IDT2MSI的單元上。當(dāng)質(zhì)量縮放可接受時，推薦用這種方法。用這種方法時質(zhì)量增量是有限的。過多的增加質(zhì)量會導(dǎo)致計算任務(wù)終止。正值：初始時間將不會小于DT2MS。單元質(zhì)量會增加或者減小以保證每一個單元的時間步都一樣。這種方法盡管不會因為過多增加質(zhì)量而導(dǎo)致計算終止，但更難以作出合理的解釋。*control_timestep卡中的參數(shù)MS1ST控制是否只是在初始化時增加一次質(zhì)量（MS1ST=1）還是任何需要維持由DT2MS所指定的時間步時都增加質(zhì)量（MS1ST=0）。

12、你可以通過在*control_termination卡片中設(shè)置參數(shù)ENDMAS來控制當(dāng)質(zhì)量增加到初始質(zhì)量一定比率時終止計算（只對自動質(zhì)量縮放有效）可變形點焊梁的質(zhì)量縮放*mat_spotweld卡的質(zhì)量縮放參數(shù)DT只影響點焊單元。如果*control_timestep卡中沒有指定質(zhì)量縮放（DT2MS=0），而且時間由可變形點焊控制，可以用參數(shù)DT來在初始化時增加慣量到點焊單元上來提高時間步達(dá)到DT指定的值。當(dāng)DT不為0時,增加到可變形點焊梁元上的質(zhì)量會輸出到d3hsp文件里oMATSUM中動量和動能不受增加到可變形點焊上的質(zhì)量的影響。GSLTAT中DOES和總的KE受增加的質(zhì)量的影響?？紤]三種

13、調(diào)用可變形點焊的質(zhì)量縮放的情況：1.當(dāng)DT2MS為負(fù)值*mat_spotweld卡DT=0時，盡管在d3hsp文件中可變形點焊質(zhì)量增量百分比不真實。下面幾個值是正確的：d3hsp中”addedspotweldmass”第一個時間步之后的”addedmass”&“percentageincrease”glstat和matsum中的addedmass”。當(dāng)DT2MS為負(fù)值且*mat_spotweld卡DT主0時，可變形點焊質(zhì)量增加不會包含在d3hsp、glstat、matsum文件中的”addedmass”里。這非常容易令人誤解。用戶必須檢查d3hsp文件的”addedspotweldmass”。

14、建議不要同時使用兩種質(zhì)量縮放標(biāo)準(zhǔn)，推薦使用第一種方法（即負(fù)的DT2MS&DT=0）。如果DT2MS=0且DT#0，初始時間步將不考慮增加點焊的質(zhì)量，但是之后每一個周期時間步都會增加10%，直到時間步達(dá)到正確的值（考慮點焊質(zhì)量增加）。glstat&matsum不包含”addedmass”的行。注意質(zhì)量增加會引起能量比率增長。EnglishVersion：Mass-scalingreferstoatechniquewherebynonphysicalmassisaddedtoastructureinordertoachievealargerexplicittimestep.Anytimeyouad

15、dnonphysicalmasstoincreasethetimestepinadynamicanalysis,youaffecttheresults(thinkofF=ma).Sometimestheeffectisinsignificantandinthosecasesaddingnonphysicalmassisjustifiable.Examplesofsuchcasesmayincludetheadditionofmasstojustafewsmallelementsinanoncriticalareaorquasi-staticsimulationswherethevelocity

16、islowandthekineticenergyisverysmallrelativetothepeakinternalenergy.Intheend,itsuptothejudgementoftheanalysttogagetheaffectofmassscaling.Youmayhavetoreduceoreliminatemassscalinginasecondruntogagethesensitivityoftheresultstotheamountofmassadded.Onecanemploymassscalinginaselectivemannerbyartificiallyin

17、creasingmaterialdensityofthepartsyouwanttomass-scale.ThismanualformofmassscalingisdoneindependentlyoftheautomaticmassscalinginvokedwithDT2MSin*control_timestep.WhenDT2MSisinputasanegativevalue,massisaddedonlytothoseelementswhosetimestepwouldotherwisebelessthanTSSFAC*|DT2MS|.Byaddingmasstotheseelemen

18、ts,theirtimestepbecomesequaltoTSSFAC*|DT2MS|.AninfinitenumberofcombinationsofTSSFandDT2MSwillgivethesameproductandthusthesametimestepbuttheaddedmasswillbedifferentforeachofthosecombinations.ThetrendisthatthesmallertheTSSF,thegreatertheaddedmass.Inreturn,stabilitymayimproveasTSSFisreduced(justasinnon

19、-mass-scaledsolutions).IfstabilityisaproblemwiththedefaultTSSFof0.9,try0.8or0.7.IfyoureduceTSSF,youcanincrease|DT2MS|proportionallysothattheproduct/timestepisunchanged.Todeterminewhereandwhenmassisautomaticallyadded,writeGLSTATandMATSUMfiles.Thesefileswillallowyoutoplotaddedmassvs.timeforthecomplete

20、modelandforindividualparts,respectively.Toproducefringeplotsofaddedmassinpartscomprisedofshellelements(DT2MSnegative),setSTSSZ=3in*database_extent_binary.Youcanthenfringetheaddedmass(perelement)usingLS-POSTbychoosingFcompMisctimestepsize.(Here,thelabel“timestepsize”isreallytheelementaddedmass.)Thedi

21、fferencebetweenusingapositiveornegativenumberforDT2MSin*control_timestepisasfollows:Negative:InitialtimestepwillnotbelessthanTSSF*-DT2MS.MassisaddedtoonlythoseelementswhosetimestepwouldotherwisebelessthanTSSF*abs(DT2MS).Whenmassscalingisappropriate,Irecommendthismethod.Theamountofmassthatcanbeaddedu

22、singthismethodislimited.Excessiveaddedmasswillcausethejobtoterminate.Positive:InitialtimestepwillnotbelessthanDT2MS.MassisaddedORTAKENAWAYfromelementssothatthetimestepofeveryelementisthesame.Thismethodishardertorationalizealthoughitisnotsubjecttoterminationfromexcessiveaddedmass.TheparameterMS1STin*

23、control_timestepcontrolswhethermassisaddedonlyonceduringinitialization(MS1ST=1)oranytimeasnecessarytomaintainthedesiredtimestepspecifiedviaDT2MS(MS1ST=0).YoucanuseENDMASin*control_terminationtostopthecalculationafteracertainamountofmasshasbeenadded(activeforautomaticmassscalingonly).Mass-scalingofde

24、formablespotweldbeams:Themass-scalingparameterin*mat_spotweld(DT)affectsonlythespotwelds.Ifnomass-scalingisinvokedin*control_timestep(DT2MS=0.)ANDthetimestepiscontrolledbythedeformablespotwelds,DTcanbeusedtoaddinertiatothespotweldsduringintializationinordertoincreasethetimesteptoavalueofDT.WhenDTisn

25、onzero,massaddedtospotweldbeamsisreportedtod3hsp.MATSUMmomentumandKEdoesNOTfactorinaddedmasstodef.spotwelds.GLSTATDOESfactorinaddedmasstototalKE(spotweld.beam.type9.mscale.initvel.k)Consider3casesofinvokingmass-scalinginamodelwithdeformablespotwelds:1.Although“percentagemassincrease”under“Deformable

26、Spotwelds:”ind3hspisboguswhenDT2MSisneg.andDTin*mat_spotweld=0,thefollowingarecorrect:“addedspotweldmass”ind3hsp“addedmass”and“percentageincrease”ind3hspAFTERthefirsttimestep“addedmass”inglstatandmatsum2.AddedspotweldmasscontrolledbyDTin*mat_spotweldisNOTINCLUDEDin“addedmass”givenind3hsp,glstat,orma

27、tsumwhenDT2MSisneg.andDTin*mat_spotweldisnonzero.Thiscanbequitemisleading.Usermustcheckfor“addedspotweldmass”ind3hsp.Recommended:Donotinvokebothmass-scalingcriteria.Neg.DT2MSwithDT=0(case1above)ispreferred.3.IfDTisnonzeroandDT2MS=0,theinitialtimestepwillNOTconsideraddedspotweldmassbutthetimestepwill

28、increaseby10%eachcycleuntilthecorrecttimestep(consideringaddedspotweldmass)isachieved.Glstatandmatsumcontainno“addedmass”lineitem.Theabovecanbeillustratedusing/j5000a_2/jday/test/weld/spotweld.beam.type9.mscale.k.Notethataddedmassmaycausetheenergyratiotorise.(See/j5000a_2/jday/test/erode/taylor.mat3

29、.noerode.mscale.k)Longruntimes長分析時間當(dāng)用顯式時間積分時，對于仿真非常小的部件而分析時間又要相當(dāng)長時沒有好的方法。質(zhì)量縮放(mass-scaling)增加了需要確認(rèn)非物理質(zhì)量的增加不會顯著影響計算結(jié)果的負(fù)擔(dān)。當(dāng)使用時間縮放(time-scaling)時也有同樣的問題。時間縮放(time-scaling)是指為了減小需要的時間步數(shù)，通過增加加載速率而縮短仿真時間。要確認(rèn)時間步不是僅由很少的小單元或者剛度大單元控制，可以通過在d3hsp文件中搜索”smallest”來顯示100個最小的時間步單元。如果只有很少的幾個單元控制時間步，可以把那些單元及鄰近區(qū)域重新reme

30、sh或者把它們變成剛體?？墒莾H運行必要長的時間是很明顯的。這意味著在一個跌落分析的情況時，給跌落物體一個初速度，把它放在離地面一個非常小的距離。沖擊之后，僅運行足夠得到需要的結(jié)果的時間。值得注意的是對于一個長時間的仿真，如果時間步數(shù)超過了50萬步，最好使用雙精度版本的LS-DYNA求解器，使截斷誤差的影響最小化。運行雙精度版本要增加30%的時間。對于長時間的分析，自動顯式/隱式轉(zhuǎn)換可能是一個選擇。使用這種方法，用戶可以指定在一個時間段內(nèi)使用隱式積分。隱式積分的優(yōu)點是時間步不由單元尺寸控制，所以可以得到大的時間步。當(dāng)然，隱式計算也非常點用cpu時間。而且，目前并不是所有的LS-DYNA的功能和材

31、料都在隱式分析中實現(xiàn)(大部分已經(jīng)實現(xiàn))。下面的FEAinformationnewsletter里討論了顯式/隱式轉(zhuǎn)換()。Seealso:mass_scaling,quasistaticEnglishVersion：Whenyoureusingexplicittimeintegration,thereisnomagiccureforlongruntimesassociatedwithsimulatingverysmallgeometriesoverrelativelylongperiodsoftime.Mass-scalingcarriesaburdenofhavingtoconfirmtha

32、ttheadditionofnonphysicalmassdoesnotsignificantlyaffecttheresults(seeattachedfile“mass_scaling”).Asimilarburdenexistswhentime-scalingisemployed.Time-scalingisatechniquewheretheloadingrateisincreasedandthusthesimulationtimeisshortenedinordertoreducetherequirednumberoftimesteps.Makesurethatyourtimeste

33、pisnotbeingcontrolledbyonlyafewsmallorstiffelementsbysearchinginthed3hspthestring“smallest”.Ifthereareonlyafewcontrollingelements,youcanremeshinthevicinityofthoseelementsorperhapsmakethemrigid.Thoughitsratherobvious,runonlyaslongasisnecessary.Thismeansinthecaseofadropsimulation,assigninganinitialvel

34、ocitytothedroppedobjectandplacingitaverysmalldistancefromthelandingsurface.Afterimpact,runonlylongenoughtogettheresultsyouneed.Beawarethatforlengthysimulationswherethenumberoftimestepsgoesabovehalfamillionorso,youdbewelladvisedtouseadoubleprecisionexecutableofLS-DYNAtominimizeerrorduetoroundoff.Runn

35、ingdoubleprecisioncarrieswithitacpupenaltyofaround30%.Automaticexplicit/implicitswitchingmaybeanoption.Usingthistechnique,thep1a1alaaattusercanspecifytimewindowsinwhichimplicittimeintegrationisusedasopposedp-4pTtoexplicittimeintegration.Anadvantageofimplicittimeintegrationisthattimestepsarenottiedto

36、elementsizeandcanthusbemuchlarger.Ofcourse,animplicittimestepisalsomuchmoreexpensiveintermsofcpu.Further,notallLS-DYNAfeaturesandmaterialsareimplementedforimplicitanalysisatthistime(thoughmostare).Explicit/implicitswitchingisdiscussedinthefollowingarchivedFEAInformationnewsletter.Seealso:mass_scalin

37、g,quasistatic.Quasi-static準(zhǔn)靜態(tài)動態(tài)松馳(Dynamicrelaxation)并不是有意為一般的準(zhǔn)靜態(tài)(quasi-static)分析設(shè)置的。它適合于當(dāng)預(yù)載只產(chǎn)生小的彈情況應(yīng)變的施加預(yù)載，或者初始化系統(tǒng)到一個預(yù)定義的幾何形狀1。但對其它更多情況并不適合。你可以通過做一個常規(guī)的顯示仿真來模擬準(zhǔn)靜態(tài)分析，通過按需要調(diào)用時間/質(zhì)量縮放(time-scaling,mass-scaling)來在可接受的時間內(nèi)得到結(jié)果，但這種方法是需要技巧地。你必須監(jiān)測系統(tǒng)動能按希望的使慣性效應(yīng)最小化?；旧蟿幽芟鄬?nèi)能應(yīng)該保持在一個較小的值。時間縮放是指加載比在準(zhǔn)靜態(tài)實驗里更快，以減少總的仿真

38、時間。關(guān)于質(zhì)量縮放更多內(nèi)容可以看”mass_scaling”一節(jié)。或者你可以嘗試用LS-DYNA運行一個隱式靜力分析?？梢钥从脩羰謨岳锏目ㄆ?control_implicit_.和AppendixM。Seealso:gravity.txt,readme.preload,mass_scaling,long_run_times,implicit.general,quick_initialization.Notel：初始化到預(yù)定義的幾何從第一次分析的最終狀態(tài)輸出一個節(jié)點位移文件。(這一部分未按原文翻譯)注意d3plot文件里不包含節(jié)點轉(zhuǎn)動信息，因此轉(zhuǎn)動輸出為0。這對初始化殼和梁單元會是個問題。LS-

39、Prepost有一個選項是輸出節(jié)點位移，在Output-NodalDisplacements里。但是這個輸出是i8,3e16格式的，但需要的是i8,3e15，所以要注意修改一下。如果你做了一個正常的動態(tài)松馳分析來得到初始狀態(tài)，一個預(yù)定義位移和轉(zhuǎn)動的drdisp.sif文件在DR階段結(jié)束時會自動創(chuàng)建。在第二次分析時，快速的初始化到第一步輸出的預(yù)定義的幾何。你需要設(shè)置卡片*control_dynamic_relaxation里的參數(shù)IDRFLG=2，而且在命令行里指定”m=”(其中指第一步創(chuàng)建的文件)。這樣在瞬態(tài)分析之前，LS-DYNA會自動做一個100步的預(yù)分析來使節(jié)點根據(jù)文件指定的數(shù)據(jù)移動到指

40、定值。EnglishVersion:Dynamicrelaxationisnotintendedforgeneralquasistaticanalysis.Itsokforapplyingpreloadwhenthepreloadproducesonlysmallelasticstrainsorforinitializingasystemtoaprescribedgeometry1butitsnotgoodformuchelse.Youcandoaquasi-staticanalysisbyrunningaregularexplicitsimulation,invokingtime-and/o

41、rmass-scalingasnecessarytocrankouttheresultsinareasonabletimeframe,butthisapproachcanbetricky.Youhavetokeepaneyeonthekineticenergyinthesystemasyouwanttominimizetheinertialeffects.Basically,thekineticenergyshouldremainsmallrelativetotheinternalenergy.(Bytime-scaling,Imeanapplyingtheloadmorequicklytha

42、ninthequasi-staticexperimentinordertoreducethesimulationtime.)Seethefile“mass_scaling”formoreonmass-scaling.Or,youcantryanimplicit,staticanalysisusingLS-DYNA.Seethecommands*control_implicit_.andAppendixMintheUsersManual.Thereareexamplesofimplicitanalysisonour“user”inthels-dyna/exampledirectory.Seeal

43、so:gravity.txt,readme.preload,mass_scaling,long_run_times,implicit.general,quick_initialization.Note1*Initializingtoaprescribedgeometry*1.Writeanodaldisplacementsfromthefinalstateofyourfirstrun.Togetthisdatainthenecessaryformat,useLS-TAURUSasfollows:ls-taurusg=d3plotexecutes1000goestofinalstatedefor

44、mwriteadescribedabovetNodalDisplacementsbuttheoutputisi8,3e16ratherthantherequiredi8,3e15andhencethesuggesteduseofLS-TAURUS.LS-TAURUSisnotavailableforWindowsPCs.ItsfreeandavailableforUnixandLinuxworkstations.Ifyoudoaregulardynamicrelaxationruntogettotheinitializedstate,aprescribeddisplacementsandrot

45、ationswillautomaticallybewrittenattheconclusionoftheDRphase(drdisp.sif).Bug#2020reportedon9/22/2004thatrigidbodynodesdonotgetinitializedaccordingtodatain“m=pres_geom_file”.Additionalexamplein/home/jday/test/cantilever/solid/typ2sol_dr_nrb.k(createsdrdisp.sif)andtyp2sol_presgeom_nrb.k(m=drdisp.sifrun

46、).Nodes11,22,33,44arenotinitializedtowhatsindrdisp.sif2.Inyoursecondrun,quicklyinitializetotheprescribedgeometrywritteninstep1.Youneedtoset:“m=”ontheexecutionIDRFLG=2in*control_dynamic_relaxationandincludelinewhere“”istheinstep1.Beforethetransientrunbegins,LS-DYNAwillautomaticallyrunaprecusoranalysi

47、sof100timestepswhereinthenodsearedisplacedaccordingtothedatain“”.5.Instability計算不穩(wěn)定一些表示計算不穩(wěn)定的消息如：“out-of-rangevelocities”速度超出范圍“negativevolumeinbrickelement”體單元負(fù)體積“terminationduetomassincrease”因質(zhì)量增加而終止用來克服顯式求解中的不穩(wěn)定的方法如下：首先（也是最重要的）是使用可獲得的最新的LS-DYNA版本。最新的執(zhí)行塊可以從上下載（注：前提是你有訪問權(quán)限）。聯(lián)系LSTC獲得user帳號的密碼。最新的BET

48、A版執(zhí)行塊可以在上找到（不需要密碼，但lstc公司對ftp訪問有IP限制）。其次是增加d3plot的輸出頻率到可以顯示出不穩(wěn)定的出現(xiàn)過程。這可以提供導(dǎo)致不穩(wěn)定性發(fā)生的線索。其它的不些解決數(shù)值不穩(wěn)定性的技巧：*試著用雙精度LS-DYNA版本運行一次*試著減小時間步（timestep）縮放系數(shù)（即使使用了質(zhì)量縮放mass-scaling）*單元類型和/或沙漏（hourglass）控制。對出現(xiàn)不穩(wěn)定的減縮體和殼單元，試著用沙漏控制type4和沙漏系數(shù)0.05?；蛘咴囍妙愋?6的殼單元，沙漏控制type8。如果殼響應(yīng)主要是彈性，設(shè)置BWC=1和PROJ=1（僅對B-T殼）。避免使用type=2體單元

49、。對體單元部件，在厚度方向最少用兩個體單元。*接觸。設(shè)置接觸的bucketsorts之間周期數(shù)為0，這樣會使用缺省的分類間隔。如果參與接觸的兩個部件的相對速度異常的大，可能需要減小bucketsort的間隔（比如減小到5,2甚至1）。如果仿真過程中有明顯的接觸穿透出現(xiàn)，轉(zhuǎn)換到使用*contact_automatic_surface_to_surface或者*contact_automatic_single_surface，并設(shè)置SOFT=1。確保幾何考慮了殼單元的厚度。如果殼非常薄，比如小于1mm，放大或者設(shè)置接觸厚度到一個更加合理的值。*避免冗余的接觸定義，也就是說不要對同樣的兩個部件定義多

50、于一個的接觸對。*查找出現(xiàn)不穩(wěn)定的部件的材料定義中的錯誤(比如誤輸入，不一致的單位系統(tǒng)等)關(guān)掉所有的*damping這些技巧是一些通用的方法，可能并不適合于所有的情況。Seealso:negative_volume_in_brick_element.tips，shooting-nodesEnglishVersion:Somemessagesthatindicateaninstabilityhasoccurred:“out-of-rangevelocities”“negativevolumeinbrickelement”“terminationduetomassincrease”Approach

51、estocombatinginstabilityofanexplicitsolution:Firstandforemost,usethelatestversion/revisionofLS-DYNAavailable.Thelatestproductionexecutablescanbedownloadedfrom.ContactLSTCforthepasswordtothis“user”.MorerecentBETAexecutablesarefoundin(nopasswordrequired).Thenextstepistowriteplotstatesfrequentlyenought

52、oseetheevolutionoftheinstability.Thisshouldoffercluesintowhatsinitiatingtheinstability.Someothergeneraltipstowardresolvingnumericalinstabilities:TryrunningadoubleprecisionLS-DYNAexecutable.*Timestep.Tryreducingthetimestepscalefactor(evenifmass-scalingisinvoked).*Elementformulationand/orhourglasscont

53、rol.Forunderintegratedsolidsorshellsthatgounstable,tryhourglasstype4withahourglasscoefficientof0.05.Or,tryshellformulation16withhourglasstype8.Ifresponseofshellsisprimarilyelastic,setBWC=1andPROJ=1(B-Tshellsonly).Avoidtype2solids.Useatleasttwosolidelementsthruthethicknessofanysolidpart.*Contact.Setn

54、umberofcyclesbetweenbucketsortstozerosothatthedefaultsortintervalwillbeused.Iftherelativevelocitybetweentwopartsincontactisexceptionallyhigh,itmaybenecessarytoreducethebucketsortinterval(forinstanceto5,2,oreven1).Ifvisiblecontactpenetrationsdevelopduringthesimulation,switchto*contact_automatic_surfa

55、ce_to_surfaceor*contact_automatic_single_surfacewithSOFTsetto1.Makesuregeometrytakesintoaccountthicknessofshells.IfshellsareVERYthin,e.g.,lessthan1mm,scaleuporsetthecontactthicknesstoamorereasonablevalue.Avoidredundantcontactdefinitions,thatis,donttreatcontactbetweenthesametwopartsusingmorethanoneco

56、ntactdefinition.*Lookformistakes(typos,inconsistentunits,etc.)inmaterialinputofpartsthatgounstable.*Turnoffall*damping.Thesetipsareofageneralnatureandmaynotbeappropriateinallsituations.Seealso:negative_volume_in_brick_element.tips,shooting-nodesNegativeVolume負(fù)體積泡沫材料的負(fù)體積(或其它軟的材料)對于承受很大變形的材料，比如說泡沫，一個單

57、元可能變得非常扭曲以至于單元的體積計算得到一個負(fù)值。這可能發(fā)生在材料還沒有達(dá)到失效標(biāo)準(zhǔn)前。對一個拉格朗日(Lagrangian)網(wǎng)格在沒有采取網(wǎng)格光滑(meshsmoothing)或者重劃分(remeshing)時能適應(yīng)多大變形有個內(nèi)在的限制。LS-DYNA中計算得到負(fù)體積(negativevolume)會導(dǎo)致計算終止，除非在*control_timestep卡里面設(shè)置ERODE選項為1,而且在*control_termination里設(shè)置DTMIN項為任何非零的值，在這種情況下，出現(xiàn)負(fù)體積的單元會被刪掉而且計算繼續(xù)進行(大多數(shù)情況)。有時即使ERODE和DTMIN換上面說的設(shè)置了,負(fù)體積可能

58、還是會導(dǎo)致因錯誤終止。有助于克服負(fù)體積的一些方法如下：*簡單的把材料應(yīng)力應(yīng)變曲線在大應(yīng)變時硬化。這種方法會非常有效。有時候修改初始網(wǎng)格來適應(yīng)特定的變形場將阻止負(fù)體積的形成。此外,負(fù)體積通常只對非常嚴(yán)重的變形情況是個問題,而且特別是僅發(fā)生在像泡沫這樣的軟的材料上面。*減小時間步縮放系數(shù)(timestepscalefactor)。缺省的0.9可能不足以防止數(shù)值不穩(wěn)定。避免用全積分的體單元(單元類型2和3),它們在包含大變形和扭曲的仿真中往往不是很穩(wěn)定。全積分單元在大變形的時候魯棒性不如單點積分單元,因為單元的一個積分點可能出現(xiàn)負(fù)的Jacobian而整個單元還維持正的體積。在計算中用全積分單元因計算

59、出現(xiàn)負(fù)的Jacobian而終止會比單元積分單元來得快。*用缺省的單元方程(單點積分體單元)和類型4或者5的沙漏(hourglass)控制(將會剛化響應(yīng))。對泡沫材料首先的沙漏方程是：如果低速沖擊type6,系數(shù)1.0;高速沖擊type2或者3。*對泡沫用四面體(tetrahedral)單元來建模，使用類型10體單元。*增加DAMP參數(shù)(foammodel57)到最大的推薦值0.5。*對包含泡沫的接觸，用*contact選項卡B來關(guān)掉shootingnodelogic。使用*contact_interior卡用partset來定義需要用contact_interior來處理的parts，在set

60、_part卡1的第5項DA4來定義contact_interior類型。缺省類型是1，推薦用于單一的壓縮。在版本970里，類型1的體單元可以設(shè)置type=2,這樣可以處理壓縮和減切混合的模式。*如果用mat_126，嘗試ELFORM=0*嘗試用EFG方程(*section_solid_EFG)。因為這個方程非常費時，所以只用在變形嚴(yán)重的地方，而且只用于六面體單元。Seealso:instabilityEnglishVersion:NegativeVolumesinFoams(orothersoftmaterials)Inmaterialsthatundergoextremelylargedef