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

1、LS_Prepost關(guān)鍵操作1.模型建立立方體模型立方體Part編號(hào)起始節(jié)點(diǎn)號(hào)：該P(yáng)art第一個(gè)節(jié)點(diǎn)號(hào)碼起始坐標(biāo)沿三個(gè)坐標(biāo)軸 方向的節(jié)點(diǎn)數(shù)粒子密度球體模型圓柱體模型圓柱SphGenSWGenBlockMSurMeshnLMeshSmooth2D meshBMeshBulkFTetMeshBestFrtCurvesSurfaceChainMTTMeshMassTrDieLineHIP201SQualPeneCkHHSShowGMorph123456 7Dp Sph gene rati on interface-CModify方向坐標(biāo)底部中心坐標(biāo)PositionDirection 10 0 Pol

2、arV：0圓柱髙圓柱半徑NY1NX境徑向節(jié)點(diǎn)數(shù)，一 般設(shè)置相同F(xiàn)ilP/o 100% 100.0NZ 1軸向節(jié)點(diǎn)數(shù)圓臺(tái)/錐模型圓臺(tái)/錐底部中心坐標(biāo)圓臺(tái)底部半徑圓臺(tái)頂部半徑方向坐標(biāo)圓臺(tái)髙徑向節(jié)點(diǎn)數(shù)，一的諂罟相同軸向節(jié)點(diǎn)數(shù)數(shù)據(jù)庫(kù)Airbag Z*Dba5e*Mat才NodeBoundry*Elem*Pa ram*Cn5trnd*Eos41 Pa 比*ompnt*Hrglass*Rgdwal*ontactInitialSectionControl*Intgrtn丈Set*Def2Rg*Intrfac*TermntDamping*LadUser坷| Keyword Input FormuAccept

3、Delete DefaultDoneUse PARAHH(Subsys: 1)Setting*DATABASE_FORMAT (0)r Ke-y w o rcnD-AHABASEEdit曲何時(shí)些I： rroRNA?DoneModelCOMMENT:1 1 1 1 1 1 1 11 1 1 1 1 ! 24 Illiiiii -*i iiii1 i i 4口 ii IIII -*IQ I I I 1 t i i i i 11117iQ 111.| 8,4 $5 DATABASE60輸出結(jié)果格式2.參數(shù)定義*DATABASE_F()RMAT計(jì)算結(jié)果輸出格式Card12345678VanableIF

4、ORMIBINARYTypeIIDefault00Remarks12IFORM Output format for D3PLOT and D3THDT files:LS-DYNA database format (default),:ANSYS database format,:Both LS-DYNA and ANSYS database formats.IBINARY Word size of the binary output files (D3PLOT, D3THDT, D3DRLF and interface files for 64 bit computer such as CRA

5、Y and NEC:default 64 bit format,:32 bit IEEE format1 256IFOR 卜 I【BINARY4D數(shù)據(jù)庫(kù)* Airbag*MatAle羊DefineNode*BoundryElem* Param*Cnstrnd*Eos*Part*ompnt*HrglassRgdwal*ontact*Initial*Section* Control*Intgrtn*Set*Def2Rg*Intrfac*TermntDamping拿 Load*Userhf3)4 56 |打DQoneAllqt_SEON_SETModel輸出二進(jìn)制結(jié)果$DATABASE OPTION

6、S$ iT* J*of j*ofJ*T T T T*T f jsJ T of ? * ofJ* T *J*T T of T 7* :fJ j* M1* of T T Ta*a。a，a，a ra，a，丄.a，a，Tr a，a，a，丄 a，$ *DATABASE_EXTEWT_BIiJilRY0Q00Remarks：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

7、size of the binary output files which are created by 64 bits computer such as CRAY and NEC二進(jìn)制*DATABASE_EXTENT_BINARY 結(jié)果輸出控制_ Q11111111 IQ 11111 12A11111113Q i 111 14口 111111154iEJ1 丨i丨 7.q丨 i BI 94341543-4294343鉗344兩3439434SMW7Edit廠 Keyw oraDATABASEVanableNEIPHNHPSMAXDJTS.IBFLGSIGFLGEPSFLGRLITLGENGF

8、LGTypeIIIIIIIIDefauli0030-1111Rgmaifcs1For the BINARY option the foHoning card applyOCard 3 is oprionab；Card 2VariableCMPFLGEEVERPBEAX11PDCOMPSHGESTSSZN31HDTIALEMATTypeIIIIIIIIDefauLt00000021RenLaiis.2VaniabieNINTSLD PKP.SEKSC1PMSSCLTHERMTypeIIFIIDefault01.000RansiksNEIPH Number of additional integr

9、ation point history variables written to the binary data- base for solid elements. The integration point data is written in the same order that it is stored in memory-each material model has its own history variables that are stored For user defined materials it is important to store the history dat

10、a that is needed for plotting before the data which is not of interest給實(shí)體單元增加附加綜合的點(diǎn)歷史變量二進(jìn)制數(shù)據(jù)庫(kù)，綜合點(diǎn)數(shù)據(jù)目的是儲(chǔ)存記憶每個(gè) 材料模型使它擁有自己的歷史變量。對(duì)于使用者而言，定義材料對(duì)于儲(chǔ)存歷史數(shù)據(jù)是很 重要的，在數(shù)據(jù)儲(chǔ)存前需要對(duì)其進(jìn)行描繪。NEIPS Number of additional integration point history variables written to the binary data- base for both shell and thick shell elemen

11、ts for each integration point, see NEIPH above給薄殼和厚殼單元增加附加綜合的點(diǎn)歷史變量二進(jìn)制數(shù)據(jù)庫(kù)，看下面的|EIPHOMAXINT Number of shell integration points written to the binary database, see also *INTE GRATION_SHELL If the default value of 3 is used then results are outputfor the outermost (top) and innermost (bottom) integratio

12、n points together with results for the n eutral axis. If MAX I NT is set to 3 and the element has 1 integration point then all three results will be the same If a value other than 3 is used then results for the first MAXINT integration points in the element will be output Note： If the element has an

13、 even number of integration points and MAX I NT is not set to 3 then you will not get mid-surface results See Remarks below. _殼單元附加綜合的點(diǎn)歷史變量二進(jìn)制數(shù)據(jù)庫(kù)，可以參看卜IVTE GRAT 1 ON SHELL.如果采用默認(rèn)值3,則結(jié)果輸出的是最遠(yuǎn)(高)和最近(下)的綜合點(diǎn)，結(jié)果包含中間軸 坐標(biāo)。如果設(shè)置為3,單元有一個(gè)綜合點(diǎn)，則三個(gè)結(jié)果將相同。如果值超過(guò)3,則在單元 中的第一個(gè)MAXINT綜合點(diǎn)將會(huì)被輸出。注意：如果單元擁有一系列綜合點(diǎn)，而MAXINT沒(méi) 有被設(shè)

14、置成3,則你將無(wú)法獲得中間表面結(jié)果值，詳見(jiàn)F面的附注STRFLG Set to 1 to dump strain tensors for solid, shell and thick shell elements for plotting by LS-PREPOST and ASCII file ELOLT. For shell and thick shell elements two tensors are written, one at the innermost and one at the outermost integration point For solid elements a

15、 single strain tensor is written.SIGFLG Flag for including stress tensor in the shell LS-DYNA database：:include (default),:excludeEPSFLG Flag for including the effective plastic strains in the shell LS-DYNA database：:include (default)t:excludeRLTFLG Flag for including stress resultants in the shell

16、LS-DYNA database：:include (default),:excludeENGFLG Flag for including shell internal energy density and thickness in the LSDYNA database：:include (default),:excludeCMPFLG Orthotropic and anisotropic material stress and strain output in local material coordinate system for solids, shells and thick sh

17、ells:global,:localIEVERP Every plot state for ud3plotn database is written to a separate file This option will limit the database to 1000 states：:more than one state can be on each plotfile,:one state only on each plotfileBEAMIP Number of beam integration points for output This option does not apply

18、 to beams that use a resultant formulation.DCOMP Data compression to eliminate rigid body data：:off (default)T no rigid body data compression,:on, rigid body data compression active,:off, no rigid body data compression, but nodal velocities and accelerations are eliminated from the database:on, rigi

19、d body data compression active and nodal velocities and accelerations are eliminated from the databaseSHGE Output shell hourglass energy density：:off (default)t no hourglass energy written,:on.STSSZ Output shell element time step, mass, or added mass：:off (default),:output time step size,:output mas

20、s, added mass, or time step size See remark 3 belowN3THDT Material energy write option for D3THDT database:off, energy is NOT written to D3THDT database,:on (default), energy is written to D3THDT databaseIALEMAT Output solid part ID list containing ale materials:on (default)NINTSLD Number of solid e

21、lement integration points written to the LS-DYNA databas巳.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_

22、SEN Flag to output the peak pressure and surface energy computed by each contact in terface into the interface force database To obtain the surface energy, FRCENG, must be sent to 1 on the control contact card When PKP_SEN=1t it is possible to identify the energies generated on the upper and lower s

23、hell surfaces, which is important in metal forming appli- cations This data is mapped after each H-adaptive remeshing:No data is written:Output the peak pressures and surface energy by contact interfaceSCLP A sealing parameter used in the computation of the peak pressure This parameter is gen erally

24、 set to un ity (the default), but it must be greater tha n 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 incompatible with other 3rd party software:(default) output temperature:output tempe

25、rature:output temperature and flux:output temperature, flux, and shell bottom and top surface temperature Remarks：1. If MAX I NT is set to 3 the n mid-surface, inn er 一 surface and outer-surface stresses are output at the center of the element to the LS-DYNA database For an even number of integratio

26、n points, 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 computing the average of these points. For MAX I NT equal to 3 LS-DYNA assumes that the data for

27、 the user defined integration rules are ordered from bottom to top even if this is not the case. If MAX I NT is not equal to 3, the n 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 stre

28、sses are first averaged2. Beam stresses are output to the LSDYNA database if and only if BEAM IP is greater thanzero. 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 val

29、ues for elastic-plastic Hughes-Liu beams： the normal stress, FT； the transverse shear stresses, rs and tr； 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

30、the beam elements of Belytschko and his coworkers, the transverse shear stress components are not used in the formulation. No 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 gl

31、obal mass scaling is used for a constant time step, the total element mass is output； however, if the mass is increased so that a minimum time step size is maintained (DT2MS is negative), the added mass is output. Also, see the control card *CONTROL TIMESTEP.*DATABASE_BINARY_D3PL0T結(jié)果輸出時(shí)間步長(zhǎng)數(shù)據(jù)庫(kù)* Airba

32、g *Dtase*Mat*Ale*Define*hodeBoundryElem*Pa ram*Cnstrnd*EosPartCompnt*Hrglass* RgdwlContactInitial“SectionControl*Intgrtn*Set*Def2Rg*Intrfac*Termnt* Damping才 LoadUser12 3 | 567 DEditDone廝 Keyv/ord Input FormPick Accept Delete Defat0 Use *PARAMET(Sub,: 1)DATABASEBINARYDDPLOT ( 0 )ModelQ11111111 iq 111

33、111111111 i11! 11 i 111 i 1111111 1 1 1 11 11 11.1 ?.Q! . . 154:32 $rmE HISTORT$94：- $翱$翱$翱$翱$霸：$沁刃$*JA7AIASE BimY D3PLDT2.000Card1234567VaiiableDT/CYCLLCDT.NRBEAMNPLTCPSETIDTypeF1IIIDefault-RemarksOptional Card tlifit only applies to rlie D3PLOT dntabaseCaid1234567SVariableIOOPTTypeIDefault0Remarks

34、DT Time interval between outputsCYCL 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 n* th convergence check interval speci- fied on the *CONTROL_DYNAMIC_RELAXATION card.NR Number of Running Restart Files, RU

35、NRSF. written in a cyclical fashion. The default -KeywoixHBwABASEPART3PROP3PLOTnumber is one, the same file is overwritten each timeLCDT 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 *D

36、ATABASE_BINARY D3PLOT or D3PART.:Discrete spring and damper elernents are added to the D3PLOT or D3PART database where they are display as beam elements The elernent global X, global Y, global Z and resultant forces are written to the database,:No discrete spring and damper elements are added to the

37、 D3PLOT or D3PART databas巳.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 exist- ing beam ele

38、ments causing 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：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 t

39、his 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 post-processing applicationsNPLTC DT=ENDTIME/NPLTC applies to D3PLOT and D3PART only. This overrides t

40、he DT specified in the first fieldPSETID SET_PART ID for D3PART only.IOOPT This option applies to the D3PLOT file only. Flag to govern behavior of the plot fre- quency load curve defined by LCDT：:At the time each plot is generated, the load curve value is added to the current time to determine the n

41、ext plot time (this is the default behavior):At the time each plot is generated, the next plot time T is computed so that T = the current time plus the load curve value at time T.:A plot is generated for each abscissa point in the load curve def initio n. The actual value of the load curve is ignore

42、d*C0NTR0L SPH畐 Keywd Input Form4JearDefaultAcceptDeleteDoneUsexPARAMHERCOffTROL.SPH( (0) )&bsys 1)SettingA1 HCBS Bono H roni HEMORY FORM SJABI HAXVn PP_ D ECOM POSITION_SCAL E_ n PP_ D ECOM POSmON_SCAL E M PP_ D ECOM POSmON_SHOW M PP_ D ECOni POSITIONS RANS riPP_IO_DINOUTONLY M PP_IO_LST 0 0巾IJSRWRU

43、SRFRCnntRNUKKPTTRDPR)000o 4.00o |o5ERKDffiKEDCYKIHTH SFPHLX0.00.00.00.00.00.00.0IGIK)RfRCEMGSKIPRWGOUTSLGSPOTSTPSPOTDELSPOTHINi 111111111 11 1 1 1 1111 -*iQ 1111111 111111i11111 1111 i i i iQ 111iiii &Qi COBTROL COUTACT1.000000 0.00010001000004.000000000100.0000.0000.0000.0000.0000.0000.000110000000

44、.000120000.0001 00000000.0000id 112345678VariableSLSFACRWPNALISLCHKSHLTHKPENOPTTHKCHGORDENENMASSTypeFFIIIIIIDefault.1none101010Card?VariableUSRSTRUSRFRCNSBCSIN 丁 ERMXPENESSTHKECDTTIEDPRJTypeIIIIFIIIDefault0010-10004.0000Card 3 is optional. The following parameters are the default values used by part

45、spKeywordXONTROL0 0bO 0* 05ISYHNSEROD RWGAP5 RWGDTH RWKST ICOVSWRADF ITHOfFo討 ovjo.o 1.00 0.006SHLEDGinautoma tic contacts. These frictional coefficie nts apply only to contact types： SINGLE_SURFACE, AUTOMATIC.GENERAL, AUTOMATIC_SINGLE_ SURFACE,AUTOMATIC NODES TO , AUTOMATIC SURFACE , and AUTOMATICO

46、NE_WAY_. , and ERODING_SINGLE_SURFACE. Also see *CONTACT and * that these default values will override the values specified for these con tact types in the*CONTACT section.Card 31234567VanableSFRICDFRICEDCVFCTHTH_SFPEN_SFTypeFFFrFFrDefault000.000000.00000Card 4 is optional. If this card is defined,

47、then C sud 3 nbove must be included A blank card mav be inserted for Card 3.Sd 412345678VariableIGNOREFRCENGSKIPR.WGOUTSEGSPOTSTPSPOTDELSPOTHINTypeIIIIIIFDefault000000lnactiv-eCard 5 is optional. If this card is defined, tlieu Cards 3 aud 4 above must be included. Blank cards mav be inserted.Card 51

48、2345678Variable1SYMNSERODRWGAPSRWGDTHRXVKSFICOVSWR.4DFITHOFFTypeIIIFFIFIDcfoult0000.1.000.0Card 6 k oprionnl. If this card is defined then Carck 3 TO 5 above must be included. Blank cards may be inserted.Card 512345678VariableSHLEDGTypeIDefault0SLSFAC Scale factor for siiding interface penalties, SL

49、SFAC：:default = 1.RWPNALScale factor for rigid wall penalties, which treat nodal points interacting with rigid walls, RWPNAL. The penalties are set so that an absolute value of unity should be optimal: however, this penalty value may be very problem dependent If rigid/deformable materials switching

50、is used, this option should be used if the switched materials are interacting with rigid walls all nodes are treated by the penalty method This is required for implicit calcula- tions Since seven (7) variables are stored for each slave node, only the nodes that may interact with the wall should be i

51、ncluded in the node listthe constraint method is used and nodal points which belong to rigid bodies are not considered rigid bodies nodes are treated by the penalty method and all other nodes are treated by the constraint methodISLCHK Initial penetration check in contact surfaces with indication of

52、initial penetration in output files (see remarks below):the default is set to 1,:no checking,:full check of initial penetration is performedSHLTHK Shell thickness considered in type surface to surface and node to surface type contact options, where options 1 and 2 below activate the new contact algo

53、rithms The thick- ness offsets are always ineluded in single surface, constraint method, and automatic surface to surface and node to surface contact types (See remarks below)：:thickness is not considered,:thickness is considered but rigid bodies are excluded,:thickness is considered including rigid

54、 bodiesPENOPT Penalty stiffness value opt ion. For default calculation of the penalty value please refer to the LS-DYNA Theory Manua1:the default is set to 1,:minimum of master segme nt and slave node (default for most con tact types),:use master segment stiffness (old way),:use slave node value,:us

55、e slave node value, area or mass weighted,:same as 4 but inversely proportional to the shell may require special scaling and is not generally recommended Options 4 and 5 can be used for metalforming calculationsTHKCHG Shell thickness changes considered in single surface contact：:no consideration (de

56、fault),:shell thickness changes are ineludedORIEN Optional automatic reorientation of contact interface segments during initialization：:default is set to 1.:active for automated (part) input only. Con tact surfaces are give n by *PART definitions:active for manual (segment) and automated (part) inpu

57、t:inactive for non-forming contact:inactive for forming contactENMASS Treatme nt of the mass of eroded n odes in con tact This option affects all contact types where nodes are removed after surrounding elements fai 1. Generally, the removal of eroded nodes makes the calculation more stable； however,

58、 in problems where erosion is impor- tant the reduction of mass will lead to incorrect results:eroding nodes are removed from the calculation.:eroding nodes of solid elernents are retained and continue to be active in contact:the eroding nodes of solid and shell elements are retained and continue to

59、 be active in contactUSRSTR Storage per contact interface for user supplied interface control subroutine, see Appe ndix F If zero, no input data is read and no in terface storage is permitted in the user subroutine This storage should be large enough to accommodate input parameters and any history i

60、n put data is available in the user supplied subroutineUSRFRC Storage per contact interface for user supplied interface friction subroutine, see Appe ndix G If zero, no in put data is read and no interface storage is permitted in the user subroutine This storage should be large enough to accommodate