CN115196036B - 一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 - Google Patents
一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 Download PDFInfo
- Publication number
- CN115196036B CN115196036B CN202210876282.6A CN202210876282A CN115196036B CN 115196036 B CN115196036 B CN 115196036B CN 202210876282 A CN202210876282 A CN 202210876282A CN 115196036 B CN115196036 B CN 115196036B
- Authority
- CN
- China
- Prior art keywords
- poisson ratio
- honeycomb
- positive
- deformation
- negative poisson
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000013461 design Methods 0.000 title claims abstract description 25
- 229910001285 shape-memory alloy Inorganic materials 0.000 claims abstract description 17
- 210000004027 cell Anatomy 0.000 claims description 113
- 210000002421 cell wall Anatomy 0.000 claims description 41
- 238000006073 displacement reaction Methods 0.000 claims description 39
- 238000005452 bending Methods 0.000 claims description 34
- 210000002287 horizontal cell Anatomy 0.000 claims description 28
- 238000004088 simulation Methods 0.000 claims description 22
- 238000012937 correction Methods 0.000 claims description 13
- 241000186321 Cellulomonas Species 0.000 claims description 10
- 230000011218 segmentation Effects 0.000 claims description 6
- 230000001413 cellular effect Effects 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims description 3
- 238000013519 translation Methods 0.000 claims description 3
- 238000003466 welding Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 abstract description 7
- 238000004364 calculation method Methods 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 229920002160 Celluloid Polymers 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 239000002520 smart material Substances 0.000 description 1
- 229920002994 synthetic fiber Polymers 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C3/00—Wings
- B64C3/26—Construction, shape, or attachment of separate skins, e.g. panels
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C3/00—Wings
- B64C3/28—Leading or trailing edges attached to primary structures, e.g. forming fixed slots
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C3/00—Wings
- B64C3/38—Adjustment of complete wings or parts thereof
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Mechanical Engineering (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明提出一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法,首先通过智能材料形状记忆合金驱动正、负泊松比蜂窝结构变形;并且通过对正、负泊松比结构的单胞设计以及相互之间的约束使得两个结构的变形之间能发生相互作用,从而实现整个蜂窝夹芯的变形,达到前缘变厚度的目的,减少了机械结构的使用;而且本发明结合机翼前缘变厚度需求,给出的结构设计机理清楚,通过计算能够得到准确的驱动参数,继而能够精确控制机翼前缘变形。
Description
技术领域
本发明属于柔性蒙皮结构技术领域,具体为一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法。
背景技术
智能变体飞行器是指一类能够实时感知外部环境,并会根据飞行任务、状态和环境等因素实时、自主地改变本体外形布局,以实现在当前影响因素综合下的最优性能飞行的新概念飞行器。机翼作为航空器的传统主要部件,其在智能变体飞行器的研究中占重要的地位。
机翼变厚度是机翼变形的一种重要方式。传统变厚度机翼主要由刚体部分,柔性蒙皮和安装在机翼内部的驱动器构成,驱动器普遍采用的是液压或电机驱动器。于是机翼变厚度的方式为驱动机械结构对机翼表面的柔性蒙皮施加载荷,通过柔性蒙皮的被动变形来改变机翼的厚度。传统机翼变厚度的方法虽然有回复力大和便于控制的优点,但会使得机翼结构重量增加和机翼内部结构设计变得复杂。随着智能材料的发展,智能材料逐渐被选择作为驱动器代替传统的液压和电机驱动器,减少了变厚度机翼的重量,但这类实现机翼变厚度的方法仍要借助机械结构,机械结构带来的问题依然存在。
将智能材料引入柔性蒙皮的设计中,使得柔性蒙皮不借助机械结构实现主动变形,是解决机械结构所带来的问题的一种方案。在这之前,有学者将智能材料形状记忆合金加入波纹式蒙皮设计中,实现了蒙皮主动地单向伸长和缩短。在此基础上,有学者将人工纤维肌肉代替形状记忆合金驱动器实现了同样的功能。但在前缘变厚度机翼中,前缘由钝变尖的过程中具有沿弦向变形大,以及尖锐前缘点附近的局部变形精度难以控制等难点,基于波纹式主动变形蒙皮对于这些问题并不能很好地解决。
发明内容
针对现有技术存在的问题,尤其是针对机翼前缘变厚度中出现的问题,本发明提出一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法。本发明通过智能材料形状记忆合金驱动正、负泊松比蜂窝结构变形;本发明通过对正、负泊松比结构的单胞设计以及相互之间的约束使得两个结构的变形之间能发生相互作用,从而实现整个蜂窝夹芯的变形,达到前缘变厚度的目的,减少了机械结构的使用;而且本发明结合机翼前缘变厚度需求,给出的结构设计机理清楚,通过计算能够得到准确的驱动参数,继而能够精确控制机翼前缘变形。
本发明的技术方案为:
一种基于正、负泊松比蜂窝夹芯主动变形的柔性蒙皮设计方法,包括以下步骤:
步骤1:将机翼初始翼型以及目标翼型各自的上翼型和下翼型均沿弦长等分为N段,计算每一段翼型的弯曲角度,要求每一段翼型的弯曲角度均不能大于设定值,如果超过设定值,则增加分段数;
步骤2:根据步骤1确定的分段数N,设计柔性蒙皮中蜂窝夹芯的蜂窝单胞结构参数:
所述蜂窝夹芯中的蜂窝单胞包括负泊松比蜂窝单胞和正泊松比蜂窝单胞;
首先确定负泊松比蜂窝单胞的结构参数:
根据步骤1确定的分段数N,以及已知的机翼弦长L,确定每一分段的弦长为L/N,则负泊松比蜂窝单胞x向的长度为L/N,并在0.5L/N<h1<L/N范围内,设定负泊松比蜂窝单胞水平胞壁的长h1,则得关系式
h1-l1 cosθ1=0.5L/N (1)
再设定负泊松比蜂窝单胞的y向等效泊松比为νyx1,得关系式:
联立(1)式和(2)式得到负泊松比蜂窝单胞的倾斜胞壁长度l1和倾斜胞壁的角度θ1,并设定负泊松比蜂窝单胞胞壁的厚度为t1,高度为b1;从而确定负泊松比蜂窝单胞的结构参数;
其次根据负泊松比蜂窝单胞的结构参数确定正泊松比蜂窝单胞的结构参数:
正泊松比蜂窝单胞x向的长度为L/N,在0<h2<0.5L/N范围内,设定正泊松比蜂窝单胞水平胞壁的长h2,得关系式
h2+l2 cosθ2=0.5L/N (3)
再由正、负泊松比蜂窝单胞在y方向上的长度相等,得关系式:
2l2 sinθ2=2l1 sinθ1 (4)
联立(3)式和(4)式求得正泊松比蜂窝单胞的倾斜胞壁长度l2和倾斜胞壁的角度θ2,并设定正泊松比蜂窝单胞胞壁的厚度为t2,高度为b2;从而确定正泊松比蜂窝单胞的结构参数;
步骤3:将步骤2设计的正、负泊松比蜂窝单胞组合成变形单元;所述变形单元由正泊松比蜂窝多胞、负泊松比蜂窝多胞和两块形状记忆合金驱动器组合而成;
其中正泊松比蜂窝多胞为正泊松比蜂窝单胞沿y方向排列组合而成的纵向多胞;负泊松比蜂窝多胞为负泊松比蜂窝单胞沿y方向排列组合而成的纵向多胞;且纵向多胞中的单胞数量n由机翼展长M除以蜂窝单胞y向的长度所得结果取整得到;
正泊松比蜂窝多胞与负泊松比蜂窝多胞沿z方向重合叠加在一起,且变形单元中每一个正泊松比蜂窝单胞的一侧中间水平胞壁端部和与其对应重合的负泊松比蜂窝单胞的同一侧中间水平胞壁端部焊接固定;两块形状记忆合金驱动器布置在正、负泊松比蜂窝多胞组合体的正面和背面,沿y向放置且位于中间位置,分别和正、负泊松比蜂窝多胞连接,多胞中每个和驱动器接触的水平胞壁的接触部分都和驱动器焊接在一起;
每个变形单元对应一个翼型分段,将N个分段对应的变形单元沿弦向依次焊接固定,就得到柔性蒙皮中的蜂窝夹芯。
进一步的,步骤1中,所述设定值为15°。
进一步的,步骤1中,对于初始翼型以及目标翼型各自的上翼型和下翼型,分别进行以下操作来计算每一段翼型弯曲角度:
首先拟合翼型曲线,之后通过曲线平移和旋转,使得翼型的后缘处于坐标系原点并在该点和坐标系x轴相切;
其次将翼型沿弦长等分,每一段翼型对应弦长相等,根据拟合的翼型曲线计算每一段翼型两端点的斜率T1和T2,T2>T1,根据公式
得到该分段的弯曲角度β。
进一步的,负泊松比蜂窝单胞以及正泊松比蜂窝单胞的胞壁厚度和高度不大于4mm,且胞壁厚度小于高度。
进一步的,泊松比蜂窝单胞胞壁厚度t2等于负泊松比蜂窝单胞胞壁厚度t1;正泊松比蜂窝单胞高度b2等于负泊松比蜂窝单胞高度b1。
进一步的,正泊松比蜂窝多胞以及负泊松比蜂窝多胞的y方向顶部和底部的水平胞壁的厚度为对应蜂窝单胞胞壁厚度的两倍。
一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮,所述柔性蒙皮具有蜂窝夹芯,所述蜂窝夹芯由上述任一方法设计得到。
一种柔性蒙皮中蜂窝夹芯主动变形的控制方法,根据公式
计算蜂窝夹芯上各个变形单元的驱动位移,并通过形状记忆合金驱动器产生驱动位移,将蜂窝夹芯从初始翼型形状变形为目标翼型形状;其中δyj为第j个翼型分段对应的变形单元的驱动位移,T是正、负泊松比蜂窝多胞的高度和,μ是修正系数,αj为机翼初始翼型第j个分段和目标翼型第j个分段的弯曲角度之差,αj=αj目标-αj初始。
进一步的,修正系数μ有限元数值分析进行数值模拟实验得到:
首先建立变形单元模型,在驱动位移小于或等于1.5mm的前提下至少选取5个不同的驱动位移δy施加给变形单元,每个驱动位移构成一组数值模拟实验,共进行s组数值模拟实验,每组实验的驱动位移δyi使得变形单元产生绕y轴的弯曲角度αi;
然后根据公式
计算得到修正系数μ。
有益效果
本发明提出了一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法。首先通过智能材料形状记忆合金驱动正、负泊松比蜂窝结构变形;并且通过对正、负泊松比结构的单胞设计以及相互之间的约束使得两个结构的变形之间能发生相互作用,从而实现整个蜂窝夹芯的变形,达到前缘变厚度的目的,减少了机械结构的使用;而且本发明结合机翼前缘变厚度需求,给出的结构设计机理清楚,通过计算能够得到准确的驱动参数,继而能够精确控制机翼前缘变形。
本发明的附加方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为圆弧两端点的斜率和圆心角的关系;
图2为负、正泊松比蜂窝单胞结构与负泊松比蜂窝纵向多胞;(a)为负泊松比蜂窝单胞结构,(b)为正泊松比蜂窝单胞结构,(c)为负泊松比蜂窝纵向多胞结构;
图3为变形单元结构;
图4为变形单元的边界条件和变形单元受载弯曲的状态;
图5为柔性蒙皮蜂窝夹芯;
图6为负、正泊松比蜂窝板;
图7为本发明对应的流程图;
图8为已知机翼的二维翼型;其中1为初始翼型,2为目标翼型;
图9为分段处理后的已知机翼二维翼型;
图10为蜂窝夹芯的初始状态数值模型;
图11为数值模拟结果和目标翼型曲线的对比图。
具体实施方式
本发明提出一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法,能够应用在变厚度机翼上,尤其能够解决前缘变厚度机翼的控制难题。
下面以前缘变厚度机翼为例,已知该前缘变厚度机翼的二维翼型弦长为L,机翼展长为M,且已知其初始状态和变形后目标状态的二维翼型,分别简称为初始翼型和目标翼型。由于本发明设计的是主动变形柔性蒙皮中的蜂窝夹芯,所以根据前缘变厚度机翼已知的初始状态和目标状态而确定的蜂窝夹芯初始状态和目标状态也是已知的。
基于上述已知参数,下面给出基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,主要设计的是柔性蒙皮中的蜂窝夹芯。
步骤1:将已知机翼的初始翼型以及目标翼型各自的上翼型和下翼型均沿弦长等分为N段,计算每一段翼型的弯曲角度,要求每一段的弯曲角度均不能大于设定值,如果超过设定值,则增加分段数。
这里对初始翼型以及目标翼型的具体处理过程相同,下面以初始翼型为例进行说明:
在数学模型中拟合出初始翼型中上翼型和下翼型的曲线,然后通过曲线平移和旋转,使得翼型的后缘为原点并在该点和x轴相切。
将初始翼型中上翼型和下翼型沿弦长等分,每一段的弦长是相等的,等分后的每一分段都可认为是光滑的圆弧,圆弧两端点为分段点,分段点的斜率通过拟合的曲线获得。假设某一分段两端的斜率分别为t1和t2(t2>t1),则通过几何关系(如图1),可获得该分段的弯曲角度,即圆弧的圆心角
对于分段数的要求是各个分段中最大的弯曲角度要小于或等于15°,若分段中最大的弯曲角度大于15°,则应该增加分段数。假设二维翼型被等分为N段,从后缘的分段开始到前缘的分段依次将每个分段编号为1~N。
步骤2:根据步骤1确定的分段数N,设计蜂窝夹芯中的蜂窝单胞的结构参数。
所述蜂窝夹芯中的蜂窝单胞包括负泊松比蜂窝单胞和正泊松比蜂窝单胞。负、正泊松比蜂窝单胞的结构如图2(a)和(b)所示,其中h1和h2分别为负泊松比蜂窝单胞以及正泊松比蜂窝单胞的水平胞壁的长、l1和l2分别为负泊松比蜂窝单胞以及正泊松比蜂窝单胞的倾斜胞壁的长、t1和t2分别为负泊松比蜂窝单胞以及正泊松比蜂窝单胞的胞壁厚度、θ1和θ2分别为负泊松比蜂窝单胞以及正泊松比蜂窝单胞的倾斜胞壁角度,b1和b2分别为负泊松比蜂窝单胞以及正泊松比蜂窝单胞的高度,即单胞沿z向的长度。
首先确定负泊松比蜂窝单胞的结构参数。
由步骤1可知,每个分段的弦长为L/N,所以这里确定负泊松比蜂窝单胞x向的长度为L/N,其中蜂窝单胞x向的长度为单胞左端和右端之间的连线长,取负泊松比蜂窝单胞水平胞壁的长为h1(0.5L/N<h1<L/N),如图2(a)可得关系式:
h1-l1 cosθ1=0.5L/N (1)
再取负泊松比蜂窝单胞的y向等效泊松比为νyx1,可得关系式:
联立(1)式和(2)式便可求得负泊松比蜂窝单胞的倾斜胞壁长度l1和倾斜胞壁的角度θ1。最后设定负泊松比蜂窝单胞胞壁的厚度为t1,高度为b1,其中取胞壁的厚度和高度要小于或等于4mm,并且t1≤b1。这样就得到了负泊松比蜂窝单胞的结构参数。
其次通过负泊松比蜂窝单胞的结构参数确定正泊松比蜂窝单胞的结构参数。同样,正泊松比蜂窝单胞x向的长度为L/N,取正比松比蜂窝单胞的水平胞壁的长为h2(0<h2<0.5L/N),可得关系式:
h2+l2 cosθ2=0.5L/N (3)
再由正、负泊松比蜂窝单胞在y方向上的长度相等,蜂窝单胞y向长度为单胞顶端和低端之间的连线长,可得关系式:
2l2 sinθ2=2l1 sinθ1 (4)
联立(3)式和(4)式便可求得正泊松比蜂窝单胞的倾斜胞壁长度l2和倾斜胞壁的角度θ2。最后设定正泊松比蜂窝单胞胞壁厚度为t2,高度为b2,这里取t2=t1,高度b2=b1。
于是,正、负泊松比蜂窝单胞的结构参数就已全部确定。
步骤3:将正、负泊松比蜂窝单胞组合成变形单元。
所述变形单元是本方法提出的概念,如图3所示,它由正泊松比蜂窝多胞、负泊松比蜂窝多胞和两块形状记忆合金驱动器组合而成。其中,负、正泊松比蜂窝多胞为沿y方向(展向)的纵向多胞。纵向多胞的结构如图2(c)所示,它由若干相同的蜂窝单胞沿y向排列组合而成,需要注意的是纵向多胞位于顶部和底部的水平胞壁的厚度为蜂窝单胞设计的胞壁厚度的两倍。
将已知机翼的展长M除以步骤2中蜂窝单胞y向的长度2(l1 sinθ1+t1)并取整,就得到了变形单元纵向多胞中的单胞数量n,并且整数n是最接近M/2(l1 sinθ1+t1)的整数。
如图3所示,变形单元中,正、负泊松比蜂窝多胞沿z方向重合叠加在一起,变形单元中每一个正泊松比蜂窝单胞的一侧中间水平胞壁端部和与其对应重合的负泊松比蜂窝单胞的同一侧水平胞壁端部焊接在一起,正、负泊松比蜂窝多胞就是通过这些地方连接在一起。两块形状记忆合金驱动器布置在正、负泊松比蜂窝多胞组合体的正面和背面,沿y向放置且位于中间位置,分别和正、负泊松比蜂窝多胞连接,多胞中每个和驱动器接触的水平胞壁的接触部分都和驱动器焊接在一起。
由此就设计得到了单个沿机翼展向的变形单元,每个变形单元对应一个翼型分段,将N个分段对应的变形单元沿弦向依次焊接固定,就得到柔性蒙皮中的蜂窝夹芯。
下面进一步对柔性蒙皮蜂窝夹芯的变形控制进行分析:
如图4(a)所示,在变形单元的左边,约束中间水平胞壁左端x、y和z向位移自由度,约束剩下的水平胞壁左端x和z向位移自由度。此时形状记忆合金驱动器产生y向的驱动位移δy时,变形单元便会绕y轴向正泊松比蜂窝多胞的一侧的弯曲(如图4(b)所示),变形单元的弯曲角度α和驱动位移δy的关系为:
其中μ是修正系数,它是一个和蜂窝单胞构型相关的量,若蜂窝单胞的构型确定,修正系数μ也就确定了。θ1和θ2分别是负、正泊松比蜂窝单胞倾斜胞壁的角度,n是变形单元多胞中单胞的个数,T是正、负泊松比多胞的高度和,即T=b1+b2。
这里,修正系数μ可以借助有限元数值分析软件进行数值模拟实验来确定。在软件中建立变形单元的模型,选用步骤2中确定的蜂窝单胞的构型。数值模拟实验中可用位移载荷代替形状记忆合金驱动器的功能。在驱动位移小于或等于1.5mm的前提下至少选取5个不同的驱动位移δy施加给变形单元,每个驱动位移构成一组数值模拟实验。若进行了s组数值模拟实验,每组实验的驱动位移δyi都会使得变形单元产生绕y轴的弯曲角度αi。
αi获取的方法为:通过数值模拟实验,变形单元会近似弯曲为一段圆弧,获取圆弧上离散点的坐标,通过坐标在数学软件中模拟出数值模拟实验结果中的圆弧,然后用步骤1中圆弧的圆心角公式就可以获得变形单元的弯曲角度αi。
最后结合式(5)和s组数值模拟实验中获得的变形单元的弯曲角度αi,可确定修正系数的值为:
这样式(5)中除了α和δy,其他参数的值都已确定。接下来涉及到式(5)的计算中,μ的值是不变的。
接下来,利用步骤1确定的各个分段的弯曲角度,来计算驱动位移,从而实现驱动控制。
在步骤1中,将已知机翼的二维翼型,包括上翼型和下翼型,分别均沿弦长等分为N段,将每个分段都用变形单元代替,变形单元的形状和编号都与初始翼型的分段相同,这样N个变形单元就组合成了已知机翼的柔性蒙皮蜂窝夹芯,如图5所示。
图5所示的是由7个变形单元组成的已知机翼初始状态的柔性蒙皮蜂窝夹芯,方框中为前缘部分的变形单元。图6所示的是图5中的蜂窝芯板是由负泊松比蜂窝板和正泊松比蜂窝板沿z轴叠加在一起的。在每两个变形单元之间,正、负泊松比蜂窝板接触地方是焊接在一起的,具体如图5中变形单元5和6之间的深色部分。
给每个变形单元施加驱动位移,使其变形成为目标翼型相应分段的形状,这样整个蜂窝芯板就变形成已知机翼目标状态的形状。
这里结合式(5)可以得到蜂窝夹芯上各个变形单元的驱动位移为(上翼型或下翼型):
其中αj为已知机翼初始翼型第j个分段和目标翼型第j个分段弯曲角度之差,即αj=αj目标-αj初始。
所计算的驱动位移为蜂窝夹芯上相应编号的变形单元上的形状记忆合金驱动器的驱动位移,通过这些驱动位移就可以将蜂窝夹芯从初始翼型形状变形为目标翼型形状。
下面详细描述本发明的实施例,所述实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
在本发明的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“长度”、“宽度”、“厚度”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”、“顺时针”、“逆时针”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
本实施例是用本发明的方法将图8中的已知机翼从初始状态变形成为目标状态,也就是将蜂窝夹芯从初始状态变形为目标状态。图8中是已知机翼的二维翼型,1是初始翼型,2是目标翼型,已知机翼的弦长为450mm,选取的变形机翼部分的翼展为330mm。
由于这里采用的已知机翼是上下对称的,所以只选取图8中上半部分即上翼型来应用本方法,下表是节选自上翼型的初始翼型和目标翼型形状的坐标。
首先设计柔性蒙皮中的蜂窝夹芯:
步骤1,将上翼型按弦长等分,计算每一分段的弯曲角度。
在数学软件中输入上表中的坐标,用二次样条差值函数拟合出初始翼型和目标翼型的曲线,然后进行旋转和平移使得翼型的后缘点为原点并且和x轴相切,如图9所示,图9中上半部分为处理后的初始翼型,下半部分为处理后的目标翼型。
将二维翼型按弦长等分为6段,则每个分段的弦长为75mm,所以分段点的坐标依次为:0、-75、-150、-225、-300、-375和-450,从后缘的分段到前缘的分段依次将每个分段编号为1~6,如图9所示。依次计算分段点的斜率然后按照公式计算各个分段的弯曲角度,按照编号填入下表:
由表中数据可知,分段中最大的弯曲角度为9.79277°,符合分段的规则。
步骤2,设计蜂窝单胞的尺寸。
由步骤1可知,负泊松比蜂窝单胞x向的长度为75mm,取单胞水平胞壁的长为50mm,再取单胞y向等效泊松比为-1,则可得关系式:
则负泊松比蜂窝单胞倾斜胞壁的长l1为25mm,倾斜胞壁的角度θ1为60°,再取单胞的胞壁厚度t1为1mm,高b1为1mm。
同样,正泊松比蜂窝单胞x向的长度为75mm,取正比松比蜂窝单胞的水平胞壁的长为30mm,再由正、负泊松比蜂窝单胞在y方向上的长度相等可得关系式:
则单胞倾斜胞壁的长l2为22.91mm,倾斜胞壁的角度θ2为70.88°,由正、负泊松比蜂窝单胞的胞壁厚度和高度相等,有t2=t1=1mm,b2=b1=1mm。
正、负泊松比蜂窝单胞的结构参数就已全部确定,如下表所示。
t/mm | l/mm | h/mm | θ° | b/mm | vyx | |
正泊松比蜂窝单胞 | 1 | 22.91 | 30 | 70.88 | 1 | 1.66 |
负泊松比蜂窝单胞 | 1 | 25 | 50 | 60 | 1 | -1 |
步骤3,将正、负泊松比蜂窝单胞组合成变形单元。
将已知机翼的展长除以步骤2中蜂窝单胞的y向的长度并取整,可得:330÷[2×(25×sin60°+1)]=7.2846,取n=7,即变形单元中的正泊松比蜂窝多胞以及负泊松比蜂窝多胞分别都是由7个单胞组合而成。
变形单元中,正、负泊松比蜂窝多胞沿z方向重合叠加在一起,变形单元中每一个正泊松比蜂窝单胞的水平胞壁和与其对应重合的负泊松比蜂窝单胞的水平胞壁在端部焊接在一起,正、负泊松比蜂窝多胞就是通过这些地方连接在一起。两块形状记忆合金驱动器布置在正、负泊松比蜂窝多胞组合体的正面和背面,沿y向放置且位于中间位置,分别和正、负泊松比蜂窝多胞连接,多胞中每个和驱动器接触的水平胞壁的接触部分都和驱动器焊接在一起。
由此就设计得到了单个沿机翼展向的变形单元,每个变形单元对应一个翼型分段,将7个分段对应的变形单元依次焊接固定,就得到柔性蒙皮中的蜂窝夹芯。
接下来,确定柔性蒙皮中蜂窝夹芯的变形控制。
步骤4、确定修正系数μ:
在ABAQUS数值模拟软件中建立变形单元的数值模型,如图3所示,模型的结构参数为步骤2中设计的蜂窝单胞的结构参数,模型的材料参数为:弹性模量为68.9GPa,泊松比为0.33。对模型分别施加0.2mm、0.4mm、0.6mm、0.8mm和1.0mm的驱动位移,进行5组数值模拟实验,五组数值模拟实验下变形单元的弯曲角度如下表:
实验组 | 1 | 2 | 3 | 4 | 5 |
驱动位移/mm | 0.2 | 0.4 | 0.6 | 0.8 | 1.0 |
弯曲角度/° | 2.72 | 5.45 | 8.15 | 10.85 | 13.55 |
所以由修正系数的计算公式可计算出修正系数的值为:
步骤5:根据步骤1中得到的每个分段的弯曲角度,计算驱动位移。
在步骤1中,已知机翼的二维翼型被等分为了6段,将每个分段都用变形单元代替,变形单元的形状和编号都与初始翼型的分段相同,如图10所示为6个变形单元组成的已知机翼初始状态中柔性蒙皮蜂窝芯的数值模型,其中未加变形单元驱动器,驱动器的功能在数值模拟实验中用位移载荷代替。
给变形单元施加驱动位移,使其变形成为目标翼型相应分段的形状。根据驱动位移和弯曲角度的关系
将步骤1中计算的弯曲角度的差值代入上式,可得各个变形单元的驱动位移列入下表,注意公式中弯曲角度的单位为弧度。
变形单元 | 1 | 2 | 3 | 4 | 5 | 6 |
弯曲角度αj/° | -0.23027 | -0.68822 | -1.25225 | -2.07094 | -3.64665 | -9.17068 |
驱动位移/mm | -0.017 | -0.051 | -0.092 | -0.152 | -0.268 | -0.675 |
所计算的驱动位移为蜂窝夹芯上相应编号的变形单元驱动器的驱动位移,通过这些驱动位移就可以将蜂窝夹芯从初始翼型形状变形为目标翼型形状。
检验变形效果:
为了检验蜂窝芯板在计算的驱动位移下是否可以从已知机翼的初始状态形状变形成为目标状态形状,下面通过数值模拟,将变形后的蜂窝芯板的形状和已知机翼目标翼型之间进行对比。
已知机翼初始状态的数值模型如图10所示,将步骤5中计算的驱动位移施加到相应的变形单元上。注意,在数值模拟中施加驱动位移时,要先后设置两个分析步,第一个分析步为负泊松比蜂窝多胞的变形,第二个为正泊松比蜂窝多胞的变形。
如图11所示为数值模拟结果模型上离散的点和目标翼型曲线的对比图,可以看到数值模拟的结果和目标翼型有很好的一致性。
此外,我们也进行了实际的变形分析试验,试验结果同样表明,采用本发明设计的基于正、负泊松比蜂窝芯主动变形的柔性蒙皮,能够准确实现机翼前缘变厚度要求。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (9)
1.一种基于正、负泊松比蜂窝夹芯主动变形的柔性蒙皮设计方法,其特征在于:包括以下步骤:
步骤1:将机翼初始翼型以及目标翼型各自的上翼型和下翼型均沿弦长等分为N段,计算每一段翼型的弯曲角度,要求每一段翼型的弯曲角度均不能大于设定值,如果超过设定值,则增加分段数;
步骤2:根据步骤1确定的分段数N,设计柔性蒙皮中蜂窝夹芯的蜂窝单胞结构参数:
所述蜂窝夹芯中的蜂窝单胞包括负泊松比蜂窝单胞和正泊松比蜂窝单胞;
首先确定负泊松比蜂窝单胞的结构参数:
根据步骤1确定的分段数N,以及已知的机翼弦长L,确定每一分段的弦长为L/N,则负泊松比蜂窝单胞x向的长度为L/N,并在0.5L/N<h1<L/N范围内,设定负泊松比蜂窝单胞水平胞壁的长h1,则得关系式
h1-l1cosθ1=0.5L/N (1)
再设定负泊松比蜂窝单胞的y向等效泊松比为νyx1,得关系式:
联立(1)式和(2)式得到负泊松比蜂窝单胞的倾斜胞壁长度l1和倾斜胞壁的角度θ1,并设定负泊松比蜂窝单胞胞壁的厚度为t1,高度为b1;从而确定负泊松比蜂窝单胞的结构参数;
其次根据负泊松比蜂窝单胞的结构参数确定正泊松比蜂窝单胞的结构参数:
正泊松比蜂窝单胞x向的长度为L/N,在0<h2<0.5L/N范围内,设定正泊松比蜂窝单胞水平胞壁的长h2,得关系式
h2+l2cosθ2=0.5L/N (3)
再由正、负泊松比蜂窝单胞在y方向上的长度相等,得关系式:
2l2sinθ2=2l1sinθ1 (4)
联立(3)式和(4)式求得正泊松比蜂窝单胞的倾斜胞壁长度l2和倾斜胞壁的角度θ2,并设定正泊松比蜂窝单胞胞壁的厚度为t2,高度为b2;从而确定正泊松比蜂窝单胞的结构参数;
步骤3:将步骤2设计的正、负泊松比蜂窝单胞组合成变形单元;所述变形单元由正泊松比蜂窝多胞、负泊松比蜂窝多胞和两块形状记忆合金驱动器组合而成;
其中正泊松比蜂窝多胞为正泊松比蜂窝单胞沿y方向排列组合而成的纵向多胞;负泊松比蜂窝多胞为负泊松比蜂窝单胞沿y方向排列组合而成的纵向多胞;且纵向多胞中的单胞数量n由机翼展长M除以蜂窝单胞y向的长度所得结果取整得到;
正泊松比蜂窝多胞与负泊松比蜂窝多胞沿z方向重合叠加在一起,且变形单元中每一个正泊松比蜂窝单胞的一侧中间水平胞壁端部和与其对应重合的负泊松比蜂窝单胞的同一侧中间水平胞壁端部焊接固定;两块形状记忆合金驱动器布置在正、负泊松比蜂窝多胞组合体的正面和背面,沿y向放置且位于中间位置,分别和正、负泊松比蜂窝多胞连接,多胞中每个和驱动器接触的水平胞壁的接触部分都和驱动器焊接在一起;
每个变形单元对应一个翼型分段,将N个分段对应的变形单元沿弦向依次焊接固定,就得到柔性蒙皮中的蜂窝夹芯。
2.根据权利要求1所述一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,其特征在于:步骤1中,所述设定值为15°。
3.根据权利要求1所述一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,其特征在于:步骤1中,对于初始翼型以及目标翼型各自的上翼型和下翼型,分别进行以下操作来计算每一段翼型弯曲角度:
首先拟合翼型曲线,之后通过曲线平移和旋转,使得翼型的后缘处于坐标系原点并在该点和坐标系x轴相切;
其次将翼型沿弦长等分,每一段翼型对应弦长相等,根据拟合的翼型曲线计算每一段翼型两端点的斜率T1和T2,T2>T1,根据公式
得到该分段的弯曲角度β。
4.根据权利要求1所述一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,其特征在于:负泊松比蜂窝单胞以及正泊松比蜂窝单胞的胞壁厚度和高度不大于4mm,且胞壁厚度小于高度。
5.根据权利要求1所述一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,其特征在于:正泊松比蜂窝单胞胞壁厚度t2等于负泊松比蜂窝单胞胞壁厚度t1;正泊松比蜂窝单胞高度b2等于负泊松比蜂窝单胞高度b1。
6.根据权利要求1所述一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮设计方法,其特征在于:正泊松比蜂窝多胞以及负泊松比蜂窝多胞的y方向顶部和底部的水平胞壁的厚度为对应蜂窝单胞胞壁厚度的两倍。
7.一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮,其特征在于:所述柔性蒙皮具有蜂窝夹芯,所述蜂窝夹芯由权利要求1~6任一所述方法设计得到。
8.一种权利要求7所述柔性蒙皮中蜂窝夹芯主动变形的控制方法,其特征在于:根据公式
计算蜂窝夹芯上各个变形单元的驱动位移,并通过形状记忆合金驱动器产生驱动位移,将蜂窝夹芯从初始翼型形状变形为目标翼型形状;其中δyj为第j个翼型分段对应的变形单元的驱动位移,T是正、负泊松比蜂窝多胞的高度和,μ是修正系数,αj为机翼初始翼型第j个分段和目标翼型第j个分段的弯曲角度之差,αj=αj目标-αj初始。
9.根据权利要求8所述控制方法,其特征在于:修正系数μ有限元数值分析进行数值模拟实验得到:
首先建立变形单元模型,在驱动位移小于或等于1.5mm的前提下至少选取5个不同的驱动位移δy施加给变形单元,每个驱动位移构成一组数值模拟实验,共进行s组数值模拟实验,每组实验的驱动位移δyi使得变形单元产生绕y轴的弯曲角度αi;
然后根据公式
计算得到修正系数μ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210876282.6A CN115196036B (zh) | 2022-07-25 | 2022-07-25 | 一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210876282.6A CN115196036B (zh) | 2022-07-25 | 2022-07-25 | 一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115196036A CN115196036A (zh) | 2022-10-18 |
CN115196036B true CN115196036B (zh) | 2024-04-26 |
Family
ID=83584085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210876282.6A Active CN115196036B (zh) | 2022-07-25 | 2022-07-25 | 一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115196036B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115985419B (zh) * | 2022-12-11 | 2024-06-11 | 西南交通大学 | 具有梯度泊松比分布特性的夹心梁蜂窝芯层结构的设计方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106800084A (zh) * | 2017-02-10 | 2017-06-06 | 哈尔滨工业大学 | 一种具有可设计泊松比的变形蒙皮结构 |
WO2020079424A1 (en) * | 2018-10-19 | 2020-04-23 | Cranfield University | Materials with structures exhibiting zero poisson's ratio |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2463320A (en) * | 2008-09-15 | 2010-03-17 | Auxetic Technologies Ltd | Multi-layer adhesive interface containing auxetic material |
EP3135949B1 (en) * | 2015-08-27 | 2019-06-19 | Airbus Operations S.L. | Deformable structure for absorption of energy from mechanical and/or acoustic impacts |
-
2022
- 2022-07-25 CN CN202210876282.6A patent/CN115196036B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106800084A (zh) * | 2017-02-10 | 2017-06-06 | 哈尔滨工业大学 | 一种具有可设计泊松比的变形蒙皮结构 |
WO2020079424A1 (en) * | 2018-10-19 | 2020-04-23 | Cranfield University | Materials with structures exhibiting zero poisson's ratio |
Non-Patent Citations (1)
Title |
---|
主动变形负泊松比蜂窝结构用于变体机翼;陶伟灏;;设备管理与维修;20171125(第17期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN115196036A (zh) | 2022-10-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Demir et al. | Design of variable stiffness composite structures using lamination parameters with fiber steering constraint | |
CN115196036B (zh) | 一种基于正、负泊松比蜂窝芯主动变形的柔性蒙皮及设计和驱动控制方法 | |
CN104866673A (zh) | 一种轴压加筋柱壳的开口补强方法 | |
CN101650756A (zh) | 一种多点成形飞机蒙皮的方法 | |
CN111709097A (zh) | 一种基于零亏格网格曲面连续变形的柔顺机构生成方法 | |
CN106885676A (zh) | 气动载荷产生的非解耦六自由度机构末端位姿误差补偿法 | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
CN110955941B (zh) | 基于向量场的复合材料结构优化设计方法及设备 | |
CN108108582A (zh) | 一种曲面件柔性轧制成形过程的数值模拟方法 | |
CN106250605B (zh) | 基于精确变形控制的薄板压电智能结构协同拓扑优化方法 | |
CN111177906A (zh) | 一种离散化模具型面精确补偿方法 | |
CN113886967A (zh) | 多巡航工况的大型飞机机翼气动弹性优化方法 | |
CN109670271A (zh) | 基于敏感点多点增量成形的大型薄壁加筋板的校形方法 | |
CN115635004A (zh) | 一种马鞍外形带筋壁板喷丸成形方法 | |
CN113779831B (zh) | 一种基于区域分解的缩聚feti工程数值方法 | |
CN113844636B (zh) | 一种ω形柔性蒙皮蜂窝结构 | |
CN113232833B (zh) | 一种形状记忆合金拉线驱动的变弯度机翼及其设计方法 | |
Zhang et al. | A morphing wing with cellular structure of non-uniform density | |
CN110543655A (zh) | 用于纳米压印设备的精密对准平台的优化设计方法及装置 | |
CN112733291B (zh) | 一种复杂曲率曲面的模具型面精准补偿方法 | |
CN102354324B (zh) | 跨音速颤振模型复合材料单梁结构设计及其刚度计算方法 | |
CN112733253B (zh) | 一种波纹板形式柔性后缘机翼结构的设计方法 | |
Zhang et al. | An intelligent method to design die profile for rubber forming of complex curved flange part | |
CN115221622A (zh) | 一种大尺寸复合材料机身壁板装配定位夹持布局优化方法 | |
CN114282410A (zh) | 一种考虑气动结构协调更新的变形机翼有限元建模方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Liang Ke Inventor after: Cheng Qian Inventor before: Cheng Qian Inventor before: Liang Ke |
|
GR01 | Patent grant | ||
GR01 | Patent grant |