AUTODYN培训二学习课程

1、2021-11-20基础培训二1、Lagrange求解器2、模型生成3、Lagrange之间的作用 第1页/共114页第一页，编辑于星期五：十一点 三十一分。2021-11-20AUTODYN 求解器类型LagrangeEulerALESPHShellBeam第2页/共114页第二页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 节点与材料一起运动 通常用于固体求解T=0.0T0.0第3页/共114页第三页，编辑于星期五：十一点 三十一分。2021-11-20 在初始计算中，通过对节点速度积分得到新节点位置； 通过新节点的坐标位置，计算新单元密度和单元应变率。Dt拉格朗日

2、求解器计 算 周 期第4页/共114页第四页，编辑于星期五：十一点 三十一分。2021-11-20 通过应变，新的应力可以计算得到； 对于各向同性材料，单元变形能分成两个独立的部分。体积改变形状改变 拉格朗日求解器计 算 周 期第5页/共114页第五页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 体积改变是由球应力张量引起； 形状改变是由于偏应力张量引起； 球应力由状态方程(EOS)求解； 偏应力是由胡克定律(弹性)和塑性屈服准则(本构关系)求解-即AUTODYN的强度模型； 再加上失效模型就构成了一个完整的材料模型。计 算 周 期第6页/共114页第六页，编辑于星期五

3、：十一点 三十一分。2021-11-20 一旦计算单元的内部应变率被决定，每一个节点的受力能被计算； 在这个阶段，任何从边界条件或者与其它物体作用的外力作为最后的节点力； 知道节点力，可以计算节点加速度； 积分节点加速度可以得到节点新的速度； 完整的一个计算周期，重复这样的周期达到要求的求解时间或循环次数。拉格朗日求解器计 算 周 期第7页/共114页第七页，编辑于星期五：十一点 三十一分。2021-11-20质量守恒能量守恒状态方程动量守恒积分积分 S S 力 S S 质量初始条件节点位移应变率，密度压力，内能 偏应力节点力节点加速度节点速度外力 (边界条件和相互作用力)拉格朗日求解器计 算

4、 周 期第8页/共114页第八页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器低 速挤 压 碰 撞第9页/共114页第九页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 结构求解器(V5和更早版本)； 2D 体单元和壳单元3D 体单元、壳单元和杆单元 非结构求解器(V6和以后版本)； 3D 体单元、壳单元和杆单元 结构和非结构求解器能用在同一个模型中； 非结构化求解器特点：求解更快 内存更小适合更多的单元类型 可以方便地输入非结构网格 结构化网格能被转换成非结构化网格。结构和非结构求解器第10页/共114页第十页，编辑于星期五：十一点 三十一分。20

5、21-11-20拉格朗日体单元求解器非结构求解器单元和面连通性数据存储满足精确要求结构求解器基于IJK 空间单元和面通过 IJK 空间定义 (与数据存储一样)数据存储直接与单元拓扑相连第11页/共114页第十一页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器 8 节点精确体积积分； 可选的 8 节点，近似单点高斯体积积分适合非结构求解器；FE 公式 (Hallquist) 对于弯曲变形，单元精度会有一些损失 使用人工粘性处理冲击问题； 使用沙漏控制保证稳定性； 五面体单元 退化六面体单元； 非结构六面体求解速度约为结构求解的两倍； 非结构求解需要的内存比结构求解少3

6、0%。六面体/五面体求解器第12页/共114页第十二页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器比较：内存(Mb) 时间(Min)加速比结构4001.741.00非结构(Exact)3400.881.97非结构(Gauss)3400.622.803.2Ghz Pentium IV Processor, Windows XP500,000 单元效率 内存 速度六面体/五面体求解器第13页/共114页第十三页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器六面体/五面体求解器结构六面体 (IJK)非结构六面体 (Exact)第14页/共11

7、4页第十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器六面体/五面体求解器Structured (IJK)UnstructuredExactUnstructuredGauss第15页/共114页第十五页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器头 盔 碰 撞第16页/共114页第十六页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器 结构 非结构 结构 非结构头 盔 碰 撞运行时间结构求解器：443分钟(36284 cycles)非结构求解器：265分钟(37013 cycles) 第17页/共114页第

8、十七页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器非结构四面体求解器SCP Tet ANP Tetstandard constant pressureAverage Nodal Pressure第18页/共114页第十八页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日体单元求解器非结构四面体求解器SCP Tet ANP Tet第19页/共114页第十九页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器 2D 采用有限差分方法； 3D 采用有限体积方法； 节点轨迹变形； 厚度是一个填充参数； 假设厚度方向为零应变； 时间步长

9、由长/宽决定，不由厚度决定。第20页/共114页第二十页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器 非结构求解比结构求解器更有效 在相同的内存条件下，非结构求解器存放的数据多于结构求解器四边形单元内存内存(Mbytes)时间时间(分钟分钟)加速比加速比结构结构5903.131.00非结构非结构5951.162.703.2Ghz Pentium IV Processor, Windows XP第21页/共114页第二十一页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器四边形壳单元结构 VS 非结构第22页/共114页第二十二页，编辑于星

10、期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器结构非结构金属成型运行时间(65,000 cycles)结构网格 (IJK)3,814 分钟非结构网格1,722分钟第23页/共114页第二十三页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器非结构主四边形壳单元金属冲压成型第24页/共114页第二十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳单元求解器四边形/三角形 壳单元非结构 (上面)、结构四边形 (中间)、结构三角形 (底部)位移时间第25页/共114页第二十五页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗

11、日壳单元求解器 用户定义层数； 每一层可以是各向同性或异性材料； 每一层指定一个厚度。多层复合壳单元外层内层第26页/共114页第二十六页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日梁单元求解器 2节点的 梁单元公式； 允许大的轴向应变。内存内存(Mbytes)时间时间(分钟分钟)加速比加速比IJK梁梁1500.441.00非结非结构梁构梁1650.251.763.2Ghz Pentium IVProcessor, Windows XP结构/非结构梁单元求解器效率 (100,000 梁单元)例子：第27页/共114页第二十七页，编辑于星期五：十一点 三十一分。2021-11-

12、20拉格朗日梁单元求解器梁单元截面形状第28页/共114页第二十八页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日梁求解器弹簧单元第29页/共114页第二十九页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日梁求解器弹簧单元速度与时间曲线第30页/共114页第三十页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日梁求解器阻尼单元第31页/共114页第三十一页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日梁求解器弹簧 & 阻尼单元联合速度与时间曲线第32页/共114页第三十二页，编辑于星期五：十一点 三十一分。2021-11-2

13、0拉格朗日梁单元 单根梁 多根梁 2D 梁构架 3D 梁构架 梁求解器集合组名 (类型，梁号单元/梁) 第33页/共114页第三十三页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日壳 & 梁求解器冲击波对结构中的影响第34页/共114页第三十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器(3D) 兼容性： 模型 非结构和结构(IJK)求解器能在一个模型中使用； 接触 在任意非结构/结构(IJK)之间使用； 连接 在任意非结构/结构(IJK)之间使用； 与3D FCT 和多物质欧拉耦合 包括厚度壳单元； 并行 可用于结构求解器和非结构求解器。第3

14、5页/共114页第三十五页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 优势 ： 更少的循环计算时间 (与 Euler/SPH 比较)； 清晰的材料界面和边界定义； 时间历史信息完整； 强度模型好； 较少的冲击扩散 ； 代码相对简单 。 壳单元优势 ： 时间步长由单元面积决定， 不是单元厚度。 梁单元优势： 时间步长由单元长度决定 ， 不是横截面积。不足 ： 单元变形导致时间步长变小； 单元变形导致网格混乱； 特殊的需要碰撞 / 滑移界面。第36页/共114页第三十六页，编辑于星期五：十一点 三十一分。2021-11-20基础培训二1、Lagrange求解器2、模型生成

15、3、Lagrange之间的作用 第37页/共114页第三十七页，编辑于星期五：十一点 三十一分。2021-11-20模型生成 结构化网格 用 (I,J,K) 指标空间； 网格在 AUTODYN 里生成； 通过 ICEM-CFD 或 TrueGrid 输入结构化网格； 完全的 Fill 使用选项 (用材料 / 初始条件) ； 边界条件用指标空间，交互式或者组施加 。 非结构化网格 使用单元/节点连通性； 目前，不能在 AUTODYN 中生成网格； 通过 LS-DYNA .k 文件和 Nastran .dat 文件输入非结构网格； 转换结构化网格为非结构化网格； 有限的使用Fill 选项； 边界条

16、件用交互式或者组施加。第38页/共114页第三十八页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 对于每一个 Part，一个结构化网格通过单元、节点和空间位置定义连通性； 网格可以看成两中空间： 指标空间（I 、J、K） 物理空间（X、Y、Z） 指标空间在节点之间定义连通性； 物理空间定义节点的空间坐标； 物理空间通常用来识别； 指标空间用来建立模型。第39页/共114页第三十九页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 指标空间 在AUTODYN里，每一个结构 Part 定义一个指标空间(i,j,k)，i,j,k是从 1 到 N 的的整数 指标空间

17、是矩形(2D)或立方体(3D)第40页/共114页第四十页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 物理空间 在AUTODYN里，每一个结构 Part 也定义一个 xyz 空间，x,y,z是实数 在物理空间里，通常网格有一个形状i = 1 第41页/共114页第四十一页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 不用区域 不是所有在指标空间定义的单元都需要定义物理空间-这就可以实现复杂几何模型划分网格不用的单元连接的节点指标空间 物理空间 第42页/共114页第四十二页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 Part 向导

18、 (predef) 允许自动划分预先确定几何的结构化网格： 2D Volume: Box，Quad，Circle，尖顶拱形 ，菱形，三角形，楔形 3D Volume: Box, Hex, Cylinder, Sphere, 尖顶拱形, Fragments/Bricks Shells: Plate, Cylinder Beams: Single, Multiple, 2D Frame, 3D Frame 网格转换 AUTODYN里创建结构化网格输入结构化网格：ICEM-CFDTrueGrid手工生成结构网格：点边面体第43页/共114页第四十三页，编辑于星期五：十一点 三十一分。2021-11-

19、20结构求解器生成 Part 向导三部曲定义几何定义网格填充材料第44页/共114页第四十四页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导-Box (2D)第45页/共114页第四十五页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导-Circle (2D)第46页/共114页第四十六页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导-2D楔形尖顶拱形四边形菱形三角形第47页/共114页第四十七页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导 - Box (3D)第48

20、页/共114页第四十八页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导 - Cylinder (3D)第49页/共114页第四十九页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器Part 向导 - 3D球六面体尖顶拱形砖块结构第50页/共114页第五十页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 一步步的生成网格： 点 线 面 体 拉伸 手工生成网格第51页/共114页第五十一页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 生成节点： 节点的 I、J 和K的值； 节点坐标。手工生成网格第52页/共114

21、页第五十二页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 生成线条： I/J/K方向； 几何比例； 直线或圆弧； 圆弧中心。手工生成网格第53页/共114页第五十三页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 生成平面： I/J/K常数； 网格类型。手工生成网格第54页/共114页第五十四页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器手工生成网格 生成体网格： 生成方向；网格比例。第55页/共114页第五十五页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器网格过渡第56页/共114页第五十六页，编辑于星期五：十一点

22、 三十一分。2021-11-20结构求解器Part 向导 - 3D圆柱面平面第57页/共114页第五十七页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 TrueGrid (XYZ Scientific) 在 TrueGrid 生成网格 输出AUTODYN-3D网格数据 (.zon) 用 “Import, TrueGrid” 下拉菜单分别读入 Part 能够定义块和壳单元 自动创建和填充多物质 Part 输入的体单元能够用于拉格朗日、ALE 和 Fill 求解器输入结构网格第58页/共114页第五十八页，编辑于星期五：十一点 三十一分。2021-11-20结构求解器 ANSY

23、S, ICEM-CFD 超强的六面体网格生成工具 直接连接 CAD CATIA，Pro/E，SDRC I-DEAS，SolidWorks，UG， 提供输入结构网格到 AUTODYN ICEM multiblock meshes (.geo 文件) 输入过程输入结构网格第59页/共114页第五十九页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 除了初始条件，有的问题还需要边界条件。 AUTODYN 中边界条件包括： 应力 速度 力 传输 流动 (Euler) 通过定义边界条件，然后通过指标空间施加，对于非结构 网格通过 Groups来施加。边界条件第60页/共114页第六

24、十页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 对于拉格朗日 Part 边界上的单元，缺省的边界条件是自由边界 (压强 =0.0)； 不用的单元也定义边界条件。拉格朗日边界条件不用单元填充单元缺省边界条件第61页/共114页第六十一页，编辑于星期五：十一点 三十一分。2021-11-20 常数t梯形t 三角形t 指数Pk = Pe-ktt用户子程序 EXSTRt拉格朗日 Part拉格朗日求解器应力边界条件(Lagrange/ALE/Shells)第62页/共114页第六十二页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 X-Velocity 常

25、数； Y-Velocity 分段表示； Z-Velocity 限制在最大和最小之间； 通过子程序 EXVEL。 R-Velocity 常数、指定中心或 X/Y/Z 轴 General Velocity在X、Y、Z (3D) 上固定常数和固定绕坐标 轴的角速度速度边界条件(Lagrange/ALE/Shells/Beams/SPH)第63页/共114页第六十三页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 节点力边界条件： 常数、分段、子程序 EXFOR 定义； X、Y、Z 和通常的方向。 单元力/单位长度边界条件： 常数、分段、子程序 EXFOR 定义； X、Y、Z 和

26、通常的方向。力边界条件(Beams)第64页/共114页第六十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日求解器 传输波通过单元面； 仅仅为正交分量传输； 可以指定边界的阻抗： 如果设置阻抗为零，使用临近单元的阻抗； 如果边界条件是近似的，不要在关心的区域设置。 传输边界条件(Lagrange/ALE/Euler)第65页/共114页第六十五页，编辑于星期五：十一点 三十一分。2021-11-20时间步长控制时间步长的稳定性0 T 0 T X/C R = C.T拉格朗日欧拉其中x = 单元尺寸，c = 声速，v = 速度第66页/共114页第六十六页，编辑于星期五：十一点

27、 三十一分。2021-11-20时间步长控制 缺省通常 OK； 如果初始步长小于零，时间步长取稳定时间步长的1/2； 如果最小时间步长小于零，时间步长取出是时间步长的1/10； 对于拉格朗日计算问题，安全系数可以增大到0.9 。时间步长选项第67页/共114页第六十七页，编辑于星期五：十一点 三十一分。2021-11-20阻尼选项 人工粘性 沙漏阻尼 静态阻尼第68页/共114页第六十八页，编辑于星期五：十一点 三十一分。2021-11-20阻尼选项 不连续 (不稳定) 二次粘性 二次粘性 + 线性粘性人工粘性的作用第69页/共114页第六十九页，编辑于星期五：十一点 三十一分。2021-11

28、-20阻尼选项 “沙漏”是一种网格畸变，它引起没有应变或体积没有改变； 四边形 (2D) /六边形 (3D) 使用单点积分产生沙漏； 一个阻尼选项用来抵制这种变形。沙漏阻尼第70页/共114页第七十页，编辑于星期五：十一点 三十一分。2021-11-20阻尼选项 静态阻尼系数 Rd ：静态阻尼 那么如果第71页/共114页第七十一页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器非结构化网格的单元/节点四面体多数四边形多数六面体第72页/共114页第七十二页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 输入选项LS-DYNA (.k) MSC.Nastr

29、an (.dat) 通过第三方软件输入： 网格 材料 边界条件 载荷 初始条件 目前不能产生非结构化网格 可以将结构化Parts (IJK)转变成非结构网格模型输入第73页/共114页第七十三页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 使用 Import 下拉菜单输入网格/模型： ICEM-CFD (ANSYS AI*Environment) TrueGrid LS-DYNA MSC Nastran 选择 Convert IJK Parts to Unstructured 菜单转换拉格朗日、壳和杆结构化 Parts 成非结构化 Parts。网格模型输入第74页/共11

30、4页第七十四页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 Retain LS-DYNA part definitions 选择该项，LS-DYNA parts / Nastran 属性将要被转变成 AUTODYN Parts 不选择，AUTODYN Parts 通过连通性创建 (如：六面体单元与四面体单元的连接) Merge duplicate materials 选择该项，使用相同模型和数据的材料将被合并 Check shell orientation 选择该项，检查所有壳单元外法向方向的一致性 (对流固耦合重要) 网格模型输入输入选项-LS-DYNA / MSC-N

31、astran第75页/共114页第七十五页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器输入后操作 材料 改变材料名字(选择) 添加材料数据(大多数材料) 可以增加刚性材料约束(如果需要的话) Parts、组件、边界条件和初始条件 重新命名(选择)网格模型输入第76页/共114页第七十六页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 整个 Parts 进行的操作 复制 填充删除 分区，平移非结构 Parts 操作测量单元/节点第77页/共114页第七十七页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 边界条件的施加在屏幕创建多变形

32、来选择节点，通过此来施加或删除边界条件非结构/结构 Parts 操作第78页/共114页第七十八页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器边界条件的施加： 通过屏幕选择节点来施加和删除边界条件。非结构/结构 Parts 操作第79页/共114页第七十九页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器 组的选择操作 用户定义节点、单元或面集合 交互式选择 通过多边形选取 通过面选取 施加/清除边界条件 节点边界条件 面边界条件 填充单元组 材料 速度 初始条件 非结构/结构 Parts 操作第80页/共114页第八十页，编辑于星期五：十一点 三十一分

33、。2021-11-20非结构求解器求 解 选 项六面体积分 有限体积公式：精确积分(缺省) 有限元公式：近似单点的高斯积分(快，精度不高)六面体沙漏控制 AD 标准 (缺省) 求解速度快，需要内存少； Flanagan-Belytschko (用于大刚体转动)。四面体压强积分 SCP (六面体网格为主) ANP (仅仅四面体网格)第81页/共114页第八十一页，编辑于星期五：十一点 三十一分。2021-11-20非结构求解器_练习拉格朗日(3D 非结构)头盔碰撞 输入 LS-DYNA网格： 头盔 头 空心圆柱 头盔用泡沫材料模型第82页/共114页第八十二页，编辑于星期五：十一点 三十一分。2

34、021-11-20基础培训二1、Lagrange求解器2、模型生成3、Lagrange之间的作用 第83页/共114页第八十三页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日连接第84页/共114页第八十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日连接 刚性连接可用于壳单元和梁单元 (具有旋转和平动自由度)UYA URA UXA UYB URB UXB 销连接UXA = UXBUYA = UYBURA URB刚性连接UXA = UXBUYA = UYBURA = URB用来指定 Part 与 Part 之间的连接方式 缺省方式为刚性连接

35、 连接处节点一致 第85页/共114页第八十五页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日连接例子梁单元和壳单元连接第86页/共114页第八十六页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 用于材料的滑移或分离： 内部破碎; 大剪切变形; 碰撞面接触; 靶板侵彻; 自身接触.接 触 作 用第87页/共114页第八十七页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 间隙的定义 每一个面段被一个接触搜索区域包围 接触搜索区域的半径就叫做间隙尺寸接触逻辑拉格朗日 壳单元 第88页/共114页第八十八页，编

36、辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 任何进入一个面断的接触搜索区域节点将受到一个反作用力，这个力的大小与该点进入接触搜索区域的深度成比例。3D 中的接触搜索区域FR R R R FF第89页/共114页第八十九页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触3D 中的接触搜索区域(分解图)第90页/共114页第九十页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 可见间隙 (2D & 3D) 作用面必须用间隙尺寸分开 可见/内部间隙内部间隙 (仅仅 2D 不推荐) 作用面可能处于接触中 第

37、91页/共114页第九十一页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 外部间隙 (2D & 3D) 最大间隙尺寸为 所有作用面中最小边的一半 内部间隙 (仅仅 2D 不推荐) 最大间隙尺寸为 所有作用面中最小单元尺寸的一半最大间隙尺寸第92页/共114页第九十二页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 如果只用点/面接触，有可能会接触失效 ； 这是就需要 edge/edge 接触 ； 在 Lag/Lag 作用菜单上选择； 增加计算时间。 Edge/Edge 接触第93页/共114页第九十三页，编辑于星期五：十一点

38、 三十一分。2021-11-20拉格朗日-拉格朗日接触Edge/Edge 接触梁 壳 拉格朗日 What are Edges?第94页/共114页第九十四页，编辑于星期五：十一点 三十一分。2021-11-20拉格朗日-拉格朗日接触 这个选项用来定义自身接触。 相应会增加计算时间。自身接触 当求解壳单元结构的时候，推荐设置 间隙尺寸为壳单元厚度的一半比较好。壳单元第95页/共114页第九十五页，编辑于星期五：十一点 三十一分。2021-11-20侵蚀 通过删除退化单元和严重变形单元，避免拉格朗日网格的重新分区； 侵蚀单元的内能被丢弃； 侵蚀单元的质量和动量可以被保留，也可以丢弃； 所有的材料模

39、型都能与失效模型结合求解； 侵蚀是一个非物理过程。 第96页/共114页第九十六页，编辑于星期五：十一点 三十一分。2021-11-20侵蚀带侵蚀作用Surface faces表面节点Part APart B第97页/共114页第九十七页，编辑于星期五：十一点 三十一分。2021-11-20侵蚀Surface faces自由节点ErodedErodedPart APart B表面节点带侵蚀作用第98页/共114页第九十八页，编辑于星期五：十一点 三十一分。2021-11-20侵蚀侵蚀应变选项Geometric Strain(from deformations)Effective PlasticStrain - EPS(from plastic material flow)InstantaneousStrain推荐推荐IncrementalStrain子程序 EXEROD 用来定义自己的侵蚀准则；对于小侵蚀情况，指定一个大侵蚀应变，用来仅仅来删除退化单元。第99页/共114页第九十九页，编辑于星期五：十一点 三十一分。202