发明名称 一种基于最大熵模型的材料疲劳寿命的预测方法
摘要 本发明提出一种基于最大熵模型的材料疲劳寿命的预测方法,具体环节为:1.获取材料在某一横幅应力为S<sub>m</sub>下的疲劳试验数据;2.计算数据的前四阶矩,分别是均值、标准差、偏度和峰度;3.结合Lagrange乘子法写出最大熵分布的表达式;4.利用无约束最小化方法求解最大熵分布的参数;5.获取给定可靠度为p,应力水平为的S<sub>m</sub>疲劳寿命N;6.获取不同应力水平S<sub>m</sub>下材料的疲劳试验数据,重复前述五个步骤,获得多个可靠度为p,应力水平不同的疲劳寿命;7.三参数幂函数拟合P‑S‑N曲线,简化材料在实际使用中的载荷谱,结合Miner累计损伤理论,对材料的疲劳寿命做出预测。此方法相对于传统基于经验的对数正态或威布尔模型,具有更高的预测精度,降低了可靠性在工程应用中的主观性。
申请公布号 CN106383959A 申请公布日期 2017.02.08
申请号 CN201610859074.X 申请日期 2016.09.23
申请人 南京航空航天大学 发明人 李洪双;闻德兵;张航;赵翔;袁骄阳;任超;顾汝佳
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 一种基于最大熵模型的材料疲劳寿命的预测方法,其特征在于以下几个步骤:步骤一:获取材料在某一横幅应力为S<sub>m</sub>下的疲劳试验数据x<sub>1</sub>,x<sub>2</sub>,x<sub>3</sub>,...x<sub>n</sub>,包括试件断裂时的寿命和指定裂纹长度下寿命数据。步骤二:根据步骤一获得的寿命数据,计算它们的前四阶矩,分别是均值、标准差、偏度和峰度,记作<img file="FSA0000134973930000011.GIF" wi="127" he="46" />ske,kur。<maths num="0001"><math><![CDATA[<mrow><mover><mi>&mu;</mi><mo>^</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>x</mi><mi>i</mi></msub></mrow>]]></math><img file="FSA0000134973930000012.GIF" wi="262" he="146" /></maths><maths num="0002"><math><![CDATA[<mrow><msup><mover><mi>s</mi><mo>^</mo></mover><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>&mu;</mi><mo>^</mo></mover><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FSA0000134973930000013.GIF" wi="481" he="144" /></maths><maths num="0003"><math><![CDATA[<mrow><mi>s</mi><mi>k</mi><mi>e</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>&mu;</mi><mo>^</mo></mover><mo>)</mo></mrow><mn>3</mn></msup><mo>/</mo><msup><mover><mi>s</mi><mo>^</mo></mover><mn>3</mn></msup></mrow>]]></math><img file="FSA0000134973930000014.GIF" wi="552" he="134" /></maths><maths num="0004"><math><![CDATA[<mrow><mi>k</mi><mi>u</mi><mi>r</mi><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>-</mo><mover><mi>&mu;</mi><mo>^</mo></mover><mo>)</mo></mrow><mn>4</mn></msup><mo>/</mo><msup><mover><mi>s</mi><mo>^</mo></mover><mn>4</mn></msup></mrow>]]></math><img file="FSA0000134973930000015.GIF" wi="637" he="141" /></maths>步骤三:由步骤二的样本前四阶统计矩,计算Lagrange乘子,并确定最大熵分布的具体表达式。最大熵分布具有如下形式的表达式<maths num="0005"><math><![CDATA[<mrow><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msup><mi>x</mi><mi>i</mi></msup><mo>)</mo></mrow></mrow>]]></math><img file="FSA0000134973930000016.GIF" wi="778" he="164" /></maths>式中:λ<sub>i</sub>,i=1,…,m为Lagrange乘子。标准化因子λ<sub>0</sub>为<maths num="0006"><math><![CDATA[<mrow><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>=</mo><mi>l</mi><mi>n</mi><mo>&lsqb;</mo><mo>&Integral;</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msup><mi>x</mi><mi>i</mi></msup><mo>)</mo></mrow><mi>d</mi><mi>x</mi><mo>&rsqb;</mo></mrow>]]></math><img file="FSA0000134973930000017.GIF" wi="669" he="177" /></maths>一般的,Lagrange乘子可以利用标准牛顿方法求解N个非线性等式来获得。对比求解N个非线性等式的方法,提出采用无约束最小化方法,该方法较传统方法效率更高、数值更稳定。步骤四:获取给定可靠度为p,应力水平为S<sub>m</sub>的疲劳寿命N。在步骤三获得最大熵分布的具体表达式后,由<maths num="0007"><math><![CDATA[<mrow><mi>p</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>N</mi></msubsup><mi>f</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mi>d</mi><mi>x</mi><mo>=</mo><msubsup><mo>&Integral;</mo><mn>0</mn><mi>N</mi></msubsup><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>&lambda;</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>&lambda;</mi><mi>i</mi></msub><msup><mi>x</mi><mi>i</mi></msup><mo>)</mo></mrow><mi>d</mi><mi>x</mi></mrow>]]></math><img file="FSA0000134973930000021.GIF" wi="946" he="135" /></maths>求得给定可靠度P下的疲劳寿命N。步骤五:由前述四个步骤,获取不同应力水平S′<sub>m</sub>下材料的疲劳试验数据,获得多个可靠度为p,应力水平不同的疲劳寿命。步骤六:由步骤五,利用三参数幂函数拟合p‑S‑N曲线,简化材料在实际使用中的载荷谱,结合Miner累计损伤理论,对材料的疲劳寿命做出预测。
地址 210018 江苏省南京市秦淮区御道街29号