TMS BOOK · ACADEMY 课题

结构 FEA 建模方法基础:单元、约束、本构与结果可信度

剔除既有课已讲内容后,七项真零覆盖(405 门 outline + 134 篇讲义均查证为 0)构成本课主体,足够 5 节 L2,与 J2-01 体量对等: (1) 结构 FEA 的几何前处理——中面抽取(等厚区自动抽取、变厚过渡与中面断裂的手工缝合、抽中面后厚度属性怎么带)、小特征清理准则(圆角/小孔/倒角/卡扣/logo 删到什么程度、删除对应力集中与刚度结果的影响边界)、几何缺陷(缝隙、自交、拓扑错误、单位与精度)在结构侧的后果;与 (2) 单元选型合并入第 1 讲,不另加讲次。 (2) 单元类型与积分阶次——壳/实体/接触/垫片单元的适用域与判据、一阶 vs 二阶单元在弯曲问题上的差别、完全积分 vs 缩减积分与沙漏控制。

J5-06 · 仿真分析 CAE / 结构/热应力/疲劳 FEA
L2 进阶
待认领listed 众筹中recruiting 已开课open · Course 售卖 一课题可并存多讲师版本

课程大纲

学完能做什么
  • 能按响应量与壁厚特征判断该抽中面用壳单元还是直接用实体,并给出小特征的删/留清单与厚度属性(含变厚区等效厚度)的回带方式
  • 能按弯曲/接触/自锁风险选单元类型与积分阶次,并用沙漏能占比与位移收敛结果反查选型是否成立
  • 能做模型的六自由度体检,区分欠约束与过约束,按适用条件正确使用对称约束、RBE2/RBE3 与惯性释放,排除凭空造出的虚假热应力
  • 能给出螺栓预紧与接触对的仿真侧设置(预紧两步法、摩擦系数上下界包络、接触收敛调法),并说清预紧力取值来自设计侧而非仿真侧
  • 能判定何时必须上几何/材料/接触非线性,并按工况选出弹塑性/蠕变/超弹的本构类别与温度相关属性的建入方式
  • 能区分应力集中与应力奇异、做结构侧的分层收敛检查,并沿壁厚做应力线性化给出合法的评定应力(Pm、Pm+Pb),选对 von Mises 或最大主应力
内容大纲
  1. 几何前处理:壳/实体判定、中面抽取与小特征清理
    • 中面抽取三件事:等厚区自动抽取可用;变厚过渡区要么按分区分别给厚度、要么就近取等效厚度(等效原则按响应量定——刚度类问题按抗弯等效、一次强度类问题按最小壁厚保守取);中面在交叉、变厚、圆角处会断裂,必须手工缝合成连续曲面,否则传力路径断开
    • 抽完中面厚度属性必须显式回带到壳单元(厚度、偏置方向、上下表面对应关系);漏带或偏置方向反了,弯曲应力符号与位置会整体错,是一类查不出来的静默错误
    • 壳还是实体的判据本体在本课:壁厚方向应力可忽略、特征跨度与壁厚之比在 10 倍以上量级、关心的是膜+弯响应时用壳;需要壁厚方向温度/应力梯度、三向应力状态、局部接触压力时用实体。J5-01「薄壁件用壳单元还是实体单元的取舍」那一条要点的判据在这里
    • 小特征清理分级判据:判据是「该特征是否在关注响应量的传力路径上、是否是评定点」——承载路径上的圆角一律不删,非承载的 logo、工艺凸台、装配标记先删;小孔、卡扣按是否传载分别处置
    • 删简的影响边界要分响应量说:删掉远小于结构特征尺寸的小特征,对整体刚度、位移、模态频率的影响通常在百分之几量级;但对局部应力是数量级差别。刚度类问题的删简自由度远大于应力类问题
    • 删圆角的代价必须当场讲透:把一个有限的应力集中变成了无界的应力奇异——那个位置从此不能作为评定点,与第 4 讲互为呼应
    • 结构侧的删简判据与 CFD 侧不重讲:响应量不同(结构看应力与刚度,CFD 看流阻与换热),判据来源随之不同。脏数模验收、破面缝合、包面、计算域抽取一律前指 J2-06,本课只讲结构侧
    • 「小特征」在 J3-03 指的是拓扑优化里的最小特征尺寸制造约束,与本课的小特征清理同形异义,不要混用
    • 几何缺陷在结构侧的后果与接收口径:缝隙(接触对与连接建不上、网格在此断开、传力路径丢失)、自交与拓扑错误(网格划分失败或生成负体积/翻转单元)、单位与精度(单位制不一致是量纲级错误,须做 mm-t-s-MPa 一类一致性自检;导出精度不足会在小圆角处生成锯齿面)
  2. 单元类型与积分阶次:剪切自锁、沙漏模式与选型反查
    • 单元类型适用域:壳(薄壁承压件、钣金支架)/实体(厚壁、压铸壳体、需壁厚方向梯度)/接触(会分离-闭合的界面)/垫片(密封垫的压缩-回弹,其本构标定与接触压力分布判据见 J5-08)/梁(长细支撑与螺栓的等效替代)
    • 一阶 vs 二阶单元在弯曲问题上的差别:一阶四面体因剪切自锁在弯曲主导问题上系统性偏刚——位移偏小、频率偏高,方向偏危险,弯曲主导问题禁用;改用二阶四面体,或用一阶六面体配缩减积分
    • 完全积分 vs 缩减积分:缩减积分省算力并缓解自锁,代价是可能出零能沙漏模式;判据是沙漏能与内能之比(见公式条),超门槛说明结果被非物理变形污染,须加密或提高沙漏控制刚度后重算
  3. 约束与刚性连接:六自由度体检与虚假热应力的排除
    • 先做六自由度体检:任何模型求解前先确认刚体位移是否被约束住、且只被约束住。欠约束表现为求解器报刚度矩阵奇异或位移量级异常大;过约束表现为约束反力与外载不平衡、应力集中出现在约束点而不是结构薄弱处
    • 固定/滑动/对称约束的合法条件:对称约束要求几何、材料、载荷、约束四者全对称——温度场不对称就不能用对称面(这一条与 J5-02 的热载荷工况直接相关);滑动约束要指明放开的是哪个方向,而不是「大概能滑」
    • 刚性连接与分布耦合分账:RBE2 把从节点刚性绑死,会在连接区局部增刚并憋死热膨胀;RBE3 只做载荷分配、不加刚度,不能拿来当承力约束。用 RBE2 顶替 RBE3 是虚假热应力最常见的来源之一
    • 过约束造出的虚假热应力是本课的责任范围:模型里六自由度全锁、对称面用错、刚性连接跨过了需要自由膨胀的方向,会凭空造出并不存在的应力。真实结构该不该留腰形孔/滑套、约束反力的物理量级(F=EA·α·ΔT 上界与 Δα 柔度串联的真值)属设计侧,前指 I4-01,两侧不重讲
    • 惯性释放的适用条件:无外部约束的自由件、载荷自平衡(如仅受自重、仅受内压的孤立件)、准静态问题才可用;有真实支撑反力的工况不能用惯性释放糊过去。J5-02 那条「不要随手锁一个点」的纪律,其方法层(何时用、怎么加、怎么校验反力残差)在本课
  4. 螺栓预紧与接触界面:两步法、摩擦包络与接触收敛
    • 螺栓预紧的仿真侧施加:在螺栓杆中段设预紧截面/预紧单元,分两步求解——第一步施加预紧力并锁定,第二步加工作载荷(含温度)观察预紧力变化。预紧力与扭矩取值(T=K·F_pre·d、VDI 2230 手算)来自设计侧 I4-01,本课不重算手册值
    • 接触对设置口径:主从面选择(网格较粗、刚度较大者作主面)、初始间隙/穿透与过盈的处理方式(几何调整 vs 逐步施加)、探测方式与容差;接触状态输出(有效接触面积、穿透量)必须作为结果的一部分被检查,而不是只看应力云图
    • 摩擦系数在仿真里是敏感参数,不是查表定值:干态金属界面通常取 0.1~0.2 量级作起点,须按表面状态与试验确认;正确做法是取上下界各算一次做包络,看结论会不会翻转
    • 接触收敛调法:先建立接触再加载的分步策略、载荷步细分与自动步长、接触刚度与穿透容差的取舍。不收敛时先查是不是初始就没接触上(悬空)或过盈过大,再去调求解参数——调参在几何/初始状态错的模型上不会收敛
  5. 非线性分界与材料本构:什么时候必须上,上哪一种
    • 三类非线性分账清楚:几何非线性(大位移/大转动/应力刚化,判据量级是位移与结构特征尺寸同量级、薄板面外位移达到板厚量级、或结构受载后刚度明显变化)、材料非线性(应力越过比例极限/屈服)、接触非线性(接触状态在加载过程中改变)
    • 何时必须上的判据:线性解算出的最大应力已越过屈服却仍按线弹性报数,结果不可用——不是「偏保守」而是位置与量值都可能错。J5-01 的误区只提到「遗漏塑性/爆破风险」,怎么判、怎么上在本课
    • 线性解是入口不是答案:先跑线弹性,看最大应力/屈服之比、最大位移/特征尺寸之比、接触面是否发生状态变化,按这三个比值决定要不要升级到非线性——这样升级是有依据的,不是凭感觉「保险起见上非线性」
    • 温度相关材料属性怎么建:E、屈服、线胀系数随温度变化,热管理件工况常跨 -40 到 150℃ 量级(制冷剂高压侧、钎焊工艺态更宽);高温段沿用室温 E 与屈服会把应力和裕度都算乐观。属性表要给到工况上下界之外,禁止靠外推撑出边界点
    • 材料取值口径不在本课:钎后态(3xxx 按钎后 O 态、6xxx 按钎后实际时效态)、压铸件本体 vs 单铸试棒、钎缝强度折减一律前指 J5-01 与 H1-01/H1-04。本课只讲怎么把已定的属性和本构类型正确建进模型
    • 弹塑性本构选用:理想弹塑性(只判是否进入塑性)、双线性硬化(估塑性区)、多段真应力-真应变曲线(算残余变形与塑性垮塌)。手册给的工程应力应变必须转成真应力-真应变才能进大变形塑性模型,直接代入会低估硬化
    • 蠕变/黏弹与超弹的选用条件:金属在高温长时保载下的预紧力衰减要蠕变本构;橡胶、密封件的大变形要超弹本构;纤维增强塑料要各向异性弹性常数。本课只给「什么条件下必须选这一类」,试验标定方法(多组曲线联合拟合、时温等效、取向场映射)与密封接触判据全部前指 J5-08
    • 钎焊冷却残余应力那一段的高温工艺态本构属工艺语境,见 H2-07,本课不重讲;本课只负责说明:残余应力作为初始状态被忽略时,后续热应力/寿命结果偏乐观
    • 本构选错的代价方向要说清:线弹性在超屈服区算出的应力偏高,对一次静强度是偏保守;但同一个模型的位移、变形形状与失效位置会算错,对密封面贴合、间隙包络、卡边剩余压缩量这类判断是偏危险的——不能拿「偏保守」一句话覆盖全部响应量
  6. 应力集中 vs 应力奇异:结构侧网格收敛的特殊性
    • 应力集中是有限值:几何过渡(圆角、开孔边、变截面)处的应力放大,K_t 定义见 J5-01;网格加密后该点应力单调趋近一个有限平台值,可以做收敛判定并作为评定点
    • 应力奇异是无限值:零圆角尖角、集中力/点约束、刚性连接的边界节点、焊趾与材料界面尖点,理论解的应力无界,加密只会单调发散、永不收敛——这类点上的「最大应力」不是物理量
    • 判别程序要落到操作:只加密关注区,配三套系统加密的网格(加密比约 1.5~2),拟合峰值应力与单元尺寸的关系(见公式条);出现平台判为集中,无平台且单调发散判为奇异。只用两套网格无法区分「慢收敛」与「弱奇异」
    • 结构侧收敛与 CFD 的「网格无关性」不是一回事:CFD 侧的 y+、边界层与 GCI 口径见 J2-01,本课不重讲。结构侧是分层收敛——位移与模态频率收敛最快,应力次之,奇异点根本不收敛;报数时必须声明「收敛的是哪个量」,不能拿位移收敛证明应力可信
    • 奇异处的四条出路,按优先级:①补回真实圆角(首选,因为奇异多半是第 1 讲的删简造成的);②改用沿壁厚的线性化路径取值(第 5 讲);③改用网格不敏感的结构应力法——方法本体在 J5-07,本课只给指路牌;④把该点排除在评定点之外并在报告里写明理由
    • 本讲是 J5-01 那条实错的方法源头:J5-01 误区「把 CFD 的网格无关性口径直接套到尖角、开孔边、焊趾上,越加密应力越高却当成物理结论」的方法层在这里,J5-02 误区「把焊趾峰值 von Mises 送进 Coffin-Manson」同源
    • 本课的「结果可信度」只指模型内部可信度——网格是否收敛、奇异是否识别、单元与约束是否设对。与试验的相关性判据、标定与 sign-off 见 J7-01;输入散差与不确定度传播见 J7-06;输入齐套、假设登记与报告评审出口见 J7-07。模型内部错着也能靠调参跟试验对上,这正是必须先守内部可信度的理由
  7. 应力评定口径:沿壁厚线性化、Pm/Pb 分解与等效应力选用
    • 为什么不能拿云图上的单元最大应力报数:它同时受网格密度、奇异成分、后处理平均方式(节点平均 vs 单元中心、是否跨材料/跨壳面平均)三重影响,换一套网格就变,不是可评定量
    • 线性化路径怎么取:起止点落在壁厚两侧表面、沿最小壁厚方向、垂直于壁面;路径不得斜穿圆角、不得跨过材料界面(如钎料层)、不得穿过多个壁。这是 J5-01 公式 σ_m ≤ σ_yield/SF_设计 里 σ_m 的唯一合法取法——J5-01 给了结论,取法在本课
    • Pm/Pb/F 三分解的物理含义:Pm 一次薄膜应力(壁厚方向平均,与外载平衡,无自限性,超限即塑性垮塌)、Pb 一次弯曲应力(线性分量)、F 峰值应力(局部、自限,含奇异成分)
    • 评定纪律两条:一次应力评定只用 Pm 与 Pm+Pb;峰值 F 不进静强度评定、只进疲劳评定。把 F 塞进静强度会把合格件误判为超限,把 Pm+Pb 当疲劳应力用会把裂纹源漏掉
    • 一次应力与二次/自限性应力口径不同:热应力属自限性应力,它控制的是循环变形与疲劳,不是一次塑性垮塌——不能拿一次应力限值直接卡热应力。耐压三档工况与安全系数判据见 J5-01,热循环侧见 J5-02
    • von Mises 还是主应力:延性金属(变形铝、钢)用 von Mises 等效应力;低塑性/脆性材料——钎料层、压铸铝(Fe 形成的针片相使延伸率塌陷,判据靠脆性侧,见 J5-01 与 H1-04)——用最大主应力;界面剥离、密封面开裂看界面法向主应力的方向与符号,等效应力会把拉压差别抹平
    • 安全系数的口径来源要在报告里点名是谁给的:主机厂/Tier1 规范或借鉴压力容器思路的限值,随产品等级浮动。本课不给任何限值数字(按现行版本与客户规范确认),本课负责的是让 SF 分母里那个应力的取法合法
    • 疲劳寿命方法本身不在本课:S-N、ε-N、Coffin-Manson、Miner 累积损伤、平均应力修正见 J5-02 与 J5-03;接头专用的名义/热点/缺口三级应力口径、热点外插取点、FAT 等级与等效结构应力法见 J5-07。本课只交付一件事——送进那些方法的应力/应变是不是合法取来的
关键公式
Pm = (1/t)·∫₀ᵗ σ(x) dx
把云图上的总应力换成可评定的一次薄膜应力,是 J5-01「σ_m ≤ σ_yield/SF_设计」里 σ_m 的合法来源;也是判断结构会不会发生塑性垮塌的量
Pb = (6/t²)·∫₀ᵗ σ(x)·(t/2 − x) dx
把壁厚方向的弯曲分量单独拿出来,使「膜应力够但弯得太狠」这类问题能被单独判出;同时把峰值 F 从评定量里剥离
E_hg / E_int ≤ 约 5%(常用工程门槛)
判断缩减积分单元有没有失效的唯一量化出口;也是决定「加密网格」还是「提高沙漏控制刚度」的依据
σ_peak(h) ≈ C·h^(−λ),λ>0 判为奇异
把「一加密应力就升高」这件事定量成「奇异 or 尚未收敛」,据此决定是补回圆角、转线性化路径,还是转 J5-07 的结构应力法
关键概念
中面抽取等效厚度壳单元厚度偏置小特征清理准则单位制一致性自检壳单元实体单元垫片单元梁单元一阶/二阶单元剪切自锁完全积分/缩减积分沙漏(零能)模式六自由度体检过约束与欠约束对称约束适用条件RBE2 刚性连接RBE3 分布耦合惯性释放螺栓预紧截面与两步求解接触对与主从面摩擦系数上下界包络几何非线性材料非线性接触非线性真应力-真应变温度相关材料属性弹塑性本构蠕变本构超弹本构应力集中应力奇异分层网格收敛应力线性化路径一次薄膜应力 Pm一次弯曲应力 Pb峰值应力 F二次/自限性应力von Mises 等效应力最大主应力模型内部可信度
推荐工具与标准
Ansys Mechanical(静力/热-结构、线性化后处理、单元库) Abaqus/Standard(非线性、接触、预紧单元) HyperMesh / ANSA(中面抽取、几何清理、小特征处理与壳网格) Simcenter Femap / NX Nastran(RBE2/RBE3 与惯性释放的经典实现口径) 求解器自带的应力线性化路径后处理功能(用于取 Pm、Pb、F) 材料属性库/供应商数据(温度相关 E、屈服与线胀系数;取值口径见 J5-01/H1-01/H1-04)
各主机厂/Tier1 的 CAE 建模与仿真报告规范(企业 BKM,无统一编号,判据与许用值随项目浮动) 弹性应力分类与沿壁厚应力线性化(Pm/Pb/F 分解)的方法思路来自压力容器规范体系(**GB/T 150.1~150.4-2024**——2011 版已于 2025-02-01 废止、2016 复审后转为推荐性;以及 ASME 锅炉及压力容器规范弹性应力分析法一路,条款与版本按现行版本确认)——车用换热器/冷板并不直接适用,本课只借其分层框架,不借限值数表 计算固体力学的验证与确认(V&V)方法框架:ASME V&V 系列中面向计算固体力学的那一分册,具体编号与现行版本按发布方现行目录确认;本课只用其「模型内部可信度先于外部相关性」这一分层思路,外部相关性判据见 J7-01 求解器自身的理论手册与单元库文档(沙漏能定义、单元积分方案、线性化算法实现口径各家不同,判据必须按所用版本手册确认)
工程案例
某电子水泵铝支架的热循环仿真报「根部应力约 380 MPa、远超屈服,判定必裂」,而同状态件在台架热冲击耐久跑完并未开裂(本例数字均为演算假设值)。复盘查出三处方法错叠在一起:①两处安装孔用 RBE2 刚性绑到固定面,把轴向热膨胀方向一起锁死,模型里凭空造出约束反力(真实结构侧本就是一端滑动设计,设计侧口径见 I4-01);②前处理为省网格把根部 R 角删掉,尖角成为应力奇异——网格从 2 mm 加密到 0.5 mm,该点应力从 180 升到 295 MPa 且无平台,λ≈0.36>0;③报数直接取节点最大应力,未做沿壁厚线性化。整改按顺序做:补回真实 R 角 → 把 RBE2 改为法向约束+RBE3 载荷分配、释放轴向 → 沿最小壁厚布线性化路径取 Pm 与 Pm+Pb。整改后评定应力落到屈服以下,与应变片实测同量级,结论由「必裂」改为「一次静强度合格,热循环疲劳另行评价(转 J5-02)」。教训是这三处错的方向并不一致:过约束与奇异把应力推高(虚警),而若当初一阶四面体网格再粗一档,刚度偏高又会把变形算小(漏判)——所以必须逐项自检,不能靠「结果偏保守」自我安慰。
动手做
交付物 · 拿一个热管理支架或冷板局部模型(要求含至少一个承载 R 角、一处螺栓连接、一处薄壁区),做一遍完整建模自检并交付「建模自检卡」一页:①决定抽中面用壳还是直接用实体,写出判据(跨厚比、关注响应量),并列出小特征删/留清单与厚度属性(含变厚区等效厚度)的回带方式;②配三套只在关注区系统加密的网格(约 2 / 1 / 0.5 mm),对关注点做分层收敛检查,用 σ_peak~h 关系判定是应力集中还是应力奇异,并写明奇异点的处置方式(补圆角/线性化/转 J5-07/排除评定);③把一处 RBE2 改成「法向约束+RBE3 载荷分配」,对比同一温度载荷下的约束反力与根部应力变化,量化「虚假热应力」有多大;④在最小壁厚处布线性化路径,给出 Pm、Pb、Pm+Pb 与该点峰值 F 四个数,并标明哪两个可用于一次静强度评定。交付物必须包含:单元类型与积分阶次的选型理由(含沙漏能占比若用缩减积分)、约束清单与六自由度体检结论、收敛判定表、可评定应力及其取法说明;不得包含任何安全系数达标结论(限值与判据归 J5-01)。
常见误区
  • 以为网格越密越准,其实尖角、点约束、刚性连接边界、焊趾处越密应力越高且永不收敛——那是应力奇异,不是精度在提升
  • 以为删圆角只是省网格量,其实把一个有限的应力集中变成了无界的应力奇异,等于把这个评定点作废;奇异问题多数在前处理阶段就已经造成了
  • 以为四面体单元通用,其实一阶四面体在弯曲主导问题上因剪切自锁而系统性偏刚——位移偏小、频率偏高,方向偏危险(不是偏保守)
  • 以为缩减积分只是省时间,其实可能出零能沙漏模式;不查沙漏能/内能占比,就不知道结果有没有被非物理变形污染
  • 以为用对称面总是划算,其实载荷或温度场不对称时用对称约束等于人为多加了约束,造出来的热应力与约束反力都是假的
  • 以为 RBE2 和 RBE3 差不多,其实 RBE2 会在连接区局部增刚并憋死热膨胀,RBE3 只分配载荷不加刚度;用 RBE2 顶替 RBE3 是虚假热应力的常见来源
  • 以为无约束件随手锁一个点「不影响整体」,其实那个点会成为集中约束奇异并改变变形形状;自平衡载荷的自由件应当用惯性释放
  • 以为线弹性算出来超屈服「只是偏保守」,其实静强度侧偏保守而变形侧偏危险——位移、变形形状与失效位置都可能错,密封贴合与间隙包络的结论会反
  • 以为摩擦系数查个表填上就行,其实它是接触结果的敏感参数;应取上下界各算一次做包络,看结论是否翻转
  • 以为拿云图上最大那个数报安全系数最保险,其实那个数含峰值 F 与奇异成分、随网格与后处理平均方式变化,不是可评定量;一次应力评定只认 Pm 与 Pm+Pb
  • 以为 von Mises 万能,其实钎料层、压铸铝这类低塑性材料的判据靠最大主应力,用 von Mises 会把脆性开裂与界面剥离风险算漏
  • 以为「结果可信」就是跟试验对上了,其实模型内部错着也能靠调参对上某个工况——本课的可信度只到模型内部,与试验的相关性判据是 J7-01 的事
相关课题
前置:体系外前置(需自备):材料力学(应力/应变状态、强度理论、梁与薄壁件受力) · 体系外前置(需自备):有限元法基础(离散与形函数、单元刚度矩阵、位移解与应力解精度差别的概念层认识) · 体系外前置(需自备):至少一次用 Ansys Mechanical / Abaqus 一类求解器跑通静力算例的操作经验——本课讲方法与判据,不讲软件菜单路径
适合:结构/仿真(CAE)工程师——J5 模块的入口必修课,读 J5-01~J5-08 任一门之前先读本课;也适合审阅仿真报告的设计工程师,用来判断「这份报告里的应力数字能不能信」。 · 时长 约 4 小时(8 讲 + 1 次建模自检实操)

需求区 · 想听众筹

0/15 人想听
登录后想听 / 点名

想听/点名均为需求登记,不涉及任何付款(意向金仅登记不收款)。满 5 人点名 → 自动向平台专家发邀约。