基于DOE敏感性分析的AISI 304亚尺寸Charpy冲击试验Johnson–Cook模型标定

《Metals》：Calibration of the Johnson–Cook Model for Sub-Size Charpy Impact Tests on AISI 304 Steel

　　

准确数值再现Charpy冲击试验需要能够表征材料在动态加载条件下行为的本构模型与损伤模型。本研究针对具有不同厚度与缺口深度的亚尺寸AISI 304不锈钢Charpy试样，开展了冲击响应的试验-数值研究。采用50 J摆锤冲击系统在3.458 m/s冲击速度下开展仪器化试验。试验获得的力、位移、速度与吸收能量历程被用作在LS-DYNA中采用Johnson–Cook本构与损伤模型建立的显式有限元模拟的参考数据。通过改变六个Johnson–Cook参数——A、B、n、D₁、D₂与D₃——相对于其基准值±20%，开展了包含64次模拟的两水平全因子试验设计（Design of Experiments，DOE）。采用一般线性模型（General Linear Model）与方差分析（Analysis of Variance，ANOVA）评估这些参数对峰值冲击力的影响。随后，敏感性结果被用于指导最具影响参数的标定，所得参数集在不同试样几何构型上进行了评估。总体而言，标定后的模型以合理精度再现了初始加载阶段、峰值力量级与能量吸收趋势。最大偏差出现在峰后与断裂阶段，包括试验观测与数值预测断裂条件之间的差异。结果表明，DOE引导的敏感性分析能够缩减参数空间并支持更系统的标定过程。然而，结果亦表明，主要基于峰值力的标定不足以确保对不同亚尺寸试样几何构型峰后行为与断裂条件的准确预测。

**研究背景与问题**　Charpy冲击试验因标准化、快速且成本相对较低而被广泛用于表征金属材料的冲击行为。传统试验结果即吸收冲击能量仅提供响应的全局度量，无法直接描述冲击过程中变形与断裂的演化，且相近的总吸收能量可能来自变形、裂纹萌生与裂纹扩展的不同贡献。仪器化Charpy试验通过记录冲锤力-时间信号可获得力-位移与吸收能量历程，从而考察屈服、最大载荷、裂纹萌生、稳定或不稳定裂纹扩展及最终断裂等特征阶段。然而，断裂阶段的识别依赖所采用的判读准则，并可能受惯性效应、试样振动及冲锤-试样动态相互作用的影响。此外，试样几何形状强烈影响冲击响应：缺口深度改变剩余韧带与局部应力状态，试样尺寸影响约束条件、变形与能量吸收，亚尺寸试样的响应不能假定按尺寸直接缩放或再现标准试样响应。当比较不同厚度与缺口深度试样时，仅凭试验响应可能无法充分解释差异机制，数值模拟因而成为研究变形、损伤发展与断裂演化的有价值补充工具。从试验观测到预测性数值模型的过渡并不简单：动态加载下大塑性变形、应变率敏感性、热效应、损伤累积、接触、惯性与断裂可能同时发生，模型可靠性不仅取决于所选本构与损伤形式，还取决于其参数的准确性与物理一致性。参数可通过直接或反演方法确定；直接法在应力-应变场可视为均匀时较直接，但动态加载下应力波传播、非均匀变形、温度演化与复杂应力状态使这些假设存疑，准静态或简化加载条件下识别的参数向冲击模拟的可移植性可能有限。反演识别将候选参数集引入数值模型并与试验测量比较、迭代调整以最小化差异，但其有效性强烈依赖目标函数选择且可能需要大量计算昂贵的模拟。因此需区分参数识别与模型标定：识别旨在确定描述材料行为内在方面的参数，而标定调整选定参数以改善与特定试验响应的吻合；标定参数集未必唯一，不同参数组合可再现相同全局量却预测不同的局部应力、变形、损伤或断裂历程，且边界条件、接触、离散化或模型形式本身的误差可能通过参数调整被无意补偿。标定目标的选择因而至关重要，仅匹配单一试验量可能对反演问题约束不足，更稳健的程序应考虑多个互补响应特征。敏感性分析是应对该问题的重要预备步骤，可在最终标定前量化单个参数及其可能交互的影响；已有研究表明本构参数影响可能强烈非线性并依赖加载条件与数值输出，识别主导参数可缩减有效标定空间、限制不必要计算量并降低通过物理影响甚微的参数补偿建模误差的风险。试验设计（DOE）提供了开展此类敏感性分析的结构化框架：不同于逐次改变单一参数，DOE按预定义试验矩阵系统组合不同参数水平，可估计各因子个体效应、识别交互并评估统计显著性，同时控制所需模拟数量；在数值应用中每种因子水平组合代表一次模拟而非物理试验，结果采用方差分析（ANOVA）等统计方法分析。DOE获得的信息随后可指导参数标定：对选定响应变量影响可忽略的参数可固定或排除，而最具影响参数的范围可细化以供进一步分析。因此DOE不应被视为标定程序本身，而是理解数值模型并定义更小、条件更好的标定问题的手段；其与反演识别或优化方法结合，为试错调整提供了系统替代方案，并建立模型参数与预测动态响应之间的透明关系。有限元建模（FEM）需表征冲锤、试样与支座之间的完整相互作用以及缺口附近高度局部化的变形与断裂，通常采用显式时间积分格式。材料形式尤为重要，因为预测的力演化与吸收冲击能量取决于应变硬化、应变率敏感性、热效应、损伤累积与断裂的表征方式。Johnson–Cook模型是一种唯象本构形式，用于描述金属材料在大应变、高应变率与高温条件下的塑性行为，考虑应变硬化、应变率敏感性与热软化，并可与损伤形式耦合以描述损伤累积与断裂，因而特别适合Charpy试验等动态冲击问题。已有大量研究采用不同显式有限元代码与本构形式再现Charpy冲击试验，包括采用ANSYS LS-DYNA与Cowper–Symonds模型模拟整体与异种结构钢、采用ANSYS Explicit Dynamics与Johnson–Cook形式研究缺口几何与试样厚度对Q&P钢冲击响应的影响等，表明不同应变率相关本构模型在参数与数值实现适当定义时可再现全局Charpy响应。Johnson–Cook本构与失效形式已被广泛用于此目的，因其对应变硬化、应变率敏感性、热软化、损伤累积与断裂提供相对简单的描述。近期研究进一步强调参数确定与验证分离的重要性：例如确定增材制造1.2709钢的Johnson–Cook流动模型本构参数及其关联失效形式参数，随后通过数值Charpy冲击试验评估其性能；或提出确定API X52管线钢Johnson–Cook变形与断裂参数的系统程序并考虑应变率与温度效应，所得模型与仪器化Charpy及落锤撕裂试验数据比较并在不同于参数确定所用温度下评估——此类验证尤为重要，因为与用于标定的试验吻合本身并不能证明模型在不同条件下保持预测能力。尽管在需显式表征不同失效机制时可能需要替代断裂形式，Johnson–Cook框架仍因相对简单地表征应变硬化、应变率效应、热软化、损伤累积与断裂而被广泛用于数值Charpy模拟。试样几何也已通过经试验验证的FEM与Johnson–Cook模型加以考察，例如采用Altair Radioss模拟不同轧制方向提取的亚尺寸S600MC试样，将Johnson–Cook塑性形式与延性失效准则结合并通过比较数值预测与试验测量吸收能量标定失效参数；或采用ABAQUS与Johnson–Cook模型研究试样宽度、侧槽深度与侧槽根部半径对压力容器钢冲击响应的影响，比较力-位移曲线、吸收能量演化、断口形貌与总冲击能量不同贡献，结果证实几何变化不仅改变最终吸收能量，也改变冲击过程中力与损伤的演化。综上，这些研究在Charpy冲击试验数值再现方面取得实质性进展，显式FEM已被用于研究不同钢种、试样几何、加载速度、冲锤构型、温度与断裂 regime，验证采用吸收冲击能量、力-位移历程、断口形貌及某些情况下附加冲击试验。然而标定策略仍高度异质：一些研究独立确定材料参数并以Charpy试验验证，另一些则直接针对Charpy响应调整一个或多个参数；本研究采用的方法属于后一类，但引入预备敏感性分析以在标定前识别对选定响应影响最强的参数，所得参数集因而未必具有相同物理意义或可移植程度。进一步局限在于吻合常主要通过最终吸收能量评估，该标量量重要但不能独立约束模拟中涉及的所有本构与损伤参数，不同参数组合可再现相近吸收能量却预测不同峰值力、位移历程、损伤萌生条件或断裂演化；纳入力-位移曲线与断口形貌的研究提供了更强验证基础，但一般未量化各模型参数的相对贡献或标定前参数间交互。因此，尽管本研究采用峰值冲击力作为敏感性分析与标定的主要响应，所得模型随后针对冲击响应的附加特征进行评估。该问题对亚尺寸试样尤为相关，因为厚度与缺口几何变化改变约束、应力三轴度、变形与断裂演化，针对单一构型标定的参数集可能再现该特定试验却不能对另一几何保持相同预测能力——这正是本研究通过考虑不同试样厚度与缺口深度所考察的可移植性问题。尽管已有若干数值Charpy模型，参数敏感性、统计设计数值试验与试验-数值标定的系统集成仍然有限，这一空白促使发展能够在多试样构型标定前识别最强支配模拟Charpy响应不同特征的Johnson–Cook本构与损伤参数的程序。总体而言，文献表明FEM可再现Charpy冲击响应的若干特征，但本构与损伤参数的标定常基于针对有限数量试验量的迭代调整，而未首先量化单个参数及其交互的相对影响；此外，针对特定试样标定的参数集应用于不同几何构型的亚尺寸试样时可能不保持相同预测能力。这些局限表明需要一种结合敏感性分析、统计评估与多试样几何评估的结构化标定程序。据此，本研究发展了一种试验-数值程序，以改善亚尺寸AISI 304试样仪器化Charpy冲击试验与LS-DYNA R16.1中显式FE模拟之间的相关性。试验方案考虑不同厚度与缺口深度的试样，从而提供考察几何影响与评估标定材料参数可移植性的基础。采用两水平全因子DOE量化六个选定Johnson–Cook本构与损伤参数的影响，各参数相对于基准值变化±20%，形成包含64次FE模拟的2⁶设计。DOE采用参考试样构型进行，选择峰值冲击力作为响应变量；随后采用一般线性模型与统计分析确定所考察因子的显著性并识别对数值响应影响最大的参数。敏感性结果随后用于指导最具影响Johnson–Cook参数针对试验响应的标定，而影响较小的参数保持基准值。所得参数集首先对参考构型评估，随后在无额外几何特定再标定情况下应用于其余亚尺寸试样。试验-数值相关性通过比较力-时间、力-位移、冲锤速度与吸收能量历程评估，特别关注峰值力、整体响应演化与断裂行为。本研究目标不是为所测AISI 304批次识别唯一的固有Johnson–Cook参数集，而是评估基于文献的参数集能否通过敏感性引导标定系统调整以再现特定Charpy冲击响应并随后迁移至不同试样几何。因此本研究的原创贡献在于将仪器化Charpy试验、显式LS-DYNA R16.1模拟、全因子DOE与统计分析集成于敏感性引导标定框架内。除改善参考构型的相关性外，该程序还使考察单一标定Johnson–Cook参数集在具有不同厚度与缺口深度的亚尺寸试样间的可迁移程度以及识别几何相关局限仍存在的响应特征成为可能。

**材料与方法**　本研究采用的试验-数值方法组织为七个主要阶段：（i）通过仪器化Charpy冲击试验对具有不同厚度与缺口深度的亚尺寸AISI 304试样进行试验表征，测得力、位移、速度与吸收能量历程用作参考数据；（ii）在LS-DYNA 2025R1中建立显式FE模型，包括几何、离散化、边界、加载与接触条件；（iii）定义Johnson–Cook本构与损伤模型并选择待考察参数，其余材料属性与模型参数保持固定；（iv）实施包含64次模拟的两水平全因子DOE，随后采用一般线性模型与方差分析（ANOVA）量化六个选定Johnson–Cook参数对峰值冲击力的影响；（v）利用参考试样构型的试验响应对选定参数进行敏感性引导标定；（vi）基于力-时间、力-位移、冲锤速度、吸收能量与变形/断裂响应开展试验-数值比较与验证；（vii）通过在无进一步几何特定调整情况下将同一参数集应用于其余试样构型评估所得标定模型，该最后阶段用于考察标定Johnson–Cook模型在不同试样几何间的可移植性与局限。AISI 304奥氏体不锈钢因广泛工程应用与充分记录的力学行为而被选用，其显著应变硬化响应与应变率敏感性使该合金特别适合评估非线性变形与断裂在动态冲击加载下如何演化。材料与试样制备遵循先前试验研究建立的一般程序，其中 machining 相关表面条件、试样清洁与处理程序被系统考察，本研究中该程序针对每个特定几何进行调整。研究考察四种亚尺寸试样构型，组合两种试样厚度与两种缺口深度；采用不同几何构型旨在提供具有不同韧带尺寸与约束条件的试验响应，并随后评估单一标定Johnson–Cook参数集能否再现不同试样几何的冲击行为。每种试样构型考虑单次试验冲击试验，该选择与本研究的目的一致，即不是再现多个名义相同试样的平均试验响应，而是评估FE模型再现单个冲击事件响应的能力；因此试验-数值比较针对每个测试试样的测得力、位移、速度与吸收能量历程进行，该途径提供模型预测能力的个案评估，尽管不能统计量化重复试验相关的试验变异性。试验Charpy冲击试验在意大利费拉拉大学工程系采用CEAST Resil Impactor（代码6957.000，CEAST，Pianezza，意大利）进行，该设备配备手动制动系统且最大冲击能量为50 J；摆锤冲击系统设计符合包括ASTM E23在内的国际冲击试验标准，摆锤从150°初始角释放，对应冲击速度3.458 m/s。冲击数据采用集成于冲击锤头的力传感器仪器化锤采集，采集系统每次冲击事件记录多达16,000个数据点，从而表征加载、变形与断裂过程中试样响应的演化；记录量包括冲击力、冲锤位移、冲锤速度与吸收能量，这些试验历程随后用作数值模型开发、标定与评估的参考响应。每次冲击试验提供力[N]、吸收能量[J]、冲锤位移[mm]、冲锤速度[m/s]的完整时间历程以及描述冲击事件的一组标量参数，即非仪器化吸收能量与锤势能、线性与常规回弹值、能量百分比与速度变化；从仪器化响应确定峰值力F_peak、峰值力处吸收能量E_peak与峰值力处位移δ_peak，以及总吸收能量E_total与总位移δ_total，这些量提供每次冲击事件主要阶段与总体严重程度的紧凑描述，并随后在相关处用于试验-数值比较。试验Charpy冲击试验采用显式FE求解器LS-DYNA 2025R1对每个试样构型进行数值再现；该求解器因其适合具有复杂接触相互作用、大变形与渐进损伤与断裂的瞬态冲击问题而被选用；每个FE模型包括试样、冲锤与支座，数值配置定义为再现相应试验装置。试样几何再现每个试验构型的尺寸，采用减缩积分三维实体单元（ELFORM 1）离散，在U型缺口区域进行局部网格细化；缺口轮廓分为两个侧向段与一个中心曲线段，每个侧向段对2 mm缺口构型（A1与A3）离散为10个单元、对5 mm缺口构型（A2与A4）离散为20个单元，而中心曲线段对所有构型离散为20个单元；该细化旨在更好地分辨断裂萌生与扩展期间预期的高变形与损伤梯度。冲锤与支座采用三维实体单元建模为刚体，在参与接触的区域引入附加网格细化；冲锤冲击区域离散为35个单元，两个支座接触区域分别离散为11与15个单元；将这些部件视为刚体使计算量集中于可变形试样，同时保持表征冲击事件所需的几何与离散化。边界条件定义为再现试验Charpy冲击配置：支座完全约束，冲锤限制为沿冲击方向的平动；冲锤施加3.458 m/s初始速度，对应试验冲击条件，其有效质量为8.36 kg，产生50 J初始冲击能量；数值分析总模拟时间为20 ms以捕捉试样完整变形与断裂响应。显式分析期间数值稳定性采用沙漏控制（IHQ = 5，QM = 0.01）与时间步缩放因子（TSSFAC）0.9维持，未应用人工质量缩放（DT2MS = 0）；材料失效由针对试样实施的Johnson–Cook损伤形式控制。在可变形试样与刚性冲锤及支座之间定义接触相互作用以再现Charpy冲击试验期间发生的力学相互作用；冲锤-试样与试样-支座界面均采用LS-DYNA *CONTACT_AUTOMATIC_SURFACE_TO_SURFACE形式与无摩擦接触（μ = 0）建模，允许切向滑动同时在相互作用表面间传递法向接触力；该建模假设与先前钢试样Charpy冲击试验数值研究一致，其中采用无摩擦接触条件表征试样、冲锤与支座之间的相互作用。未采用侵蚀接触形式，因为接触仅用于描述试样、冲锤与支座之间的相互作用；材料失效与单元删除由Johnson–Cook损伤模型直接处理，而新产生的断裂表面或分离碎片之间的接触未被视为对整体冲击响应重要。材料分离改由根据Johnson–Cook损伤准则的单元失效控制；相同接触形式与参数对所有试样构型以及整个DOE与标定分析保持，从而使数值响应的变化可归因于所考察的材料参数而非接触定义的变化。Johnson–Cook本构与损伤形式通过*MAT_015/MAT_JOHNSON_COOK在LS-DYNA中实施，该模型将等效应力描述为塑性应变、应变率与温度的函数，其中各术语分别考虑应变硬化、应变率与热软化对材料响应的贡献；材料失效采用Johnson–Cook损伤形式描述，其中失效时的等效塑性应变表示为应力三轴度的函数，D₁至D₃为Johnson–Cook损伤参数，具体地D₁、D₂与D₃支配断裂应变对应力状态的依赖，而D₄与D₅分别考虑其应变率与温度依赖；损伤根据方程在塑性变形期间增量累积，材料失效假定在累积损伤变量达到临界值时发生。AISI 304的Johnson–Cook参数代表性文献范围从先前涉及本构与断裂行为的研究汇编；具体地，本构参数主要源自Jiang等报道的比较文献分析，而损伤参数范围由Seo等与Frontán等的断裂表征研究建立；这些范围表明本研究模型采用的基准本构参数集直接源自Tatarintsev等，处于文献报道的变异性之内；同时若干Johnson–Cook参数报道的相对宽范围突显其识别相关的实质变异性以及 consequently 定义适合模型标定的参数空间的复杂性，此类变异性部分反映各研究考虑的试验条件、材料状态、应变率与温度范围、试样几何与识别程序的差异。因此应注意这些文献范围不应必然被解释为对本文所考察特定条件同样合理的标定区间；更有针对性的文献数据选择可通过排除在明显不同于本文Charpy冲击试验条件下获得的值而实质性缩减可接受参数空间，例如从 machining 相关试验或从强烈改性或纳米结构材料状态识别的参数集可能较不能代表块体AISI 304在室温冲击加载下的响应；在不进行逐源详细比较的情况下，本研究在选择用于后续DOE与标定分析的基准参数集时进行了该附加文献筛选。数值响应从LS-DYNA生成的时间历程输出评估：接触力历程从界面力输出（RCFORC）获得，冲锤运动学响应通过刚体输出（RBDOUT）监测，全局与材料能量历程分别从GLSTAT与MATSUM获得；这些时间历程量以1 × 10^?5 s（0.01 ms）间隔获得；选择这些数值输出以提供与试验信息直接可比的响应量，包括冲击力与能量响应；此外，试样在整个冲击事件期间的变形与断裂演化从时间相关FE构型评估。选择DOE方法量化选定Johnson–Cook参数对数值Charpy冲击响应的影响；六个Johnson–Cook参数A、B、n、D₁、D₂与D₃被选为敏感性分析的可变因子，而C、m、D₄与D₅在整个DOE期间保持其采用的基准值；这些选定参数共同支配材料的初始屈服、应变硬化行为与应力三轴度相关失效响应，因而代表敏感性分析中考虑的主要本构与损伤项。关于未选参数，应变率系数C与热软化指数m保留其基于文献的基准值，因为试验方案未包括冲击速度或温度的独立变化； consequently，尽管应变率与热效应可能在冲击期间对局部材料响应有贡献，可用试验数据未提供足够独立信息用于其可靠单独识别；将这些参数作为附加DOE因子将因此增加参数空间维数与非唯一或补偿参数组合的可能性，而不提供相应试验基础用于其独立识别。支配失效应变应变率与温度依赖的损伤参数D₄与D₅也被排除在DOE之外，其采用基准值设为零，使失效应变形式中相应的速率与温度相关贡献在模型中 inactive；文献中也提出了基于Charpy的程序确定应变率相关损伤参数D₄。应强调将C与m保持基准值并不意味着在本构形式中忽略应变率敏感性或热软化，这些贡献仍由采用的Johnson–Cook参数支配，尽管其个体敏感性未被考察；相反，将D₄与D₅设为零从失效应变形式中移除显式应变率与温度依赖；采用该简化是因为试验方案在单一冲击速度与相同名义温度条件下进行，无适合识别这些附加损伤参数的独立测量； consequently，应变率与温度对缺口附近局部损伤与断裂演化的可能耦合效应未被完全表征，代表本文数值模型的局限。然而缺口附近局部应变率分布未被独立量化，峰值冲击力与吸收能量对C变化的敏感性未被评估，因此采用的C值在局部应变率范围内的适用性仍是本研究的附加局限。此外，塑性变形相关的绝热温升未在数值模拟中考虑；尽管热软化参数m保留其基于文献的值，温度演化未被显式建模； consequently，可能局部热软化效应，特别是在塑性变形集中的缺口附近，未被捕捉；该简化可能影响冲击后期阶段预测的流变应力、力-位移响应与损伤演化；因此，缺少热力耦合代表本文模型的附加局限。每个选定因子在两个水平评估，对应相对于其基准值的?20%与+20%；采用两水平全因子设计（2⁶），产生64个唯一参数组合，每个组合进行一次数值模拟；基准值仅报告供参考，未作为附加因子水平纳入因子设计。比较采用的Johnson–Cook基准值及其±20%区间与代表性文献范围（归一化于0与1之间），多数参数落在报道范围内：A位于较低部分；B、m与D₁位于较高部分，而n、C与D₂占据中间位置；相反，D₁落在代表性文献范围之外（0.04 vs. 0.69–1.21）；然而D₁、D₂与D₃被联合识别并共同定义断裂应变对应力三轴度的依赖，其个体值可能因此随用于识别的三轴度范围与加载条件而实质性变化；原始D₁值因而与关联的D₂与D₃值一起保留以保持损伤参数集的一致性。在采用的Johnson–Cook失效方程中，D₁是断裂应变的加性贡献，而D₂与D₃确定其应力三轴度依赖；因此，个体基准D₁值与文献区间之间的差异本身并不确立组合失效应变函数的有效性或无效性；然而报道区间不能独立证实D₁ = 0.04：该值并非从本文AISI 304批次的断裂试验识别；所选三元组因而应被视为临时的、耦合的参数化，与独立识别的断裂应变-三轴度曲线的比较将需要用于材料特定验证。DOE采用A1试样作为参考情况；对64个参数组合中的每一个，从相应数值模拟提取峰值冲击力（F_peak）并用作响应变量；所得敏感性分析因而具体针对A1构型与所考察±20%参数域。全2⁶因子设计采用一般线性模型（GLM）分析，包括六个所考察Johnson–Cook参数的主效应及其两因子交互；主效应量化F_peak在各参数低水平与高水平之间的平均变化，而交互项识别一个参数的效应是否依赖另一个参数。统计显著性用于导出缩减的层级模型，随后应用方差分析（ANOVA）量化保留项的相对贡献；主效应图用于识别参数影响的方向与量级；所得敏感性排序随后用于指导而非直接确定Johnson–Cook参数的标定。所得敏感性排序用于优先后续标定而非直接确定标定参数值；对F_peak显示最强影响的参数因此在标定期间被给予主要考虑。在基于DOE的敏感性分析之后，进行第二阶段以针对A1参考试样的试验响应标定Johnson–Cook模型；DOE采用的±20%区间未保留为标定约束，因为其仅被引入以定义敏感性分析的受控域；在标定期间，最具影响参数因而被允许相对于其基准值变化至±30%。参数调整由敏感性排序与DOE分析识别的效应方向共同指导；随后将附加数值模拟与试验A1响应比较以评估修改的参数组合是否改善试验-数值吻合；因此，DOE用于指导和限制标定搜索而非直接确定标定参数值。最终参数集从所考察组合中基于A1的数值与试验响应之间的吻合选择；一旦建立，同一参数集在无进一步几何特定调整情况下应用于其余试样构型以评估其可移植性。本方法的一个关键方法学方面是固定模型输入与标定期间视为可调参数之间的区分：表中报告的物理与弹性属性在整个数值分析中保持不变，而Johnson–Cook本构与损伤参数从文献数据初始化而非从本文试验所用材料批次直接识别；将标定限制于非线性本构与损伤描述缩减反演问题的维数并限制通过无关模型量的补偿变化获得相似全局响应的可能性；采用全场测量与有限元分析的联合试验-数值程序也已被用于金属材料Johnson–Cook参数确定。敏感性分析限于根据可用试验条件与采用基准值选择的六个Johnson–Cook参数；该限制通过将分析集中于预期影响选定响应的参数并减少潜在补偿参数组合数量而改善标定问题的可解释性。峰值冲击力被选为主要响应变量，因为它提供明确定义的标量量，可从试验与数值响应一致提取并在所有DOE模拟间比较；其使用也避免在敏感性分析阶段引入与峰后行为、损伤萌生或断裂进展解释相关的附加准则；该选择为因子设计的统计分析提供实用且可再现的基础；然而该简化带来固有局限：对F_peak的敏感性未必反映对冲击动力学其他方面（即吸收能量、峰后软化、变形或断裂演化）的敏感性；不同参数组合因而可再现相似峰值力却产生不同后续响应。同理，F_peak被用作敏感性分析与标定指导的主要量，而标定模型的更广泛评估考虑完整力、位移、速度、能量与断裂响应。然而敏感性结果必须针对所选响应变量解释：本研究中DOE具体量化参数对峰值冲击力的影响；对F_peak显示有限敏感性的参数仍可能影响响应其他方面，包括峰后软化、损伤演化、吸收能量、冲锤运动学或断裂进展；因此DOE用于指导标定，而标定模型的后续评估考虑附加试验-数值响应量与变形/断裂行为。据此，所提程序应被视为敏感性引导标定策略而非唯一固有Johnson–Cook参数集的识别；此外，标定参数不应被解释为从专门力学表征试验（如拉伸或高应变率试验）获得的材料参数的替代；所得参数集最后通过在无几何特定再标定情况下应用于其余试样构型评估。

**研究结果**　3.1 试验Charpy冲击响应：从A1–A4试样仪器化Charpy冲击试验获得的主要特征量包括峰值力、峰值力处能量、总吸收能量、峰值力处位移与总位移。试验结果揭示四种试样构型之间的显著差异，不仅在于峰值力，还在于力、位移与吸收能量之间的关系。最大峰值力为试样A3（9751.220 N），而最高总吸收能量为A4（49.638 J）；值得注意的是，A3与A4表现出相对相似的总吸收能量（46.997与49.638 J），尽管其峰值力显著不同（9751.220与4593.948 N），这表明可比水平的总能量吸收可与显著不同的力与变形历程相关联。试样A3特征为最高峰值力与最小峰值力处位移（1.990 mm），连同相对较小的总位移（4.943 mm），该组合表明响应由在相对有限变形范围内达到的高载荷水平支配；相反，A4达到低得多的峰值力但表现出最高峰值力处能量（20.614 J）、最高总吸收能量与相当大的总位移（12.731 mm），表明更变形密集的能量吸收响应。A1显示中间行为，峰值力5320 N与总吸收能量35.753 J；约56%的最终吸收能量在达到最大力时已累积（E_peak/E_total ≈ 0.56），为所考察构型中最高比例；相比之下，相应比率对A2约为0.33、对A3约为0.32、对A4约为0.42，表明能量吸收在力峰值前后分布方式的实质差异。试样A2表现出最低峰值力（2132.587 N）与最低总吸收能量（20.131 J），同时显示最大总位移（20.355 mm）；低力与能量水平与相对大位移的共存进一步表明整体冲击响应不能仅由峰值力或吸收能量表征。总体而言，这些试验结果突显Charpy响应对试样构型的强烈依赖，如预期，并表明力、位移与吸收能量提供互补信息；具有可比总能量或相似力水平的试样间观察到的实质差异支持在评估后续数值模拟与解释标定Johnson–Cook模型的有效性与可移植性时使用多个响应量。3.2 冲击动力学：数值模型预测的A1试样渐进变形与代表性情况（即基准配置）的相应力-时间响应。在阶段（a），冲击力迅速增加而变形开始在U型缺口周围局部化；在阶段（b），试样接近断裂萌生，力达到其最大值附近且变形在剩余韧带中强烈集中；阶段（c）对应峰后 regime，裂纹扩展明确建立且试样逐渐失去承载能力；最后在阶段（d），裂纹扩展推进且力已实质性下降，试样接近完全断裂。3.3 数值输出分析：展示当单个参数围绕其基准值变化时A1参考试样数值力-时间响应可能如何变化；在每种情况下，仅一个参数修改±20%，而所有其余模型参数保持固定；该单因子逐一表示仅旨在说明单个参数对模拟响应形状与量级的特征效应，不应被解释为全因子DOE的替代，其中所有参数组合与交互同时考虑；本构参数A、B与n主要影响力加载阶段与力峰值的量级与位置，而损伤参数预期对峰后响应、损伤演化与承载能力丧失施加强烈影响；这些效应的实际量级以及参数间可能的交互通过全2⁶因子设计与后续GLM与ANOVA分析量化。在±20% DOE域内，提供最高峰值力的参数组合产生N（即模拟n. 28），仍比试验A1目标5320.3 N低约4.2%；以下敏感性趋势将验证参数的进一步变化是否会使数值响应向试验值移动；这为将最具影响参数的标定范围扩展至±30%提供了基础。3.4 基于DOE的敏感性分析：试验设计分析表明峰值冲击力表现出对Johnson–Cook模型若干本构与失效参数的显著依赖。缩减层级模型达到R²为95.31%、调整R²为94.82%、预测R²为94.09%，表明高度表征观测响应变异性的能力。方差分析（ANOVA）表明A、B、n、D₂与D₃参数对峰值冲击力具有统计显著效应（α = 0.05）；B×n交互也显著（α = 0.05）；如标准化效应Pareto图所示，所有保留项均超过对应α = 0.05的显著性阈值。此外，表中呈现显著项对响应的百分比贡献；B×n交互也统计显著（α = 0.05），表明这两个参数对峰值冲击力的效应相互关联；相反D₁在所考察范围内未表现出统计显著效应，因而从缩减层级模型中排除。具体地，n参数对峰值冲击力变异性表现出最大百分比贡献（28.11%），其次为B（27.13%）与A（24.58%），识别它们为所考察范围内对响应影响最大的三个参数；失效参数D₃与D₂显示较低贡献分别为9.27%与5.23%，尽管两者均保持对峰值冲击力的统计显著效应；B×n交互表现出1.00%的贡献与F值12.12（α = 0.05），表明这两个参数之间的显著交互。主效应与交互图总结显著Johnson–Cook参数及B×n交互对峰值冲击力的主效应；在主效应与交互图中，垂直轴代表DOE内预测的平均峰值冲击力，对其他所考察因子水平平均，因此不应被解释为试验目标值；增加A与B产生响应的显著增加，而增加n导致显著下降；损伤参数D₂与D₃也显示正效应，尽管不如本构参数显著；峰值冲击力在D₂两个水平下均随D₃增加而增加，但该增加的幅度在D₂较低水平更大，如交互图中非平行趋势所示；这证实GLM/ANOVA分析识别的统计显著B×n交互；总体而言，这些结果提供参数效应的相对重要性与方向，从而指导Johnson–Cook模型的后续标定。上述ANOVA排序具体适用于所考察±20%参数范围内A1试样的F_peak；然而当考虑其他响应量时敏感性层级可能不同；特别地，总吸收能量依赖完整力-位移历程，包括峰后阶段，并可能因此对支配损伤与断裂的参数更敏感；为考察该方面，同一64次运行DOE数据集进一步分析考虑能量相关响应。为进一步考察所选响应变量对敏感性排序的影响，使用同一64次运行全因子DOE结果进行补充分析；除峰值冲击力（F_peak）外，考虑两个能量相关量：峰值力处吸收能量（E_peak）与总吸收能量（E_total）；对每个响应，六个Johnson–Cook参数的主效应从?20%与+20%因子水平处平均响应之间的差异计算；其相对贡献随后针对相应响应的总变异性评估；该附加评估旨在识别取决于所选标定目标应优先考虑哪些参数，而无需进一步数值模拟。结果表明Johnson–Cook参数的相对重要性随所选响应变量显著变化；对峰值冲击力，本构参数A、B与n占主导，合计约占响应总变异性的87.8%，这证实先前ANOVA分析识别的一般敏感性层级；然而当考虑能量相关响应时观察到不同行为：对E_peak，损伤参数D₂与D₃约占74.4%的总变异性，而其合计贡献对E_total增加至约81.9%；相反A与n的影响变得显著更小；参数B对两个能量相关响应保持中等贡献。结果还表明D₁在所考察±20%范围内对F_peak、E_peak与E_total影响有限；因此，其调整不能仅由这些全局响应量的显著敏感性证明，尽管其对局部损伤萌生与断裂演化的可能影响不能在没有专门分析的情况下排除。这些发现为某些试样构型峰后响应中观察到的差异提供可能解释：峰值力主要受本构参数A、B与n控制，而能量吸收更强烈受损伤参数D₂与D₃影响； consequently，主要集中于F_peak的标定可能不足以约束整个能量吸收过程；然而局部断裂萌生与扩展的敏感性在本分析中未被直接评估，将需要附加断裂特定响应量。尽管如此，当前结果代表补充敏感性评估而非新的多目标标定，因为最终参数集使用原始标定程序获得。3.5 Johnson–Cook参数的敏感性引导标定：由识别的敏感性趋势指导，通过沿预期改善与试验A1响应吻合的方向渐进修改参数进行标定；增加A与B并减小n渐进改善峰值冲击力预测，最佳吻合分别对+30%、+30%与?30%变化获得。因为这些三个参数在所考察DOE范围内显示可比敏感性与平行趋势，故意采用相等相对变化幅度，符号由每个参数效应方向确定；±30%变化定义为相对于基准值的最大允许偏差以约束标定过程并避免开放式参数搜索，其中相关参数间补偿变化可导致渐进不同但未必更可识别的参数组合。损伤参数D₁、D₂与D₃同时调整+20%以探索其对损伤演化与峰后行为的联合影响；尽管D₁在所考察DOE范围内未表现出对峰值冲击力的统计显著影响，该结果特定于所选响应变量且未必延伸至其他损伤相关量；的确，D₁直接贡献Johnson–Cook损伤形式中的失效等效塑性应变，并可能因此影响断裂萌生与后续材料分离；其调整因而为探索性而非由峰值力敏感性驱动；然而由于三个损伤参数同时修改且未对峰后或断裂相关响应进行单独敏感性分析，D₁对标定响应的个体贡献不能从当前结果建立；最终标定参数值报告于表8。更具体地，D₁在DOE内对F_peak的可忽略效应不证明使用其+20%调整来校正该响应；调整为探索性并与D₂与D₃同时实施；可用结果不能识别归因于D₁的单独改善；据此，D₁ = 0.048不应被解释为唯一标定或试验验证的损伤参数，保留基准值仍是合理的替代方案，待独立损伤响应证据。3.6 试验-数值比较与模型可移植性：标定后试验-数值相关性对所有试样（A1–A4）通过比较（a）力-时间、（b）力-位移与（c）冲锤速度与吸收能量响应评估。A1：数值模型令人满意地再现试验响应的主要特征，包括初始加载阶段与峰值力；峰后阶段出现中等差异，数值力下降更快；冲锤速度与吸收能量也良好再现，仅在冲击事件末期有小偏差。A2：标定参数集与试验响应提供良好吻合，特别是峰值力与整体力演化；观察到最终吸收能量的中等低估与冲锤速度的略大数值下降；总体而言，基于A1的标定保持合理可迁移至A2几何。A3：数值响应在可可靠测量的冲击事件部分与试验数据良好吻合；在较高力水平，试验力传感器超过其测量范围，后续力信号因而不再能被认为代表实际试样响应；试验与数值力历程之间的表观差异应因此主要归因于该测量局限而非数值模型的缺陷；模拟表明试样不断裂并继续经历变形与卸载，提供对未能试验捕捉的响应部分的物理一致重构；冲锤速度与吸收能量历程中观察到的吻合进一步支持该解释。A4：数值与试验响应在峰值力附近前合理吻合；峰值后数值力比试验响应下降稍快，表明预测峰后损伤演化的差异；然而模拟正确预测无试样断裂，与试验观察一致；吸收能量与冲锤速度响应也显示令人满意吻合，支持标定模型再现整体冲击响应与A4几何观察到的非断裂条件的能力。下降分支的差异，特别是对A2与A4，突显采用单响应标定的局限：标定在A1上使用F_peak作为主要目标进行，同一参数集在无进一步调整情况下迁移至其他几何；匹配F_peak不决定峰后斜率、后续位移期间吸收的功量或局部损伤发展的时刻与模式；剩余差异因此不应被解释为断裂演化已独立验证的证据；相反它们识别未来多响应标定需要包括的附加约束，如总吸收能量、峰后曲线误差、特征事件处位移与直接观察的断裂特征。网格离散对力-时间响应的影响：采用四种有限元配置（A1a、A1b、A1c与A1d）进行网格离散敏感性分析，考虑峰值冲击力作为主要评估参数；考察名义单元尺寸1.2与0.5 mm，连同用于离散U型缺口直边与圆角段的单元数量变化；结果表明配置A1b与A1c表现出相似数值响应，峰值力分别为5.381,86与5.386,71 N，对应相对差异仅0.09%；该吻合表明峰值力对这两个配置间离散差异的低敏感性；然而配置A1b提供与试验峰值力5308.83 N最接近的吻合，产生1.38%相对误差，相比A1c的1.47%；此外，A1b沿直边采用10个单元、沿缺口圆角段采用20个单元，提供比仅在直边区域采用2个单元的配置A1c与A1d更详细的直边区域几何表示；该改善的局部分辨率提供捕捉与缺口几何相关变形梯度与应力集中的潜在优势；相反，配置A1d表现出3.18%试验偏差，而A1a导致5.62%误差，表明峰值力预测较低吻合；因此，配置A1b被选为所评估替代中最适合的离散，因为其实现相对于试验峰值力的最低误差同时提供沿缺口直边更详细的局部离散；尽管A1b与A1c之间的紧密吻合表明峰值力稳定性，建立严格网格收敛将需要附加系统网格细化。3.7 模型性能定量评估：为补充试验与数值响应历程的定性比较，表10总结四种试样构型的总吸收能量、峰值冲击力与最终冲锤速度的试验与数值；对总吸收能量与峰值力计算绝对相对误差（|Δ|），而最终速度值直接报告以提供冲击后冲锤响应的附加信息。尽管标定专门针对A1参考试样的峰值冲击力进行（精度99.70%），数值模型在所考察构型间以总体良好精度预测其他关键响应量；对试样A1、A2与A4，峰值冲击力以0.30%至3.40%的误差再现；总吸收能量的吻合也令人满意，A1–A4误差分别为11.11%、5.00%、0.04%与2.38%；该吻合尤其相关，因为吸收能量与冲锤速度均未被直接用作标定目标。最终冲锤速度提供冲击后行为的附加信息；对A1与A2，数值模型紧密再现试验最终速度；对A3与A4，试验与数值最终速度均为负，表明摆锤回弹与无试样断裂；然而最终速度的定量差异对A3（23.88%）比对A4（7.88%）更显著。特别地，A3构型需要特别注意，因为试验力传感器在冲击试验期间达到其测量极限； consequently，高力水平处试验力-时间响应不可靠，所报道16.73%的峰值力误差不能被认为模型精度的有效度量；同理，A3不能用于完整力演化或损伤萌生时间的可靠定量验证；然而吸收能量、最终冲锤速度与试验观察到的无断裂仍可用于评估全局冲击响应；因此A3结果应被视为模型的部分验证而非力与损伤演化的完整验证。3.8 变形与断裂行为的数值预测：比较A1–A4试样预测的最终变形配置；数值结果表明不同厚度与缺口深度组合导致显著不同的变形与损伤模式；在A1与A2中，变形在缺口处强烈局部化并演化为裂纹通过剩余韧带扩展，导致试样断裂；相反，A3与A4表现出显著弯曲与缺口周围局部化损伤而无完整韧带分离，与试验观察到的非断裂条件一致；A3显示最有限的整体变形，而A4经历实质性更大弯曲变形同时仍保持试样连续性；这些结果证实试样几何不仅影响力与吸收能量的量级，也影响变形模式与断裂与非断裂行为之间的转变。标定模型使用同一Johnson–Cook参数集再现断裂与非断裂构型的能力提供其跨所考察试样几何可移植性的进一步证据。

**讨论与结论**　本研究发展了一种敏感性引导程序，用于针对AISI 304钢的仪器化亚尺寸Charpy冲击试验标定Johnson–Cook本构与损伤模型。主要假设、方法学成果与一般发现可总结如下：方法学假设与范围——标定基于单一参考构型A1，选择峰值冲击力作为主要响应变量；通过将A、B、n、D₁、D₂与D₃围绕其基准值变化±20%进行两水平全因子DOE，而其余材料属性与Johnson–Cook参数保持固定；DOE范围用于表征参数敏感性且未被施加为后续标定的约束；所得参数集因而应被解释为所考察冲击条件的标定集而非固有材料常数的唯一识别。敏感性分析识别受限的主导参数集——峰值冲击力主要受n、B与A支配，贡献分别为28.11%、27.13%与24.58%；D₃与D₂显示较小但统计显著的贡献，而D₁在所考察DOE范围内未显著影响峰值力；还识别了统计显著的B×n交互。DOE提供参数效应的相对重要性与方向：增加A与B增加预测峰值冲击力，而增加n降低它；D₂与D₃产生较弱正效应；这些趋势用于指导标定而非直接确定最终参数值。标定被故意约束以限制参数补偿——因为A、B与n表现出可比敏感性与近似线性主效应趋势，对这些参数采用相同相对变化幅度，符号由其效应方向确定；标定极限设为相对于基准值的±30%以避免开放式调整过程，其中补偿参数变化可导致日益不同但非唯一的参数组合；最终标定导致A = +30、B = +30与n = ?30，而D₁、D₂与D₃修改+20%以改善损伤演化与峰后响应的表征。可移植性实质好于单响应、单试样标定可能预期的——尽管标定仅使用A1的峰值冲击力，所得参数集再现所有四种构型的峰值力，误差约0.30%至3.40%；总吸收能量虽未直接用作标定目标，也以令人满意精度预测，A1–A4误差约11.11%、5.00%、0.04%与2.38%；未明确纳入标定的量所获得的吻合支持标定参数集跨所考察几何的可移植性。模型还再现主要冲击后行为——对A3与A4，试验与数值响应均表明摆锤回弹与无试样断裂；在A3中，试验力传感器在高力水平超过其测量范围； consequently，数值模拟提供对未能可靠测量的冲击响应部分的物理一致重构；冲锤速度与吸收能量的吻合进一步支持该解释。剩余差异主要集中于详细峰后响应而非整体冲击行为——力软化、损伤演化与最终冲锤运动学的差异表明仅准确再现峰值力不能完全约束动态损伤响应；峰值力因而对敏感性分析与标定指导是有效标量量，但需要附加响应量以更完整识别断裂行为。重要的是，参数排序是响应特定的：对F_peak证明的A、B与n主导不能推广至总吸收能量或断裂度量；初始D₁ = 0.04尽管在耦合D₁–D₂–D₃失效律内使用，却落在表2引用的个体文献区间之外，且其值及其后续+20%调整均未被峰值力DOE独立验证；未来分析应在采用多目标标定前量化对E_total与损伤相关响应变量的敏感性。总体而言，结果表明结合全因子DOE、统计敏感性分析与针对性数值标定为纯试错参数调整提供透明且计算上可管理的替代方案；一个特别相关的结果是，对单一试样单一响应量进行的标定能够以总体良好精度再现若干独立响应量与不同试样几何。未来发展应将程序扩展至多目标标定，同时纳入力演化、吸收能量、冲锤运动学、变形与断裂相关量以进一步改善参数可识别性与模型稳健性。