J5-08 · 仿真分析 CAE / 结构/热应力/疲劳 FEA
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 能把模流解出的纤维取向张量场经均匀化映射成结构网格上逐单元的各向异性弹性常数，并分段说明闭合近似、两套网格插值、皮层—芯层分层三处各自引入什么误差、方向偏保守还是偏危险
- 能区分恒载荷下的蠕变与恒压缩量下的应力松弛并选对本构，用多温度点数据拼主曲线与移位因子，给出目标寿命末期的剩余预紧力/密封力，并说清为什么单指数式与 CS 反推系统性偏乐观
- 能按单轴拉伸+平面剪切+等双轴三类试验联合拟合完成超弹标定，写明拟合适用应变区间并做整区间材料稳定性检查，且能当场拒收只用硬度换算或只用单轴数据拟出的材料卡
- 能输出接触压力沿密封线的分布曲线，按最低点判据在冷/热两态与装配态/寿命末期各判一次，并判定何时必须上挤出边界的非线性校核
- 能建立卡边/扣压件的成形＋服役两步分析，把成形残余应力与回弹后的实际几何继承进服役工况，算出寿命末期剩余接触比压并与开漏门槛比
- 能写齐一份非金属件仿真报告必须声明的材料卡溯源项（试验类型与温度点、介质与老化态、拟合区间、取向场来源与映射方式、本构类型及其失效边界），并说明缺项为什么使该报告不具备被采信的资格
内容大纲
- 增强塑料的各向异性本构与模流取向场向结构网格的映射
- 入口料由 H4-02 给定并不再重选：牌号、玻纤含量（GF30/35/50）、等时应力-应变曲线、DAM/调湿态两套数据是材料参数侧的输入；H4-02 §5 已写死「选完牌号不等于已经拿到长期预紧力预测」，本课补的正是这一步
- 等效各向同性材料卡为什么不成立：GF 增强件流向与横流向的模量、强度、CTE 与收缩率均不同，取平均值会同时高估横流向承载、低估纵向刚度，两个方向的错法相反，加安全系数补不上；「必须分方向取值」这条设计判据由 G6-04 交付，本课交付把它算出来的方法
- 映射链分三段、报告要分段声明：模流（Moldflow/Moldex3D）解出纤维取向张量场 → 均匀化算子把纤维与基体折成逐单元各向异性弹性常数 → 按材料主方向映射到结构网格。三段各有独立误差源，只报最后一个数等于没做
- 取向张量的信息损失必须当场讲透：二阶取向张量只保留取向分布的二阶矩，均匀化实际需要四阶量，工程实现靠闭合近似补出；取向各向异性越强，该近似误差越大，不是可以略过的实现细节
- 两套网格不重合是常态：模流网格按充填导向（常为中面/双层面或体网格）、结构网格按应力导向，密度与拓扑都不同，映射须走空间插值；插值方式、源网格是否覆盖全部目标单元、映射后是否做过质量抽查，三项都要写进报告
- 皮层—芯层取向反转不能抹平：注塑件壁厚方向存在取向反向层，壳单元只取一个厚度平均取向会同时算错面内刚度与弯曲刚度；关心弯曲响应或密封面翘曲时须按厚度分层建材料，或改用实体单元
- 与工艺侧划界：玻纤取向致收缩不均→翘曲与密封面平面度的成因、浇口与模温对策归 G8-03、H4-02 与 G4-03 §5.1，本课不重讲成因，只接住「取向场」这个中间产物往结构侧走
- 熔接线不在这条映射链里：熔接线的强度折减系数取法归 H4-02（旋转件上的落位准则归 G6-04）；映射得到的各向异性常数不含熔接线折减，两者必须分开叠加，不能指望取向场自动体现
- 蠕变/应力松弛本构与时温等效：长期预紧力与残余密封力的定量预测
- 先分清在算哪一件事：恒载荷下变形持续增长＝蠕变（塑料螺栓 boss 在恒定预紧下继续变形），恒变形下应力持续衰减＝应力松弛（垫片/O 圈压缩量锁定后反力衰减）。密封力问题绝大多数是后者，拿蠕变曲线回答松弛问题是最常见的错配
- 本构选用判据：小应变、响应仍近似线性的区间用线性黏弹（松弛模量以 Prony 级数表述）；应力或温度高到进入非线性蠕变区时须换蠕变律。判据是应变量级与是否在同一应力水平上验证过，不能靠求解器默认材料卡
- 时温等效是把等不起的时间换成等得起的温度：多温度点短时试验沿 log(t) 用移位因子平移拼成主曲线，再外推到常温长期。WLF 与 Arrhenius 两条移位路线适用区间不同（玻璃化转变附近的物理松弛 vs 远离转变的化学老化），常数必须由本材料多温度点数据拟出，⛔ 不许套用文献「通用值」
- 外推有硬边界：外推区间超出实测时间窗口若干数量级、或跨过玻璃化/结晶/化学老化的机理切换点时，主曲线整体作废——这与 I6-02「加速不能改失效机理」是同一条纪律的仿真侧表述
- 与 I6-02 的分工：I6-02 给试验侧（CS 与 CSR 的口径差、ISO 3384/ASTM D6147 与 ISO 815/ASTM D395 的取数、Arrhenius 加速因子与 Weibull 判读），本课接过它的 CSR 剩余反力曲线当标定数据，输出件级的接触压力场随时间演化；I6-02 那条自标定性的 F_seal(t)≈F₀·exp(−t/τ) 明写「单指数会低估长期衰减」，本讲的多项松弛谱正是它留出的定量替代
- ⛔ 不许拿 CS 反推的剩余密封力当标定数据：G11-03 与 I6-02 都已判定 σ(t)≈σ₀·(1−CS(t)) 系统性乐观（CS 是卸载回弹后的几何量），本构标定只认恒压缩量下连续记录反力的 CSR 曲线
- 螺栓 boss 与金属嵌件周边的预紧力松弛是塑料件的典型场景：H4-02 §4 只到定性，本课把它算成三段账——初始预紧 → 目标寿命末期剩余预紧 → 是否仍高于密封所需夹紧力；金属件高温保载该不该上蠕变本构这条选用条件在 J5-06，本课不重讲选用，只讲标定与外推
- 温度史必须用真实工况谱而不是单一恒温：热循环下松弛与部分回弹交替，按最高温恒定保载算会高估衰减、按平均温算会低估；工况谱的提取与当量归 I6-01，本课只负责把它建成载荷步序列
- 超弹本构的试验标定：三类试验联合拟合与材料卡验收
- 本讲就是 G11-03 讲义里那条被直接指名的流程——它已写死「按 J5-08 的标定流程把试验数据补齐后再算」，交付物是一张可被评审采信的超弹材料卡，而不是一次拟合操作
- 三类试验对应三种独立的变形状态，缺一不可：单轴拉伸、平面剪切（纯剪的工程实现）、等双轴拉伸。只用单轴数据拟出的参数在双轴与受约束压缩状态下外推失真——而 O 圈压在沟槽里正是接近受约束多轴的状态
- ⛔ 只用硬度换算得到的材料卡不可采信：硬度—模量换算只给一个单点线性刚度，拟不出应变硬化段与多轴响应；G11-03 已写明「没有这套材料试验数据时，不要套一个默认本构就当结果可信」
- 本构族取舍是拟合能力与外推稳定性的权衡：Mooney-Rivlin 一类低阶多项式参数少、区外稳；Ogden 一类拟合能力强但高阶参数在数据区外容易给出非物理响应。选型判据是「服役应变范围内拟合好 + 范围外仍单调稳定」，不是比谁的拟合残差小
- 拟合必须限定应变区间并做稳定性检查：把适用应变区间写进材料卡（超出即不适用），并在整个区间上检查参数组是否满足材料稳定性条件——不稳定的参数组会在接触压力分布上给出局部反常峰值，判读时容易被当成真实应力集中
- 长期问题要在超弹之上再配黏弹/松弛参数：G11-03 原话已点明「长期问题还要配黏弹/应力松弛参数」；瞬时超弹拟合只回答装配态接触压力，寿命末期的接触压力必须接第 2 讲的松弛谱
- 介质、温度与老化态是标定的前提条件而不是事后修正项：溶胀后的橡胶与干态不是同一材料；相容与溶胀结论归 H4-01/G11-03，本课只要求标定试样的介质、温度与老化状态必须与被评价工况一致，并逐项写进材料卡
- 不可压缩性怎么建：弹性体近似不可压缩，须用杂交/混合单元一类处置体积锁定；单元与积分方案的一般选型归 J5-06，本课只加一条专用纪律——用普通位移单元算高压缩率密封件会因体积锁定把接触压力算高，方向偏乐观
- 密封接触仿真：沿线分布、最低点判据、冷热两态与挤出边界
- 判据不是平均值：G4-03 §6.5 已写死「内漏发生在分布的最低点，不在平均值上」，P_contact=F_seal/A_contact 只是名义平均值。本课交付的是那条分布曲线本身——沿密封线逐点输出接触压力，取最低点判
- 冷热两态必须分别校核：异种材料 CTE 失配使密封面压紧力随温度变（机理由 G4-03 §5.4 给出），一态合格不等于两态合格；温度场来自热分析，热-结构耦合方式与温度载荷施加归 J5-02，本课只接住温度场做接触求解
- 输出项要成套而不是只有一张云图：接触压力沿线分布曲线、最低点位置与数值、沿线极差、有效接触宽度/面积、穿透量与接触状态、密封线是否出现分离段。接触对设置口径（主从面、初始间隙与过盈、摩擦系数上下界包络、收敛调法）归 J5-06，本课不重讲收敛，只规定输出项与判读
- 挤出边界的触发条件由 G11-03 给出——间隙 × 压力一旦越过挤出边界就必须上非线性（超弹）FEA，线弹性会把大变形下的接触应力与挤出量算得离谱；本课接住这条触发条件，交付的是校核怎么做：大变形、自接触、挡圈接触与低压侧间隙的建模
- 与 I3-04 划死界：I3-04 管设计侧的压缩量沿线分布（紧固点间距、法兰刚度、跨中张开量），本课管仿真侧的接触压力沿密封线分布与最低点判据；同一条密封线两门各给一半，引用时不得互相顶替
- 垫片单元与实体建模的取舍：J5-06 已把垫片单元列入单元适用域，本课给判据——只要关心的是沿线分布的最低点与局部分离，就不能用一维垫片压缩-回弹响应曲线代替实体+超弹/黏弹本构；垫片单元适合的是整体刚度与螺栓载荷分配那一类问题
- 寿命末期态才是真判据：装配态最低点合格只是起点，须把第 2 讲的松弛谱与第 3 讲的超弹卡合起来算到目标寿命末期再判一次——只算装配态是密封仿真最常见的漏项，同口径在 G2-11 的密封力链三段账里已立
- 沿线不均比总量不足更难查：中位值合格而沿线极差超限时，会先在极差最大处渗漏；报数必须给分布与极差，只给一个平均数的报告不可判
- 卡边/扣压件的成形＋服役两步分析
- 为什么必须两步：卡边是一次性塑性成形，成形后件里留着残余应力与回弹后的实际几何；直接从名义 CAD 加服役载荷算，等于假设装配态无应力且几何理想，会把剩余接触比压算高——方向偏危险
- 第一步（成形）要算的是大变形塑性、回弹后的实际压缩量与主片爪/法兰的残余应力状态；欠卷渗漏与过卷开裂的双侧工艺窗口、结构量清单（爪高、爪数、爪距、爪根圆角、主片槽形）归 G2-11，本课不重定窗口，只把它当输入工况
- 第二步（服役）要继承的是上一步的应力/应变状态与变形后网格，再叠内压、热循环与介质老化下的松弛；继承方式（结果传递路径、是否重划网格、状态变量怎么带）必须写进报告——重划网格而不传状态变量，等于把第一步白算
- 与 J5-01 划死界（本课自设的硬边界）：J5-01 管卡边件的几何建模与线弹性本构下的一次耐压静强度，其三个设计输入（卷边几何、装配压缩量、剩余接触比压门槛）由 G2-11 给；本课管两步分析与长期松弛，两门不重叠也不许互相顶替，本课不给任何静强度达标结论
- 判据落在剩余接触比压上而不是应力上：终点是「热循环＋介质老化后剩余接触比压是否仍高于开漏门槛」，门槛值来自 G2-11/H4-01 的设计侧，本课只负责把这个数算出来
- 试验侧对标引 I6-02 §3（压接率、回弹、蠕变松弛）与 L1-03（耐压/泄漏）；仿真给的是上台架前的预测，对标不成立时的排查顺序是先查材料卡（第 2、3 讲）与成形状态继承，再去调接触参数——顺序反了就会用调参把错误的材料卡掩盖掉
- 扣压（软管金属套径向压缩锁死胶管）与卷边（主片爪压紧橡胶垫）机理不同，辨异 H3-02 已讲完、本课不重复；但两者都走同一条「成形＋服役」两步流程，差别只在第一步的载荷与失效判据
- 塑料水室侧的初始状态源不止卡边：注塑残余内应力与翘曲会改变法兰面的初始平面度，是这条链上与卡边并列的第二个初始状态源（成因归 H4-02/G8-03）；报告必须声明有没有把它建进来，没建就要说明这会让最低点结果偏乐观
关键公式
a_ij = ∮ p_i·p_j·ψ(p) dΩ，且 a_11+a_22+a_33 = 1
把「纤维朝哪儿排」变成能进均匀化算子的场量，是模流→结构映射链的第一个中间产物，也是判断「该不该分方向取值」的量化入口
E(t) = E_∞ + Σ_{i=1..N} E_i·exp(−t/τ_i)
把 CSR 实测的剩余反力衰减曲线变成求解器能用的材料卡，用来算「目标寿命末期还剩多少密封力/预紧力」
log a_T = −C₁·(T−T₀)/(C₂+T−T₀)，折合时间 t_ref = t / a_T
把等不起的十年换成等得起的高温短时试验，是长期预紧力与残余密封力预测唯一的取数路径
Mooney-Rivlin：W = C₁₀(Ī₁−3) + C₀₁(Ī₂−3) + (1/D₁)(J−1)²；Ogden：W = Σ_p (μ_p/α_p)·(λ̄₁^{α_p} + λ̄₂^{α_p} + λ̄₃^{α_p} − 3)
弹性体材料卡的拟合对象——密封件装配态接触压力算得对不对，首先取决于这组参数是怎么标出来的
p_c,min(T_state, t) = min_{s∈密封线} p_c(s, T_state, t) ≥ p_system + 裕量，须在 {冷态, 热态} × {t=0, t=t_target} 四种组合上分别成立
把密封仿真那一堆云图收成一个能判合格与否的数，是本课交付物的验收口
关键概念
纤维取向张量闭合近似均匀化（micromechanics homogenization）逐单元各向异性材料卡皮层—芯层取向反转模流网格→结构网格映射蠕变与应力松弛的区分线性黏弹Prony 级数松弛谱松弛主曲线时温等效WLF 移位因子折合时间CSR 剩余反力曲线超弹本构Mooney-RivlinOgden三类试验联合拟合（单轴/平面剪切/等双轴）拟合适用应变区间材料稳定性检查体积锁定与杂交单元接触压力沿密封线分布最低点判据沿线极差冷热两态校核有效接触宽度挤出边界的非线性校核垫片单元适用域成形＋服役两步分析残余应力继承与状态变量传递回弹后实际压缩量剩余接触比压寿命末期态材料卡溯源声明
推荐工具与标准
Moldflow / Moldex3D（充填分析与纤维取向张量场输出，作为映射链的源） Digimat 一类微观力学均匀化与材料卡映射工具（纤维+基体折算成逐单元各向异性弹性常数、向结构网格映射；实现口径与插值方式按所用版本理论手册确认） Abaqus / Standard、ANSYS Mechanical（超弹/黏弹/蠕变本构、接触求解、成形＋服役两步分析与状态量传递） 求解器自带的材料参数拟合模块（超弹多组曲线联合拟合、Prony 级数拟合与稳定性检查） 超弹标定用试验能力：拉力试验机 + 单轴拉伸/平面剪切/等双轴三类试样与夹具（三类缺一则材料卡不完整） 压缩应力松弛 CSR 夹具 + 力传感器 + 恒温老化箱（恒压缩量下连续记录剩余反力；取数口径与标准归 I6-02） DMA 动态力学分析仪（多温度点数据、主曲线拼接与移位因子确定） HyperMesh / ANSA（结构网格与密封线路径定义、沿线结果提取） Python（numpy/scipy）或 MATLAB（主曲线移位与拟合、沿密封线分布与极差的后处理）
橡胶压缩应力松弛（CSR）与压缩永久变形（CS）的试验方法：站内已由 I6-02 给出 ISO 3384 / ASTM D6147（CSR，恒压缩量下测剩余反力）与 ISO 815 / ASTM D395（CS，卸载回弹后的几何量）；本课把 CSR 曲线当作松弛本构的标定输入，CS 数据不得用于标定。版本按现行版本确认 超弹标定所需的单轴拉伸、平面剪切、等双轴试验方法标准：站内尚无可指名来源，⛔ 本大纲不给编号；实施前按现行版本与试验室能力确认，并把所依据的方法与试样尺寸写进材料卡 塑料长期力学性能与设计数据（等时应力-应变曲线一类）的取数口径：站内由 H4-02 的标准条给出，本课不另立，版本按现行版本确认 各主机厂 / Tier1 的 CAE 建模与仿真报告规范（企业 BKM，无统一编号；材料卡溯源要求、许用值与判据随项目与产品等级浮动，按现行客户规范确认） 纤维取向均匀化与取向场映射：无公开统一标准，闭合近似形式、均匀化算子与插值实现各家不同，判据必须按所用软件版本的理论手册确认
工程案例
某塑铝散热器的塑料水室—铝集流管卡边密封，在热循环耐久后期出现沿周局部渗漏，而交付的仿真报告按线弹性＋等效各向同性材料卡算出「装配态接触比压裕量充足、判定合格」（本例数字均为演算假设值）。复盘查出三处方法错叠在一起：① 材料卡用等效各向同性，未做取向场映射——法兰面沿周刚度被抹平，跨中的挠曲与压缩量损失算不出来；② 只算了装配态，且松弛参数是由 CS 数据反推的（σ(t)≈σ₀·(1−CS(t))），系统性乐观，寿命末期那一段账根本没建；③ 从名义 CAD 直接加服役载荷，没做成形分析，回弹后的实际压缩量比名义值小、爪根残余应力也未继承。整改按本课链条顺序做：先按第 1 讲把模流取向张量场经均匀化映射成逐单元各向异性材料卡（分皮层—芯层两层）；再按第 2 讲用 CSR 实测曲线拟 Prony 松弛谱、用多温度点主曲线外推到目标寿命；最后按第 5 讲补成形＋服役两步分析并传递状态变量。重算后装配态最低接触比压 1.8 MPa 仍合格，但寿命末期热态跌到 0.25 MPa、沿线极差 0.9 MPa，最低点落在极差最大的爪位——与台架实际渗漏位置一致，结论由「合格」改为「不合格，须加大爪高或改垫圈截面后重算」。真正的教训在方向上：这三处错的方向一致，全部偏乐观，叠加后是单向偏危险——不像 J5-06/J5-07 案例里正负相消还留个「看着刚好合格」的迹象，这里连一处会给出预警的地方都没有，所以材料卡溯源声明必须逐项检查，不能靠结果看起来合理来兜底。
动手做
交付物 · 取一件塑料水室—铝集流管卡边件（含一段完整密封线、至少一处紧固/卷边爪位、给定模流结果与一组材料试验数据：三类超弹试验曲线 + 三个温度点的 CSR 松弛曲线），走完「取向映射 → 成形残余应力 → 热循环蠕变松弛 → 剩余接触比压是否守住开漏门槛」全链，交付一页《非金属件仿真评价卡》，须含六项且逐项可检查：① 由模流取向张量场经均匀化生成逐单元各向异性材料卡，写明闭合近似的处置、模流网格→结构网格的插值方式、是否按皮层—芯层分层，并给一张映射后材料主方向的抽查图；② 用三个温度点数据拼松弛主曲线、给出移位因子与折合时间，拟出 Prony 参数（报支数、τ 覆盖的时间跨度、E_∞/E(0)），并说明为什么不用 CS 反推、不用单指数式；③ 用单轴+平面剪切+等双轴三类曲线联合拟合超弹参数，报所选本构族与理由、适用应变区间与整区间稳定性检查结论，并另附一次「只用单轴数据拟合」的对照，量化两者在等双轴状态下的偏差；④ 做成形＋服役两步分析，写明第一步的回弹后实际压缩量与爪根残余应力状态、第二步的继承方式（是否重划网格、状态变量怎么传），并给出「不做两步、直接从名义 CAD 起算」的对照结果差多少；⑤ 输出接触压力沿密封线的分布曲线，在 {冷态, 热态} × {装配态, 目标寿命末期} 四种组合上分别给最低点数值与位置、沿线极差，按 p_c,min ≥ p_system+裕量 逐个判定，并指出哪一种组合最先失守；⑥ 在卡片末尾写齐材料卡溯源声明五项——试验类型与温度点、试样的介质与老化态、拟合适用区间、取向场来源与映射方式、本构类型及其失效边界。交付物不得包含：任何静强度或耐压安全系数达标结论（限值与判据归 J5-01 与客户规范）、任何 ESC 的定量判据、任何未声明拟合区间的材料参数、以及只报平均接触压力而不报沿线分布与极差的结论。
常见误区
- 以为 GF 增强件用一张等效各向同性材料卡再加个安全系数就够了——其实流向与横流向差一档，平均值同时高估横流向承载、低估纵向刚度，两个方向的错法相反，安全系数那点倍数补不上
- 以为把模流取向场导进结构模型就算完成映射——其实二阶取向张量要靠闭合近似补出四阶量、两套网格还要空间插值、壳件还要保留皮层—芯层分层，三段各有独立误差；报告不分段声明等于没做
- 以为壳单元取一个厚度平均取向省事又保守——其实皮层与芯层取向相反，平均后面内刚度与弯曲刚度会同时算错，密封面的翘曲与贴合结论可能整体反过来
- 以为蠕变曲线可以拿来回答密封力衰减——其实那是恒载荷下的变形增长，密封问题是恒压缩量下的反力衰减（应力松弛），两者是不同试验、不同本构
- 以为用 CS（压缩永久变形）反推剩余密封力更省事——其实 CS 是卸载回弹后的几何量、系统性乐观，G11-03 与 I6-02 都已判死；本构标定只认 CSR 的剩余反力曲线
- 以为时温等效可以随便往外推——其实主曲线只在同一失效机理内成立，跨过玻璃化、结晶或化学老化的切换点即作废，这与 I6-02「加速不能改失效机理」是同一条纪律
- 以为超弹参数用硬度换算或只用单轴拉伸拟合就能用——其实 O 圈在沟槽内接近受约束多轴状态，单一变形状态拟出的参数在那里外推失真；必须三类试验联合拟合并写明适用应变区间
- 以为 Ogden 阶数越高越准——其实高阶参数在数据区外容易给出非物理响应，选型判据是「服役区间内拟合好、区间外仍单调稳定」，不是比谁的残差小
- 以为拟合完就可以直接用——其实还要在整个服役应变区间上做材料稳定性检查；不稳定的参数组会在接触压力分布上造出局部反常峰值，判读时容易被当成真实应力集中
- 以为密封仿真报一个平均接触压力就够了——其实内漏发生在分布的最低点（G4-03 §6.5 已写死），必须给沿线分布曲线、最低点位置与沿线极差，只给中位值的报告不可判
- 以为装配态接触压力合格就算密封设计做过了——其实还要接松弛谱算到寿命末期，并在冷、热两态各判一次；只算装配态是密封仿真最常见的漏项
- 以为卡边件可以从名义 CAD 直接加服役载荷——其实那等于假设装配态无应力且几何理想，回弹后的实际压缩量比名义小，剩余接触比压会被算高，方向偏危险
- 以为成形分析算完就自动带进服役工况——其实重划网格而不传状态变量等于把第一步白算；继承方式必须显式写进报告
- 以为用普通位移单元算高压缩率弹性体只是慢一点——其实近似不可压缩材料会体积锁定，把接触压力算高，方向偏乐观
- 以为 ESC（环境应力开裂）可以顺手在仿真里给个判据——其实工程上没有普遍接受的可算判据，写成仿真方法只会给错误预期；本课只作边界提示（机理与选材侧见 H4-02，PP 与吹塑件的 ESCR 见 H4-04）
- 以为求解器自带的默认材料卡拿来就能用——其实介质、温度、老化态与拟合区间任一不匹配这张卡就不适用；材料卡溯源声明缺项的报告不具备被评审采信的资格
相关课题