CN103324782B - 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 - Google Patents
一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 Download PDFInfo
- Publication number
- CN103324782B CN103324782B CN201310207200.XA CN201310207200A CN103324782B CN 103324782 B CN103324782 B CN 103324782B CN 201310207200 A CN201310207200 A CN 201310207200A CN 103324782 B CN103324782 B CN 103324782B
- Authority
- CN
- China
- Prior art keywords
- compound substance
- load
- substance pressurized
- bounded
- pressurized
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,该方法首先将复合材料分散性参数区间量化,并在不确定性参数的区间内配点,基于复合材料受压蒙皮结构特点,建立复合材料受压蒙皮配点型区间有限元分析模型,并建立复合材料受压蒙皮稳定性及承载能力的两层不确定性优化模型,依据不确定性优化模型的设计变量、约束条件、目标函数对复合材料受压蒙皮稳定性及承载能力的不确定性优化模型进行求解,得到复合材料受压蒙皮稳定性及承载能力的优化方案。
Description
技术领域
本发明主要是适用于复合材料蒙皮结构稳定性及承载能力的不确定性优化及方案的制定,具体的涉及一种适合复合材料受压蒙皮稳定性及承载能力的不确定性优化模型和不确定性传播分析方法。
背景技术
复合材料蒙皮由高强度、低密度的纤维材料与基体组成,与传统的金属蒙皮材料相比具有更高的比强度和比模量等优点,因此在现代航空航天、汽车等工业领域已得到广泛关注和应用。有资料显示,美国波音公司设计生产的B787全机结构复合材料覆盖率达到50%,其机翼部件复合材料覆盖率更是高达90%以上。
传统的复合材料蒙皮优化设计方法是通过选择确定性的材料参数及结构参数,建立优化模型,寻求满足约束条件的,并使性能指标最优的确定性优化方案。然而,复合材料蒙皮结构中存在着大量不确定的优化参数,如复合材料的几何特征、材料特征、外部属性等。传统的复合材料蒙皮优化设计方法往往没有考虑不确定因素的影响,其数学上的最优方案对未来的情况而言可能并非最优,甚至可能需要进行大量的补偿投入,造成损失和浪费。为此考虑了复合材料分散性的不确定性优化设计方法受到学术界和工程界的高度重视。
为充分利用不确定性参数的信息寻求性能指标最优的不确定性优化方案,国内外学者做了大量工作,主要集中在两个方面:一方面是将不确定性参数定量化;另一方面是改善不确定传播分析方法,主要是发展随机有限单元法、区间有限单元法,如Taylor展开随机有限单元法,Taylor展开区间有限单元法等。这些研究在一定程度上丰富了不确定性优化设计方法,但是仅仅利用了不确定性参数的均值、方差、中心值和区间半径等信息,采用Taylor展开法区间运算具有扩张性质,单纯基于区间自然扩张原理的区间运算求解含有区间数的不确定性问题,常会因得到过宽区间界限而失去实际应用意义。
由于以往的不确定性传播分析方法未能充分利用不确定性参数的更多信息(例如除不确定性参数的均值、方差、中心值和区间半径之外的信息),导致响应估计的区间界限过宽而过于保守,使得基于此类不确定性传播分析方法的不确定性优化设计不能获得很好的优化方案。
发明内容
本发明的目的在于提供一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法。为实现上述目的,本发明提供的复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其具体实现步骤是:
第一步:将复合材料受压蒙皮中的区间不确定性参数向量处理成为一元区间不确定性参数向量,区间不确定性参数向量表示为:
式中,为区间不确定性参数向量α的中值;为区间不确定性参数向量α的半径;e∈IIn,IIn定义为所有元素值包含在[-1,1]内的n维向量集合。符号“×”定义为两个向量各对应元素相乘的算子,乘积仍然为维数相同的向量。
一元区间不确定性参数向量表示为:
αi=αc+Δα×Xi
式中,Xi=(0,…,x,…,0)T,角标i代表αi中第i个分量为区间不确定性参数。由此可见,一个n维的区间不确定性参数向量通过处理变为n个一元区间不确定性参数向量。
第二步:在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集,配点原则是采用高斯积分点在区间内配点,区间内的Gauss积分点记为xk,表示为:
式中,xk为区间内配置的第k个高斯积分点,q为区间内配点个数。
第三步:建立复合材料受压蒙皮稳定性及承载能力的不确定性优化模型,其基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(F,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮的稳定性及承载能力为优化目标maxFc=f(x1,x2...,xm,X1,X2,…,Xn),xi为设计变量;Xi为区间不确定性参数。
第四步:将步骤三建立的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型分解为两个子模型,第一个子模型的基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮结构的屈曲临界载荷为优化目标,xi为复合材料受压蒙皮铺层角,即设计变量;Xi为复合材料受压蒙皮区间不确定性参数;和分别为第i个设计变量的下界和上界;m为设计变量的个数。
第二个子模型的基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(P,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮结构的压力临界载荷为优化目标,xi为复合材料受压蒙皮铺层角,即设计变量;Xi为复合材料受压蒙皮区间不确定性参数;εmax为复合材料受压蒙皮结构中最大应变或最大应力;ε0为复合材料受压蒙皮结构的最大许用应变或最大许用应力;和分别为第i个设计变量的下界和上界;m为设计变量的个数。
第五步:求解第一个子模型,得到复合材料受压蒙皮结构屈曲临界载荷上、下界值并给出复合材料受压蒙皮结构屈曲临界载荷的中值和复合材料受压蒙皮结构屈曲临界载荷不确定区间半径对复合材料受压蒙皮结构屈曲临界载荷求解时采用最佳平方逼近多项式逼近结构响应函数,引入r阶第一类Chebyshev多项式,其正交多项式系{Tj(x)}和最佳平方逼近函数为:
Tj(x)=cos(jarccosx),-1≤x≤1
式中,j(0≤j≤r)为非负整数,aj为逼近函数展开式系数。
第六步:求解第二个子模型,得到满足许用应变或许用应力约束下复合材料受压蒙皮结构最大压力载荷上、下界值(pL,pU),并给出复合材料受压蒙皮结构最大压力载荷的中值和复合材料受压蒙皮结构最大压力载荷不确定区间半径
第七步:将步骤五和步骤六求得的解相比较,取两个中值中最小的载荷作为复合材料受压蒙皮稳定性及承载能力的性能指标,性能指标Fc是复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的优化目标;
第八步:采用序列二次规划法寻找步骤三建立的基于配点型区间方法的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的最优解,得到复合材料受压蒙皮稳定性及承载能力的性能指标Fc最大时复合材料的铺层方案,此方案为复合材料受压蒙皮稳定性及承载能力的最终优化方案。
本发明与现有技术相比的优点在于:本发明提供了一种适合复合材料受压蒙皮稳定性及承载能力不确定性优化的新方法,与随机不确定性优化方法和基于Taylor展开的区间不确定性优化方法相比,该方法不需要知道复合材料受压蒙皮结构不确定性参数的概率分布形式,只要知道复合材料受压蒙皮结构不确定性参数的上下界限就可以方便的求解并得到复合材料受压蒙皮结构稳定性及承载能力的估计区间界限。通过在复合材料受压蒙皮结构不确定性参数区间内配置高斯点,避免了Taylor方法带来的区间信息丢失,充分利用了除Taylor展开点以外的其它区间内的信息,有效的缩小了复合材料受压蒙皮稳定性及承载能力估计区间界限,得到更“紧”的区间界限,在工程实际使用中更为方便和可靠,对复合材料蒙皮结构稳定性及承载能力的不确定性优化及方案的制定有很好的直接应用价值。
附图说明
图1是本发明不确定性优化求解流程;
图2是本发明配点型区间分析流程;
图3是本发明实例中复合材料受压蒙皮结构三维模型;
图4是本发明实例中复合材料受压蒙皮结构试验件部分试件图纸;
其中a为试验件主视图;b为试验件A向剖视图;c为试验件筋条放大图;
图5是本发明实例中复合材料受压蒙皮结构试验件应变片贴片位置部分图纸;
图6是本发明实例中复合材料受压蒙皮结构各层合板优化前后铺层角示意图;
其中(a)是筋条下翼使用的23层层合板,(b)是筋条上翼和筋条腹板使用的32层层合板,(c)是壁板使用的35层层合板。
具体实施方式
如图1和图2所示,本发明提出了一种基于配点型区间分析的复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其具体实现步骤是:
(1)将复合材料受压蒙皮中的区间不确定性参数向量处理成为一元区间不确定性参数向量,区间不确定性参数向量表示为:
式中,为区间不确定性参数向量α的中值;为区间不确定性参数向量α的半径;e∈IIn,IIn定义为所有元素值包含在[-1,1]内的n维向量集合。符号“×”定义为两个向量各对应元素相乘的算子,乘积仍然为维数相同的向量。
取e中的第i(1≤i≤n)个元素为x,其它元素为0,记为:
式中,x∈[-1,1]。由式(1)和式(2)得到一元区间不确定性参数向量为:
αi=αc+Δα×Xi(3)
式中,αi为一元区间不确定性参数向量,角标i代表αi中第i个分量为区间不确定性参数。由此可见,一个n维的区间不确定性参数向量通过处理变为n个一元区间不确定性参数向量。
(2)在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集,配点原则是采用高斯积分点在区间内配点,在[-1,1]上配置q个Gauss积分点,记为xk。Gauss积分点xk,k=1,2,…,q为Tq(x)的零点,表示为:
式中,xk为区间内配置的第k个高斯积分点,q为区间内配点个数。
(3)建立复合材料受压蒙皮稳定性及承载能力的不确定性优化模型,其基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(F,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮的稳定性及承载能力为优化目标maxFc=f(x1,x2...,xm,X1,X2,…,Xn),xi为设计变量;Xi为区间不确定性参数。
(4)将步骤(3)建立的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型分解为两个子模型,第一个子模型的基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮结构的屈曲临界载荷为优化目标,xi为复合材料受压蒙皮铺层角,即设计变量;Xi为复合材料受压蒙皮区间不确定性参数;和分别为第i个设计变量的下界和上界;m为设计变量的个数。
第二个子模型的基本形式为:
式中,以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(P,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮结构的压力临界载荷为优化目标,xi为复合材料受压蒙皮铺层角,即设计变量;Xi为复合材料受压蒙皮区间不确定性参数;εmax为复合材料受压蒙皮结构中最大应变或最大应力;ε0为复合材料受压蒙皮结构的最大许用应变或最大许用应力;和分别为第i个设计变量的下界和上界;m为设计变量的个数。
(5)求解第一个子模型,得到复合材料受压蒙皮结构屈曲临界载荷上、下界值并给出复合材料受压蒙皮结构屈曲临界载荷的中值和复合材料受压蒙皮结构屈曲临界载荷不确定区间半径对复合材料受压蒙皮结构屈曲临界载荷求解时采用最佳平方逼近多项式逼近结构响应函数。引入r阶第一类Chebyshev多项式,其正交多项式系{Tj(x)}和最佳平方逼近函数为:
Tj(x)=cos(jarccosx),-1≤x≤1(8)
式中,j(0≤j≤r)为非负整数,aj为逼近函数展开式系数。
由Gauss积分点求出多项式系数并代入式(9),可进一步得到:
简记为:
式中:
T(x)=[1T1(x)T2(x)…Tr(x)]T(14)
先考虑如何求解的最值,对式(11)关于x求导并令导数为零,得:
求解式(15)的根,并联合和根据连续函数在闭区间上的最值定理,可得一元逼近函数的最小值点和最大值点,分别记为和
重复以上过程,直到i遍历完1~n时,得到具有n个元素的最值点向量,记为:
由式(2)可知,Xmin和Xmax是n维向量集合IIn中的向量。进而通过式(1)、式(3)、式(10)及式(16)可以得到结构响应的近似区间估计为按下式计算:
式中,和分别为第i个一元区间不确定性参数向量在区间[-1,1]内的最小值点和最大值点;Xmin和Xmax分别是由最小值点和最大值点构成的最值点向量;αc为区间不确定性参数向量的中值向量;Δα为区间不确定性参数向量的区间半径向量。
(6)求解第二个子模型,得到满足许用应变或许用应力约束下复合材料受压蒙皮结构最大压力载荷上、下界值(pL,pU),并给出复合材料受压蒙皮结构最大压力载荷的中值和复合材料受压蒙皮结构最大压力载荷不确定区间半径对复合材料受压蒙皮结构静力受压分析时,同样采用最佳平方逼近多项式逼近结构响应函数;
(7)将步骤(5)和步骤(6)求得的解相比较,取两个中值中最小的载荷作为复合材料受压蒙皮稳定性及承载能力的性能指标,性能指标Fc是复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的优化目标;
(8)采用序列二次规划法寻找步骤三建立的基于配点型区间方法的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的最优解,得到复合材料受压蒙皮稳定性及承载能力的性能指标Fc最大时复合材料的铺层方案,此方案为复合材料受压蒙皮稳定性及承载能力的最终优化方案。
实施例:
1.结构参数及模型分析介绍
为了更充分的了解该发明的特点及其对工程实际的适用性,本发明采用一段飞行器复合材料受压蒙皮结构件进行不确定性优化设计和试验验证。飞行器复合材料受压蒙皮结构为复合材料加筋板见图3,蒙皮结构由4部分构成:筋条上翼、筋条下翼、筋条腹板和壁板。各部分的复合材料铺层情况见表1和图6。复合材料单层板的材料参数见表2。由于蒙皮结构在实际使用中主要承受沿筋条方向的载荷,对飞行器复合材料受压蒙皮结构稳定性及承载能力进行有限元分析时,将蒙皮结构筋条方向的一端施加位移约束,在另一端施加沿筋条方向的压力,对结构依次进行两次分析:第一次是屈曲分析,得到复合材料受压蒙皮的屈曲临界载荷;第二次是受压静力分析,得到复合材料受压蒙皮所能承受的最大压力载荷,其条件是在受压过程中,蒙皮结构的最大应变不能超过3000应变。将得到的屈曲临界载荷和压力载荷进行比较,取二者较小的载荷作为衡量复合材料受压蒙皮稳定性及承载能力。
2.优化模型框架
飞行器复合材料受压蒙皮稳定性及承载能力的不确定性优化模型,其基本形式为:
将上述不确定性优化模型分解为两个子模型,第一个子模型的基本形式为:
该子模型以配点型区间方法对蒙皮结构进行不确定性屈曲分析,得到屈曲临界载荷的上下界值,再以屈曲临界载荷的中心值最大且半径最小为目标函数进行优化。
第二个子模型的基本形式为:
该子模型以配点型区间方法对蒙皮结构进行不确定性受压静力分析,得到满足约束条件的最大压力载荷上下界值,再以压力载荷的中心值最大为目标函数进行优化。
其中参数意义如说明书中所述。
3.优化结果与分析
基于配点型区间方法的不确定性优化模型能够解决以区间形式存在的不确定性,从而为含有区间不确定性参数的复合材料蒙皮结构设计提供最优的规划方案。表3给出了本发明的优化结果及其与其它设计方法的比较。初始设计是工程中常见的复合材料铺层角度,一般取30°、45°、60°等。传统设计采用确定性的优化设计方法,但未考虑不确定因素的影响。从几种优化结果的比较中得出:采用传统设计与本发明设计设计的复合材料蒙皮结构稳定性及承载能力优于传统设计的设计值,这是因为后两者都采用了优化设计的思想。传统设计未考虑不确定因素的影响,导致结构在实际情况中存在小于设计值的潜在危险。本发明设计考虑了复合材料参数的区间不确定性,得到的蒙皮稳定性及承载能力优于初始设计,略次于传统设计,但是其设计中不存在潜在危险,从工程实际出发,考虑了不确定因素的本发明设计方法要优于前两种设计方法。
4.试验结果比较
根据本发明不确定性优化设计方法得到的蒙皮铺层角度,生产了复合材料蒙皮结构试验件,试验件部分图纸见图4。将复合材料蒙皮结构试验件置于力学试验机上进行压力试验,部分应变片贴片位置见图5。表4给出了试验时应变片测量的部分试验数据。经过数据处理得到了试验测得的复合材料蒙皮稳定性及承载能力。表3给出了复合材料蒙皮结构试验件的稳定性及承载能力试验结果和数值计算结果。从表中可知试验结果与本发明设计计算结果吻合较好,说明了本发明设计方法的准确性,同时也体现了本发明设计方法对复合材料蒙皮结构稳定性及承载能力优化设计及方案制定具有很好的直接应用价值。
表1
表2
表3
表4
本发明未详细阐述部分属于本领域技术人员的公知技术。
Claims (6)
1.一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于实现步骤如下:
第一步:将复合材料受压蒙皮中的区间不确定性参数向量处理成为一元区间不确定性参数向量,区间不确定性参数向量是以区间形式表示的不确定性参数向量,一个n维的区间不确定性参数向量变为n个一元区间不确定性参数向量,n维的区间不确定性参数向量中每一维都是区间不确定性参数,一元区间不确定性参数向量中只有其中一维是区间不确定性参数,其它n-1维是确定性参数;
第二步:在一元区间不确定性参数向量的区间内配点,生成区间不确定性参数向量的区间配点集;
第三步:以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用应变或许用应力为约束条件,以复合材料受压蒙皮稳定性及承载能力为优化目标,建立基于配点型区间方法的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型;
第四步:将步骤三建立的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型分解为两个子模型:第一个子模型为复合材料受压蒙皮结构屈曲分析模型;第二个子模型为复合材料受压蒙皮结构静强度受压分析模型,根据区间配点方案求解响应的界值估计;
第五步:求解第一个子模型,得到复合材料受压蒙皮结构屈曲临界载荷上、下界值并给出复合材料受压蒙皮结构屈曲临界载荷的中值和复合材料受压蒙皮结构屈曲临界载荷不确定区间半径
第六步:求解第二个子模型,得到满足许用应变或许用应力约束下复合材料受压蒙皮结构最大压力载荷上、下界值(pL,pU),并给出复合材料受压蒙皮结构最大压力载荷的中值和复合材料受压蒙皮结构最大压力载荷不确定区间半径
第七步:将步骤五和步骤六求得的解相比较,取两个中值中最小的载荷作为复合材料受压蒙皮稳定性及承载能力的性能指标,性能指标Fc是复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的优化目标;
第八步:采用序列二次规划法寻找步骤三建立的基于配点型区间方法的复合材料受压蒙皮稳定性及承载能力的不确定性优化模型的最优解,得到复合材料受压蒙皮稳定性及承载能力的性能指标Fc最大时复合材料的铺层方案,此方案为复合材料受压蒙皮稳定性及承载能力的最终优化方案。
2.根据权利要求1所述的一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于:步骤一中区间不确定性参数向量表示为:
式中,为区间不确定性参数向量α的中值;为区间不确定性参数向量α的半径;e∈IIn,IIn定义为所有元素值包含在[-1,1]内的n维向量集合;符号“×”定义为两个向量各对应元素相乘的算子,乘积仍然为维数相同的向量;
一元区间不确定性参数向量表示为:
αi=αc+△α×Xi
式中,Xi=(0,…,x,…,0)T,角标i代表αi中第i个分量为区间不确定性参数,由此可见,一个n维的区间不确定性参数向量通过处理变为n个一元区间不确定性参数向量。
3.根据权利要求1所述的一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于:所述步骤二中的在一元区间不确定性参数向量的区间内配点的原则是采用高斯积分点在区间内配点,区间内的Gauss积分点记为xk,表示为:
式中,xk为区间内配置的第k个高斯积分点,q为区间内配点个数。
4.根据权利要求1所述的一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于:所述步骤三中建立的不确定性优化模型是以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(F,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮的稳定性及承载能力为优化目标maxFc=f(x1,x2...,xm,X1,X2,…,Xn),不确定性优化模型为:
式中,xi为设计变量;Xi为区间不确定性参数,εmax为复合材料受压蒙皮结构中最大应变或最大应力;ε0为复合材料受压蒙皮结构的最大许用应变或最大许用应力,m为设计变量的个数。
5.根据权利要求1所述的一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于:所述步骤四中分解的两个优化子模型,第一个优化子模型是以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮结构的屈曲临界载荷为优化目标,建立的无条件约束优化模型为
式中,xi为复合材料受压蒙皮铺层角,即设计变量;Xi为复合材料受压蒙皮区间不确定性参数;和分别为第i个设计变量的下界和上界;m为设计变量的个数;
第二个优化子模型是以复合材料受压蒙皮铺层角为设计变量,以复合材料受压蒙皮的许用压应变或许用压应力为约束条件εmax(P,x1,x2,…,xm,X1,X2,…,Xn)≤ε0,以复合材料受压蒙皮结构的压力临界载荷为优化目标,建立的优化模型为
式中,εmax为复合材料受压蒙皮结构中最大应变或最大应力;ε0为复合材料受压蒙皮结构的最大许用应变或最大许用应力。
6.根据权利要求1所述的一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法,其特征在于:所述步骤五中对复合材料受压蒙皮结构屈曲分析模型求解时以及步骤六中对复合材料受压蒙皮结构静强度受压分析模型求解时,采用最佳平方逼近多项式逼近结构响应函数,表示为:
式中,Tj(x)=cos(jarccosx),-1≤x≤1,0≤j≤r,为正交多项式系;为配点处的结构响应;(x)为r阶第一类Chebyshev多项式,角标i表示针对第i个一元区间不确定性参数向量,采用最佳平方逼近多项式逼近结构响应函数;xk为区间内配置的第k个高斯积分点,q为区间内配点个数;
求解x∈[-1,1]的最小值点和最大值点,分别记为和重复以上过程,直到i遍历完1~n时,就能得到具有n个元素的最值点向量,记为和将Xmin和Xmax分别代入结构响应函数中,得到结构响应的近似区间估计为按下式计算
式中,和分别为第i个一元区间不确定性参数向量在区间[-1,1]内的最小值点和最大值点;Xmin和Xmax分别是由最小值点和最大值点构成的最值点向量;αc为区间不确定性参数向量的中值向量;△α为区间不确定性参数向量的区间半径向量;响应的中值由下式给出
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310207200.XA CN103324782B (zh) | 2013-05-29 | 2013-05-29 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310207200.XA CN103324782B (zh) | 2013-05-29 | 2013-05-29 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103324782A CN103324782A (zh) | 2013-09-25 |
CN103324782B true CN103324782B (zh) | 2016-01-20 |
Family
ID=49193524
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310207200.XA Expired - Fee Related CN103324782B (zh) | 2013-05-29 | 2013-05-29 | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103324782B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678786B (zh) * | 2013-11-29 | 2016-06-29 | 北京宇航系统工程研究所 | 一种运载火箭蒙皮桁条结构的有限元分析方法 |
CN105760619B (zh) * | 2016-03-10 | 2018-06-29 | 北京航空航天大学 | 一种考虑疲劳可靠性的复合材料铺层优化设计方法 |
CN105808884B (zh) * | 2016-03-30 | 2021-05-25 | 北京航空航天大学 | 一种基于分形理论的有界不确定性平面裂纹应力强度因子上下界的预测方法 |
CN106021654B (zh) * | 2016-05-09 | 2018-04-03 | 北京航空航天大学 | 一种考虑铺层相容性的复合材料铺层方案库建库方法 |
CN105893716B (zh) * | 2016-06-02 | 2018-06-29 | 北京航空航天大学 | 一种基于分形理论的结构断裂非概率可靠性分析方法 |
CN105956336B (zh) * | 2016-06-12 | 2019-01-18 | 北京航空航天大学 | 一种面向静动态混合不确定性的拟建结构非概率可靠性优化设计方法 |
CN107451335A (zh) * | 2017-07-06 | 2017-12-08 | 华东建筑设计研究院有限公司 | 一种引入蒙皮效应的铝合金网格结构整体稳定分析方法 |
CN108846149B (zh) * | 2018-04-20 | 2020-02-07 | 北京航空航天大学 | 一种基于多源不确定性的结构分布式动态载荷识别的方法 |
CN110008530B (zh) * | 2019-03-15 | 2020-04-24 | 东南大学 | 一种空间柔性复合材料分布式概率建模方法 |
CN111948044B (zh) * | 2020-07-24 | 2021-11-19 | 中国飞机强度研究所 | 一种拉压剪复合载荷作用的加筋壁板失效预测方法 |
CN112434381B (zh) * | 2020-12-16 | 2024-06-04 | 大连理工大学 | 一种重型运载火箭栅格加筋结构在缩比模型中的刚度等效方法 |
CN114291250B (zh) * | 2021-12-20 | 2023-11-03 | 北京机电工程研究所 | 一种剪切变后掠翼面及其设计方法 |
CN114486518B (zh) * | 2021-12-31 | 2024-06-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种结构复合材料选用效果评估方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101939217A (zh) * | 2008-02-08 | 2011-01-05 | 空中客车营运有限公司 | 用于制造纤维复合材料构件的方法、纤维复合材料构件以及飞机的纤维复合材料机身部件 |
CN102060101A (zh) * | 2010-12-21 | 2011-05-18 | 南京航空航天大学 | 一种用于变体机翼的蒙皮 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ES2384349B1 (es) * | 2009-12-30 | 2013-05-16 | Airbus Operations, S.L. | Cuaderna de fuselaje de aeronave en material compuesto con costillas estabilizadoras. |
-
2013
- 2013-05-29 CN CN201310207200.XA patent/CN103324782B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101939217A (zh) * | 2008-02-08 | 2011-01-05 | 空中客车营运有限公司 | 用于制造纤维复合材料构件的方法、纤维复合材料构件以及飞机的纤维复合材料机身部件 |
CN102060101A (zh) * | 2010-12-21 | 2011-05-18 | 南京航空航天大学 | 一种用于变体机翼的蒙皮 |
Non-Patent Citations (2)
Title |
---|
《Eigenvalue bounds of structures with uncertain-but-bounded parameters》;Zhiping Qiu 等;《JOURNAL OF SOUND AND VIBRATION》;20041014;第297-312页 * |
配点型区间有限元法;邱志平;《力学学报》;20110531;第43卷(第3期);第496-504页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103324782A (zh) | 2013-09-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103324782B (zh) | 一种复合材料受压蒙皮稳定性及承载能力的不确定性优化方法 | |
Bazilevs et al. | Fluid–structure interaction modeling for fatigue-damage prediction in full-scale wind-turbine blades | |
Vescovini et al. | A variable-kinematic model for variable stiffness plates: Vibration and buckling analysis | |
Karama et al. | A new theory for laminated composite plates | |
Egan et al. | Modelling bearing failure in countersunk composite joints under quasi-static loading using 3D explicit finite element analysis | |
WO2019047529A1 (zh) | 一种纤维增强复合材料动态剪切本构模型的构建方法 | |
Vescovini et al. | Efficient post-buckling analysis of variable-stiffness plates using a perturbation approach | |
Nachtane et al. | Numerical investigation of damage progressive in composite tidal turbine for renewable marine energy | |
Yang et al. | Prediction of post-buckling and ultimate compressive strength of composite plates by semi-analytical methods | |
Bekemeyer et al. | Rapid gust response simulation of large civil aircraft using computational fluid dynamics | |
Ye et al. | Higher order semi-analytical solution for bending of angle-ply composite laminated cylindrical shells based on three-dimensional theory of elasticity | |
CN109507040B (zh) | 一种蜂窝夹层结构面板压缩应力评估方法 | |
Loseille et al. | Anisotropic adaptive simulations in aerodynamics | |
Cheng et al. | Modelling pre-fatigue, low-velocity impact and post-impact fatigue behaviours of composite helicopter tail structures under multipoint coordinated loading spectrum | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
McElroy | An enriched shell element for delamination simulation in composite laminates | |
Pandya et al. | CFD Computations for a Generic High-Lift Configuration Using TetrUSS | |
Chitale et al. | Boundary layer adaptivity for incompressible turbulent flows | |
Xiaohui et al. | Buckling of soft-core sandwich plates with angle-ply face sheets by means of a C 0 finite element formulation | |
Fang et al. | X-FEM co-simulation of delamination and matrix cracking in fiber metal laminated structures under fatigue loading | |
Fallah et al. | Computational modelling of dynamic delamination in morphing composite blades and wings | |
Mohammad Panahiha et al. | Buckling analysis of composite grid-stiffened cylindrical shells using a generalized equivalent single layer theory | |
Xue et al. | The application of reliability based optimization of tophat stiffened composite panels subject to bi-directional buckling loads | |
Mehmood et al. | Failure analysis of single and double bolted joints for composite laminates | |
Chen et al. | Modeling vibration behavior of delaminated composite laminates using meshfree method in Hamilton system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160120 |