TMS BOOK · ACADEMY 讲义

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

大纲的完整展开版——讲师授课蓝本 / 学员自学材料

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

课程代码 J5-06 · 板块 J 仿真分析 CAE / 结构/热应力/疲劳 FEA 时长 约 4.0 小时(8 讲) 前置 体系外前置(需自备):材料力学(应力/应变状态、强度理论、梁与薄壁件受力);体系外前置(需自备):有限元法基础(离散与形函数、单元刚度矩阵、位移解与应力解精度差别的概念层认识);体系外前置(需自备):至少一次用 Ansys Mechanical / Abaqus 一类求解器跑通静力算例的操作经验——本课讲方法与判据,不讲软件菜单路径 本讲义定位 讲师授课蓝本 / 学员自学讲义,是 J5-06 大纲的完整展开版


按任务导航

你手上的事 直接看 还要配
审一份仿真报告,判那个应力数能不能采信 第 8 讲:五类评定量各进哪一档、报数必附的取法说明 第 7 讲(那点是集中还是奇异)· 第 3 讲(约束造没造出应力)
判一个「必裂」的热应力结论是真的还是模型造的 第 3 讲:六自由度体检、RBE2/RBE3 分账、只改约束算两次的对照法 第 5 讲(该不该上非线性)· 第 1 讲(圆角是不是被删了)
定一个零件用壳还是用实体,并列小特征删留清单 第 1 讲:按响应量入口的六档判别、小特征处置的唯一判据 第 2 讲(单元类型与积分阶次跟着变)
写一份能过评审的网格收敛举证 第 7 讲:三套系统加密网格、集中与奇异的二分判据、分层收敛怎么报 第 1 讲(奇异多半是前处理删出来的)· 第 8 讲(收敛后怎么取数)
给一处螺栓连接与接触界面定仿真侧设置 第 4 讲:预紧两步法、接触四要素与必查输出、摩擦上下界包络 第 6 讲(温度相关属性怎么建)

本课不覆盖:许用值、安全系数、耐压工况与材料取值口径(→ J5-01 换热器/冷板结构强度与耐压 FEA H1-01 热管理用铝合金牌号、性能与选材 H1-04 铸造铝合金(压铸为主)的材料特性与热管理壳体选材)· 预紧力/扭矩取值与腰形孔、滑套的设计决定(→ I4-01 支架、悬置与减振结构设计)· 热疲劳寿命方法(→ J5-02 热应力、热变形与热疲劳仿真 J5-03 振动、模态与随机振动疲劳)· 接头三级应力与等效结构应力法(→ J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N)· 非金属本构标定与密封接触判据(→ J5-08 非金属件的结构仿真:各向异性本构、蠕变松弛、超弹标定与密封接触)· CFD 侧删简与 y+/GCI(→ J2-06 整车/机舱 CFD 前处理:从脏数模到可算网格(数模验收·删简·包面·计算域抽取) J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+))· 钎焊残余应力工艺态本构(→ H2-07 钎焊变形、残余应力与尺寸控制)· 试验相关性、不确定度传播与报告评审出口(→ J7-01 仿真—试验相关性验证(Correlation) J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度 J7-07 仿真任务的输入、交付与模型资产管理)· 求解器菜单、单元编号与网格生成算法(→ 求解器现行版本手册)· 屈曲失稳的方法本体(本课程体系暂无对应课 → 企业 CAE 规范与稳定性专著);以及任何一个可直接取用的许用值、安全系数、收敛阈值百分数与角点理论奇异指数表——由客户/企业规范、所引标准现行版本与求解器手册给出,⛔ 本课一个数都不给(课内算例数字全为教学假设值,⛔ 禁止照抄)。


引言:算出来了,看着也像那么回事

一份结构 FEA 报告交出去的东西很朴素:一张应力云图、一个最大值、一句「够不够」。而拿到它的人——你自己,或者下游的设计与项目——要判的只有一件事:这个数能不能拿来做决定。是照它加厚一处壁、改一版支架、推翻一个方案,还是照它签字放行。

麻烦在于,一个方法上错得离谱的模型,和一个可信的模型,交出来的东西在纸面上长得几乎一样:都有云图,都有连续的应力分布,最大值都落在一个「看起来就该应力大」的位置。分辨它们的信息不在那个数字上,在它的来路上——它是从哪条路径、哪一档分解、哪个等效口径取出来的,那个位置又还算不算一个可以评定的点。

本课只做一件事,这件事叫模型内部可信度:单元与约束设没设对、非线性与本构该不该上、网格收敛的是哪个量、那个点是不是奇异、评定应力取得合不合法。它⛔ 不包含与试验对不对得上(→ J7-01 仿真—试验相关性验证(Correlation))、输入散差怎么传播(→ J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度)、输入齐套与报告评审出口(→ J7-07 仿真任务的输入、交付与模型资产管理)。次序不能颠倒,理由很直接:模型内部错着,也能靠调参跟某一个工况的试验对上——而那样对上的模型,换一个工况就会把你带到沟里去。

钉子①:可信度是一条九环的链,而九环里只有一环会自己报错

一次结构 FEA 从几何走到报数,中间是九个环节:① 几何理想化(维度选择、中面、小特征)· ② 单元类型与积分阶次 · ③ 网格 · ④ 约束与连接 · ⑤ 载荷与边界 · ⑥ 材料与本构 · ⑦ 求解与收敛设置 · ⑧ 网格收敛判定 · ⑨ 后处理取法与报数。⚠ 这条九环链是本课归纳的责任链,⛔ 不是任何标准的原文条款。

在每一环下面问同一句话——「这一环错了,谁会告诉我?」——答案很难看:只有第 ⑦ 环会(不收敛、刚度矩阵奇异、负体积或翻转单元)。其余八环的失效形态一律是:算出来了,看着也像那么回事。

模型可信度责任链的九环流程图:横向一条链,九个方框依次为几何理想化、单元与积分阶次、网格、约束与连接、载荷与边界、材料与本构、求解与收敛设置、网格收敛判定、后处理取法与报数;每个方框正下方一行小字回答同一个问句「这一环错了,谁会告诉我?」,只有第七环求解与收敛设置写「求解器会报:不收敛、刚度矩阵奇异、负体积」并用主蓝描边,其余八环一律写「没有人」并用灰底;链条下方一条横带写「八环的失效形态一律是:算出来了、看着也像那么回事」;链条右端一个出口框写「一个应力数」,出口框旁的小框列出报数三件套:怎么取来的、收敛的是哪个量、还算不算合法评定点。九环的划分是本课归纳的责任链,不是任何标准的原文条款;框的大小与间距不代表工作量或重要性。本图为定性示意,不得据图读取任何数值。
图1 这张图要读出三件事:一,九环里只有第 ⑦ 环会自己报错,可信度⛔ 不能靠「跑通了」来证明;二,「网格」只是九环之一,⛔ 不是可信度的同义词;三,出口处那个应力数必须带着三件套一起走——怎么取来的、收敛的是哪个量、还算不算合法评定点。九环划分为本课归纳,⛔ 不是任何标准的原文条款。本图为定性示意,⛔ 不得据图读取任何数值。

这条链解释了一件常见的事:工程师对 FEA 质量的直觉几乎全部压在「网格够不够细」上——因为那是唯一一处「做得越多越安心」的地方,而且它是九环里少数你能一眼看出「做没做」的。可网格只是第 ③ 环。本课那个典型案例的三处方法错,没有一处在网格上:过约束在第 ④ 环、删掉承载圆角在第 ①、直接取节点最大值报数在第 ⑨;而复盘时第一件做的事恰恰是加密网格,加密之后那个点的应力不降反升,看起来「更准了」,实际是把一个根本不会收敛的点喂得更饱。

⇒ 由此推出本课贯穿到底的一条报数纪律:一个应力数不能单独出门,它必须带着三件套一起走——怎么取来的 · 收敛的是哪个量 · 还算不算合法评定点。三件里缺一件,这个数就还不构成工程结论:它可能是对的,但你没有理由相信它

钉子②:这些错的方向并不一致——有的把应力推高(虚警),有的把变形算小(漏判)

意识到八环静默之后,下一个念头通常是「那我就往保守里做」。这句话在本课不成立,因为「保守」是分响应量的,而同一个模型里各响应量的错误方向可以相反。

把应力推高的那一侧:用刚性连接把该自由膨胀的方向绑死,凭空造出约束反力与热应力;删掉承载路径上的 R 角,把一个有限的应力集中变成无界的应力奇异——那个点的读数从此由「你停在哪一档网格」决定,而不是由零件决定。后果是虚警:好件被判必裂,代价是白改一轮设计、白排一轮验证。

把变形算小的那一侧:一阶四面体在弯曲主导问题上因剪切自锁而系统性偏刚,位移偏小、频率偏高;缩减积分出了零能沙漏模式却没查沙漏能占比,结果被非物理变形污染;线弹性算到越过屈服仍按线弹性报数,位移、变形形状与失效位置都可能错。后果是漏判,而且漏掉的恰恰是最贵的那一类判断——密封面贴合量、间隙包络、卡边剩余压缩量,这些结论会直接反过来。

⇒ 所以「反正结果偏保守」这句话,对应力侧成立,对位移与变形形状侧不成立,⛔ 不能用一句话覆盖全部响应量。说「保守」之前先补一个宾语:对哪个响应量保守?而你这次要交的结论挂在哪一个响应量上,决定了哪个方向的错才是致命的。

反钉子:线性化不是「把过高的应力合法地降下来」,它是把一个数拆开归类

学完「云图上的单元最大值不是可评定量」之后,最容易顺出来的下一句是:「那就沿壁厚做一次应力线性化——这就是把过高的应力合法地降下来的办法。」 这句话在纸面上几乎无懈可击,而它是本课要拆掉的那一句。

它对应的那个具体动作是这样的:某个点数值太大、报告不好交代 ⇒ 沿壁厚拉一条路径做线性化 ⇒ 线性化后的一次应力落到屈服以下 ⇒ 写一句「一次静强度合格」结案。而当这个点的圆角正是前处理时自己删掉的、它已经是一个应力奇异时,这个动作看上去把问题解决了,实际上只是把病灶盖住了。

第一处要拆开:数变小不是因为它被「降下来」了,是因为剥掉的那部分本来就不进这一档。 线性化做的是分解与归类:把路径上的总应力拆成沿壁厚平均的一次总体薄膜应力 Pm、沿壁厚的线性分量一次弯曲应力 Pb,以及剩下的峰值应力 F。一次静强度只用 Pm 与 Pm+Pb,因为它判的是整壁的塑性垮塌,而 F 是局部的、自限的。教学算例(教学假设值,不对应任何平台,⛔ 禁止照抄取用):壁厚 t=3 mm 的路径上,两表面应力 220.0 MPa 与 40.0 MPa 呈线性分布,起点表面的云图读数 260.0 MPa ⇒ Pm=(220.0+40.0)/2=130.0 MPa、Pb=(220.0−40.0)/2=90.0 MPa、Pm+Pb=220.0 MPa、F=260.0−220.0=40.0 MPa。被剥掉的那 40.0 MPa 没有消失,它换一条评定链走——疲劳裂纹恰恰从表面峰值处起,把 F 丢掉等于把裂纹源丢掉。完整的规则是「F 换一条链走」,⛔ 不是「F 可以不管」。

第二处要拆开:在奇异点上做线性化,⛔ 不恢复这个点的局部应力与疲劳评定资格。 奇异处有四条出路,而且有优先级:① 补回真实圆角(首选——奇异多半是前处理删出来的,补回去等于消除病因,代价是重做网格);② 改用沿壁厚的线性化路径取值(给得出合法的一次应力评定,但那个点从此不能再作局部应力与疲劳评定点);③ 改用网格不敏感的结构应力法(方法本体 → J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N,本课只给指路牌);④ 把该点排除在评定点之外,并在报告里写明理由与影响范围。直接跳到 ② 并宣布合格,就是这一节点名的那个做反动作——报告看起来合格了,而删简的真实代价被掩盖了,代价还在零件上。

⇒ 凡是用线性化取出来的数,报告里必须跟着两句:这条路径穿的是哪一段壁(沿最小壁厚、垂直壁面,⛔ 不斜穿圆角、⛔ 不跨材料界面),以及这个点是集中还是奇异;若是奇异,明写一句「本点⛔ 不再作局部应力与疲劳评定点」。这两句写不出来,那个 Pm 就还不是一个可评定量。

全课八讲怎么排

第 1 讲几何理想化:维度选择有六档而不是壳与实体的二选一,中面在三处必然断裂、要手工缝合,厚度与偏置方向必须显式回带;小特征的判据只有一条——它在不在关注响应量的传力路径上、是不是评定点,⛔ 不是尺寸大小。第 2 讲单元与积分阶次:剪切自锁与零能沙漏两个机理,以及用沙漏能占比反查选型成不成立(分母取域必须声明)。第 3 讲约束与连接:六自由度体检、刚性连接与分布耦合分账、对称约束的四条件(温度场不对称就不能切对称面)、惯性释放的三条件与反力残差校验——虚假热应力全部长在这一讲。第 4 讲螺栓预紧与接触界面:预紧两步法的时序、接触对四要素与必须被检查的输出、摩擦系数取上下界做包络。第 5 讲非线性分界:用线性解的三个比值决定要不要升级,并说清这三个比值不覆盖哪一整类——薄壁受压失稳的特征恰恰是应力还远低于屈服、位移也不大,而结构已经垮了。第 6 讲材料本构:工程应力应变必须转成真应力应变、弹塑性三档各能回答什么、温度相关属性要覆盖到工况上下界之外且⛔ 不许靠外推撑边界点。第 7 讲集中还是奇异:三套系统加密网格给出二分判据,以及结构侧是分层收敛——报「收敛」必须说清收敛的是哪个量。第 8 讲评定口径:线性化路径的合法取法、五类评定量各进哪一档、von Mises 还是最大主应力。

⇒ 本课最后交出去的不是一个数,是一页建模自检卡(本课归纳):壳与实体的判定及小特征删留清单 · 收敛判定表(写明收敛的是哪个量) · 一处刚性连接改成「法向约束+分布耦合载荷分配」前后的对照(那个差值就是这个模型的虚假热应力) · 一组合法取来的评定应力及其取法说明。⛔ 卡上不写任何「安全系数满足要求」的结论——那不是本课交付的东西。

全课统一的四条口径

放在这里一次说清,后面各讲不再重复:

  • 符号撞名先钉住:本课与传热类课程共用一批字母而含义完全不同——t 是壁厚(保载时长写 t_hold)· h 是单元特征尺寸,⛔ 不是对流换热系数 · λ 是奇异指数,⛔ 不是导热系数 · Q 是二次应力(热量写 Q_heat)· F 是峰值应力,力一律带下标(预紧力 F_pre、约束反力 F_react)· μ 是摩擦系数,⛔ 不是动力粘度。完整的一张表在第 1 讲开头,全课不改。
  • 单位制一律 mm-t-s-MPa:长度 mm、质量 t、时间 s、力 N、应力 MPa、能量 mJ、密度 t/mm³。换算本身是恒等式,但它是量纲级错误的高发点——只要有一个属性沿用了 SI 数值,整份结果错的是数量级,而⛔ 不报错。
  • 本课出现的全部算例数字都是教学假设值,不对应任何平台、任何牌号、任何供应商产品,⛔ 禁止照抄取用。它们的算术逐点可复算(纯代数与恒等式),可复算的是关系,⛔ 不是数值本身。
  • 两个短语禁止裸写:①「网格已收敛」——必须写成「哪一个量、在哪个加密序列上、相邻两套的相对变化是多少、且这个变化在衰减」;②「等效应力」——必须写明是 von Mises 还是最大主应力口径,第 8 讲整讲的结论都挂在这个区分上。

读完这一节,你今天就能用的东西

拿到任何一个应力数(别人报的,或你自己算的),逐个问这四句,答不上来的数不能进结论:

  1. 怎么取来的?哪条路径、哪一档分解、哪个等效口径(回答「等效应力」不算答案)。
  2. 收敛的是哪个量?位移收敛与峰值应力收敛是两码事,拿前者证明后者,是拿一个量的结论替另一个量背书。
  3. 这个点还算不算合法评定点?是有限的应力集中,还是加密只会越来越大的应力奇异;若是奇异,四条出路走的是哪一条、代价写没写。
  4. 这次要交的结论挂在哪个响应量上?应力、位移、变形形状、失效位置——同一个模型里它们的错误方向可以相反,「偏保守」不带宾语就是一句空话。

这四句就是钉子①、钉子②与反钉子的现场形态。后面八讲要做的,是把每一句背后的做法讲到能照着执行:第 1~4 讲与第 6 讲把「来路」上的五个环节逐个讲清(第 1 问的前半);第 7 讲管第 2、3 两问;第 8 讲管第 1 问的最后一段——取法本身,以及第 3 问判为奇异之后怎么收场;第 4 问那笔方向账由第 2、5、6 三讲各出一段。

第 1 讲 几何理想化:责任链的第一环,也是最不会喊的那一环

一份结构 FEA 报告递到手上,你真正要判的只有一句话:这个应力数能不能用。本课把这句话拆成一条有九个环节的责任链,而第 1 讲讲的是第一环——几何被简化成了什么样子:用壳还是用实体、中面怎么抽、抽完的厚度属性怎么带回去、哪些小特征被删掉了、几何缺陷与单位制有没有被验过。

把它排在最前面,不是因为它最难,而是因为它同时占了两条性质:最不会喊,而且不可逆。删掉一个承载圆角、抽错一张中面、把一条另一套单位制的密度拖进模型——求解器全都照常算完,云图照常出,颜色照常漂亮;等到第 7 讲发现某个点怎么加密都不收敛,回头要重做的往往是整套网格。⇒ 本讲每一节其实都在回答同一个问题的一个分支:这一步做错了,谁会告诉我?答「没有人」的那几步,就是必须靠人主动去核的地方。

本讲交付什么:一张全课通用的撞名符号表与一套单位制;一张射程表(哪些东西本课只给去向);一张九环自检卡;维度选择的六档判据与各自选错的方向;中面抽取的三处必然断裂与缝合判据;厚度与偏置属性的回带三问;变厚过渡区等效厚度的两档取法与取错的量化代价;小特征的唯一判据与五种处置;删圆角的代价与删简的影响边界;几何缺陷与单位制这两类「输入本身就错了」的自检。

本讲不交付什么:任何前处理器的菜单路径与自动清理参数;网格尺寸怎么定、网格质量指标怎么看(收敛判定本身在第 7 讲);CFD 侧的数模验收、破面缝合、包面与计算域抽取(去向 J2-06 整车/机舱 CFD 前处理:从脏数模到可算网格(数模验收·删简·包面·计算域抽取));材料属性取哪个数(本课只管怎么把已定的属性正确建进模型,取值口径去向 J5-01 换热器/冷板结构强度与耐压 FEA · H1-01 热管理用铝合金牌号、性能与选材 · H1-04 铸造铝合金(压铸为主)的材料特性与热管理壳体选材)。

路线:先钉住符号与单位制(1.1)、划清射程边界(1.2)、摆出九环自检卡(1.3);然后进几何理想化本身——维度选择(1.4~1.5)、中面与厚度属性(1.6~1.8)、小特征的判据与处置(1.9~1.10)、删圆角与删简的代价(1.11~1.13);最后两节是「输入本身就错了」的那两类(1.14~1.15)。

1.1 先钉住九组撞名符号与一个单位制——本讲开头一次钉死,全课不改

本课同时用到固体力学、有限元实现与热管理三套符号,撞名极其密集。下面这张表在本讲开头一次钉死,全课不改,后面各讲不再重复声明。

本课符号 在本课的含义 你多半会带进来的另一个含义 那个含义在本课怎么写
t 壁厚(也是应力线性化路径的长度),mm 时间 t_hold(保载 / 蠕变时长,h 或 s)
h 单元特征尺寸(网格尺寸),mm 对流换热系数 本课全课不出现换热系数
F 峰值应力,MPa F_pre 预紧力 · F_n 法向力 · F_react 约束反力(N)
Q 二次应力(自限性,热应力属这一类),MPa 热量 / 换热量 Q_heat(W)
E 弹性模量,MPa 能量 E_hg 沙漏能 · E_int 内能(mJ)
λ 奇异指数,无量纲 导热系数 本课全课不出现导热系数
μ 摩擦系数,无量纲 动力粘度 本课全课不出现粘度
T 温度,°C(温差写 ΔT,单位 K) 拧紧扭矩 T_tq(N·m)
K 拧紧系数,无量纲 应力集中系数 K_t

为什么这张表值得占一节而不是一个脚注:本课读者是热管理背景的结构 / 仿真工程师。h、λ、Q、μ 这四个字母在你的肌肉记忆里全是传热量(换热系数、导热系数、热量、粘度),而在结构评定语境里它们是四个完全不同的东西。更要命的是 F:同一段文字里既有「峰值应力 F」又有「预紧力 F」,两者量纲不同,而读者可以毫无障碍地把它们代进同一个式子。

三条配套纪律:

  • ⛔ 全课禁止裸写 F 表示力——力一律带下标(F_pre / F_n / F_react)。
  • ⛔ 全课禁止裸写「等效应力」四个字,必须写明是 von Mises 等效应力(σ_eq_vM)还是最大主应力(σ_1)。第 8 讲整讲的结论都挂在这个区分上。
  • Pm 与 PL 是两个量:一次总体薄膜应力与一次局部薄膜应力。两者都是一次应力,但成立的范围不同,⛔ 不得混成一个(第 8 讲展开)。

单位制:全课统一取 mm-t-s-MPa 制(长度 mm、质量 t、时间 s、力 N、应力 MPa、能量 mJ、密度 t/mm³)。它为什么值得和符号表并列、以及它错了会怎样,放在本讲最后一节(1.15)。

⚠ 本课全部教学算例都是教学假设值:不对应任何平台、任何牌号、任何供应商产品,⛔ 禁止照抄取用(本讲用到的是算例 D,另引用第 7 讲的算例 A 与 B 作对照)。它们的算术可以逐步复算(纯代数与恒等式)——可复算的是关系,⛔ 不是数值本身。

今天就能做的一件事。打开你手上正在写的那份仿真报告,检索 F、Q、h、λ 四个字母,看同一份文件里它们有没有被用在两个含义上。这是一次零成本的自检,而报告发出去之后,读它的人不会有机会来问你「这个 F 是哪个 F」——他会直接代进式子。

1.2 本课的射程:八件事只给去向,⛔ 不给内容

是什么。本课交付的每一个判据式,分母或输入里都有一项来自别处。下面这张表把它们一次列清,本课一律只给去向、⛔ 不给内容

你在本课会遇到、而本课不给的东西 去向
许用值、安全系数、三档耐压工况 J5-01 换热器/冷板结构强度与耐压 FEA 与客户 / 企业规范
材料取值口径(钎后态、压铸本体与单铸试棒之别、钎缝折减) J5-01 换热器/冷板结构强度与耐压 FEA · H1-01 热管理用铝合金牌号、性能与选材 · H1-04 铸造铝合金(压铸为主)的材料特性与热管理壳体选材
螺栓预紧力与拧紧扭矩的取值;真实结构该不该留腰形孔 / 滑套 I4-01 支架、悬置与减振结构设计
热应力工况与热疲劳寿命方法 J5-02 热应力、热变形与热疲劳仿真
接头的三级应力口径与等效结构应力法 J5-07 钎缝/焊缝接头的疲劳评价方法:三级应力口径、等效结构应力法与自建 S-N
非金属件的各向异性 / 蠕变 / 超弹标定与密封接触判据 J5-08 非金属件的结构仿真:各向异性本构、蠕变松弛、超弹标定与密封接触
CFD 侧的删简、包面、计算域抽取与网格无关性口径 J2-06 整车/机舱 CFD 前处理:从脏数模到可算网格(数模验收·删简·包面·计算域抽取) · J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+)
钎焊残余应力的工艺态本构 H2-07 钎焊变形、残余应力与尺寸控制

另有三层属于模型外部可信度,同样只给去向:与试验的相关性判据与放行 ⇒ J7-01 仿真—试验相关性验证(Correlation);输入散差与不确定度传播 ⇒ J7-06 仿真不确定度量化(UQ)与稳健/可靠性设计:从确定性单点到带分布的裕度;输入齐套、假设登记与报告评审出口 ⇒ J7-07 仿真任务的输入、交付与模型资产管理

本课的「结果可信度」只指模型内部可信度——单元与约束设没设对、非线性与本构该不该上、网格收没收敛、奇异识没识别出来、评定应力是不是合法取出来的。⚠ 为什么这个边界要划在这里,而不是干脆一路讲到跟试验对上:因为模型内部错着,也能靠调参跟某个工况的试验对上。一个过约束的模型配一个偏软的材料常数,完全可能在某个测点上跟应变片吻合得很好——而它在别的工况、别的测点上会随机地错。⇒ 内部可信度是外部相关性的前置条件,顺序反过来就变成了拟合。

为什么必须先写边界再讲内容。因为不先说清「从哪取、取不到怎么办」,读者会默认把取不到的那一项当成不存在——而一个被默认掉的输入不会触发任何检查:它不在假设清单里,不在评审问题单上,也不在报告的任何一页上。

易错点。① 把「去向某课」读成「那门课会给我一个数」——去向给的是方法与口径,平台相关量在那边同样不给数;② 把边界外的东西顺手补一个「常见值」进报告,而不标注它是谁给的;③ 反过来,因为「本课不给限值」就连评定应力的取法也不写——取法恰恰是本课的交付物,本课负责的是让安全系数分母里那个应力取得合法

今天就能做的一件事。在报告的输入表上给每一行加一列「谁给的」,只填三种值:本项目规范 / 设计侧 / 我自己定的。第三种如果不是空的,那几行就是这份报告真正的风险所在——它们没有出处,也就没有人会来复核。

1.3 责任链九环:逐环问一句「这一环错了,谁会告诉我?」

是什么。一次结构 FEA 拆成九环,每一环下面挂同一个问句:

内容 这一环错了,谁会告诉我?
① 几何理想化 维度选择、中面、小特征 ⛔ 没有人
② 单元类型与积分阶次 阶次、完全 / 缩减积分 ⛔ 没有人
③ 网格 尺寸、过渡、形状质量 ⛔ 没有人(质量指标要人去看)
④ 约束与连接 约束、刚性连接、接触 ⛔ 没有人
⑤ 载荷与边界 含温度场的施加方式 ⛔ 没有人
⑥ 材料与本构 属性、本构类别、温度相关性 ⛔ 没有人
⑦ 求解与收敛设置 步长、收敛准则、求解器设置 会报:不收敛 / 刚度矩阵奇异 / 负体积
⑧ 网格收敛判定 分层收敛、奇异识别 ⛔ 没有人(不做也不会报)
⑨ 后处理取法与报数 平均方式、路径、分类 ⛔ 没有人

九环里只有第 ⑦ 环会自己报错。其余八环的失效形态一律是同一句:算出来了,看着也像那么回事。

⚠ 这九环的划分是本课归纳的责任链,⛔ 不是任何标准的原文条款,也⛔ 不是本课程体系既有的分环——引用时要把这条来路一并写上。本课覆盖 ①②④⑥⑧⑨ 六环,并对 ③⑤⑦ 各给一句自检问句。

为什么这是全课的骨架。因为多数人对「FEA 质量」的印象几乎全部压在「网格够不够细」这一环上,而本课典型案例里的三处错,没有一处在网格上:刚性连接过约束在第 ④ 环,删掉根部圆角在第 ①,直接取节点最大值报数在第 ⑨。更值得记住的是复盘时人们的第一反应——先去加密网格。加密之后那一点的应力从 180.0 升到 295.0 MPa(教学假设值,见 1.11),看上去「更准了」,实际上是把一个奇异点喂得更饱。

易错点。① 把九环当成一份流程顺序背下来,而不是当成九个「谁报错」的问句——它的价值全在那一列答案上;② 认为第 ⑦ 环收敛了就说明前六环没问题——收敛只说明这组方程被解出来了,⛔ 不说明解的是对的那组方程;③ 把「自检」理解成再跑一遍——同一个错模型跑两遍还是那个错,而两次结果一致会让人更放心。

今天就能做的一件事。把这九行抄成一张卡贴在报告模板的第一页,交付前逐环填一句「我怎么知道这一环没错」。⚠ 判据不是填满九行,是那八行里有没有出现「因为跑通了」「因为收敛了」这类回答——出现了,就说明那一环实际上没有被检查过。

1.4 壳还是实体:三条判据同时问,任何一条不成立就往实体走

是什么。三条判据,同时成立才用壳

  1. 壁厚方向的应力可以忽略——你关心的不是壁厚方向的梯度;
  2. 特征跨度与壁厚之比在 10 倍以上量级(行业典型判据,⚠ 它是量级,⛔ 不是一个分界开关,8 与 12 之间不存在一个开关);
  3. 你要读的是膜 + 弯响应

反过来,只要出现下面任何一样,就往实体走:需要壁厚方向的温度或应力梯度、三向应力状态、局部接触压力分布。⚠ 三条里第 ① 与第 ③ 是定性的,而且优先级更高:跨厚比过关、但你要看壁厚方向的梯度时,仍然用实体。

为什么判据要挂在响应量上,而不是零件类型上。同一个钣金支架,算模态与整体刚度时用壳完全够;算螺栓孔周的局部承压时必须实体。零件没变,问题变了,答案就变。这是本课反复出现的一条主线:几乎所有建模选择的入口都是「我要读的是哪个量」,⛔ 不是「这个件长什么样」。

工程量级。跨厚比 10 倍以上量级用壳(行业典型区间,⛔ 不作可承诺值);需要看壁厚方向梯度时一律实体,且壁厚方向至少要有能表达那个梯度形状的层数——层数由所解的梯度形状定,本课⛔ 不给固定值。可执行的判据是:把壁厚方向的层数加一倍重算一次,看你要读的那个梯度形状变不变;不变,说明层数够了。

易错点。① 把判据记成「薄就用壳」——薄而要读接触压力时,壳给不出;② 忽略三向应力状态(厚壁、约束密集区、压装配合处的应力三轴度高,壳的平面应力假设在那里直接失效);③ 用壳算完之后再去读「壁厚方向的应力分布」——那张分布是假设出来的,不是算出来的:壳的本构里根本没有那个方向的自由度,你读到的是后处理按线性假设铺出来的结果。

今天就能做的一件事。在选单元之前,先把「我这次要报的量」写成一行字(不是「做个强度分析」,而是「根部最小壁厚处的一次薄膜应力」或者「盖板的一阶模态频率」)。这一行写不出来,后面每一个建模选择都没有判据可依——而它多半也说明这次分析的任务书本身就是空的。

1.5 维度选择有六档,⛔ 不是「壳还是实体」两档——被整类漏掉的是轴对称、平面应力 / 应变与子模型

是什么。几何理想化的维度选择,全集是六档,另加一支性质不同的接口:

什么时候用它 选错的方向
① 三维实体 要壁厚方向梯度、三向应力状态、局部接触压力 网格量与算力代价
② 壳 / 中面 膜 + 弯响应,跨厚比在 10 倍以上量级 壁厚方向的分布是假设出来的
③ 梁 长细支撑,只要轴力与弯矩;螺栓的等效替代 局部承压与孔周应力读不出来
轴对称 ⚠ 本课归纳补上 储液件、管接头、压装衬套这类回转体 载荷与约束必须同样轴对称,否则整个假设不成立
平面应力 / 平面应变 ⚠ 本课归纳补上 薄板面内问题取平面应力;长条截面取平面应变 方向是系统性的:平面应变把面外约束当成无限刚,应力偏高;平面应力反之
子模型 ⚠ 本课归纳补上 关注区要加密,而整模型算不起 切割边界取得不对,局部解跟着错
壳—实体过渡(接口,⛔ 不是维度)⚠ 本课归纳补上 一个模型里两种维度并存时的接缝 直接共节点会形成一条假铰

⚠ 后四支是本课归纳补上的,⛔ 不是本课程体系既有的分档,也⛔ 不是任何标准既有的分类。

为什么必须补。只沿「壳还是实体」这条二选一线索去枚举,漏掉的不是「少一条」,是整类不出现——清单看上去仍然是完整的,因为没有任何空位提示你少了什么。其中漏掉子模型的后果最实在:读者会默认「关注区加密」等于「整模型加密」,于是第 7 讲那套三套网格的收敛检查在一个真实零件上算不起,然后这件事就不做了。而子模型(用粗网格全局解的位移场当局部模型的边界条件)正是它在工程上唯一做得起的做法。

漏掉平面应力 / 平面应变则是方向性的错。这两个选项在菜单上像同一个「二维简化」开关,实际是两个方向相反的假设:一个把面外完全放开,一个把面外锁成无限刚。

工程量级。子模型可以把关注区的网格尺寸降一到两个数量级而总自由度不变;轴对称件用轴对称单元可以把网格量降一到两个数量级——两句都是量级判断,具体倍数取决于零件,⛔ 不作可承诺值。轴对称的代价是载荷与约束必须同样轴对称(⚠ 与第 3 讲对称约束的条件同源:几何对称只是其中一条)。

易错点。① 子模型的切割边界离关注区太近——全局解的位移场在那里本身就不准,切割面要落在应力梯度平缓处;② 壳与实体直接共节点——壳有转动自由度、实体没有,接缝处会形成一条假铰,传力路径被静默削弱;③ 轴对称模型上加了非轴对称载荷还继续用。

今天就能做的一件事。拿你手上这个零件,把六档挨个问一遍「这一档为什么不适用」。⚠ 判据不是「哪一档最合适」,是逐档说得出排除理由——说不出理由的那一档,往往就是你根本没想到的那条路,而它多半正是你觉得「这个模型算不动」的出口。

几何理想化维度选择图。上半部横排六个方框,各带文字标签与一行适用场景:三维实体、壳与中面、梁、轴对称、平面应力与平面应变、子模型;后三个框右上角标注为本课归纳补上的分档。下半部是一棵判别树,入口写「我要读的是哪个响应量」,第一问「需要壁厚方向的梯度、三向应力或局部接触压力吗」,答是走向实体;答否进第二问「特征跨度与壁厚之比在 10 倍以上量级吗」,答是走向壳,答否回到实体。主干旁三条支路各带标签:几何是回转体且载荷与约束同样轴对称,走向轴对称;长细杆件只要轴力与弯矩,走向梁;关注区要加密而整模型算不起,走向子模型。树右一个警示框写:壳与实体直接共节点会形成一条假铰,转动自由度与平动自由度不匹配。本图为定性示意,图中 10 倍是量级判据而不是分界开关,不得据图读取任何数值。
图2 这张图要读出三件事:一,选择有六档,⛔ 不是壳与实体的二选一;二,判别的入口是「我要读哪个响应量」,⛔ 不是零件叫什么名字;三,子模型是「关注区加密」这件事在真实零件上唯一算得起的那条路——漏掉它,第 7 讲的收敛检查就变成一件做不了的事。★ 右上角标「本课归纳补上」的四支,⛔ 不是本课程体系既有的分档。⛔ 图中 10 倍为量级判据、⛔ 不是分界开关;本图为定性示意,⛔ 不得据图读取任何数值。

1.6 中面抽取:断裂只发生在三个可预期的位置,而缝合的判据是节点共用,⛔ 不是渲染看着连着

是什么。三件事:

等厚区自动抽取可用,直接采信; ② 变厚过渡区两条路——按分区分别赋厚度(首选),或就近取等效厚度(等效原则按响应量定,见 1.8); ③ 中面在三处会断裂:交叉(T 型 / 十字加强筋)、变厚、圆角。断裂处必须手工缝合成连续曲面。

为什么缝合不是美观问题。中面断开的地方传力路径就断了——两片壳没有共用节点,载荷传不过去,而求解器不会报错:它只是算出一个「这一支不受力」的结果。这个结果在云图上完全正常,只是应力分布跟你预期的不一样,而人对「跟预期不一样」的第一反应通常是去改预期,不是去查拓扑。

工程量级。抽中面在实际零件上是人工工时的大头,一个带加强筋的钣金或压铸支架通常有十几处断裂要缝合(量级观察,⛔ 不是指标)。

易错点。① 只看渲染出来的中面「连着」就采信——渲染连着与拓扑连着是两回事,判据是查节点是否共用;② 变厚过渡区直接交给工具自动抽——工具会在过渡处生成一段斜的中面并给一个自己算的厚度,而那个厚度是几何平均,不是响应量意义上的等效(下面两节就是治这一条的);③ 缝合时把两片壳硬拉到同一位置,制造出一段几何偏差——弯曲臂长跟着变了,而这个变化不会有任何提示。

今天就能做的一件事。抽完中面后,把模型的自由边显示出来看一遍。⚠ 判据不是「自由边多不多」,是每一条自由边你都说得出它为什么应该自由(真实的开口边界、对称面边、外缘)——说不出来的那几条,就是没缝上的地方。这一步的成本是几分钟,漏掉它的代价是一整支结构不受力。

中面断裂与缝合对照图。左列三个剖面小图,各带文字标签,画出中面在三处断开的形态:T 型与十字加强筋交叉处,三片中面互不相交并留下三个自由边;变厚过渡处,两段中面错开半个厚度差;圆角处,中面在拐角断成两截。右列三个对应小图画出缝合后的连续曲面,缝合线加粗并标注节点共用。每一对之间一个箭头,标注缝合前这一支不受力而求解器不报错,缝合后传力路径连通。底部一行提示:渲染上看着连着不等于拓扑上连着,判据是节点是否共用。本图为定性示意,剖面为示意比例,不得据图读取任何数值。
图3 这张图要读出三件事:一,中面断裂只发生在三个可预期的位置(交叉、变厚、圆角),⛔ 不是随机的——所以它是可以逐处去查的;二,断裂处的后果是传力路径丢失而无人报错,它不表现为算不出来,而表现为「这一支不受力」;三,缝没缝上的判据是节点是否共用,⛔ 不是渲染看起来连着。⛔ 本图为定性示意、不含任何尺寸与数值,⛔ 不得据图读取任何数值。

1.7 厚度与偏置方向必须显式回带——偏置反了,位移看不出差别,而「哪一面受拉」整体反过来

是什么。壳单元只是一张面。厚度、偏置方向(中面 / 上表面 / 下表面)、上下表面的对应关系三样全部是后加的属性,必须逐分区显式赋回去。

为什么它是「静默错误」的典型。两种失败,隐蔽程度还不一样:

  • 漏带厚度:求解器会用一个默认值算下去(不少前处理器默认 1 mm),云图照常出。这一档还算有救——默认值离谱时,结果会离谱到被人看出来。
  • 偏置方向反了:更隐蔽。模型的总刚度几乎不变,位移看不出差别,而壳的上下表面被互换,弯曲应力的符号、以及它出现在哪一侧,全部反过来。而结构评定恰恰关心「哪一面受拉」——裂纹从受拉面起。

工程量级。偏置量是半个壁厚,对一块 2 mm 板是 1 mm。这个量对整体刚度可以忽略,对装配间隙、接触面位置与弯曲臂长不能忽略——它正好落在「刚度上无所谓、几何上有所谓」这条缝里,而人们习惯用刚度去判断一个几何量重不重要。

易错点。① 只检查「厚度赋了没有」,不检查「赋对了没有」——分区错位比漏赋更常见,而漏赋还有机会被离谱的结果暴露,分区错位不会;② 靠肉眼看渲染出来的实体化显示判断偏置方向——必须核单元法向;③ 用壳算完之后去读「表面应力」而不说清读的是上表面还是下表面。⚠ 第 ③ 条与第 8 讲那条「Pb 的符号随线性化路径起点翻转」是同一类错:一个带方向的量在报数时没有声明方向,于是两个人可以从同一个模型报出符号相反的结论,而两人都没算错。

今天就能做的一件事。交付前逐分区做一次三问自检:这一区的厚度是多少?单元法向指向哪一面?我报的表面应力读的是哪一面?三问缺任何一个,弯曲侧的结论就没有方向——而没有方向的弯曲结论,在评定里等于没有结论。

壳单元厚度属性与偏置方向示意图。中央一条粗线代表中面,上下各一条虚线代表上表面与下表面,三者之间标出偏置量为半个壁厚并带标签。左侧画偏置方向正确时的弯曲应力分布,上表面标受拉面、下表面标受压面;右侧画偏置方向反了之后的同一张图,上下两个标签互换并以警示色标出,旁边一条注解线写:模型总刚度几乎不变,位移看不出差别。底部一个方框列出三件必须显式回带的属性,各带文字标签:厚度、偏置方向、上下表面对应关系。图上另标出单元法向箭头。本图为定性示意,应力分布只画符号与位置关系、不含量值刻度,不得据图读取任何应力值或厚度值。
图4 这张图要读出三件事:一,偏置反了会让「哪一面受拉」整体反过来,而裂纹从受拉面起;二,这个错在位移与总刚度上几乎不留痕迹 ⇒ ⛔ 不能靠看位移发现它,只能去核单元法向;三,厚度、偏置方向、上下表面对应关系三件都要显式回带,漏哪一件都不报错。⛔ 本图为定性示意,⛔ 不得据图读取任何应力值或厚度值。

1.8 变厚过渡区的等效厚度按响应量分两档,⛔ 不得混用

是什么。两档,⛔ 不得混用:

  • 刚度 / 位移 / 模态类问题 ⇒ 按抗弯等效取。抗弯刚度正比于 t³,于是按 t³ 做面积加权:t_eq = (Σ f_i · t_i³)^(1/3),其中 f_i 是各段厚度所占的面积份额。⚠ 这个等效式是本课归纳的——由「抗弯刚度正比于 t³」与面积加权直接得出,⛔ 不是任何标准既有的公式;报告里用它,要把这条来路一并写上。
  • 一次强度类问题 ⇒ 按最小壁厚保守取。理由:一次薄膜应力与外载平衡,它随真正承载的那个壁厚变,⛔ 不随一个加权平均变。

⚠ 连带一条符号纪律:⛔ 不得裸写 t_eq 而不说是哪一档——两档给出的是两个不同的数(下面就是),一个不带档位标记的 t_eq 在报告里等于没定义。

为什么两档不能混。因为两者的方向相反:t³ 加权会被厚的那一段拉上去,所以抗弯等效厚度总是大于最小壁厚;拿它去算一次薄膜应力,就会把 Pm 算小。

工程量级(教学算例 D)。⚠ 教学假设值:不对应任何平台、任何牌号、任何供应商产品,⛔ 禁止照抄取用。过渡区一半 2 mm、一半 3 mm,面积份额各 0.5。逐步复算:

  • 抗弯等效:t_eq = (0.5×2³ + 0.5×3³)^(1/3) = (0.5×8 + 0.5×27)^(1/3) = 17.5^(1/3) = 2.5962 mm
  • 算术平均:2.5000 mm
  • 最小壁厚:2.0000 mm

⇒ 拿抗弯等效厚度去做一次强度评定,会把 Pm 低估 29.81%(2.5962/2.0000 − 1)。

⚠ 这个 29.81% 不是「差一点」。它直接落在安全系数的分母上,而分子(许用值)由别处给(本课不给)——于是报告里凭空多出约 1.3 倍的裕度,而它长得跟真的裕度一模一样:中间没有任何一个量会显得反常。

易错点。① 用算术平均 2.5000 mm——它既不是抗弯等效、也不是保守取法,两头不靠;② 同一块过渡区被两个响应量共用(既要算模态又要评强度)——这时必须分区给真实厚度,⛔ 不许取任何等效值,因为一个值不可能同时满足两个方向相反的等效原则;③ 把这个等效式当成标准公式去引用。

今天就能做的一件事。把模型里所有用了「等效厚度」的分区列一张表,每行填三样:这一区服务于哪个响应量 · 用的是哪一档等效 · 该区的最小壁厚是多少。第二列与第一列对不上的行,直接改成分区给真实厚度——这比事后解释那 29.81% 从哪来便宜得多。

变厚过渡区三种厚度取法的对照条形图。三条并列的横条各带文字标签与数值:最小壁厚 2.0000 毫米,标注为一次强度类的保守取法;算术平均 2.5000 毫米,标注为两头不靠;抗弯等效 2.5962 毫米,标注为刚度、位移与模态类的取法。右侧一列标出各自相对最小壁厚的比值 1.0000、1.2500、1.2981,抗弯等效那一行另标注:用它做一次强度评定会把一次薄膜应力低估 29.81%。图上另设一个小框写出等效式 t_eq 等于面积份额加权的厚度立方和再开立方根,并给出本例代入过程:0.5 乘 2 的立方加 0.5 乘 3 的立方等于 17.5,开立方根得 2.5962,并标注该等效式为本课归纳、不是任何标准既有。全图顶部标注:教学假设值,过渡区一半 2 毫米、一半 3 毫米,不对应任何平台,禁止照抄取用。
图5 这张图要读出三件事:一,同一块过渡区,三种取法给出三个不同的数,⛔ 不存在一个通用的等效厚度;二,取哪一个由响应量决定——刚度类用抗弯等效,一次强度类用最小壁厚;三,取错的方向是把一次薄膜应力低估 29.81%,⛔ 不是「差一点」,而它会以「裕度」的面目出现在报告里。⚠ 图中三个厚度值为教学假设值、不对应任何平台、⛔ 禁止照抄取用;等效式为本课归纳,⛔ 不是任何标准既有。

1.9 小特征清理的判据只有一条,而且⛔ 不挂在尺寸上

是什么。分级判据只有一句:这个特征在不在你关注的那个响应量的传力路径上、是不是评定点。由它直接得出三条处置:

  • 承载路径上的圆角,一律不删
  • 非承载的 logo、工艺凸台、装配标记,先删
  • 小孔、卡扣,按是否传载分别处置。

为什么判据不能挂在尺寸上。一个 R1 的根部承载圆角,与一个 R1 的装饰倒角,尺寸完全一样:前者删掉会作废整个评定点(1.11),后者删掉一点代价都没有。尺寸这条线索天然分不开它们——而自动清理工具用的恰恰是尺寸。

工程量级。删简对整体刚度、位移、模态频率的影响通常在百分之几量级,对局部应力是数量级差别(下一层展开在 1.12)。

易错点。① 用「小于某个尺寸就删」的自动清理批量跑一遍——它会连承载 R 角一起删掉,而清理日志里那一条与 logo 那一条长得一模一样;② 只看几何大小,不看它属于哪个响应量:同一个零件做模态时可以大删,做应力评定时几乎不能删;③ 删完不记账——删了什么必须成表,否则第 7 讲发现某点不收敛时,你无从回溯它是哪一次删简造成的,只能整套重做。

今天就能做的一件事。把自动清理那一步的日志导出来,逐条打勾:这一条删的东西在不在传力路径上。⚠ 判据不是「删了多少个」,是有没有一条你说不出它属于哪一类——说不出的那几条,就是这次删简真正的风险。这张勾完的表同时就是 1.10 要的那本账。

1.10 处置有五种,⛔ 不是「删 / 不删」两种——承载路径上的螺栓、焊点、卡扣全落在后三类

是什么。小特征的处置全集是五种:

处置 用在什么上
① 删除 非承载的 logo、工艺凸台、装配标记
② 保留并局部加密 承载路径上的圆角、评定点
等效替代 ⚠ 本课归纳补上 螺纹 → 光杆;圆角 → 倒角;翅片 / 芯体 → 等效体
用连接单元替代 ⚠ 本课归纳补上 螺栓 → 梁 + 分布耦合;焊点 → 专用连接单元
子模型局部恢复 ⚠ 本课归纳补上 全局模型里删掉,关注区在子模型里原样建回来

⚠ 后三种是本课归纳补上的,⛔ 不是本课程体系既有的分类。

为什么必须补后三类。因为承载路径上的螺栓、焊点、卡扣,既不能原样建,也不能删:原样建(实体 + 螺纹 + 接触)会把自由度与接触对数量推到求解不可用;删掉则传力路径直接断了。一个只有「删 / 不删」两档的清单,会把它们逼进其中一档,而两档都是错的——这就是「漏掉的那一类不会以少一条的形式暴露」的典型:清单本身没有空位,它只是把不属于任何一档的东西硬塞进最近的那一档。

工程量级。一个带十几处螺栓与卡扣的支架,若全部原样实体建模并带接触,自由度与接触对数量会把求解时间推到不可用(量级判断,⛔ 不作可承诺值)。

易错点。① 把连接单元当成「简化」而不记它引入的刚度假设——刚性连接式的螺栓替代会在孔周局部增刚,第 3 讲专门讲这一条;② 等效替代之后不回头核对等效前后的传力路径是否一致(等效的是刚度还是载荷分配?两者的下游结论不同);③ 子模型的切割边界落在应力梯度陡的地方。

今天就能做的一件事。把 1.9 那张删简账再加一列:这一条用的是五种里的哪一种。所有落在第 ③④⑤ 类的行,再补一句「它引入了什么假设」——这一列就是评审时最先被问到的那一列,也是你唯一有机会提前写好的那一列。

小特征处置决策图。顶端入口框写「这个小特征怎么办」,下接唯一一道菱形判据,标注为「它在不在关注响应量的传力路径上、是不是评定点」,框旁一行小字写:判据不是它有多小。两条出边分别通向五个处置框,各带文字标签与一行例子:删除,例子是 logo、工艺凸台、装配标记;保留并局部加密,例子是承载 R 角与评定点;等效替代,例子是螺纹换光杆、圆角换倒角、芯体换等效体;用连接单元替代,例子是螺栓换梁加分布耦合、焊点换专用连接单元;子模型局部恢复。后三个框右上角标注为本课归纳补上,并用一条横带把它们框在一起,横带上写:承载路径上的螺栓、焊点、卡扣全落在这三类,既不能原样建也不能删。底部一个记账框写:删了什么必须成表,第 7 讲发现奇异时要靠它回溯。本图为定性示意,不含任何尺寸门槛,不得据图读取任何数值。
图6 这张图要读出三件事:一,处置有五种,⛔ 不是删与不删两种;二,判据只有一条,而且挂在传力路径上,⛔ 不挂在尺寸上;三,后三类正是被整类漏掉的那一支——而螺栓与焊点正好全在里面,它们既不能原样建也不能删。★ 底部那个记账框不是流程装饰:删简账是第 7 讲回溯奇异来源时唯一的凭据。⛔ 本图不含任何尺寸门槛,⛔ 不得据图读取任何数值。

1.11 删圆角的代价:把一个有限的应力集中,换成了一个无界的应力奇异

是什么。两种性质,⛔ 不是同一件事的两个数值:

  • 真实圆角处的应力是有限值:网格加密后,那一点的应力单调趋近一个平台值,可以做收敛判定,可以作评定点(应力集中系数 K_t 的定义在 J5-01 换热器/冷板结构强度与耐压 FEA,本课只引用、不重定义)。
  • 删成零圆角尖角之后,理论解的应力无界:加密只会单调发散、永不收敛。该点的「最大应力」不再是一个物理量——它等于多少,完全取决于你停在哪一档网格。

为什么必须在第 1 讲讲、而不是留到第 7 讲。因为奇异问题绝大多数在前处理阶段就已经造成了。第 7 讲能做的只是识别与善后,而那里四条出路里优先级最高的一条就是「补回真实 R 角」——那等于回到第 1 讲重做一遍,连带整套网格。

工程量级。⚠ 教学假设值:不对应任何平台、任何牌号、任何供应商产品,⛔ 禁止照抄取用(完整的判定程序与双对数拟合在第 7 讲)。同一个模型的两种做法,各配三套系统加密的网格 h = 2 / 1 / 0.5 mm,读那一点的峰值应力 σ_peak:

  • 删了 R 角:σ_peak = 180.0 / 230.4 / 295.0 MPa。相邻比 230.4/180 = 1.2800、295/230.4 = 1.2804——几乎恒定;由首末两点拟合 λ = ln(295/180)/ln(2/0.5) = 0.3564 > 0⛔ 无平台
  • 补回真实 R 角:σ_peak = 168.0 / 176.0 / 178.0 MPa。相邻变化 (176−168)/168 = +4.762% 与 (178−176)/176 = +1.136%——迅速衰减;末段拟合 λ = ln(178/176)/ln2 = 0.0163 ≈ 0出现平台

⚠ 两组数必须成对看:单看任何一条序列,都分不出「发散」与「收敛得慢」。判据是相邻比衰不衰减,⛔ 不是某个百分比阈值。

易错点。① 认为「反正加密到某个尺寸就够了」——奇异点上没有「够了」这回事;② 删了圆角又在那一点报应力,且因为数值很大而判「必裂」——这是虚警,本课典型案例里的第二处错就是它;③ 删了圆角、改用沿壁厚线性化取值,就以为那次删简的代价没了。⛔ 这是本课最容易被读反的一个动作:线性化的作用是把性质不同的应力分量分开(哪一部分与外载平衡、超限即塑性垮塌,哪一部分是局部的、自限的),⛔ 不是把数值降下来。它作废的是这个位置作为评定点的资格(局部应力与疲劳),而一次应力合格并不把这个资格救回来——奇异点上的峰值本来就不是物理量,对一个不是物理量的东西做分解,得到的不会因此变成物理量。正确动作是按第 7 讲那四条出路的优先级走,首选补回真实 R 角

今天就能做的一件事。把你模型里所有被删成尖角的位置标出来,然后对照评定点清单划掉重叠的那些行——划掉的每一行都是一个「报告里还在报数、而那个数不是物理量」的位置。⚠ 这一步的产物不是「哪里应力高」,是一份收缩过的评定点清单;先有它,第 7 讲的收敛检查才知道该在哪几个点上做。

删圆角前后的应力分布对照图。左右两栏。左栏上方画带真实 R 角的几何剖面并标注真实 R 角,下方一条应力沿表面分布的定性曲线,峰值处画一条水平虚线并标注:有限平台值,可作评定点。右栏上方画删成零圆角的尖角剖面并标注零圆角尖角,下方同一位置的分布曲线在尖角处向上不封顶,箭头旁标注:无界,加密即发散,该点不再是评定点。两栏之间一个由左指向右的粗箭头标注删圆角,箭头下方写:代价不是省了网格,是作废了一个评定点。底部一行前指提示:量化对照见第 7 讲的三套网格序列。本图为定性示意,纵轴不设刻度、两条曲线不代表任何平台的实际应力,不得据图读取任何数值。
图7 这张图要读出三件事:一,两条曲线的差别不在高矮,在右边那条没有平台;二,删圆角改变的是这一点的性质,⛔ 不是它的数值大小——所以它不能靠换一种取值方式补救;三,一旦成为奇异点,它就退出评定点清单,而这一步没有任何东西会提醒你。⛔ 本图为定性示意,⛔ 纵轴不设刻度,⛔ 不得据图读取任何数值。

1.12 删简的影响边界必须分响应量说:刚度类是百分之几量级,应力类是数量级

是什么。删掉远小于结构特征尺寸的小特征,影响按响应量分开算:

响应量 删简的影响 结论
整体刚度与位移 百分之几量级 删简自由度大
模态频率 百分之几量级 删简自由度大
沿壁厚的膜应力 视删掉的东西是否承载 要看传力路径
该处的局部应力 数量级差别(极端情形:有限的集中 ⇒ 无界的奇异) 删简自由度小

同一次删简,在模态分析里几乎免费,在应力评定里可能直接作废一个评定点。

为什么这条必须写在删简之前而不是之后。因为删简是不可逆的前处理动作:等到第 7 讲发现某点不收敛再回头补几何,往往整套网格都要重来。而做删简那一刻,人手上通常只有一个模糊的目标(「做个分析」),还没有确定要读哪个量——于是默认按最省事的那一档删。

工程量级。上面两句本身就是量级判断(行业典型观察),⛔ 不是可对外承诺的误差。⚠ 这一点要说清楚:把「影响百分之几」写进报告当成精度保证,是把一句量级观察升格成了一个承诺。

易错点。拿模态模型的删简程度去做应力模型——本课点名的高发形态,因为两个模型往往是同一份几何派生的,而派生时没有人重新问一遍判据;② 反过来,为了应力评定一点都不敢删,网格量爆炸后被迫整体粗化网格——在第 ③ 环上把省下来的还回去了;③ 把「百分之几」当成可承诺的精度对外报出。

今天就能做的一件事。给每一份删简后的几何打上它服务的响应量标签(「模态用」「应力评定用」),并规定标签不同的模型不共用几何。⚠ 判据不是「能不能共用」(多数时候技术上能),是共用之后谁来重新走一遍删简判据——答不上来这个问题,就分成两份。

删简影响按响应量分层的带状图。纵向按响应量分四行,各带文字标签:整体刚度与位移、模态频率、沿壁厚的膜应力、该处的局部应力。前三行的影响带画在一个窄区并标注文字:百分之几量级。第四行的影响带跨越若干个刻度并标注:数量级差别,行尾另标注:极端情形是有限的应力集中变成无界的应力奇异。横轴只标影响量级并以对数刻度示意,不标任何具体百分数。右侧一列结论标签,前三行写刚度类问题的删简自由度大,第四行写应力类问题的删简自由度小。本图为定性示意,横轴为量级示意且无刻度值,不得据图读取任何百分数。
图8 这张图要读出三件事:一,同一次删简在不同响应量上差了若干个量级——所以「这次删简可不可以」这个问题在没有响应量之前无解;二,拿模态模型的删简程度去做应力模型,是本课点名的高发形态;三,第四行的极端情形会直接作废一个评定点(1.11)。⛔ 横轴为量级示意、⛔ 无刻度值,两句量级判断是行业典型观察、⛔ 不是可承诺的精度,⛔ 不得据图读取任何百分数。

1.13 结构侧的删简与 CFD 侧不重讲;而「小特征」这个词在另一门课里是另一个意思

是什么,两件事。

其一,结构侧与 CFD 侧的删简判据来源不同。两边都做「删简」,动作名称一模一样(清理、缝合、简化),但删什么、留什么由各自的响应量决定:结构侧问的是「在不在传力路径上、是不是评定点」,CFD 侧问的是「影不影响流道截面、壁面剪应力与换热面积」。同一个特征在两侧的处置可以完全相反——一个非承载的导流小筋,结构侧可以删,CFD 侧删掉会改变流场分配;一个承载 R 角,结构侧必须留,CFD 侧无所谓。⇒ 脏数模验收、破面缝合、包面、计算域抽取一律去向 J2-06 整车/机舱 CFD 前处理:从脏数模到可算网格(数模验收·删简·包面·计算域抽取),网格无关性口径去向 J2-01 CFD 基础、湍流模型与网格划分(Fluent/STAR-CCM+),本课只讲结构侧。

⚠ 其中最要紧的一条:⛔ 不要拿 CFD 侧的包面 / 水密模型直接当结构模型的输入——包面取的是外表面,壁厚信息已经没有了,而本课几乎每一条判据(跨厚比、线性化路径、Pm 的取法)都要用到壁厚。

其二,「小特征」这个词在拓扑优化里指的是另一件事。本课的「小特征」是已有几何上要不要删掉的局部细节J3-03 冷板/流道拓扑优化 拓扑优化里的「最小特征尺寸」是一条制造约束,限定优化结果里允许出现的最细杆件与最薄壁,防止优化出造不出来的形状。⚠ 两者的动作方向恰好相反:本课是「从已有几何里拿掉」,那边是「不许生成小于某尺寸的东西」。串起来读会得出「优化时已经设了最小特征尺寸,所以结构模型里就没有小特征要清理了」这种完全不成立的结论。

易错点。① 两侧共用一份删简清单而不标注它服务于哪个响应量;② 在同一份报告里两处含义混写而不加限定;③ 把拓扑优化输出的几何直接拿来做应力评定——优化输出的过渡处往往是数值上的尖角,正是 1.11 讲的那类奇异来源。

1.14 几何缺陷的四条后果链:四条里只有一条会自己喊

是什么。四类缺陷,各自在结构侧的后果链:

缺陷 结构侧的后果链 谁会告诉我?
缝隙 接触对与连接建不上 → 网格在此断开 → 传力路径丢失 ⛔ 没有人——只算出「这一支不受力」
自交 / 拓扑错误 网格划分失败,或生成负体积 / 翻转单元 会报(落在第 ⑦ 环,响亮,好处理)
单位制不一致 量纲级错误 → 结果错的是数量级 ⛔ 没有人(下一节专治)
导出精度不足 小圆角处生成锯齿面 → 每个折点都是一个人造奇异点 ⛔ 没有人

为什么把它们放在本讲末尾。因为这四类与本讲前面所有内容属于同一层——输入本身就错了;而四条里只有第二条会自己喊,其余三条落在九环里「没有人会告诉你」的那一档。

工程量级。导出的弦高容差要与关注区的最小圆角半径匹配:容差取到远小于那个 R 角半径的一档。⚠ 具体值按你所用的前处理器与零件定,本课⛔ 不给数;可执行的判据是把容差再收紧一档、重新导一次,看关注区那个圆角上的峰值应力变不变——变了,说明你之前读到的是锯齿而不是几何。

易错点。① 只看「模型打开没报错」就当几何合格;② 缝隙用「自动缝合」一键处理而不看它把哪些面拉到了一起(它可能把两个本该分开的件缝成一体,于是接触界面整个消失);③ 把锯齿面上取到的峰值应力当成真实的应力集中去评定——它是一个人造的奇异点,加密只会让它更高,而你会以为自己找到了一个危险位置。

今天就能做的一件事。几何进网格之前跑一遍三查:查自由边(缝隙)· 查自交与拓扑报错 · 查关注区那个圆角上有几个面片。⚠ 第三条最容易被跳过而它最便宜:一个 R 角上只有两三个面片,说明导出精度不够,那里之后无论加密到什么程度,读到的都是锯齿的峰值。

几何缺陷在结构侧的四条后果链图。四条并列的横向后果链,每条起点是一类缺陷、终点是它在结构侧的后果,链上每一步各带文字标签。第一条:缝隙,到接触对与连接建不上,到网格在此断开,到传力路径丢失,终点框标注求解器不报错并用灰底。第二条:自交与拓扑错误,到网格划分失败或生成负体积与翻转单元,终点框标注求解器会报并用描边突出。第三条:单位制不一致,到量纲级错误,到结果错的是数量级,终点框标注求解器不报错并用灰底。第四条:导出精度不足,到小圆角处生成锯齿面,到每个折点都是人造奇异,终点框标注求解器不报错并用灰底。底部一行写:四条里只有第二条会自己喊。本图为定性示意,不含任何容差数值,不得据图读取任何数值。
图9 这张图要读出三件事:一,四类几何缺陷里只有一类会被求解器抓到,另外三类要靠人主动去核;二,缝隙的后果是「有一支不受力」而不是「算不出来」——它以一个正常的结果出现;三,锯齿面把一个几何精度问题变成了第 7 讲的奇异问题,而在云图上它跟真实的应力集中长得一样。⛔ 本图不含任何弦高容差与缝隙容差数值(那属前处理器与项目规范),⛔ 不得据图读取任何数值。

1.15 单位制一致性自检是量纲级的:mm-t-s-MPa 制下,铝的密度写作 2.7e−9 t/mm³

是什么。结构求解器多数不带单位:数值的含义由你选定的那套一致单位制决定。本课全课取 mm-t-s-MPa 制:

本制单位 与常见写法的关系
长度 mm
质量 t 1 t = 1000 kg
时间 s
N
应力 / 弹性模量 MPa 1 MPa = 1 N/mm²
能量 mJ 1 mJ = 1e−3 J
密度 t/mm³ 铝的 2.7 g/cm³ ⇒ 2.7e−9 t/mm³

⚠ 上表的换算是恒等式,⛔ 不是经验值——每一行你都可以自己验一遍,而这正是它能当自检用的原因。⚠ 但要分清两件事:铝的密度是物性值,可以给;而某个牌号在某个状态下的弹性模量 E 与屈服强度 σ_y 属材料取值口径,⛔ 不在本课(去向 J5-01 换热器/冷板结构强度与耐压 FEA · H1-01 热管理用铝合金牌号、性能与选材 · H1-04 铸造铝合金(压铸为主)的材料特性与热管理壳体选材)——本课只管把已定的属性正确地建进模型。

为什么它的后果最不成比例。因为只要有一个属性沿用了另一套制的数值,错的就是数量级——而云图照常出、收敛照常好、报告照常写得出来。它和几何缺陷属于同一层(都是「输入本身就错了」),却比几何缺陷更难被发现:几何至少还能看,而一个数字看不出它属于哪套制。

易错点。① 只自检长度与应力,忘了密度(影响自重与模态)与线胀系数(影响热应力——而热应力恰恰是本模块读者最常算的那一类);② 从材料库里直接拖一条另一套制的材料进 mm 制模型;③ 把能量单位当成 J 去比第 2 讲那个沙漏能。⚠ 第 ③ 条要说清方向,否则会被读成「所以沙漏能判据也不可信」:能量比 E_hg/E_int 是无量纲的,那个比值不受单位制影响;受影响的是能量的绝对值读数会差 1000 倍——拿绝对值去判断的人会判反,拿比值判断的人不受影响。

⚠ 还有一句关于危险分布的:真正危险的不是那些差 10⁹、10⁶ 的整数量级偏差——那些一眼可辨(如果你去看的话)。危险的是只差一两个量级、看起来仍然像一个工程数的那些组合:它们不会触发任何人的怀疑。

今天就能做的一件事。求解之前做一次「一物一验」:挑三个不是从模板继承的输入(比如密度、线胀系数、某个集中力),各自亲手算一遍它在本制下应该是多少,再跟模型里的值比。⚠ 判据不是「看着像不像」,是亲手换算一遍;这三个数对了,其余从同一套模板来的通常也对——而这三个恰恰是模板保护不到的那几个。

后面还有 7 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做

会员专属

后续为会员深水区内容——四库数据与深度拆解。

查看会员方案