LSDYNA常见问题集锦文档格式.docx

资源描述

LSDYNA常见问题集锦文档格式.docx

《LSDYNA常见问题集锦文档格式.docx》由会员分享，可在线阅读，更多相关《LSDYNA常见问题集锦文档格式.docx（35页珍藏版）》请在冰豆网上搜索。

LSDYNA常见问题集锦文档格式.docx

5LS-DYNA求解有时为什么有负的滑动能

这是由于在建立模型时PART与PART之间有初始穿透，尤其是壳单元模型时很容易发生，应当防止这种情况的出现，否那么容易在有初始穿透的地方产生塑性铰，原因是程序在求解的开场阶段给与穿透相应的接触力消除穿透，使材料发生局部塑性变形。

解决方法见2002年11月的应用技巧。

6在DYNA中如何考虑材料失效

在LS-DYNA的材料库中，能考虑失效的材料其失效模式往往比较单一，或者是应力失效，或是应变失效，如果材料本身较为复杂，在破坏过程中可能涉及多种失效模式，能否在一种材料中同时定义多种实效模式？

答：

可以。

LS-DYNA材料库中提供了专门定义失效准那么的命令，即*mat_add_erosion，利用该命令，可以同时定义压力、主应力、等效应力、主应变、临界应力以及应力脉冲六种失效准那么，在加载过程中满足任何一种失效准那么都会使材料发生破坏。

7在LS-DYNA中能否施加跟随力和跟随力矩？

能，对于一些应用，施加的载荷相对与坐标系不仅大小变化，而且方向变化，此时按照通常的施加力方向〔X、Y、Z〕不能满足要求，在LS-DYNA中，可以方便的施加跟随力和跟随力矩，在关键字*LOAD_NODE_OPTION中，对DOF选择4和8就可以施加跟随力和跟随力矩。

8如果在工程上遇到壳的厚度是坐标位置的函数时，这样的壳单元模型如何建立？

我们常用到等厚度的壳单元，如果在工程上遇到壳的厚度是坐标位置的函数时，这样的壳单元模型如何建立？

要用到RTHICK命令。

/PREP7

ET,1,63

RECT,,10,,10

ESHAPE,2

ESIZE,,20

AMESH,1

EPLO

MXNODE=NDINQR（0,14）

*DIM,THICK,,MXNODE

*DO,NODE,1,MXNODE

*IF,NSEL（NODE）,EQ,1,THEN

THICK（node）=0.5+0.2*NX（NODE）+0.02*NY（NODE）**2

*ENDIF

*ENDDO

NODE=

MXNODE=

RTHICK,THICK

（1）,1,2,3,4

/ESHAPE,1.0

/USER,1

/DIST,1,7

/VIEW,1,-0.75,-0.28,0.6

/ANG,1,-1

/FOC,1,5.3,5.3,0.27

9不同自由度的单元在merge时，ANSYS如何处理自由度的不匹配问题

ANSYS可以将在给定容差〔tolerance〕内的节点通过merge可以合并到一起，它可以是不同的单元类型和不同的自由度〔如：

壳或梁〔6自由度〕与块单元〔3自由度〕，ANSYS在处理共节点的自由度关系使，类似于藕合自由度〔CoupleDOFs〕而非约束方程，即只把一样的自由度联系起来，这样就可以方便的处理如柱销、铰链之类的约束问题。

10如何方便地建立压力容器的有限元模型?

在建立压力容器模型时，充分考虑模型的对称性，灵活利用ANSYS软件的工作平面和坐标系，利用合理的分网工具，可以得到漂亮的有限元模型，如下面的命令流所示：

/prep7

et,1,45

mp,ex,1,2e11

mp,nuxy,1,0.3

cylind,0.89,0.8,0,-1.7125,90,270,

wpoff,0,0,-1.7125

wprot,0,90,

cylind,0.47/2,0.37/2,-1.5,0,90,180,

vovlap,all

vsel,s,loc,x,-0.11,0

vdel,all,,,1

vsel,s,loc,y,0.3,0.5

vsel,s,loc,y,-0.3,0.3

vsbw,all

afillt,21,11,0.1

al,33,50,5

al,37,53,7

alls,all

va,15,13,25,24,11

kgen,2,35,,,,-0.2,,,0

lstr,35,15

adrag,54,,,,,,21

vsba,2,16

lsel,s,loc,y,1.5

lsel,u,length,,0,0.06

lesize,all,,,16,

vsel,s,loc,y,0.9,1.5

vsweep,all

vsel,s,loc,y,0.89,0.99

vsel,s,loc,y,0.8,0.89

asel,s,loc,z,

asel,r,loc,y,0,1

asel,a,loc,y,0

accat,all

vsel,s,loc,y,0,0.6

vmesh,all

vsel,s,loc,y,0,-0.6

WPCSYS,-1,0

CSWPLA,11,0,1,1,

VSYMM,z,all,

nummrg,all

numcmp,all

CSWPLA,12,0,1,1,

wpoff,0,0,-3.425

wprot,0,-90

sphere,0.89,0.8,90,180,

csys,0

vsel,s,loc,z,-8.0,-6.85

vsweep,all

11ANSYS是否具有混合分网的功能？

ANSYS具有混合网格剖分的功能。

例如两个粘在一起的面，可以对一个面进展三角形划分，再对另一个面进展四边形划分。

过程见以下命令：

et,1,42

rect,,1,,1

rect,1,2,,1

aglue,all

mshape,0,2d

amesh,1

mshape,1,2d

amesh,3

12在ANSYS中怎样给面施加一个非零的法向位移约束?

给面施加非零的法向位移约束的过程如下:

1）在面上施加一个对称约束条件（DA,2,SYMM）

2）将实体模型上的载荷传递到有限元模型（SBCTRAN）

3）选择需要施加约束条件的面（ASEL,S,,,2）

4）选择附在面上的节点（NSLA,S,1）

5）创立节点组元（CM,AREA2_N,NODE）

6）删除面上的对称约束条件（DADELE,2,SYMM）

7）删除节点上的零位移约束（DDELE,AREA2_N,UY）

8）在节点组元上施加一个非零的法向位移约束（D,AREA2_N,UY,.05）

9）图示节点验证约束是否正确（NPLOT）

13如何得到径向和周向的计算结果?

在圆周对称构造中，如圆环构造承受圆周均布压力。

要得到周向及径向位移，可在后处理/POST1中，通过菜单GeneralPostproc>

OptionsforOutp>

Rsys>

Globalcylindric或命令Rsys,1将结果坐标系转为极坐标，那么X方向位移即为径向位移，Y向位移即为周向位移。

14如何正确理解ANSYS的节点坐标系

节点坐标系用以确定节点的每个自由度的方向，每个节点都有其自己的坐标系，在缺省状态下，不管用户在什么坐标系下建立的有限元模型，节点坐标系都是与总体笛卡尔坐标系平行。

有限元分析中的很多相关量都是在节点坐标系下解释的，这些量包括：

输入数据：

1自由度常数

2力

3主自由度

4耦合节点

5约束方程等

输出数据：

1节点自由度结果

2节点载荷

3反作用载荷等

但实际情况是，在很多分析中，自由度的方向并不总是与总体笛卡尔坐标系平行，比方有时需要用柱坐标系、有时需要用球坐标系等等，这些情况下，可以利用ANSYS的“旋转节点坐标系〞的功能来实现节点坐标系的变化，使其变换到我们需要的坐标系下。

具体操作可参见ANSYS联机帮助手册中的“分析过程指导手册->

建模与分网指南->

坐标系->

节点坐标系〞中说明的步骤实现。

15如何考虑构造分析中的重力

在构造分析中，如何模拟构造自重和设备重量是一个经常遇到的问题，对于构造自重有两点要注意：

1．在材料性质中输入密度，如果不输入密度，那么将不会产生重力效果。

2．因为ANSYS将重力以惯性力的方式施加，所以在输入加速度时，其方向应与实际的方向相反。

对于构造上的设备重量可以用MASS21单元来模拟，该单元为一个空间“点〞单元。

设备重量可通过单元实常数来输入。

下面附上一个小例子〔设重力方向向下〕。

et,2,21

r,2,10,10,10

mp,ex,1,2e5

mp,dens,1,1

rect,,10,,1

esize,.5

amesh,all

type,2

real,2

e,node（5,1,0）

fini

/solu

dk,1,all

dk,2,uy,

acel,,10

solve

/post1

plnsol,u,sum,2

/SOLU

ANTYPE,MODAL

MODOPT,subspa,2,,,2,ON

MXPAND,2,,,YES

SOLVE

FINISH

/post1

set,1,1

etabl,kene,kene

ssum

*get,keneval1,ssum,,item,kene

*get,freqval1,mode,1,freq

eigen1=（2*3.14159*freqval1）**2

pmass1=2*keneval1/eigen1

set,1,2

ssum

*get,keneval2,ssum,,item,kene

*get,freqval2,mode,2,freq

eigen2=（2*3.14159*freqval2）**2

pmass2=2*keneval2/eigen2

finish

16对于具有高度不规那么横截面的3D模型什么是最正确网格划分方法？

在横截面上自由划分四边形网格，然后在体内扫掠成六面体单元。

在扫掠前可对四边形网格加密〔如需要〕。

确认加密后生产的单元保持四边形以保证扫掠成六面体单元。

〔ANSYS5.6/FLOTRAN支持锲形单元，所以无此要求。

〕

17在交互方式下如何施加任意矢量方向的外表载荷?

假设需在实体外表上施加任意方向的外表载荷，可通过在实体外表生成外表效应单元〔比方SURF154单元〕的方法来完成。

施加面载荷时，可施加在外表效应单元上，这样可以任意控制面力的方向。

加载过程中，选定外表效应单元，对话框中LKEY取值不同，那么所加外表载荷的方向不同。

〔请仔细看一看surf154的单元手册〕。

比方：

LKEY=1〔缺省〕，载荷垂直于外表；

LKEY=2，载荷为+X切向；

LKEY=3，载荷为+Y切向；

LKEY=4，载荷垂直于外表；

LKEY=5，那么可输入任意矢量方向的载荷。

特别地：

LKEY=5，VALUE项为均布压力值

VAL2、VAL3、VAL4三项的值确定矢量的方向。

18LS-DYNA94版后〔95和96〕在爆炸及流固耦合方面的功能增强

在LS-DYNA中，处理爆炸和流固耦合单元一般采用ALE列式和Euler列式〔也可采用Lagrange〕，从而克制单元严重畸变引起的数值计算困难，并实现流体-固体耦合的动态分析。

ALE列式先执行一个或几个Lagrange时步计算，此时单元网格随材料流动而产生变形，然后执行ALE时步计算：

〔1〕保持变形后的物体边界条件，对内部单元进展重分网格，网格的拓扑关系保持不变，称为SmoothStep；

〔2〕将变形网格中的单元变量〔密度、能量、应力张量等〕和节点速度矢量输运到重分后的新网格中，称为AdvectionStep。

用户可以选择ALE时步的开场和终止时间，以及其频率。

Euler列式那么是材料在一个固定的网格中流动，在LS-DYNA中只要将有关实体单元标志Euler算法，并选择输运（advection）算法。

LS-DYNA还可将Euler网格与全Lagrange有限元网格方便地耦合，以处理流体与构造在各种复杂载荷条件下的相互作用问题，并在95和96版中得到了极大的增强。

19ANSYS坐标系总结

工作平面〔WorkingPlane〕

工作平面是创立几何模型的参考（X,Y）平面，在前处理器中用来建模（几何和网格）

总体坐标系

在每开场进展一个新的ANSYS分析时，已经有三个坐标系预先定义了。

它们位于模型的总体原点。

三种类型为:

CS,0:

总体笛卡尔坐标系

CS,1:

总体柱坐标系

CS,2:

总体球坐标系

数据库中节点坐标总是以总体笛卡尔坐标系，无论节点是在什么坐标系中创立的。

局部坐标系

局部坐标系是用户定义的坐标系。

局部坐标系可以通过菜单路径Workplane>

LocalCS>

CreateLC来创立。

激活的坐标系是分析中特定时间的参考系。

缺省为总体笛卡尔坐标系。

当创立了一个新的坐标系时，新坐标系变为激活坐标系。

这说明后面的激活坐标系的命令。

菜单中激活坐标系的路径Workplane>

ChangeactiveCSto>

。

节点坐标系

每一个节点都有一个附着的坐标系。

节点坐标系缺省总是笛卡尔坐标系并与总体笛卡尔坐标系平行。

节点力和节点边界条件〔约束〕指的是节点坐标系的方向。

时间历程后处理器/POST26中的结果数据是在节点坐标系下表达的。

而通用后处理器/POST1中的结果是按结果坐标系进展表达的。

例如:

模型中任意位置的一个圆，要施加径向约束。

首先需要在圆的中心创立一个柱坐标系并分配一个坐标系号码（例如CS,11）。

这个局部坐标系现在成为激活的坐标系。

然后选择圆上的所有节点。

通过使用"

Prep7>

Move/Modify>

RotateNodalCStoactiveCS"

选择节点的节点坐标系的朝向将沿着激活坐标系的方向。

未选择节点保持不变。

节点坐标系的显示通过菜单路径Pltctrls>

Symbols>

NodalCS。

这些节点坐标系的X方向现在沿径向。

约束这些选择节点的X方向，就是施加的径向约束。

注意：

节点坐标系总是笛卡尔坐标系。

可以将节点坐标系旋转到一个局部柱坐标下。

这种情况下，节点坐标系的X方向指向径向，Y方向是周向〔theta）。

可是当施加theta方向非零位移时，ANSYS总是定义它为一个笛卡尔Y位移而不是一个转动〔Y位移不是theta位移〕。

单元坐标系

单元坐标系确定材料属性的方向〔例如，复合材料的铺层方向〕。

对后处理也是很有用的，诸如提取梁和壳单元的膜力。

单元坐标系的朝向在单元类型的描述中可以找到。

结果坐标系

/Post1通用后处理器中（位移,应力，支座反力〕在结果坐标系中报告，缺省平行于总体笛卡尔坐标系。

这意味着缺省情况位移，应力和支座反力按照总体笛卡尔在坐标系表达。

无论节点和单元坐标系如何设定。

要恢复径向和环向应力，结果坐标系必须旋转到适当的坐标系下。

这可以通过菜单路径Post1>

Optionsforoutput实现。

/POST26时间历程后处理器中的结果总是以节点坐标系表达。

显示坐标系

显示坐标系对列表圆柱和球节点坐标非常有用（例如,径向，周向坐标）。

建议不要激活这个坐标系进展显示。

屏幕上的坐标系是笛卡尔坐标系。

显示坐标系为柱坐标系，圆弧将显示为直线。

这可能引起混乱。

因此在以非笛卡尔坐标系列表节点坐标之后将显示坐标系恢复到总体笛卡尔坐标系。

20显式隐式分析转换的考前须知

运用ANSYS/LS-DYNA进展隐、显式分析时，由于隐、显式分析过程中所涉及的对象一般会有所不同，ANSYS/LS-DYNA使用手册中只介绍了一种方法，即下面所述的第一种。

实际上，共有两种方法可以选择使用：

1、将隐、显式分析过程中涉及到的所有对象都在隐式分析前建好模型，把隐式分析不需要的对象的所有节点自由度都约束住，进展隐式求解，转换单元类型，进入显式求解阶段，将显式part的约束去除，执行动力松弛求解以便对相应part进展应力初始化，并按照需要施加新的边界和载荷条件，进展显式分析。

2、在隐式分析时只对隐式分析涉及的对象建模，而不考虑显式分析需要的part，完成隐式分析后，单元类型转换完成后，通过定义新的单元类型和材料，创立显式分析所需的模型，生成新的part列表，选择所有节点，读入隐式求解结果文件进展动力松弛求解，对相应part进展应力初始化，施加必要的约束和载荷条件，执行显式求解。

实际上，动力松弛过程是执行一次稳态或是准静态分析，目的就是将隐式分析的结果中的位移、温度结果作为体载荷施加到相关节点上，实现相应部件的应力初始化，作为后续分析的初始条件。

需要注意的是，LS-DYNA中无高阶单元，所以在进展隐式求解时要选择缩减积分的低阶单元。

如果隐式分析使用高阶单元，那么程序无法自动转换单元类型，需要手动转换。

上面所述的是利用ANSYS作为隐式求解器时的操作方法。

我们知道，近几年来，LSTC公司不断加强LS-DYNA程序本身的隐式分析能力，所以我们也可以利用LS-DYNA本身的隐式求解器来完成隐式分析，也根本有两种方法：

1、进展隐式分析时，涉及的关键字主要有：

*control_implicit_solver

*control_implicit_general

*control_implicit_solution

*control_implicit_auto

*control_implicit_dynamics

等。

在这些命令中，设置隐式求解的求解方法〔波前、迭代〕、时间步长等控制参数。

在dyna的输入文件中参加以下命令，

*interface_springback_nike3D

在该关键字中，声明需要进展应力初始化的part号，完成隐式求解后，生成一个nikin文件，包含了相关part的应力应变信息。

在后续的显式分析中，在inputdeck中参加以下命令，

*include

nikin

程序就会自动将存在应力、应变的相关part导入，进展显式分析。

2、另外，可以LS-DYNA的动力松弛方式来对某一构件进展应力初始化。