




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
第三章疲劳应用统计学基础3.1疲劳数据的分散性3.2正态分布3.3威布尔分布3.4二元线性回归分析3.5S-N曲线和P-S-N曲线的拟合返回主目录1第三章疲劳应用统计学基础3.1疲劳数据的分散性1)实验:7075-T6铝R=-1,恒幅45678X=lgN105099.9Pf100
7075-T6铝合金对数疲劳寿命分布2430.117030999051Sinclair和Dolan,1953.应力水平越低,寿命越长,分散性越大。152207MPa下57件,寿命:2×106108次;240MPa下29件,寿命:7×1054×106次275MPa下34件,寿命:1×1058×105次310MPa下29件,寿命:4×1041×105次430MPa下25件,寿命:1.5×1042×104次。分散性:共174件NS(MPa)400300200104101010105678klgN2015105678+207MPa共57件寿命分布直方图100102倍对数正态分部3材质不均匀,加工质量,加载误差,试验环境等。原因:裂纹、缺口件的疲劳破坏局限在裂纹或缺口高应力局部,上述因素影响较小。光滑件寿命分散>缺口件>裂纹扩展寿命
给定应力水平下,寿命小于N的概率pf?
存活率为ps(如99%)的疲劳寿命?问题疲劳寿命常用对数正态分布、威布尔分布描述。50f(x)mx=X正态概率密度曲线3.2正态分布对数疲劳寿命lgN常常是服从正态分布的。令X=lgN,X即服从正态分布。一、正态分布的密度函数和分布函数密度函数:(-<x<)是均值;f(x)关于x=对称为标准差,是非负的。6越小,f()越大,曲线越瘦,X的分散性越小。故标准差反映X的分散性。(1)f(x)0;随机变量X取值的可能性非负。在x=处,f(x)最大,且:f(x=)=密度函数性质:(无论分布形式如何)(2);所有取值的总可能为1。0f(x)mx=X正态概率密度曲线7二、标准正态分布令,即有:注意dx=du,由密度函数变换公式可得到标准正态分布密度函数为:
(-<u<)U0-uuf(u)标准正态分布密度函数
u服从均值=0、标准差=1的正态分布。标准正态分布函数则为:9u<0或(u)<0.5,利用(-u)=1-(u)的关系求解。注意有:
(0)=0.5;
(-u)=1-(u);
Pr(a<u<b)=(b)-(a)u(u)关系,还可用近似表达式表达,如:u0(u)0.5且由,还有:
F(x)=Pr(Xx)=Pr(Uu)=(u)故求正态分布函数F(x),只需求得(u)即可。U0-uuf(u)标准正态分布密度函数(-u)F(u)1-FU0abf(u)(a)(b)-(a)10分布参数估计:设在某si下,样本含n个疲劳寿命数据xi=lgNi;破坏概率为p的对数疲劳寿命xp为:三、给定疲劳寿命下的破坏概率估计则样本均值为:样本方差s2为:标准差s是偏差(xi-)2的度量,反映分散性大小。只有(n-1)个偏差独立。up可由p确定。存活概率R=1-p。11若=95%,意味着100个样本估计的xp中,有95个小于xp(g)。即有95%的把握认为估计量小于真值。四、置信水平估计量Np=+ups,若大于真值+up,偏于危险。置信度:估计量小于真值的概率。破坏率p,置信度的对数寿命写为:单侧容限系数k:有表可查若u=0,有k=up,则xp(g)=+ups;=50%。13五、正态概率纸问题:X是否服从正态分布?已知:xF(x)关系:非线性
xu关系:线性F(x)=(u)u:一一对应能否作出xF(x)呈线性关系的坐标纸?先画x-u坐标,即若随机变量X服从正态分布,则有线性关系;再按u-(u)关系,依据u标定F(x),则线性关系不变。若X服从正态分布,F(x)-x在概率纸上呈线性。x正态概率纸up0321-1-2-3p1000.010.1150103070909999.914利用正态概率纸检验随机变量X是否服从正态分布,需xiF(xi)数据描点,由其是否线性作出判断。F(xi)是对数寿命X小于xi的概率,即破坏概率。其均秩估计量为:F(xi)=pi=i/(n+1)无论X服从何种分布,此式均适用。序号iN/10i3x=lgNiiin+1123456789101241341351381401471541601661812.09342.12712.13032.13992.14612.16732.18752.20412.22012.25770.09090.18180.27270.36360.45450.54550.63640.72730.81820.9091例3.1之xiF(xi)数据如表所列,可在正态概率纸上描点,观察是否呈线性,判断X是否服从正态分布。15分析计算框图:疲劳试验R、S给定样本数据n个N排序i破坏率F(Ni)=i/(n+1)概率纸上描点[x=lgNi,F(Ni)]是否正态分布线性?估计分布参数,s(计算或图解法)x给定破坏概率pf下的疲劳寿命?寿命N对应的pf?up=(lgN-)/s;pf=F(up)xxp=lgNp=+upsx17寿命有大于零的下限,正态分布不能反映。3.3威布尔分布Weibull1951一、密度函数和分布函数1.密度函数定义为:
(NN0)下限N0,最小寿命;
尺度参数Na,反映数据的分散性;
形状参数b;反映f(N)曲线形状。三个参数f(N)N-NN-N00ab=1b=3.5-4b=2Weibull密度函数曲线0指数Reyleigh正态分布18N=N0,F(N0)=0,即寿命小于N0的概率为零;N=Na,F(Na)=1-1/e=0.632,Na称特征寿命参数。2.分布函数:F(N)--寿命小于等于N的概率。令x=(N-N0)/(Na-N0),则有dN=(Na-N0)dx,可得:注意F(N)=F(x),故得Weibull分布函数F(N)为:19能否作出威布尔概率纸?N-F(N),非线性关系;lglg[1-F(N)]-1-lg(N-N0),线性lglg[1-F(N)]-1-F(N),一一对应二、分布参数的图解估计二个问题:N是否服从威布尔分布?如何确定其分布参数?0.90.50.1F(N)lglg[1-F(N)]-1F(N)lglg[1-F(N)]-10.010-0.521-1.339-2.360F(N)lglg[1-F(N)]-1—对应值结论:可作威布尔概率纸。若N服从威布尔分布,概率纸上lg(N-N0)-F(N)应有线性关系。威布尔概率纸lg(N-N)0456lglg[1-F(N)]-10-0.5-1.0-1.5-2.00.90.50.10.050.02F(N)21解:1)Ni排序,估计F(Ni)2)估计下限:0N0N1例3.2二组疲劳寿命数据如表。判断其是否服从威布尔分布并估计分布参数。iN(10)5N(10)5例A例B123456782.03.75.08.011.513.020.023.54.05.06.07.38.09.010.613.0F(N)=1n+10.1110.2220.3330.4440.5560.6670.7780.889.1.5120.90.50.1F(N)威布尔概率的应用0.050.02N-N0(10)6ABB'0.632B':N0=N1/2=2×105A、B:N0=022注意F(N)=0.9时,lglg[1-F(N90)]-1=0,有:Na对应的破坏概率为63.2%。3)估计分布参数Na和b。如A:N90-N0=23.5×105,Na-N0=11.5×105。有:b=0.3622/lg(23.5/11.5)=1.17A组:Na-N0=11.5×105,因为N0=0,Na=11.5×105;B‘组:Na-N0=6.8×105,N0=2×105,Na=8.8×105。.1.5120.90.50.1F(N)N-N0(10)6AB'0.63223框图:疲劳试验R、S给定样本数据n个N排序i破坏率F(Ni)=i/(n+1)取N0,概率纸上描点[x=lg(Ni-N0
),F(Ni)]是否Weibull分布线性?调整N0估计分布参数N0,b;给定破坏概率pf=F(N)下的疲劳寿命N?寿命N对应的pf?25习题:3-2,3-5b)(取N0=50千周)再见再见再见再见再见!谢谢!第一次课完
请继续第二次课
返回主目录26回归分析的主要任务是:确定回归方程的形式及回归系数;检验回归方程的可用性;利用回归方程进行预测和统计推断。设X、Y间存在着相关关系。X=x时,Y的数学期望E(Y/X=x)是x的函数,即:E(Y/X=x)=f(x)E通常未知,一般只能通过样本求其估计量:=f(x)
称为Y对X的回归方程。
~y若回归方程是线性的,有=A+Bx;常数A、B是待定的回归系数。
~y29XY0散点图分散带二、最小二乘法拟合回归方程获取数据样本(xi,yi)n对描点作散点图回归方程形式回归系数是否存在相关关系回归方程估计量与观测值yi之偏差平方和为:最小二乘法
~yQ是A、B的函数,Q最小的条件为:;由此给出方程组30正规方程组为:þýüå=å+åå=å+iiiiiiyxxBxAyxBnA2解得:
ïïþïïýü-å--å=å-ååå-å=yXxYyXxxxnyxyxnBiiiiiiiii2222-=å-ååå-åå=XBYxxnxxyxAiiiiiii22)()())(()(式中,n为样本数据点数,、分别为变量X、Y的样本均值,且=xi/n;=yi/nX_Y_X_Y_注意,均值点(、)落在回归直线上。Y_X_31三、相关系数及相关关系的检验相关系数r定义为:若令:有:LLxyxx=/r==LLLxyxxyy/BLLxxyy[/]/12å-ååå-å=xxnyxyxnBiiiiii22)(32偏差平方和为:QABxyYBXBxyiiii=+-=-+-åå()[()]__22å__=---[()()]YyBXxii2____=----+-å[()()()()]YyBXxYyBXxiiii2222___=---+-ååå()()()YyBXxBXxiii222222=---=-åå()()__YyBXxLBLiiyyxx2222上式二端除以Lyy,即得:注意:Lyy>Q>0,故相关系数r£133XY0r~1完全相关XY0r~-1完全相关XY0r~0完全不相关当时,有Q0,数据点基本在回归直线上,变量X、Y相关密切;®1r,Q,数据点越分散,相关越差;若0,X、Y完全不相关。相关系数r与B同号,r>0,则B>0,正相关;
r<0,B<0,负相关。rr相关系数的几何意义:XY00<r<1正相关XY0-1<r<0负相关34回归方程能否反映随机变量间的相关关系?是相关系数起码值,可查表。与样本容量n有关,n越大,越小。
与置信水平有关,=1-越大,越大。
是显著性水平,或纳伪概率。rarara相关性检验条件为:rra表3-4相关系数的起码值
n-20.050.01n-20.050.01n-20.050.0110.9971.00050.7540.874100.5760.708200.4230.537300.3490.449400.3040.393ra35四、利用回归方程进行统计推断对应于任一x0,y0正态分布,n大时,分布参数可估计为:
==A+Bx
=s=[Q/(n-2)]1/2=[(Lyy-B2Lxx)/(n-2)]1/2y~YX0x0y0~y~=A+Bx利用回归方程进行统计推断y~y0=+3sy~y0=-3s直线y=+ups所对应的概率为p=Pr(Yy)。如,up=3时,p=99.87%,故y落在y=+3s之下的概率为99.87%,上限。y~y~up=-3时,p=0.13%,y=-3s为0.13%的下限。y在y=+3s间的概率为p=99.74%。up=0时,p=50%,故y==A+Bx对应概率50%y~y~y~36获取样本数据(xi,yi)共n对下面通过一例题,进一步了解其分析步骤。五、
二元线性回归分析的基本方法:作散点图回归方程的形式y=A+Bx~最小二乘法确定回归系数A、B相关系数r相关性检验
rra用回归方程进行预测和统计推断37例3.3表中为某材料在R=0.1下的疲劳试验结果,试估计其S-N曲线。解:S-N曲线为SmN=C;取对数后有:lgS=lgC/m-(1/m)lgN;令y=lgS,x=lgN,回归方程可写为:y=A+Bx其中:A=lgC/m,B=-(1/m)试验数据 计算
Sa(MPa)N 199 94124166 146656141.2 298263120.2 981070
21.60638.7478117.300119.161347.1351yi=lgSai2.29892.22012.14982.0799xi=lgNi
4.97375.16635.47465.9917xi224.737726.690729.971235.9005yi25.28494.92884.62164.3260xiyi
11.434011.469711.769312.462138得到:x=xi/n=21.6063/4=5.40158;
y=2.18718Lxx=
=117.3001-21.60632/4=0.59205Lyy=0.02636;Lxy=-0.12166xxnii-åå()/22回归系数:B=Lxy/Lxx=-0.2054m=-1/B=4.87A=y-Bx=3.2974C=lg-1(Am)=1.12×1016
相关系数:r==-0.975(负相关);取=0.05时,注意n-2=2,查表知=0.950,故有r,满足相关性检验条件。LxyLLxxyy/rara39破坏率为1%时,up=-2.326,即有:y=A+Bx-2.326s=3.2362-0.2054x
破坏率为1%的S-N曲线为:
(p=0.01)回归方程可用,且S-N曲线为:
×
(p=0.50)SNa4871611210..=注意:由于数据点少(n=4),若取=0.01,则r<ra;回归方程是不可接受的。要估计破坏率为1%的S-N曲线,需先计算样本剩余标准差s。此处有:
s=[(Lyy-B2Lxx)/(n-2)]1/2=0.026340例3.4试用最小二乘法进行回归分析,估计例3.2中B组数据的分布参数。N(10)5F(N)=1n+1i12345674.05.06.07.38.09.010.60.1110.2220.3330.4440.5560.6670.778813.00.889解:1)设x=lgN正态分布有:
回归方程为:
Y=A+BXxxus=+对照可知,回归方程中变量为:
Y=x=lgN;X=u=-1[F(N)]系数:A=x
;B=s41列表计算,得到回归系数为:B=Lxy/Lxx=0.2096;A=Y-BX=0.8674;故正态分布参数:
s=B=0.2096;
x=A=0.8674YX0.6021-1.220正态分布0.69900.77820.86330.90310.95421.02531.1139-0.765-0.431-0.1400.1400.4310.7651.220若估计寿命为N=3×105时的破坏概率,则有:
=3.1%puNxsp==-=-1862FFF()(lg)(.)_f相关系数为:
r=-0.9976ra
=0.834(=0.01,n=8)422)设寿命N服从威布尔分布,有:回归方程写为:
Y=A+BX时,有:
Y=lglg[1-F(N)]-1;X=lg(N-N0);
系数:A=lglge-blg(Na-N0);B=b设N0=2×105,列表计算后得到回归系数为:B=Lxy/Lxx=1.7196,A=Y-BX=-1.7985Weibull分布YX-1.29120.3010-0.9620-0.7543-0.5931-0.4532-0.3213-0.1849-0.02030.47710.60210.72430.77820.84510.93451.041443故威布尔分布参数:
b=B=1.7196,Na=lg-1[(lglge-A)/b]+N0=8.84×105。注意:对于本例,威布尔分布给出比正态分布更好的拟合精度,即更大的r值。相关系数为:r=-0.9988
r
ra=0.834(=0.01,n=8)估计寿命N=3×105时的破坏概率pf,则有:pf==3.6%FN()exp[(.)].=---13288417196443.5S-N曲线和P-S-N曲线的拟合实验得到:Ly12铝合金板材,在Smax为199、166、141.2、120.2Mpa四种应力水平下的疲劳试验结果x=lgN,循环应力比R=0.1S-N曲线和P-S-N曲线拟合计算实例试用最小二乘法拟合S-N曲线和P-S-N曲线。45序号 Smax(MPa)破坏率存活率
i199166141.2120.2Pi=i/n+1Ri=1-pi14.9145.0935.3255.7210.09090.909124.9145.1275.3605.8510.18180.818234.9295.1305.4355.8590.27270.727344.9645.1405.4415.9380.36360.636454.9645.1465.4706.0120.45450.545564.9825.1675.4716.0150.54550.454574.9825.1885.5016.0820.63640.363684.9965.2045.5496.1360.72730.272795.0295.2205.5826.1380.81820.1818105.0635.2485.6126.1650.90910.0909464.55.05.56.06.5X=lgN99.999907050301010.1Ps100对数疲劳寿命分布表中数据在正态概率纸上描点结果如图。四种应力水平下的xps数据,均呈线性,即x=lgN,服从正态分布。Smax=199Smax=166Smax=141.2Smax=120.247各应力水平下的xup拟合结果Si
x=A+Bup=lgN i(MPa)AB r lgSPs=50%Ps=99.9%up=0up=-3.091199.04.97370.05660.9752.29894.97374.7988 2166.05.16630.05710.9882.22015.16634.9899 3141.25.47460.10840.9892.14985.47465.1396 4120.25.99170.17220.9732.07995.9917
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年份第二季度数据资产质押借款保证合同安全审计附件
- 2019-2025年期货从业资格之期货基础知识模考预测题库(夺冠系列)
- 2025租房合同模板CC
- 2025家居定制家具购销合同范本模板
- 2025冰箱供货合同范本
- 2025年中外合作经营合同示范文本
- 2025房屋买卖居间合同范本
- 2025建筑外墙涂料施工及景观绿化不锈钢围栏工程合同
- 养牛入股合同样本
- 机构职能体系 司法责任制
- 起诉保险公司的诉讼书范本
- 人教部编版六年级下册语文【选择题】专项复习训练真题100题(附答案解析)
- H3C新员工文化培训报到指引(201607期)
- 《功和机械能》 单元作业设计
- 《辅酶q10》教学讲解课件
- 第十章痰液检查课件
- 《融媒体实务》教学课件(全)
- 牛津译林版六年级下册英语期中检测试卷+答案
- 重庆农艺师考试(种植业卷)
- 散文阅读理解文中重要句子的含意公开课一等奖市优质课赛课获奖课件
- 2023学年完整公开课版《认识洗衣机》
评论
0/150
提交评论