J5-07 · 仿真分析 CAE / 结构/热应力/疲劳 FEA
待认领listed
→
众筹中recruiting
→
已开课open · Course 售卖
一课题可并存多讲师版本
课程大纲
学完能做什么
- 能按板厚、接头形式与可用曲线三条判据,判定某接头该走名义/热点/缺口哪一级应力口径,并说清「口径与 S-N 必须成对」在两个方向上分别错成什么样
- 能对 mm 级板厚焊接接头完成 0.4t/1.0t 表面外插取得热点应力,并用三问判据判定该接头是否落在 FAT 等级体系的适用边界之内(薄壁异质钎缝一律判否)
- 能用节点力法沿焊线算出膜+弯分解与等效结构应力,并用三套网格实测验证其网格不敏感性,说明何时该方法只能作排序量而非寿命输入
- 能建立合法的子模型——切割边界位置、位移驱动、母模型网格充分性检查——并指出子模型结果不可信的三种情形
- 能用台架试件数据回归自建接头 S-N,给出目标存活率与置信水平下的设计下限曲线,并据此算出标注了取值寿命的接头/HAZ 折减系数
- 能把接头寿命评价按 Miner 合账并写齐四项声明(应力口径、曲线来源与标定口径、存活率/置信度、应力比与平均应力处理),使结果具备被评审采信的资格
内容大纲
- 三级应力口径的分层与「口径—曲线成对」纪律
- 三级口径不是三个精度档,而是三种分账方式:名义应力只含总体几何的宏观应力、不含焊缝局部效应;热点应力含结构不连续(板边、开孔、搭接错位)引起的放大、但排除焊趾局部缺口效应;缺口应力把焊趾/填角的实际圆角也算进应力。同一接头三种口径算出的应力可以差数倍,而它们各自都是对的——差别在于剩下那部分集中被记进了 S-N 曲线里。
- 本课第一条硬判据:每级口径必须配同口径标定的 S-N 曲线,混用即错,且两个方向的错法后果相反——把热点应力送进名义 S-N,结构不连续被计入两遍(双重计入),寿命偏保守到件永远不合格;把名义应力送进热点/缺口 S-N,集中被漏掉一遍,寿命偏乐观,方向偏危险。
- 口径的选择判据只有两条:该级集中能不能被网格与建模方式分辨、以及有没有同口径的 S-N 曲线可用。没有对应曲线时,再精细的应力也落不成寿命——这条判据直接决定了第 2 讲的 FAT 适用边界与第 5 讲的自建 S-N 是同一个问题的两个出口。
- 报数纪律(贯穿全课的验收项):任何接头疲劳结果必须同时声明四件事——① 应力口径;② S-N 曲线来源及其标定口径;③ 存活率/置信度;④ 平均应力(应力比 R)的处理方式。缺一项,这个寿命数字连被评审判对错的资格都没有,这正是 J5-02、J5-03 两门误区里「寿命数字没有意义」的根源。
- 母材段与接头段必须分段评价:母材段可用母材 S-N(J5-03 已立此纪律),焊缝/钎缝段换接头口径与接头曲线;同一模型上两段用两套曲线,交界处按较严的一侧取,不得为了「统一」把整件套一条曲线。
- 本课接住 J5-01 的一条死链:J5-01「疲劳压力循环工况转入 J5-02/J5-03 校核」实为死链——J5-02 通篇热循环、J5-03 通篇随机振动,两门均无压力循环内容。压力循环下接头的仿真侧评价从本课入,载荷谱侧归 I6-01,试验侧归 L1-03/L1-04。
- 边界:S-N/ε-N/Coffin-Manson/Miner 累积损伤/平均应力修正等疲劳方法本体不在本课,见 J5-02(热循环)、J5-03(随机振动)、I6-01(损伤累积与折减纪律)、G2-07(换热器部件设计侧算例);本课只解决「送进这些方法的应力幅按哪级口径取、配哪条曲线」。
- 边界:应力奇异的一般识别、沿壁厚线性化与 Pm/Pb/F 分解、单元与约束选型见 J5-06;本课从「该点已判为奇异、一次静强度评定已完成」处接手,只做接头处的专用评价路径。
- 热点应力的表面外插取点与 FAT 等级体系的适用边界
- 表面外插的做法:在焊趾前方沿板面按壁厚 t 的固定比例布两点(常用 0.4t 与 1.0t),线性外插回焊趾位置得到热点应力(公式见公式条)。要讲透的是它的性质——取点位置是「避开焊趾局部缺口、保留结构不连续」的工程约定,外插值不是该点应力的真值,所以它只在配套曲线体系内有意义。
- 取点落地的三条约束:外插区网格尺寸必须与 t 挂钩(常取 t 或 t/2 量级)、须取单元表面应力而非跨面节点平均、外插方向沿板面且垂直于焊趾线。网格一变取点位置随之变,同一模型换套网格就换个热点应力值——所以网格规则必须与 S-N 曲线一起冻结、一起写进报告。
- 热点要先分类再取点:焊趾落在板面上的一类按 t 的比例取点;焊趾落在板边/端部的一类须按绝对距离取点。把后者按 t 取点会取到一个物理上不存在的位置,算出来的数没有对应曲线可配。
- FAT 等级体系是什么:按接头几何形式给出一族参考疲劳强度等级(以指定循环次数下的应力幅表征),用固定斜率延伸成设计曲线,并绑定一个存活率约定。它的价值就一句——不做试验也能有一条曲线可用。
- 适用边界(本课写死,不许含糊):对 FSW/激光焊冷板、钣金支架这类 mm 级板厚的焊接接头,0.4t/1.0t 取点落得进去、接头形式能在分级表里找到对位,FAT 分级可借用(工艺侧与缺陷谱见 H3-01);对 0.2~0.5 mm 壁厚、钎料与母材异质的炉中钎焊铝换热器接头,0.4t 只有零点几毫米,取点落进焊趾影响区内部、外插失去意义,且分级表里没有铝钎焊接头的标定依据。
- 由此定本课的取用纪律:借三级口径的分层框架,不借 FAT 数表。薄壁钎焊接头一律走第 3 讲的等效结构应力法配第 5 讲的自建 S-N;本课不给任何 FAT 等级数值。
- 「能不能用 FAT」的三问判据,逐条答是才可引:① 板厚是否 mm 级、外插取点是否落在母材段;② 接头形式能否在分级表里找到几何对位(而不是「差不多像」);③ 被评价件的母材/焊材族与该分级的标定试件族是否一致、是否同质。三问有一问否,引 FAT 数表等于拿别人的试件替自己的件背书。
- 标准归属只作指路:三级应力口径与 FAT 分级出自国际焊接学会(IIW)焊接接头疲劳评定推荐做法一路文件,具体编号与现行版本按发布方现行目录确认;本课只借其分层框架。
- 网格不敏感的等效结构应力法:节点力反算膜+弯分解
- 方法本体:沿焊线在焊趾截面上取出单元节点力与节点力矩,按截面平衡反算该截面的膜应力与弯曲应力,得到结构应力 σ_s = σ_m + σ_b。节点力是单元积分出的合力、天然满足平衡,不随网格加密而漂移——这就是「网格不敏感」的来源,也是它与「取节点应力再外插」的根本区别。
- 为什么它是奇异位置唯一能落地的定量出口:焊趾是应力奇异点,节点应力永不收敛(识别程序见 J5-06),但截面合力收敛。J5-06 那条「奇异处改用结构应力/网格不敏感法」的指路牌所指的方法本体就在本讲。
- 等效结构应力参数在做什么:把结构应力按壁厚与载荷弯曲比(膜/弯占比)作无量纲化,使不同板厚、不同弯曲比的试验数据能归并到一条主曲线上。工程价值是一条主曲线覆盖多种接头;代价是必须接受它内含的裂纹扩展假设——这个代价要写进报告,不能当成免费的。
- 落地纪律四条:焊线必须在网格里显式定义成一条节点路径(不能只挑一个最大点报数);沿焊线逐节点算、沿线取包络并给出分布;单元排布须垂直焊线(斜排会把切向力算进法向);壳单元与实体单元的节点力提取方式不同,实体须沿壁厚成组求合力。
- 失效边界(比方法更重要):它评的是焊趾/焊根处贯穿壁厚型裂纹扩展主导的失效。对钎缝这类钎料层本身是异质薄层、失效可能是界面剥离或钎料层内部开裂的情形,「把截面当均质壁厚做膜+弯分解」这个前提不成立——此时 σ_s 只能作方案间的排序量,绝对寿命必须靠第 5 讲的自建 S-N 回归到自己的试件上。
- 与静强度侧线性化的分工要讲清:Pm/Pb 线性化(J5-06)与结构应力都在做膜+弯分解,区别是一个用路径上的应力积分、一个用截面节点力;应力积分在奇异点附近随网格漂移,节点力不漂移。这就是「一次静强度评定用线性化、接头疲劳评定用节点力法」的分工依据,不是两派习惯之争。
- 报数要求:交付 σ_s 沿焊线的分布曲线、包络位置、膜/弯占比,以及三套网格下 σ_s 的偏差——网格不敏感性必须被验证过,不能被假定;三套网格差得远时先查焊线定义与单元排布,而不是去调 S-N。
- 子模型与局部加密:把焊趾级网格合法地放进大模型
- 为什么必须用:整芯体/整总成模型不可能在每条焊线上都布焊趾级网格(芯体本就按等效简化建模,口径见 J5-01)。做法是粗网格全局模型先跑出位移场,再在关注接头处切出细网格子模型,用全局位移驱动其切割边界。
- 切割边界三条判据:① 边界须离开应力集中区足够远——按「离焊趾若干倍壁厚以上、且落在全局模型应力已回到名义水平的区域」取,用全局应力云图判而不是凭经验定;② 边界不得跨过接触面、间隙或材料界面;③ 子模型内不得改变刚度,加圆角、改壁厚、换连接方式后它已不是同一结构,位移驱动随即失效。
- 位移驱动而非力驱动的理论依据:位移是分层收敛里收敛最快的量(J5-06),所以粗网格的位移可以信、粗网格的应力不可以信。反过来说,若全局模型的局部刚度本身错了(单元偏刚、连接用 RBE2 憋死了热膨胀),子模型算得再细也是把错误高分辨率地复制一遍。
- 母模型网格充分性判据:切割边界附近要有足够单元描述位移梯度。可操作的检查是用两套疏密不同的全局网格分别驱动同一个子模型、比子模型结果差异;差异大说明母模型不够,该加密的是母模型而不是子模型。
- 不切子模型、直接局部加密的替代路线及其代价:省一次数据传递,但会引入过渡带的畸变单元。判据是加密比按 2 倍量级逐级过渡,且过渡带不得落在评定点或焊线上。
- 与第 2 讲的耦合约束(最容易被当成优点犯的错):一旦采用 0.4t/1.0t 外插,子模型网格尺寸就必须按 t 规则设定,不能「能加多密就加多密」——超出取点规则的精细化在这里是违规,不是精度提升。走第 3 讲节点力法时无此约束,但仍受单元排布垂直焊线的约束。
- 术语防混:本课的子模型是 FEA 空间子模型;J1-01/J1-04 讲的「子模型」是一维系统模型里的子系统模块,J4-01 的「网格加密」是 CFD 语境,三者同形异义,检索时不要互相当证据。
- 钎缝 S-N 数据缺失时的工程处置:自建接头 S-N 与折减系数定法
- 起点判据:查不到同口径曲线就必须自建。铝钎焊接头没有可引的 FAT 数表(第 2 讲已判死),I6-01 只立「缺口、钎焊/焊接热影响区、残余应力必须折减」这条纪律而不给数——那个数只能自己标,做法就是本讲。
- 试件与载荷设计:试件必须带真实工艺态(同炉钎焊、同钎料与包覆层、同装配间隙与填角形貌),加载方式的膜/弯占比要与服役载荷一致,否则第 3 讲的弯曲比无量纲化无从对应。内压脉动型接头的数据从 L1-03 的耐压/压力循环台架取,振动与热循环型从 L1-04 取;试验执行、当量循环数与 DVP 编排归那两门,本课只提「要什么试件、要多少件、要哪几个应力级」。
- 应力比必须随曲线一起声明:内压脉动是 R≈0 的脉动载荷(平均应力≈应力幅,J5-01 已把「直接套对称循环 S-N」列为实错)。自建曲线要写明标定时的 R;用到不同 R 的工况上才需要平均应力修正,修正式本体在 I6-01。拿 R≈0 的数据当对称循环曲线用,系统性偏乐观。
- 数据处理:在双对数坐标下按 Basquin 形式回归(以应力幅为自变量还是以寿命为自变量结果不同,须声明选哪种),交付斜率、截距与残差分散度;分散度用对数寿命标准差表征。钎焊/焊接接头的分散通常显著大于母材,这一点直接决定下一条的下移量。
- 存活率与置信度下限曲线(本讲核心交付):设计不能用中值曲线(存活率 50%),须按目标存活率与置信水平取单侧容许限,把中值曲线在寿命轴上下移得到设计曲线。样本量少时下移量急剧变大——这就是「多做几个试件比多算几遍模型划算」的定量理由。样本量由目标置信度与可接受的下移量反推,不按凑整数定。
- 批间分散必须单独立项:钎焊接头寿命分散有相当份额来自批次(炉温均匀性、装载方式、钎剂涂覆量,机理见 H2-04 与 H2-05),单批试件拟合出的标准差会低估真实分散、设计曲线随之偏危险。判据是至少跨两个生产批取样,并把批间差作为独立分散项,不并进随机误差。
- 标定样本的合格性前提:本课假定标定试件与被评价件都是合格件——虚焊、漏焊、母材溶蚀这类工艺缺陷件不得进 S-N 样本(缺陷判别与形貌反查见 H2-05,检测与判读见 H2-06)。把带缺陷件混进样本,会同时抬高分散度并压低曲线截距,等于把一个可整改的工艺问题标定成不可改的材料属性。
- 接头折减系数怎么落成一个能填进计算书的数:在同一目标寿命处,取母材同口径 S-N 的应力幅比接头设计曲线的应力幅(公式见公式条)。这正是 I6-01 那条纪律缺的定法。关键限制——两条曲线斜率不同,折减系数随目标寿命变化,报数必须标注是在哪个寿命量级下取的;不标寿命的折减系数是没有定义的。
- 曲线作为资产的纪律:自建曲线须连同试件图、工艺态记录、加载方式与 R、拟合参数、分散度、存活率/置信度与样本量一并登记;工艺态一变(换钎料、改间隙、动炉温窗口)该曲线即失效、须重标,不能因为「上个项目标过」就继续引用。
- 边界:熔焊接头的静强度折减取值与 HAZ 软化带位置见 H3-06,本课只用其结论作输入;本讲折减的是疲劳曲线,两者不可互相代用。
关键公式
σ_hs = 1.67·σ_0.4t − 0.67·σ_1.0t
把焊趾处不可评定的奇异应力换成一个可以去配热点 S-N 的量,是三级口径中间那一级的唯一落地取法
σ_s = σ_m + σ_b,其中 σ_m = F_n/t、σ_b = 6·M/t²
在奇异位置给出一个网格不敏感的定量应力,是 J5-06 那条「奇异处改用结构应力/网格不敏感法」指路牌所指的方法本体
lg N_design = (a − m·lg Δσ) − k(n,P,C)·s_lgN
把台架寿命数据变成能用于设计的下限曲线——设计不许用中值曲线,这一步也是接头折减系数的唯一定法出口
k_joint = Δσ_母材(N_target) / Δσ_接头设计曲线(N_target)
把 I6-01「缺口、钎焊/焊接热影响区必须折减」这条只立纪律不给数的要求,落成一个能填进耐久计算书的数
关键概念
名义应力热点应力缺口应力三级应力口径口径—曲线配套(成对不可混用)表面外插(0.4t/1.0t)热点类型(板面型/板边型)FAT 等级结构应力节点力法膜应力+弯曲应力分解等效结构应力载荷弯曲比主 S-N 曲线网格不敏感性FEA 子模型切割边界位移驱动母模型网格充分性Basquin 回归对数寿命标准差存活率置信度下限曲线(单侧容许限)中值曲线批间分散接头折减系数应力比 R焊趾钎缝填角接头 S-N
推荐工具与标准
Ansys Mechanical(节点力提取、子模型 submodeling、线性化后处理) Abaqus/Standard(节点力输出、submodeling 的位移驱动实现、接触与非线性) 求解器内的网格不敏感结构应力模块(如 Verity 一路的商用实现;名称、可用性与无量纲化参数口径按所用求解器版本确认) nCode DesignLife / FEMFAT(接头疲劳后处理、焊缝求解器与自定义 S-N 曲线导入;本课只用其自定义曲线通道,不使用其内置 FAT 数表) HyperMesh / ANSA(焊线定义、沿焊线垂直单元排布、外插区按 t 规则布网、子模型切割) Python(numpy/scipy)或 MATLAB、nCode GlyphWorks(Basquin 回归、对数寿命标准差、单侧容许限与设计下限曲线计算) L1-03 压力循环/耐压台架与 L1-04 振动/热冲击台架(自建接头 S-N 的试件数据来源;试验执行归那两门) 金相/断口检测能力(用于确认标定试件为合格件、核对失效位置与预测位置是否一致;方法归 H2-05/H2-06)
焊接结构疲劳评定的三级应力口径(名义/热点/缺口)、表面外插取点规则与 FAT 等级体系:出自国际焊接学会(IIW)焊接接头与构件疲劳设计推荐做法一路文件,具体编号与现行版本按发布方现行目录确认。本课只借其口径分层框架,不引用其 FAT 数表(薄壁异质钎缝无标定依据) 等效结构应力(网格不敏感结构应力)法的评定路线在压力容器规范体系的弹性应力分析/疲劳评定一路中已有采纳,具体规范册次与现行版本按现行版本确认;车用薄壁钎焊换热器并不直接适用,本课只借方法本体 疲劳试验数据的统计处理与存活率/置信度下限曲线的取法:金属材料疲劳试验数据统计分析方法一类标准,发布方与现行版本按现行目录确认;本课不给容许限系数数表 试验侧依据随所选台架体系确定,见 L1-03(耐压/压力循环/爆破/检漏)与 L1-04(振动/热冲击/温度循环/综合耐久)所列标准体系,本课不另立 各主机厂/Tier1 的接头疲劳评价与仿真报告企业规范(企业 BKM,无统一编号;许用应力、目标存活率与折减系数随项目与产品等级浮动,按现行客户规范确认) 求解器侧的节点力输出、结构应力模块与子模型实现口径各家不同(焊线定义方式、沿壁厚成组规则、无量纲化参数取法均有差异),判据必须按所用软件版本的理论手册确认
工程案例
某电池冷板集流管—扁管接头在压力脉冲台架跑到数十万循环量级出现渗漏,而交付的仿真报告按母材 S-N 给出「寿命大于 1e7 循环、裕度充足」(本例数字均为演算假设值)。复盘查出三处方法错叠在一起:① 应力取的是焊趾节点最大应力,却去配母材 S-N——口径与曲线双错;② 网格从 1 mm 加密到 0.3 mm 后,同一点应力由约 180 升到约 300 MPa 且无平台,预测寿命从 1e7 掉到 3e5,证明用的一直是奇异点应力而不是物理量;③ 内压脉动是 R≈0 的脉动载荷,却直接套了对称循环曲线,未做平均应力处理。整改按本课链条顺序做:先按 J5-06 的 σ_peak~h 判据确认焊趾为奇异点,改用节点力法沿焊线算等效结构应力(包络处 σ_s≈200 MPa、弯曲比约 0.5);再用同炉钎焊、跨两个生产批的 12 件试件在 L1-03 压力循环台架取三个应力级,回归得 m≈4、s_lgN≈0.25,按 P95/C90 取设计下限曲线;最后按 Miner 合账,预测寿命落到 2e5~5e5 循环,与台架实测同量级,并算出 1e6 处 k_joint≈2.5、1e5 处约 2.0 回写进 I6-01 口径的耐久计算书。真正的教训在方向上:三处错的方向并不一致——奇异应力偏高是虚警,母材曲线与对称循环假设各偏乐观一档;两边抵消后报告「看着刚好合格」,这种互相抵消比单向错更危险,因为它连不合理的迹象都不留。
动手做
交付物 · 取一个集流管—扁管接头(或冷板封道焊接头)的模型与一组给定台架数据(12 组寿命点、跨两个生产批、3 个应力级),走完「奇异识别→等效结构应力→自建 S-N→Miner 合账」全链,交付一页《接头疲劳评价卡》,须含六项且逐项可检查:① 用三套网格的峰值应力与 σ_peak~h 拟合结果确认焊趾为应力奇异点(写出拟合指数与判定结论);② 建一个合法子模型——写明切割边界离焊趾多少倍壁厚、依据全局模型哪张云图判定、位移驱动方式,以及用两套疏密不同的全局网格驱动同一子模型的差异检查结论;③ 沿焊线用节点力法给出 σ_m / σ_b / σ_s 的分布与包络值、弯曲比,并附三套网格下 σ_s 的偏差表以证明网格不敏感性(不是声称,是给数);④ 用三问判据回答本接头能否引 FAT 数表(板厚、接头形式对位、母材同质性),并写明因此走哪条路线;⑤ 用给定数据回归自建 S-N,报 m、s_lgN、批间差是否单列,给出 P95/C90 设计下限曲线,并算出 1e5 与 1e6 两个寿命点的 k_joint(各自标注取值寿命);⑥ 按 Miner 合账给出寿命结论,并在卡片末尾写齐四项声明(应力口径/曲线来源与标定口径/存活率与置信度/应力比 R 与平均应力处理)。交付物不得包含:任何 FAT 等级数值的引用(本接头已判不适用)、以中值曲线作判据的结论、未声明口径的应力数字,以及安全系数达标结论(限值归 J5-01 与客户规范)。
常见误区
- 以为应力算得越精细越准,其实焊趾处网格越细峰值越高且永不收敛,精细化只是在放大一个非物理量——奇异点上唯一的定量出口是节点力法的结构应力(奇异识别本体见 J5-06)
- 以为热点应力就是「更准的应力」,其实它是带取点约定的工程量:换取点方案、换网格尺寸就换一个数,必须与配套 S-N 和网格规则一起冻结才有意义
- 以为 FAT 等级表是通用查表工具,其实它按接头形式分级、标定在 mm 级板厚的同质焊接接头上;0.2~0.5 mm 壁厚的异质钎缝在表里没有对位等级、没有标定依据,照查等于拿别人的试件替自己的件背书
- 以为把焊趾应力配母材 S-N「再加个安全系数」就稳了,其实口径与曲线不成对时误差是数量级级的,安全系数那点倍数补不上——这正是 J5-02、J5-03 已点名的断点
- 以为口径混用只会偏保守,其实两个方向后果相反:热点应力送名义曲线是双重计入(偏保守到件永不合格),名义应力送热点曲线是漏计集中(偏乐观,偏危险)
- 以为子模型越细越好,其实其合法性来自母模型位移的收敛与切割边界的位置;母模型刚度错了或边界切在应力集中区内,细网格只是把错误高分辨率地复制一遍
- 以为切割边界「离焊缝远一点就行」,其实判据要落到全局应力云图上——边界必须落在应力已回到名义水平的区域,且不得跨过接触面、间隙或材料界面
- 以为拟合出一条 S-N 就能用,其实设计必须用下限曲线;直接用中值曲线等于按 50% 存活率设计,目标寿命前有一半的件就该坏
- 以为多算几遍模型能省试验,其实小样本时置信下限的下移量急剧增大——寿命裕度被样本量吃掉的份额,远大于被模型精度吃掉的份额
- 以为一个接头折减系数可以全寿命区通用,其实母材与接头两条曲线斜率不同,折减系数随目标寿命变化;不标注取值寿命的折减系数是没有定义的
- 以为拟合出的分散度就是材料分散,其实钎焊接头的分散有相当份额来自批次(炉温均匀性、装载、钎剂涂覆量,见 H2-04/H2-05);单批取样低估分散,设计曲线随之偏危险
- 以为把带缺陷试件一起拟合「更贴近实物」,其实那会同时抬高分散度并压低曲线截距,等于把一个可整改的工艺问题标定成不可改的材料属性(缺陷判别见 H2-05)
- 以为内压脉动可以直接套对称循环 S-N,其实它是 R≈0 的脉动载荷、平均应力≈应力幅,不做处理会系统性偏乐观(J5-01 已把这条列为实错)
- 以为等效结构应力法对所有接头通用,其实它评的是贯穿壁厚型裂纹扩展主导的失效;钎料层异质薄层的界面剥离型失效不满足均质截面假设,此时 σ_s 只能作排序量
- 以为寿命算不准就是网格问题,其实四项声明(口径/曲线来源/存活率/平均应力处理)缺一项,这个数字连被评审判对错的资格都没有
相关课题