LS-DYNA常见问题集锦教学文案

一、影响穿透的一些因素解释I.接触厚度接触厚度定义的是一个参数——当接触体/面相互穿透的距离大于接触厚度时，程序将不计算这个接触，即认为没有接触了。

什么是接触厚度与距离？在自动接触中，接触厚度是一个默认值，大概是面厚度的几倍，在普通接触中，接触厚度无穷大。

II.壳厚度和接触厚度1. 壳厚度：影响刚度和单元质量；2. 接触厚度：①决定解除中的厚度偏移量；②并不影响刚度或壳体质量；③默认接触厚度等于壳厚度；④可以在*CONTACT 或*PART_CONTACAT 中直接缩放接触厚度；⑤在穿透节点被释放之前影响最大允许穿透深度。

III.运动速度对穿透的影响如果物体相对运动速度过大，在一个时间步长中所走过的距离会远超过一个单元的尺寸，若缩小时间步长，即缩小在一个时间步长内所走过的距离和单元尺寸的差异，基础检查可以正常进行，若初速度过高，会搜索不到接触，计算会出现问题。

IV.非对称接触算法中，主从面的定义原则①粗网格表面定义为主面，细网格表面为从面；②主从面相关材料刚度相差悬殊，材料刚度大的一面为主面；③平直或者凹面为主面，凸面为从面。

V.接触刚度的影响穿透可以认为是一种虚拟穿透，如果设定的穿透刚度（fkn）值，就可以减小这种穿透，但却不可避免。

如果fkn 值过大，会使到那元刚度病态，而不能求解。

二、穿透的可能解决方案I.接触方面：1. 修改接触类型，尝试自动接触类型：①STS（面面接触），当一个体的表面穿透另外一个体的表面是创建②SS（单面接触），当一个体的表面自身接触或者接触另一个体的表面时创建2. 接触定义存在问题：①增加接触刚度因子②改变接触面的主从设置，将刚体设置为主面，同时使用单向接触③修改关键字CONTROL_CONTACT中RWPNAL=23. 接触穿透距离超过了接触厚度，从而不再计算接触；4. 如果两个接触体的材料属性和网格差别较大，可以修改SOFT值为1 或者2.5. 接触群组设置不直接用PART，将可能接触的地方设置为segment；6. 修改摩擦系数：Fs和Fd通常设置为相同的值，避免额外的噪声产生。

LS‐DYNA求解器
支持单机和并行计算 支持Windows、Unix、Linux等系统 有32位、64位版本 有单精度、双精度版本
LS‐DYNA求解器
2011/11/9
13
2011/11/9
14
湖南大学 汽车车身先进设计制造国家重点实验室
LS‐DYNA基础
湖南大学 汽车车身先进设计制造国家重点实验室
2011/11/9
汽车工业：碰撞分析、气囊设计 、约束系统匹配、部件加工 航空航天：鸟撞、飞机结构冲击动力分析、 碰撞，坠毁、冲 击爆炸及动态载荷、火箭级间分离模拟分析、宇宙垃圾碰撞、 特种复合材料设计 制造业：冲压、锻造、铸造、切割 建筑业：地震安全、混凝土结构、爆破拆除、公路桥梁设计 国防工业：内弹道和终点弹道、装甲和反装甲系统、穿甲弹 与破甲弹设计、战斗部结构设计、空气，水与土壤中爆炸 电子领域：跌落分析、包装设计、热分析、电子封装 生物医学
9
体育器材 2011/11/9
10
湖南大学 汽车车身先进设计制造国家重点实验室
LS‐DYNA基础 LS‐DYNA基本分析流程
几何模型修正 几何模型读取
湖南大学 汽车车身先进设计制造国家重点实验室
LS‐DYNA基础
前处理软件
FEMB LS‐PREPOST HYPERMESH ANSYS PATRAN FEMAP
沙漏的产生
有限元计算的数值积分采用高斯积分方法
2011/11/9
23
2011/11/9
24
4
2011/11/9
湖南大学 汽车车身先进设计制造国家重点实验室
Hourglass
湖南大学 汽车车身先进设计制造国家重点实验室

LS-DYNA常见问题集锦教学文案1 如何处理LS-DYNA中的退化单元？在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。

该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。

实用标准文案
精彩文档 错误名称
原因及修改措施
输入格式不正确
找到错误位置，修改数据格式
浮点溢出，核心代码崩溃 (1) 材料和单元特性定义错误，如零密度、零厚度等
(2) 单元定义错误
(3) 各种曲线定义有问题，如求解时间比定义曲线的
时间长，导致数据外推失效
(4) 单位不一致
(5) 时间步太大，导致求解不稳定
程序一开始就
崩溃
(1) 单位不一致
(2) 边界条件、荷载曲线定义错误
(3) 某一材料未定义
节点速度无限
大
(1) 应用LS-PREPOST 检查最后一帧动画，找出崩
溃的单元和PART
(2) PART 的材料定义有问题，如某个参数设置过大
或与该PART 相关的接触定义有问题，检查接触
的设置
负体积
(1) 材料参数设置错误，选择合适材料模式
(2) 沙漏模式的变形积累，改为全积分单元
(3) 太高的局部接触力，需调整间隙、降低接触刚度
或降低时间步。

1 如何处理LS-DYNA中的退化单元？在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。

该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。

/viewthread.php?tid=68546LS-DYNA在1976年由美国劳伦斯·利沃莫尔国家实验室（Lawrence Livermore National Laboratory）J.O.Hallquist博士主持开发，时间积分采用中心差分格式，当时主要用于求解三维非弹性结构在高速碰撞、爆炸冲击下的大变形动力响应，是北约组织武器结构设计的分析工具。

LS-DYNA的源程序曾在北约的局域网PubicDomain公开发行，因此在广泛传播到世界各地的研究机构和大学。

从理论和算法而言，LS-DYNA是目前所有的显式求解程序的鼻祖和理论基础。

1988年，J.O.Hallquist创建利沃莫尔软件技术公司（Livermore Software Technology Corporation），LS-DYNA开始商业化进程，总体来看，到目前为止在单元技术、材料模式、接触算法以及多场耦合方面获得非常大的进步。

以下为LS-DYNA初学者常见的问题:  一、LS-DYNA与市面上其它的前处理软件兼容性如何？解答：由于LS-DYNA是全球使用率最高的结构分析软件，因此其它的前处理软件与LS-DYNA 是完全兼容的。

在此要强调的是：LS-DYNA的官方前处理程序为FEMB，因为FEMB是专门为LS-DYNA量身订作的前处理程序，有许多设定条件及架构逻辑是其它前处理软件所难望其项背的，为了避免在学习LS-DYNA的过程及操作上产生困扰，强烈建议使用者采用原厂出品的FEMB来做为LS-DYNA的前处理工具，使用者必定更能体会LS-DYNA 直觉式的设定与强大的分析能力。

.二、LS-DYNA似乎很重视「Contact Algorithm」，这是为什么？解答：是的，LS-DYNA很早以前就已经发展「接触算法」，这是因为基础力学所分析的对像均只考虑「力的受体」，故输入条件皆为外力量值。

LYDYNA能量平衡GLSTAT(参见*database_glstat)文件中报告的总能量是下面几种能量的和：内能internal energy动能kinetic energy接触(滑移)能contact(sliding) energy沙漏能houglass energy系统阻尼能system damping energy刚性墙能量rigidwall energyGLSTAT 中报告的弹簧阻尼能”Spring and damper energy”是离散单元(discrete elements)、安全带单元(seatbelt elements)内能及和铰链刚度相关的内能(*constrained_joint_stiffness…)之和。

而内能”InternalEnergy”包含弹簧阻尼能”Spring and damper energy”和所有其它单元的内能。

因此弹簧阻尼能”Spring anddamper energy”是内能”Internal energy”的子集。

由SMP 5434a 版输出到glstat 文件中的铰链内能”joint internal energy”跟*constrained_joing_stiffness 不相关。

它似乎与*constrained_joint_revolute(_spherical,etc)的罚值刚度相关连。

这是SMP 5434a 之前版本都存在的缺失的能量项，对MPP 5434a 也一样。

这种现象在用拉格朗日乘子(Lagrange Multiplier)方程时不会出现。

与*constrained_joint_stiffness 相关的能量出现在jntforc 文件中，也包含在glstat 文件中的弹簧和阻尼能和内能中。

回想弹簧阻尼能”spring and damper energy”，不管是从铰链刚度还是从离散单元而来，总是包含在内能里面。

在MATSUM 文件中能量值是按一个part 一个part 的输出的(参见*database_matsum)。

在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1, N2, N3, N4, N4, N4, N4, N4和N1, N2, N3, N4，N5, N5, N6, N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID 下划分网格，通过*CONTROL_SOLI关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1.采用默认B-T算法，同时利用*control_shell 控制字设置参数BWC= 1,激活翘曲刚度选项；2.采用含有翘曲刚度控制的单元算法,第10号算法。

该算法是针对单元翘曲而开发的算法,处理这种情况能够很好的保证求解的精度。

除了上述方法外,在计算时要注意控制沙漏,确保求解稳定。

3在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign 命令和重启动技术；以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4 关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOFSOLUTIN -Translation ，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。

一、如何得知L‎S-DYNA为‎一套完整的‎软件？<TOP>解答：LS-DYNA一‎直以来都是‎独立软件，原因很简单‎：(1)美国原厂L‎iverm‎o re Softw‎a re Techn‎o logy‎Corp.(简称LST‎C)授权「势流科技」为台湾地区‎L S-DYNA软‎件的独家代‎理商。

(2)LS-DYNA拥‎有LSTC‎官方的前、后处理程序‎：FEMB及‎L S-Pre/Post，无需另外增‎购。

LS-DYNA强‎调「One Model‎、One Code、Multi‎-resul‎t s」，单一核心程‎序即可进行‎I mpli‎c it、Expli‎c it及混‎合Impl‎icit/Expli‎c it的求‎解方式，可以为您解‎决静力、准静力、掉落、冲击、模态分析、振动、预应力/预应变、流固耦合、流场、热应力....等等的力学‎分析，使用者不必‎再额外添购‎模块程序，徒增预算的‎浪费与软件‎学习的负担‎。

二、LS-DYNA与‎市面上其它‎的前处理软‎件兼容性如‎何？<TOP>解答：由于LS-DYNA是‎全球使用率‎最高的结构‎分析软件，因此其它的‎前处理软件‎与L S-DYNA是‎完全兼容的‎。

在此要强调‎的是：LS-DYNA的‎官方前处理‎程序为FE‎MB，因为FEM‎B是专门为‎L S-DYNA量‎身订作的前‎处理程序，有许多设定‎条件及架构‎逻辑是其它‎前处理软件‎所难望其项‎背的，为了避免在‎学习LS-DYNA的‎过程及操作‎上产生困扰‎，强烈建议使‎用者采用原‎厂出品的F‎EMB来做‎为LS-DYNA的‎前处理工具‎，使用者必定‎更能体会L‎S-DYNA直‎觉式的设定‎与强大的分‎析能力。

.三、LS-DYNA似‎乎很重视「Conta‎c t Algor‎i thm」，这是为什幺‎？<TOP>解答：是的，LS-DYNA很‎早以前就已‎经发展「接触算法」，这是因为基‎础力学所分‎析的对像均‎只考虑「力的受体」，故输入条件‎皆为外力量‎值。

$ then imparts an elastic stiffness to prevent further rotation)
$ - Define between rigid body A (part 1) and rigid body B (part 2)
$ - Define a local coordinate system along the revolute axis
铰链
铰链连接刚体 6个基本类型的铰链定义 »球型铰，转动铰 ，万向铰 »平面铰，平动铰 »圆柱铰 节点对(1,2), (3,4)和 (5,6)应该一致，除了圆柱铰和平动铰 铰链仅应用于刚体 关键字
»*CONSTRAINED_JOINT_option »*CONSTRAINED_JOINT_STIFFNESS
另外的铰链 常速度 齿轮 联轴 滑轮 镙杆 齿轮齿条 转动电机 平动电机
可变形开关 材料可以在变形体和刚体之间转换 »模拟具有长时间不变形运动过程的冲击事件时的一种高 效方式（比如，汽车翻滚） *DEFORMABLE_TO_RIGID »在开始时把变形体转换为刚体 »然后在重启动中part转换可以指定 – *RIGID_DEFORMABLE_R2D – *RIGID_DEFORMABLE_D2R *DEFORMABLE_TO_RIGID_AUTOMATIC »part 转换可以自动发生 »激活的时间 »激和的接触
BUILT-IN HYBRID III 假人
HYBRID III 假人
•用*COMPONENT_HYBRID3
关键字生成一个42 个自由度的 Hybrid III 假人
•与标准的FE假人相比无很大的 开销
•一种混合的公式包含可变形的 夹克衫，脸和盆骨

LS-DYNA是一种流体-固体耦合的计算机仿真软件。

在工程实践中，流固耦合固体不正常变形是一个常见的问题，特别是在汽车碰撞、航天航空以及地震工程中。

本文将对LS-DYNA流固耦合中出现的固体不正常变形问题进行探讨和分析。

流固耦合是指流体和固体之间相互作用的物理过程。

LS-DYNA软件利用数值方法模拟这一物理过程，为工程实践提供了可靠的仿真工具。

然而，在实际应用中，由于参数设置、模型建立等因素的影响，时常会出现固体不正常变形的问题，给工程设计和分析带来困扰。

1. 前期模型建立在使用LS-DYNA进行流固耦合仿真之前，首先需要建立合理的模型。

这涉及到对流体和固体的边界条件、材料性质、网格划分等参数的设定。

一些常见的问题包括网格划分不均匀、流体和固体之间的边界条件设置不合理等。

2. 模拟过程中的参数选择在仿真过程中，一些关键的参数选择对于流固耦合固体不正常变形起着至关重要的作用。

时间步长的选择、材料参数的设定等都会直接影响到仿真结果的准确性。

如果这些参数选择不当，很容易导致固体出现不正常的变形。

3. 材料模型的使用LS-DYNA提供了多种材料模型供用户选择，每种模型都有其适用的范围和局限性。

在进行流固耦合仿真时，对于复杂的材料，如橡胶、泡沫塑料等，如何选择合适的材料模型也是一个关键的问题。

如果材料模型选择不当，很容易导致固体的不正常变形。

4. 网格划分及网格质量对于流固耦合仿真，网格的划分和质量对于仿真结果起着至关重要的作用。

不合理的网格划分和低质量的网格都会导致固体不正常变形的问题。

在进行流固耦合仿真时，需要对网格划分和网格质量进行充分的考虑和优化。

在LS-DYNA流固耦合中，固体不正常变形是一个影响仿真结果准确性的重要问题。

为了避免出现固体不正常变形的问题，需要在模型建立、参数选择、材料模型使用、网格划分和网格质量等方面进行全面的考虑和优化。

对于流固耦合中出现的固体不正常变形问题，需要进行详细的分析和调试，以找到问题的根源并加以解决。

1 如何处理LS-DYNA中的退化单元？在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。

该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。

LS-DYNA 常见问题汇总 2.0
资料整理：yuminhust2005 职 务：A2:LS-DYNA 版主 资料来源：网络和自己的总结 更新时间：2008-10-12
LS-DYNA 常见问题汇总 2.0
yuminhust2005
2008-10...................................................................................................................................................1 2. 单位制度 .................................................................................................................................................................5 3. 质量缩放 .................................................................................................................................................................6 4. 长分析时间 .............................................................................................................................................................7 5. 准静态 .....................................................................................................................................................................8 6. 计算不稳定 .............................................................................................................................................................9 7. 负体积 ...................................................................................................................................................................10 8. 能量平衡 ...............................................................................................................................................................11 9. 沙漏控制 ...............................................................................................................................................................13 10. 阻尼 ...................................................................................................................................................................14

错误名称
原因及修改措施
输入格式不正确
找到错误位置，修改数据格式
浮点溢出，核心代码崩溃
（1）材料和单元特性定义错误，如密度、零厚度
（2）单元定义错误
（3）各种曲线定义有问题，如求解时间比定义曲线的时间解不稳定
程序一开始就崩溃
（1）单位不一致
（2）边界条件、载荷曲线定义错误
（3）某一材料未定义
节点速度无限大
（1）应用LS-PREPOST检查最后一帧动画，找出崩溃的单元和PART
（2）PART的材料定义有问题，如某个参数设置过大或与该PART相关的接触定义有问题，检查接触的设置
负体积
（1）材料参数设置错误，选择合适材料模式
（2）沙漏模式的变形积累，改为全积分单元
（3）太高的局部接触力，需调整间隙、降低接触刚度或降低时间步

几乎看了本版块所有关于转动的讨论，现在终于对刚体的转动彻底的了解了，现在把前一段时间的心得拿出来，希望能对和我一样正在为这个问题迷茫的人有所帮助，以下两个文件是用ansys命令流写的，因为关于转动的k文件在“关于刚体转动的问题”这个主题里面已经讨论的很清除了，相信有很多人都是用ansys做前处理和后处理的，所以即使看懂了k文件，也未必就知道用ansys的前处理怎么来实现。

关于刚体转动，主要要注意以下几点：1. 如果是绕自己的质心转动，那很简单，不用设转动惯量和转动轴，因为ls-dyna默认就是绕着质心转动的。

2.如果是绕任意轴转动，那么就必须设新的质心center、惯性张量和转动轴。

这主要是用edlcs和edipart两个命令来实现，前一个设转动轴，后一个设center和惯性张量，关于惯性张量可以用solidworks来计算（wswu在那个主题曾提到这个方法），当然惯性张量也可以自己估算,如果质量体相对坐标轴比较远就可以把质量体当作集中质量，比如计算IXY=m*x*y(x和y是把质量中心投影到xy平面后的那个点距x和y的距离）。

有人曾提到惯性张量可以随便给，当然随便给可能能算出结果来，但是那个结果可能并不是你所想要的结果，所以，如果是针对课题，这些东西不妨花点时间，弄出一个确定的值。

其实也没有什么太多要注意的，以上两点都做好了，你肯定就ok了。

我这两个文件借助了iamaxelong提出的用刚体带动弹性体的方法，因为iamaxelong老大那个ansys命令流贴在网页上被弄了一些乱码在里面，复制过来之后要做大量的修改，而且分析所需要的时间太长，所以我这里整了一个特别简单的例子，我的电脑分析只需要1分钟，可以帮助你快速看到分析的结果，相信即使比我还菜的人也能很快了解刚体转动的问题啦。

/PREP7!定义单元和材料ET,1,SOLID164MP,DENS,1,7800MP,EX,1,2.1e11MP,NUXY,1,0.3EDMP,RIGI,2,7,4MP,DENS,2,7800MP,EX,2,2.1e11MP,NUXY,2,0.3!画刚体和弹性体的实体模型BLOCK,-0.1,0.1,0.4,0.6,0,0.05,BLOCK,-0.05,0.05,0.45,0.55,0.05,0.1,KWPAVE, 9wpro,,90.000000,VSBW, 1KWPAVE, 10VSBW, 3VGLUE,all!确定实体模型上线条的网格密度LESIZE,all,0.025, , , , , , ,1lesize,11,,,2lesize,12,,,2lesize,31,,,2lesize,32,,,2lesize,9,,,2lesize,10,,,2lesize,39,,,2lesize,40,,,2lesize,21,,,2lesize,22,,,2lesize,23,,,2lesize,24,,,2numcmp,volu!划分单元vsweep,2vsweep,3vsweep,4mat,2vsweep,1!建立partedpart,create!定义载荷数组*DIM,time,ARRAY,2,1,1, , ,*SET,TIME(2,1,1) , 10*DIM,wz,ARRAY,2,1,1, , ,*SET,WZ(1,1,1) , 5*SET,WZ(2,1,1) , 5!加载EDLOAD,ADD,RBOZ,0, 2,TIME,WZ, 0, , , , ,!定义分析选项FINISH/SOLTIME,0.5,EDRST,30,EDHTIME,50,EDDUMP,1,EDENERGY,1,1,1,1/STATUS,SOLUSOLVE!输出动画FINISH/POST1SET,LAST/DSCALE,1,1.0/REPLOTPLDI, ,ANTIME,50,0.1, ,1,2,0,0FINISH/PREP7!定义单元和材料ET,1,SOLID164MP,DENS,1,7800MP,EX,1,2.1e11MP,NUXY,1,0.3EDMP,RIGI,2,7,4MP,DENS,2,7800MP,EX,2,2.1e11MP,NUXY,2,0.3!画刚体和弹性体的实体模型BLOCK,-0.1,0.1,0.4,0.6,0,0.05, BLOCK,-0.05,0.05,0.45,0.55,0.05,0.1,KWPAVE, 9wpro,,90.000000,VSBW, 1KWPAVE, 10VSBW, 3VGLUE,all!确定实体模型上线条的网格密度LESIZE,all,0.025, , , , , , ,1 lesize,11,,,2lesize,12,,,2lesize,31,,,2lesize,32,,,2lesize,9,,,2lesize,10,,,2lesize,39,,,2lesize,40,,,2lesize,21,,,2lesize,22,,,2lesize,23,,,2lesize,24,,,2numcmp,volu!划分单元vsweep,2vsweep,3vsweep,4mat,2vsweep,1!建立partedpart,create!定义载荷数组*DIM,time,ARRAY,2,1,1, , ,*SET,TIME(2,1,1) , 10*DIM,wz,ARRAY,2,1,1, , ,*SET,WZ(1,1,1) , 5*SET,WZ(2,1,1) , 5!自定义质心和旋转轴!质心定义在原点，模型本来的质心不在原点*DIM,center,ARRAY,3,1,1, , ,!定义惯性张量，反应了实体的质量相对自定义质心的分布*DIM,II,ARRAY,6,1,1, , ,*SET,II(1,1,1) , 1.1*SET,II(4,1,1) , 0.02475*SET,II(5,1,1) , 0.165*SET,II(6,1,1) , 1.1edipart,2,add,center,4.4,,II!定义旋转轴EDLCS,ADD,11,0,0,0,0,0,10,0,0,0,!加载EDLOAD,ADD,RBOZ,11, 2,TIME,WZ, 0, , , , ,!定义分析选项FINISH/SOLTIME,0.5,EDRST,30,EDHTIME,50,EDDUMP,1,EDENERGY,1,1,1,1/STATUS,SOLUSOLVE!输出动画FINISH/POST1SET,LAST/DSCALE,1,1.0/REPLOTPLDI, ,ANTIME,50,0.1, ,1,2,0,0FINISH用关键字 *INITIAL_VELOCITY_GENERATION 可以实现。

lsdyna常见问题汇总LYDYNA能量平衡GLSTAT(参见*database_glstat)⽂件中报告的总能量是下⾯⼏种能量的和：内能internal energy动能kinetic energy接触(滑移)能contact(sliding) energy沙漏能houglass energy系统阻尼能system damping energy刚性墙能量rigidwall energyGLSTAT 中报告的弹簧阻尼能”Spring and damper energy”是离散单元(discrete elements)、安全带单元(seatbelt elements)内能及和铰链刚度相关的内能(*constrained_joint_stiffness…)之和。

⽽内能”Internal Energy”包含弹簧阻尼能”Spring and damper energy”和所有其它单元的内能。

因此弹簧阻尼能”Spring and damper energy”是内能”Internal energy”的⼦集。

由SMP 5434a 版输出到glstat ⽂件中的铰链内能”joint internal energy”跟*constrained_joing_stiffness 不相关。

它似乎与*constrained_joint_revolute(_spherical,etc)的罚值刚度相关连。

这是SMP 5434a 之前版本都存在的缺失的能量项，对MPP 5434a 也⼀样。

这种现象在⽤拉格朗⽇乘⼦(Lagrange Multiplier)⽅程时不会出现。

与*constrained_joint_stiffness 相关的能量出现在jntforc ⽂件中，也包含在glstat ⽂件中的弹簧和阻尼能和内能中。

回想弹簧阻尼能”spring and damper energy”，不管是从铰链刚度还是从离散单元⽽来，总是包含在内能⾥⾯。

1 如何处理LS-DYNA中的退化单元？在网格划分过程中，我们常遇到退化单元，如果不对它进行一定的处理，可能会对求解产生不稳定的影响。

在LS-DYNA中，同一Part ID 下既有四面体，五面体和六面体，则四面体，五面体既为退化单元，节点排列分别为N1，N2，N3，N4，N4，N4，N4，N4和N1，N2，N3，N4，N5，N5，N6，N6。

这样退化四面体单元中节点4有5倍于节点1-3的质量，而引起求解的困难。

其实在LS-DYNA的单元公式中，类型10和15分别为四面体和五面体单元，比退化单元更稳定。

所以为网格划分的方便起见，我们还是在同一Part ID下划分网格，通过*CONTROL_SOLID关键字来自动把退化单元处理成类型10和15的四面体和五面体单元。

2 LS-DYNA中对于单元过度翘曲的情况有何处理方法有两种方法：1. 采用默认B-T算法，同时利用*control_shell控制字设置参数BWC＝1，激活翘曲刚度选项；2. 采用含有翘曲刚度控制的单元算法，第10号算法。

该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。

除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。

3 在ANSYS计算过程中结果文件大于8GB时计算自动中断，如何解决这个问题？解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。

特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。

具体操作方法和过程参见下列命令流文件的演示。

4关于梁、壳单元应力结果输出的说明问题：怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。