用于飞机铝合金结构性能预测的状态近场动力学方法
【专利摘要】一种力学数值模拟【技术领域】的用于预测飞机铝合金结构性能的状态近场动力学方法，首先将铝合金结构离散成一系列空间域内的物质点；其次根据近场动力学状态线性化理论推导出节点刚度矩阵，并在此基础上改进Kilic针对近场动力学键理论提出动力松弛法，求解计算所需参量；然后将铝合金结构总的位移载荷分为若干增量步，对每个增量步采用改进的动力松弛法迭代计算；结合相对和绝对准则改进不平衡力准则，选用该准则判定系统的平衡与否，若平衡则进入下一增量步，直至加载加速，否则进入下一轮迭代，最终得到铝合金结构的性能数据。通过本发明能够实现运用近场动力学状态理论的常规材料模型求解复杂平衡路径准静态问题。
【专利说明】用于飞机铝合金结构性能预测的状态近场动力学方法
【技术领域】
[0001]本发明涉及的是一种力学数值模拟【技术领域】的方法，具体是一种基于近场动力学状态理论常规材料模型的准静态数值模拟方法。
【背景技术】
[0002]近场动力学(Peridynamics,简称PD)是近期发展起来的一种多尺度力学方法，能有效解决宏、微观力学领域中的不连续问题。它的优势源于理论本身的特点:它将物体离散成一系列空间域内的物质点，基于非局部作用思想建立物质点间的作用关系，通过空间积分方程描述物质点的运动，可采用统一的模型和求解体系描述原子尺度到宏观尺度的力学行为，不连续现象自然产生，突破了连续性假设和空间微分方程在不连续问题上出现的求解瓶颈。ro方法兴起于2000年，目前存在两个理论分支一近场动力学键理论和近场动力学状态理论，其中:近场动力学状态理论突破了近场动力学键理论对传统材料建模的局限，更具先进性。按物质点间作用力方向的不同，可将近场动力学状态理论中的材料模型分为两类:常规材料和非常规材料，根据Seleson P, Parks M L.“On the role ofthe influence function in the peridynamic theory，，([J].Journal for MultiscaleComputational Engineering, 2011, 9(6):689 - 706.)，其中:常规材料模型与近场动力学键理论中的材料模型在一定条件下具有一致性。近场动力学状态理论提出于2007年，目前关于常规材料模型准静态问题模拟的数值方法研究很少，Mitchell J A.在“Anonlocal, ordinary, state - based plasticity model for peridynamics，，( [R].SandiaNational Laboratory Report, SAND2011 - 3166, Albuquerque, New Mexico, 2011.)中只提到了一种类似于混合法的数值求解思想，但没有给出具体实施方法及应用实例，且该方法对一些特殊问题的求解会失效，如极值型失稳问题。Kilic B, Madenci E.在“An adaptivedynamic relaxtion method for quas1- static simulation using the peridynamictheory”([J].Theoretical and Applied Fracture Mechanics.2010, 53:194 -204.)中提出的动力松弛法能解决此类问题，但该方法仅限于近场动力学键理论。
【发明内容】
本发明针对现有技术存在的上述不足，提出一种用于飞机铝合金结构性能预测的状态近场动力学方法，通过本发明能够实现运用近场动力学状态理论的常规材料模型求解复杂平衡路径准静态问题。
本发明是通过以下技术方案实现的，本发明首先将结构离散成一系列空间域内的物质点；其次根据近场动力学状态线性化理论推导出节点刚度矩阵，并在此基础上改进Kilic针对近场动力学键理论提出动力松弛法，求解计算所需参量；然后将结构总的位移载荷分为若干增量步，对每个增量步采用改进的动力松弛法迭代计算，增量步开始时根据节点的空间位置选取初始位移场，使结构的空间构型接近新的平衡状态以提高算法的收敛速度；结合相对和绝对准则改进不平衡力准则，选用该准则判定系统的平衡与否，若平衡则进入下一增量步，直至加载结束，否则进入下一轮迭代；最终得到铝合金结构的性能数据。
[0003]本发明具体包括以下步骤:
[0004]步骤一、将结构离散成一系列空间域内的物质点:将物体离散成均匀的正方体晶格，节点(积分计算点)置于正方体晶格的中心。
[0005]步骤一中，当采用薄板铝合金结构时，将其离散成板状晶格以减少节点数量提高计算效率。
[0006]所述的薄板铝合金结构中的板状晶格在其他方向的尺寸一致，以使节点键族分布均匀。
[0007]所述的正方体晶格在各方向的晶格数量不少于2，以避免节点刚度矩阵和某方向相关的元素都为O导致的计算失败。
[0008]步骤二、根据近场动力学状态线性化理论推导出节点刚度矩阵，具体步骤包括:
[0009]设结构离散后，物质点X附近有Nx个节点qi (1≤i≤Nx),满足| q「x |〈2 δ，根据PD状态线性化理论，导出节点X的增量平衡方程为:
【权利要求】
1.一种用于预测飞机铝合金结构性能的状态近场动力学方法，其特征在于，首先将铝合金结构离散成一系列空间域内的物质点；其次根据近场动力学状态线性化理论推导出节点刚度矩阵，并在此基础上改进K1l1C针对近场动力学键理论提出动力松弛法，求解计算所需参量；然后将铝合金结构总的位移载荷分为若干增量步，对每个增量步采用改进的动力松弛法迭代计算；结合相对和绝对准则改进不平衡力准则，选用该准则判定系统的平衡与否，若平衡则进入下一增量步，直至加载加速，否则进入下一轮迭代，最终得到铝合金结构的性能数据。
2.根据权利要求1所述的方法，其特征是，所述方法具体包括以下步骤: 步骤一、将物体离散成均匀的正方体晶格，节点即积分计算点置于正方体晶格的中心； 步骤二、物质点X附近有Nx个节点q1 (1≤1≤Nx)，满足1 Q1-X 1〈2 δ，根据PD状态线性化理论，导出节点X的增量平衡方程为:
3.根据权利要求1或2所述的方法，其特征是，当采用薄板铝合金结构时，将其离散成板状晶格以减少节点数量提高计算效率。
4.根据权利要求3所述的方法，其特征是，所述的薄板铝合金结构中的正方体晶格在其他方向的尺寸一致，以使节点键族分布均匀。
5.根据权利要求2所述的方法，其特征是，所述的正方体晶格在各方向的晶格数量不少于2，以避免节点刚度矩阵和某方向相关的元素都为O导致的计算失败。
6.根据权利要求1或2所述的方法，其特征是，所述的节点在不同自由度上可采用不同的相对或绝对平衡容差。
【文档编号】G06F17/50GK103942381SQ201410150670
【公开日】2014年7月23日 申请日期:2014年4月15日 优先权日:2014年4月15日
【发明者】孙璐妍, 余音, 汪海 申请人:上海交通大学