J5-08 非金属件的结构仿真：各向异性本构、蠕变松弛、超弹标定与密封接触

项	内容
课程代码	J5-08
所属板块	J 仿真分析 CAE ／ J5 结构、热应力与疲劳 FEA
难度	L3 专家
时长	4.0~4.5 小时（含 1 次全链实操）
适合对象	结构／仿真（CAE）工程师中做非金属件与密封件仿真的人——把 J5-01 换热器/冷板结构强度与耐压 FEA（耐压）、J5-02 热应力、热变形与热疲劳仿真（热循环）做完，却仍回答不出「这个塑料件三年后还剩多少预紧力、这条密封线寿命末期还封不封得住」的，本课是那条断链的接口；也适合密封／塑料件设计工程师与评审人，用来判断一份非金属件仿真报告的材料卡有没有资格被采信。⛔ 不面向只做金属件线弹性静强度校核的岗位。
前置	J5-06 结构 FEA 建模方法基础：单元、约束、本构与结果可信度 结构 FEA 建模方法基础：单元、约束、本构与结果可信度。⚠ 另有四项本课程体系之外、需自备的基础：① 连续介质力学与非线性有限元的概念层认识（大变形、几何／材料非线性、接触求解）；② 高分子材料黏弹性的概念层认识（蠕变与应力松弛的区分、玻璃化转变、时温等效）；③ 曲线拟合与最小二乘的概念层认识（第 2 讲拼主曲线、第 3 讲多组曲线联合拟合直接要用）；④ 至少一次用 Abaqus／ANSYS 一类求解器跑通过含接触的非线性算例——本课讲标定、映射与判据，⛔ 不讲软件菜单路径。
本讲义定位	本课交付的是非金属／弹性体件的「本构—标定—判据」链：怎么把材料的真实行为变成一张能被评审采信的材料卡，怎么把它算到寿命末期，怎么把结果收成一个能判合格与否的数。⛔ 它不交付任何静强度或耐压的达标结论（那是 J5-01 换热器/冷板结构强度与耐压 FEA），⛔ 不交付密封结构的设计准则（那是 G11-03 O 型圈、密封面与静/动密封设计／I3-02 静密封/动密封结构设计与选型），⛔ 不交付有限元建模的通用方法（那是 J5-06 结构 FEA 建模方法基础：单元、约束、本构与结果可信度），⛔ 不交付寿命评价（那是 J5-02 热应力、热变形与热疲劳仿真／J5-03 振动、模态与随机振动疲劳／J5-07 钎缝/焊缝接头的疲劳评价方法：三级应力口径、等效结构应力法与自建 S-N）。

引言

有一类仿真报告，每一页都是绿的，而件在台架上漏了。

它通常长这样：塑料水室—铝集流管的卡边密封，材料卡是「PA66-GF30，E = 某个数，ν = 0.35」，本构线弹性，几何取名义 CAD，载荷加内压，结果是接触比压平均值远高于门槛，结论「合格」。热循环耐久跑到后期，沿周某几个爪位开始渗漏。复盘的时候会发现，这份报告里没有一处是「算错了」——每一步在它自己的假设下都成立。问题在于假设本身：GF 增强件不是各向同性的、橡胶垫圈的反力不是不随时间变的、卡边件的初始状态不是无应力的。

这三处假设有一个共同点，也是本课最重要的一句话：它们的误差方向一致，全部偏乐观。等效各向同性抹平了法兰面沿周刚度的弱点、不建松弛等于假设密封力永不衰减、从名义 CAD 起算等于假设回弹后的压缩量就是设计值。三者叠加之后是单向偏危险——不像 J5-06 结构 FEA 建模方法基础：单元、约束、本构与结果可信度／J5-07 钎缝/焊缝接头的疲劳评价方法：三级应力口径、等效结构应力法与自建 S-N 里常见的那种正负相消、最后留下一个「看着刚好合格」的可疑迹象，这里连一处会给出预警的地方都没有。所以这门课的验收口不在结果上，在材料卡的溯源声明上。

这门课在整条链上的位置

金属件的结构仿真有一条成熟的路：材料卡是各向同性线弹性、常数从手册取、应力线性化后与许用值比。非金属件把这条路的每一环都拆掉了：

• 材料卡不是常数表——增强塑料的模量随方向变，而方向由注塑工艺决定，不由设计者决定；
• 材料卡不是时间无关的——塑料与弹性体在恒定变形下反力会一直衰减，寿命末期与装配态是两个不同的力学状态；
• 材料卡不是随手能查的——超弹参数只能由试验数据拟，而且必须是三类试验，不是一类；
• 判据不是应力——密封失效的判据是接触压力沿密封线的分布最低点，它甚至不是一个标量。

本课的五块内容正好对应这四处断裂，第六讲把它们收成一份可被评审采信的报告。

一句话的边界

本课只在「本构类别已经选定」之后接手。「什么条件下必须选蠕变／超弹／各向异性本构」这条选用条件由 J5-06 结构 FEA 建模方法基础：单元、约束、本构与结果可信度 给出，本课不重讲；本课交付的是选定之后的标定、映射与判据。这条界线在后面每一讲都会再点一次，因为它是最容易被越过的一条——一份报告写「我们用了 Ogden 本构」，评审真正要问的从来不是「为什么用 Ogden」，而是「这组参数是拿什么数据拟的、适用到多大应变」。

---

第 1 讲　增强塑料的各向异性本构与取向场映射

1.1 入口在哪里：材料参数由谁给，本课补哪一步

增强塑料件的仿真，材料参数侧的入口是给定的，⛔ 本课不重选：牌号、玻纤含量（GF30／GF35／GF50 这一族）、等时应力—应变曲线、干态（DAM）与调湿态两套数据，全部归 H4-02 塑料件（尼龙/PPS）注塑与耐冷媒相容。那门课已经把话说到位：选完牌号不等于已经拿到长期预紧力预测。本课补的正是从「有牌号数据」到「有一张能算的材料卡」这一步。

⚠ 请注意这一步的性质：它不是查表，是一次计算。手册上的模量是一个（或两个）数，而件上每一个单元的模量都不同，因为纤维在每个位置的朝向不同。这一节的全部内容，就是把「纤维朝哪儿排」这件工艺事实变成结构模型里逐单元的材料常数。

1.2 为什么等效各向同性材料卡不成立，而且加安全系数补不上

先把最常见的那条捷径钉死。工程上很自然的做法是：从手册取一个模量（通常是流向值，或者流向与横流向的某种平均），配上泊松比，当各向同性材料卡用，然后在结论上加一个安全系数兜底。

这条路结构上就走不通，理由不是「精度差一点」：

GF 增强件的流向与横流向在四个量上都不同——弹性模量、强度、线胀系数、成型收缩率。取平均值之后会同时发生两件事：

• 高估横流向承载——横流向的模量与强度都低于平均值，模型给的刚度和许用都偏高；
• 低估纵向刚度——流向的模量高于平均值，模型给的刚度偏低。

两个方向的错法方向相反。这就是安全系数补不上的原因：安全系数是一个标量倍数，它只能整体地把结论往一个方向挪；而这里的误差是方向相关的、正负各半的，任何一个标量都无法同时修正两个符号相反的偏差。你把系数放大到能盖住横流向那一侧的偏差时，纵向那一侧已经保守到设计做不出来了。

⚠ 划一条界：「GF 增强件必须分方向取值」这条设计判据由 G6-04 冷却风扇总成的结构可靠性、耐久与失效分析 交付（连同随之而来的选材与浇口决策），本课交付的是把它算出来的方法。两者不重叠：那门课回答「要不要分方向」，本课回答「分方向的那组常数从哪儿来」。

1.3 映射链分三段，报告必须分段声明

从「注塑工艺」到「结构模型里的材料卡」，中间有三段，各有独立的误差源：

段	输入 → 输出	典型工具	这一段引入什么误差
① 充填分析	工艺条件 → 纤维取向张量场（在模流网格上）	Moldflow／Moldex3D 一类	取向预测模型本身的误差；浇口位置与工艺参数若与实际不符，整个场都错
② 均匀化	取向张量 + 纤维／基体性能 → 逐单元各向异性弹性常数	Digimat 一类微观力学均匀化工具	闭合近似的误差（见 1.4）；纤维长度分布与界面结合的假设
③ 网格映射	模流网格上的材料常数 → 结构网格上的材料常数	均匀化工具或前处理器的映射模块	空间插值误差；源网格覆盖不全的区域怎么处理

这三段各自独立、误差不共享。因此一条硬要求：报告要分段声明——写明每一段用了什么工具、什么设置、做过什么检查。⛔ 只报最后一个数（「我们做了各向异性映射」）等于没做，因为读报告的人无法判断这个数值得几分信任。

⚠ 这一条不是格式洁癖。三段里任何一段出错，最终材料卡都会以「看起来完全正常的数」的形式呈现——各向异性材料卡不会因为映射错了就报错，它只是给出一组错的常数，然后一路算到底。这条链上没有任何一处会自己说话。

1.4 取向张量的信息损失：闭合近似不是可以略过的实现细节

这一小节是本讲的技术核心，也是最容易被当成「软件内部的事」跳过的一处。

纤维在某一点的排布，完整描述是一个取向分布函数 ψ(p)——p 是单根纤维方向的单位矢量，ψ 给出各个方向上的份额。工程上不直接用它，因为它是一个连续分布，存不下也算不动。于是取它的矩：

a_ij = ∮ p_i·p_j·ψ(p) dΩ，且 a₁₁ + a₂₂ + a₃₃ = 1

这就是二阶纤维取向张量。三个对角元 a₁₁、a₂₂、a₃₃ 是三个坐标方向上的取向份额，无量纲，恒和为 1。它把一个连续分布压成了 6 个独立分量（对称张量），是模流软件实际输出的那个量。

问题在这里：把纤维与基体折算成等效刚度的均匀化算子，数学上需要的是四阶取向张量，而手上只有二阶。这个缺口靠闭合近似补——用二阶量按某种假设「造」出四阶量。常见的闭合近似有若干族，各家软件默认值不同。

⚠ 由此得到一条必须写进报告的判据：取向各向异性越强，闭合近似的误差越大。在取向高度对齐的区域（比如薄壁远端的表层，a₁₁ 接近 1）和取向接近随机的区域（三个对角元接近 1/3），闭合近似的表现不一样；而这两种区域在同一个件上并存。⇒ ⛔ 不能把它当成「求解器内部实现、跟我无关」——它是你选的，而且它改变结果。

一个演算例子（数字均为教学假设值，⛔ 不对应任何牌号或平台）：某注塑法兰面表层单元的模流结果给 a₁₁ = 0.72、a₂₂ = 0.20、a₃₃ = 0.08（三者和为 1）。这说明该处主取向沿流动方向、取向度相当高——流向与横流向的模量差在这里不可能用一张等效各向同性卡蒙混过去。同一件的芯层单元若给 a₁₁ = 0.30、a₂₂ = 0.60，说明该处已经发生取向反转，这正是下一小节要讲的。

1.5 两套网格不重合是常态，插值必须写进报告

模流网格与结构网格几乎从不重合，因为它们是按不同目的划的：

• 模流网格按充填导向——常见形态是中面／双层面网格，或壁厚方向布若干层的体网格；密度取决于要不要分辨熔接线与流动前沿；
• 结构网格按应力导向——在圆角、孔边、接触面加密，在远场放粗。

两者的密度、拓扑、甚至维度都不同。所以第 ③ 段必然是一次空间插值。三件事必须写进报告，缺一项这一段就不可复核：

1. 插值方式——最近点、加权平均、还是别的；
2. 源网格是否覆盖全部目标单元——结构网格上总有一些单元落在模流网格覆盖之外（倒角细节、后加的加强筋、简化差异），这些单元拿到了什么材料常数？软件通常有一个「未覆盖时的缺省」，那个缺省是什么、有多少单元用了它，必须报出来；
3. 映射后有没有做质量抽查——最起码的一次抽查是：在若干代表性位置把结构网格上的材料主方向画出来，看它是不是与流动方向一致、有没有出现明显不合理的跳变。

⚠ 第 2 条最容易漏，而它的后果最难发现：一批单元静默地退回了缺省的各向同性材料卡，而这批单元恰好常常在几何过渡处——也就是应力最高的地方。

1.6 皮层—芯层取向反转不能抹平

注塑件在壁厚方向上不是均质的。熔体在型腔里流动时，靠近型壁的一层受剪切主导，纤维沿流动方向排；中心层受拉伸流主导，纤维转向横流向。结果是壁厚方向上存在一个取向反转的分层结构（工程上常概括为「皮层—芯层」）。

这对建模的直接后果是：用壳单元、只取一个厚度平均取向，会同时算错面内刚度与弯曲刚度。

原因不难看：

• 面内刚度大致由厚度方向上的取向平均决定，所以取平均取向对它的伤害相对小；
• 弯曲刚度由取向的厚度分布加权决定，离中面越远的层权重越大——而离中面最远的正是取向与芯层相反的皮层。把两层平掉，弯曲刚度的方向性会被整体抹掉。

⇒ 判据落在你要读什么上：

你关心的量	建模要求
整体传力、远场面内应力	壳单元 + 厚度平均取向，通常可接受；但须在报告里声明这个简化
弯曲响应（法兰挠曲、盖板变形）	必须按厚度分层建材料（分层壳），或改用实体单元逐层赋材料
密封面翘曲与贴合	同上；且这正是本课第 4 讲那条接触压力分布曲线的上游

⚠ 这条判据与 J5-06 结构 FEA 建模方法基础：单元、约束、本构与结果可信度 的单元选型判据不是同一条，⛔ 别互相顶替：那门课问的是「壳还是实体」（按跨厚比与要读的量），本课加的是一条材料侧的附加条件——即使跨厚比允许用壳，只要关心弯曲或密封面贴合，就得把厚度方向的材料分层建出来。

1.7 与工艺侧、与熔接线的两条界

界一：翘曲成因不在本课。玻纤取向导致收缩不均、进而造成翘曲与密封面平面度超差——这条成因链，连同浇口布置与模温对策，归 G8-03 集成模块的流道注塑/压铸工艺与密封、H4-02 塑料件（尼龙/PPS）注塑与耐冷媒相容 与 G4-03 多通阀（四通/八通）结构与密封设计。本课不重讲成因，只接住「取向场」这一个中间产物往结构侧走。⚠ 请留意这个分工的形状：同一份模流结果，工艺侧拿它去调浇口，结构侧拿它去生成材料卡，两条用途互不替代。

界二：熔接线不在这条映射链里。熔接线（两股熔体前锋汇合处）的强度是折减的，折减系数的取法归 H4-02 塑料件（尼龙/PPS）注塑与耐冷媒相容（旋转件上熔接线该落在哪儿的准则归 G6-04 冷却风扇总成的结构可靠性、耐久与失效分析）。这里有一条必须写死的纪律：

⛔ 映射得到的各向异性弹性常数里，不含熔接线折减。两者必须分开叠加，⛔ 不能指望取向场自动体现。

原因是它们描述的是两件不同的事：取向张量描述纤维朝哪儿排，而熔接线的弱化主要来自汇合面处纤维不跨越界面与分子链缠结不足——这在取向张量上表现为「取向平行于熔接线」，但取向张量不携带界面结合强度的信息。⇒ 一个只做了取向映射的模型，会在熔接线处给出一个看起来正常的刚度，而实际强度显著低于该处的计算值。这一条要在报告里显式写：熔接线位置在哪、折减怎么单独处理、还是判为本次不评价。

本讲小结：增强塑料的材料卡是算出来的不是查出来的。等效各向同性卡的两个方向错法相反，⛔ 安全系数补不上。映射链三段（充填 → 均匀化 → 网格映射），三段误差独立，报告必须分段声明。二阶取向张量只是分布的二阶矩，均匀化实际要四阶量、靠闭合近似补出，取向越对齐误差越大，⛔ 它不是实现细节。两套网格不重合是常态，插值方式、覆盖率、映射后抽查三件缺一不可复核。皮层—芯层取向反转不能抹平，判据是「你要读弯曲吗」而⛔ 不是「壁薄不薄」。翘曲成因与熔接线折减各有属主，⛔ 都不在这条映射链里。