复合材料参数多阶段识别方法技术

本发明专利技术公开了一种复合材料确定性参数多阶段识别方法，包括以下步骤：步骤1，在有限元分析软件中建立复合材料的初始有限元模型，并计算其静力学位移响应；步骤2，根据复合材料结构静力学位移响应与复合材料弹性参数之间的力学关系，计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度；步骤3，在步骤1和2的基础上，根据相对灵敏度分析结果使用多阶段方法识别复合材料参数；步骤4，根据试验数据识别复合材料全部弹性参数。本发明专利技术使用三点弯位移数据，多阶段模型修正方法识别复合材料全部弹性参数。多阶段识别方法同时解决了识别灵敏度数量级相差较大的参数，导致修正过程中灵敏度矩阵病态的问题。

【技术实现步骤摘要】
复合材料参数多阶段识别方法
本专利技术涉及复合材料结构动力学参数识别
，具体涉及一种采用基于静力学位移响应来识别复合材料全部弹性参数的多阶段识别方法。
技术介绍
复合材料结构具有高比强度、高比刚度、耐疲劳等优异性能，广泛应用航空航天领域。但由于复合材料结构的不确定性，其材料参数也会随着生产工艺，内部结构不同而发生改变。准确的复合材料的弹性参数是研究复合材料力学性能的基础。复合材料等效参数为正交各向异性，具有9个独立的参数。传统的测量复合材料参数方法通过拉伸，剪切测量应变。再通过应力应变关系计算其刚度。传统的模型修正方法是通过计算参数对响应的灵敏度，选取灵敏度较大的参数进行识别，放弃了其他的参数。或同时识别所有的参数，对识别结果保留一定的误差。显然这样建立的复合材料模型不会准确，对于结构后续实际工况研究也加入了一定的误差。难以捕捉其实际载荷下的力学性能。
技术实现思路
本专利技术的目的在于提供一种复合材料确定性参数多阶段识别方法，以解决现有技术存在的复合材料仿真模拟的精度差的问题。为实现上述目的，本专利技术采用的技术方案为：一种复合材料确定性参数多阶段识别方法，包括以下步骤：步骤1，在有限元分析软件中建立复合材料的初始有限元模型，并通过有限元分析软件计算其静力学位移响应；步骤2，根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系，计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度；步骤3，在步骤1和2的基础上，研究复合材料参数多阶段预测方法，将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集，基于计算静力学位移参数与实验静力学参数的残差建立目标函数，并对待修正参数进行迭代优化，优先识别灵敏度大的参数集合，直至精度达到要求且所有的参数集合均收敛；并进行仿真验证研究；步骤4，根据试验数据识别复合材料全部弹性参数。所述步骤1包括以下步骤：步骤1.1，在有限元软件中，根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型；其中，9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33，三个泊松比ν12、ν23、ν31，与三个平面内的剪切弹性模量G12、G23、G31；步骤1.2，根据步骤1.1得到的初始有限元模型，通过有限元分析软件对其进行静力学分析，得到复合材料的静力学位移响应。所述步骤2包括以下步骤：步骤2.1，由复合材料力学知识得到材料刚度矩阵C如下式表示：上式中，E1、E2、E3分别表示主轴方向1,2,3的弹性模量；νij表示单轴应力σ在j方向作用时引起i方向应变的泊松耦合系数，即泊松比；G12、G23、G31分别表示1-2,2-3,3-1平面的剪切弹性模量；刚度矩阵C中各元素与独立弹性参数之间的关系如下所示：C11＝(1-ν23ν32)·E11·B,C12＝(ν12+ν13ν32)·E22·BC22＝(1-ν13ν31)·E22·B,C23＝(ν23+ν21ν13)·E33·BC33＝(1-ν12ν21)·E33·B,C31＝(ν31+ν21ν32)·E11·B其中上式中，Ckl分别为6×6刚度矩阵中对应位置上的元素，其中，0≤k≤6，0≤l≤6；步骤2.2，研究复合材料结构静力学响应与复合材料弹性参数之间的力学关系，并进行相对灵敏度分析，计算公式为：其中，fe为静力试验参数，通过三点弯静力学试验测得；Δf为静力学位移响应与实验参数的残差；p为需要识别的弹性参数；S为相对灵敏度，各个弹性参数对于静力学响应的影响程度，计算公式如下表示：其中f抽象地代表静力学响应，fi为仿真模型有限元分析出的静力学位移响应，由于刚度矩阵C为对称矩阵故只需要求静力学相响应fi对刚度矩阵中元素C11、C12、C13、C22、C23、C333、C443、C553、C66的偏导，为刚度矩阵中各元素对待求的独立弹性参数的偏导数。所述步骤3包括以下步骤：步骤3.1，根据相对灵敏度分析结果，将相差小于十倍的参数灵敏度划分为相同的集合；则复合材料的9个独立弹性参数划分为n个子集；步骤3.2，构建以计算响应参数与实验响应参数的加权残差为目标函数；即在参数p的取值范围p1≤p≤p2内找到一个pa，其中p的取值范围为9个独立弹性参数可能值的上下限；使目标函数J(p)为试验与计算参数的加权残差取极小值；e为特征值的残差，λe,λa(p)分别为试验与计算的参数，加权矩阵W为反映各特征值残差相对权重的对角阵；步骤3.3，采用分布迭代修正的模式，对步骤3.1中n个集合分别识别，优先识别灵敏度大的参数集合，直至精度达到要求且所有的参数集合均收敛；从而识别复合材料全部弹性参数。所述步骤4的根据试验数据识别复合材料全部弹性参数包括以下步骤包括以下步骤：步骤4.1，开展试验验证研究，采用复合材料板结构三点弯试验位移数据，来识别复合材料的全部的弹性参数。有益效果：本专利技术解决了繁琐的传统复合材料参数测量方法，使用三点弯位移数据，多阶段模型修正方法识别复合材料全部弹性参数。多阶段识别方法同时解决了识别灵敏度数量级相差较大的参数，导致修正过程中灵敏度矩阵病态的问题。对于建立准确的复合材料动力学模型具有指导作用，具有工程应用价值。附图说明图1是复合材料确定性参数多阶段识别方法具体步骤的流程图；图2是编织复合材料板仿真模型；图3是三点弯节点位移对9个独立弹性参数的相对灵敏度绝对值柱状图；图4是多阶段迭代曲线图。具体实施方式下面结合附图对本专利技术做更进一步的解释。如图1所示，本专利技术的一种复合材料确定性参数多阶段识别方法，包括以下步骤：一种复合材料确定性参数多阶段识别方法，其特征在于：包括以下步骤：步骤1，在有限元分析软件中建立复合材料的初始有限元模型，并通过有限元分析软件计算其静力学位移响应；具体步骤为：步骤1.1，在有限元软件中，根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型；其中，9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33，三个泊松比ν12、ν23、ν31，与三个平面内的剪切弹性模量G12、G23、G31；步骤1.2，根据步骤1.1得到的初始有限元模型，通过有限元分析软件对其进行静力学分析，得到复合材料的静力学位移响应。步骤2，根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系，计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度；具体步骤为：步骤2.1，由复合材料力学知识得到材料刚度矩阵C如下式表示：上式中，E1、E2、E3分别表示主轴方向1,2,3的弹性模量；νij表示单轴应力σ在j方向作用时引起i方向应变的泊松耦合系数，即泊松比；G12、G23、G31分别表示1-2,2-3,3-1平面的剪切弹性模量；刚度矩阵C中各元素与独立弹性参数之间的关系如下所示：C11＝(1-ν23ν32)·E11·B,C12＝(ν12+ν13ν32)·E22·BC22＝(1-ν13ν31)·E22·B,C23＝(ν23+ν21ν13)·E33·BC33＝(1-ν12ν21)·E33·B,C31＝(ν31+ν21ν32)·E11·B其中上式中，Ckl分别为6×6刚度矩阵中对应位置上的元素，其本文档来自技高网...

【技术保护点】
1.一种复合材料确定性参数多阶段识别方法，其特征在于：包括以下步骤：步骤1，在有限元分析软件中建立复合材料的初始有限元模型，并通过有限元分析软件计算其静力学位移响应；步骤2，根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系，计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度；步骤3，在步骤1和2的基础上，研究复合材料参数多阶段预测方法，将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集，基于计算静力学位移参数与实验静力学参数的残差建立目标函数，并对待修正参数进行迭代优化，优先识别灵敏度大的参数集合，直至精度达到要求且所有的参数集合均收敛；并进行仿真验证研究；步骤4，根据试验数据识别复合材料全部弹性参数。

【技术特征摘要】
1.一种复合材料确定性参数多阶段识别方法，其特征在于：包括以下步骤：步骤1，在有限元分析软件中建立复合材料的初始有限元模型，并通过有限元分析软件计算其静力学位移响应；步骤2，根据复合材料的结构静力学位移响应与复合材料的独立弹性参数之间的力学关系，计算复合材料的独立弹性参数对于静力学位移响应的相对灵敏度；步骤3，在步骤1和2的基础上，研究复合材料参数多阶段预测方法，将复合材料全部材料参数集合根据相对灵敏度大小划分成n个子集，基于计算静力学位移参数与实验静力学参数的残差建立目标函数，并对待修正参数进行迭代优化，优先识别灵敏度大的参数集合，直至精度达到要求且所有的参数集合均收敛；并进行仿真验证研究；步骤4，根据试验数据识别复合材料全部弹性参数。2.根据权利要求1所述的复合材料确定性参数多阶段识别方法，其特征在于：所述步骤1包括以下步骤：步骤1.1，在有限元软件中，根据复合材料的结构的几何尺寸、正交各向异性9个独立弹性参数、边界条件及载荷建立复合材料的初始有限元模型；其中，9个独立弹性参数包括三个材料主轴方向的弹性模量E11、E22、E33，三个泊松比ν12、ν23、ν31，与三个平面内的剪切弹性模量G12、G23、G31；步骤1.2，根据步骤1.1得到的初始有限元模型，通过有限元分析软件对其进行静力学分析，得到复合材料的静力学位移响应。3.根据权利要求1所述的复合材料确定性参数多阶段识别方法，其特征在于：所述步骤2包括以下步骤：步骤2.1，由复合材料力学知识得到材料刚度矩阵C如下式表示：上式中，E1、E2、E3分别表示主轴方向1,2,3的弹性模量；νij表示单轴应力σ在j方向作用时引起i方向应变的泊松耦合系数，即泊松比；G12、G23、G31分别表示1-2,2-3,3-1平面的剪切弹性模量；刚度矩阵C中各元素与独立弹性参数之间的关系如下所示：C11＝(1-ν23ν32)·E11·B,C12＝(ν12+v13ν32)·E22·BC22＝(1-v13ν31)·E22·B,C23＝(ν23+ν21...

【专利技术属性】
技术研发人员：姜东，朱东辉，费庆国，曹芝腑，陈素芳，朱锐，
申请(专利权)人：南京林业大学，
类型：发明
国别省市：江苏,32