CN104091022A - 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法 - Google Patents

一种速度滑移工况下静压主轴系统流固耦合仿真分析方法 Download PDF

Info

Publication number
CN104091022A
CN104091022A CN201410328613.8A CN201410328613A CN104091022A CN 104091022 A CN104091022 A CN 104091022A CN 201410328613 A CN201410328613 A CN 201410328613A CN 104091022 A CN104091022 A CN 104091022A
Authority
CN
China
Prior art keywords
mrow
bearing
mfrac
finite element
rotor
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.)
Pending
Application number
CN201410328613.8A
Other languages
English (en)
Inventor
陈东菊
边艳华
周帅
范晋伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201410328613.8A priority Critical patent/CN104091022A/zh
Publication of CN104091022A publication Critical patent/CN104091022A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Magnetic Bearings And Hydrostatic Bearings (AREA)

Abstract

本发明涉及一种速度滑移工况下静压主轴系统流固耦合仿真分析方法,该方法包括:选择Navier速度滑移模型模拟轴承流固界面速度滑移现象,建立考虑速度滑移情况下静压轴承油膜压力分布的有限差分模型,求解油膜的压力分布;借助ansys有限元软件建立转子和轴瓦及轴承座的有限元模型,应用APDL参数化有限元分析技术对转子和轴瓦及轴承座有限元模型求解;提取转子和轴瓦及轴承座有限元模型相关节点的变形值,并将数据导入到油膜厚度的计算中,实现对油膜压力的重新求解;重复上述过程直至转子和轴瓦及轴承座的变形达到稳定。本发明考虑了固流界面速度滑移的影响,实现了静压轴承-转子系统的流体-结构之间的双向耦合。

Description

一种速度滑移工况下静压主轴系统流固耦合仿真分析方法
技术领域
本发明流体动力学计算领域,涉及一种适用于静压轴承-转子系统流体-结构之间的双向耦合仿真方法。
背景技术
静压轴承以支撑精度高等优点,广泛应用于机械制造和仪器制造等领域。在静压轴承中,润滑介质在轴承间隙中的流动已经进入到了微尺度流动的研究领域。然而流体在微尺度下流动时将会表现出异于宏观尺度的流动的性质,流体在微尺度下会表现出来大表面积体积比、梯度参数效应及界面效应等问题。将微尺度的性质引入到静压轴承的研究中,从而使油膜性能的研究更加准确。
静压轴承性能分析在很大程度上依赖于仿真方法,仿真方法的优劣直接决定静压轴承分析的结果的准确性。但传统的仿真方法是对静压轴承工作环境做出在一定程度上合理的假设,传统的方法忽略了在静压轴承流固界面的速度滑移现象,而且没有考虑轴承和转子之间的相互作用,没有对轴承和转子之间进行流固耦合分析,速度滑移和流固界面的相互耦合作用在一定程度上影响着油膜行为,油膜行为对决定静压轴承性能起着关键作用。在以往的一些研究中单纯研究速度滑移对轴承承载力的影响,静压轴承受滑移的影响承载力会在一定程度上得到提高,单纯研究流固界面流体和固体的相互耦合作用对静压轴承的影响,轴承的承载力在一定程度上会有所降低,而在静压轴承的实际工况作用下,两种作用都会对静压轴承产生一定的影响,而两种影响作用不能简单的相互抵消,最终的影响结果和很多因素有关,例如轴瓦及转子材料、油膜温度等,这必然会使传统的仿真方法的结果和实际工作结果之间产生一定的误差。
发明内容
针对现有技术中存在的上述问题,本发明提供一种考虑速度滑移影响下静压轴承-转子系统的流体-结构耦合仿真方法,考虑了固-流界面速度滑移的影响,具有计算简洁、通用性强、可人工控制耦合迭代过程等优点,实现了真正意义上的静压轴承-转子系统的流体-结构之间的双向耦合。
一种速度滑移工况下静压主轴系统流固耦合仿真分析方法,包括以下步骤:
(1)根据静压轴承的工况选择Navier速度滑移模型,模拟轴承-转子系统中流固界面上出现的速度滑移现象。
Navier模型假设流体的滑移速度vs与流体在表面的剪切力成正比,即:
v s = b ∂ u 1 ∂ z | wall - - - ( 1 )
式中,u1是流体沿固体表面的流动速度,Z是沿界面法向方向坐标,wall表示流体和固体相互作用的界面位置。b称之为滑移系数,是虚构固体表面与实际界面的距离,如图1所示。
(2)建立考虑速度滑移情况下静压轴承油膜压力分布的有限差分模型,根据油膜结构划分合理网格。
考虑速度滑移条件下的油膜压力分布符合修正后的N-S方程为:
∂ ∂ x { [ - h 3 12 η - hb bx 2 η ( h - h * ) - hh * b ax 2 η ] ∂ p ∂ x + hu 2 } + ∂ ∂ z { ( - h 3 12 η - h 2 b az 2 η ) ∂ p ∂ z } = 0 - - - ( 2 )
式中,u为旋转主轴表面的线速度;bax、baz分别为轴瓦表面周向、轴向滑移系数;bbx为轴颈表面周向滑移系数;h*为油膜周向方向剪切力为零时的油膜厚度;p为油膜压力;η为润滑油的动力粘度。
采用有限差分方法将方程(2)线性化处理,得到如下线性化方程:
P i , j = 1 2 A i , j + 2 C i , j [ ( A i , j + B i , j ) P i + 1 , j + ( A i , j - B i , j ) P i - 1 , j + C i , j P i , j + 1 + C i , j P i , j - 1 - D i , j P i , j ] - - - ( 3 )
其中,
C i , j = β 2 H i , j 3 ( 1 + 6 b 0 ) Δλ 2
β = 2 R L
式中,b0为滑移系数,R为轴承半径,L为轴承长度,Δλ为轴承轴向方向网格长度,为轴承周向方向网格长度,(i,j)为油膜位置坐标,Hi,j表示(i,j)处无量纲油膜厚度。
应用超松弛迭代法将线性化方程(3)在matlab软件平台上进行数值求解,即得考虑速度滑移影响下油膜压力分布。
(3)用matlab语言编写程序将计算所得的压力分布数据导出,将数据文件存储成txt格式,便于将此计算数据在步骤(5)中作为压力边界条件加载到转子和轴瓦及轴承座有限元模型中。
(4)借助ansys有限元软件建立转子和轴瓦及轴承座的有限元模型,控制有限元模型的网格划分形式与油膜的网格划分一致,以保证数据交换过程的顺利完成。
在转子模型建立过程中,采用自下而上的建模方法,先创建关键点,依次构成线和面,用lsize控制转子有限元模型的网格划分,将面绕轴线旋转出周向规则的网格,然后在转子模型的表面建立表面效应单元,用于加载与表面相互垂直的压力载荷,定义转子有限元模型的材料属性,转子的有限元模型建立完成。
在轴瓦及轴承座模型的建立中,静压腔的深度与轴承的尺寸相比很小,所以采用有限元方法对轴瓦建模时忽略轴瓦表面的静压腔,有限元模型的建立过程与转子模型的建立过程相同,但需根据实际情况分别对轴瓦及轴承座定义有限元模型的材料属性。
(5)应用APDL参数化有限元分析技术编写相应程序实现matlab与ansys之间的数据交换,编写APDL程序将matlab数据读入ansys,将油膜压力分布的matlab计算结果作为边界条件加载到转子和轴瓦及轴承座有限元模型中,实现对转子和轴瓦及轴承座有限元模型的求解。
(6)应用APDL参数化有限元分析技术提取转子和轴瓦及轴承座有限元模型相关节点的变形值。有限元模型在某节点的变形可以分解为X,Y,Z三个方向,节点在X,Y方向上的变形将影响到油膜厚度的变化。
(7)将ansys导出的节点变形数据导入到matlab油膜厚度的计算中,实现对油膜压力分布的重新求解。
(8)重复步骤(3)~(7)直至转子和轴瓦及轴承座的变形达到稳定,实现真正意义上静压轴承-转子系统的流体-结构之间的双向耦合。
与现有技术相比,本发明具有以下优点:
本发明考虑了传统设计中忽略的固液界面速度滑移现象,将油膜分析尺度引入到微尺度的研究领域;并将流体和固体分析过程进行耦合,使得分析结果更加准确,更加有利于指导静压轴承应用于实践中。应用本发明仿真方法和传统仿真方法在不忽略速度滑移和不进行流固耦合分析的情况下计算某静压轴承的承载力,并与实验方法得到的承载力进行比较,本发明仿真方法计算得到承载力误差为2.97%,而传统仿真方法得到的误差分别为10.89%和38.6%。实验表明,本发明仿真方法的精度明显高于传统仿真方法。
附图说明
图1为Navier滑移模型示意图;
图2为本发明所涉及方法的流程图。
具体实施方式
本发明所述方法由安装在计算机上的软件程序实现。所述计算机上安装由美国TheMathWorks公司出品的商业数学软件MATLAB和由美国ANSYS公司开发Mechanical APDLProduct Lanucher软件。
本发明所述方法的的流程图如图2所示,具体包括以下步骤:
步骤1,根据静压轴承的工况选择Navier速度滑移模型,模拟轴承-转子系统中流固界面上出现的速度滑移现象。
步骤2,建立考虑速度滑移情况下静压轴承油膜压力分布的有限差分模型,根据油膜结构划分合理网格。
步骤3,用matlab语言编写程序将计算所得的压力分布数据导出,将数据文件存储成txt格式,便于将此计算数据作为压力边界条件加载到转子和轴瓦及轴承座有限元模型中。
步骤4,借助ansys有限元软件建立转子和轴瓦及轴承座的有限元模型,控制有限元模型的网格划分形式与油膜的网格划分一致,以保证数据交换过程的顺利完成。
步骤5,应用APDL参数化有限元分析技术编写相应程序实现matlab与ansys之间的数据交换,编写APDL程序将matlab数据读入ansys,将油膜压力分布的matlab计算结果作为边界条件加载到转子和轴瓦及轴承座有限元模型中,实现对转子和轴瓦及轴承座有限元模型的求解。
步骤6,应用APDL参数化有限元分析技术提取转子和轴瓦及轴承座有限元模型相关节点的变形值。有限元模型在某节点的变形可以分解为X,Y,Z三个方向,节点在X,Y方向上的变形将影响到油膜厚度的变化。
步骤7,将ansys导出的节点变形数据导入到matlab油膜厚度的计算中,实现对油膜压力分布的重新求解。
步骤8,重复步骤3~7,直至转子和轴瓦及轴承座的变形达到稳定,实现了真正意义上静压轴承-转子系统的流体-结构之间的双向耦合。
表1给出了三种不同条件下通过仿真计算得到的一定结构的静压轴承承载力,并列出通过一定的实验方法得到的静压轴承承载力数据。三种条件分别为:无滑移并不考虑流固界面的相互耦合;考虑速度滑移但不进行流固界面的相互耦合分析;考虑速度滑移并进行流固界面的相互耦合分析。表1同时还给出了实验方法得出的静压轴承的承载力,即采用测微仪检测主轴偏心量,并通过砝码对轴承施加静力外载荷,当主轴的偏心率为0.1时静压轴承的承载力。传统的仿真方法一般忽略静压轴承流固界面的速度滑移,同时也不会进行流固耦合分析,本发明仿真方法将两者都作为影响静压轴承性能的因素并考虑在仿真过程中。从表1中可以看出实验测得的静压轴承的承载力为1010N,通过传统仿真方法在不忽略速度滑移和不进行流固耦合分析的情况下轴承的承载力为1120N,考虑速度滑移但不对流固界面的相互耦合作用进行分析时,静压轴承的承载力为1400N,而采用本发明仿真方法得到的静压轴承的承载力为980N。通过计算可得三种不同条件下得到的承载力相对实验结果的误差分别为10.89%、38.6%和2.97%,因此,本发明仿真方法的准确度明显高于传统方法。
表1 不同条件下静压轴承的承载力

Claims (3)

1.一种速度滑移工况下静压主轴系统流固耦合仿真分析方法,其特征在于包括以下步骤:
步骤1,根据静压轴承的工况选择Navier速度滑移模型,模拟轴承-转子系统中流固界面上出现的速度滑移现象;
Navier模型假设流体的滑移速度vs与流体在表面的剪切力成正比,即:
v s = b ∂ u 1 ∂ z | wall - - - ( 1 )
式中,u1是流体沿固体表面的流动速度,Z是沿界面法向方向坐标,wall表示流体和固体相互作用的界面位置;b称之为滑移系数,是虚构固体表面与实际界面的距离;
步骤2,建立考虑速度滑移情况下静压轴承油膜压力分布的有限差分模型,根据油膜结构划分合理网格;
考虑速度滑移条件下的油膜压力分布符合修正后的N-S方程为:
∂ ∂ x { [ - h 3 12 η - hb bx 2 η ( h - h * ) - hh * b ax 2 η ] ∂ p ∂ x + hu 2 } + ∂ ∂ z { ( - h 3 12 η - h 2 b az 2 η ) ∂ p ∂ z } = 0 - - - ( 2 )
式中,u为旋转主轴表面的线速度;bax、baz分别为轴瓦表面周向、轴向滑移系数;bbx为轴颈表面周向滑移系数;h*为油膜周向方向剪切力为零时的油膜厚度;p为油膜压力;η为润滑油的动力粘度;
采用有限差分方法将方程(2)线性化处理,得到如下线性化方程:
P i , j = 1 2 A i , j + 2 C i , j [ ( A i , j + B i , j ) P i + 1 , j + ( A i , j - B i , j ) P i - 1 , j + C i , j P i , j + 1 + C i , j P i , j - 1 - D i , j P i , j ] - - - ( 3 )
其中,
C i , j = β 2 H i , j 3 ( 1 + 6 b 0 ) Δλ 2
β = 2 R L
式中,b0为滑移系数,R为轴承半径,L为轴承长度,Δλ为轴承轴向方向网格长度,为轴承周向方向网格长度,(i,j)为油膜位置坐标,Hi,j表示(i,j)处无量纲油膜厚度;应用超松弛迭代法将线性化方程(3)在matlab软件平台上进行数值求解,即得考虑速度滑移影响下油膜压力分布;
步骤3,用matlab语言编写程序将计算所得的压力分布数据导出,将数据文件存储成txt格式,便于将此计算数据在步骤5中作为压力边界条件加载到转子和轴瓦及轴承座有限元模型中;
步骤4,借助ansys有限元软件建立转子和轴瓦及轴承座的有限元模型,控制有限元模型的网格划分形式与油膜的网格划分一致,以保证数据交换过程的顺利完成;
步骤5,应用APDL参数化有限元分析技术编写相应程序实现matlab与ansys之间的数据交换,编写APDL程序将matlab数据读入ansys,将油膜压力分布的matlab计算结果作为边界条件加载到转子和轴瓦及轴承座有限元模型中,实现对转子和轴瓦及轴承座有限元模型的求解;
步骤6,应用APDL参数化有限元分析技术提取转子和轴瓦及轴承座有限元模型相关节点的变形值;有限元模型在某节点的变形可以分解为X、Y、Z三个方向,节点在X、Y方向上的变形将影响到油膜厚度的变化;
步骤7,将ansys导出的节点变形数据导入到matlab油膜厚度的计算中,实现对油膜压力分布的重新求解;
步骤8,重复步骤3~7,直至转子和轴瓦及轴承座的变形达到稳定,实现真正意义上静压轴承-转子系统的流体-结构之间的双向耦合。
2.根据权利要求1所述的一种速度滑移工况下静压主轴系统流固耦合仿真分析方法,其特征在于,在步骤4所述转子模型建立过程中,采用自下而上的建模方法,先创建关键点,依次构成线和面,用lsize控制转子有限元模型的网格划分,将面绕轴线旋转出周向规则的网格,然后在转子模型的表面建立表面效应单元,用于加载与表面相互垂直的压力载荷,定义转子有限元模型的材料属性,转子的有限元模型建立完成。
3.根据权利要求1所述的一种速度滑移工况下静压主轴系统流固耦合仿真分析方法,其特征在于,在步骤4所述轴瓦及轴承座模型的建立中,静压腔的深度与轴承的尺寸相比很小,所以采用有限元方法对轴瓦建模时忽略轴瓦表面的静压腔,有限元模型的建立过程与转子模型的建立过程相同,但需根据实际情况分别对轴瓦及轴承座定义有限元模型的材料属性。
CN201410328613.8A 2014-07-11 2014-07-11 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法 Pending CN104091022A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410328613.8A CN104091022A (zh) 2014-07-11 2014-07-11 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410328613.8A CN104091022A (zh) 2014-07-11 2014-07-11 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法

Publications (1)

Publication Number Publication Date
CN104091022A true CN104091022A (zh) 2014-10-08

Family

ID=51638738

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410328613.8A Pending CN104091022A (zh) 2014-07-11 2014-07-11 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法

Country Status (1)

Country Link
CN (1) CN104091022A (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095583A (zh) * 2015-07-26 2015-11-25 北京工业大学 一种微尺度下的静压主轴模态分析方法
CN105972081A (zh) * 2016-06-15 2016-09-28 北京工业大学 一种微尺度下空气静压轴承性能优化方法
CN106383930A (zh) * 2016-08-31 2017-02-08 中国石油大学(华东) 一种尾轴承‑转子系统多重流固耦合计算方法
CN107665286A (zh) * 2017-10-31 2018-02-06 华晨汽车集团控股有限公司 一种关于汽车发动机轴承动力学分析方法
CN109117522A (zh) * 2018-07-24 2019-01-01 北京工业大学 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法
CN110287576A (zh) * 2019-06-20 2019-09-27 哈尔滨理工大学 一种基于Matlab的液体静压推力轴承润滑油膜三维温度场显示方法
CN112069710A (zh) * 2020-09-07 2020-12-11 电子科技大学 用于气体静压主轴自激振动的预测方法
CN113027921A (zh) * 2021-02-09 2021-06-25 太原重工股份有限公司 获取静动压轴承油膜压力分布的方法和装置
CN113609734A (zh) * 2021-08-03 2021-11-05 国能大渡河检修安装有限公司 一种滑动轴承水轮机仿真油膜厚度可视化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102880765A (zh) * 2012-10-18 2013-01-16 哈尔滨工业大学 用于气体静压主轴动力学预测的仿真方法
CN102880766A (zh) * 2012-10-18 2013-01-16 哈尔滨工业大学 基于参数化设计的液体静压主轴制作方法
US8543365B1 (en) * 2004-10-25 2013-09-24 Novatorque, Inc. Computer-readable medium, a method and an apparatus for designing and simulating electrodynamic machines implementing conical and cylindrical magnets

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8543365B1 (en) * 2004-10-25 2013-09-24 Novatorque, Inc. Computer-readable medium, a method and an apparatus for designing and simulating electrodynamic machines implementing conical and cylindrical magnets
CN102880765A (zh) * 2012-10-18 2013-01-16 哈尔滨工业大学 用于气体静压主轴动力学预测的仿真方法
CN102880766A (zh) * 2012-10-18 2013-01-16 哈尔滨工业大学 基于参数化设计的液体静压主轴制作方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
LI W L, CHU H M, CHEN M D.: "The partially wetted bearing—extended Reynolds equation", 《TRIBOLOGY INTERNATIONAL》 *
何芝仙,桂长林: "计入轴瓦弹性变形的滑动轴承润滑分析的快速近似算法", 《润滑与密封》 *
刘赵淼,王国斌,申峰: "基于Navier 滑移的油膜缝隙微流动特性数值分析", 《机械工程学报》 *
姚方,马希直: "辊压机液体静压轴承的设计及性能研究", 《润滑与密封》 *
边艳华,陈东菊,范晋伟: "微尺度下液体静压轴承性能研究", 《润滑与密封》 *
陈东菊,边艳华,范晋伟,周帅: "基于速度滑移的液体静压轴承油膜微流动模型", 《中国科技论文在线》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105095583B (zh) * 2015-07-26 2017-12-15 北京工业大学 一种微尺度下的静压主轴模态分析方法
CN105095583A (zh) * 2015-07-26 2015-11-25 北京工业大学 一种微尺度下的静压主轴模态分析方法
CN105972081A (zh) * 2016-06-15 2016-09-28 北京工业大学 一种微尺度下空气静压轴承性能优化方法
CN106383930B (zh) * 2016-08-31 2019-06-28 中国石油大学(华东) 一种尾轴承-转子系统多重流固耦合计算方法
CN106383930A (zh) * 2016-08-31 2017-02-08 中国石油大学(华东) 一种尾轴承‑转子系统多重流固耦合计算方法
CN107665286A (zh) * 2017-10-31 2018-02-06 华晨汽车集团控股有限公司 一种关于汽车发动机轴承动力学分析方法
CN109117522A (zh) * 2018-07-24 2019-01-01 北京工业大学 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法
CN109117522B (zh) * 2018-07-24 2022-11-22 北京工业大学 一种基于matlab-ansys软件的静压溜板流固耦合工作性能的计算方法
CN110287576A (zh) * 2019-06-20 2019-09-27 哈尔滨理工大学 一种基于Matlab的液体静压推力轴承润滑油膜三维温度场显示方法
CN110287576B (zh) * 2019-06-20 2023-04-07 哈尔滨理工大学 一种基于Matlab的液体静压推力轴承润滑油膜三维温度场显示方法
CN112069710A (zh) * 2020-09-07 2020-12-11 电子科技大学 用于气体静压主轴自激振动的预测方法
CN113027921A (zh) * 2021-02-09 2021-06-25 太原重工股份有限公司 获取静动压轴承油膜压力分布的方法和装置
CN113027921B (zh) * 2021-02-09 2022-10-11 太原重工股份有限公司 获取静动压轴承油膜压力分布的方法和装置
CN113609734A (zh) * 2021-08-03 2021-11-05 国能大渡河检修安装有限公司 一种滑动轴承水轮机仿真油膜厚度可视化方法

Similar Documents

Publication Publication Date Title
CN104091022A (zh) 一种速度滑移工况下静压主轴系统流固耦合仿真分析方法
Gao et al. A FEM based modeling method for analyzing the static performance of aerostatic thrust bearings considering the fluid-structure interaction
Tiwari et al. A ghost fluid Lattice Boltzmann method for complex geometries
Zhang et al. A Kriging model for dynamics of mechanical systems with revolute joint clearances
Frosina et al. A modeling approach to study the fluid-dynamic forces acting on the spool of a flow control valve
Song et al. Development and validation of a three-dimensional computational fluid dynamics analysis for journal bearings considering cavitation and conjugate heat transfer
Huang et al. Theoretical and experimental investigation on tip forces and temperature distributions of the brush seal coupled aerodynamic force
Lopez et al. A recommended correction to the kT− kL− ω transition-sensitive eddy-viscosity model
Zhang et al. Determination method of scaling laws based on least square method and applied to rectangular thin plates and rotor-bearing systems
CN104239625A (zh) 一种基于修正流体运动方程线性迭代的稳态求解方法
Ding et al. Exact and efficient isogeometric reanalysis of accurate shape and boundary modifications
Meyer et al. The path-independent M integral implies the creep closure of englacial and subglacial channels
Feng et al. An efficient method to predict the chatter stability of titanium alloy thin-walled workpieces during high-speed milling by considering varying dynamic parameters
Filiz et al. A model for bending, torsional, and axial vibrations of micro-and macro-drills including actual drill geometry—part I: model development and numerical solution
Jhon Ramos et al. Stability analysis of rotating systems supported by textured journal bearings
Venkatesan et al. Finite element computations of viscoelastic two‐phase flows using local projection stabilization
Akyildiz et al. Exact solution of the startup electroosmotic flow of generalized Maxwell fluids in triangular microducts
Esipov et al. Direct numerical simulation of the Segre–Silberberg effect using immersed boundary method
DENG Computational Analysis of Turbulence and Thermal Characteristics in Fluid Film Thrust Bearings
Deng et al. Efficient numerical simulation of injection mold filling with the lattice Boltzmann method
Pan et al. Estimation of dynamic behaviors of hydraulic forging press machine in slow-motion manufacturing process
Wang et al. Parametric study of fluid–solid interaction for single-particle dissipative particle dynamics model
Helgason et al. Optimization Using Arbitrary Lagrangian–Eulerian Formulation of the Navier–Stokes Equations
Lei et al. Numerical calculation of oil film for ship stern bearing based on matrix method
Felici et al. Analysis of vortex induced vibration of a thermowell by high fidelity FSI numerical analysis based on RBF structural modes embedding

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20141008

RJ01 Rejection of invention patent application after publication