CN103678775B - 一种输送管动强度分析方法 - Google Patents
一种输送管动强度分析方法 Download PDFInfo
- Publication number
- CN103678775B CN103678775B CN201310577106.3A CN201310577106A CN103678775B CN 103678775 B CN103678775 B CN 103678775B CN 201310577106 A CN201310577106 A CN 201310577106A CN 103678775 B CN103678775 B CN 103678775B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- delivery pipe
- mfrac
- msup
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Pipeline Systems (AREA)
Abstract
一种输送管动强度分析方法，本发明由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分组成，首先通过有限元软件Abaqus的static general模块、Frequency模块、Random reaponse模块依次实现输送管的静强度分析、模态分析和随机响应分析，然后采用Mises应力谱密度估算方法计算出Mises应力谱密度，基于随机振动的功率谱密度法，得到输送管的定量疲劳损伤后进行评估。本发明能够评价输送管动强度是否满足要求，在管路产品设计时提前揭示风险，提高了输送管动强度评价效率。
Description
技术领域
本发明涉及一种适用于输送管的动强度分析方法，属于运载火箭力学环境技术领域。
背景技术
疲劳失效是弹（箭）体结构件最常见的破坏模式之一，随着在研运载型号推力的提升，输送管特别是发动机舱内的输送管振动量级相应增加，载荷环境比较恶劣，动强度成为输送管是否满足使用要求的关键指标。
传统的疲劳寿命估算方法是通过试验获得危险点的应力—时间历程,采用不同的计数方法,将不规则应力时间历程转化为不同应力幅下的循环次数,再利用疲劳累积损伤理论估算疲劳损伤及寿命。但对于输送管，由于缺乏试验条件，无法获得危险点的应力时间历程，因此只能采用基于频域的统计方法来估算疲劳寿命。
目前基于频域的统计方法仅针对振动载荷下的动强度进行分析，未考虑输送管承受内压、温度及位移等静载荷对动强度的影响；仅针对实体结构的动强度进行分析，未考虑输送管内部携带推进剂重量的影响；采用应力分量功率谱密度作为参量进行疲劳寿命估算，则输送管处于多轴应力状态下时，忽略了其他应力分量的影响。
发明内容
本发明的技术解决问题是：克服现有技术的不足，提供一种输送管动强度分析方法，本方法通过采集频率域内的主要信息-功率谱密度，根据频率域内载荷过程的统计量分析完成输送管路的强度分析，本发明解决了管路疲劳寿命的评估标准问题，填补了管路设计领域在该专业方面的空白。
本发明的技术解决方案是：
一种输送管动强度分析方法包括步骤如下：
（1）首先将推进剂重量等效处理后，在有限元软件Abaqus中对输送管进行建模，所述的输送管道建模包括输送管的几何属性和材料属性；
（2）判断输送管模型几何和材料的正确性，若不正确返回步骤（1），若正确则对模型施加内压、位移和温度载荷，利用有限元软件Abaqus（采用static general模块）进行静强度分析；
（3）判断输送管静强度是否满足要求，若不满足要求，改进输送管道结构返回步骤（1），若满足要求，则利用有限元软件Abaqus（采用Frequency模块）进行模态分析；
（4）利用有限元软件Abaqus（采用Random response模块），结合步骤（3）中分析得到的模态以及对模型施加随机激励（例如：加速度功率谱密度）进行随机响应分析得到的应力分量均方根及功率谱密度；
（5）对步骤（4）中分析得到的应力分量均方根及功率谱密度，根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度；
（6）根据步骤（5）中计算得到的Mises应力功率谱密度，由Dirlik公式计算出输送管的疲劳损伤；
（7）判断输送管疲劳损伤是否满足要求，若不满足要求，改进输送管道结构后返回步骤（1），若满足要求，则结束。
所述步骤（1）中推进剂重量等效处理包括如下步骤：
输送管等效密度ρ通过公式得到：
其中，D为输送管路内径，δ为输送管路壁厚，ρ1为输送管材密度，ρ2为推进剂密度。
所述步骤（3）中判断输送管静强度是否满足要求包括如下步骤：
（a）输送管静强度安全因子η通过如下公式得到：
其中，σ为静强度分析得到的输送管Mises应力，σs为管路材料在设计温度下的屈服强度；
（b）若η≥1，则输送管静强度满足要求，反之若η<1不满足要求。
所述步骤（5）中Mises应力功率谱密度估算方法如下：
输送管Mises应力谱密度G(σM)通过如下公式得到：
其中，G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度。
所述步骤（6）中疲劳损伤计算方法如下：
输送管疲劳损伤DL通过如下公式得到：
其中，C、m为材料疲劳曲线参数，σb为管路材料在设计温度下的强度极限，T为输送管振动时间，p(S)为Mises应力概率密度，E[p]为单位之间内的峰值数。
所述单位之间内的峰值数E[p]通过如下公式得到：
其中，f1为随机激励的频率下限，f2为随机激励的频率上限。
所述Mises应力概率密度p(S)通过公式：
其中，
所述步骤（7）中判断输送管疲劳损伤是否满足要求方法如下：
若输送管疲劳损伤DL≥0.1，则输送管疲劳损伤不满足要求，反之若DL<0.1满足要求。
本发明与现有技术相比的有益效果是：
（1）本发明给出了由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分依次组成的动强度分析流程，采用有限元软件Abaqus的static general模块、Frequency模块、Random reaponse模块实现输送管的静强度分析、模态分析和随机响应分析；通过该流程，内压、温度及位移等静载荷对动强度的影响将通过静强度分析向模态分析、随机响应分析和疲劳损伤计算评估传递，避免了单纯振动载荷下动强度分析的片面性；
（2）本发明对推进剂重量进行等效，可以在动强度分析中引入推进剂重量的影响。等效处理的好处是，既灵活的考虑了推进剂重量，又避免了需要对推进剂进行额外的建模、计算工作，使计算效率得到了保证。
（3）本发明给出了Mises应力谱密度与各应力分量应力谱密度的换算关系，采用Mises应力谱密度作为疲劳寿命估算参量，较之单一应力分量谱密度，充分考虑了输送管多轴应力状态的影响，可以更全面、更精确的评估输送管的疲劳寿命。
附图说明
图1为本发明方法流程图。
具体实施方式
下面结合附图对本发明的具体实施方式进行进一步的详细描述。
如图1所示，本发明一种输送管动强度分析方法，由静强度分析、模态分析、随机响应分析和疲劳损伤计算评估六部分组成，具体步骤如下：
（1）首先将推进剂重量等效处理后，在有限元软件Abaqus中对输送管进行建模，所述的输送管道建模包括输送管的几何属性和材料属性；
推进剂重量等效处理包括如下步骤：
输送管等效密度ρ通过公式得到：
其中，D为输送管路内径，δ为输送管路壁厚，ρ1为输送管材密度，ρ2为推进剂密度。
（2）判断输送管模型几何属性和材料属性的正确性（几何属性和材料属性的要求根据用户实际的需求相一致），若不正确返回步骤（1），若正确则对模型施加内压、位移和温度载荷，利用有限元软件Abaqus（采用有限元软件Abaqus中的static general模块）进行静强度分析；
（3）判断输送管静强度是否满足要求，若不满足要求，改进输送管道结构返回步骤（1），若满足要求，则利用有限元软件Abaqus（采用有限元软件Abaqus中的Frequency模块，截断频率阶数建议取值为20，小于等于20阶，采用subspace求解器；大于20阶，采用Lanczos求解器）进行模态分析；
判断输送管静强度是否满足要求包括如下步骤：
（a）输送管静强度安全因子η通过如下公式得到：
其中，σ为静强度分析得到的输送管Mises应力，σs为管路材料在设计温度下的屈服强度；
（b）若η≥1，则输送管静强度满足要求，反之若η<1不满足要求。
（4）利用有限元软件Abaqus（采用有限元软件Abaqus中的Random response模块，阻尼系数建议取值为0.03），结合步骤（3）中分析得到的模态以及对模型施加随机激励（例如：加速度功率谱密度）进行输送管道的随机响应分析得到应力分量均方根及功率谱密度；
（5）对步骤（4）中分析得到的应力分量均方根及功率谱密度，根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度；
Mises应力功率谱密度估算方法如下：
输送管Mises应力谱密度G(σM)通过如下公式得到：
其中，G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度。
（6）根据步骤（5）中计算得到的Mises应力功率谱密度，由Dirlik公式计算出输送管的疲劳损伤；
疲劳损伤计算方法如下：
输送管疲劳损伤DL通过如下公式得到：
其中，C、m为材料疲劳曲线参数，σb为管路材料在设计温度下的强度极限，T为输送管振动时间，p(S)为Mises应力概率密度，E[p]为单位之间内的峰值数。
单位之间内的峰值数E[p]通过如下公式得到：
其中，f1为随机激励的频率下限，f2为随机激励的频率上限。
Mises应力概率密度p(S)通过公式：
其中，
（7）判断输送管疲劳损伤是否满足要求，若不满足要求，改进输送管道结构后返回步骤（1），若满足要求，则结束。
判断输送管疲劳损伤是否满足要求方法如下：
若输送管疲劳损伤DL≥0.1，则输送管疲劳损伤不满足要求，反之若DL<0.1满足要求。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。
Claims (1)
1.一种输送管动强度分析方法，其特征在于步骤如下：
(1)首先将推进剂重量等效处理后，在有限元软件Abaqus中对输送管进行建模，所述的输送管道建模包括输送管的几何属性和材料属性；
(2)判断输送管模型几何和材料的正确性，若不正确返回步骤(1)，若正确则对模型施加内压、位移和温度载荷，利用有限元软件Abaqus进行静强度分析；
(3)判断输送管静强度是否满足要求，若不满足要求，改进输送管道结构返回步骤(1)，若满足要求，则利用有限元软件Abaqus进行模态分析；
(4)利用有限元软件Abaqus，并结合步骤(3)中分析得到的模态以及对模型施加随机激励进行随机响应分析得到的应力分量均方根及功率谱密度；
(5)对步骤(4)中分析得到的应力分量均方根及功率谱密度，根据Mises应力功率谱密度估算方法计算Mises应力功率谱密度；
(6)根据步骤(5)中计算得到的Mises应力功率谱密度，由Dirlik公式计算出输送管的疲劳损伤；
(7)判断输送管疲劳损伤是否满足要求，若不满足要求，改进输送管道结构后返回步骤(1)，若满足要求，则结束；
所述步骤(1)中推进剂重量等效处理包括如下步骤：
输送管等效密度ρ通过公式得到：
<mrow>
<mi>&rho;</mi>
<mo>=</mo>
<msub>
<mi>&rho;</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<mfrac>
<mi>D</mi>
<mrow>
<mn>4</mn>
<mi>&delta;</mi>
</mrow>
</mfrac>
<msub>
<mi>&rho;</mi>
<mn>2</mn>
</msub>
</mrow>
其中，D为输送管路内径，δ为输送管路壁厚，ρ1为输送管材密度，ρ2为推进剂密度；
所述步骤(3)中判断输送管静强度是否满足要求包括如下步骤：
(a)输送管静强度安全因子η通过如下公式得到：
<mrow>
<mi>&eta;</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>1.67</mn>
<msub>
<mi>&sigma;</mi>
<mi>s</mi>
</msub>
</mrow>
<mi>&sigma;</mi>
</mfrac>
</mrow>
其中，σ为静强度分析得到的输送管Mises应力，σs为管路材料在设计温度下的屈服强度；
(b)若η≥1，则输送管静强度满足要求，反之若η＜1不满足要求；
所述步骤(5)中Mises应力功率谱密度估算方法如下：
输送管Mises应力谱密度G(σM)通过如下公式得到：
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mi>M</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>22</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msqrt>
<mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>11</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>22</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
<mo>+</mo>
<mn>3</mn>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>&sigma;</mi>
<mn>12</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
其中，G(σ11)、G(σ22)、G(σ12)为随机响应分析得到的输送管各应力分量的谱密度；
所述步骤(6)中疲劳损伤计算方法如下：
输送管疲劳损伤DL通过如下公式得到：
<mrow>
<msub>
<mi>D</mi>
<mi>L</mi>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>p</mi>
<mo>&rsqb;</mo>
<mi>T</mi>
</mrow>
<mi>C</mi>
</mfrac>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>&sigma;</mi>
<mi>b</mi>
</msub>
<mrow>
<msub>
<mi>&sigma;</mi>
<mi>b</mi>
</msub>
<mo>-</mo>
<mi>&sigma;</mi>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
<mi>m</mi>
</msup>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&infin;</mi>
</msubsup>
<msup>
<mi>S</mi>
<mi>m</mi>
</msup>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>S</mi>
</mrow>
其中，C、m为材料疲劳曲线参数，σb为管路材料在设计温度下的强度极限，T为输送管振动时间，p(S)为Mises应力概率密度，E[p]为单位之间内的峰值数；
所述单位之间内的峰值数E[p]通过如下公式得到：
<mrow>
<mi>E</mi>
<mo>&lsqb;</mo>
<mi>p</mi>
<mo>&rsqb;</mo>
<mo>=</mo>
<msup>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<msub>
<mi>M</mi>
<mn>4</mn>
</msub>
<msub>
<mi>M</mi>
<mn>2</mn>
</msub>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
</msup>
</mrow>
其中，n＝2,4，f1为随机激励的频率下限，f2为随机激励的频率上限；
所述Mises应力概率密度p(S)通过公式：
<mrow>
<mi>p</mi>
<mrow>
<mo>(</mo>
<mi>S</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mfrac>
<msub>
<mi>D</mi>
<mn>1</mn>
</msub>
<mi>q</mi>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<mi>z</mi>
<mi>q</mi>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<mfrac>
<mrow>
<msub>
<mi>D</mi>
<mn>2</mn>
</msub>
<mi>z</mi>
</mrow>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mfrac>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
<mrow>
<mn>2</mn>
<msup>
<mi>R</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
</mrow>
</msup>
<mo>+</mo>
<msub>
<mi>D</mi>
<mn>3</mn>
</msub>
<msup>
<mi>ze</mi>
<mrow>
<mo>-</mo>
<mfrac>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
<mn>2</mn>
</mfrac>
</mrow>
</msup>
</mrow>
<mrow>
<mn>2</mn>
<msqrt>
<msub>
<mi>M</mi>
<mn>0</mn>
</msub>
</msqrt>
</mrow>
</mfrac>
</mrow>
其中， D3＝1-D1-D2，
所述步骤(7)中判断输送管疲劳损伤是否满足要求方法如下：
若输送管疲劳损伤DL≥0.1，则输送管疲劳损伤不满足要求，反之若DL＜0.1满足要求；
所述的σ为静强度分析得到的输送管Mises应力；
所述的σ11、σ22、σ12为随机响应分析得到的输送管各应力分量。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310577106.3A CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310577106.3A CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103678775A CN103678775A (zh) | 2014-03-26 |
| CN103678775B true CN103678775B (zh) | 2017-08-29 |
Family
ID=50316311
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310577106.3A Active CN103678775B (zh) | 2013-11-18 | 2013-11-18 | 一种输送管动强度分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103678775B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105651478A (zh) * | 2015-12-15 | 2016-06-08 | 西安交通大学青岛研究院 | 一种基于振动信号测试零部件疲劳寿命的分析方法 |
| CN110114157B (zh) * | 2016-12-20 | 2020-11-27 | 杰富意钢铁株式会社 | 钢管的变形性能评价方法、钢管的制造方法 |
| CN116432335A (zh) * | 2023-03-03 | 2023-07-14 | 江阴市南方管件制造有限公司 | 一种基于有限元分析的化工用非标管件的制备工艺 |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101509855B (zh) * | 2009-03-17 | 2012-01-25 | 华东理工大学 | 一种高温管道损伤及寿命的在线预测方法 |
| CN102567567B (zh) * | 2011-11-15 | 2014-05-28 | 北京宇航系统工程研究所 | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 |
-
2013
- 2013-11-18 CN CN201310577106.3A patent/CN103678775B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN103678775A (zh) | 2014-03-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Szolwinski et al. | Mechanics of fretting fatigue crack formation | |
| Li et al. | Deep learning regression-based stratified probabilistic combined cycle fatigue damage evaluation for turbine bladed disks | |
| EP2983099B1 (en) | Method for determining reduction factor of axial load bearing capacity of a cylindrical shell structure in a rocket | |
| CN104200122B (zh) | 复杂焊接结构随机振动疲劳寿命预测方法 | |
| Mirzaei et al. | Experimental and numerical analysis of dynamic rupture of steel pipes under internal high-speed moving pressures | |
| CN107145641B (zh) | 叶片振动疲劳概率寿命预测方法 | |
| CN103926152A (zh) | 一种高温多轴谱载下低周蠕变-疲劳寿命评估方法 | |
| CN102567567A (zh) | 一种基于有限元分析的管路随机振动疲劳寿命分析方法 | |
| CN103678775B (zh) | 一种输送管动强度分析方法 | |
| CN105260574A (zh) | 一种基于临界面法疲劳破坏准则的高周多轴疲劳寿命预测方法 | |
| Bäcker et al. | Analysis of the deep rolling process on turbine blades using the FEM/BEM-coupling | |
| CN107977516A (zh) | 一种考虑多轴载荷非比例度的缺口件局部应力应变确定方法 | |
| Balcı et al. | Free vibration analysis of a laminated composite beam with various boundary conditions | |
| Zheng et al. | Autofrettage and shakedown analysis of strain-hardening cylinders under thermo-mechanical loadings | |
| Xu et al. | Assessment of fatigue life of remanufactured impeller based on FEA | |
| CN113919078A (zh) | 一种尺寸效应下耦合应力梯度的缺口结构疲劳分析方法 | |
| CN103837343B (zh) | 基于振动疲劳耦合分析的船艇轴系寿命预测方法 | |
| CN103674521A (zh) | 一种基于分工况设计的补偿器验证方法 | |
| Arai et al. | Investigation on method of elasto-plastic analysis for piping system (Benchmark analysis) | |
| Slima et al. | Fatigue analysis-based numerical design of stamping tools made of cast iron | |
| CN107748817A (zh) | 一种考虑非比例附加强化的高温多轴本构关系确定方法 | |
| Wang et al. | The Reliability Analysis for Pistons on Fracture Mechanics | |
| Zhu et al. | Fatigue Life Analysis of Liquid-Filled Pipeline Under the Combined Excitation of Random Vibration and Fluid-Structure Interaction | |
| CN116429362B (zh) | 一种用于发动机管路结构的疲劳试验方法 | |
| CN115358113B (zh) | 一种航空发动机用聚四氟乙烯软管脉冲寿命计算方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |