LS-DYNA常见问题集锦.doc

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号算法。该算法是针对单元翘曲而开发的算法，处理这种情况能够很好的保证求解的精度。 除了上述方法外，在计算时要注意控制沙漏，确保求解稳定。解决超大结果文件的方案:1. 将不同时间段内的结果分别写入一序列的结果记录文件；2. 使用/assign命令和重启动技术；3. ANSYS采用向指定结果记录文件追加当前计算结果数据方式使用/assign指定的文件，所以要求指定的结果记录文件都是新创建的文件，否则造成结果文件记录内容重复或混乱。特别是，反复运行相同分析命令流时，在重复运行命令流文件之前一定要删除以前生成的结果文件序列。具体操作方法和过程参见下列命令流文件的演示。问题： 怎样显示梁单元径向和轴向的应力分布图（我作的梁单元结果只有变形图DOF SOLUTIN –Translation，但是没有stress等值线图，只有一种颜色）和壳单元厚度方向的应力、变形图（我们只能显示一层应力、变形，不知道是上下表层或中间层的结果）。 解答：如果想显示梁单元的应力等值线图，请打开实际形状显示功能（PLotCtrl-Style-Size and Shape-/ESHAPE选为ON），然后即可绘制。注意梁单元（如BEAM188，BEAM189）的应力结果是在单元坐标系中显示的，即SXX为轴向正应力，SXY，SXZ为截面剪应力，没有其他应力分量。另外，缺省情况下，只输出SXX，如果想观察SXY，SXZ，请将BEAM188或189的KEYOPT（4）选为Include both（以这两个单元为例，其他单元可能不同，请看帮助文件，推荐使用BEAM188，BEAM189，这是功能最强的梁单元）。 至于壳的应力显示也类似，请打开实际形状显示功能，即可如同在实体上一样显示结果，您可以很清楚地看出不同位置、高度的应力值。当然如果你只想画出顶部、中部或底部的应力图也可以，以shell63为例，首先需关闭powergraphics（Toolbar上点POWRGRPH，选择OFF），然后进入General PostProc-Option for outp-SHELL中选择位置即可。 这是由于在建立模型时PART与PART之间有初始穿透，尤其是壳单元模型时很容易发生，应当避免这种情况的出现，否则容易在有初始穿透的地方产生塑性铰，原因是程序在求解的开始阶段给与穿透相应的接触力消除穿透，使材料发生局部塑性变形。解决方法见2002年11月的应用技巧。 /PREP7ET,1,63RECT,,10,,10ESHAPE,2ESIZE,,20AMESH,1EPLOMXNODE = NDINQR(0,14)*DIM,THICK,,MXNODE*DO,NODE,1,MXNODE*IF,NSEL(NODE),EQ,1,THENTHICK(node) = 0.5 + 0.2*NX(NODE) + 0.02*NY(NODE)**2*ENDIF*ENDDONODE = 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 EPLO 9不同自由度的单元在merge时，ANSYS如何处理自由度的不匹配问题 ANSYS可以将在给定容差（tolerance）内的节点通过merge可以合并到一起，它可以是不同的单元类型和不同的自由度（如：壳或梁（6自由度）与块单元（3自由度），ANSYS在处理共节点的自由度关系使，类似于藕合自由度（Couple DOFs）而非约束方程，即只把相同的自由度联系起来，这样就可