在LS-Prepost 建立sph模型的详细过程

LS_Prepost

关键操作

1.

模型建立

立方体模型

起始坐标

立方体

Part

编号

沿三个坐标轴

方向的节点数

起始节点号：

该

Part

第一个节点号码

粒子密度

1

球体模型

球心坐标

球体

沿三个坐标轴

方向的节点数

球体半径

2

圆柱体模型

圆柱

方向坐标

底部中心

坐标

圆柱高

圆柱半径

径向节点数，

一

般设置相同

轴向节点数

3

圆台

/

锥模型

圆台

/

锥

方向坐标

底部中心

坐标

圆台高

圆台底部半径

径向节点数，

一

般设置相同

圆台顶部半径

轴向节点数

4

2.

参数定义

*DATABASE_FO RMAT

计算结果输出格式

数据库

输出结

果格式

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 computer such as CRAY and NEC.

EQ.0: default 64 bit format,

5

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 createdby

64 bits computer such as CRAY and NEC

.

二进制

*DATABASE_EXTENT_BINARY

结果输出控制

数据库

输出二进

制结果

6

NEIPH

Number of additional integration point history variables written to thebinarydata- base

for

solid

elements.

The

integration

point

data

is

writtenin

the

same

order

that

it

is

stored

in

memory-each

material

model

has

itsown

history

variables

that

are

stored.

For

user

defined

materials it isimportant to store the history data that is needed for plotting before thedatawhich

is not of interest.

给实体单元增加附加综合的点历史变量二进制数据库，综合点数据目的是储存记忆每个

材料模型使它拥有自己的历史变量。对于使用者而言，定义材料对于储存历史数据是很

< br>重要的，在数据储存前需要对其进行描绘。

NEIPS

Number of additional integration point history variables written to thebinarydata- base

for both shell and thick shell elements for eachintegration point, see NEIPH above.

给薄壳和厚壳单元增加附加综合的点历史变量二进制数据库，看下面的

NEI PH

。

MAXINT

Number

of

shell

integration

points

written

to

the

binary

database,

seealso

*INTEGRATION_SHELL.

If

the

default

value

of

3

is

used

thenresults

are

output

for

the

outermost (top) and innermost (bottom)integration points together with results for the neutral

axis. If MAXINTis set to 3 and the element has 1 integration point then all three resultswill be

7

the same. If a value other than 3 is used then results for the firstMAXINT integration points in

the element will be output. Note: If theelement has an even number of integration points and

MAXINT is notset to 3 then you will not get mid-surface results. See Remarks below.

壳单元附加综合的点历史变量二进制数据库，可以参看

*INTEGRATI ON_SHELL.

如果采用默认值

3

，则结果输出的是最远（高）和最近（下）的综合点，结果包含中间轴

坐标。如果设置 为

3

，单元有一个综合点，则三个结果将相同。如果值超过

3

，则在单元

中的第一个

< br>MAXINT

综合点将会被输出。

注意：

如果单元拥有一系列综合点，

而

MAXINT

没有被设置成

3

，则你将无法获得中间表面结果 值，详见下面的附注。

STRFLG

Set to 1 to dump strain tensors for solid, shell and thick shell elementsfor plotting by

LS-PREPOST

and

ASCII

file

ELOUT.

For

shell

andthick

shell

elements

two

tensors

are

written, one at the innermost andone at the outermost integration point. For solid elements a

single straintensor 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-DYNAdatabase:

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

LSDYNAdatabase:

EQ.1: include (default),

EQ.2: exclude.

CMPFLG

Orthotropic and anisotropic material stress and strain output in localmaterial

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. Thisoption 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 notapply to beams

that use a resultant formulation.

DCOMP

Data compression to eliminate rigid body data:

EQ.1: off (default), no rigid body data compression,

EQ.2: on, rigid body data compression active,

EQ.3: off, no rigid body data compression, but nodal velocities andaccelerations are

eliminated from the database.

EQ.4: on, rigid body data compression active and nodal velocitiesand accelerations are

eliminated from the database.

8

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 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-DYNAdatabase. The

default

value

is

1.

For

solids

with

multiple

integrationpoints

NINTSLD

may

be

set

to

tly,

no

other

values

forNINTSLD

are

allowed.

For

solids

with

multiple

integration

points, anaverage value is output if NINTSLD is set to 1.

PKP_SEN

Flag

to

output

the

peak

pressure

and

surface

energy

computed

by

eachcontact

interface into the interface force database. To obtain the surfaceenergy, FRCENG, must be sent

to

1

on

the

control

contact

card.

WhenPKP_SEN=1,

it

is

possible

to

identify

the

energies

generated

on

theupper

and

lower

shell

surfaces,

which

is

important

in

metal

formingappli-

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 contactinterface

SCLP

A

scaling

parameter

used

in

the

computation

of

the

peak

pressure.

Thisparameter

is

generally set to unity (the default), but it must be greaterthan 0.

MSSCL

Output

nodal

information

related

to

mass

scaling

into

the

D3PLOTdatabase.

This

option can be activated if and only if DT2MS < 0.0, seecontrol card *CONTROL_TIMESTEP.

This option is available startingwith 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 incompatible with other 3

rd

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 areoutput at

the center of the element to the LS-DYNA database. For an even number ofintegration points,

the

points

closest

to

the

center

are

averaged

to

obtain

the

midsurfacevalues.

If

multiple

integration points are used in the shell plane, the stresses at the centerof the element are found

by computing the average of these points. For MAXINT equalto 3 LS-DYNA assumes that the

9

data for the user defined integration rules are orderedfrom bottom to top even if this is not the

case. If MAXINT is not equal to 3, then thestresses at the center of the element are output in

the order that they are stored for theselected integration rule. If multiple points are used inplane

the stresses are firstaveraged.

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 theintegration

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 transverseshear stresses,

?

rs

and

?

tr

; the effective plastic strain, and the axial strain which islogarithmic. For beams that are not

elastic-plastic,

the

first

history

variable,

if

any,

isoutput

instead

of

the

plastic

strain.

For

the

beam elements of Belytschko and his coworkers,the transverse shear stress components are not

used in the formulation. No datais output for the Belytschko-Schwer resultant beam.

3. If mass scaling is active, the output of the time step size reveals little information aboutthe

calculation.

If

global

mass

scaling

is

used

for

a

constant

time

step,

the

total

elementmass

is

output;

however,

if

the

mass

is

increased

so

that

a

minimum

time

step

size

ismaintained

(DT2MS

is

negative),

the

added

mass

is

output.

Also,

see

the

control

card*CONTROL_TIMESTEP.

*DATABASE_BINARY_D3PLOT

结果输出时间步长< /p>

数据库

输出时间

步长

10

DT

Time interval between outputs.

CYCL

Output interval in time steps (a time step is a cycle). For the D3DRFLfile a positive

number 'n' will cause plot dumps to be written at everyn'th convergence check intervalspeci-

fied on the *CONTROL_DYNAMIC_RELAXATION card.

NR

Number of Running Restart Files, RUNRSF, written in a cyclicalfashion. The default

number is one, i.e. the same file is overwritten eachtime.

LCDT

Optional load curve ID specifying time interval between dumps. Thisoption is only

available for the D3PLOT, D3PART, D3THDT andINTFOR files.

BEAM

Option flag for *DATABASE_BINARY_D3PLOT or D3PART.

EQ.0: Discrete spring and damper elements are added to theD3PLOT or D3PART

database where they are display as beamelements. The element global X, global Y, global Z

and resultantforces are written to the database,

EQ.1: No discrete spring and damper elements are added to theD3PLOT or D3PART

database. This option is useful whentranslating old LS-DYNA input decks to KEYWORD

input. Inolder input decks there is no requirement that beam and springelements have unique

ID's, and beam elements may be created forthe spring and dampers with identical ID's to exist-

ing beamelements causing a fatal error. Contact interfaces which are basedon part IDs ofseat-

belt elements will not be properly generated ifthis option is .2: Discrete spring and

damper elements are added to theD3PLOT or D3PART database where they are displayed as

beamelements (similar to option 0). In this option the element resultantforce is written to its

first database position allowing beam axialforces and spring resultant forces to be plotted at the

same can be useful during some post- processing applications.

NPLTC

DT=ENDTIME/NPLTC applies to D3PLOT and D3PART only. Thisoverrides the

11

DT specified in the first field.

PSETID

SET_PART ID for D3PART only.

IOOPT

This option applies to the D3PLOT file only. Flag to govern behavior ofthe plot fre-

quency load curve defined by LCDT:

EQ.1: At the time each plot is generated, the load curve value isadded to the current time

to determine the next plot time (this is thedefault behavior).

EQ.2: At the time each plot is generated, the next plot time T iscomputed so that T = the

current time plus the load curve value attime T.

EQ.3: A plot is generated for each abscissa point in the load curvedefinition. The actual

value of the load curve is ignored

*CONTROL_SPH

SPH

12