笔记摘录（隐式分析专题）
1、显式与隐式区别
0、隐式、显式算法区别
隐式求解法不考虑惯性效应[C]和[M]。对于线性问题，无条件稳定，可以用大的时间步。对于非线性问题，通过一系列线性逼近（Newton-Raphson）来求解；要求转置非线性刚度矩阵[K]，收敛时候需要小的时间步，对于高度非线性问题无法保证收敛。因此，隐式求解一般用于线性分析和非线性结构静动力分析，包括结构固有频率和振型计算。 ansys使用的Newmark时间积分法即为隐式求解法；
由于两者解题的出发点，所以一般来说显式用于求解和时间相关的动力学问题。隐式用来求解和时间无关的静力学问题。但也不是绝对的；
隐式存在收敛问题；显式不存在收敛问题，但求解时间要长一些；而且可能存在沙漏能问题；显式算法基于动力学方程，因此无需迭代；而静态隐式算法基于虚功原理，一般需要迭代计算；
显式算法 最大优点是有较好的稳定性。 动态显式算法采用动力学方程的一些差分格式（如广泛使用的中心差分法、线性加速度法、Newmark法和wilson法等），不用直接求解切线刚度，不需要进行平衡迭代，计算速度快，时间步长只要取的足够小，一般不存在收敛性问题。因此需要的内存也比隐式算法要少。并且数值计算过程可以很容易地进行并行计算，程序编制也相对简单。但显式算法要求质量矩阵为对角矩阵，而且只有在单元级计算尽可能少时速度优势才能发挥, 因而往往采用减缩积分方法，容易激发沙漏模式，影响应力和应变的计算精度；
1、LS-DYNA显式求解：
采用中心差分方法进行显式时间积分
l方程非耦合，可以直接求解(显式)
l但需要常小的时间步保持稳定状态
l不需要求解刚度矩阵
l适合冲击、穿透等高频非线性动力响应问题
2 LS-DYNA隐式求解：
采用Newmark隐式时间积分
l对于线性问题，无条件稳定
l可以采用大的时间步
l对于非线性问题，需要一系列线性逼近（Newton-Raphson）叠代求解
l需要求解刚度矩阵
l适合静力问题、低频动力问题及特征值分析。
3、LS-DYNA在进行Implicit分析时，有没有什幺该注意的地方？
解答：您这个问题好象有点「蛮广泛」...
在这边先跟您说明一下大原则：Implicit method分析主要就是在求解矩阵式，给定的力学条件不外乎是「Force Prescribed」与「Displacement Prescribed」，一般最常犯的毛病就是边界条件给定不足，造成model发生「刚体运动」的情形，您可以将这种情形想象成物体浮在空中一样，宜尽量避免这种情形产生。 在此提供您一个判断的方式：您可以先利用LS-DYNA eigenvalue an alysis的功能先进行模态分析，藉后处理程序检查有无刚体模态，以判断Implicit model的束制条件是否充足。
4、显式隐式分析转换的注意事项
运用ANSYS/LS-DYNA进行隐、显式分析时，由于隐、显式分析过程中所涉及的对象一般会有所不同，ANSYS/LS-DYNA使用手册中只介绍了一种方法，即下面所述的第一种。实际上，共有两种方法可以选择使用：
1、 将隐、显式分析过程中涉及到的所有对象都在隐式分析前建好模型，把隐式分析不需要的对象的所有节点自由度都约束住，进行隐式求解，转换单元类型，进入显式求解阶段，将显式part的约束去除，执行动力松弛求解以便对相应part进行应力初始化，并按照需要施加新的边界和载荷条件，进行显式分析。
2、 在隐式分析时只对隐式分析涉及的对象建模，而不考虑显式分析需要的part，完成隐式分析后，单元类型转换完成后，通过定义新的单元类型和材料，创建显式分析所需的模型，生成新的part列表，选择所有节点，读入隐式求解结果文件进行动力松弛求解，对相应part进行应力初始化，施加必要的约束和载荷条件，执行显式求解。
实际上，动力松弛过程是执行一次稳态或是准静态分析，目的就是将隐式分析的结果中的位移、温度结果作为体载荷施加到相关节点上，实现相应部件的应力初始化，作为后续分析的初始条件。需要注意的是，LS-DYNA中无高阶单元，所以在进行隐式求解时要选择缩减积分的低阶单元。如果隐式分析使用高阶单元，则程序无法自动转换单元类型，需要手动转换。
1、进行隐式分析时，涉及的关键字主要有：
*control_implicit_solver
*control_implicit_general
*control_implicit_solution
*control_implicit_auto
*control_implicit_dynamics
2、另外，可以LS-DYNA的动力松弛方式来对某一构件进行应力初始化。
相关的关键字为：
*DEFINE_CURVE
将此卡片的SIDR参数设置为1即可启动动力松弛分析。
*CONTROL_DYNAMIC_RELAXATION
此卡片在随后的显式分析中用来进行应力初始化操作。
*LOAD_BODY_RX（RY、RZ）等
运行后收敛的结果即为初始化应力，同时生成动力松弛文件drdisp.sif，该文件与drelax文件结构、用法完全一致，只是精度上较差。
建议：使用ANSYS作为隐式求解器，因为它的隐式功能和计算精度都优于LS-DYNA。
2、LS-DYNA隐式求解
作为著名的显式求解程序只能求解瞬态动力问题，对于时间历程较长的静、动力问题， LS-DYNA的显式中心差分法有它的局限性，而一些与瞬态动力分析紧密相关的问题要求LS-DYNA也能够求解，如：
Ø冲压成型后的回弹计算【显式联合隐士】
Ø应力初始化
Ø冲击后常时间低频动力响应
Ø静力分析
Ø特征值分析
Ø频率响应分析
1、隐式的激活-隐式/显式_GENERAL
LS-DYNA中隐式分析的激活及相关关键字
在LS-DYNA中，缺省的求解是显式求解，可以通过下面的关键字来激活隐式求解：
*CONTROL_IMPLICIT_GENERAL
$ imflag dt0 iefs nsteps b igso
1 0.01 0 0 0
其中参数
imflag=1激活全隐式求解
imflag=0为缺省的显式求解。
imflag=2为显式求解后无缝进行隐式回弹求解。
该关键字对于所有隐式求解分析来说都是必需的。
与隐式求解相关的其它关键字：
*CONTROL_IMPLICIT_LINEAR
(v960版本改为*CONTROL_IMPLICIT_SOLVER)
*CONTROL_IMPLICIT_NONLINEAR
(在v960版后改为*CONTROL_IMPLICIT_SOLUTION)
*CONTROL_IMPLICIT_AUTO
*CONTROL_IMPLICIT_STABILIZATION
*CONTROL_IMPLICIT_DYNAMICS
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_IMPLICIT_BUCKLE（v970）
*CONTROL_IMPLICIT_MODES（v970）
使用*CONTROL_IMPLICIT_GENERAL激活隐式求解后，还需要与上面的相关关键字进行联合，从而对不同的问题设置不同的关键字进行求解。
2、隐式关键字2-线性-刚度矩阵_SOLVER
*CONTROL_IMPLICIT_LINEAR
(v960版后改为*CONTROL_IMPLICIT_SOLVER)
*CONTROL_IMPLICIT_SOLVER
$ lsolver prntfl g negeig
4 0 0
$
该关键字用来设置线性方程求解器的相关参数，主要用来求解下列线性方程：
要注意的是该关键字不是激活线性分析，该求解器的功能主要是对刚度矩阵进行转置计算（线性分析还是非线性分析由下面的关键字指定）。
3、隐式关键字3-非线性/线性_SOLUTION
*CONTROL_IMPLICIT_NONLINEAR
(在v960版后改为*CONTROL_IMPLICIT_SOLUTION)
*CONTROL_IMPLICIT_SOLUTION
$ nlsolvr ilimit maxref dctol ectol rctol lstol
2 0 0 0.0 0.0 0 0
$ dnorm divflag inistif nlprint
0 0 0 0
该关键字是用来指定是线性分析还是非线性分析，若是非线性分析，则指定非线性叠代求解器类型、控制平衡叠代搜索次数和各种容限。
若设置nlsolvr=1，则为线性分析，实际上相当于叠代次数为1的非线性分析；
若设置nlprint=1，则屏幕上输出如下，在求解过程中，也可以使用CTRL+ nlprint来输出。
4、隐式关键字4-时间步AUTO
*CONTROL_IMPLICIT_AUTO
*CONTROL_IMPLICIT_AUTO
$ iauto iteopt itewin dtmin dtmax
1 0 0 0.0 0.0
该关键字控制激活自动步长控制，缺省是固定时间步长，当平衡叠代失效后，可以试着改为自动时间步长。
5、隐式关键字5-冲压回弹-人工粘度STABILIZATION
*CONTROL_IMPLICIT_STABILIZATION
该关键字仅用于多步回弹分析，施加人工粘度，提高回弹分析的收敛性，具体应用见冲压成型分析章节。
6、隐式关键字6-静力/动力DYNAMICS
*CONTROL_IMPLICIT_DYNAMICS
*CONTROL_IMPLICIT_DYNAMICS
$ imass gamma beta
0 0.0 0.0
该关键字激活是静力还是动力分析的选项。
Imass=0为静力分析
Imass=1低频动力分析（采用Newmark时间积分）
Imass=2低频动力分析(采用模态综合法）
7、隐式关键字7-特征值分析EIGENVALUE
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_IMPLICIT_EIGENVALUE
$ neig
15
该关键字激活特征值分析，可输出各阶模态，neig=15表示输出前15阶模态。
8、隐式单元-LS-DYNA隐式单元公式SECTION_OPTION
一般情况下，不建议使用LS-DYNA缺省的显式单点积分单元来求解隐式问题， LS-DYNA开发有专门的各种隐式单元应用于隐式分析中，下面具体介绍各种单元：
使用*SECTION_OPTION来定义。
8-1 隐式壳单元公式：
1．S/R Hughes-Liu单元：6号单元公式
l2 x2选择性缩减积分
l每个节点6个自由度(dx, dy, dz, rx, ry, rz)
l在*CONTROL_SHELL中设置统一的节点法向
l可以退化为三角形单元公式
l CPU开销很大
2．快速壳：16号单元公式
l 每个节点6个自由度(dx, dy, dz, rx, ry, rz)
l面内2 x2积分点，厚方向建议5个积分点
l不能退化为三角形
l利用沙漏控制8增加翘曲刚度，改善收敛性
l非常适合用于回弹计算
3．线性壳：18号单元公式
l每个节点6个自由度(dx, dy, dz, rx, ry, rz)
l面内2 x2积分点
l由于是线弹性，厚度方向2个积分点足够
l包含三角形壳公式
l仅应用于线性问题和特征值分析
l必须使双精度求解器
8-2 隐式体单元公式：
1．积分S/R体单元：2号单元公式
l每个节点3个自由度(dx, dy, dz)
l2 x2 x 2选择性缩减积分
lCPU开销大
2．线性体单元：18号单元公式
l每个节点3个自由度(dx, dy, dz)
l2 x2 x 2积分点
l仅应用于线性问题和特征值分析
l必须使用双精度求解器
8-3 隐式梁单元公式：
l几号？就WB默认的1号即可；
l每个节点6个自由度(dx, dy, dz, rx, ry, rz)
l长度方向一积分点
l截面积分
l横向剪切修正因子
l可移动的中面方式
ls，t方向由向节点n3控制
8-4 隐式平面应变单元：13号单元公式：
8-5 隐式轴对称单元：15号单元公式
l特殊条件下的4节点壳单元
lX-Y平面建模
lY轴为15号单元的对称轴
l2个自由度(dx,dy)
l建议约束自由度DOF (dz, rx, ry, rz)
l建议设置NIP=4
l非常快，很好的收敛特性
9、隐式LS-DYNA隐式求解材料：
LS-DYNA开发的材料主要应用于显式分析，并不全部支持隐式功能（但在970中大多数材料都支持隐式功能）
不同的单元公式材料支持的程度也不一样，下面为960的材料支持程度：
l3D Solid单元
1,2,3,4,5,6,7,9,12,13,18,20,24,26,30,49,57,63,91,92,103,104,106,115
lShell单元
1,2,3,4,6,9,18,20,24,36,37,49,91,92,103,104,106,116,117,118
lBeam单元
1,3,4,6,9,18,20,24,100
l2D Solid单元
1,2,3,4,5,6,7,9,12,13,18,20,24,26,57,63
10、LS-DYNA隐式接触界面：
LS-DYNA能用于隐式求解的接触界面如下：
l所有的隐式接触除了TIED方式外都使用罚函数方式
l节点发射逻辑在隐式分析中自动失效（可选卡片B中参数SNLOG=1）
l建议使用重定向法向矢量（ORIEN=1或2）
l自动接触类型在隐式求解中常失效，不建议使用
其他边界条件、约束条件、载荷、输出文件等与显式分析大同小异，这里不再详述。
11、隐式求解非线性收敛问题：
应用LS-DYNA进行隐式求解常遇到的问题就是收敛问题，如：
l到达平衡叠代步数—不满足位移和能量容限，放弃该子步
l发散—不平衡力增加（而不是减小），重新构造刚度矩阵，继续叠代
l能量爆炸—能量相比初始值增加百万倍，放弃该子步
l负的特征值— 这是在使用线性方程求解器在转置刚度矩阵时出现的错误，缺省状态下忽略，试着继续运算
l线性搜索子步尺寸为零— 新的位移不降低不平衡力；在没有到达最后叠代容限情况下重新构造刚度矩阵，否则放弃该子步
解决收敛问题的步骤：
l决定终止的原因（检查错误文件）
l激活输出开关得到更多的信息（NLPRINT=1）
l在叠代过程中使用D3ITER检查几何体变形情况
l仔细坚持关键字输入文件
上述收敛问题的解决方法：
12、几种诊断工具：
l使用NLPRINT输出开关（NLPRINT=1）或〈CTRL+C〉NLPRINT得到诊断信息。
l负的特征值开关可以使程序在遇到负的特征值时立即用更小的时间步去计算。
l〈CTRL+C〉CONV可以采用力收敛容限
l〈CTRL+C〉ITER创建每次叠代的D3ITER数据库文件
13、到达平衡叠代步数错误信息：
l观察位移和能量规范的收敛过程（查看D3HSP文件，或激活NLPRINT）
l如果发现收敛容限几乎要满足：
— 增加更多的叠代步数
—稍微放宽收敛容限
— 使用〈CTRL+C〉CONV采用力收敛容限
l如果收敛过程非常缓慢
— 改为全牛顿叠代方法（ILIMIT=1）
—检查是否有不正确的材料特性
— 增加时间步长
l如果收敛过程不稳定
— 降低接触刚度（10倍量级缩放因子）
— 降低时间步尺寸
14、发散和能量爆炸错误信息：
l一些PART模型可能变形太快
l使用〈CTRL+C〉ITER创建每次叠代的D3ITER数据库文件
l在叠代过程中观察模型的变形情况
l检查是否有过大的载荷
l检查是否有不正确屈服应力或硬化特性
l检查接触定义，是否有过大的渗透
l降低时间步尺寸
15、负的特征值错误信息：
出现负的特征值有四种可能性：
1.不正确的材料特性：
l负杨氏模量，负的应力应变斜率等
l检查D3HSP文件确保LS-DYNA正确读入输入关键字文件
2．单元法向翻转：
l在初始网格中有极度变形的单元
l接触初始化过程中由于节点的重新分布导致网格变形
l由于过度的变形导致单元扭曲
3．刚体模式
l检查边界条件和约束条件
l关注梁单元绕轴的自由旋转
4．模型中的某些特征在隐式分析中并不支持
l当使用非法的特征时，警告信息并不总是出现
l参考户手册仔细检查关键字文件
16、特征值（模态）分析EIGENVALUE
不同的应用需要的关键字不一样，如做特征值分析，需要3个关键字：
如图所示，对一个顶端固支的薄板进行特征值分析，输出前15阶模态：使用
LS-DYNA的隐式求解进行分析，并和ANSYS得到的结果进行对比。
知道应用LS-DYNA做特征值分析，需要3个关键字：
*CONTROL_IMPLICIT_GENERAL
*CONTROL_IMPLICIT_EIGENVALUE
*CONTROL_TERMINATION【貌似不需要】
需要使用双精度的求解器，单元公式用18号壳单元公式（见下面关键字文件），求解完后自动生成两个文件：D3EIGV和EIGOUT
D3EIGV：与D3POLT文件一样的二进制输出文件，包含每一模态下的振型，可以用LS-PREPOST显示，下面为LS-PREPOST显示的几阶模态振型：
EIGOUT：ASCII文件，是所有输出模态的文档总结。
第六阶
模态关键字：
*KEYWORD
$
*TITLE
Implicit Example 6 - Eigenvalue An alysis
*CONTROL_IMPLICIT_EIGENVALUE
$控制为模态特征值求解，输出前8阶频率及振型
8 0.0 0 0.0 0 0.0 0 0.0
0 0 0 0 0
*CONTROL_IMPLICIT_GENERAL
$控制为隐式求解
1 0.1 0 0 0 0 0 0
*MAT_ELASTIC
$材料
1 7.85E-9 210000.0 0.3 0.0 0.0 0.0
*SECTION_SHELL
$18号壳单元适用于特征值问题
1 18 0.0 0 0.0 0.0 0 0
1.0 1.0 1.0 1.0 0.0 0.0 0.0 0
*PART
$PR_PART_COL 1 78005561
Sheet
1 1 1 0 0 0 0 0
*NODE
$前处理生成
*ELEMENT_SHELL
*BOUNDARY_SPC_SET
$边界-固定-前处理生成
1 0 1 1 1 1 1 1
*SET_NODE_LIST
$双精度呢？
*END
模态后处理：
Lsprepost中，使用plot方式，打开D3EIGV文件即可；Fcomp中选择位移，既可以查看频率及振型（通过动画按钮）；
Eigout文件好像就自动导入了，就可以查看频率相关信息了；
17、悬臂梁模态分析-WB下前处理
WB的DM中建立梁的截面、线等模型
（wb-lsdyna）中建立网格、约束、载荷（不加载不能导出K文件，导出后再删除载荷即可）；
K文件插入模态相关关键字即可；
关键字
可输出模态应力；
*KEYWORD
*TITLE
$Created with ANSYS Workbench v17.0
$
$Units: mm, mg, ms, mN, K
$
$*DATABASE_FORMAT
$ 1IFORM 2IBINARY
$ 0
$
*NODE
*ELEMENT_BEAM
*CONTROL_IMPLICIT_EIGENVALUE
$控制为模态特征值求解，输出前8阶频率及振型，最后一个1可以输出应力；
8 0.0 0 0.0 0 0.0 0 0.0
0 0 0 0 0/1
*CONTROL_IMPLICIT_GENERAL
$控制为隐式求解
1 0.1 0 0 0 0 0 0
*SECTION_BEAM
$WB建立
$ 1SECID 2ELFORM 3SHRF 4QR/IRID 5CST 6SCOOR 7NSM
1 1 1 2 1
$ 1TS1 2TS2 3TT1 4TT2 5NSLOC 6NTLOC
6 6 0 0 0 0
*MAT_ELASTIC
$ 1MID 2RO 3E 4PR 5DA 6DB 7K
1 7.85 2e+008 0.3
*PART
$ HEADING
Line Body
$ 1PID 2SECID 3MID 4EOSID 5HGID 6GRAV 7ADPORT 8TMID
1 1 1 0 0
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ LOAD DEFINITIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DEFINE_CURVE
$没用到
$ 1LCID 2SIDR 3SFA 4SFO 5OFFA 6OFFO 7DATTYP
1
$ 1A 2O
0 1e+011
1 1e+011
10 1e+011
*DEFINE_CURVE
$ 1LCID 2SIDR 3SFA 4SFO 5OFFA 6OFFO 7DATTYP
2
$ 1A 2O
0 100000
1 100000
10 100000
$
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ CONTROL OPTIONS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$没用到
$*CONTROL_TERMINATION
$ 1ENDTIM 2ENDCYC 3DTMIN 4ENDENG 5ENDMAS
$ 1 10000000 0.01 10 0
$*CONTROL_TIMESTEP
$ 0 0.9 0 0 0 1 1 0
$*CONTROL_HOURGLASS
$ 1IHQ 2QH
$1为沙漏公式，不同单元类型建议使用不同单元公式，如shell-2号壳建议用4或5号；实体就1号；
$ 1 0.1
$*CONTROL_BULK_VISCOSITY
$ 1Q1 2Q2 3TYPE
$ 1.5 0.06 -2
$*DAMPING_GLOBAL
$ 0 0
$*CONTROL_ENERGY
$ 1HGEN 2RWEN 3SLNTEN 4RYLEN
$ 2 1 2 2
$*CONTROL_ACCURACY
$ 1OSU 2INN 3PIDOSU
$ 1 4
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$ LIST SETS $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$约束
*SET_NODE_LIST
$ 1SID 2DA1 3DA2 4DA3 5DA4
1
$ 1NID1 2NID2 3NID3 4NID4 5NID5 6NID6 7NID7 8NID8
1
*BOUNDARY_SPC_SET
$ 1NSID 2CID 3DOFX 4DOFY 5DOFZ 6DOFRX 7DOFRY 8DOFRZ
1 0 1 1 1 1 1 1
*END
后处理：
类似壳模态后处理，可以显示梁截面；
18、壳静力分析-WB下做-加入相关关键字即可
关键字如下：
那些可以前处理加？
*KEYWORD
*TITLE
Implicit Example 1 - Simple Cantilever Beam
$ =============
*CONTROL_BULK_VISCOSITY
0.0 0.0 -2 0
*CONTROL_IMPLICIT_GENERAL
1 0.5 0 0 0 0 0 0
*CONTROL_IMPLICIT_SOLUTION
2 0 0 0.0 0.0 0.0 0.0 0.0
*CONTROL_IMPLICIT_SOLVER
4 0 0 0 0 0.0 0 0.0
*CONTROL_TERMINATION
$时间
1.0 0 0.0 0.0 0.0
$
*DATABASE_BINARY_D3PLOT
1.0E-2 0 0 0
*DATABASE_BINARY_D3THDT
1.0E-3 0
*DATABASE_HISTORY_NODE
$输出节点历程
93
*MAT_PLASTIC_KINEMATIC_TITLE
Steel
1 7.85E-9 210000.0 0.3 200.0 1000.0 0.0
0.0 0.0 0.0 0.0
*SECTION_SHELL_TITLE
1mm
1 16 0.0 0 0.0 0.0 0 0
1.0 1.0 1.0 1.0 0.0 0.0 0.0 0
*PART
Cantilever Beam
1 1 1 0 0 0 0 0
*NODE
*ELEMENT_SHELL
$载荷
*DEFINE_CURVE
$ Cross-reference summary for Load-curve 1
$ -----------------------------------------------
$
$ Loading definition <No label>: Nodal force vs time
$ X axis : Time (Units: Time)
$ Y axis : Nodal force (Units: Force)
$
$ Usage: Transient a nalysis
$
1 0 0.0 0.0 0.0 0.0 0
0.0 0.0
0.89999998 0.77800000
1.1000000 0.77800000
$
*LOAD_NODE_POINT
$93号节点施加曲线定义的载荷；
93 3 1 -1.0 0
$边界
*BOUNDARY_SPC_SET
1 0 1 1 1 1 1 1
$
*SET_NODE_LIST_TITLE
SPC
1 0.0 0.0 0.0 0.0
1 32 63 94 125
$
$
*CONSTRAINED_NODAL_RIGID_BODY
$1为part号或刚性体号，0为坐标系，2为节点set号
1 0 2 0 0 0 0
$
*SET_NODE_LIST
2 0.0 0.0 0.0 0.0
124 155 93 62 31*end
*end
19、梁静力分析-WB下做-加入相关关键字即可
关键字：
*KEYWORD
$
*TITLE
$*DATABASE_FORMAT
$ 1IFORM 2IBINARY
$ 0
*node
*element
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$时间控制、隐式控制、非线性控制等
*CONTROL_BULK_VISCOSITY
$人工体积粘性
0.0 0.0 -2 0
*CONTROL_IMPLICIT_GENERAL
$0为显示分析，1为隐式分析,0.5为隐式分析的时间步
1 0.5 0 0 0 0 0 0
*CONTROL_IMPLICIT_SOLUTION
$1为线性分析，2为非线性分析BFGS updata
2 0 0 0.0 0.0 0.0 0.0 0.0
*CONTROL_IMPLICIT_SOLVER
$4为默认的SMP并行求解矩阵
4 0 0 0 0 0.0 0 0.0
*CONTROL_TERMINATION
$结束时间
1.0 0 0.0 0.0 0.0
$*CONTROL_TIMESTEP
$ 0 0.9 0 0 0 1 1 0
$
$
$ ==============
$ DATABASE cards
$ ==============
$输出
*DATABASE_BINARY_D3PLOT
1.0E-2 0 0 0
*DATABASE_BINARY_D3THDT
1.0E-2 0
$
*DATABASE_HISTORY_NODE
$93号节点定义为输出到history及ascii文件；
2
*DATABASE_EXTENT_BINARY
$没这个输出就没应力！
0 0 3 1 0 0 0 0
0 0 4 0 0 0
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*SECTION_BEAM
$单元、材料、part
$ 1SECID 2ELFORM 3SHRF 4QR/IRID 5CST 6SCOOR 7NSM
1 1 1 2 1
$ 1TS1 2TS2 3TT1 4TT2 5NSLOC 6NTLOC
6 6 0 0 0 0
$
*MAT_ELASTIC
$ 1MID 2RO 3E 4PR 5DA 6DB 7K
1 7.85 2e+008 0.3
$
*PART
$ HEADING
Line Body
$ 1PID 2SECID 3MID 4EOSID 5HGID 6GRAV 7ADPORT 8TMID
1 1 1 0 0
*DEFINE_CURVE
$ 1LCID 2SIDR 3SFA 4SFO 5OFFA 6OFFO 7DATTYP
1
$ 1A 2O
0 1e+011
1 1e+011
10 1e+011
*DEFINE_CURVE
$载荷曲线用lid2；WB中定义；
$ 1LCID 2SIDR 3SFA 4SFO 5OFFA 6OFFO 7DATTYP
2
$ 1A 2O
0 100000
1 100000
10 100000
*CONTROL_ACCURACY
$1为打开双精度，4为实体壳体都支持；
$ 1OSU 2INN 3PIDOSU
1 4
*SET_NODE_LIST
$约束
$ 1SID 2DA1 3DA2 4DA3 5DA4
1
$ 1NID1 2NID2 3NID3 4NID4 5NID5 6NID6 7NID7 8NID8
1
$
*BOUNDARY_SPC_SET
$ 1NSID 2CID 3DOFX 4DOFY 5DOFZ 6DOFRX 7DOFRY 8DOFRZ
1 0 1 1 1 1 1 1
*LOAD_NODE_POINT
$节点载荷，第一个2为节点号2，如果用节点组，就是节点组号
$ 1ID 2DOF 3LCID 4SF 5CID 6M1 7M2 8M3
2 2 2 1.000 0
*END
后处理：
可以看位移
Beam下可以看st方向的受力、力矩情况、应力等；
History面板下可以看2号节点的位移曲线
