《LSDYNA中文教程.docx》由会员分享，可在线阅读，更多相关《LSDYNA中文教程.docx（15页珍藏版）》请在三一办公上搜索。
1、LSDYNA中文教程第二部分 ANSYS/LS- DYNA程序的使用方法 1 概述 ANSYS/LS-DYNA 程序系统是将非线性动力分析程序LS-DYNA显式积分部分与ANSYS程序的前处理PREP7和后处理POST1、POST26连接成一体。这样既能充分运用LS-DYNA程序强大的非线性动力分析功能，又能很好地利用ANSYS程序完善的前后处理功能来建立有限元模型与观察计算结果，它们之间的关系如下。 ANSYS前处理PREP7 Jobname.DB Jobname.k ANSYS/LS-DYNA Jobname.RST d3plot Jobname.HIS d3thdt ANSYS/LS-D
2、YNA程序系统的求解步骤为： ANSYS后处理POST1,POST26 后处理LS-TAURUS 1.1 前处理Preprocessor建模 1.设置Preference(Main Menu:Preference)选项 置Structural LS-DYNA explicit。这样，以后显示的菜单完全被过滤成ANSYS/LS-DYNA的输入选项。再定义一种显式单元类型，即可激活LS-DYNA求解。 GUI: Main MenuPreferences a.选择Structural. b.选择LS-DYNA Explicit. c.OK. 2.定义单元类型Element Type和Option和实
3、常数Real Constant。 3.定义材料性质Material Properties。 4.建立结构实体模型Modeling。 5.进行有限元网格剖分Meshing。 6.定义接触界面Contact。 1.2 加载和求解Solution 1.约束、加载和给定初始速度。 2.设置求解过程的控制参数。 3.选择输出文件和输出时间间隔。 4.求解Solve。 1.3 后处理POST1和POST26，也可连接LSTC公司的后处理程序LS-TAURUS。 在各程序模块之间传递数据的文件有： (1) ANSYS数据文件 数据库文件(Database File)Jobname.DB 二进制文件 图形数据
4、文件(Results File)Jobname.RST 二进制文件 时间历程数据文件 Jobname.HIS 二进制文件 输出文件(Output File) Jobname.OUT ASCII文件 命令文件(Log File) Jobname.LOG ASCII文件 (2) LS-DYNA数据文件 输入数据文件(Iuput file) Jobname.K ASCII文件 重起动文件(Dump File) D3DUMP 随机文件 图形数据文件(Plot File) D3PLOT 随机文件 时间历程文件 D3THDT 随机文件 由于ANSYS前处理程序还不能满足LS-DYNA程序系统的全部功能，用
5、户可以生成LS-DYNA的输入数据文件Jobname.K，经过编辑、修改后，再直接调用LS-DYNA程序求解，其计算结果图形数据文件仍然可以连接ANSYS后处理程序POST1和POST26以及LS-DYNA的后处理程序LS-TAURUS观察计算结果。 ANSYS/LS-DYNA可以与ANSYS结构分析程序之间传递几何数据和结果数据来执行隐式显式和显式隐式分析，如跌落试验、回弹计算等。 使用ANSYS/LS-DYNA时，建议用户使用程序提供的缺省设置。多数情况，这些设置适合求解问题。 下面详细介绍ANSYS/LS-DYNA程序的使用方法，由于ANSYS前处理程序PREP7的实体建模和网格剖分在A
6、NSYS培训教材中详细介绍，这里不再重复。 2 单元 160-168九种单元 2.1 概述 在ANSYS/LS-DYNA程序的显式动态分析中可以使用如下类型单元： LINK160 BEAM161 杆单元 轴向力 梁单元 两个端点 PLANE162 二维实体单元 SHELL163 SOLID164 薄壳单元 薄壳 三维实体单元 三维一般情况 弹簧阻尼单元 弹簧 质量单元 一个点 COMBI165 MASS166 LINK167 仅位伸单元 索 SOLID168 四面体实体单元 这些单元都是三维的，采用线性位移函数，每种单元都有多种算法可供用户选择。其高斯积分点数，缺省时实体单元是单点积分，薄壳单
7、元是面内单点积分、沿壳厚多点积分。线性位移函数和单积分点的显式单元能够很好地用于大变形和材料失效等非线性问题。近年来新版本LS-DYNA的单元也可以采用多点积分，供用户选择。各种单元能够选用的材料模型是有限制的，详见下一节“材料性质”。 (1) SOLID164实体单元 它是8节点六面体单元，可以采用某些节点重复来退化成6节点楔形单元，5节点锥形单元。在受弯曲时退化单元比较僵硬，要尽量避免使用。实体单元算法有二种：(1)KEYOPT(1)=1，缺省算法，采用单点积分和沙漏控制。它能节省机时并在大变形条件下增加可靠性。(2)KEYOPT(1)=2，采用2X2X2多点高斯积分。它没有零能模式，不需
8、要沙漏控制。对于某种材料，如泡沫材料，计算结果会好些。但计算机时将有显著增加。 (2) SHELL163薄壳单元 它是4节点四边形单元，或3节点三角形单元。有11种算法，用KEYOPT(1)值来定义不同算法。和实体单元一样，积分点数会显著影响机时耗费，对于一般课题建议采用面内单点积分。以下概述这十一种算法。 a. 4节点四边形薄壳单元 Belytschko-Tsay(KEYOPT(1)=2)缺省算法 采用面内单点积分，计算速度快，建议在大多数分析中使用，单元过度翘曲时不要使用。 Belytschko-Wong-Chiang(KEOPT(1)=10) 比Belytschko-Tsay算法慢1/4
9、，采用面内单点积分，用于翘曲情况，一般可得到正确结果。 Belytschko-Leriathan(KEYOPT(1)=8) 采用面内单点积分，比Belytschko-Tsay算法慢2/5，自动含物理上的沙漏控制。 Hughes-Liu(KEYOPT(1)=1) 采用面内单点积分，比Belytschko-Tsay算法慢250%。 Fast(Co-Rotational) Hughes-Liu 采用面内单点积分，比Belytschko-Tsay算法慢150%。 S/R Hughes-Liu(KEYOPT(1)=6) 采用面内X积分点，没有沙漏，比Belytschko-Tsay算法慢20倍，如果在分析
10、中遇到沙漏麻烦的话，建议使用这种算法。 S/R co-rotational Hughes-Liu(KEYOPT(1)=7) 采用面内X积分点，没有沙漏，比Belytschko-Tsay算法慢8.8倍，如果在分析中遇到沙漏麻烦的话，建议使用这种算法。 b. 薄膜单元算法 Belytschko-Tsay(KEYOPT(1)=5)薄膜单元 采用单个积分点，计算速度快，建议在大多数薄膜分析中使用，可很好地用于纤维织品。 Fully integrated Belytschko-Tsay(KEYOPT(1)=9)薄膜单元 采用2X2积分点，无沙漏控制，比Belytschko-Tsay单积分点薄膜单元要显著
11、地慢。 c. 3节点三角形薄壳单元 C0三角形壳单元 面内单点积分，根据Mindlin-Reissner薄板理论导出，相当僵硬，建议不要用它做整体网格剖分。 BCIZ三角形壳单元 采用面内单点积分，根据Kirchhoff薄板理论导出，比C0三角形壳单元的计算速度慢。 退化的四边形单元在横向剪切时会发生锁死，应改用三角形壳单元。在同一种材料中，只要把单元分类标记置1，就可混合使用四边形单元和三角形单元。 所有壳单元算法，沿壳厚方向用户可任选2-5个高斯积分点，对于弹性材料沿壳厚方向2个积分点已足够，但对于塑性材料，至少要3个或更多的积分点。 薄壳单元的输出应力，尽管后处理中采用顶面Top和底面B
12、ottom术语，实际上它们是最外面积分点处的应力，而不是薄壳表面上的应力，因此，在分析计算结果时应特别注意。 (3) 梁单元和杆单元 a. BEAM161梁单元 采用单元端部二个节点定义，并以第3个节点对单元主轴面定向，它有二种算法： Hughes-Liu梁单元缺省算法 它是一个很方便的退化单元，可以用梁单元中间跨度横截面上的一组积分点来模拟矩形或圆形横截面。用户还可以自定义一个横截面积分法则来模拟任意截面形状，梁单元与实体单元、梁单元一样，网格必须足够细分以保证计算精度。由于积分点位于单元跨度中间横截面上，计算时只在单元中心检验材料屈服，沿单元长度为不变力矩。因此，必须在单元中间截面处产生全
13、塑性力矩，而不是在单元节点处，悬臂梁模型将在一个稍高的外力作用下才会屈服。 Belytschko-Schwer梁单元 单元的内力矩沿其长度方向线性分布。有“正确”的弹性应力和在单元端部检验屈服。与Hughes-Liu梁单元一样，单元质量集中在节点上，而在动态分析问题中，质量分布正确很重要，必须细分网格。 b. LINK160杆单元 它与Belytschko-Schwer梁单元很相似，但只能承受轴向加载，而不能承受弯曲力矩，用于桁架系统。适用的材料模型有弹性和随动硬化弹塑性。 c. LINK167索单元 用单元端部二个节点定义，是仅能拉伸的杆单元，可以模拟索，它与弹簧单元类似。由用户直接输入力与
14、变形的关系式。 (4) 离散单元和质量单元 a. COMBI165弹簧阻尼单元 用单元端部二个节点定义。 弹簧单元在改变单元长度时产生沿单元轴向的力。如：正力时在节点1处沿轴的正方向，而对节点2是沿轴的负方向。缺省时，单元轴就是从节点1到节点2。当单元旋转时，力作用线的方向也将随之旋转，可以模拟弹性、弹塑性和非线性弹性的弹簧性质。 阻尼单元可以模拟线性粘性阻尼和非线性粘性阻尼。 扭转弹簧阻尼单元也可使用，它由KEYOPT(1)选项来选择。它的力位移关系可以认为是力矩转角关系。旋转弹簧单元只影响它们节点的转动自由度。 COMBI165单元也可以与其它显式单元混合使用，由于COMBI165单元没有
15、质量，在分析中一个节点不能只与COMBI165一种类型单元连接，必须在该节点上加质量单元MASS166。 b. MASS166质量单元 质量单元由一个节点和一个质量值定义。采用质量单元可以简化部分结构，以减少动力分析所需单元数目，因而减少求解机时。 每个实体单元、壳单元、梁单元、杆单元的质量都平均分配到单元的节点上。在壳单元和梁单元情况下，每个节点上还将附加一个转动惯量：只采用一个单值，它的作用就是让质量围绕节点球形分布。 2.2 定义单元类型和实常数 在ANSYS/LS-DYNA的前处理时，通常采用图形用户界面操作来选择和输入各种参数，也可采用命令输入。 用ET命令定义单元类型 用R命令来定
16、义实常数 具体执行过程如下： 第一步：定义单元类型 GUI：Main Menu：Preprocessorelement TypeAdd/Edit/Delete 1. Add一个单元类型 2. 选LS-DYNA Explicit单元族 3. 采用薄壳单元SHELL163 4. OK接受这个选择 指定S/R Co-rotational Hughes-Liu面内多积分点改进型单元算法，以消除某种沙漏模态。 5. 选择Options 6. 选择S/R Co-rotational算法 7. OK 8. Close 第二步：定义实常数 GUI：Main Menu：PreprocessorReal Cons
17、tants 1. Add 2. OK，定义SHELL163单元的Real Constant Set 3. 如果单元是均匀壳厚，那么只要定义节点1处的壳厚，本例输入T1=0.1。剪切因子SHRF和沿壳厚积分点数NIP取缺省值。NIP的缺省值为3。 4. OK 5. Close 3 材料 LS-DYNA程序可以使用近100种材料模型，但在ANSYS/LS-DYNA前处理PREP7中直接输入的材料模型约30种，其他材料模型可以修改LS-DYNA程序输入数据文件Jobname.K来加入。ANSYS/LS-DYNA程序可以直接输入的材料模型为： 3.1 线弹性模型 材料类型1：各向同性O, 1H, 1B
18、, 1I, 1T, 2, 3 材料类型2：正交各向异性0, 2, 3 材料类型2：各向异性0 3.2 非线性弹性模型 材料类型7：Blatz-Ko橡胶0, 2 材料类型27：Mooney-Rivlin橡胶0, 2 材料类型6：粘弹性0, 1H 3.3 弹塑性模型 材料类型3：双线性随动硬化0, 1H, 1I, 1T, 2, 3，与应变率无关。 材料类型3：双线性各向同性硬化0, 1H, 1I, 1T, 2, 3 ，与应变率无关。 材料类型3：随动塑性0, 1H, 1I, 1T, 2, 3 ，与应变率相关、失效各向同性、随动和组合硬化。 材料类型18：幂指数硬化塑性0, 1H, 2, 3 ，与应
19、变率相关。 材料类型19：与应变率相关各向同性硬化0, 2, 3 材料类型64：与应变率相关的幂指数硬化塑性0, 2, 3 材料类型36：Barlat和Lian发展的平面应力状态各向异性弹塑性，用于壳单元2 材料类型33：由Barlat、Lege和BremXX年发展，用于模压过程的各向异性弹塑性0, 2, 3 材料类型24：分段线性弹塑性0, 1H, 2, 3 材料类型37：横向各向异性弹塑性 ，用于板压成型2, 3 材料类型10：弹塑性流体动力，0 1. 泡沫模型 材料类型53：低密度、闭合的多孔聚氨酯泡沫0 材料类型57：低密度氨基 酸 酯泡沫0 材料类型62：粘性泡沫0 材料类型63：可
20、压扁泡沫0 材料类型26：正交异性可压扁Honeycomb蜂窝结构 2. 复合材料 材料类型22：考虑损伤的复合材料0, 2, 3 3. 混凝土模型 材料类型72：考虑损伤的混凝土0 4. 需要状态方程的材料模型 材料类型15：John/Cook与应变率、温度相关的塑性0, 2 材料类型9：无偏应力流体动力模型0 材料类型65：Zerilli-Armstrong模型，用于金属成型过程和高速碰撞，其应力与应变、应变率及温度有关，0, 2 材料类型51：Bamman模型，用于金属成型过程，其塑性与应变率和温度相关，0，2，3 这四种材料类型，均需附加状态方程，可选用(1)Linear Polyno
21、mial，(2)Gruneisen和(3)Tabulated状态方程。 5. 其他 材料类型20：刚性材料0, 1H, 1B, 1T, 2, 3 材料类型71：索，1D 材料类型1：弹性流体,0 在方括号 中的数字符号标志该材料模型可以用于何种单元类型： 0 实体单元 1H Hughes-Liu梁单元 1B Belytschko Resultant梁单元 1I Belytschko Integrated Solid and Tubular Beams 1T 杆单元 1D 离散单元 薄壳单元 厚壳单元 用图形用户界面输入材料数据的过程： 第一步：增加一个新的材料模型 特殊的气囊单元 Prepro
22、cessor: Material PropsDefine MAT ModelAdd 第二步：选择所需材料模型 检查输入材料号是否正确 先选择材料模型所属的材料族 再选择所需材料模型 选择OK键 第三步：输入所需值 几乎所有材料模型，都需输入 DENS, EX和NUXY 材料应力应变数据应该是真 应力真应变格式 待所有数据已经输入后，选择 OK键 下面介绍各种材料模型的数据输入方法： 第一类 线弹性材料 这一材料族有三种材料模型： Elastic，Orthotropic和Anisotropic。 1. Elastic(Isotropic)各向同性弹性 多数工程材料都是各向同性材料，在弹性变形范围
23、内采用此材料模型，只要输入DENS、EX和NUXY。缺省值EX=30E6, NUXY=0.3。 2. Orthotropic正交各向异性弹性 通用正交各向异性弹性材料由9个独立常数和质量密度定义。 横向正交各向异性弹性材料由5个独立常数定义。 正交各向异性弹性材料定义与一个给定的坐标系统号有关。 (1) 选择Linear ElasticOrthotropic (2) 输入所需材料常数和质量密度 (3) 输入Coordinate System ID，它由EDLCS命令定义，或 GUI：PreprocessorLS-DYNA Options-Local CSCreate Local CS 3. A
24、nisotropic各向异性弹性 各向异性弹性材料由21个独立的材料常数和质量密度定义。仅可用于实体单元。各个材料常数在本构矩阵中位置如下式所示。 (1) 选择Linear ElasticAnisotropic (2) 输入所需材料常数和质量密度 (3) 输入由EDLCS命令定义的Coordinate System ID 第二类 非线性弹性材料 这一材料族有三种材料模型：Platz-Ko, Mooney-Rivlin和Viscoelastic。 1. Blatz-Ko接近不可压缩橡胶 Blatz-Ko橡胶是由Blatz和Ko定义的超弹性橡胶模型，该模型使用第二类Piola-Kirchoff应力
25、。 11Sij=GCij-V1-2ndijV G剪切模量；V相对体积；n泊松比；Cij右柯西格林应变张量；dijkronecker delta。 当剪切模量作为仅有的材料性质定义时就能运行这种材料选项。泊松比NUXY自动取0.463。 2.Mooney-Rivlin不可压缩橡胶 Mooney-Rivlin橡胶是不可压缩橡胶模型。它与ANSYS的Mooney-Rivlin二参数模型很相似。 通过C10、C01和n来定义应变能量密度函数： W=C10(l1-3)+C01(l2-3)+C(C=C10+C012 1l32-1)+D(l3-1)2D=C10(5n-2)+C01(11n-5)2(1-2n)
26、 输入数据DENS、NUXY和Mooney-Rivlin常数C10与01。 3.粘弹性 Viscoelastic 由Herrmann和Peterson提出的线粘弹性材料模型。该模型采用偏量特性： (t)eijsij=2f(t-t)dt0t t式中剪切松弛模量为： f(t)=G+(G0-G)e-bt 在模型中，由弹性体积模量K和体积V计算压力：P=KlnV；需要参数G、G0、K和b来定义线粘弹性材料模型。 输入数据： G0短期限的弹性剪切模量 G长期限的弹性剪切模量 K弹性体积模量 Yb衰减常数 DENS质量密度 第三类 弹塑性材料 这一材料族有十种材料模型，其中2种是与应变率无关的各向同性材料模型，5种是与应变率相关的各向同性材料模型和3种是与应变率相关的各向异性材料模型。 (A) 与应变率无关的各向同性材料 二种： a.双线性随动硬化Bilinear Kinematic b.双线性各向同性硬化Bilinear Isotropic 这二种材料模型的应力应变曲线用二条直线表达，其不同在于硬化假定，随动硬化或各向同性硬化，这类与应变率无关的材料模型，可用于多种工程金属材料。 与应变率相关的各向同性材料
