chapter7- machine elements-2.ppt

等效正态分布方法 将非标准正态分布在设计点处转换成一个等效的正态分布，然后利用一次二阶矩的迭代法求解。 适用于任何分布类型。 也适用于极限状态方程中有多个随机变量的情况。 例 题 一拉杆受外力作用，若外力的均值和标准差分别为 杆横截面积的均值和标准差分别为 求拉应力的均值和标准差。 解： 用蒙特卡洛法求应力、强度干涉解 从应力分布中随机选取一个应力值样本，从强度分布中随机选取一个强度值样本，将两个样本相比较，如果应力大于强度，则零件失效；反之，则零件可靠。每一次随机模拟相当于对一个随机抽取的零件进行一次试验，通过大量重复的随机抽样及比较，可得到零件的总失效数，从而可以求得零件的失效概率和可靠度近似值。抽样系数越多，模拟的精度越高。 统计试验法 随机模拟法 Monte Carlo 仿真算例 工作寿命 lgn 均值 MPa 标准差 MPa 4.30 4.40 4.50 4.70 4.70 4.80 4.90 5.00 5.10 5.20 5.30 685 661 638 617 596 578 562 546 530 514 499 14.0 13.1 12.6 13.3 13.0 12.3 13.0 13.8 14.4 14.8 15.0 钢轴受弯矩作用，其最大应力服从正态分布，分布参数为 钢轴的强度也服从正态分布，数据列于表中。求钢轴运转105次时，此轴的可靠度。 10 Simulations R S Y 543.4234 361.0918 182.3316 556.0159 310.0448 245.9711 537.8812 384.1888 153.6925 576.1280 390.9098 185.2183 544.1177 331.5361 212.5817 547.5723 428.3039 119.2684 560.7214 428.2314 132.4900 546.8181 377.4420 169.3761 544.6801 392.5499 152.1301 534.5136 386.2301 148.2835 100 and 10000 simulations Matlab 代码 clear all close all % Define input variables S_mean=379; S_std=41.4; R_mean=546; R_std=13.8; %Step 1 -Sampling on random variables randn(state,0); N=10000; S_sample=normrnd(S_mean,S_std,N,1); R_sample=normrnd(R_mean,R_std,N,1); % Step 2 - Experimentation g=R_sample-S_sample; Matlab 代码 % Step 3 - statistical analysis Nf=sum(g0); Pf=Nf/N; R=1-Pf; % plot plot(S_sample,R_sample,.); xlabel(S); ylabel(R); text(300,400,°2è???); text(400,350,ê§D§??); text(400,300,[Number of failures=,num2str(Nf)]); text(400,250,[Number of simulations=,num2str(N)]); text(400,200,[probability of failures=,num2str(Pf)]); hold on ezplot(R,[200,600]); 用等效正态分布方法计算零件的可靠度 功能函数（极限状态方程） 功能函数中含有非正态变量时，可以用一个与原函数等效的正态分布替代。要求满足以下2个条件： 原函数值F(xi*)与当量正态函数值F’(xi*)相等 原概率密度值f(xi*)与当量正态分布概率密度值f’(xi*)相等 原分布 FXi(xi*),fXi(xi*) 等效正态分布F’Xi(xi*),f’Xi(xi*) O 等效正态分布的均值和标准差分别为： 之后求解可靠指标的过程与前面随机变量服从正态分布的过程相同。 原密度函数 原分布函数 主 要 内 容 机械可靠性设计概述 应力——强度干涉理论 机械零件静强度可靠性设计 机械零件疲劳强度可靠性设计 机械零件静强度可靠性设计 主要是把应力、强度分布和可靠度在概率的意义下联合起来，构成一种设计计算的依据。 基本步骤： 根据零部件的功能、复杂程度、重要程度、使用条件、生产难易程度、相似产品失效的历史数据以及研制水平等确定可靠性指标。 明确零部件失效模式，如屈服、失稳、断裂、过量变