CN105760639B - 一种液体热膨胀引起的密闭工作容器应力分析方法 - Google Patents
一种液体热膨胀引起的密闭工作容器应力分析方法 Download PDFInfo
- Publication number
- CN105760639B CN105760639B CN201610320221.6A CN201610320221A CN105760639B CN 105760639 B CN105760639 B CN 105760639B CN 201610320221 A CN201610320221 A CN 201610320221A CN 105760639 B CN105760639 B CN 105760639B
- Authority
- CN
- China
- Prior art keywords
- working container
- liquid
- airtight
- volume
- deformation
- 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.)
- Expired - Fee Related
Links
- 239000007788 liquid Substances 0.000 title claims abstract description 43
- 238000004458 analytical method Methods 0.000 title claims abstract description 25
- 239000007787 solid Substances 0.000 claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 10
- 238000012805 post-processing Methods 0.000 claims abstract description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 230000003993 interaction Effects 0.000 claims description 3
- 238000012986 modification Methods 0.000 abstract description 5
- 230000004048 modification Effects 0.000 abstract description 5
- 238000013461 design Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 238000001816 cooling Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
一种液体热膨胀引起的密闭工作容器应力分析方法，包括以下步骤：建立密闭工作容器的模型；划分网格；去除容器中的液体，在密闭工作容器内表面施加均匀的法向压强；利用有限元求解得内表面上节点变形前、后的节点坐标；构造封闭的密闭工作容器的内表面并生成内表面实体模型，得到等效液体变形前、后的体积；计算液体积增量百分比，得到此时液体的温度；按有限元分析后处理方法，得到该温度下密闭工作容器各处的应力；依次对密闭工作容器内表面施加不同压强，重复上述步骤可以得到液体温度上升到不同值时的应力。对密闭工作容器进行应力分析以确定密闭工作容器各处的应力，最终确定密闭工作容器容易失效的位置，从而为结构修改提供依据。
Description
技术领域
本发明涉及一种液体热膨胀引起的密闭工作容器应力分析方法，属于应力分析技术领域。
背景技术
有限元分析是一种现代计算方法，它采用数学近似的方法对真实物理系统进行模拟，利用简单而又相互作用的元素，以有限数量的未知量去逼近无限未知量的真实系统。利用有限元分析对结构件进行详细的力学分析，可以获得构件在实际使用状态下尽可能真实的受力分析，从而可以在设计阶段对可能出现的各种问题进行安全评判，并进行相应的设计参数修改。
传统的流固耦合应力分析方法，流体区域的网格随固体边界的运动而动态更新，网格运算和网格质量控制过程异常复杂。为此，本发明提出了等效应力分析方法，提高网格划分质量，分析密闭容器热膨胀应力简单、有效。
发明内容
本发明的内容是提出一种液体热膨胀引起的密闭工作容器应力分析方法，对密闭工作容器进行应力分析以确定密闭工作容器各处的应力，最终确定密闭工作容器容易失效的位置，从而为结构修改提供依据。
为实现上述目的，本发明提供一种液体热膨胀引起的密闭工作容器应力分析方法，其特征在于，包括以下步骤：
S1建立密闭工作容器的三维模型；
S2划分网格；将简化后的密闭工作容器三维模型导入到有限元分析软件中，进行网格的划分，然后导出密闭工作容器内表面上所有节点的坐标；
S3去除密闭工作容器中的液体，在密闭工作容器内表面施加均匀的法向压强来模拟液体与密闭工作容器之间的相互作用；
S4在密闭工作容器内表面施加压强P1，利用有限元求解后导出内表面上所有节点的位移并叠加到变形前的内表面节点坐标上，从而得到内表面上节点变形后的节点坐标；
S5根据变形前的节点坐标和变形后的节点坐标，在软件中构造密闭工作容器的内表面并生成内表面实体模型，密闭工作容器内表面实体模型的体积等效为液体的体积，对内表面实体模型变形前和变形后的体积进行计算得到等效液体变形前的体积V0和变形后的体积V1；
S6将V0、V1、P1和K的值代入下式，得到液体热膨胀在不受约束情况下的体积增量百分比每个温度都对有与之相对应的液体体积增量百分比所以根据得此时的液体温度T1，
式中K为液体体积弹性模量，单位为Pa，体积弹性模量在不同温度下变化不大，视为固定值；
S7按有限元分析后处理方法，得到密闭工作容器各处的应力，即液体温度上升到T1时的应力；
S8依次对密闭工作容器内表面施加压强P2、P3、…、Pn，每次重复步骤S3～S7，即得到液体温度上升到不同值时的应力。
与现有技术相比，本发明具有如下有益效果：
(1)本发明适用于任何密闭工作容器由工作液热膨胀引起的应力分析，同时该方法也适用于密闭工作液热膨胀引起的电池冷却箱应力分析。
(2)本发明分析液体热膨胀引起的密闭工作容器应力较流固耦合热膨胀应力分析简单、实用，能够缩短修改结构的周期，降低研发的成本。
附图说明
图1为本发明中一种液体热膨胀引起的密闭工作容器应力分析方法的实现流程图
具体实施方式
以下结合附图和实施例对本发明作进一步详细说明。
S1建立密闭工作容器的三维模型。
S2划分网格。将简化后的密闭工作容器三维模型导入到有限元分析软件中，进行网格的划分，然后导出密闭工作容器内表面上所有节点的坐标。
S3去除密闭工作容器中的液体，在密闭工作容器内表面施加均匀的法向压强来模拟液体与密闭工作容器之间的相互作用。
S4在密闭工作容器内表面依次施加0Mpa、0.2Mpa、0.4Mpa、0.6Mpa、0.8Mpa、1.0Mpa的压强，利用有限元求解后导出内表面上所有节点的位移并叠加到变形前的内表面节点坐标上，从而得到内表面上节点变形后的节点坐标。
S5根据变形前的节点坐标和变形后的节点坐标，在软件中构造密闭工作容器的内表面并生成内表面实体模型，密闭工作容器内表面实体模型的体积可等效为液体的体积，对内表面实体模型变形前和变形后的体积进行计算可得到等效液体变形前的体积2269897mm3和变形后的体积2275160mm3、2280373mm3、2285896mm3、2290604mm3、2295397mm3。
S6工作液在不受约束的情况下，设20℃时体积为V0，即液体变形前的体积2269897mm3，若温度增加到T时，液体体积相对20℃的体积增量百分比每个温度都对有与之相对应的液体体积增量百分比将V0、V1、P1和K的值带入下式，得到液体热膨胀在不受约束情况下的体积增量百分比所以根据可得此时的液体温度T1。
式中K为液体体积弹性模量，单位为Pa，体积弹性模量在不同温度下变化不大，可视为固定值。因此，可以得到不同压力下冷却液温度为20℃、26℃、31℃、35℃、40℃、44℃。
S7按有限元分析一般后处理方法，得到密闭工作容器各处的应力，不同温度下易失效位置的最大应力分别为0Mpa、422Mpa、835Mpa、1263Mpa、1655Mpa、2060Mpa，根据密闭容器的材料可以得到其抗拉强度，由此可判断液体温度达到多少时密闭容器容易发生局部断裂，最终为密闭容器的结构设计和修改提供依据。本专利所采用的一种液体热膨胀引起的密闭工作容器应力分析方法得到的结果与一般流固耦合的热膨胀应力分析相比，误差在2％以内；但其在整个分析过程中能够节省30％的时间，能够缩短修改结构的周期，降低结构设计和修改的成本。
Claims (1)
1.一种液体热膨胀引起的密闭工作容器应力分析方法，其特征在于，包括以下步骤：
S1建立密闭工作容器的三维模型；
S2划分网格；将简化后的密闭工作容器三维模型导入到有限元分析软件中，进行网格的划分，然后导出密闭工作容器内表面上所有节点的坐标；
S3去除密闭工作容器中的液体，在密闭工作容器内表面施加均匀的法向压强来模拟液体与密闭工作容器之间的相互作用；
S4在密闭工作容器内表面施加压强P1，利用有限元求解后导出内表面上所有节点的位移并叠加到变形前的内表面节点坐标上，从而得到内表面上节点变形后的节点坐标；
S5根据变形前的节点坐标和变形后的节点坐标，在软件中构造密闭工作容器的内表面并生成内表面实体模型，密闭工作容器内表面实体模型的体积等效为液体的体积，对内表面实体模型变形前和变形后的体积进行计算得到等效液体变形前的体积V0和变形后的体积V1；
S6将V0、V1、P1和K的值代入下式，得到液体热膨胀在不受约束情况下的体积增量百分比每个温度都对有与之相对应的液体体积增量百分比所以根据得此时的液体温度T1，
式中K为液体体积弹性模量，单位为Pa，体积弹性模量在不同温度下变化不大，视为固定值；
S7按有限元分析后处理方法，得到密闭工作容器各处的应力，即液体温度上升到T1时的应力；
S8依次对密闭工作容器内表面施加压强P2、P3、…、Pn，每次重复步骤S3～S7，即得到液体温度上升到不同值时的应力。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610320221.6A CN105760639B (zh) | 2016-05-15 | 2016-05-15 | 一种液体热膨胀引起的密闭工作容器应力分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610320221.6A CN105760639B (zh) | 2016-05-15 | 2016-05-15 | 一种液体热膨胀引起的密闭工作容器应力分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105760639A CN105760639A (zh) | 2016-07-13 |
| CN105760639B true CN105760639B (zh) | 2018-12-18 |
Family
ID=56324064
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610320221.6A Expired - Fee Related CN105760639B (zh) | 2016-05-15 | 2016-05-15 | 一种液体热膨胀引起的密闭工作容器应力分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105760639B (zh) |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105043647A (zh) * | 2015-04-23 | 2015-11-11 | 成都柏森松传感技术有限公司 | 一种密闭容器内液体压力监测系统 |
| CN105352816A (zh) * | 2015-11-26 | 2016-02-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料结构失效预测分析方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2526501A4 (en) * | 2010-01-21 | 2017-09-20 | Firehole Technologies | Automated method to determine composite material constituent properties |
-
2016
- 2016-05-15 CN CN201610320221.6A patent/CN105760639B/zh not_active Expired - Fee Related
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105043647A (zh) * | 2015-04-23 | 2015-11-11 | 成都柏森松传感技术有限公司 | 一种密闭容器内液体压力监测系统 |
| CN105352816A (zh) * | 2015-11-26 | 2016-02-24 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种复合材料结构失效预测分析方法 |
Non-Patent Citations (1)
| Title |
|---|
| 基于ANSYS的某密闭容器静态应力分析;牛碧凯 等;《机械工程与自动化》;20150430(第2期);89-90,93 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105760639A (zh) | 2016-07-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107391885B (zh) | 基于有限体积法的剪切滑移动网格方法 | |
| CN106485032A (zh) | 一种计算叶盘振动可靠性的双重极值响应面法 | |
| CN105335568B (zh) | 一种基于有限元技术考虑热膨胀的超塑成形模具设计方法 | |
| CN101767185B (zh) | 一种基于定量设置反变形量的设计铸件模型的方法 | |
| CN107563066A (zh) | 一种基于数值模拟的回弹过程模具补偿方法 | |
| CN115422654B (zh) | 基于cfd/csd技术的跨/超声速飞行器高效高精度颤振时域分析方法 | |
| CN105631125A (zh) | 基于降阶模型的气动-热-结构耦合分析方法 | |
| CN106873568A (zh) | 基于h无穷鲁棒未知输入观测器的传感器故障诊断方法 | |
| CN106909719A (zh) | 集成多元线性回归算法 | |
| CN110826208B (zh) | 一种气动数值模拟加速收敛方法 | |
| CN105760639B (zh) | 一种液体热膨胀引起的密闭工作容器应力分析方法 | |
| CN107563071A (zh) | 一种综合环境应力的仿真方法 | |
| CN105224726B (zh) | 结构网格动网格技术用于非结构网格流场求解器的方法 | |
| CN104699901A (zh) | 基于差量采样解的GappyPOD翼型反设计方法 | |
| Yao et al. | Superconvergence analysis of two-grid FEM for Maxwell’s equations with a thermal effect | |
| CN105956347A (zh) | 纤维织物增强陶瓷基复合材料力学行为的仿真预测方法 | |
| CN118734457A (zh) | 一种基于多保真流形神经算子的复合材料构件固化状态场预测方法 | |
| CN106407547A (zh) | 针对各向异性材料铸造残余应力的数值模拟方法 | |
| CN105631061B (zh) | 一种复杂点式连接结构的动力学建模方法 | |
| CN107944116B (zh) | 一种针对时变结构的瞬态能量响应高效预示方法 | |
| CN110633556A (zh) | 一种陶瓷基复合材料流固耦合响应计算方法 | |
| CN102411656A (zh) | 一种模拟接触的有限元建模方法 | |
| CN105787999A (zh) | 一种火箭总体骨架曲面搭建系统及方法 | |
| CN104239696B (zh) | 一种约束状态下构件焊接变形的预测方法 | |
| Cook et al. | Applications of the CESE method in LS-DYNA |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20181218 Termination date: 20210515 |