飞机结构静强度计算

飞机强度计算方法飞机结构静强度计算3.1飞机结构静强度与结构可靠性计算机翼和机身的强度估算结构可靠性返回

结构静强度计算方法有多种，但结构静强度计算仍是结构设计的基础，主要体现在下列三个阶段。

飞机总体设计中的结构布局和结构形式的确定对结构连接部位、开口区、复合材料铺层等细节进行设计计算结构静强度校核阶段结构有限元分析结构优化设计1、机翼和机身的强度估算

一般采用有限元方法，但在结构初步设计和结构强度分析时，常采用薄壁结构力学方法。具体的公式和简化方法可参见设计手册，不一一讲解。2、结构有限元分析

MSC/NASTRAN3、结构优化设计4、结构可靠性4.1结构可靠性概念

可靠性是指结构在规定条件下和规定时间内，完成规定功能的能力。结构可靠性定义的要素是三个“规定”(“规定条件”、“规定时间”、“规定功能”）结构在规定的条件下和规定的时间内，完成规定功能的概率称为可靠度。结构在规定的条件下和规定的时间内，丧失规定功能的概率称为不可靠度或失效概率。作为飞机结构的可靠性问题，从定义上可以理解为：“结构在规定的使用载荷/环境工作下及规定的时间内，为防止各种失效或有碍正常工作功能的损伤，应保持其必要的强刚度、抗疲劳断裂以及耐久性能力。”可靠度则应是这用能力的概率度量。例如：结构静强度可靠性是指结构元件或结构系统的强度大于工作应力的概率；结构安全寿命可靠性是指结构的裂纹形成寿命小于使用寿命的概率；结构损伤容限可靠性则一方面指结构剩余强度大于工作应力的概率，另一方面指结构在规定的未修使用期内，裂纹扩展小于裂纹容限的概率。其它可靠度度量方法：

结构的失效概率F(t)，指结构在t时刻之前破坏的概率；

失效率λ(t)，指在t时刻以前未发生破坏的条件下，在t时刻的条件破坏概率密度；

平均无故障时间MTTF(MeanTimeToFailure)，指从开始使用到发生故障的工作时间的期望值。4.1结构可靠性概念4.2结构安全余量方程进行结构元件可靠性分析时，需要建立起元件设计变量与元件能力表征量间的分析关系，这类似于确定性分析设计中的工程破坏判据，但可靠性分析是建立在随机变量的分析基础之上。这个概率型的联系设计变量与结构元件固有性能表征量间的破坏判据，通常称为元件的安全余量方程（功能函数）。讨论结构元件的静强度可靠性时，可初步认为只有两个随机变量，即元件的强度R和元件的内力S。元件的强度由于材料的强度特性、元件尺寸等不确定因素呈随机性；而元件所承受的内力，由于作用载荷的随机性以及元件尺寸在结构系统中所处的位置等不确定因素显然是随机变量。如果元件能够承受载荷，则安全余量方程为4.2结构安全余量方程可靠度定义为元件能可靠承载的概率，可以表示为则元件的失效概率可以表示为4.3应力强度干涉模型R，SfOfSfRμSμR干涉区可靠度一般表达式4.3应力强度干涉模型应当指出应力强度干涉模型揭示了概率设计的本质。从干涉模型可以看到，就统计数据观点而言，任何一个设计通常存在着失效概率，即可靠度小于1，而我们设计能够做到的仅仅是将失效概率限制在一个可以接受的限度之内，该观点在常规设计的安全系数法中是不明确的。可靠性设计的这一重要特征客观地反映了产品设计和运行的真实情况，同时还定量地给出了产品在使用中的失效概率或可靠度，因而收到重视与发展。4.4可靠性指标当应力和强度均为正态分布时，有可靠性指标由此可以看出，在分析线性安全余量方程且变量间服从正态分布的可靠性概率时，可靠性指标与失效概率一样，可表征可靠性程度。对于非线性安全余量、变量不服从正态分布的情况，可将非线性安全余量在设计验算点处作近似线性展开，并将非正态分布变量转换成正态分布变量。因此，可靠性指标β在可靠性分析中具有重要的实际意义。4.4可靠性指标β00.51.01.52.02.53.04.05.0Pf0.50.30850.15870.06680.02280.00620.00143.27×10-53×10-7Pr0.50.69150.84130.97720.97720.99380.99860.99996730.99999974.4可靠性指标例如某构件强度和所受应力均服从正态分布，具体数据如下：则以上讨论的为线性安全余量，且变量服从正态分布。4.5可靠性指标（均值一次二阶矩法）以上讨论的为线性安全余量，当安全余量为非线性时，将安全余量方程在各变量均值点处进行泰勒展开，仅取展开项中的线性项（一次项），忽略高次项，则有这样，安全余量成为线性函数，当各变量相互独立时，其均值和方差如下则可靠性指标为4.5可靠性指标（均值一次二阶矩法）算例：某受拉铝杆，已知材料强度均值为μR=360N/mm2，标准差为σR=20N/mm2；杆的直径d的均值μd=10mm，标准差为σd=0.04mm；所受拉力P的均值μP=20000N，标准差σP=600N。求该拉杆的可靠性指标。解：安全余量为则4.5可靠性指标（均值一次二阶矩法）在上例中若安全余量取为采用同样方法求得的可靠性指标为从计算结果可以看出，取不同的安全余量，用均值一次二阶矩方法求得结果是不同的，因此需要改进。最常用的方法为改进的一次二阶矩方法（验算点法、JC法）。但由于一次二阶矩方法有计算方便简单的特点，应用较广泛，对于初步估算较好。4.6可靠性指标（验算点法、JC法）a）随机变量为正态分布情况

Hasofer和Lind建议根据临界破坏面而不是安全余量方程定义失效模式的可靠度指标

。对于同一物理问题，根据H-L算法计算得到的可靠度指标

，不会由于选择不同形式的等价安全余量方程而发生变化。H-L方法的计算程序为将随机变量

进行正则化处理相应的可靠度指标定义为4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）按泰勒级数展开并取一次项有上式两端同除以得4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）b）随机变量为非正态分布情况4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法）4.6可靠性指标（验算点法、JC法、算例）取极限状态方程为4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）M=pi*d.^2.*s/4-p;p=100*10^3;us=290*10^6;ds=25*10^6;ud=3*10^(-2);dd=3*10^(-3);fori=1:10%;dM_d=pi*d*s/2;%偏导数dM_s=pi*d^2/4;K=sqrt((dM_d*dd)^2+(dM_s*ds)^2);a_d=dM_d*dd/K;%敏度a_s=dM_s*ds/K;g=pi*d^2*s/4-p;uM=g+dM_d*(ud-d)+dM_s*(us-s);%在设计点处的均值dM=a_d*dM_d*dd+a_s*dM_s*ds;%在设计点处的标准差bat(i)=uM/dM;%可靠性指标d1(i)=ud-bat(i)*a_d*dd;s1(i)=us-bat(i)*a_s*ds;g1(i)=pi*d1(i)^2*s1(i)/4-p;d=d1(i);%下一个迭代点s=s1(i);end4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.6可靠性指标（验算点法、JC法、算例）4.7可靠度与安全系数传统意义上的静强度设计安全系数法从概率观点上可以理解为概率变量（强度与应力）的均值化设计。那么，它所获得的可靠性效益何如呢？设一拉杆的设计安全系数为n＝1.5，设计为满应力设计，且假定强度和应力均服从正态分布，其变异系数VR＝σR/μR=0.1,VS＝σS/μS=0.15,求拉杆的强度可靠度Pr?4.7可靠度与安全系数若安全系数n＝1，则β＝0，Pr＝0.5；