在LS-Prepost建立sph模型的詳細(xì)過(guò)程

1、LS Prepost關(guān)鍵操作1 模型建立SphGenSWGenBlockNSurMeshnLlleshSmoothZDmeshBMeshBulkFTetri eshBestFitCurvesSurfaceChainMTTMeshMasslrDie LineH1P201SQualP eneCknHSShowGMrph1 I 2 I 3 I 4 I 5 I 6D圓柱體模型圓柱底部中心坐標(biāo)SphGenSWGenBkcknSurMeshnLMeshSmooth2Dme&hBMeshB-ulkFTetMeshBestFrtCurvesSurfaceChainMTTMeshMassTrDielineH1P

2、201SQualPenetkHHSShowGMorph1234CPositionPSph generation interface CreateMcxJify圓柱半徑圓臺(tái)/錐模型5v：oDYr 0D2y1Len;10 Polar10P1D; 2MXNID: 801Direction方向坐標(biāo)圓柱高M(jìn)YlDen: 7-6Fil%100cknSurMeshnLMeshSmoothZDmeshBHeshBulkFTetNeshBestFrtCurvesSurfaceChainMTTMeshHassTrDieLineHIP201SQualP enenP Sph generation interface

3、CreateModify00圓臺(tái)底部半徑dl: 5DirectkmDDYDZLen:方向坐標(biāo)圓臺(tái)咼圓臺(tái)頂部半徑2: 8PIO： 1O PolarNIO: 1NX WNV 10NZ 20en: 78100-0軸向節(jié)點(diǎn)數(shù)數(shù)據(jù)庫(kù)2參數(shù)定義*DATABASE_FORMAT計(jì)算結(jié)果輸出格式爾 Keyword Input Formar* Airbag*Dba5e*HatAle* Define*Node*Boijndry*Elem* Pa ram*Cii5tmd*Eos*Part*Compnt*Hrglass*RgdwQl*CoiitactSectionControl*Intgrtn*Set*Ief2Rg*

4、TermntSam ping*User1 I 3 J 4 I 5 I 6 I 7 I D -Keywora*IWCTABA5EEdit0 Use *PARAMnIFORHCOMMENT:IBINARYHEAllDoneModelAccept Delete Default Done-DATABASE_FORHAT (0(5uEsys: 1) SettingQ, , 1, 11】11,1, 2,Q, I, 31, I 4A, 1115A, i】&A ii, ZQi, 3輸岀結(jié)果格式ijfsrSStSOR4$5 DATABASE 7OIMA7i063VaiubleIFORMffiDWEYTypeII

5、De&ult0012Card2EQ、EQ、EQ、IBINARYIFORM Output format for D3PLOT and D3THDT files0: LS-DYNA database format (default),I: ANSYS database format,2: Both LS-DYNA and ANSYS database formatsWord size of the binary output files (D3PLOT. D3THDT. D3DRLF andinterface files for 64 bit puter such as CRAY and NEC、

6、EQ、0: default 64 bit format,EQ、1: 32 bit IEEE formatRemarks:1、This option is not available for every platform Check LS-DYNA Banner upon executionof the program2、By using this option one can reduce the size of the binary output files which are created by 64 bits puter such as CRAY and NEC、二進(jìn)制*DATABAS

7、E_EXTENT_BINARY 結(jié)果輸出控制數(shù)據(jù)庫(kù)* Airbag*fbBS*Mat* Define*Uode*Bouiidry*Elem* Pa ram*Ciistmd*Eos*Part*Coinpnt*Hrgla*Rgdwal*ContactInitial*SectKnControl*Intgrtn*Set*Def2Rg*Intrfac*TermntPam pingRLoad*User Keyword Input FormNewIDAddAcceptDeleteDefaultDofusePARAMETER(Subsys 1Settv*DArABA5t_tXIWT_BinAKY (0)_3

8、M I 5-K&ywor3*O5GABASEEditDoneAllModelBINARY94342 W343 9434i94軸94346dPHHOPSHAXmTSTRFIGSIGFLGEPSFLGRLTFL6OKHG003011巾h巾 -1CHPFLClEVERPBEAHtPOCOHPSHCE5TSSZN3TrtTWLEHAT1|o0*1*I11b巾 JNINTSLDPKPSOSOPUNUSEDHssaTHERHWTOOTfWOOUTlTqt|io0(l-jsTRESS|srRS5 JCOHHHH:xqB,q.34一赳Q，5Q 飛Q加&$DATABASE OPTIONS$住住住住住住住住住住住

9、住住住住住住住 住 住 住 住 住住住住住住住住住住住住住住住住住住住住住住 住 住 住 住 住bbbbpnbbbbb$*SATAASE EXTENT 町NKRY00000,VajiablfiNEIPHNEPSlAHNTSTEFLGSIGFLGEPSFLGRLTTLGO For solid elements a single strain tensor is written Flag for including stress tensor in the shell LS-DYNA database:I: include (default),2: exclude、Flag for includi

10、ng the effective plastic strains in the shell LS-DYNA database: 1: include (default),2: exclude、SIGFLGEQ、 EQ、EPSFLGEQ、EQ、RLTFLG Rag for including stress resultants in the shell LS-DYNA database:EQ、I: include (default),EQ、2: exclude、ENGFLG Flag for including shell internal energy density and thicknes

11、s in the LSDYNA database:EQ、I: include (default),EQ、2: exclude、CIVIPFLG Orthotropic and anisotropic material stress and strain output in local material coordinate system for solids, shells and thick shells EQ、0: globakEQ、1: local AlEVERP Every plot state for “d3plot” database is written to a separat

12、e file a This option will limit the database to 1000 states:EQ、0: more than one state can be on each plotfile,EQ、I: one state only on each plotfileBEAMIP Number of beam integration points for out put、This op tion does not apply to beams that use a resultant formulation DP Data pression to eliminate

13、rigid body data:EQ、I: off (default), no rigid body data pression,EQ、2: on. rigid body data pression active,EQ、3: off. no rigid body data pression, but nodal velocities and accelerations are eliminated from the database、EQ、4: on, rigid body data pression active and nodal velocities and accelerations

14、are eliminated from the database、SHGE Output shell hourglass energy density:EQ、I: off (default), no hourglass energy written,EQ、2: on、STSSZ Output shell element time stepnass, or added mass:EQ、1: off (default),EQ、2: out put time ste psize,See remark 3 below、EQ、3: output mass, added mass, or time ste

15、p size、N3THDT Material energy write option for D3THDT databaseEQ、I: off. energy is NOT written to D3THDT database,EQ、2: on (default), energy is written to D3THDT database lALEMAT Output solid part ID list containing ale materials EQ、1: on (default)NINTSLD Number of solid element integration points w

16、ritten to the LS-DYNA database The default value is 1、 For solids with multiple integration points NINTSLD may be set to 8、 Currently, no other values for NINTSLD are allowed、 For solids with multiple integration points, an average value is output if NINTSLD is set to 1、PKP_SEN Flag to output the pe

17、ak pressure and surface energy puled by each contact interface into the interface force database To obtain the surface energy, FRCENG, must be sent to 1 on the control contact card、 When PKP_SEN=1, it is possible to identify the energies generated on the upper and lower shell surfaces, which is impo

18、rtant in metal forming appli- cations、 This data is mapped after each H-adaptive remeshingEQ、0: No data is writtenEQ、I: Output the peak pressures and surface energy by contact interfaceSCLP A scaling parameter used in the putation of the peak pressure、 This parameter isgenerally set to unity (the de

19、fault), but it must be greater than 0、MSSCL Output nodal information related to mass scaling into the D3PLOT database This option can be activated if and only if DT2MS 0) may make the database inpatible with other 3rd party software、EQ、EQ、EQ、EQ、 Remarks:0: (default) output ternperatureI: output temp

20、erature2: output temperature and flux3: output temperature, flux, and shell bottom and top surface temperature1、If MAXINT is set to 3 then inid-surface, inner-surface and outer-surface stresses are output at the center of the element to the LS-DYNA database、 For an even number of integration points,

21、 the points closest to the center are averaged to obtain the midsurface values If multiple integration points are used in the shell plane, the stresses at the center of the element are found by puling the average of these points、 For MAXINT equal to 3 LS-DYNA assumes that the data for the user defin

22、ed integration rules are ordered from bottom to top even if this is not the case、 If MAXINT is not equal to 3. then the stresses at the center of the element are output in the order that they are stored for the selected integration rule、 If multiple points are used in plane the stresses are first av

23、eraged、2、Beam stresses are output to the LS-DYNA database if and only if BEAMIP is greater than zero In this latter case the data that is output is written in the same order that the integration points are defined、 The data at each integration point consists of the followingfive values for elastic-p

24、lastic Hughes-Liu beams: the normal stress, frr; the transverse shear stresses, frs and ftr： the effective plastic strain, and the axial strain which is logarithmic、 For beams that are not elastic-plastic, the first history variable, if any, is output instead of the plastic strain For the beam eleme

25、nts of Belytschko and his coworkers, the transverse shear stress ponents are not used in the formulationNo data is output for the Belytschko-Schwer resultant beam、3、If mass scaling is active, the output of the time step size reveals little information about the calculation If global mass scaling is

26、used for a constant time step, the total element mass is output; however, if the mass is increased so that a ininiinuin time step size is maintained (DT2MS is negative), the added mass is output、 Also, see the control card CONTROL _TIMESTEP、數(shù)據(jù)庫(kù)*DATABASE_BINARY_D3PLOT結(jié)果輸出時(shí)間步長(zhǎng)* Airbag1 *Dtase* Define*

27、Mat*Ale*Node*foundry*Elem* Pa ram*Cnstmd*Eos*Part*CompntHrglass*RgdvralContactInitial*SectionControlTntgrtnSet*Def2Rg*IntrfacTerrnnt* Damping*LoadUser1訂3)4丨5 16 1 7 1 DP KeyworJM詞ABASEEditDoneAllModel2斷 Keyword Input Form0 Use PARAMFTDII8PTIXDTPick Accept Delete Defai*DAIABASt_BlNARY_D3PLOT ( 0 )BEA

28、MNPLTC0,1 八 II I1 口 I,11i i2(C：I I勺Q(Subsys: 1)P5ET1D)3PART)3PROP53532 $riHE HIST3Kf$5i；33 轉(zhuǎn)$報(bào)$；$擲畀$畀$鶉$9$ $M,35 *5A7AIASEjnAJlY D3PLOT93：3 2,000輸岀時(shí)間步長(zhǎng)VariableDT/CYCLLCDT/NRBEAMNPLTCPSETIDTypeF1IIIDefault-Remarks23567$4CardOptional Card ihal ouh npp lies to the DSP LOT databaseVariablelOOPTType1Defa

29、ult0RemarksCard24578Time interval between outputs、DTCYCL Output interval in time steps (a time step is a cycle) For the D3DRFL file a positive number n will cause plot dumps to be written at every nth convergence check interval specifled on the *CONTROL_DYNAMIC_RELAXATION card、NR Number of Running R

30、estart Files, RUNRSF, written in a cyclical fashion、The default number is one, i、e、 the same overwritten each time、LCDT Optional load curve ID specifying time interval between dumps、 This option is only available for the D3PLOT. D3PART, D3THDT and INTFOR files、BEAM Option flag for *DATABASE_BINARY_D

31、3PL0T or D3PART、EQ、0: Discrete spring and dumper elements are added to the D3PLOT or D3PART database where they are display as beam elements The element global X, global Y, global Z and resultant forces are written to the database,EQ、I: No discrete spring and damper elements are added to the D3PLOT

32、or D3PART database A This option is useful when translating old LS-DYNA input decks to KEYWORD input、In older input decks there is no requirement that beam and spring elements have unique IDs, and beam elements may be created for the spring and dampers with identical IDs to existing beam elements ca

33、using a fatal error Contact interfaces which are based on part IDs of seat- belt elements will not be properly generated if this option is used、EQ、2: Discrete spring and damper elements are added to the D3PLOT or D3PART database where they are displayed as beam elements (similar to option 0)、 In thi

34、s option the element resultant force is written to its first database position allowing beam axial forces and spring resultant forces to be plotted at the same time、 This can be useful during some postprocessing applications NPLTC DT二ENDTIME/NPLTC applies to D3PLOT and D3PART only、This overrides the

35、 DT specified in the first field、PSETID SET.PART ID for D3PART only、 lOOPT This option applies to the D3PLOT、Flag to govern behavior of the plot fre- quency load curve defined by LCDT:EQ、I: At the time each plot is generated, the load curve value is added to the current time to determine the next pl

36、ot time (this is the default behavior)、EQ、2: At the time each plot is generated, the next plot time T is puted so that T = the current time plus the load curve value at time T、EQ、3: A plot is generated for each abscissa point in the load curve definition The actual value of the load curve is ignored

37、control sphXCBSBOWDTIDEMMEMORYFORMSTART1AXVTypeIIFIIIFFDefauli101 e20oeae1500OO1nCSIE4|CtIntgrtn*Set*Intrfac*Tc：riTint*Donipinq*Lgd*Uxr1 S 16 1 7 1 D-Key W OJ/CpIfTRO L Wrt 1Done 1ModelAlly Kewad Input FormdevDefaJtAcceptMeteDoneWUsexpAWHETER(Sabsysl)SeUng*C(HlTROt_SPH(0)BonoDIHEHORYFORMSTMTfWWI LE+

38、200.010E-H5ripR MRR MPR MPR ripR MPP MPR MPP ripR MPPDECONPOsmon scale, ECOMPOSmON_SCALE D ECOM POSmON_SHOW DtCOMPOSIT lOri_TRAMS IO_DINOUTOriLy IO_LSTC_REDUCE lO_MOD3rUMP IO_ hlODlIMP IO_MOFULLIO SWAPBVTESHON toe Al OUTPUT PARALLEL RORE_FLUID REMESHIOG RIGID SHELL SOLIDSPOTWELD BEAMDfRIVHE(WHEW；11.

39、 X.Q . 111.11 2|Q 1111111 3|Q 1111111 4|Q . . 5.0 . . 6.0 . . 7.0 i i 23 C0NTR0L SPH 0 02425OLOOOOE+20 01500. OOOL OOOOE+15VanableCOXTDERTV環(huán)IDefiul!00NCBS Number of cycles between particle sorting 粒子分類搜索循環(huán)次數(shù)BOXID SPH approximations are puted inside a specified BOX、 When a particle has gone outside t

40、he BOX. it is deactivated This will save putational time by eliminating particles that no longer interact with the structure、指定BOX內(nèi)得SPH粒子參與汁算。、“1某個(gè)SPH粒子位于BOX Z外時(shí)，該粒子失效，通過(guò) 消除某些不再與結(jié)構(gòu)發(fā)生作用得粒子,可以節(jié)省計(jì)算時(shí)間。DT Death time、Detenniiies when the SPH calculations are stopped粒子失效時(shí)間IDIM Space dimension for SPH parti

41、cles: SPH 粒子得空間維蛟3 for 3D problems參數(shù)為3時(shí)，表示3維問(wèn)題。2 for 2D plane stiain problems 參數(shù)為 2 時(shí)，表示 2 維問(wèn)題。-2 for 2D axisyninietric problems參數(shù)為-2時(shí)，表示軸對(duì)稱問(wèn)題。When a value is not specified LS-DYNA determines the space dimension automatically by checking the use of 3D, 2D or 2D asisyninietric elements打這個(gè)值無(wú)法自動(dòng)指定LS-DY

42、NA得維數(shù)，程序通過(guò)核對(duì)使用3維、2維或軸對(duì)稱得 單元來(lái)確定空間維數(shù)。MEMORY Defines the initial number of neighbors per particle、This variable is just for memory allocation of arrays during the initialization phase During the calculation, some particles can request more neighbors and LS-DYNA will automatically adapt the size of that

43、 variable Default value should apply for most applicaUons、定義每個(gè)粒子得初始相鄰粒子得數(shù)量，該變量只就是在初始化階段調(diào)整內(nèi)存分配。在計(jì)算 中，如果某些粒子需要更多得相鄰粒fXS-DYNA會(huì)自動(dòng)調(diào)整該變ft設(shè)置。默認(rèn)值適用于 大部分問(wèn)題。FORM Particle approximation theory:粒子近似理論EQ、0: default fbnnulation, 默認(rèn)公式EQ、 1: reniomialization approximationifi歸-近似化START Start time for particle approximation、 Particle approxiinations will be put