《Materials & Design》：A density-dependent Drucker–Prager cap model for green powder compacts: identification, LS-DYNA user-material implementation and benchmark validation
编辑推荐：
摘要（据稿件开篇要点）译文：每个帽盖（cap）常数均为密度的函数，且每条密度律的形式均经选择确定。本构定律必须采用当前密度；若采用塑性密度，模具应力（die stress）代价最高可达32%。工业零件的生坯密度（green density）可在无需重新标定的条
摘要（据稿件开篇要点）译文：每个帽盖（cap）常数均为密度的函数，且每条密度律的形式均经选择确定。本构定律必须采用当前密度；若采用塑性密度，模具应力（die stress）代价最高可达32%。工业零件的生坯密度（green density）可在无需重新标定的条件下复现。在锥面–帽盖交线（cone–cap junction）以下，剪切不能致密化：这是帽盖模型内置的极限。
研究人员围绕生坯粉末压坯的密度相关Drucker–Prager帽盖模型（DPC，Drucker–Prager Cap model，一种由剪切包络与硬化帽盖描述粉末塑性致密化的本构模型）开展研究。冷压金属粉末零件的质量由密度分布决定，而密度分布受摩擦、工具运动学与本构响应共同控制；已有轮组对比表明，在材料与过程数据一致时，有限元代码选择并非主要误差源，误差更多归因于粉末表征及硬化曲线拟合。现有商业显式代码中的帽盖模型把各常数作为单一数值冻结，不能在密度范围内更新；同时，帽盖族把致密化仅赋予帽盖支，导致低围压剪切不能致密化的问题未被系统量化；此外，固定帽盖端点和低压力剪切包络的实验约束较弱。研究人员因此建立可辨识的密度相关DPC：每个常数写成相对密度（relative density）函数，函数形式逐参数从候选族中选择；将其实现为LS-DYNA用户材料子程序（UMAT，user-defined material subroutine），并与两种内置帽盖卡比较；最后用独立工业压坯基准检验无需重标定的可迁移性。该研究的意义在于，把文献中已报道的密度相关性落实为本构实现，明确冻结常数与连续帽盖形状各自造成的误差，并给出帽盖模型在剪切致密化上的构造性压力阈值。
方法上，研究人员采用两组已发表数据：扩散合金化铁粉Distaloy AE与WC–Co硬质合金，二者由同一协议表征。辨识试验包括径向与轴向巴西试验、闭模压缩（含径向应力测量）、等静压（isostatic compaction）及卸载–重载；将力–位移记录转换为屈服面应力点后，以最小二乘拟合DPC参数。每个参数用幂、指数、线性、二次、偏置指数或偏置幂律等候选律回归，以加权调整决定系数（adjusted coefficient of determination）选择，并用AICc（corrected Akaike information criterion）、BIC（Bayesian information criterion）、留一交叉验证（LOO，leave-one-out cross-validation）及重采样检验稳定性。UMAT以总体积应变推进当前相对密度，更新各常数，并执行弹性预测与回映（return mapping）；随后在四种加载程序与三种网格上做材料点验证，并用Dienet Case Study 2进行无重标定的工业基准计算。
结果部分按原文小标题概括如下。4.1 校准模型参数：对两种粉末各七个密度点进行参数辨识与密度律选择，发现多数强度与帽盖参数可由偏置幂律高精度复现，摩擦角可作为常数；硬质合金帽盖偏心率近似常数，而铁基帽盖偏心率随密度显著变化，说明函数形式不能预先固定；经典两参数硬化律对静水截距的重建存在形状限制，模量残差较其他参数更大。4.2 重建屈服面：通过比较密度相关用户材料与冻结的两种内置卡，得出静水极限和模量随致密化大幅增长；冻结常数造成的误差随密度变号，连续帽盖形状造成近似固定的乘性偏差，且该偏差由摩擦角和帽盖偏心率控制，不能靠调整偏心率完全消除。4.3 材料点验证：在静水压缩、闭模压缩、恒围压简单剪切与压–剪–卸载组合程序中，用户材料复现闭模轴向应力和侧压比（lateral stress ratio），而两种内置卡均出现明显偏差，连续帽盖返回的侧压比显著偏低，冻结椭圆趋势与实测相反；循环后内置卡保留残余剪应力，密度相关模型近于无应力回归；组合路径表明差异主要来自冻结常数，低围压下流法则起决定作用。4.4 独立闭模基准验证：用Dienet Case Study 2的带法兰衬套工业压坯，在不调整任何参数条件下复现轮毂及下部截面的生坯密度，下部截面优于基准公开参考模拟；上冲头和外下冲头载荷较七家独立帽盖实现更接近测量，但仍低于测量，作者将残差归于冲头运动学与摩擦等输入不确定性，并强调密度场才是该基准对本构验证的主要依据。
讨论部分指出，最大误差来源之一是密度定义：已发表律以承载状态下测得的密度为自变量，其中含弹性体积变化；若按塑性体积应变推进密度，会低估闭模应力，最高达32%，因此必须以总体积应变推进当前相对密度。评估还表明，内置卡的误差可分解为冻结常数误差和帽盖形状误差：前者随密度变号，只在局部可修；后者为固定乘性偏差。侧压比因直接关联壁摩擦、脱模力和回弹，是最有实验意义的判别量。局限性包括：相对密度是唯一状态变量，无织构和各向异性记忆；弹性按各向同性等效处理且为次弹性；不描述断裂和残余孔洞演化；静水端点与剪切支实验约束弱；模型为等温；帽盖构造导致低于帽盖–锥面交线压力时剪切不致密化，限制其超出压坯主导路径的应用。
研究结论部分译文：1. 研究人员针对两个公开数据集辨识了修正DPC模型，每个常数均为相对密度函数，且各律形式经选择确定；两种材料保留形式不同，且均需候选族原本缺少的偏置幂律，该形式可使参数在填充密度处归零。2. 驱动本构律的密度必须是当前相对密度，并由总体积应变推进，而非帽盖实现惯用的塑性密度；采用塑性密度会低估闭模应力，最高达32%，该定义应与密度律本身同等明确。3. 按此定义实现的密度相关用户材料在LS-DYNA中闭合了源数据回路：多种加载程序下复现仪器化模具的轴向应力、侧压比和等静压截距律，单单元、125单元与千单元网格结果一致；将同一卡不经修改用于独立工业零件，可复现截面生坯密度，且上部与外部下冲头载荷较基准七家实现更接近测量，但仍低于测量；该残差为基准各实现共有，并由输入不确定性界定。4. 推导了辨识常数向两种内置卡的闭式映射；映射到*MAT_025是对分段曲面的等价重写，而映射到*MAT_159仅保持剪切支、静水截距和交线，连续帽盖在二者之间明显高于母椭圆，故文献DPC标定不宜直接转移到连续帽盖卡。5. 只要剪切支采用非关联流动，保留分段子午线的角点无需顶点处理；帽盖椭圆在交线处体积流动消失，与锥面零体积流动相容。6. 两种内置帽盖卡在材料点评估中以相同方式失败，原因均为常数冻结：冻结常数误差随密度变号，连续帽盖侧压比偏低，冻结椭圆趋势错误，压–剪–释放循环后保留峰值模具应力十分之一至三分之一的残余剪应力。7. 帽盖模型存在剪切致密化压力阈值：致密化仅赋给帽盖，且帽盖体积流动在其顶点消失，因此低于帽盖–锥面交线压力时，无论累积剪切多大均不产生致密化；高于该压力时，可达密度受交线硬化至施加围压的状态限制。该阈值源于模型构造而非特定标定，并界定了闭模标定模型可合法应用的范围。
