CN102902859A - 基于相变模型的动静压滑动轴承流固耦合设计方法 - Google Patents
基于相变模型的动静压滑动轴承流固耦合设计方法 Download PDFInfo
- Publication number
- CN102902859A CN102902859A CN2012103811235A CN201210381123A CN102902859A CN 102902859 A CN102902859 A CN 102902859A CN 2012103811235 A CN2012103811235 A CN 2012103811235A CN 201210381123 A CN201210381123 A CN 201210381123A CN 102902859 A CN102902859 A CN 102902859A
- Authority
- CN
- China
- Prior art keywords
- model
- fluid
- solid
- grid
- calculation
- 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.)
- Granted
Links
- 239000007787 solid Substances 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000010168 coupling process Methods 0.000 title claims abstract description 35
- 230000008878 coupling Effects 0.000 title claims abstract description 34
- 238000005859 coupling reaction Methods 0.000 title claims abstract description 34
- 230000008859 change Effects 0.000 title abstract description 7
- 238000004364 calculation method Methods 0.000 claims abstract description 54
- 239000012530 fluid Substances 0.000 claims abstract description 52
- 230000008676 import Effects 0.000 claims abstract description 14
- 230000008569 process Effects 0.000 claims abstract description 8
- 238000005457 optimization Methods 0.000 claims abstract description 6
- 241000682719 Adina Species 0.000 claims description 9
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000002131 composite material Substances 0.000 claims description 6
- 238000006073 displacement reaction Methods 0.000 claims description 6
- 230000010354 integration Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims description 3
- 238000005059 solid analysis Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 230000002950 deficient Effects 0.000 claims 1
- 230000007480 spreading Effects 0.000 claims 1
- 230000009466 transformation Effects 0.000 claims 1
- 230000003068 static effect Effects 0.000 abstract description 11
- 239000003921 oil Substances 0.000 description 18
- 230000002706 hydrostatic effect Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000012805 post-processing Methods 0.000 description 2
- 230000008439 repair process Effects 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000005489 elastic deformation Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- FKLFBQCQQYDUAM-UHFFFAOYSA-N fenpiclonil Chemical compound ClC1=CC=CC(C=2C(=CNC=2)C#N)=C1Cl FKLFBQCQQYDUAM-UHFFFAOYSA-N 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000010687 lubricating oil Substances 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
Images
Landscapes
- Sliding-Contact Bearings (AREA)
- Magnetic Bearings And Hydrostatic Bearings (AREA)
Abstract
一种基于相变模型的动静压滑动轴承流固耦合设计方法，首先根据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体油膜以及固体的几何模型和网格计算模型。之后将网格导入ansa软件进行网格优化。然后将流体和固体计算模型分别导入CAE软件ADINA中相应的流体与结构计算模块，设定边界条件及求解控制参数后，进行流固耦合数值计算。在数值计算过程中加入相变公式，计算运转过程中产生的气泡问题。根据计算结果分析油膜及轴承的性能，之后调整腔结构及几何模型参数，得到满意的计算结果，从而为轴承结构设计提供依据。
Description
技术领域
本发明涉及一种动静压滑动轴承的设计方法，特别涉及一种利用流固耦合数值计算技术实现动静压滑动轴承性能参数优化设计方法。
背景技术
动静压滑动轴承系统能够在一定的工况下保持耦合系统的稳定运行，从而使主轴可以稳定的运转，达到承压、润滑的目的。尤其对于速度高，弹性变形大的高速动静压轴承，考虑耦合系统影响的流固耦合方法更能够真实反应轴承的运转工况。
目前大部分研究润滑油膜使用计算流体动力学（CFD）方法和雷诺方程法对滑动轴承进行计算，但是这些方法不考虑耦合系统的影响，不考虑轴颈和轴承对油膜压力的影响，计算出来的各项参数不能真实的反应轴承的运转工况，对轴承结构的优化设计有一定的影响，因此数值计算使用流固耦合（FSI）的方法是至关重要的。这种方法能够在转速和载荷的作用下自动确定偏心位置，而不需要事先给定偏心位置进行求解，因此，FSI方法更能真实的反应轴承的运转情况。
发明内容
本发明的目的在于克服传统方法无法真实反应轴承的运转状态，无法获取真实的轴承性能指标，提供了一种基于相变模型的动静压滑动轴承流固耦合设计方法。
为达到上述目的，本发明采用的技术方案是：
1）首先根据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体油膜以及固体的几何模型和网格计算模型；
2）之后将网格计算模型导入ansa软件进行网格优化；
3）然后将流体油膜以及固体的网格计算模型分别导入CAE软件ADINA中相应的流体与结构计算模块，设定边界条件，进行流固耦合数值计算，在数值计算过程中加入相变公式，计算运转过程中产生的气泡问题；
4）根据计算结果分析油膜及轴承的性能，之后调整腔结构及几何模型参数，得到所要设计的动静压滑动轴承。
本发明包括以下步骤：
第一步，几何模型的建立
据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体域几何模型即流体油膜几何模型以及固体域几何模型；
第二步，模型离散
在gambit软件中对建立的流体和固体几何模型进行几何离散，划分网格得到计算网格模型；
第三步，网格计算模型的优化
将计算网格模型导入ansa软件中，利用ansa软件的几何处理和网格处理功能，对流体和固体网格进行优化，对几何模型和网格模型存在的缺陷进行修补；
第四步，模型导入
将流体及固体的网格计算模型分别导入CAE计算软件ADINA相应的流体与结构计算模块，并分别设置流体、固体的材料参数，流体参数为不可压缩粘性流体，设置湍流模型；
第五步，边界条件设置
模型导入ADINA之后，在ADINA软件中设置流体域和固体域的边界条件；
在流体域中，设置7个边界条件，分别是两端进油孔的进油压力边界、前后圆柱端面为出油口边界、施加在内表面上的转速边界、内表面的流固耦合边界、内表面的滑移边界、外表面的流固耦合边界、整个流体域的相变边界；
在固体域中，设置5个边界条件，分别是轴瓦外边面的固定边界、轴瓦内表面的流固耦合边界、轴颈外表面的流固耦合边界、在轴颈上的动载载荷边界、轴颈的约束边界；
第六步，流固耦合数值计算
设置完边界条件之后，在ADINA中根据已知的边界条件和模型控制参数进行流固耦合计算；
设置计算过程中的模型控制参数，流固耦合采用迭代耦合分析方法，控制方程求解采用完全牛顿迭代方法，设置位移与压力收敛判据，流体控制方程求解采用有限体积方法，格式选用二阶空间离散，时间积分采用二阶composite积分格式，求解采用simple算法，设置变量收敛判据，固体分析采用隐式动力分析方法，有限元方法离散，二阶精度Bathe composite时间积分格式，完全牛顿迭代法计算，设置位移收敛判据；
第七步，后处理
计算完成后统计油膜压力值，轴瓦应力值，绘制轴心轨迹曲线以及气穴面积曲线图，并与设计要求进行对比：若计算结果与设计要求偏差在设计精度要求允许范围内，则完成计算；否则返回第一步,修改模型参数重新计算直至满足设计要求，并记录已完成计算的模型参数及计算结果，为修改模型参数提供依据。
本发明首先根据实际轴承尺寸要求，采用gambit软件建立两腔结构的动静压滑动轴承的流体油膜以及固体的几何模型和网格计算模型。之后将网格导入ansa软件进行网格优化。然后将流体和固体计算模型分别导入CAE软件ADINA中相应的流体与结构计算模块，设定边界条件及求解控制参数后，进行流固耦合数值计算。在数值计算过程中加入相变公式，计算运转过程中产生的气泡问题。根据计算结果分析油膜及轴承的性能，之后调整腔结构及几何模型参数，得到满意的计算结果，从而为轴承结构设计提供依据。
附图说明
图1是本发明技术路线图；
图2是本发明动静压滑动轴承流体域计算几何模型图（一半）；
图3是本发明动静压滑动轴承固体域计算几何模型图（一半）；
图4是本发明在gambit中建立的网格计算模型（一半）；
图5是本发明在ansa中的优化网格模型（一半）；
图6是本发明动静压滑动轴承流固耦合导入ADINA中的流体域计算网格模型图（一半）；
图7是本发明动静压滑动轴承流固耦合导入ADINA中的固体域计算网格模型图（一半）；
具体实施方式
下面结合附图对本发明作进一步详细说明。
参见图1，本发明的过程如下：
1）首先根据动静压滑动轴承设计要求，采用gambit软件进行结构设计，建立动静压滑动轴承的流体区域几何模型及固体域几何模型；
2）之后在gambit软件中继续进行计算网格模型设计，建立网格单元；
3）将网格计算模型导入ansa中进行网格优化设计，对网格的缺陷和不足进行优化设计；
4）然后将优化过的网格计算模型导入CAE软件ADINA中，设定与实际工况相同的边界条件；
5）设置好所有的边界条件之后，输入求解控制参数，进行流固耦合数值计算；
6）根据计算结果调整几何模型参数，若符合设计要求完成设计，若不符合设计要求转步骤1）重新对结构设计进行优化，修改几何模型。
参见图2-图7，设计方法的具体过程如下：
第一步，几何模型的建立
参见图2、3，根据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体域几何模型以及固体域几何模型；
第二步，模型离散
参见图4，在gambit软件中对建立的几何模型进行几何离散，划分网格得到计算网格模型；
第三步，网格计算模型的优化
参见图5，将计算网格模型导入ansa软件中，利用ansa软件的几何处理和网格处理功能，对流体和固体网格进行优化，对几何模型和网格模型存在的缺陷进行修补；
第四步，模型导入
参见图6、7，将流体及固体的网格计算模型分别导入CAE计算软件ADINA相应的流体与结构计算模块，并分别设置流体、固体的材料参数，流体参数为不可压缩粘性流体，设置湍流模型；
第五步，边界条件设置
模型导入ADINA之后，在ADINA软件中设置流体域和固体域的边界条件；
在流体域中，设置7个边界条件，分别是两端进油孔的进油压力边界、前后圆柱端面为出油口边界、施加在内表面上的转速边界、内表面的流固耦合边界、内表面的滑移边界、外表面的流固耦合边界、整个流体域的相变边界；
在固体域中，设置5个边界条件，分别是轴瓦外边面的固定边界、轴瓦内表面的流固耦合边界、轴颈外表面的流固耦合边界、在轴颈上的动载载荷边界、轴颈的约束边界；
第六步，流固耦合数值计算
设置完边界条件之后，在ADINA中根据已知的边界条件和模型控制参数进行流固耦合计算；
设置计算过程中的模型控制参数，流固耦合采用迭代耦合分析方法，控制方程求解采用完全牛顿迭代方法，设置位移与压力收敛判据，流体控制方程求解采用有限体积方法，格式选用二阶空间离散，时间积分采用二阶composite积分格式，求解采用simple算法，设置变量收敛判据，固体分析采用隐式动力分析方法，有限元方法离散，二阶精度Bathe composite时间积分格式，完全牛顿迭代法计算，设置位移收敛判据；
第七步，后处理
计算完成后统计油膜压力值，轴瓦应力值，绘制轴心轨迹曲线以及气穴面积曲线图，并与设计要求进行对比：若计算结果与设计要求偏差在设计精度要求允许范围内，则完成计算；否则返回第一步,修改模型参数重新计算直至满足设计要求，并记录已完成计算的模型参数及计算结果，为修改模型参数提供依据。
Claims (2)
1.一种基于相变模型的动静压滑动轴承流固耦合设计方法，其特征在于：
1）首先根据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体油膜以及固体的几何模型和网格计算模型；
2）之后将网格计算模型导入ansa软件进行网格优化；
3）然后将流体油膜以及固体的网格计算模型分别导入CAE软件ADINA中相应的流体与结构计算模块，设定边界条件，进行流固耦合数值计算，在数值计算过程中加入相变公式，计算运转过程中产生的气泡问题；
4）根据计算结果分析油膜及轴承的性能，之后调整腔结构及几何模型参数，得到所要设计的动静压滑动轴承。
2.根据权利要求1所述的基于相变模型的动静压滑动轴承流固耦合设计方法，其特征在于包括以下步骤：
第一步，几何模型的建立
据实际轴承尺寸要求，采用gambit软件建立含油腔结构的动静压滑动轴承的流体域几何模型即流体油膜几何模型以及固体域几何模型；
第二步，模型离散
在gambit软件中对建立的流体和固体几何模型进行几何离散，划分网格得到计算网格模型；
第三步，网格计算模型的优化
将计算网格模型导入ansa软件中，利用ansa软件的几何处理和网格处理功能，对流体和固体网格进行优化，对几何模型和网格模型存在的缺陷进行修补；
第四步，模型导入
将流体及固体的网格计算模型分别导入CAE计算软件ADINA相应的流体与结构计算模块，并分别设置流体、固体的材料参数，流体参数为不可压缩粘性流体，设置湍流模型；
第五步，边界条件设置
模型导入ADINA之后，在ADINA软件中设置流体域和固体域的边界条件；
在流体域中，设置7个边界条件，分别是两端进油孔的进油压力边界、前后圆柱端面为出油口边界、施加在内表面上的转速边界、内表面的流固耦合边界、内表面的滑移边界、外表面的流固耦合边界、整个流体域的相变边界；
在固体域中，设置5个边界条件，分别是轴瓦外边面的固定边界、轴瓦内表面的流固耦合边界、轴颈外表面的流固耦合边界、在轴颈上的动载载荷边界、轴颈的约束边界；
第六步，流固耦合数值计算
设置完边界条件之后，在ADINA中根据已知的边界条件和模型控制参数进行流固耦合计算；
设置计算过程中的模型控制参数，流固耦合采用迭代耦合分析方法，控制方程求解采用完全牛顿迭代方法，设置位移与压力收敛判据，流体控制方程求解采用有限体积方法，格式选用二阶空间离散，时间积分采用二阶composite积分格式，求解采用simple算法，设置变量收敛判据，固体分析采用隐式动力分析方法，有限元方法离散，二阶精度Bathe composite时间积分格式，完全牛顿迭代法计算，设置位移收敛判据；
第七步，后处理
计算完成后统计油膜压力值，轴瓦应力值，绘制轴心轨迹曲线以及气穴面积曲线图，并与设计要求进行对比：若计算结果与设计要求偏差在设计精度要求允许范围内，则完成计算；否则返回第一步,修改模型参数重新计算直至满足设计要求，并记录已完成计算的模型参数及计算结果，为修改模型参数提供依据。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210381123.5A CN102902859B (zh) | 2012-10-09 | 2012-10-09 | 基于相变模型的动静压滑动轴承流固耦合设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201210381123.5A CN102902859B (zh) | 2012-10-09 | 2012-10-09 | 基于相变模型的动静压滑动轴承流固耦合设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN102902859A true CN102902859A (zh) | 2013-01-30 |
| CN102902859B CN102902859B (zh) | 2014-11-05 |
Family
ID=47575089
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201210381123.5A Expired - Fee Related CN102902859B (zh) | 2012-10-09 | 2012-10-09 | 基于相变模型的动静压滑动轴承流固耦合设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN102902859B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103235857A (zh) * | 2013-04-28 | 2013-08-07 | 昆明学院 | 中空轴式静压轴承多场耦合仿真分析方法 |
| CN104794356A (zh) * | 2015-04-28 | 2015-07-22 | 河海大学 | 一种高耸升船机塔柱结构抗风分析方法 |
| CN106383930A (zh) * | 2016-08-31 | 2017-02-08 | 中国石油大学(华东) | 一种尾轴承‑转子系统多重流固耦合计算方法 |
| CN108319737A (zh) * | 2017-01-17 | 2018-07-24 | 沈阳工业大学 | 高速列车铝合金齿轮箱的流场温度场耦合仿真分析方法 |
| CN110705181A (zh) * | 2019-10-13 | 2020-01-17 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
| CN111222291A (zh) * | 2019-12-30 | 2020-06-02 | 太原重工股份有限公司 | 油膜流体域离散化方法 |
| CN111666644A (zh) * | 2020-06-18 | 2020-09-15 | 哈尔滨工业大学 | 一种参数化cfd气体静压轴承性能分析方法 |
| CN111751220A (zh) * | 2020-06-22 | 2020-10-09 | 哈尔滨工业大学 | 一种考虑流固耦合的气体静压轴承性能计算方法 |
| CN111967107A (zh) * | 2020-08-21 | 2020-11-20 | 西安交通大学 | 内反馈动静压滑动轴承油膜压力场的nurbs等几何求解方法 |
| CN114169259A (zh) * | 2021-11-15 | 2022-03-11 | 广东核电合营有限公司 | 主泵静压轴封在极限工况运行可靠性的计算评价方法 |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201858256U (zh) * | 2010-11-15 | 2011-06-08 | 南宁广发重工集团有限公司 | 滑动轴承 |
-
2012
- 2012-10-09 CN CN201210381123.5A patent/CN102902859B/zh not_active Expired - Fee Related
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN201858256U (zh) * | 2010-11-15 | 2011-06-08 | 南宁广发重工集团有限公司 | 滑动轴承 |
Non-Patent Citations (3)
| Title |
|---|
| 孟晶 等: "基于FLUENT的液体动静压轴承油膜特性的分析", 《液压与气动》 * |
| 林起崟 等: "高速滑动轴承流固耦合传热及流场分析", 《润滑与密封》 * |
| 牛犇 等: "水润滑轴承螺旋密封性能的仿真分析", 《润滑与密封》 * |
Cited By (15)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103235857A (zh) * | 2013-04-28 | 2013-08-07 | 昆明学院 | 中空轴式静压轴承多场耦合仿真分析方法 |
| CN104794356A (zh) * | 2015-04-28 | 2015-07-22 | 河海大学 | 一种高耸升船机塔柱结构抗风分析方法 |
| CN104794356B (zh) * | 2015-04-28 | 2017-07-18 | 河海大学 | 一种高耸升船机塔柱结构抗风分析方法 |
| CN106383930A (zh) * | 2016-08-31 | 2017-02-08 | 中国石油大学(华东) | 一种尾轴承‑转子系统多重流固耦合计算方法 |
| CN106383930B (zh) * | 2016-08-31 | 2019-06-28 | 中国石油大学(华东) | 一种尾轴承-转子系统多重流固耦合计算方法 |
| CN108319737B (zh) * | 2017-01-17 | 2021-04-20 | 沈阳工业大学 | 高速列车铝合金齿轮箱的流场温度场耦合仿真分析方法 |
| CN108319737A (zh) * | 2017-01-17 | 2018-07-24 | 沈阳工业大学 | 高速列车铝合金齿轮箱的流场温度场耦合仿真分析方法 |
| CN110705181A (zh) * | 2019-10-13 | 2020-01-17 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
| CN110705181B (zh) * | 2019-10-13 | 2022-06-24 | 重庆交通大学 | 基于卷积长短时记忆循环神经网络的滚动轴承剩余寿命预测方法 |
| CN111222291A (zh) * | 2019-12-30 | 2020-06-02 | 太原重工股份有限公司 | 油膜流体域离散化方法 |
| CN111222291B (zh) * | 2019-12-30 | 2023-03-10 | 太原重工股份有限公司 | 油膜流体域离散化方法 |
| CN111666644A (zh) * | 2020-06-18 | 2020-09-15 | 哈尔滨工业大学 | 一种参数化cfd气体静压轴承性能分析方法 |
| CN111751220A (zh) * | 2020-06-22 | 2020-10-09 | 哈尔滨工业大学 | 一种考虑流固耦合的气体静压轴承性能计算方法 |
| CN111967107A (zh) * | 2020-08-21 | 2020-11-20 | 西安交通大学 | 内反馈动静压滑动轴承油膜压力场的nurbs等几何求解方法 |
| CN114169259A (zh) * | 2021-11-15 | 2022-03-11 | 广东核电合营有限公司 | 主泵静压轴封在极限工况运行可靠性的计算评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN102902859B (zh) | 2014-11-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102902859B (zh) | 基于相变模型的动静压滑动轴承流固耦合设计方法 | |
| CN103226633B (zh) | 一种静压油垫流场特性仿真优化方法 | |
| CN103366067B (zh) | 仿真预测液力变矩器性能的方法 | |
| CN104298869B (zh) | 一种弹性水翼的流固耦合特性数值预测方法 | |
| Castilla et al. | Three-dimensional numerical simulation of an external gear pump with decompression slot and meshing contact point | |
| CN111967107B (zh) | 内反馈动静压滑动轴承油膜压力场的nurbs等几何求解方法 | |
| CN105022900B (zh) | 基于热固耦合分析的重型数控立车静压转台结构优化方法 | |
| CN103186698B (zh) | 一种重型机床静压转台动、静态性能仿真优化方法 | |
| CN102141084B (zh) | 基于膜厚可变的静压推力轴承间隙油膜温度与厚度关系数值模拟方法 | |
| CN102799705A (zh) | 负载敏感控制轴向柱塞泵虚拟样机的建立方法 | |
| CN107526914B (zh) | 基于结构化动网格的可倾瓦滑动轴承变流域流场计算方法 | |
| CN102880757A (zh) | 一种基于流固耦合数值计算的微灌压力调节器设计方法 | |
| CN104091021A (zh) | 一种计算倾斜时油垫承载力的方法 | |
| CN105095536B (zh) | 一种考虑表面形貌的静压油垫流场特性数值仿真方法 | |
| CN102693339A (zh) | 一种斯特林制冷机曲柄连杆机构综合分析方法 | |
| CN104091004B (zh) | 一种基于改进的粒子群算法的定量式静压转台优化设计方法 | |
| CN113343462B (zh) | 基于高阶等几何的多油腔动静压滑动轴承油膜特性仿真方法 | |
| Tang et al. | A local and fast interpolation method for mesh deformation | |
| CN111666644A (zh) | 一种参数化cfd气体静压轴承性能分析方法 | |
| CN113158369B (zh) | 静压推力轴承油垫封油边的油膜流量模拟监测方法 | |
| Cheng | Finite element analysis and structural optimization of the box on the ANSYS workbench | |
| Hui et al. | Investigation on numerical computational methodology for the flow characteristics of aviation piston pump | |
| CN111222291B (zh) | 油膜流体域离散化方法 | |
| Zhou et al. | A well-balanced stable grp scheme for shallow water equations for adaptive unstructured triangular meshes | |
| Li et al. | Analysis of turbine fluid dynamics based on FSI |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141105 Termination date: 20181009 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |