法律状态公告日
法律状态信息
法律状态
2019-08-13
授权
授权
2018-03-30
实质审查的生效 IPC(主分类):G06F17/50 申请日:20171107
实质审查的生效
2018-03-06
公开
公开
技术领域
本发明涉及一种手征蜂窝结构材料等效弹性模量计算方法,特别涉及一种周期性手征蜂窝结构材料等效弹性模量预测方法。
背景技术
应用于自适应机翼设计中周期性手征蜂窝结构材料韧带长度与厚度比大,具有柔性特征,在分析计算柔性手征蜂窝结构的等效弹性常数时需考虑韧带的几何非线性变形。
文献“A.Spadoni,M.Ruzzene.Elasto-static micropolar behavior of achiral auxetic lattice[J].Journal of the Mechanics and Physics of Solids,60(2012),p156-171”公开了一种六韧带手性蜂窝结构的等效弹性矩阵的计算方法。该方法把结构简化为梁计算出结构的应变能,运用微极理论推导出等效本构方程,进而得出结构的各弹性常数。但是,此方法是基于微极理论,只适用于线弹性结构,无法反应结构在大变形情况下的非线性等效弹性模量等。
发明内容
为了克服现有手征蜂窝结构材料等效弹性模量计算方法实用性差的不足,本发明提供一种周期性手征蜂窝结构材料等效弹性模量预测方法。该方法取手征性蜂窝结构一个单胞为研究对象,把单胞简化为一系列梁和刚性节圆组成的结构,模拟结构单轴拉伸加载,运用梁的大变形理论计算单胞每个梁的变形,再列出结构整体力和力矩的平衡方程、结构的变形要求方程组,用牛顿迭代法求解方程组可得到结构受力及变形状态,进而得到结构的等效弹性模量。此方法基于结构大变形条件,可以预测周期性多孔材料在发生大变形下的非线性等效等弹性模量。
本发明解决其技术问题所采用的技术方案:一种周期性手征蜂窝结构材料等效弹性模量预测方法,其特点是包括以下步骤:
步骤一、手征蜂窝简化为半径为r的圆环及与其相切的长为L的韧带组成,韧带的有效梁变形部分长度为l:
l=L-2l1
其中
步骤二、取手征性蜂窝结构单胞,由圆环和六条长为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点处切线。
根据梁的弹性弯曲理论,AB的变形微分方程为:
式中,I表示梁的惯性矩,s表示梁上任一点的位置。
化简得到:
式中,C为常数。已知当
化简得微分表达式:
式中,S=s/l(0≤S≤0.5),
引入中间量η,令
两边积分得到无量纲力表达式:
式中,参数
然后求AB的投影长度。建立局部坐标系XBY,以力F1方向为X轴,其垂直方向为Y轴。沿力方向及垂直方向投影无量纲长度分别为:
并进一步推导计算,得:
式中,E(m1)表示第二类完全椭圆积分,E(η1,m1)表示第二类不完全椭圆积分。
AB在全局坐标系下的投影长度为:
x1=[X1cos(α1-β)+Y1cos(α1-β-π/2)]·l
y1=[X1sin(α1-β)+Y1sin(α1-β-π/2)]·l
A、B点坐标为:
xA=xB+x1
yA=yB+y1
(b)CD边变形分析。
C0D0为变形前位置,曲线CD为变形后状态,C’D为D点处切线。
坐标系xOy为全局坐标系。以力F2方向为X轴,垂直方向为Y轴,建立局部坐标系XDY。
根据梁的弹性弯曲理论,CD的变形微分方程为:
化简得到无量纲力表达式:
式中,参数
求CD在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为:
CD在全局坐标系下的投影长度为:
x2=[X2cos(α2-β+π/3)+Y2cos(α2-β-π/6)]·l
y2=[X2sin(α2-β+π/3)+Y2sin(α2-β-π/6)]·l
求C、D点坐标为:
xC=xD+x2
yC=yD+y2
(c)EF边变形分析。
E0F0为变形前位置,曲线EF为变形后状态,E’F为F点处切线。
坐标系xOy为全局坐标系。以F为原点,以力F3方向为X轴,垂直方向为Y轴,建立局部坐标系XFY。
根据梁的弹性弯曲理论,EF的变形微分方程为:
化简得无量纲力表达式:
式中,参数
求EF在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为:
EF在全局坐标系下的投影长度为:
x3=[X3cos(α3-β-π/3)+Y3cos(α3-β-5π/6)]·l
y3=[X3sin(α3-β-π/3)+Y3sin(α3-β-5π/6)]·l
求E、F点坐标为:
xE=xF+x3
yE=yF+y3
步骤五、根据单胞处于平衡状态及单胞变形为中心对称形式,列出力Fi及韧带位置满足的平衡方程组,其未知量为αi、γi。
结构变形后C、E两点关于横向轴线上下对称,变形后仍然满足:
xC-xE=0
yC+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
根据前文中的分析,以上两式改写为:
对于手性蜂窝单轴拉伸时,给定一个转角
步骤六、利用牛顿迭代法求解非线性方程组,得出αi、γi各量;
非线性方程组通过Newton-Raphson迭代法求得。为了保证迭代的收敛性,
步骤七、计算结构等效应力、等效应变,得到等效杨氏模量。
求得变形后A点的横坐标后,A点的位移即是手性蜂窝结构的变形,等效应变为:
单轴拉伸下AB、CD、EF三条韧带合力在横向投影即为单胞的等效外力,等效应力为:
改写为:
最后得等效拉伸杨氏模量为:
由此得到六韧带手征蜂窝结构等效杨氏模量计算表达式。
本发明的有益效果是:该方法取手征性蜂窝结构一个单胞为研究对象,把单胞简化为一系列梁和刚性节圆组成的结构,模拟结构单轴拉伸加载,运用梁的大变形理论计算单胞每个梁的变形,再列出结构整体力和力矩的平衡方程、结构的变形要求方程组,用牛顿迭代法求解方程组可得到结构受力及变形状态,进而得到结构的等效弹性模量。此方法基于结构大变形条件,可以预测周期性多孔材料在发生大变形下的非线性等效等弹性模量。
考虑周期性手征蜂窝结构材料的柔性特征,采用本发明方法预测其等效弹性模量时考虑了韧带的几何非线性变形,在手征性蜂窝结构等效应变0~0.2范围内,计算出的等效应力应变曲线与仿真结果很吻合,等效弹性模量与仿真结果误差在10%以内,预测模型在更接近实际变形情况,预测结果更可靠,计算精度更高。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明周期性手征蜂窝结构材料等效弹性模量预测方法中手征蜂窝示意图。
图2是本发明方法实施例中手征蜂窝单胞受力图。
图3是本发明方法实施例中手征蜂窝AB边受拉变形图。
图4是本发明方法实施例中手征蜂窝CD边受拉变形图。
图5是本发明方法实施例中手征蜂窝EF边受拉变形图。
具体实施方式
参照图1-5。本发明周期性手征蜂窝结构材料等效弹性模量预测方法具体步骤如下:
步骤一、手征蜂窝可简化为半径为r的圆环及与其相切的长为L的韧带组成,韧带的有效梁变形部分长度为l:
l=L-2l1
其中
步骤二、取手征性蜂窝结构单胞,由圆环和六条长为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点处切线。
根据梁的弹性弯曲理论,AB的变形微分方程为:
式中,I表示梁的惯性矩,s表示梁上任一点的位置。
化简得到:
式中,C为一常数。已知当
化简得微分表达式:
式中,S=s/l(0≤S≤0.5),
引入中间量η,令
两边积分可得到无量纲力表达式:
式中,参数
然后求AB的投影长度。建立局部坐标系XBY,以力F1方向为X轴,其垂直方向为Y轴。沿力方向及垂直方向投影无量纲长度分别为:
并进一步推导计算,可得:
式中,E(m1)表示第二类完全椭圆积分,E(η1,m1)表示第二类不完全椭圆积分。
AB在全局坐标系下的投影长度为:
x1=[X1cos(α1-β)+Y1cos(α1-β-π/2)]·l
y1=[X1sin(α1-β)+Y1sin(α1-β-π/2)]·l
A、B点坐标为:
xA=xB+x1
yA=yB+y1
(b)CD边变形分析。
C0D0为变形前位置,曲线CD为变形后状态,C’D为D点处切线。
坐标系xOy为全局坐标系。以力F2方向为X轴,垂直方向为Y轴,建立局部坐标系XDY。
根据梁的弹性弯曲理论,CD的变形微分方程为:
化简可得到无量纲力表达式:
式中,参数
求CD在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为:
CD在全局坐标系下的投影长度为:
x2=[X2cos(α2-β+π/3)+Y2cos(α2-β-π/6)]·l
y2=[X2sin(α2-β+π/3)+Y2sin(α2-β-π/6)]·l
求C、D点坐标为:
xC=xD+x2
yC=yD+y2
(c)EF边变形分析。
E0F0为变形前位置,曲线EF为变形后状态,E’F为F点处切线。
坐标系xOy为全局坐标系。以F为原点,以力F3方向为X轴,垂直方向为Y轴,建立局部坐标系XFY。
根据梁的弹性弯曲理论,EF的变形微分方程为:
化简可得无量纲力表达式:
式中,参数
求EF在局部坐标系下投影长度。沿力方向及垂直方向投影无量纲长度分别为:
EF在全局坐标系下的投影长度为:
x3=[X3cos(α3-β-π/3)+Y3cos(α3-β-5π/6)]·l
y3=[X3sin(α3-β-π/3)+Y3sin(α3-β-5π/6)]·l
求E、F点坐标为:
xE=xF+x3
yE=yF+y3
步骤五、根据单胞处于平衡状态及单胞变形为中心对称形式,可列出力Fi及韧带位置满足的平衡方程组,其未知量为αi、γi。
结构变形后C、E两点关于横向轴线上下对称,变形后仍然满足:
xC-xE=0
yC+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
根据前文中的分析,以上两式可改写为:
对于手性蜂窝单轴拉伸时,给定一个转角
步骤六、利用牛顿迭代法求解非线性方程组,得出αi、γi各量;
非线性方程组通过Newton-Raphson迭代法求得。为了保证迭代的收敛性,
步骤七、计算结构等效应力、等效应变,得到等效杨氏模量。
求得变形后A点的横坐标后,A点的位移即是手性蜂窝结构的变形,等效应变为:
单轴拉伸下AB、CD、EF三条韧带合力在横向投影即为单胞的等效外力,等效应力为:
可改写为:
最后可得等效拉伸杨氏模量为:
由此得到六韧带手征蜂窝结构等效杨氏模量计算表达式,该式考虑了韧带的几何非线性变形,体现了周期性手征手蜂窝结构材料的柔性特征。
机译: 广义等效等效转换厚度和等效弹性模量的多层基底沉降或沉降/应力分析方法
机译: 广义等效折算厚度等效弹性模量的多层地基桩基沉陷应力分析方法
机译: 木材组合材料弯曲强度极限和弹性模量的预测方法