在ls-prepost 建立sph模型的详细过程.docx

LS_Prepost关键操作
1. 模型建立
立方体模型
沿三个坐标轴方向得节点数
立方体
起始节点号:该Part第一个节点号码
粒子密度
Part编号
起始坐标

球体模型
球体半径
球心坐标
球体
沿三个坐标轴方向得节点数

圆柱体模型
轴向节点数
径向节点数,一般设置相同
圆柱高
圆柱半径
底部中心坐标
方向坐标
圆柱

圆台/锥模型
圆台/锥
圆台顶部半径
圆台底部半径
轴向节点数
径向节点数,一般设置相同
圆台高
底部中心坐标
方向坐标

2. 参数定义
*DATABASE_FORMAT计算结果输出格式
数据库

输出结果格式

IFORM Output format for D3PLOT and D3THDT files
EQ、0: LS-DYNA database format (default),
EQ、1: ANSYS database format,
EQ、2: Both LS-DYNA and ANSYS database formats、
IBINARY Word size of the binary output files (D3PLOT, D3THDT, D3DRLF and
interface files for 64 bit puter such as CRAY and NEC、
EQ、0: default 64 bit format,
EQ、1: 32 bit IEEE format
Remarks:
1、 This option is not available for every platform、 Check LS-DYNA Banner upon executionof the program
2、 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、二进制
*DATABASE_EXTENT_BINARY结果输出控制
数据库

输出二进制结果

NEIPH Number of additional integration 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 data that is needed for plotting before the data which is not of interest、
给实体单元增加附加综合得点历史变量二进制数据库,综合点数据目得就是储存记忆每个材料模型使它拥有自己得历史变量。对于使用者而言,定义材料对于储存历史数据就是很重要得,在数据储存前需要对其进行描绘。
NEIPS Number of additional integration point history variables written to the binary data- base for both shell and thick shell elements for each integration point, see NEIPH above、
给薄壳与厚壳单元增加附加综合得点历史变量二进制数据库,瞧下面得NEIPH。
MAXINT 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 output for the outermost (top) and innermost (bottom) integration points together with results for the neutral axis、 If MAXINT 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 even number of integration points and MAXINT is not set to 3 then you will not get mid-surface results、 See Remarks below、
壳单元附加综合得点历史变量二进制数据库,可以参瞧*INTE GRATION_SHELL、

如果采用默认值3,则结果输出得就是最远(高)与最近(下)得综合点,结果包含中间轴坐标。如果设置为3,单元有一个综合点,则三个结果将相同。如果值超过3,则在单元中得第一个MAXINT综合点将会被输出。注意:如果单元拥有一系列综合点,而MAXINT没有被设置成3,则您将无法获得中间表面结果值,详见下面得附注。
STRFLG Set to 1 to dump strain tensors for solid, shell and thick shell elements for plotting by LS-PREPOST and ASCII 、 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 single strain tensor is written、
SIGFLG Flag for including stress tensor in the shell LS-DYNA database:
EQ、1: include (default),
EQ、2: exclude、
EPSFLG Flag for including the effective plastic strains in the shell LS-DYNA database:
EQ、1: include (default),
EQ、2: exclude、
RLTFLG Flag for including stress resultants in the shell LS-DYNA database:
EQ、1: include (default),
EQ、2: exclude、
ENGFLG Flag for including shell internal energy density and thickness in the LSDYNA database:
EQ、1: include (default),
EQ、2: exclude、
CMPFLG Orthotropic and anisotropic material stress and strain output in local material coordinate system for solids, shells and thick shells、
EQ、0: global,
EQ、1: local、
IEVERP Every plot state for “d3plot” database is written to a separate file、 This option will limit the database to 1000 states:
EQ、0: more than one state can be on each plotfile,
EQ、1: one state only on each plotfile、
BEAMIP Number of beam integration points for output、 This option does not apply to beams that use a resultant formulation、
DP Data pression to eliminate rigid body data:
EQ、1: 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 are eliminated from the database、
SHGE Output shell hourglass energy density:
EQ、1: off (default), no hourglass energy written,
EQ、2: on、
STSSZ Output shell element time step, mass, or added mass:
EQ、1: off (default),
EQ、2: output time step size,
EQ、3: output mass, added mass, or time step size、 See remark 3 below、
N3THDT Material energy write option for D3THDT database
EQ、1: off, energy is NOT written to D3THDT database,
EQ、2: on (default), energy is written to D3THDT database、
IALEMAT Output solid part ID list containing ale materials、
EQ、1: on (default)
NINTSLD Number of solid element integration points written 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 peak pressure and surface energy puted 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 important in metal forming appli- cations、 This data is mapped after each H-adaptive remeshing、
EQ、0: No data is written
EQ、1: Output the peak pressures and surface energy by contact interface
SCLP A scaling parameter used in the putation of the peak pressure、 This parameter is generally set to unity (the default), 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、0, see control card *CONTROL_TIMESTEP、 This option is available starting with the second release of Version 971、
EQ、0: No data is written
EQ、1: Output incremental nodal mass
EQ、2: Output percentage increase in nodal mass
THERM Output of thermal data to d3plot、 The use of this option (THERM>0) may make the database inpatible with other 3rd party software、
EQ、0: (default) output temperature
EQ、1: output temperature
EQ、2: output temperature and flux
EQ、3: output temperature, flux, and shell bottom and top surface temperature
Remarks:
1、 If MAXINT is set to 3 then mid-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, 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 puting the average of these points、 For MAXINT equal to 3 LS-DYNA assumes that the data for the user defined 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 averaged、
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-plastic Hughes-Liu beams: the normal stress, órr; 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 the beam elements of Belytschko and his coworkers, the transverse shear stress ponents 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 global 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_D3PLOT结果输出时间步长
输出时间步长

DT Time interval between outputs、
CYCL 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, 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_D3PLOT or D3PART、
EQ、0: Discrete spring and damper 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、1: No discrete spring and damper elements are added to the D3PLOT or D3PART database、 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 ID's, and beam elements may be created for the spring and dampers with identical ID's to exist- ing beam elements 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、 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 this 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 applications、
NPLTC DT=ENDTIME/NPLTC applies to D3PLOT and D3PART only、 This overrides the DT specified in the first field、
PSETID SET_PART ID for D3PART only、
IOOPT This option applies to the D3PLOT 、 Flag to govern behavior of the plot fre- quency load curve defined by LCDT:
EQ、1: At the time each plot is generated, the load curve value is added to the current time to determine the next plot 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
*CONTROL_SPH

SPH

NCBS Number of cycles between particle sorting 粒子分类搜索循环次数
BOXID SPH approximations are puted inside a specified BOX、 When a particle has gone outside the BOX, it is deactivated、 This will save putational time by eliminating particles that no longer interact with the structure、
指定BOX内得SPH粒子参与计算。当某个SPH粒子位于BOX之外时,该粒子失效,通过消除某些不再与结构发生作用得粒子,可以节省计算时间。
DT Death time、 Determines when the SPH calculations are stopped、 粒子失效时间
IDIM Space dimension for SPH particles: SPH粒子得空间维数
3 for 3D problems 参数为3时,表示3维问题。
2 for 2D plane strain problems 参数为2时,表示2维问题。
-2 for 2D axisymmetric problems参数为-2时,表示轴对称问题。
When a value is not specified LS-DYNA determines the space dimension automatically by checking the use of 3D, 2D or 2D asisymmetric elements、
当这个值无法自动指定LS-DYNA得维数,程序通过核对使用3维、2维或轴对称得单元来确定空间维数。
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 variable、 Default value should apply for most applications、
定义每个粒子得初始相邻粒子得数量,该变量只就是在初始化阶段调整内存分配。在计算中,如果某些粒子需要更多得相邻粒子,LS-DYNA会自动调整该变量设置。默认值适用于大部分问题。
FORM Particle approximation theory: 粒子近似理论
EQ、 0: default formulation, 默认公式
EQ、 1: remormalization approximation 重归-近似化
START Start time for particle approximation、 Particle approximations will be puted when time of the analysis has reached the value defined in START、 粒子近似开始时间。当分析时间达到所设定得值时,粒子将开始计算。
MAXV Maximum value for velocity for the SPH particles、 Particles with a velocity greater than MAXV are deactivated SPH粒子速度最大值,如果速度超过该值,质点将会失效。
CONT Defines the putation of the particle approximation between two different SPH parts: 两个不同得SPH组之间粒子近似得计算设置。
EQ、 0: Particle approximation is defined (default)参数设为0 计算粒子近似(默认值)
EQ、 1: Particle approximation is not puted、 Two different SPH materials will not interact with each other and penetration is allowed、 参数设为1 不计算粒子近似,两种不同得SPH材料不会发生相互作用,允许相互穿透。
DERIV Time integration type for the smoothing length:光滑长度得时间积分类型。

*CONTROL_CONTACT
CONTACT

Card 3 is optional、 The following parameters are the default values used by parts in
automatic contacts、 These frictional coefficients apply only to contact types:
SINGLE_SURFACE, AUTOMATIC_GENERAL, AUTOMATIC_SINGLE_ SURFACE,
AUTOMATIC_NODES_TO_、、、, AUTOMATIC_SURFACE_、、、, and AUTOMATIC_
ONE_WAY_、、、、, and ERODING_SINGLE_SURFACE、 Also see *CONTACT and *PART、Note that these default values will override the values specified for these contact types in the*CONTACT section、

SLSFAC Scale factor for sliding interface penalties, SLSFAC:
EQ、0: default = 、1、
RWPNAL Scale 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 is used, this option should be used if the switched materials are interacting with rigid walls、 LT、0、0: 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 included in the node list、
EQ、0、0: the constraint method is used and nodal points which belong to rigid bodies are not considered、 GT、0、0: rigid bodies nodes are treated by the penalty method and all other nodes are treated by the constraint method、
ISLCHK Initial penetration check in contact surfaces with indication of initial penetration in output files (see remarks below):
EQ、0: the default is set to 1,
EQ、1: no checking,
EQ、2: full check of initial penetration is performed、
SHLTHK 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 algorithms、 The thick- ness offsets are always included in single surface, constraint method, and automat