一种蜂窝芯弯曲过程中芯格变形评估方法技术

本发明专利技术提出一种针对夹层结构的蜂窝芯弯曲过程中芯格变形评估方法，通过提取变形后蜂窝芯表面芯格轮廓并对芯格顶点进行识别，将每个芯格假设为由m个端点组成的m边形单元，通过对比m边形单元变形前后节点位置，确定每个节点位移；构建m边形单元形函数，由节点位移求得单元形函数对应系数进而确定单元对应的变形梯度；计算单元对应的应变张量，用单元对应的主应变及单元面积对芯格变形进行综合评估，实现了在复杂自由曲面外形夹层结构制造过程中判定芯格完整性。

【技术实现步骤摘要】
一种蜂窝芯弯曲过程中芯格变形评估方法
本专利技术涉及一种蜂窝芯弯曲过程中芯格变形评估方法，属于复合材料夹层结构制造领域。
技术介绍
蜂窝夹层结构是由两块高强度上、下面板和中间一层厚而轻的蜂窝芯夹层组成，具有轻质、高比强度、高比刚度及高能量吸收效率等优异的特性，广泛运用于航空航天等领域，如机翼前缘U型件、襟翼、副翼、扰流板方向舵、升降舵以及垂、平尾壁板、发动机短舱、整流罩等结构。由于蜂窝制造简单、质量较轻同时具有良好的韧性和抗冲击能力，因而在航空夹层结构中运用比例最大。然而，采用蜂窝夹层的复合材料构件大多具有复杂自由曲面外形，在蜂窝芯与型面贴合过程中，芯格由于弯曲产生严重的拉伸、挤压变形，使得芯格发生坍塌，严重影响蜂窝芯的力学特性。如何快速准确对弯曲后蜂窝芯格变形程度进行定量化评估是保证蜂窝芯格完整性及夹层结构优良力学特性的前提。
技术实现思路
为解决现有技术存在的问题，实现在复杂自由曲面外形夹层结构制造过程中判定芯格完整性，本专利技术提出一种蜂窝芯弯曲过程中芯格变形评估方法，其基本原理为：对弯曲后蜂窝芯试件非贴模表面芯格轮廓进行提取，获得蜂窝芯格变形后二维轮廓图，遍历单个芯格轮廓，提取芯格m个端点坐标，对比芯格变形前后端点位置，计算芯格各个端点对应的位移；利用芯格m个端点构造一个平面m边形单元并构建用于描述m边形单元位移场的二次多项式函数；由单元位移场及单元变形前后节点坐标关系计算变形梯度进而求得单元对应的应变张量，用单元对应的主应变及单元面积对芯格变性进行综合评估。本专利技术的技术方案为：所述一种蜂窝芯弯曲过程中芯格变形评估方法，包括以下步骤：步骤1：对弯曲后蜂窝芯芯格轮廓进行提取，得到芯格变形后二维轮廓图；对芯格变形后二维轮廓图进行处理后提取芯格轮廓边界线li，i＝1,2,3,…,n，n为芯格数量；遍历芯格轮廓线li，得到芯格m条边界线两两相交的m个端点，提取各个端点在预先设定的全局笛卡尔坐标系下对应的坐标(xi,j,yi,j)，j＝1,2,3,…,m；设蜂窝芯格变形前为标准的正m边形，以标号为1的芯格端点为原点建立局部笛卡尔坐标系，在局部坐标系下，变形前芯格端点坐标为(Xi,j,Yi,j)，j＝1,2,3,…,m，变形后端点坐标为(x′i,j,y′i,j)，通过坐标变换得到；根据局部坐标系下的变形前后芯格端点坐标计算各个端点对应的相对位移其中ui,j为第i个芯格端点j在局部坐标系x轴方向的相对位移，vi,j为第i个芯格中端点j在y轴方向的相对位移，且作为原点的端点相对位移值为0；步骤2：利用包含[1XYX2XY]项的二次不完全多项式表征芯格平面m边形位移场[uv]＝C[ab]，其中u、v分别为芯格在x、y方向的相对位移场，a、b分别为x、y方向的位移场对应的二次多项式系数矩阵，C为芯格端点坐标矩阵：分别求得x、y方向的位移场对应的二次多项式系数矩阵a＝C-1u，b＝C-1v；步骤3：对于蜂窝芯某一芯格，变形后的各个位置坐标是变形前坐标的函数，坐标关系满足x＝X+D，其中根据得到变形梯度F；进一步根据变形梯度F得到Green应变张量E在弯曲曲率方向和垂直弯曲曲率方向的值，以及得到Cauchy-Green应变张量中表征单元面积变化的第三不变量I3来对芯格变型进行综合评估；其中Green应变张量I为单位张量，C＝FTF为Cauchy-Green应变张量，I3＝det(C)。进一步的，步骤1中对芯格变形后二维轮廓图进行处理的过程为：去噪、增强和二值化处理。进一步的，步骤1中芯格m个端点按照某一时针顺序进行排序。有益效果本专利技术中提出的蜂窝芯弯曲过程中芯格变形评估方法通过提取变形后蜂窝芯表面芯格轮廓并对芯格顶点进行识别，将每个芯格假设为由m个端点组成的m边形单元，通过对比m边形单元变形前后节点位置，确定每个节点位移；构建m边形单元形函数，由节点位移求得单元形函数对应系数进而确定单元对应的变形梯度；计算单元对应的应变张量，用单元对应的主应变及单元面积对芯格变形进行综合评估；能够快速准确对弯曲后蜂窝芯格变形程度进行定量化评估，为保证蜂窝芯格完整性及夹层结构优良力学特性提供了前提。本专利技术的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本专利技术的实践了解到。附图说明本专利技术的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：图1弯曲后蜂窝芯试件图2提取后的芯格表面形状图3变形后芯格端点提取流程图4局部坐标系下芯格示意图图5相对坐标系下变形前后芯格示意图图6变形后蜂窝芯格面积变化示意图图7变形后蜂窝芯格应变统计图具体实施方式下面详细描述本专利技术的实施例，所述实施例是示例性的，旨在用于解释本专利技术，而不能理解为对本专利技术的限制。本实施例中以图2所示的六边形蜂窝芯(m＝6，因此芯格边界及端点数均为6)弯曲后表面芯格形状轮廓图为例，结合附图，说明一种蜂窝芯弯曲过程中芯格变形评估方法的具体过程。步骤1：对弯曲后蜂窝芯的芯格形状轮廓图进行数学形态学处理实现图像的去噪及增强，将图片转换为只有0/1两种像素值的二值图；根据蜂窝芯格边界像素的邻接关系对芯格轮廓线进行提取得到芯格轮廓线li，i＝1,2,3,…,n，n为芯格数量。遍历芯格轮廓线li所有像素点，根据芯格端点为芯格6条边界线两两相交的几何关系，得到芯格6个端点对应的坐标(xi,j,yi,j)(j＝1,2,3,…,6)，局部蜂窝芯格端点坐标提取流程如图3所示。以芯格最底端端点(y值最小，即min(yi,j)对应的端点)作为第一个点并按顺时针顺序对芯格的6个端点进行排序，并以第一个点为原点建立一个局部笛卡尔坐标系o′x′y′，如图4所示，在该坐标系下对芯格6个端点对应的全局坐标(xi,j,yi,j)按进行坐标变换得到局部坐标系下6个端点对应局部坐标(x′i,j,y′i,j)(j＝1,2,3,…,6)，其中端点1的局部坐标为(0,0)。假设蜂窝芯格变形前均为标准的正六边形，在该局部坐标系下对比变形前后芯格形状(变形前芯格形状为123456，变形后芯格形状为1′2′3′4′5′6′)，如图5所示，芯格变形前端点局部坐标为(Xi,j,Yi,j)(j＝1,2,3,…,6)。根据局部坐标系下芯格变形前后端点坐标，计算变形后芯格各个端点对应的相对位移(中ui,j为第i个芯格端点j在局部坐标系x轴方向的相对位移，vi,j为第i个芯格中端点j在y轴方向的相对位移，且作为原点的端点相对位移值为0；)，其中步骤2：将蜂窝芯格6个端点作为节点构造平面六边形单元，以构建的平面六边形单元的变形来表征芯格的变形，变形后平面六边形单元6个节点位移即为步骤1中芯格6个端点对应的位移其中节点1的位移值恒为0。利用包含[1XYX2XY]项的二次不完全多项式表征平面六边形位移场[uv]＝C[ab]，u、本文档来自技高网...

【技术保护点】
1.一种蜂窝芯弯曲过程中芯格变形评估方法，其特征在于，包括以下步骤：/n步骤1：对弯曲后蜂窝芯芯格轮廓进行提取，得到芯格变形后二维轮廓图；对芯格变形后二维轮廓图进行处理后提取芯格轮廓边界线l

【技术特征摘要】
1.一种蜂窝芯弯曲过程中芯格变形评估方法，其特征在于，包括以下步骤：
步骤1：对弯曲后蜂窝芯芯格轮廓进行提取，得到芯格变形后二维轮廓图；对芯格变形后二维轮廓图进行处理后提取芯格轮廓边界线li，i＝1,2,3,…,n，n为芯格数量；
遍历芯格轮廓线li，得到芯格m条边界线两两相交的m个端点，提取各个端点在预先设定的全局笛卡尔坐标系下对应的坐标(xi,j,yi,j)，j＝1,2,3,…,m；
设蜂窝芯格变形前为标准的正m边形，以标号为1的芯格端点为原点建立局部笛卡尔坐标系，在局部坐标系下，变形前芯格端点坐标为(Xi,j,Yi,j)，j＝1,2,3,…,m，变形后端点坐标为(x′i,j,y′i,j)，通过坐标变换得到；
根据局部坐标系下的变形前后芯格端点坐标计算各个端点对应的相对位移其中ui,j为第i个芯格端点j在局部坐标系x轴方向的相对位移，vi,j为第i个芯格中端点j在y轴方向的相对位移，且作为原点的端点相对位移值为0；
步骤2：利用包含[1XYX2XY]项的二次不完全多项式表征芯格平面m边形位移场[uv]＝C[ab]，其中u、...

【专利技术属性】
技术研发人员：赵志勇，李玉军，刘闯，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61