成员变更粒子子系统的构型动量与应力

《International Journal of Engineering Science》：Configurational momentum and stress in particle subsystems with changing membership

　　

利用跳变-输运恒等式推导了开放粒子子系统的平衡方程。当标准力非零时，为构造内禀做功需要引入构型力。研究恢复了构型动量与Eshelby张量的形式。

论文解读

研究背景与问题

构型力学(configurational mechanics)关注与材料结构变化功共轭(work-conjugate)的力，区别于与标准运动对应的标准力(standard force)。Eshelby在弹性奇异性、晶格缺陷和能量-动量张量(energy–momentum tensor)研究中奠定了构型力学基础；后续研究将构型力视为具有自身平衡与做功的基本相互作用；Anderson等已在不依赖自然材料构型的空间连续统理论中得到构型动量(configurational momentum)和Eshelby关系。然而，对于粒子系统中成员通过接纳(admission)或删除(deletion)变化的开放粒子子系统(open particle subsystem)，尚缺少从对力(pair force)框架出发的统一推导。若仅用标准力描述成员变化，做功的内禀性(intrinsicity)会强制成员加权标准力为零，说明标准力不足以表征材料结构变化，因此需要引入构型力。该研究旨在建立开放粒子子系统的构型动量与构型应力理论。

研究内容与主要结论

研究人员以零-壹成员函数(zero–one membership function)指定子系统成员，以成员跳变表示接纳与删除。由于成员值只在有限粒子位置给定，无法唯一确定二维分离面，因此引入辅助光滑标量成员场(auxiliary smooth scalar membership field)Θ，在不决定实际成员归属的前提下提供局部几何、力与能量权重以及成员变化率。利用跳变-输运恒等式(jump-transport identity)将封闭子系统(closed subsystem)的质量、线动量、角动量、动能平衡推广到开放子系统；对构型动量单独假设封闭子系统平衡后再推广。通过要求成员转移做功(membership-transfer working)对辅助水平集(level set)观测速度的切向变化不变，并比较独立的成员率与局部法向变化，得到逐粒子构型动量关系。再通过机械能不平衡(mechanical energy imbalance)和对相互作用能识别成员张力，最终构造标准应力与构型应力，恢复空间Eshelby关系的形式。该研究为含成员变化的粒子系统提供了无歧义且内禀的构型力学框架，并澄清了离散边际能量与连续统自由能之间需额外对应。

关键研究方法

研究方法主要包括：一是跳变-输运恒等式，将固定成员期间量的变化与成员跳变时的转移量分离；二是辅助光滑成员场，构造水平集运动学并利用切向速度自由性检验做功内禀性；三是标架无差异性(frame indifference)导出对力作用-反作用与共线条件及封闭子系统平衡；四是在做功中引入逐粒子构型力与构型动量转移；五是通过机械能不平衡和非耗散/耗散成员转移确定成员张力，耗散情形采用对称半正定迁移率矩阵(mobility matrix)；六是沿连接粒子位置的线段积分负成员场梯度构造标准应力，并通过对辅助梯度局部变分表征构型应力。研究为纯理论推导，无样本队列。

研究结果

开放子系统的标准平衡

将跳变-输运恒等式与封闭子系统的质量、线动量、角动量、动能平衡结合，得到开放子系统的对应积分平衡。成员跳变处仅由被接纳或删除粒子所携带的量贡献，且假设不存在集中于跳变的额外标准冲量、角冲量或能量供给。

构型动量平衡

引入比构型动量(specific configurational momentum)、构型接触力(configurational contact force)合力以及内部和外部构型体力(configurational body force)合力，对封闭子系统假设构型动量平衡，再用同一跳变-输运恒等式得到开放子系统构型动量平衡。内部构型体力合力在线动量平衡中没有对应项。

做功内禀性与逐粒子构型关系

对只含标准力贡献的做功施加切向不变性，并比较独立选取的局部法向，得到所有成员加权标准力必须为零。加入逐粒子构型力和构型动量转移后，组合力沿局部法向，并通过成员率和法向的独立变化得到逐粒子关系u_i=?v_i以及预Eshelby关系(pre-Eshelby relation)f_i^Θ+c_i^Θ=π_i^Θh_i，其中π_i^Θ为成员张力(membership tension)。

成员张力的能量识别

引入加权对相互作用能(weighted pair-interaction energy)Ψ，其成员导数?Ψ/?θ_i为边际对相互作用能(marginal pair-interaction energy)。机械能不平衡被分解为固定成员值下的内部耗散(internal dissipation)与成员转移贡献。非耗散成员转移要求转移贡献对独立成员率消失，得到π_i^Θ=?Ψ/?θ_i?k_i，即边际对相互作用能减去粒子动能。耗散成员转移则由对称半正定迁移率矩阵θ?_i=Σ_jM_ijφ_j给出，保证成员转移耗散非负。

应力构造与连续统恢复

利用作用-反作用和对力共线条件，沿连接粒子位置的线段积分得到标准力合力的精确表示及余项误差界。赋予代表体积(representative volume)ΔV_i后构造对称标准应力T_i。在非耗散成员转移下，通过局部改变辅助梯度表征构型应力C_i=(?Ψ/?θ_i?k_i)ΔV_i^?11?T_i，并证明其对称。引入光滑连续场后，在粒子位置成立u=?v和C=ρ(ψ?|v|²/2)1?T，形式与空间Eshelby关系(spatial Eshelby relation)一致。

讨论与结论

讨论指出，构型力应作为原始量处理，成员跳变处不存在额外冲量或能量供给；标准力不足以描述材料结构变化，因此必须引入构型力。耗散成员转移需要状态相关的共轭系数与迁移率矩阵。研究结论可概括为：在空间表示中恢复了Anderson等的Eshelby关系形式，但当前ψ插值的是粒子的比边际对相互作用能，而非子系统单位被代表质量的相互作用能；将ψ解释为连续统比自由能(specific free energy)仍需额外建立离散与连续统能量之间的对应。此外，多体相互作用以及将成员赋予键、面或缺陷芯等结构特征，是进一步扩展该框架的方向。