发明名称 一种提高基于不完全数据飞机结构可靠性分析精度的方法
摘要 本发明提供一种提高基于不完全数据飞机结构可靠性分析精度的方法,对传统方法确定安全寿命的精度与基于不完全数据的飞机结构可靠性分析方法确定安全寿命的精度进行比较,确定精度修正系数k以及最少试件数n,使利用基于不完全数据飞机结构可靠性分析方法确定安全寿命的精度与传统方法确定安全寿命的精度相同。本发明解决了基于不完全数据飞机结构可靠性分析方法的精度问题,对飞机结构进行可靠性分析时,确保基于不完全数据飞机结构可靠性分析方法的可信度与传统安全寿命确定方法的可信度基本相同。
申请公布号 CN105320826A 申请公布日期 2016.02.10
申请号 CN201410306803.X 申请日期 2014.07.01
申请人 中国人民解放军空军工程大学 发明人 何宇廷;高潮;崔荣洪;杜金强;安涛
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 中国人民解放军空军专利服务中心 11035 代理人 孙艳
主权项 一种提高基于不完全数据飞机结构可靠性分析精度的方法,包括以下步骤:步骤1:传统飞机结构安全寿命分析方法的精度分析确定飞机结构寿命母体的安全寿命Np;根据失效数据确定安全寿命估计值Np′;确定Np′个数n<sub>1</sub>,Np′小于Np的个数m<sub>1</sub>,计算存活率P<sub>1</sub>:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>P</mi><mn>1</mn></msub><mo>=</mo><mfrac><msub><mi>m</mi><mn>1</mn></msub><msub><mi>n</mi><mn>1</mn></msub></mfrac></mrow>]]></math><img file="FSA0000105748870000011.GIF" wi="169" he="130" /></maths>可用存活率P<sub>1</sub>表示采用传统方法确定飞机结构安全寿命的精度;步骤2:基于不完全数据飞机结构可靠性分析方法的精度分析根据基于不完全数据的飞机结构可靠性分析方法确定安全寿命Np<sub>d</sub>″;确定Np<sub>d</sub>″个数n<sub>1</sub>,Np<sub>d</sub>″小于Np的个数m<sub>2</sub>,计算存活率P<sub>2</sub>:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>P</mi><mn>2</mn></msub><mo>=</mo><mfrac><msub><mi>m</mi><mn>2</mn></msub><msub><mi>n</mi><mn>2</mn></msub></mfrac></mrow>]]></math><img file="FSA0000105748870000012.GIF" wi="181" he="130" /></maths>存活率P<sub>2</sub>表示基于不完全数据飞机结构可靠性分析方法确定飞机结构安全寿命的精度;步骤3:确定精度修正系数κ步骤(1):生成随机数随机产生j组服从对数正态分布/双参数威布尔分布的随机数,用来表示飞机结构的寿命特征,每组随机数的样本容量为l;j组随机数则代表有j组飞机结构失效试验数据,而每组随机数中的l个随机数代表每组有l个飞机结构失效试验数据;步骤(2):确定每组随机数的Np、Np<sub>d</sub>″和Np′①根据l个失效试验数据确定理论上的飞机结构寿命母体安全寿命Np;②随机抽取r个试验数据作为部分失效试验数据,看作为飞机结构疲劳/耐久性试验结果,利用传统方法估计确定飞机结构安全寿命值Np′;③将②中的样本失效数据作为不完全数据中的失效数据,无失效数据为传统安全寿命Np′,根据基于不完全数据的飞机结构可靠性分析方法估计确定飞机结构安全寿命值Np<sub>d</sub>″;步骤(3):确定每组随机数的P<sub>i,1</sub>和P<sub>i,2</sub>根据第i组数据确定Np′的个数n<sub>i,1</sub>:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><msubsup><mi>C</mi><mi>l</mi><mi>r</mi></msubsup></mrow>]]></math><img file="FSA0000105748870000021.GIF" wi="226" he="83" /></maths>统计Np′小于等于Np的个数m<sub>i,1</sub>,确定Np′的存活率P<sub>i,1</sub>:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>m</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mn>1</mn></mrow></msub></mfrac></mrow>]]></math><img file="FSA0000105748870000022.GIF" wi="199" he="135" /></maths>根据第i组数据确定Np<sub>d</sub>″的个数n<sub>i,2</sub>:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>=</mo><msubsup><mi>C</mi><mi>l</mi><mi>r</mi></msubsup></mrow>]]></math><img file="FSA0000105748870000023.GIF" wi="233" he="82" /></maths>统计Np<sub>d</sub>″小于等于Np的个数m<sub>i,2</sub>,确定Np<sub>d</sub>″的存活率P<sub>i,2</sub>:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>P</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub><mo>=</mo><mfrac><msub><mi>m</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub><msub><mi>n</mi><mrow><mi>i</mi><mo>,</mo><mn>2</mn></mrow></msub></mfrac></mrow>]]></math><img file="FSA0000105748870000024.GIF" wi="210" he="135" /></maths>根据上述方法确定每组随机数的P<sub>i,1</sub>和P<sub>i,2</sub>,其中1≤i≤j;步骤(4):确定精度修正系数k对第i组数据<img file="FSA0000105748870000025.GIF" wi="81" he="90" />个Np<sub>d</sub>″都同除以一个精度修正系数k<sub>i</sub>进行修正,统计修正后的Np<sub>d</sub>″小于等于Np的个数,确定P<sub>i,2</sub>,使修正后的P<sub>i,2</sub>与P<sub>i,1</sub>相同;确定每组数据Np<sub>d</sub>″的精度修正系数k<sub>i</sub>;对j组数据Np<sub>d</sub>″的精度修正系数k<sub>i</sub>取平均值可作为Np<sub>d</sub>″的精度修正系数k:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>l</mi></munderover><msub><mi>k</mi><mi>i</mi></msub></mrow><mi>j</mi></mfrac></mrow>]]></math><img file="FSA0000105748870000031.GIF" wi="197" he="223" /></maths>步骤4:基于不完全数据飞机结构可靠性分析方法的飞机结构安全寿命的修正利用步骤3中确定的精度修正系数k对步骤2中的飞机结构安全寿命值Np<sub>d</sub>″进行修正;对机结构安全寿命Np<sub>d</sub>″除以精度修正系数k就可以得到修正后的安全寿命值Np<sub>d</sub>″;步骤5:确定最少试件个数根据极大似然估计得到的<img file="FSA0000105748870000032.GIF" wi="51" he="97" />作为子样平均值<img file="FSA0000105748870000033.GIF" wi="72" he="77" />对数疲劳寿命标准差σ<sub>0</sub>作为子样标准差s,计算变异系数<img file="FSA0000105748870000034.GIF" wi="129" he="77" />确定最少试件数n;步骤6:验证精度修正系数k及最少试件数n根据试验数据与随机数确定飞机结构寿命母体的安全寿命Np、安全寿命Np′以及安全寿命Np<sub>d</sub>″;利用精度修正系数k对飞机结构安全寿命Np<sub>d</sub>″进行修正,Np′的存活率P<sub>1</sub>和修正后Np<sub>d</sub>″的存活率P<sub>2</sub>;比较P<sub>1</sub>与P<sub>2</sub>的大小,当P<sub>1</sub>与P<sub>2</sub>基本上相同时,说明上述所确定的精度修正系数k及最少试件数n的可行性。
地址 710051 陕西省西安市长乐东路甲字一号