★ FEATURED ARTICLE 简介LS-DYNA教程第二天课件聚焦显式动力学分析中的加载、边界条件与刚体定义面向已具备有限元基础、希望从隐式静力分析转向显式求解的工程师与高年级学生。内容从显式载荷必须为时间函数这一前提切入讲清TIMEFORCE与时间、载荷数组的配合方式说明D命令只作零位移约束说明CP与CE在大变形中仅对平动和转动有效等易错点并逐一展开阻尼类型选择、EDVEL初始速度、EDMP/RIGID刚体设定、EDDRELAX动力松弛与隐式—显式顺序求解各自的适用场合。一般加载流程按创建组元或Part、定义数组参数、EDLOAD施加载荷三步给出并提示压力作用于单元组元、刚体载荷作用于Part的差异。整包仅1个PPT文件约3.07MB章节按主题A至H划分并配有例题便于课堂讲授或自学对照。目前已有299人学习。1. 从隐式切到显式F 和 SF 为什么在 LS-DYNA 里直接失效把一份隐式静力的 deck 直接丢给 LS-DYNA 显式求解器最常见的翻车点不是网格也不是材料而是加载这一步F 和 SF 加上的力没有跟时间绑在一起求解器读到的是一个常量瞬时加载在显式里等于给结构一记硬冲击响应曲线能难看到让人怀疑人生。LS-DYNA 教程第二天的内容正好压在这个转换点上——加载、边界条件、刚体、初始速度、阻尼、动力松弛这几件事在显式动力学里得换一套逻辑没有时间载荷步所有载荷先写成“时间数组 载荷数组”的一对曲线再由 EDLOAD 挂到节点组元、单元组元或 Part 上。这份材料适合已经能独立建模、但第一次把模型交给显式求解器的人也适合做了一段时间却总在“加载到底对不对”上反复返工的老手。2. 载荷时间曲线与 EDLOAD组元、数组、加载三件套显式求解器读载荷不是读一个数而是读一条对时间的函数。所以整条加载链路上真正决定成败的是三件事载荷挂在谁身上、时间与载荷怎么配对、EDLOAD 的参数位填得对不对。这三件任何一件出偏差结果都不会报一个漂亮的错误而是安静地给出错误响应排查成本极高。我一般按“先建载体、再定曲线、最后加载、加载完必验”的顺序走一步不跳。2.1 组元和 Part 的分工节点组元吃 FX单元组元吃 PRESS刚体吃 RBFXANSYS/LS-DYNA 里能当载荷载体的对象只有三类节点组元、单元组元、Part。这个划分是求解器按自由度归属做的硬约束不是建议。力和位移类载荷作用在节点自由度上压力必须知道具体是哪个单元面刚体载荷必须整块处理——刚体内部的节点自由度在求解阶段已经被凝聚掉往单个节点上施加力没有意义。载荷载体建法可用标签典型用途节点组元NSEL 选节点后 CMFX/FY/FZ、UX/UY/UZ、VX/VY/VZ、AX/AY/AZ、ACLX/ACLY/ACLZ、MX/MY/MZ、ROTX/ROTY/ROTZ、OMGX/OMGY/OMGZ、TEMP局部冲击力、把面力折算成节点力、初速度与初温度单元组元ESEL 选单元后 CMPRESS面压力、爆轰压力、分布气动载荷Part由 Part 号或刚体标识RBFX/RBFY/RBFZ、RBUX/RBUY/RBUZ、RBVX/RBVY/RBVZ、RBMX/RBMY/RBMZ、RBRX/RBRY/RBRZ、RBOX/RBOY/RBOZ刚体整体受力、刚体整体初速度或角速度建组元的动作本身很简单坑在命名和选择集管理上。模型一大几十个组元名字全是COMP1、COMP2加载时手动选错一个跑完才发现力加到了不该加的面只能重来。我习惯用“部位 类型”的前缀命名比如IMPACT_NODES、PRESS_ELEM、PART_HEAD看名字就知道该挂什么载荷。/PREP7 ! 冲击面节点组元先按整体坐标粗选再在结果集内二次筛选 NSEL,S,LOC,Y,0 NSEL,R,LOC,X,-30,30 CM,IMPACT_NODES,NODE NSEL,ALL ! 压力必须挂在单元组元上挂到节点组元上求解器读不到面 ESEL,S,MAT,,2 ESEL,R,CENT,X,0,50 CM,PRESS_ELEM,ELEM ESEL,ALL这段里NSEL,S是从全体里选NSEL,R是在当前选择集内再筛两者顺序反了结果完全不同这是 APDL 的老规矩但在加载场景下特别容易出错因为选错之后 CM 依然能成功建立只是内容不对。CM操作的对象永远是当前选择集所以NSEL,ALL/ESEL,ALL一定要跟在后面复位否则后续建模操作会莫名其妙只作用于一部分实体。压力挂在单元组元时面号按单元类型的面编号约定填SOLID164 是 1 到 6SHELL163 按壳的面定义这个号填错的表现是压力方向反了或者压到了背面位移结果会给出完全相反的符号。2.2 用 *DIM 定义时间与载荷数组终止时间要留出余量显式动力学里没有任何“把某个力一直加着”的写法力要么是常数曲线两个点时间 0 和求解结束都给同一个值要么是一条真正的时程。数组必须成对存在一个时间数组、一个载荷数组长度一致、下标一一对应。三个细节决定这条曲线能不能用。第一时间必须严格单调递增不允许出现两个相同的时间点。重复时间点会让曲线在该处的斜率变成无穷大求解器在插值时行为不可预测表现出来就是某个时刻突然出现一个尖峰。第二数组下标从 1 开始这是 APDL 的语法不是 0。写TM(0)不会报错会静默做出奇怪的事情。第三也是最容易忘的一条载荷数组的终止时间要延长到求解结束时间之后。一旦时间推进到曲线最后一个时间点载荷会被迅速置零等于给结构来了一个突然卸载尾部响应会出现一段明显不真实的振荡。后续如果还要重启动这个尾巴会直接被带进下一段分析。! 求解总时长 0.02 s载荷曲线故意延到 0.03 s *DIM,T_LOAD,ARRAY,4 *DIM,F_Z,ARRAY,4 T_LOAD(1) 0.0, 1.0E-3, 5.0E-3, 3.0E-2 F_Z(1) 0.0, 1200.0, 1200.0, 1200.0这里前两个点构成一段 1 ms 的斜坡用来避免阶跃加载激起的初始振荡平台段拉到 5 ms 之后就不再变化但数组一直延续到 30 ms比 20 ms 的求解时间多出 10 ms 的余量。1200这个量级不是随便给的它应该来自静力预估——如果加载后的等效静位移远超预期说明坡升太快或者量级本身有问题。做完曲线先别急着跑看一眼时间跨度是不是覆盖了整个求解窗口这是成本最低的一次自检。2.3 EDLOAD 的参数位、EDCURVE 复用与开始/终结时间EDLOAD 的参数位是这套流程里最容易填错的一段因为它一个位置在不同载荷类型下含义不同。关键位有两个一个是载荷标签FX、PRESS、RBFX 这一类决定了这条曲线作用到哪类自由度另一个是坐标/面号位对压力载荷它是单元面号对其他类型它可以指向一个由 EDLCS 定义的局部坐标系。这两个位置填错命令照样执行结果就是力加在奇怪的坐标系里或者压力压在错误面上。/SOLU ! 直接给数组时间数组 载荷数组成对传入 EDLOAD,ADD,FZ,0,IMPACT_NODES,T_LOAD,F_Z ! 同一条曲线要挂到多个组元上时先建曲线号再引用避免重复定义数组 EDCURVE,ADD,7,T_LOAD,F_Z EDLOAD,ADD,FZ,0,IMPACT_NODES, , , ,7 EDLOAD,ADD,FZ,0,NODETOP, , , ,7 ! 压力挂单元组元第三位填单元面号 EDLOAD,ADD,PRESS,2,PRESS_ELEM,T_LOAD,P_PRESSEDCURVE那两行的价值在于复用同一个时程要加在三个不同的组元上时不需要定义三份数组一条曲线编号到处引用就行。曲线引用支持载荷缩放系数但不支持时间缩放需要拉伸时间轴只能回去改数组本身。这也是为什么我习惯把时间数组单独命名缩放需求一来改一个地方就够。参数顺序每个版本可能略有差异最稳的校准办法是在 GUI 里手动走一遍加载菜单然后去工作目录翻.log文件里面会记下完整的等价命令行照着抄比自己猜参数位可靠得多。对于多种载荷叠加的工况与其在数组里给不加载的时间段赋零值不如直接用 EDLOAD 的开始时间和终结时间参数控制激活窗口。赋零值在曲线上会形成一个台阶斜率跳变同样会激起不必要的响应给窗口的写法在数值上干净得多。最后一个参数是分析阶段默认Transient only做动力松弛时选Dynamic relax两段都要就选Trans and Dynam——这个选项的坑放在第 5 章细说。2.4 加载后的验证EDPL 看曲线EDFPLOT 看符号加载命令执行完不等于加载正确。我见过太多次命令没有任何报错、跑出来的结果却完全不对的情况原因都是组元选错或者曲线引用错了编号。所以加载后必须走两步验证而且顺序不能反。第一步用 EDPL 把每条曲线画出来看形状、看幅值、看终止时刻。这一步能抓住绝大多数数组层面的错误长度不一致、时间没递增、终止时间压线、幅值少了三个零。曲线画出来那一刻很多问题肉眼就能看出来。第二步用 EDFPLOT 打开载荷符号显示。符号只会标注载荷施加在哪个组元或哪个 Part 上并不会告诉你方向和数值而且符号永远沿坐标系正向绘制不代表实际受力方向。所以符号只用来确认“有没有加、加在谁身上”方向和大小必须回去看曲线。EDPL,7 ! 画出 7 号载荷曲线核对形状与终止时刻 EDFPLOT,ON ! 打开载荷符号显示 EPLOT ! 重绘观察符号落在哪些组元上 EDLOAD,LIST ! 列表核对载荷类型、目标组元、时间与数值EDLOAD,LIST的输出值得逐行看一遍它把载荷类型、作用的组元名、时间数组和载荷数组都列出来是唯一能把“我以为加了什么”和“实际加了什么”直接对上的地方。跑之前养成看这个列表的习惯能省下大量跑完再回查的时间。符号显示在 GUI 里的路径是 PlotCtrls 下的 Symbols 菜单进去勾上 LS-DYNA 载荷符号即可记得每次重新 EPLOT 才会刷新。3. 边界条件不是载荷D 命令、CP/CE 与 EDCNSTR 的适用边界初学者最容易混的一对概念就是“约束”和“载荷”。在隐式里两者界限还算清楚到了显式里因为时间相关性这一条界限被放大了凡是和时间无关的东西都不能用来加载。约束的作用是把某些自由度钉死或者绑在一起它本身不输入能量所以它和时间无关也没有关系。这条分界线一旦想清楚很多“为什么这条命令没反应”的问题就自动消失了。3.1 D 命令只做永久零位移约束不能用来加载D 命令在 LS-DYNA 里依然可用但只限于一种语义把指定节点的某个自由度约束为永久的零位移。它不跟时间走所以不能拿来做强迫位移加载也不能用来做随时间变化的位移边界。有人试图用 D 命令施加一个位移来模拟冲头下压结果模型完全不动就是因为这条命令根本没有载荷通道。约束的施加对象应该尽量收敛到最小的必要集合。显式求解的稳定时间步由最小单元尺寸决定而约束过多会让局部刚度突变在约束附近产生应力集中这些位置的应力值往往不可信。我的习惯是在远离关心区域的地方布约束关心区域内的约束一律用刚体或接触去替代。想做的事该用的东西不该用的东西固定底面D 命令施加零位移给底面加一个恒定为零的力曲线让冲头匀速下压给刚体 Part 一个初速度或 RB 位移曲线D 命令施加位移让两个面始终贴合接触或 EDCNSTR 绑定CP 耦合节点施加分布压力EDLOAD 压力挂单元组元SF 命令表格里最后一行值得单独说SF 这类隐式加载命令在显式分析里不可用因为它的语义是“在这个面上加一个不变的载荷”而显式需要的是时程。这条差异不是版本问题是求解器架构决定的。3.2 CP 和 CE 在显式里只对 U 和 ROT 有效耦合CP和约束方程CE在显式分析里是被保留的但适用范围被严格限制在位移的平动和转动自由度上。也就是说用 CP 把一组节点的 UX 绑在一起是可以的想通过它传递力或者其他物理量就不行。这个限制很多文档一句话带过实际用起来出问题的概率却不低。真正麻烦的是大变形场景。CP 和 CE 的约束关系是在建模阶段按初始几何定义的节点之间的相对位置一旦在大变形中发生显著变化约束方程描述的几何关系就和实际变形脱节了。表现出来是约束区域附近出现异常的应力集中或者变形模式明显不符合物理直觉但求解器不会给任何警告。这类问题排查起来极耗时间因为模型能跑完、能量也守恒、就是结果不对。所以我的判断标准很简单变形量级在弹性小变形范围内CP 和 CE 可以放心用一旦涉及大变形、大转动、材料失效后的碎片运动就不要在大变形区堆 CP 和 CE换别的路子。3.3 大变形下用 EDCNSTR 收口EDCNSTR 是显式分析里专门给约束准备的一组工具典型用途是把一组节点绑成一个整体或者做节点刚性连接、绑定断开这类行为。它和 CP 的区别在于约束是在求解过程中按当前构型处理的不会因为几何变化而失真。! 把螺栓孔周围一圈节点绑成一个刚性整体替代在该处堆叠 CP CM,NODEBOLT,NODE EDCNSTR,ADD,NRB,1,NODEBOLT ! 跨部件的连接优先考虑绑定而不是把两边节点一个个 CP 起来 EDCNSTR,ADD,TD,1,TIE_SETNRB这类节点刚性体的用法在螺栓连接、局部加固区特别顺手它把一组节点凝聚成一个整体自由度数量直接降下来对计算时间的帮助很明显。TD之类的绑定类型适合两个部件之间的面连接比手动 CP 稳定得多。参数位的具体含义和可选类型以本机命令手册为准各版本之间有些差异建议先用一个小模型试一次确认约束方向和范围符合预期再上大模型。约束加完之后一定要跑一个极短的验证算例看约束区域的接触力或者支反力是否在合理量级这一步能提前发现方向反了或者范围选大了的问题。4. 初始速度、刚体与阻尼三处最容易被忽略的省时设置模型建好、载荷曲线也写完了剩下三件事看起来都是配角实际对计算时间和结果稳定性的影响比很多主流程设置都大。初始速度能避免渐增加载带来的初始振荡刚体能砍掉大量自由度阻尼能压住不该出现的数值振荡。这三样用对了同一个模型的计算时间能差出几倍。4.1 EDVEL 的两条路径VGEN 与 VELO 怎么选EDVEL 只干一件事给节点或节点组元施加初始速度。它有两条路径一条是 VGEN一条是 VELO选择依据主要看你要的是平动、绕轴转动还是直接往节点上写转动自由度。VGEN 定义的是相对总体坐标系的平动速度外加绕任意轴的转动。它需要给出旋转轴的基点坐标和方向向量所以适合模拟一个整体在空间中绕某条轴线旋转的工况比如一个绕固定轴翻转的构件。VGEN 也能对已经成组元的对象施加速度这在批量处理上很方便。VELO 则是把初始平动或转动速度直接写到节点自由度上。正因为是直接写自由度它只对具有对应转动自由度的单元有效BEAM161 和 SHELL163 是典型代表。对实体单元用 VELO 施加转动不会有任何效果因为 SOLID164 的节点根本没有转动自由度可写。命令施加对象转动能力重启动限制EDVEL,VGEN节点或节点组元可绕任意轴需给轴点与轴方向不能用于简单重启动或完全重启动EDVEL,VELO节点或节点组元直接写节点转动自由度仅限 BEAM161/SHELL163不能用于完全重启动EDPVELPart 或 Part 集合同 VGEN 语义按 Part 批量施加受 Part 定义约束! VGEN整体沿 X 平动同时绕 Z 轴转动 ! 参数位依次为组元名、平动速度分量、旋转轴基点、轴方向、角速度 EDVEL,VGEN,NODEWHEEL,30.0,0.0,0.0, 0.0,0.0,0.0, 0.0,0.0,1.0, 120.0 ! VELO只对带转动自由度的单元有意义 EDVEL,VELO,NODEEDGE,0.0,0.0,-5.0, 0.0,0.0,0.0上面第一段里30.0,0.0,0.0是平动速度分量后面三组坐标一个定义轴基点、一个定义轴方向最后一个数是角速度。这套参数位容易在中间几个零上出错因为连续的零看不出位置建议每个数后面跟注释或者干脆用 GUI 生成一次再抄日志。第二段的VELO版本只给了平动分量转动分量留零因为一旦给非零转动而单元类型不对命令不会报错速度会悄悄丢失。重启动这条限制要特别记住如果打算做完全重启动VGEN 和 VELO 都不能用。常见现象是重启动之后模型初始速度凭空消失整个结构静止在那里跑完才发现第一步就没动起来。需要重启动的场景初速度要么在第一次计算里就一次性给对要么用 EDPVEL 按 Part 处理。4.2 EDPVEL 与 EDMP,RIGID刚体载荷那一组 RB 标签刚体的价值在于减少模型自由度。一个刚体无论网格多密内部自由度在求解时都被凝聚掉了计算量只跟它的质量和惯性张量有关。对于变形不是重点的部件——夹具、冲头、底座、配重块——设成刚体是最直接的一笔时间收益。用 EDMP 的 RIGID 选项把材料或者 Part 标成刚体之后这个 Part 上的载荷就不能再用 FX 这类节点标签了必须换成 RB 系列RBFX 是刚体整体受力RBUX 是刚体整体位移RBVX 是整体速度RBMX 是整体力矩RBRX 是整体转角RBOX 是整体角速度。这套标签操作的是刚体作为整体的运动不是某个节点。刚体服从 F ma所以想给它一个平动加速度直接加等效的力就行不需要专门找加速度标签。/PREP7 EDMP,RIGID,1,0,0 ! 材料 1 的整体运动按刚体处理 /SOLU EDPVEL,ADD,PART_HEAD,0.0,0.0,-20.0 ! 按 Part 给整体初速度 EDLOAD,ADD,RBFY,0,PART_HEAD,T_LOAD,F_Y ! 刚体载荷必须用 RB 系列标签EDMP,RIGID后面跟的材料号必须和实际材料定义对上改成刚体之后该材料的本构参数基本不再起作用但密度必须给对否则质量和惯性矩全错动力学响应跟着错。EDPVEL是给 Part 施加初速度的入口它内部其实就是给该 Part 的所有节点建一个集合并施加所以用法和 EDVEL 类似。RB 标签那条加载命令如果错写成 FX命令可能执行成功但完全不生效因为刚体的节点自由度已经被凝聚没有一个节点能接收这个力。这种情况没有任何报错只有结果里的刚体纹丝不动是最难查的一类问题。4.3 阻尼的量级与施加位置阻尼的作用是减小结构对载荷的动力响应本质上是耗散能量。显式动力分析里常见的阻尼作用位置有两层全局阻尼和材料/Part 级阻尼。全局阻尼通常用 EDDAMP 设定质量阻尼系数和刚度阻尼系数质量阻尼主要压制低频振荡刚度阻尼主要压制高频振荡材料级阻尼则跟具体本构绑定作用范围更局部。量级是最难拿捏的部分。给太小等于没给振荡照样在给太大能量被过度耗散结果看起来“特别稳”但那是假的稳——结构的真实动力响应被吃掉了。我的经验是从一个很小的系数起步跑一个短算例看能量曲线里阻尼耗散能占总能量的比例超过百分之几就要往回退。判据不是“看起来平滑”而是阻尼耗散能在总能量里的占比合理动能和内能的转换关系还符合物理预期。! 全局阻尼质量阻尼与刚度阻尼系数分开给先从小量级试 EDDAMP,0.05,0.002这两个系数的物理意义不同不能互相替代。第一个偏质量项压低频第二个偏刚度项压高频。什么时候需要动阻尼两种情况最常见一是初始加载阶段有明显的高频振荡需要在起步阶段压一下二是做动力松弛靠阻尼把动能逐步耗到零让系统落到静力平衡位置。后者是阻尼最重要的用法第 6 章展开。5. 排查清单六个高频翻车点与固定验证顺序加载这块的问题有个共同特点命令几乎不会报错模型照样跑完错误全部藏在结果里。这就需要一套固定的排查顺序把问题在小成本阶段拦下来。下面六条是我自己踩过、也帮别人查过的记录按现象、原因、解决三段写每条都对应一个具体的命令行为。5.1 现象、原因、解决六条实战记录第一条尾部剧烈振荡。现象是位移时程曲线在求解末尾突然出现一段大振幅振荡能量曲线里的动能末尾翘起来。原因是载荷曲线的最后一个时间点正好压在求解结束时刻载荷被瞬间置零等于给结构一个突然卸载。解决是把载荷数组的终止时间延长到求解结束之后一般多留 20% 到 50% 的时间窗口如果后面还要重启动延长得更多一些。第二条同一位置反复赋零值。现象是某个不该有响应的时段出现了小幅度振荡。原因是在载荷数组里用一串零值表示“这段时间不加载”零值段和加载段之间形成了台阶斜率跳变在数值上等效于一次冲击。解决是不用赋零值改用 EDLOAD 的开始时间和终结时间参数控制激活窗口曲线本身只保留真正有效的那一段。第三条压力加上去没反应。现象是压力载荷命令执行成功模型不动。原因是压力挂到了节点组元上而 PRESS 标签要求载体是单元组元压力需要知道具体作用在哪个单元面上。解决是回去用 ESEL 建单元组元并在 EDLOAD 的坐标位填对单元面号SOLID164 的面号是 1 到 6填反了压力方向就反了。第四条热载荷互相覆盖。现象是温度场施加完之后之前设的温度全被替换掉了。原因是 LDREAD 会覆盖 TUNIF/BFUNIF 施加的温度而且瞬态热载荷和静态热载荷在同一个模型里不能同时存在。解决是先确定这次分析到底是静态预载还是瞬态热载二选一隐式-显式顺序求解走静态路线冲击热载走瞬态路线两条不要混。第五条重启动后初速度消失。现象是重启动算出来的模型从头到尾不动。原因是 VGEN 和 VELO 都不支持完全重启动初速度在重启动时不会带过去。解决是重启动场景改用 EDPVEL 按 Part 施加或者在第一次计算里就把初速度固定下来后续不再依赖它。第六条动力松弛之后出现莫名振动。现象是松弛阶段好好的进入瞬态阶段后结构开始剧烈振动。原因是分析阶段选项用了“Trans and Dynam”并且只有一条曲线松弛结束后时间被重置为零预加载荷被移走然后重新渐变加载这个重新加载的过程本身就是一次激励。解决是用两条曲线一条专门给松弛阶段做从零渐变另一条从松弛结束点接着给瞬态载荷分析阶段分开指定。5.2 一套固定的验证顺序上面这些问题的共同点是都不报错所以只能靠主动验证。我把验证固化成四步每一步成本都不高但能拦掉大部分低级错误。第一步EDLOAD,LIST逐行核对。载荷类型对不对、目标组元是不是想要的那个、时间数组和载荷数组的名字有没有引错。这一步是唯一能把意图和实际配置直接对上的地方。第二步EDPL画曲线。看形状、看幅值量级、看终止时刻有没有压线、看有没有意外的时间点重复。第三步EDFPLOT,ON打开符号显示再重绘。确认符号落在预期的组元或 Part 上符号数量对不对。符号只标位置不标方向方向要回去看曲线和坐标系定义。第四步跑一个极短的算例比如只跑几十个时间步然后去看能量汇总和支反力。动能是不是从零平滑上升、内能有没有异常跳变、约束位置的支反力是不是在预估量级。这一步花的时间不到一分钟能提前发现量级错位、方向反了、约束过刚这几类问题。做完这四步再上完整算例比跑完几小时再回头查要划算得多。6. 动力松弛与热载荷用两条曲线把预载做干净前面所有加载都是瞬态语义但有一类工况需要先让系统落到一个静力平衡状态再开始动力学计算——预应力结构、装配预紧、重力场下的初始沉降都属于这一类。EDDRELAX 提供三档模式OFF是纯瞬态也是绝大多数分析该用的默认值DYNA让显式求解器通过逐步增加阻尼把动能耗到零自己找到静力平衡ANSYS则要求先做一次隐式静力分析把变形和应力读进来做应力初始化也就是常说的隐式-显式顺序求解。DYNA这一档的关键全在载荷曲线上。松弛阶段需要一条从零开始渐变到稳定值的曲线让阻尼有足够时间把动能吃掉瞬态阶段需要另一条从松弛终点接续的曲线。如果只给一条曲线又选了“Trans and Dynam”松弛结束后时间会被重置为零之前预加的载荷被移走然后重新渐变加载一次这个二次加载本身就是一次激励模型会在这个环节莫名其妙地振动起来。我一般会老老实实写两条曲线一条标Dynamic relax一条标Transient only两段的起始幅值对齐接缝处就不会有跳变。热载荷的处理逻辑和结构载荷不太一样因为它是两套互斥的机制。静态热载荷走 TUNIF/BFUNIF 施加均匀温度、走 LDREAD 读取热分析结果做非均匀温度场后者假设温度不随时间变化并且会覆盖前者瞬态热载荷走 EDLOAD 的 TEMP 标签挂在节点组元上。这两套不能混用同一个模型里也不能同时存在瞬态和静态热载荷。结构载荷倒是可以和热载荷共存静态瞬态都行。材料这一侧有个容易忽略的约束热载要求材料带温度相关性支持的是双线性各向同性模型最多定义六组温度数据实际温度不能超出这个范围。模拟热弹性材料时把屈服应力设得足够大让它始终不屈服。单元支持方面PLANE162、SHELL163、SOLID164 都可以但 PLANE162 不支持动力松弛文件只能配瞬态热载或冲击热载做瞬态热-显式分析时把 ALPX 设为 0只保留与温度相关的力学性能稳定性会好很多。另外别忘了实际施加的温度是输入值减去参考温度 TREF默认 TREF 为 0参考温度设错整条热应变曲线就整体平移。从那以后我每次配 EDLOAD不管模型多简单都强制走一遍“列表、曲线、符号、短算例”这四步宁可多花几分钟也不愿意第二天早上看到一个跑完但完全不对的结果。希望帮到你。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
