1.一种用于结构变形分析的微势近场动力学方法,其特征在于:所述该方法包括以下步骤:步骤1:建立实体模型,确定相应材料属性;
步骤2:将结构离散成一系列空间域内的物质点;
步骤3:采用键长变化的二阶变形度量物质点的变形;
步骤4:结构离散的物质点间采用微键势能函数;
步骤5:引入物质点间的微键势能函数,得到的微键势能函数在物质点近场范围积分来获得该物质点的非局部弹性应变能密度;
步骤6:通过非局部弹性应变能密度的变分推导物质点间的键力模型;
步骤7:设置计算初始条件,施加力或位移边界条件;
步骤8:提交静动力计算;
步骤9:由变形结果开展应变分析。
2.根据权利要求1所述的微势近场动力学方法,其特征在于:所述步骤3中采用键长变化的二阶变形度量物质点的变形,二阶变形度量为:
2 2
Δ=|Y<ξ〉|‑|ξ|,
其中式中ξ是参考构型下键的长度,Y<ξ〉是当前构型下键的变形状态,变形状态与未变形的键长有关,遵循Cauchy‑Born法则Y<ξ>=Fξ,其中F是变形梯度;
物质点键长的二阶变形度量可以分解为Δ=Δv+Δd, 和Δd=2eijξiξj其中下标i和j表示物质点编号, 和
3.根据权利要求1所述的微势近场动力学方法,其特征在于:所述步骤4中结构离散的物质点间采用微键势能函数,微键势能函数为:其中Δv表示体积变形,Δd表示切向变形,ω1<ξ>和ω〈2ξ〉表示影响函数,ω3<ξ>=ω1<ξ>+ω2<ξ>。
4.根据权利要求1所述的微势近场动力学方法,其特征在于:所述步骤5中引入物质点间的微键势能函数,得到的微键势能函数在物质点近场范围积分来获得该物质点的非局部弹性应变能密度,非局部弹性应变能密度为:其中H表示物质点的近场,Vξ表示物质点的体积。
5.根据权利要求1所述的微势近场动力学方法,其特征在于:所述步骤6中通过非局部弹性应变能密度的变分推导物质点间的键力模型,键力模型为:δW=ξH[(2ω1<ξ>Δv+ω3<ξ)Δd)δΔv+(2ω2<ξ)Δd+ω3<ξ)Δv)δΔd]dVξ,经过交换积分变量可以得到微键势函数的键力模型其中F1=(ω〈1β〉‑ω〈2β〉)Δ,F2=ω2<β>Δ+ω1<β>Δv+ω2<β>Δd,A为物质点变形前的面积,Sβ为物质点变形后的面积,影响函数具体为其中ν为材料的泊松比,δ为物质点
的近场半径尺度。
6.根据权利要求1所述的微势近场动力学方法,其特征在于:所述步骤8中提交静动力计算通过更新粒子位移,更新粒子位移为:其中n表示计算步数,n=0表示初始时刻,c为阻尼系数 b为体积力,Λ为对角密度矩阵,取值为 ξmin为物质点间的最小键长,Vδ为物质点近场范围体积,α≥1为安全系数,取a=2,等效刚度矩阵K由物质点间的相互作用力对相对位移求偏导得到:其中,Vj为粒子编号为j的体积,更新粒子的位移为