发明名称 一种基于寿命理论计算的作动器产品加速因子确定方法
摘要 本发明一种基于寿命理论计算的作动器产品加速因子确定方法,其具体步骤如下:1.主机理分析,确定作动器产品薄弱环节及其对应的主机理;2.寿命计算模型选择,根据主机理类型选择对应的寿命计算模型;3.加速因子综合确定,根据确定的应力类型与范围计算不同应力水平下产品各薄弱环节对应主机理的加速因子,根据加速因子标准差最小原则确定加速应力水平,再根据加速因子取小原则确定加速因子。本发明可以直接应用于产品,通过计算主机理分析所得到的产品各薄弱环节的加速因子确定产品的加速因子,为产品层的加速寿命试验方法研究提供理论支撑,具有思路明确、可操作性强等优点。
申请公布号 CN104156600B 申请公布日期 2017.04.05
申请号 CN201410396411.7 申请日期 2014.08.13
申请人 北京航空航天大学 发明人 陈云霞;井海龙;康锐;廖寻
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 北京慧泉知识产权代理有限公司 11232 代理人 王顺荣;唐爱华
主权项 一种基于寿命理论计算的作动器产品加速因子确定方法,其特征在于:该方法具体步骤如下:步骤一:主机理分析;确定作动器产品的薄弱环节及其对应的耗损型失效机理,主机理分析根据给定的载荷谱或任务剖面,结合作动器产品的组成、结构、原理,进行故障模式、机理和影响分析即FMMEA,在FMMEA的基础上根据各机理对应的故障模式的严酷度及发生频度,综合确定作动器产品的薄弱环节及其对应的主机理;步骤二:寿命计算模型选择;根据步骤一中主机理分析结果和对应的寿命计算模型,代入公式计算作动器产品各薄弱环节的寿命;主要包括:a.疲劳类机理:作动器产品疲劳类机理运用名义应力法结合线性累积损伤原理式(1)~式(7)计算薄弱环节的寿命;<maths num="0001"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>m</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mrow><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>+</mo><msub><mi>&sigma;</mi><mi>min</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000011.GIF" wi="1078" he="111" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>a</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>&sigma;</mi><mi>max</mi></msub><mo>-</mo><msub><mi>&sigma;</mi><mi>min</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000012.GIF" wi="1078" he="116" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>e</mi></msub><mo>=</mo><mfrac><msub><mi>&sigma;</mi><mi>a</mi></msub><mrow><mn>1</mn><mo>-</mo><mfrac><msub><mi>&sigma;</mi><mi>m</mi></msub><msub><mi>&sigma;</mi><mi>b</mi></msub></mfrac></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000013.GIF" wi="1022" he="191" /></maths><maths num="0004"><math><![CDATA[<mrow><msup><mrow><mo>(</mo><mfrac><msub><mi>&sigma;</mi><mi>e</mi></msub><msub><mi>&sigma;</mi><mrow><mo>-</mo><mn>1</mn><mi>A</mi></mrow></msub></mfrac><mo>)</mo></mrow><mfrac><mn>1</mn><mi>b</mi></mfrac></msup><mo>=</mo><mfrac><msub><mi>N</mi><mi>i</mi></msub><msub><mi>N</mi><mn>0</mn></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000014.GIF" wi="1046" he="150" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mrow><mo>-</mo><mn>1</mn><mi>A</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><mfrac><mrow><msub><mi>&epsiv;&alpha;&beta;C</mi><mi>L</mi></msub></mrow><msub><mi>K</mi><mi>f</mi></msub></mfrac><mo>)</mo></mrow><msub><mi>&sigma;</mi><mrow><mo>-</mo><mn>1</mn></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000015.GIF" wi="1102" he="151" /></maths><maths num="0006"><math><![CDATA[<mrow><mfrac><msub><mi>n</mi><mi>i</mi></msub><msub><mi>N</mi><mi>i</mi></msub></mfrac><mo>=</mo><msub><mi>D</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000016.GIF" wi="990" he="134" /></maths><maths num="0007"><math><![CDATA[<mrow><mi>L</mi><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>&Sigma;D</mi><mi>i</mi></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000017.GIF" wi="1012" he="151" /></maths>其中:σ<sub>max</sub>为最大应力;σ<sub>min</sub>为最小应力;σ<sub>m</sub>为应力均值;σ<sub>a</sub>为应力幅值;σ<sub>b</sub>材料强度极限;σ<sub>e</sub>平均应力修正后的等效疲劳强度;σ<sub>‑1A</sub>零件疲劳极限;ε尺寸系数;α可靠性系数;β表面质量系数;K<sub>f</sub>疲劳缺口系数;σ<sub>‑1</sub>材料疲劳极限;N<sub>i</sub>各级载荷下的循环次数;N<sub>0</sub>应力幅为疲劳极限时的循环次数;b疲劳强度指数;n<sub>i</sub>各级载荷下的工作循环数;D<sub>i</sub>各级载荷下的损伤;L疲劳寿命安全系数;b.磨损类机理:作动器产品磨损类机理为粘着磨损,运用Archard阿恰德粘着磨损原理式(8)~式(9)计算薄弱环节的寿命;<maths num="0008"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><mi>K</mi><mfrac><msub><mi>F</mi><mi>N</mi></msub><mi>H</mi></mfrac><msub><mi>L</mi><mi>m</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>8</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000021.GIF" wi="926" he="115" /></maths>logK=5logμ‑2.27                  (9)其中:V磨损体积;K磨损系数;μ摩擦系数;F<sub>N</sub>法向载荷;H摩擦副材料硬度;L<sub>m</sub>磨损行程;公式(8)中若摩擦副为直接接触则法向载荷F<sub>N</sub>直接带入对应的正压力值,若摩擦副为间隙配合,其法向载荷F<sub>N</sub>用微凸体载荷Wa代替:<maths num="0009"><math><![CDATA[<mrow><msub><mi>W</mi><mi>a</mi></msub><mo>=</mo><mfrac><mn>1</mn><mn>30</mn></mfrac><mi>A</mi><mi>H</mi><mi>&psi;</mi><mo>&lsqb;</mo><msub><mi>F</mi><mfrac><mn>3</mn><mn>2</mn></mfrac></msub><mrow><mo>(</mo><mfrac><mi>h</mi><mi>&sigma;</mi></mfrac><mo>)</mo></mrow><mo>-</mo><msub><mi>F</mi><mfrac><mn>3</mn><mn>2</mn></mfrac></msub><mrow><mo>(</mo><mfrac><mi>h</mi><mi>&sigma;</mi></mfrac><mo>+</mo><mfrac><mn>1</mn><msup><mi>&psi;</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>&rsqb;</mo><mo>+</mo><mfrac><mn>1</mn><mn>10</mn></mfrac><msub><mi>&pi;AHF</mi><mn>1</mn></msub><mrow><mo>(</mo><mfrac><mi>h</mi><mi>&sigma;</mi></mfrac><mo>+</mo><mfrac><mn>1</mn><msup><mi>&psi;</mi><mn>2</mn></msup></mfrac><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000022.GIF" wi="1398" he="133" /></maths><maths num="0010"><math><![CDATA[<mrow><msub><mi>F</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>u</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Integral;</mo><mi>u</mi><mi>&infin;</mi></munderover><msup><mrow><mo>(</mo><mi>t</mi><mo>-</mo><mi>u</mi><mo>)</mo></mrow><mi>n</mi></msup><msup><mi>e</mi><mrow><mo>-</mo><mfrac><msup><mi>t</mi><mn>2</mn></msup><mn>2</mn></mfrac></mrow></msup><mi>d</mi><mi>t</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000023.GIF" wi="1054" he="147" /></maths>其中:W<sub>a</sub>微凸体载荷;A名义接触面积;ψ塑性指数;h/σ膜厚比;t橡胶密封件热氧老化寿命;c.老化类机理:作动器产品老化类机理为橡胶固定密封件的老化,通过老化特征量与老化时间的关系并联立老化速度常数与温度间的Arrhenius关系建立其老化寿命计算模型;1‑y=A<sub>1</sub>exp{‑kt<sup>α</sup>}                      (12)<maths num="0011"><math><![CDATA[<mrow><mi>k</mi><mo>=</mo><msub><mi>A</mi><mn>2</mn></msub><mi>exp</mi><mo>{</mo><mo>-</mo><mfrac><mi>E</mi><mrow><mi>R</mi><mi>T</mi></mrow></mfrac><mo>}</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>13</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000024.GIF" wi="1030" he="119" /></maths>其中:t橡胶密封件热氧老化寿命;y压缩永久变形;A<sub>1</sub>零时刻压缩永久变形常数;k老化速度常数;α老化反应时间指数;A<sub>2</sub>频率因子;E橡胶材料的活化能;R气体常数;T老化反应温度;步骤三:加速因子综合确定:基于寿命计算模型结合产品工作极限与试验设备加载能力确定加速应力类型与范围,并带入公式计算不同应力水平下产品各薄弱环节的加速因子,最后根据标准差最小原则与加速因子取小原则综合确定产品的加速因子;主要包括:a.应力类型与范围确定;基于产品主机理分析结果,结合步骤二中所选择的寿命计算模型以及产品工作极限与试验设备加载能力确定加速应力类型与范围;b.不同应力水平下加速因子计算;在所确定的应力类型与范围内带入公式分别计算不同应力水平下产品各薄弱环节所对应主机理的加速因子A<sub>f</sub>,加速因子按下式计算:L<sub>D0</sub>常规载荷谱累积损伤达到D时的寿命;L<sub>Dr</sub>加速载荷谱累积损伤达到D时的寿命;<maths num="0012"><math><![CDATA[<mrow><msub><mi>A</mi><mi>f</mi></msub><mo>=</mo><mfrac><msub><mi>L</mi><mrow><mi>D</mi><mn>0</mn></mrow></msub><msub><mi>L</mi><mrow><mi>D</mi><mi>r</mi></mrow></msub></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000031.GIF" wi="1010" he="143" /></maths>c.加速应力水平确定;分别计算所选加速应力水平下产品各薄弱环节所对应主机理的加速因子标准差σ<sub>AF</sub>:<maths num="0013"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mrow><mi>A</mi><mi>F</mi></mrow></msub><mo>=</mo><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>&lsqb;</mo><msub><mi>A</mi><mrow><mi>f</mi><mi>x</mi><mi>y</mi></mrow></msub><mo>-</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>x</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msub><mi>A</mi><mrow><mi>f</mi><mi>x</mi><mi>y</mi></mrow></msub><mo>)</mo></mrow><mo>&rsqb;</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001112493320000032.GIF" wi="1253" he="147" /></maths>其中n为步骤一中确定的产品主机理数量;σ<sub>AF</sub>每级加速应力水平下的加速因子标准差;A<sub>fxy</sub>产品第x个主机理在第y个应力水平下的加速因子;取σ<sub>AF</sub>最小的加速应力水平作为加速试验的加速应力水平;d.加速因子综合确定:在确定的加速应力水平下根据取小原则取该应力水平下产品各薄弱环节所对应主机理中加速因子的最小值为产品的加速因子。
地址 100191 北京市海淀区学院路37号