[发明专利]一种用于结构变形分析的微势近场动力学方法

说明书

本发明公开了一种用于结构变形分析的微势近场动力学方法，该方法包括以下步骤：建立实体模型，确定相应材料属性；将结构离散成一系列空间域内的物质点；采用键长变化的二阶变形度量物质点的变形；结构离散的物质点间采用微键势能函数；引入物质点间的微键势能函数，得到的微键势能函数在物质点近场范围积分来获得该物质点的非局部弹性应变能密度；通过非局部弹性应变能密度的变分推导物质点间的键力模型；设置计算初始条件，施加力或位移边界条件；提交静动力计算；由变形结果开展应变分析。本发明通过引入物质点微键势能函数，推导得出全新的物质点键力模型，键力模型不包含任何未知数，不需要简化模型校对核键参数，简化步骤。

技术领域

本发明属于近场动力学本构推导和数值模拟技术领域，具体涉及一种用于结构变形分析的微势近场动力学方法。

背景技术

近场动力学(Peridynamics，简称PD)是近期发展起来的一种多尺度力学方法，能有效解决宏、微观力学领域中的不连续问题。近场的概念来自长程力作用，不同于局部理论即连续介质力学理论，只有在几何空间上相邻的物质点之间才会发生相互作用；非局部理论就像分子动力学，物质点之间能够在一定距离范围内发生相互作用，不需要在几何空间上接触，即物质点间存在长程力。目前存在两个理论分支——近场动力学键理论和近场动力学状态理论。

近场动力学状态理论提出于2007年，目前关于常规材料模型建立物质点间势函数推导键力的相关研究很少，Jincheng Fan在A micro-potential based Peridynamicmethod for deformation and fracturing in solids:A two-dimensionalformulation.([J].Computer Methods in Applied Mechanics and Engineering 360(2020):112751.)中将Xu-Needleman势函数引入，但该模型对键参数的确定需要通过简化获得。

中国专利数据库中(申请号为：201910679765.5，申请公布号：CN 110457790A)用于结构变形分析的近场动力学非连续伽辽金有限元方法，所述该方法包括如下步骤：步骤一，建立实体模型，确定各材料区域并赋予相应材料属性；步骤二，采用标准有限单元法的网格生成算法划分单元网格；步骤三，采用结点复制算法、生成匹配近场动力学的非连续伽辽金有限元这种新方法的非连续单元；步骤四，初始时刻计算、并在后续计算步中实时更新单元结点与高斯点坐标；步骤五，引入权函数构造键型近场动力学的弱形式方程，引入等参单元插值形函数以便用结点变量值表征任意点变量，根据高斯点受其近场域范围内的其它高斯点的非局部作用、以及高斯点间相对位移计算单元结点等效载荷，计算应力边界和体力的结点等效外载荷，计算结点质量矩阵，给出整个构型的受力变形计算方法；步骤六，设置初始条件，采用拉格朗日乘子法施加体力、应力和位移边界条件；步骤七，提交静动力计算；步骤八，由变形结果开展应力、应变分析。

对比上述专利，尚有以下不足：通过非连续伽辽金静动力有限元解法完善近场动力学的数值求解体，但是非连续伽辽金静动力有限元解法建立的模型参数运算步骤繁琐，且局限在原有的材料变形本构模型。

发明内容

本发明针对上述的不足之处提供一种不需要对模型参数进行简化的微势近场动力学方法。

本发明目的是这样实现的：一种用于结构变形分析的微势近场动力学方法，其特征在于：所述该方法包括以下步骤：

步骤1：建立实体模型，确定相应材料属性；

步骤2：将结构离散成一系列空间域内的物质点；

步骤3：采用键长变化的二阶变形度量物质点的变形；

步骤4：结构离散的物质点间采用微键势能函数；

步骤5：引入物质点间的微键势能函数，得到的微键势能函数在物质点近场范围积分来获得该物质点的非局部弹性应变能密度；