发明名称 油气弹簧不等厚环形阀片最大周向应力的计算方法
摘要 本发明涉及油气弹簧不等厚环形阀片最大周向应力的计算方法,属于油气悬架技术领域。先前对于油气弹簧不等厚环形阀片最大周向应力一直没有可靠的计算方法,大都利用有限元软件对给定压力下最大周向应力进行仿真,不能满足油气弹簧现代化CAD设计的要求。本发明提供的油气弹簧不等厚环形阀片最大周向应力的计算方法,可根据油气弹簧不等厚环形阀片的结构参数和材料性能参数,对最大周向应力进行精确、快速的计算。通过与ANSYS仿真结果比较可知,该计算方法是准确、可靠的,为实现油气弹簧现代化CAD设计和不等厚环形阀片应力强度计算,提供了可靠的油气弹簧不等厚环形阀片最大周向应力计算方法。
申请公布号 CN105279319A 申请公布日期 2016.01.27
申请号 CN201510646317.7 申请日期 2015.10.08
申请人 山东理工大学 发明人 周长城;赵雷雷;潘礼军;王凤娟;邵明磊;于曰伟
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 油气弹簧不等厚环形阀片最大周向应力的计算方法,其具体计算步骤如下:(1)确定不等厚环形阀片最大周向应力系数构成的常数项X<sub>1</sub>和X<sub>2</sub>:根据不等厚环形阀片的弹性模量E,泊松比μ,变厚度半径r<sub>t</sub>,有效内圆半径r<sub>a</sub>,外圆半径r<sub>b</sub>,建立不等厚环形阀片最大周向应力系数构成的常数项的特征方程,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>X</mi><mn>1</mn></msub><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup><mo>-</mo><mfrac><mrow><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mrow><mn>4</mn><mi>E</mi></mrow></mfrac><mrow><mo>(</mo><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><msub><mi>lnr</mi><mi>a</mi></msub><mo>+</mo><mn>2</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>+</mo><msub><mi>X</mi><mn>2</mn></msub><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000816726750000011.GIF" wi="884" he="119" /></maths><maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mfrac><mrow><mn>85</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup></mrow><mn>96</mn></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>Er</mi><mi>t</mi><mn>3</mn></msubsup><mrow><mo>(</mo><mn>11</mn><mo>+</mo><mn>3</mn><msqrt><mn>17</mn></msqrt><mo>)</mo></mrow></mrow><mrow><mn>72</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>r</mi><mi>b</mi><mrow><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn><mo>-</mo><mn>5</mn><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mi>Y</mi><mn>2</mn></msub><mo>+</mo><mfrac><mrow><msubsup><mi>Er</mi><mi>t</mi><mn>3</mn></msubsup><mrow><mo>(</mo><mn>11</mn><mo>-</mo><mn>3</mn><msqrt><mn>17</mn></msqrt><mo>)</mo></mrow></mrow><mrow><mn>72</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>r</mi><mi>b</mi><mrow><mo>-</mo><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn><mo>-</mo><mn>5</mn><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mi>Y</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000816726750000012.GIF" wi="1166" he="150" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>X</mi><mn>1</mn></msub><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup><mo>+</mo><msub><mi>X</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mrow><mn>4</mn><msub><mi>Er</mi><mi>t</mi></msub></mrow></mfrac><mrow><mo>(</mo><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup><mo>-</mo><mn>4</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><msub><mi>lnr</mi><mi>t</mi></msub><mo>+</mo><mn>2</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>-</mo><msub><mi>Y</mi><mn>1</mn></msub><msubsup><mi>r</mi><mi>t</mi><mrow><mn>5</mn><mo>/</mo><mn>2</mn><mo>-</mo><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn></mrow></msubsup><mo>-</mo><msub><mi>Y</mi><mn>2</mn></msub><msubsup><mi>r</mi><mi>t</mi><mrow><mn>5</mn><mo>/</mo><mn>2</mn><mo>+</mo><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn></mrow></msubsup><mo>+</mo><mfrac><mrow><mn>3</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup></mrow><mrow><mn>8</mn><mi>E</mi></mrow></mfrac><mrow><mo>(</mo><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup><mo>-</mo><mn>8</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><mo>)</mo></mrow><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000816726750000013.GIF" wi="1686" he="132" /></maths><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mfrac><mrow><mn>19</mn><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup></mrow><mn>96</mn></mfrac><mo>-</mo><mfrac><mrow><mn>13</mn><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup></mrow><mn>12</mn></mfrac><mo>+</mo><mfrac><mi>E</mi><mrow><mn>9</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msub><mi>X</mi><mn>1</mn></msub><mo>-</mo><mfrac><mi>E</mi><mrow><mn>18</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow><msup><mi>r</mi><mn>2</mn></msup></mrow></mfrac><msub><mi>X</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><mn>5</mn><msubsup><mi>r</mi><mi>t</mi><mn>2</mn></msubsup></mrow><mn>24</mn></mfrac><mo>+</mo><mfrac><mrow><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><msub><mi>lnr</mi><mi>t</mi></msub></mrow><mn>3</mn></mfrac><mo>+</mo><mfrac><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><mn>12</mn></mfrac></mrow></mtd></mtr><mtr><mtd><mrow><mo>-</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mn>11</mn><mo>+</mo><mn>3</mn><msqrt><mn>17</mn></msqrt><mo>)</mo></mrow></mrow><mrow><mn>72</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>r</mi><mi>t</mi><mrow><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn><mo>+</mo><mn>1</mn><mo>/</mo><mn>2</mn></mrow></msubsup><msub><mi>Y</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><mi>E</mi><mrow><mo>(</mo><mn>11</mn><mo>-</mo><mn>3</mn><msqrt><mn>17</mn></msqrt><mo>)</mo></mrow></mrow><mrow><mn>72</mn><mrow><mo>(</mo><mn>1</mn><mo>-</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow></mfrac><msubsup><mi>r</mi><mi>t</mi><mrow><mo>-</mo><msqrt><mn>17</mn></msqrt><mo>/</mo><mn>2</mn><mo>+</mo><mn>1</mn></mrow></msubsup><msub><mi>Y</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></mrow></mtd></mtr></mtable><mo>;</mo></mrow>]]></math><img file="FDA0000816726750000014.GIF" wi="1318" he="295" /></maths>利用Matlab程序,求解上述关于X<sub>1</sub>、X<sub>2</sub>、Y<sub>1</sub>和Y<sub>2</sub>的四个方程组成的方程组,求得不等厚环形阀片最大周向应力系数构成的常数项X<sub>1</sub>和X<sub>2</sub>;(2)确定不等厚环形阀片的最大周向应力系数G<sub>σθmax</sub>:根据油气弹簧不等厚环形阀片的弹性模量E,有效内圆半径r<sub>a</sub>,外圆半径r<sub>b</sub>,以及步骤(1)求得的不等厚环形阀片最大周向应力系数构成的常数项X<sub>1</sub>和X<sub>2</sub>,确定不等厚环形阀片的最大周向应力系数G<sub>σθmax</sub>,即:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>G</mi><mrow><mi>&sigma;</mi><mi>&theta;</mi><mi>max</mi></mrow></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><msub><mi>lnr</mi><mi>a</mi></msub></mrow><mn>2</mn></mfrac><mo>-</mo><mfrac><mrow><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><msub><mi>r</mi><mi>a</mi></msub><msub><mi>lnr</mi><mi>a</mi></msub></mrow><mn>2</mn></mfrac><mo>+</mo><mfrac><mrow><msub><mi>r</mi><mi>a</mi></msub><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup></mrow><mn>4</mn></mfrac><mo>-</mo><mfrac><mrow><mn>3</mn><msub><mi>EX</mi><mn>1</mn></msub><msub><mi>r</mi><mi>a</mi></msub></mrow><mn>16</mn></mfrac><mo>+</mo><mfrac><mrow><mn>3</mn><msub><mi>EX</mi><mn>2</mn></msub></mrow><mrow><mn>16</mn><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup></mrow></mfrac><mo>-</mo><mfrac><mrow><mn>3</mn><msub><mi>EX</mi><mn>2</mn></msub></mrow><mrow><mn>16</mn><msub><mi>r</mi><mi>a</mi></msub></mrow></mfrac><mo>+</mo><mfrac><mrow><mn>3</mn><msubsup><mi>r</mi><mi>a</mi><mn>2</mn></msubsup></mrow><mn>8</mn></mfrac><mo>+</mo><mfrac><msubsup><mi>r</mi><mi>a</mi><mn>3</mn></msubsup><mn>8</mn></mfrac><mo>-</mo><mfrac><msubsup><mi>r</mi><mi>b</mi><mn>2</mn></msubsup><mn>4</mn></mfrac><mo>-</mo><mfrac><mrow><mn>3</mn><msub><mi>EX</mi><mn>1</mn></msub></mrow><mn>16</mn></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000816726750000015.GIF" wi="1525" he="133" /></maths>(3)计算不等厚环形阀片的最大周向应力σ<sub>θmax</sub>:根据不等厚环形阀片的等厚度部分的厚度h<sub>0</sub>,所受均布压力p,步骤(2)求得的G<sub>σθmax</sub>,计算不等厚环形阀片的最大周向应力σ<sub>θmax</sub>,即:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mrow><mi>&theta;</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>=</mo><msub><mi>G</mi><mrow><mi>&sigma;</mi><mi>&theta;</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mfrac><mi>p</mi><msubsup><mi>h</mi><mn>0</mn><mn>2</mn></msubsup></mfrac><mo>.</mo></mrow>]]></math><img file="FDA0000816726750000016.GIF" wi="342" he="110" /></maths>
地址 255086 山东省淄博市高新技术产业开发区高创园A座313室