dyna第七章培训课件.ppt

《dyna第七章培训课件.ppt》由会员分享，可在线阅读，更多相关《dyna第七章培训课件.ppt（29页珍藏版）》请在三一办公上搜索。

1、求解和模拟控制,第7章,March 7, 2002Inventory #0016307-2,求解和模拟控制本章目的,本章学习的是求解和模拟控制。 讨论用于优化和控制显式分析的工具。主要内容:基本求解控制二进制输出文件 ASCII输出文件Cpu 时间控制质量缩放模拟控制小时间步显示自适应剖分编辑LS-DYNA输入文件 LS-DYNA 的批处理求解和模拟控制练习,March 7, 2002Inventory #0016307-3,在很多方面, 显式分析中设定的求解控制参数同隐式运算中遇到的参数相似。在显式求解中设定的基本参数有：1. 时间:模拟的物理过程的真实时间 TIMESolution: Ti

2、me Controls - Solution Time,实际的求解时间很短，通常在几毫秒内。,求解和模拟控制A. 基本求解控制,March 7, 2002Inventory #0016307-4,2.梁和壳单元的积分点数 EDINTSolution: Output Controls - Integ Pt Storage壳单元至少要有3-5个积分点才能捕捉到塑性效应3. 写入到 .rst 和 .his 文件的时间步数 EDRST, EDHTIME . 或者定义输出频率：Solution: Output Controls - File Output Freq -Number of Steps.,.

3、rst 文件记录了可供通用后处理器使用的整个模型的结果。一般典型的输出步数为 10-100 个 (缺省为100)来减少磁盘存储空间.his 文件记录了用于时间历程后处理器的模型的一个子组的结果。一般典型的输出步数为 1000-100,000个(缺省为 1000）重启动的输出频率也可以输入(通过EDDUMP 命令).,求解和模拟控制. 基本求解控制,March 7, 2002Inventory #0016307-5,文件选项包括 ADD,DELETE和LIST.可输出 ANSYS (.rst and .his) 或LS-TAURUS (d3thdt, d3plot) 或两者兼之。,求解和模拟控制

4、B. 二进制输出文件的控制,由于ANSYS/LS-DYNA 还带有 LS- POST 后处理器, 故程序除了可以写 ANSYS 的 .his 和 .rst 二进制结果文件以外，还可以写 LS-DYNA 结果文件 d3plot 和 d3thdt. 用EDOPT命令控制程序将输出何种二进制文件。Solution: Output Controls - Output File Type,March 7, 2002Inventory #0016307-6,求解和模拟控制C. ASCII 输出文件的控制,除了 LS-POST 二进制文件, 用户还可以输出一系列包含特定信息的 ASCII 文件.文件描叙如下

5、:GLSTAT全局模型数据 (全局统计) BNDOUT边界条件力与能量RWFORC刚性墙力DEFORC离散单元力MATSUM材料能量 (基于PART)NCFORC节点界面力RCFORC合成界面力DEFGEO变形几何数据SPCFORC单点约束力SWFORC节点约束反力 (点焊)RBDOUT刚体数据GCEOUT几何接触实体SLEOUT滑移面能量JNTFORC接点力NODOUT 节点数据ELOUT单元数据,March 7, 2002Inventory #0016307-7,EDOUT命令控制ASCII文件的输出：Solution: Output Controls - ASCII Output.,仅仅

6、选择结果文件应该生成哪个ASCII 文件(如 MATSUM)还可以选择下面的选项:ALL 写所有的 ASCII 输出文件LIST 列出所有的时间历程输出定义DELE 删除所有的 ASCII 输出定义,应该I 文件要写的节点与单元:应该注意到ACSII文件包含的信息只是针对模型的一个小子集而言，EDHIST命令控制ACSII文件要写的节点和单元：应该注意到 ASCII文件中包含的信息只是针对模型Solution: Output Controls - Select Component,包含在ASCII 输出文件中的信息是相对于节点和单元组元而定义的。ASCII文件的输出间隔通过 EDHTIME命令

7、控制。对于不同的组元允许多次定义。,求解和模拟控制. ASCII 输出文件的控制,March 7, 2002Inventory #0016307-8,ANSYS/LS-DYNA中, 经常使用三种高级求解控制1. CPU 控制: 在 ANSYS/LS-DYNA 分析中说明CPU限制。2. 质量放缩: 调节单元质量来增加时间步长。3. 子循环: 调整模型以减少CPU时间（不推荐）。,CPU 时间限制 (单位为秒) 用来终止分析，零值 (缺省) 表示无时间限制,CPU 控制l用EDCPU命令制定CPU时间：Solution: Analysis Options-CPU Limit,求解和模拟控制D.

8、CPU 时间控制,March 7, 2002Inventory #0016307-9,说明所用的最小时间步尺寸和缩放系数。,求解和模拟控制E. 质量缩放,质量缩放是通过调整每个单元的密度来调整单元达到合适的时间步长。以前讲过，时间步长依赖于 EX, NUXY 和 DENS 以及单元尺寸。使用质量缩放, 可以根据单元的大小调整任何单元的密度，从而达到合适的时间步。用 EDCTS 命令定义质量缩放:Solution: Time Controls-Time Step Ctrls.,March 7, 2002Inventory #0016307-10,求解和模拟控制. 质量缩放,有两种质量缩放方案：通

9、过调整单元密度，使得所有单元都具有相同的时间步长，只用于惯性效应不重要时.EDCTS, DTMS 这里DTMS为正的必要的时间步尺寸值 (乘以 1.111)质量缩放只用于 指定时步Dt的单元EDCTS, DTMS这里 DTMS 为负的必要的时间步尺寸值 (乘以 1.111),March 7, 2002Inventory #0016307-11,element 1 2 3,l1,l2,l3,求解和模拟控制. 质量缩放,质量缩放控制最小时间步长:,调整密度以实现用户定义的时间步长: :,March 7, 2002Inventory #0016307-12,求解和模拟控制. 质量缩放,质量缩放例题

10、(EDCTS)汽车碰撞模型140 parts42981 nodes1580 bricks60 beams35170 shellsTermination time 150 ms,100 个最小单元时步 (参见 LS-DYNA 输出文件 d3hsp): element time-step shell 151018 0.44612E-06 shell 150894 0.46867E-06 shell 52321 0.48682E-06 shell 51321 0.48682E-06 shell 16923 0.52225E-06 shell 16458 0.52225E-06. shell 1524

11、83 0.70112E-06 shell 92708 0.70113E-06 shell 92308 0.70114E-06 shell 38547 0.70223E-06 shell 38047 0.70223E-06,March 7, 2002Inventory #0016307-13,求解和模拟控制. 质量缩放,无质量缩放:初始时间步长取决于最小单元Dt = 4.46115E-07 sec有质量缩放:希望时间步长为 6.534E-07 秒，用EDCTS命令定义负的111%的时间步长。EDCTS, -7.26E-07Main Menu Preprocessor Output Ctrls M

12、ass scaling .初始时间步长Dt = 6.534E-07 secCPU 时间减少到68%质量误差：物理质量 1.26 吨增加质量 0.000027 吨 (27 grams)质量误差 0.002%质心坐标也发生改变。,March 7, 2002Inventory #0016307-14,求解和模拟控制F. 模拟控制,状态开关控制允许用户中断求解过程并检查求解状态。如下使用状态开关控制：在UNIX平台上的ANSYS输出窗口或在NT平台上独立的LS-DYNA输出窗口使用CTRL_C。 CTRL_C中断求解过程并在ANSYS窗口等待下一步输入。输入sw1终止求解并生成重启动文件。输入sw2得

13、到实际状态统计情况，LS-DYNA 继续进行。输入sw3写出当前时刻的重启动文件， LS-DYNA 继续运行。输入sw4写出当前时刻的结果数据组，LS-DYNA 继续运行。,March 7, 2002Inventory #0016307-15,求解和模拟控制. 模拟控制,通常第一次估计的CPU时间过高，用CTRL_C中断求解过程，输入sw2得到求解统计情况。,LS-DYNA 求解器 将所有的重要的信息（如错误、警告、失效单元、接触情况）都写入ANSYS输出窗口（NT上的独立窗口）及d3hsp文件。,March 7, 2002Inventory #0016307-16,求解和模拟控制. 模拟控制

14、,LS-DYNA 求解器也将数据写入messag 文件. 这个文件包含了时间步信息、错误、警告和一些写入d3hsp文件内容的一个简短摘要。 对于ANSYS检查出的错误和警告，用户可以通过 messag 文件了解详细情况 :,“Solution is done!” 代表没有警告和错误。这表示运算已经成功。,March 7, 2002Inventory #0016307-17,求解和模拟控制G. 小时间步长的查看,LS-DYNA 求解器可以根据每个单元的大小和密度自动计算出每个单元的最小时间步长。LS-DYNA 实际采用的时间步长是这些值中最小的值乘以一个系数(通常为0.9）。 EDTP 命令可以

15、用在LS-DYNA求解之前显示具有最小步长的单元。 这个信息允许用户评价网格质量，可进行适当的修改（如重新划分网格或运用质量缩放）。,EDTP, OPTION, VALUE1, VALUE2 其中:OPTION = 1, 2, or 3 描述如下:1 = 画出具有最小时间步长的单元2 = #1 + 此时间步值单元列表3 = #2 + 其它单元透明显示VALUE1 = 画/列出 最小描述的限制 (红色单元表示)VALUE2 = 透明 ( 0 = 不透明, 1 = max., 0.9 = default level),March 7, 2002Inventory #0016307-18,求解和模拟

16、控制. 小时间步长的目测,EDTP 在求解之前映射出时间步。模型中最小的单元将控制CPU的时间。具有最小时间步的单元以红色绘制（中间时间步单元为黄色) 具有透明选项和列出时间步大小的功能。 根据估计的时间步大小进行划分网格/质量缩放。Solution: Time ControlsTime Step Predic.,March 7, 2002Inventory #0016307-19,求解和模拟控制H. 自适应网格划分,在分析求解中，shell163 单元可以重新自动划分网格以保持变形误差均匀。 -自适应网格是基于PART来定义 (EDADAPT 命令)EDCADAPT 控制自适应网格生成的频率

17、、准则、起始和终止时间。每种网格具有独自的 POST1 ，求解过程中生成 (Jobname.rs01, Jobname.rs02, ) 和POST26 (Jobname.hi01, Jobname.hi02, ) 等多个后处理文件自适应剖分特别适用于冲压或是钣金加工应用中，这些应用大多产生较大的塑性变形。,March 7, 2002Inventory #0016307-20,求解和模拟控制. 自适应网格划分,EDADAPT, PART, Key ! 设置 Key=On 激活该part的自适应划分EDCADAPT, FREQ, TOL, OPT, MAXLVL, BTIME, DTIME, LC

18、ID, ADPSIZE, ADPASS, IREFLG, ADPENE, ADPTH, MAXELFREQ = 自适应剖分细划时间（实际时间）间隔TOL = 基于OPT=1(初始）或2（变形后）网格的自适应角度。 MAXLVL = 网格细划最大级别BTIME/DTIME = 模型中进行剖分的开始时间和终止时间LCID = 确定重新剖分时间间隔的数据曲线ADPSIZE = 基于单元边长的最小自适应单元尺寸ADPASS = 单向或双向自适应IREFLG = 均匀细化级别ADPENE = 当达到或穿透模具表面时开始自适应的标志ADPTH =开始自适应的绝对壳厚度值MAXEL =中止自适应的最大单元数

19、EDCADAPT 命令总体控制所有的进行自适应的Part。,March 7, 2002Inventory #0016307-21,求解和模拟控制 . 自适应网格划分,通过对发生大扭曲的壳单元进行网格重划, ANSYS/LS-DYNA 确保得到更精确的结果。定义自适应重划分为两步1.选择一个需要自适应划分的PART。,Solution: Analysis Options Adaptive Meshing Apply to Part,March 7, 2002Inventory #0016307-22,求解和模拟控制. 自适应网格划分,2. 指定合适的自适应剖分控制,Solution: Analy

20、sis Options Adaptive Meshing Global Settings,March 7, 2002Inventory #0016307-23,求解和模拟控制. 自适应网格划分,一旦定义了自适应重划，LS-DYNA则根据 EDCADAPT 中的参数在求解过程中重划网格。不能直接动画显示不同的结果文件，但可以用FILE 和/SEG 命令编辑宏来动画显示。,RefinedMesh,March 7, 2002Inventory #0016307-24,求解和模拟控制I. 编辑 LS-DYNA 输入文件,ANSYS/LS-DYNA支持大多数 LS-DYNA 功能，通过 GUI可以很容易

21、进入。但是，有一些其它的LS-DYNA功能,无法通过ANSYS GUI实现，例如 材料模型: : Fabric, unified creep, Drucker-Prager单元: Air bags, seat belts, explosives约束: 刚体局部坐标系虽然这些功能不能直接使用，但用户可以直接编辑LS-DYNA输入文件来实现这些功能 。但是, 这样就不能用ANSYS进行后处理，但可以用 LS-POST 来验证求解结果。 编辑 LS-DYNA 输入文件一般有5步。,March 7, 2002Inventory #0016307-25,求解和模拟控制. 编辑 LS-DYNA 输入文件,

22、STEP 1执行 EDWRITE 命令当 ANSYS/LS-DYNA完成所有的建模后, 用 EDWRITE 命令来代替SOLVE 命令 。Solution: Write Jobname.K,EDWRITE 命令写出 Jobname.K 文件但不启动 ANSYS/LS-DYNA 求解过程。用 EDWRITE 命令还可以创建 .rst 和 .his 文件文件头。 EDWRITE 是一执行命令，可立即写出指定的 LS-DYNA 输入文件。,March 7, 2002Inventory #0016307-26,求解和模拟控制. 编辑 LS-DYNA 输入文件,STEP 2: 退出 ANSYS/LS-D

23、YNA 程序Utility Menu: File-Exit,STEP 3: 编辑输入文件 Jobname.k用标准的vi或文本编辑器,根据 LS-DYNA所需的关键字格式增加信息到 Jobname.K 文件中。,March 7, 2002Inventory #0016307-27,求解和模拟控制. 编辑 LS-DYNA 输入文件,STEP 4: 运行 LS-DYNA 求解器在 Jobname.RST 和 Jobname.HIS 所在的目录中直接执行 LS-DYNA 输入文件（它将直接激活求解器）。 LS-DYNA运行的结果直接追加到ANSYS结果文件中执行的LS-DYNA求解器的命令是:/an

24、sys60/bin/ls960 i=Jobname.K pr= ANSYS product variable (see license guide)此外， m=drelax 用于显-隐式求解.同时, MEMORY=# (单位：words) 可扩展内存 (see EDSTART). 还有, R=D3dumpnn 用于重启动。STEP 5: 重新进入 ANSYS/LS-DYNA 程序在LS-DYNA求解后，重新进入 ANSYS/LS-DYNA 程序 。可以在通用后处理器（jobname.rst)和时间历程后处理器(jobname.his)中查看结果NOTE: 在此过程中不能改变节点和单元，并且不能

25、编辑Jobname.K 文件 。,March 7, 2002Inventory #0016307-28,求解和模拟控制J. LS-DYNA 批处理,在ANSYS/LS-DYNA6.0 版本中有一个LS-DYNA launcher 的测试 版，称为 windyna.exe目录为 ansys60binintel directory 与ansys.exe 同一路径。在Pc和UNIX系统中的 ANSYS/LS-DYNA 6.1都支持该launcher。 LS-DYNA launcher 使得递交.K输入文件到LS-DYNA求解器中变得更容易。,March 7, 2002Inventory #0016307-29,求解和模拟控制K. 练习题,练习题内容：练习7： 受轴向载荷的梁的屈服分析请参考习题集,