周期性手征蜂窝结构材料等效弹性模量预测方法技术

本发明专利技术公开了一种周期性手征蜂窝结构材料等效弹性模量预测方法，用于解决现有手征蜂窝结构材料等效弹性模量计算方法实用性差的技术问题。技术方案是取手征性蜂窝结构一个单胞为研究对象，把单胞简化为一系列梁和刚性节圆组成的结构，模拟结构单轴拉伸加载，运用梁的大变形理论计算单胞每个梁的变形，再列出结构整体力和力矩的平衡方程、结构的变形要求方程组，用牛顿迭代法求解方程组可得到结构受力及变形状态，进而得到结构的等效弹性模量。此方法基于结构大变形条件，可以预测周期性多孔材料在发生大变形下的非线性等效等弹性模量。

A method for predicting the equivalent elastic modulus of periodic chiral honeycomb structures

The invention discloses a prediction method for the equivalent elastic modulus of a periodic chiral honeycomb structure material, which is used to solve the technical problem of the practical calculation method of the equivalent elastic modulus of the existing honeycomb structure material. The technical scheme is take the chirality of the honeycomb structure of a single cell as the research object, the single cell is simplified as a series of circular beams and structures, structural simulation of uniaxial tensile loading, deformation calculation of single cell of each beam using large deformation theory of beam deformation equations, equilibrium equations, and then lists the structure requirements the whole structure of force and moment, with Newton iterative method for solving equations can be obtained by structure stress and deformation, and then obtain the equivalent elastic modulus of the structure. Based on the large deformation conditions, the method can predict the nonlinear equivalent elastic modulus of periodic porous material under large deformation.

【技术实现步骤摘要】
周期性手征蜂窝结构材料等效弹性模量预测方法
本专利技术涉及一种手征蜂窝结构材料等效弹性模量计算方法，特别涉及一种周期性手征蜂窝结构材料等效弹性模量预测方法。
技术介绍
应用于自适应机翼设计中周期性手征蜂窝结构材料韧带长度与厚度比大，具有柔性特征，在分析计算柔性手征蜂窝结构的等效弹性常数时需考虑韧带的几何非线性变形。文献“A.Spadoni,M.Ruzzene.Elasto-staticmicropolarbehaviorofachiralauxeticlattice[J].JournaloftheMechanicsandPhysicsofSolids,60(2012),p156-171”公开了一种六韧带手性蜂窝结构的等效弹性矩阵的计算方法。该方法把结构简化为梁计算出结构的应变能，运用微极理论推导出等效本构方程，进而得出结构的各弹性常数。但是，此方法是基于微极理论，只适用于线弹性结构，无法反应结构在大变形情况下的非线性等效弹性模量等。
技术实现思路
为了克服现有手征蜂窝结构材料等效弹性模量计算方法实用性差的不足，本专利技术提供一种周期性手征蜂窝结构材料等效弹性模量预测方法。该方法取手征性蜂窝结构一个单胞为研究对象，把单胞简化为一系列梁和刚性节圆组成的结构，模拟结构单轴拉伸加载，运用梁的大变形理论计算单胞每个梁的变形，再列出结构整体力和力矩的平衡方程、结构的变形要求方程组，用牛顿迭代法求解方程组可得到结构受力及变形状态，进而得到结构的等效弹性模量。此方法基于结构大变形条件，可以预测周期性多孔材料在发生大变形下的非线性等效等弹性模量。本专利技术解决其技术问题所采用的技术方案：一种周期性手征蜂窝结构材料等效弹性模量预测方法，其特点是包括以下步骤：步骤一、手征蜂窝简化为半径为r的圆环及与其相切的长为L的韧带组成，韧带的有效梁变形部分长度为l：l＝L-2l1其中α＝0.9。式中l1表示韧带与节圆相切部分长度，α为系数。步骤二、取手征性蜂窝结构单胞，由圆环和六条长为L/2的韧带组成。圆环的圆心距为R，圆环半径为r，圆环和韧带的壁厚为t，结构垂直于纸面的厚度为b，韧带与圆心连线夹角为β，材料本身杨氏模量Es＝1.6GPa，泊松比v＝0.3。步骤三、对周期性手征蜂窝结构进行单轴加载模拟，根据结构的变形机理，假设单胞第i条韧带的受力状态：受力大小Fi，力的方向角αi，并设为未知量，i取值范围为1、2、3。步骤四、根据结构的变形机理，设计未知量第i条韧带的转角γi，计算每条韧带的变形，得到每条韧带变形后的形状及位置；(a)AB边变形分析。A0B0为变形前位置，曲线AB为变形后状态，A’B为B点处切线。为节圆转角，α1为A点受力F1与A0B0的夹角，γ为AB上任一点处的转角。根据梁的弹性弯曲理论，AB的变形微分方程为：式中，I表示梁的惯性矩，s表示梁上任一点的位置。化简得到：式中，C为常数。已知当时，γ＝γ1转角达到最大，将这一条件代入上式得将C代入上式整理得：化简得微分表达式：式中，S＝s/l(0≤S≤0.5)，分别表示无量纲位置和无量纲力。引入中间量η，令得：两边积分得到无量纲力表达式：式中，参数F(m1)表示第一类完全椭圆积分，F(η1,m1)表示第一类不完全椭圆积分。然后求AB的投影长度。建立局部坐标系XBY，以力F1方向为X轴，其垂直方向为Y轴。沿力方向及垂直方向投影无量纲长度分别为：并进一步推导计算，得：式中，E(m1)表示第二类完全椭圆积分，E(η1,m1)表示第二类不完全椭圆积分。AB在全局坐标系下的投影长度为：x1＝[X1cos(α1-β)+Y1cos(α1-β-π/2)]·ly1＝[X1sin(α1-β)+Y1sin(α1-β-π/2)]·lA、B点坐标为：xA＝xB+x1yA＝yB+y1(b)CD边变形分析。C0D0为变形前位置，曲线CD为变形后状态，C’D为D点处切线。为节圆转角，α2为C点受力F2与初始位置C0D0的夹角，γ为CD上任一点处的转角。坐标系xOy为全局坐标系。以力F2方向为X轴，垂直方向为Y轴，建立局部坐标系XDY。根据梁的弹性弯曲理论，CD的变形微分方程为：化简得到无量纲力表达式：式中，参数求CD在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为：CD在全局坐标系下的投影长度为：x2＝[X2cos(α2-β+π/3)+Y2cos(α2-β-π/6)]·ly2＝[X2sin(α2-β+π/3)+Y2sin(α2-β-π/6)]·l求C、D点坐标为：xC＝xD+x2yC＝yD+y2(c)EF边变形分析。E0F0为变形前位置，曲线EF为变形后状态，E’F为F点处切线。为节圆转角，α3为E点受力F3与初始位置E0F0的夹角，γ为EF上任一点处的转角。坐标系xOy为全局坐标系。以F为原点，以力F3方向为X轴，垂直方向为Y轴，建立局部坐标系XFY。根据梁的弹性弯曲理论，EF的变形微分方程为：化简得无量纲力表达式：式中，参数求EF在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为：EF在全局坐标系下的投影长度为：x3＝[X3cos(α3-β-π/3)+Y3cos(α3-β-5π/6)]·ly3＝[X3sin(α3-β-π/3)+Y3sin(α3-β-5π/6)]·l求E、F点坐标为：xE＝xF+x3yE＝yF+y3步骤五、根据单胞处于平衡状态及单胞变形为中心对称形式，列出力Fi及韧带位置满足的平衡方程组，其未知量为αi、γi。结构变形后C、E两点关于横向轴线上下对称，变形后仍然满足：xC-xE＝0yC+yE＝0变形前横向节圆连线水平，变形后已然满足，即：yA＝0变形前C、E点横坐标为A点横坐标的1/2，变形后仍然满足，即：xA-2xE＝0横向单轴拉伸时，手性蜂窝结构横向排与排之间的等效作用力为零，得：F3sin(α3-β-π/3)-F2sin(α2-β+π/3)＝0单轴拉伸时，假设单胞上中心对称的两条韧带受力大小相等方向相反，因此单胞自然满足静力平衡，同时单胞满足力矩平衡，即：2F1xAsinω1+2F2(xCsinω2-yCcosω2)+2F3(xEsinω3-yEcosω3)＝0式中，定义参数为：ω1＝α1-βω2＝α2-β+π/3ω3＝α3-β-π/3根据前文中的分析，以上两式改写为：对于手性蜂窝单轴拉伸时，给定一个转角值，6个基本变量：α1、γ1、α2、γ2、α3、γ3；并联立以上六个方程组成非线性方程组。步骤六、利用牛顿迭代法求解非线性方程组，得出αi、γi各量；非线性方程组通过Newton-Raphson迭代法求得。为了保证迭代的收敛性，为分步加载。在迭代过程中，第一步迭代各变量的初值取结构线性变形结果的近似值，第k步迭代的初值取第k-1步的结果。至此，得到各个载荷步下的6个基本变量的解。步骤七、计算结构等效应力、等效应变，得到等效杨氏模量。求得变形后A点的横坐标后，A点的位移即是手性蜂窝结构的变形，等效应变为：单轴拉伸下AB、CD、EF三条韧带合力在横向投影即为单胞的等效外力，等效应力为：改写为：最后得等效拉伸杨氏模量为：由此得到六韧带手征蜂窝结构等效杨氏模量计算表达式。本专利技术的有益效果是：该方法取手征性蜂窝结构一个单胞为研究对象，把单胞简化为一系列梁和刚性节圆本文档来自技高网...

【技术保护点】
一种周期性手征蜂窝结构材料等效弹性模量预测方法，其特征在于包括以下步骤：步骤一、手征蜂窝简化为半径为r的圆环及与其相切的长为L的韧带组成，韧带的有效梁变形部分长度为l：l＝L‑2l1其中

【技术特征摘要】
1.一种周期性手征蜂窝结构材料等效弹性模量预测方法，其特征在于包括以下步骤：步骤一、手征蜂窝简化为半径为r的圆环及与其相切的长为L的韧带组成，韧带的有效梁变形部分长度为l：l＝L-2l1其中α＝0.9；式中l1表示韧带与节圆相切部分长度，α为系数；步骤二、取手征性蜂窝结构单胞，由圆环和六条长为L/2的韧带组成；圆环的圆心距为R，圆环半径为r，圆环和韧带的壁厚为t，结构垂直于纸面的厚度为b，韧带与圆心连线夹角为β，材料本身杨氏模量Es＝1.6GPa，泊松比v＝0.3；步骤三、对周期性手征蜂窝结构进行单轴加载模拟，根据结构的变形机理，假设单胞第i条韧带的受力状态：受力大小Fi，力的方向角αi，并设为未知量，i取值范围为1、2、3；步骤四、根据结构的变形机理，设计未知量第i条韧带的转角γi，计算每条韧带的变形，得到每条韧带变形后的形状及位置；(a)AB边变形分析；A0B0为变形前位置，曲线AB为变形后状态，A’B为B点处切线；为节圆转角，α1为A点受力F1与A0B0的夹角，γ为AB上任一点处的转角；根据梁的弹性弯曲理论，AB的变形微分方程为：式中，I表示梁的惯性矩，s表示梁上任一点的位置；化简得到：式中，C为常数；已知当时，γ＝γ1转角达到最大，将这一条件代入上式得将C代入上式整理得：化简得微分表达式：式中，S＝s/l(0≤S≤0.5)，分别表示无量纲位置和无量纲力；引入中间量η，令得：两边积分得到无量纲力表达式：式中，参数F(m1)表示第一类完全椭圆积分，F(η1,m1)表示第一类不完全椭圆积分；然后求AB的投影长度；建立局部坐标系XBY，以力F1方向为X轴，其垂直方向为Y轴；沿力方向及垂直方向投影无量纲长度分别为：并进一步推导计算，得：式中，E(m1)表示第二类完全椭圆积分，E(η1,m1)表示第二类不完全椭圆积分；AB在全局坐标系下的投影长度为：x1＝[X1cos(α1-β)+Y1cos(α1-β-π/2)]·ly1＝[X1sin(α1-β)+Y1sin(α1-β-π/2)]·lA、B点坐标为：xA＝xB+x1yA＝yB+y1(b)CD边变形分析；C0D0为变形前位置，曲线CD为变形后状态，C’D为D点处切线；为节圆转角，α2为C点受力F2与初始位置C0D0的夹角，γ为CD上任一点处的转角；坐标系xOy为全局坐标系；以力F2方向为X轴，垂直方向为Y轴，建立局部坐标系XDY；根据梁的弹性弯曲理论，CD的变形微分方程为：化简得到无量纲力表达式：式中，参数求CD在局部坐标系下投影长度；沿力方向及垂直方向投影无量纲长度分别为：CD在全局坐标系下的投影长度为：x2＝[X2cos(α2-β+π/3)+Y2cos(α2-β-π/6)]·ly2＝[X2sin(α2-β+π/3)+Y2sin(α2-β-π/6)]·l求C、D点坐标为：xC＝xD+x2yC＝yD+y2(c)EF边变形分析；E0F0为变形前位置，曲线EF为变形后状态，E’F为F点处切线；为节圆转角，α3为E点受力F3与初始位置E0F0的夹角，γ为EF上任一点处的转角；坐标系xOy为全局坐标系；以F为原点，以力F3方向为X轴，垂直方向为Y轴，建立局部坐标系XFY；根据梁的弹性弯曲理论，EF的变形微分方程为：化简得无量纲力表达式：

【专利技术属性】
技术研发人员：邱克鹏，王智，张卫红，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西,61