CN107506529A - 一种复合材料加筋壁板轴压稳定性计算方法 - Google Patents

一种复合材料加筋壁板轴压稳定性计算方法 Download PDF

Info

Publication number
CN107506529A
CN107506529A CN201710645288.1A CN201710645288A CN107506529A CN 107506529 A CN107506529 A CN 107506529A CN 201710645288 A CN201710645288 A CN 201710645288A CN 107506529 A CN107506529 A CN 107506529A
Authority
CN
China
Prior art keywords
mrow
msub
mfrac
msup
eta
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.)
Granted
Application number
CN201710645288.1A
Other languages
English (en)
Other versions
CN107506529B (zh
Inventor
苏雁飞
赵占军
惠红军
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201710645288.1A priority Critical patent/CN107506529B/zh
Publication of CN107506529A publication Critical patent/CN107506529A/zh
Application granted granted Critical
Publication of CN107506529B publication Critical patent/CN107506529B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及一种复合材料加筋壁板轴压稳定性计算方法,首先确定复合材料加筋壁板的计算模型,由于所述复合材料加筋壁板由蒙皮和多个均匀排布的筋条组成,故以相邻两个筋条之间的典型单元进行分析,并建立坐标系;其次确定筋条对蒙皮的支持系数,设一筋条对蒙皮的支持系数为kL,另一筋条对蒙皮的支持系数为kR,则当kL=kR=0时,筋条对蒙皮的支持状态为简支,当kL=kR=∞时,筋条对蒙皮的支持状态为固支,但筋条对蒙皮的支持状态介于简支与固支之间,故建立筋条对蒙皮支持系数k的表达式;最后确定复合材料加筋壁板轴压稳定性临界屈曲载荷。本发明的复合材料加筋壁板轴压载荷计算方法时,计算精度与试验结果非常接近,计算精度高、计算简便。

Description

一种复合材料加筋壁板轴压稳定性计算方法
技术领域
本发明属于航空结构强度分析领域,尤其涉及一种复合材料加筋壁板轴压稳定性计算方法。
背景技术
现有的复合材料设计手册及飞机设计参考资料均没有考虑筋条对复合材料加筋板轴压稳定性的影响,也没有给出筋条对蒙皮弹性支持(支持状态介于简支和固支之间)时,加筋壁板轴压载荷下的稳定性计算公式。通过对大量的试验数据分析整理,应用目前现有的工程算法得到的初始屈曲载荷与试验值的偏差比较大,有的计算误差甚至高达30%~40%。究其原因,应用工程计算公式时,轴压稳定性的边界条件按照四边简支考虑,导致计算结果太过保守;然而如果按固支考虑,又偏于危险,不利于飞机安全。事实上,加筋条对蒙皮的支持介于简支和固支之间,属于弹性支持的范畴。为了提高分析精度,提高复合材料结构使用效率,迫切需要对工程分析方法进行修正,得到精确的计算复合材料加筋板轴压初始屈曲载荷的公式。
发明内容
本发明的目的是提供一种复合材料加筋壁板轴压稳定性计算方法,用于解决目前在计算复合材料加筋壁板的轴压载荷方法中,筋条对蒙皮的影响仅有简支和固支两种极端状态的弊端,使得对其的计算结果偏差较大的问题。
为达到上述目的,本发明采用的技术方案是:一种复合材料加筋壁板轴压稳定性计算方法,首先确定复合材料加筋壁板的计算模型,由于所述复合材料加筋壁板由蒙皮和多个均匀排布的筋条组成,故以相邻两个筋条之间的典型单元进行分析,并建立坐标系;其次确定筋条对蒙皮的支持系数,设一筋条对蒙皮的支持系数为kL,另一筋条对蒙皮的支持系数为kR,则当kL=kR=0时,筋条对蒙皮的支持状态为简支,当kL=kR=∞时,筋条对蒙皮的支持状态为固支,但筋条对蒙皮的支持状态介于简支与固支之间,故建立筋条对蒙皮支持系数k的表达式;最后确定复合材料加筋壁板轴压稳定性临界屈曲载荷。
进一步地,本发明中的复合材料加筋壁板设其长度为a,宽度为B,筋条间距为b,受到单位长度的均布轴压载荷Nx作用,因此所述坐标系建立准则为:以复合材料加筋壁板的长度方向为x方向,宽度方向为y方向,z向符合右手坐标系。
进一步地,在本发明中,建立的筋条对蒙皮支持系数k的表达式具体为:
式中:Dij(i,j=1,2,6)为复合材料层压板的弯曲刚度矩阵,角标skin指代蒙皮,stringer指代筋条。
进一步地,在本发明中,确定复合材料加筋壁板轴压稳定性临界屈曲载荷的过程具体为:
对于建立的典型单元,由于筋条对蒙皮的支持状态介于简支与固支之间,即属于弹性支持,因此设定z方向的位移函数w(x,y)表达式为:
式中:β为常数;m为复合材料加筋壁板沿着x方向的屈曲半波数; 为常数项,其由位移边界条件确定;
其中位移边界条件为:
根据位移边界条件最终确定位移函数w(x,y)表达式为:
之后,复合材料加筋壁板的弹性应变能记为Ue,沿着板弹性支持边的弹性应变能记为UΓ,沿着加载方向的载荷Nx所做的外力功记为V,则复合材料加筋壁板的总势能∏为:Π=Ue+UΓ-V;
根据最小势能原理,则有:δΠ=δUe+δUΓ-δV=0
其中:
将位移函数w(x,y)表达式带入最小势能原理,可得屈曲载荷Nx为:
可得临界屈曲载荷表达式为:
式中,各系数表达式如下:
本发明的复合材料加筋壁板轴压稳定性计算方法首次提出了一种涉及筋条支持效应的复合材料加筋壁板轴压稳定性计算方法,解决了现有参考资料在计算复合材料加筋板轴压屈曲载荷时,认为筋条的支持状态只有简支和固支的两种极端状态,因而并没有给出当筋条对蒙皮的支持介于简支和固支之间(即筋条对蒙皮弹性支持)时的屈曲载荷计算公式的问题,而本发明则通过公式推导了筋条对蒙皮弹性支持时,复合材料加筋板轴压稳定性计算公式,进而提出的弹性支持方式与实际结构支持状态更接近。本发明的复合材料加筋壁板轴压稳定性计算方法解决了目前工程算法在计算复合材料的加筋板稳定性时,所得计算结果精度更高、计算过程简便。
附图说明
此处的附图被并入说明书中并构成本说明书的一部分,示出了符合本发明的实施例,并与说明书一起用于解释本发明的原理。
图1为本发明中的复合材料加筋壁板结构示意图。
图2为本发明中的典型单元结构图。
图3为本发明一实施例的复合材料筋条截面图。
图4为采用本发明的计算方法与试验结果对比图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。
本发明的复合材料加筋壁板轴压稳定性计算方法,包括以下步骤:
步骤一、确定复合材料加筋壁板计算模型
如图1所示的复合材料加筋壁板,复合材料加筋壁板由蒙皮和筋条构成,其中,蒙皮及筋条均为复合材料制成。复合材料加筋壁板的整体长度(也就是筋条的长度)为a,宽度为B,筋条间距为b,受到单位长度的均布轴压载荷Nx作用,轴压载荷Nx的作用方向与筋条的轴线平行。取两个相邻筋条之间的典型单元进行分析,如图2所示。定义复合材料加壁板的长度方向(也就是筋条轴线方向)为x方向,宽度方向为y方向,z向(垂直于xy平面)符合右手坐标系。左侧的筋条对蒙皮的支持系数记为kL,右侧筋条对蒙皮的支持系数记为kR
步骤二、确定筋条对蒙皮的支持系数k
当kL=kR=0时,为筋条对蒙皮的支持状态为简支,kL=kR=∞时,为筋条对蒙皮的支持为固支。但实际上,筋条对蒙皮的支持状态介于简支与固支之间,因此建立筋条对蒙皮的支持系数k,其支持系数k的表达式为:
式中:Dij(i,j=1,2,6)为复合材料层压板的弯曲刚度矩阵,角标skin指代蒙皮,stringer指代筋条。
步骤三、确定复合材料加筋壁板轴压稳定性临界屈曲载荷
对于图2中的典型单元,由于筋条对蒙皮的支持属于弹性支持(筋条对蒙皮的支持状态介于简支与固支之间的状态称为弹性支持),设定z方向的位移函数w(x,y)表达式为:
上式中,β为常数;m为板沿着x方向的屈曲半波数;为常数项,可由位移边界条件确定。
位移边界条件如下:
最后根据边界条件可以最终确定位移函数w(x,y)表达式为:
其中,Dij(i,j=1,2,6)为复合材料层压板的弯曲刚度矩阵。
复合材料加筋壁板的弹性应变能记为Ue,沿着板弹性支持边的弹性应变能记为UΓ,沿着加载方向外力Nx所做的外力功记为V,板的总势能∏为:
Π=Ue+UΓ-V
根据最小势能原理,则有:δΠ=δUe+δUΓ-δV=0
其中:
之后,将位移函数w(x,y)表达式带入最小势能原理,可得屈曲载荷Nx为:
可得临界屈曲载荷表达式为:
上式中,各系数表达式如下:
通过上述计算及公式,即可求得临界屈曲载荷
为了使本发明的复合材料加筋壁板轴压稳定性计算方法更加明了,下面以某一具体数据对本发明的复合材料加筋壁板轴压稳定性计算方法做更进一步地详细说明。
已知:某飞机复合材料加筋壁板结构,其中蒙皮所用材料为CF3031/3238A织物,长桁(与本发明中的筋条为同一作用)所用材料为CF3031/3238A织物和CCF300/323A单向带,蒙皮的铺层方式为[(±45)/(0/90)/(±45)/(0/90)/(±45)],长桁铺层方式为:[(±45)/02/(±45)/02/(±45)]s,受到轴压载荷Nx作用。复合材料加筋壁板的长度为a=500mm,长桁间距为b=200mm,长桁剖面(截面)为T型材,长桁的高度为20mm,型材顶面宽度为25mm,长桁的剖面形状及尺寸见图3所示,在本实施例中所用复材的材料性能见表1。
表1材料特性表
材料牌号 纵向拉伸模量MPa 横向拉伸模量MPa 纵横剪切模量MPa 泊松比
3238A/CCF300 125000 8280 3700 0.3
3238A/CF3031 59200 58000 3770 0.054
首先需要确定计算的模型,设筋条对蒙皮的支持系数为k,由于两侧的筋条方式一致,故kL=KR,且a=500mm,b=200mm;
然后确定筋条对蒙皮的支持系数k
根据蒙皮以及长桁的铺层方式和材料性能,得到长桁、蒙皮的弯曲刚度系数为见表2:
表2弯曲刚度系数表
结合图1和表1中数据,得到T型长桁对蒙皮的支持系数k为:
最后确定复合材料加筋壁板轴压稳定性临界载荷不过首先需求得各中间参数,结果为:
进一步求得临界屈曲载荷为
现有的复合材料设计手册及飞机设计参考资料在计算复合材料加筋板轴压屈曲载荷时,仅有简支边界、固支边界下的工程计算公式,并没有给出当筋条对蒙皮的支持介于简支和固支之间(即筋条对蒙皮的弹性支持)时的屈曲载荷计算公式。本发明的复合材料加筋壁板计算方法填补了这一空白,在计算复合材料加筋板轴压屈曲载荷时加入筋条对蒙皮弹性支持,求得复合材料加筋板轴压稳定性计算公式,所得结果与试验值吻合很好,实施效果可参考下图4。在图4中可见,当筋条对蒙皮的支持仅为简支时,计算结果偏小(较为保守,余量较大),当当筋条对蒙皮的支持仅为固支时,计算结果偏大(较为危险,余量较小),若取在简支与固支之间的平均值,没有相关依据,且计算偏差也较大,但采用本发明的复合材料加筋壁板轴压载荷计算方法时,计算精度与试验结果非常接近,计算精度高、计算简便。
本发明首次提出了一种涉及筋条支持效应的复合材料加筋壁板轴压稳定性计算方法,提出解决了目前工程算法在计算加筋板的稳定性时,认为筋条的支持状态只有简支和固支的两种极端状态的弊端,提出的弹性支持方式与实际结构支持状态更接近,所得计算结果精度更高,该复合材料加筋壁板轴压载荷计算方法在目前的技术领域是一重大突破。
以上所述,仅为本发明的最优具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。

Claims (4)

1.一种复合材料加筋壁板轴压稳定性计算方法,其特征在于:
首先确定复合材料加筋壁板的计算模型,由于所述复合材料加筋壁板由蒙皮和多个均匀排布的筋条组成,故以相邻两个筋条之间的典型单元进行分析,并建立坐标系;
其次确定筋条对蒙皮的支持系数,设一筋条对蒙皮的支持系数为kL,另一筋条对蒙皮的支持系数为kR,则当kL=kR=0时,筋条对蒙皮的支持状态为简支,当kL=kR=∞时,筋条对蒙皮的支持状态为固支,但筋条对蒙皮的支持状态介于简支与固支之间,故建立筋条对蒙皮支持系数k的表达式;
最后确定复合材料加筋壁板轴压稳定性临界屈曲载荷。
2.根据权利要求1的复合材料加筋壁板轴压稳定性计算方法,其特征在于,复合材料加筋壁板其长度为a,宽度为B,筋条间距为b,受到单位长度的均布轴压载荷Nx作用,因此所述坐标系建立准则为:以复合材料加筋壁板的长度方向为x方向,宽度方向为y方向,z向符合右手坐标系。
3.根据权利要求2的复合材料加筋壁板轴压稳定性计算方法,其特征在于,建立的筋条对蒙皮支持系数k的表达式具体为:
<mrow> <mi>k</mi> <mo>=</mo> <mfrac> <mrow> <munder> <mo>&amp;Sigma;</mo> <mi>i</mi> </munder> <msub> <mrow> <mo>&amp;lsqb;</mo> <mfrac> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mn>11</mn> </msub> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>D</mi> <mn>66</mn> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>D</mi> <mn>12</mn> </msub> <msub> <mi>D</mi> <mn>26</mn> </msub> <msub> <mi>D</mi> <mn>16</mn> </msub> <mo>-</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <msubsup> <mi>D</mi> <mn>16</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>D</mi> <mn>66</mn> </msub> <msubsup> <mi>D</mi> <mn>12</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>D</mi> <mn>11</mn> </msub> <msubsup> <mi>D</mi> <mn>26</mn> <mn>2</mn> </msubsup> <mo>)</mo> <mo>&amp;CenterDot;</mo> <mn>4</mn> <mi>b</mi> </mrow> <mrow> <msub> <mi>D</mi> <mn>11</mn> </msub> <msub> <mi>D</mi> <mn>22</mn> </msub> <mo>-</mo> <msubsup> <mi>D</mi> <mn>12</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mi>i</mi> <mo>-</mo> <mi>s</mi> <mi>t</mi> <mi>r</mi> <mi>i</mi> <mi>n</mi> <mi>g</mi> <mi>e</mi> <mi>r</mi> </mrow> </msub> </mrow> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <msub> <mrow> <mo>(</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>s</mi> <mi>k</mi> <mi>i</mi> <mi>n</mi> </mrow> </msub> </mrow> </mfrac> </mrow>
式中:Dij(i,j=1,2,6)为复合材料层压板的弯曲刚度矩阵,角标skin指代蒙皮,stringer指代筋条。
4.根据权利要求3的复合材料加筋壁板轴压稳定性计算方法,其特征在于,确定复合材料加筋壁板轴压稳定性临界屈曲载荷的过程具体为:
对于建立的典型单元,由于筋条对蒙皮的支持状态介于简支与固支之间,即属于弹性支持,因此设定z方向的位移函数w(x,y)表达式为:
式中:β为常数;m为复合材料加筋壁板沿着x方向的屈曲半波数; 为常数项,其由位移边界条件确定;
其中位移边界条件为:
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>w</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>w</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>b</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>M</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mn>0</mn> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>y</mi> <mo>=</mo> <mn>0</mn> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mi>L</mi> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>y</mi> <mo>=</mo> <mn>0</mn> </mrow> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>M</mi> <mi>y</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>b</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>y</mi> <mo>=</mo> <mi>b</mi> </mrow> </msub> <mo>=</mo> <mo>-</mo> <msub> <mi>k</mi> <mi>R</mi> </msub> <msub> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mi>y</mi> <mo>=</mo> <mi>b</mi> </mrow> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced>
根据位移边界条件最终确定位移函数w(x,y)表达式为:
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>w</mi> <mrow> <mo>(</mo> <mrow> <mi>x</mi> <mo>,</mo> <mi>y</mi> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <mo>{</mo> <mfrac> <mi>y</mi> <mi>b</mi> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>k</mi> <mi>L</mi> </msub> <mi>b</mi> </mrow> <mrow> <mn>2</mn> <msub> <mi>D</mi> <mn>22</mn> </msub> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>y</mi> <mi>b</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mfrac> <mrow> <mn>12</mn> <msubsup> <mi>D</mi> <mn>22</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <mrow> <mo>(</mo> <mrow> <mn>5</mn> <msub> <mi>k</mi> <mi>L</mi> </msub> <mo>-</mo> <mn>3</mn> <msub> <mi>k</mi> <mi>R</mi> </msub> </mrow> <mo>)</mo> </mrow> <mi>b</mi> <mo>-</mo> <msup> <mi>b</mi> <mn>2</mn> </msup> <mo>&amp;CenterDot;</mo> <msub> <mi>k</mi> <mi>L</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>k</mi> <mi>R</mi> </msub> </mrow> <mrow> <mn>6</mn> <msubsup> <mi>D</mi> <mn>22</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>bk</mi> <mi>R</mi> </msub> <msub> <mi>D</mi> <mn>22</mn> </msub> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mi>y</mi> <mi>b</mi> </mfrac> <mo>)</mo> </mrow> <mn>3</mn> </msup> <mo>+</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mn>12</mn> <msubsup> <mi>D</mi> <mn>22</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <mrow> <mo>(</mo> <mrow> <mn>4</mn> <msub> <mi>k</mi> <mi>L</mi> </msub> <mo>-</mo> <mn>4</mn> <msub> <mi>k</mi> <mi>R</mi> </msub> </mrow> <mo>)</mo> </mrow> <mi>b</mi> <mo>-</mo> <msup> <mi>b</mi> <mn>2</mn> </msup> <mo>&amp;CenterDot;</mo> <msub> <mi>k</mi> <mi>L</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>k</mi> <mi>R</mi> </msub> </mrow> <mrow> <mn>12</mn> <msubsup> <mi>D</mi> <mn>22</mn> <mn>2</mn> </msubsup> <mo>-</mo> <mn>2</mn> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>bk</mi> <mi>R</mi> </msub> </mrow> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mi>y</mi> <mi>b</mi> </mfrac> <mo>)</mo> </mrow> <mn>4</mn> </msup> <mo>}</mo> <mo>&amp;times;</mo> <mi>&amp;beta;</mi> <mo>&amp;CenterDot;</mo> <mi>sin</mi> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> <mi>x</mi> </mrow> <mi>a</mi> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
之后,复合材料加筋壁板的弹性应变能记为Ue,沿着板弹性支持边的弹性应变能记为UΓ,沿着加载方向的载荷Nx所做的外力功记为V,则复合材料加筋壁板的总势能∏为:Π=Ue+UΓ-V;
根据最小势能原理,则有:δΠ=δUe+δUΓ-δV=0
其中:
<mrow> <msub> <mi>U</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <mi>&amp;Omega;</mi> </munder> <mo>{</mo> <msub> <mi>D</mi> <mn>11</mn> </msub> <mo>&amp;CenterDot;</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <mo>&amp;CenterDot;</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <mn>2</mn> <msub> <mi>D</mi> <mn>12</mn> </msub> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>x</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </mfrac> <mo>+</mo> <mn>4</mn> <msub> <mi>D</mi> <mn>66</mn> </msub> <mo>&amp;CenterDot;</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>}</mo> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow>
<mrow> <msub> <mi>U</mi> <mi>&amp;Gamma;</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munder> <mo>&amp;Integral;</mo> <mi>&amp;Gamma;</mi> </munder> <msub> <mi>k</mi> <mi>L</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>y</mi> <mo>=</mo> <mn>0</mn> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>d</mi> <mi>&amp;Gamma;</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <munder> <mo>&amp;Integral;</mo> <mi>&amp;Gamma;</mi> </munder> <msub> <mi>k</mi> <mi>R</mi> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>y</mi> </mrow> </mfrac> <msub> <mo>|</mo> <mrow> <mi>y</mi> <mo>=</mo> <mi>b</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>d</mi> <mi>&amp;Gamma;</mi> </mrow>
<mrow> <mi>V</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>N</mi> <mi>x</mi> </msub> <munder> <mrow> <mo>&amp;Integral;</mo> <mo>&amp;Integral;</mo> </mrow> <mi>&amp;Omega;</mi> </munder> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mi>w</mi> </mrow> <mrow> <mo>&amp;part;</mo> <mi>x</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mi>d</mi> <mi>x</mi> <mi>d</mi> <mi>y</mi> </mrow>
将位移函数w(x,y)表达式带入最小势能原理,可得屈曲载荷Nx为:
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>N</mi> <mi>x</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mfrac> <mn>1</mn> <mn>4</mn> </mfrac> <msub> <mi>D</mi> <mn>11</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mn>4</mn> </msup> <msub> <mi>b&amp;eta;</mi> <mn>1</mn> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mn>4</mn> </mfrac> <msub> <mi>D</mi> <mn>22</mn> </msub> <mo>&amp;CenterDot;</mo> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>3</mn> </msup> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>D</mi> <mn>12</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <mo>&amp;CenterDot;</mo> <msub> <mi>&amp;eta;</mi> <mn>3</mn> </msub> <mo>+</mo> <msub> <mi>D</mi> <mn>66</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mfrac> <mn>1</mn> <mi>b</mi> </mfrac> <msub> <mi>&amp;eta;</mi> <mn>4</mn> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mrow> <mn>4</mn> <msup> <mi>b</mi> <mn>2</mn> </msup> </mrow> </mfrac> <msub> <mi>&amp;eta;</mi> <mn>5</mn> </msub> </mrow> <mrow> <mfrac> <mn>1</mn> <mn>4</mn> </mfrac> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>b&amp;eta;</mi> <mn>6</mn> </msub> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msub> <mi>D</mi> <mn>11</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>+</mo> <msub> <mi>D</mi> <mn>22</mn> </msub> <msup> <mrow> <mo>(</mo> <mfrac> <mi>a</mi> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>4</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>-</mo> <mn>2</mn> <msub> <mi>D</mi> <mn>12</mn> </msub> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>3</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>+</mo> <mn>4</mn> <msub> <mi>D</mi> <mn>66</mn> </msub> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>4</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>+</mo> <msup> <mrow> <mo>(</mo> <mfrac> <mi>a</mi> <mrow> <mi>m</mi> <mi>&amp;pi;</mi> </mrow> </mfrac> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>3</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>5</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
可得临界屈曲载荷表达式为:
<mrow> <msubsup> <mi>N</mi> <mi>x</mi> <mrow> <mi>c</mi> <mi>r</mi> </mrow> </msubsup> <mo>=</mo> <mfrac> <msqrt> <mrow> <msub> <mi>D</mi> <mn>11</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mrow> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>b&amp;eta;</mi> <mn>5</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </msqrt> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mrow> </mfrac> <mo>+</mo> <mfrac> <mrow> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <msqrt> <mrow> <msub> <mi>D</mi> <mn>11</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> </mrow> </msqrt> </mrow> <mrow> <msup> <mi>b</mi> <mn>2</mn> </msup> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> <msqrt> <mrow> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>b&amp;eta;</mi> <mn>5</mn> </msub> </mrow> </msqrt> </mrow> </mfrac> <mo>-</mo> <mn>2</mn> <msub> <mi>D</mi> <mn>12</mn> </msub> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>3</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>+</mo> <mn>4</mn> <msub> <mi>D</mi> <mn>66</mn> </msub> <mfrac> <mn>1</mn> <msup> <mi>b</mi> <mn>2</mn> </msup> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>4</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>6</mn> </msub> </mfrac> <mo>+</mo> <mfrac> <mrow> <msqrt> <mrow> <msub> <mi>D</mi> <mn>11</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> </mrow> </msqrt> <mo>&amp;CenterDot;</mo> </mrow> <msqrt> <mrow> <msub> <mi>D</mi> <mn>22</mn> </msub> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>b&amp;eta;</mi> <mn>5</mn> </msub> </mrow> </msqrt> </mfrac> <mfrac> <msub> <mi>&amp;eta;</mi> <mn>5</mn> </msub> <mrow> <msub> <mi>b&amp;eta;</mi> <mn>6</mn> </msub> </mrow> </mfrac> </mrow>
式中,各系数表达式如下:
CN201710645288.1A 2017-08-01 2017-08-01 一种复合材料加筋壁板轴压稳定性计算方法 Active CN107506529B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710645288.1A CN107506529B (zh) 2017-08-01 2017-08-01 一种复合材料加筋壁板轴压稳定性计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710645288.1A CN107506529B (zh) 2017-08-01 2017-08-01 一种复合材料加筋壁板轴压稳定性计算方法

Publications (2)

Publication Number Publication Date
CN107506529A true CN107506529A (zh) 2017-12-22
CN107506529B CN107506529B (zh) 2020-09-18

Family

ID=60688993

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710645288.1A Active CN107506529B (zh) 2017-08-01 2017-08-01 一种复合材料加筋壁板轴压稳定性计算方法

Country Status (1)

Country Link
CN (1) CN107506529B (zh)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108256214A (zh) * 2018-01-16 2018-07-06 滨州学院 一种航空正交加筋板刚度计算方法及装置
CN108710736A (zh) * 2018-05-10 2018-10-26 北京航空航天大学 一种考虑应力约束的宏微一体结构拓扑优化方法
CN109635448A (zh) * 2018-12-14 2019-04-16 西安石油大学 一种回转装置临界涡动速度的计算方法
CN111143941A (zh) * 2019-12-24 2020-05-12 中国航空工业集团公司西安飞机设计研究所 一种复合材料加筋壁板轴压承载能力计算方法
CN111274670A (zh) * 2019-12-31 2020-06-12 中国航空工业集团公司沈阳飞机设计研究所 一种钛合金型材组合壁板轴压承载能力试验设计方法
CN111553031A (zh) * 2020-04-24 2020-08-18 中国飞机强度研究所 一种整体加筋板极限载荷计算方法
CN111948043A (zh) * 2020-07-24 2020-11-17 中国飞机强度研究所 一种拉压剪复合载荷作用的加筋壁板屈曲预测方法
CN112560182A (zh) * 2020-12-17 2021-03-26 中航沈飞民用飞机有限责任公司 一种基于稳定性特征的民机复合材料加筋壁板长桁刚度确定方法
CN112711798A (zh) * 2020-12-29 2021-04-27 中国航空工业集团公司西安飞机设计研究所 一种四边固支复合材料层压板轴压稳定性计算方法
CN115031590A (zh) * 2022-05-19 2022-09-09 上海交通大学 一种承担接头载荷的加筋壁板结构稳定性优化方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103293056A (zh) * 2013-04-19 2013-09-11 空军工程大学 加筋板结构轴向压缩稳定性试验夹具及方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103293056A (zh) * 2013-04-19 2013-09-11 空军工程大学 加筋板结构轴向压缩稳定性试验夹具及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PIZHONG QIAO 等: "Explicit local buckling analysis and design of fiber-reinforced plastic composite structural shapes", 《COMPOSITE STRUCTURES》 *
苏雁飞 等: "轴压载荷下复合材料薄壁加筋板后屈曲承载能力分析与试验验证", 《航空科学技术》 *
陈金睿 等: "考虑筋条扭转弹性支持的轴压复合材料加筋板局部屈曲分析方法", 《南京航空航天大学学报》 *

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108256214A (zh) * 2018-01-16 2018-07-06 滨州学院 一种航空正交加筋板刚度计算方法及装置
CN108710736A (zh) * 2018-05-10 2018-10-26 北京航空航天大学 一种考虑应力约束的宏微一体结构拓扑优化方法
CN109635448A (zh) * 2018-12-14 2019-04-16 西安石油大学 一种回转装置临界涡动速度的计算方法
CN111143941B (zh) * 2019-12-24 2022-10-11 中国航空工业集团公司西安飞机设计研究所 一种复合材料加筋壁板轴压承载能力计算方法
CN111143941A (zh) * 2019-12-24 2020-05-12 中国航空工业集团公司西安飞机设计研究所 一种复合材料加筋壁板轴压承载能力计算方法
CN111274670A (zh) * 2019-12-31 2020-06-12 中国航空工业集团公司沈阳飞机设计研究所 一种钛合金型材组合壁板轴压承载能力试验设计方法
CN111274670B (zh) * 2019-12-31 2023-08-08 中国航空工业集团公司沈阳飞机设计研究所 一种钛合金型材组合壁板轴压承载能力试验设计方法
CN111553031A (zh) * 2020-04-24 2020-08-18 中国飞机强度研究所 一种整体加筋板极限载荷计算方法
CN111948043A (zh) * 2020-07-24 2020-11-17 中国飞机强度研究所 一种拉压剪复合载荷作用的加筋壁板屈曲预测方法
CN112560182A (zh) * 2020-12-17 2021-03-26 中航沈飞民用飞机有限责任公司 一种基于稳定性特征的民机复合材料加筋壁板长桁刚度确定方法
CN112711798B (zh) * 2020-12-29 2022-10-11 中国航空工业集团公司西安飞机设计研究所 一种四边固支复合材料层压板轴压稳定性计算方法
CN112711798A (zh) * 2020-12-29 2021-04-27 中国航空工业集团公司西安飞机设计研究所 一种四边固支复合材料层压板轴压稳定性计算方法
CN115031590A (zh) * 2022-05-19 2022-09-09 上海交通大学 一种承担接头载荷的加筋壁板结构稳定性优化方法

Also Published As

Publication number Publication date
CN107506529B (zh) 2020-09-18

Similar Documents

Publication Publication Date Title
CN107506529A (zh) 一种复合材料加筋壁板轴压稳定性计算方法
CN106156449A (zh) 一种复合材料机翼壁板优化设计方法
CN106055764B (zh) 基于三维壳有限元-梁模型的风力机叶片位移计算方法
CN106758738B (zh) 弹性转动约束边界钢-混组合梁高腹板纵向加劲肋设置的方法
Patton et al. Numerical modeling of lean duplex stainless steel hollow columns of square, L-, T-, and+-shaped cross sections under pure axial compression
CN107423468A (zh) 一种基于等效方法的prseus结构分析方法
Kwan Shear lag in shear/core walls
CN107085631A (zh) 一种基于舱门细节模型的强度校核方法
CN105844033A (zh) 船舶总纵极限强度简化逐步破坏分析方法
CN109101692B (zh) 基于最大应力准则的复合材料层合板极限载荷计算方法
CN110852012A (zh) 轴压加筋板极限强度预测方法
Ali et al. Finite element analysis of cold-formed steel connections
Polilov et al. Estimating the effect of misorientation of fibers on stiffness and strength of profiled composite elements
CN111914351B (zh) 一种机身结构加筋壁板总体稳定性计算方法
CN103397711A (zh) 多管格构式屈曲约束支撑
CN106940736A (zh) 一种考虑侧向载荷作用的纵骨多跨失稳的极限载荷确定方法
CN106055731A (zh) 一种复合材料层合板非概率可靠性优化方法
Tosh et al. Fibre steering for a composite C-beam
CN102519715B (zh) 一种核电站支架强度简化计算方法
CN107908870A (zh) 一种体内预应力胶合木梁抗弯承载能力的分析方法
CN111539107A (zh) 一种格栅增强复合材料圆柱壳体临界轴压载荷快速预测方法
CN115630458A (zh) 一种基于弹性地基梁理论的现浇连续梁的方法及其应用
CN105501462B (zh) 一种机翼结构设计方法
Jáger et al. Experimental based numerical modelling of girders with trapezoidally corrugated web subjected to combined loading
CN105631162B (zh) 一种开口加筋矩形层合薄板的基频预报方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
CB03 Change of inventor or designer information

Inventor after: Su Yanfei

Inventor after: Zhao Zhanwen

Inventor after: Hui Hongjun

Inventor before: Su Yanfei

Inventor before: Zhao Zhanjun

Inventor before: Hui Hongjun

CB03 Change of inventor or designer information
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant