一种基于虚功理论的膜式空气弹簧建模方法技术

本发明专利技术提供一种基于虚功理论的膜式空气弹簧建模方法，涉及车辆悬架设计分析领域，包括以下步骤：S1：定义虚位移并得到虚功表达式；S2：求解空气弹簧空气体积相对于气压的变化率；S3：建立气固耦合模型；S4：定义膜式空气弹簧各部分间接触类型与摩擦类型；S5：定义膜式空气弹簧材料类型。本发明专利技术以建立一个精确的空气弹簧模型，能够用于模拟和分析膜式空气弹簧在不同配置和参数下的特性为目的，基于虚功理论同时考虑体积变化、气固耦合以及气体和固体单元边界定义来建立膜式空气弹簧数值模型。单元边界定义来建立膜式空气弹簧数值模型。单元边界定义来建立膜式空气弹簧数值模型。

【技术实现步骤摘要】
一种基于虚功理论的膜式空气弹簧建模方法


[0001]本专利技术涉及车辆悬架控制
，尤其涉及一种膜式空气弹簧的建模方法。

技术介绍

[0002]随着人们对车辆性能要求的不断提高，汽车悬架系统得到了快速的发展，因空气悬架可提高车辆的乘坐舒适性、道路保持能力和控制稳定性，成为了最为流行的悬架类型之一。空气弹簧作为空气悬架系统的主要部件，也成为了学术界和工业界研究的热点。
[0003]通常，空气弹簧主要分为两类，一种称为膜式空气弹簧，其中摩擦壁在压缩过程中在活塞上滚动，另一种是波纹管型空气弹簧，其橡胶壁在压缩冲程中膨胀和弯曲。膜式空气弹簧由于其突出的弹性特性而广泛应用于重型车辆的悬架系统。由于弹性力的非线性和气囊内气体的不确定性，在膜式空气弹簧的早期开发阶段存在特殊的困难，为了解决这一问题，首先应该建立一个精确的空气弹簧模型，以便该模型能够用于模拟和分析膜式空气弹簧在不同配置和参数下的特性。然而现有方法都不能提供有效的方法来解决膜式空气弹簧的建模问题。
[0004]综上所述，急需专利技术一种基于虚功理论的膜式空气弹簧建模方法。

技术实现思路

[0005]本专利技术提供一种基于虚功理论的膜式空气弹簧建模方法，旨在用有限元方法研究非传统材料对膜式空气弹簧垂直特性和物理尺寸的影响，建立能够用于模拟和分析膜式空气弹簧在不同配置和参数下特性的模型。
[0006]为达到以上目的，本专利技术采取的技术方案是：
[0007]一种基于虚功理论的膜式空气弹簧建模方法，所述设计方法包括以下步骤：
[0008]S1、分析膜式空气弹簧压缩与拉伸过程，定义虚位移并得到虚功表达式；
[0009]S2、求解空气弹簧空气体积相对于气压的变化率；
[0010]S3、同时考虑空气单元和气囊外壳的相互作用，建立气固耦合模型；
[0011]S4、定义膜式空气弹簧各部分间接触类型与摩擦类型；
[0012]S5、定义膜式空气弹簧材料类型，完成膜式空气弹簧建模。
[0013]优选地，所述S1中空气弹簧气囊内瞬时空气体积可定义为：
[0014][0015]其中p是空气弹簧的内部空气压力，θ是空气温度。假设粒子的数量为m，则虚功可以表示为：
[0016][0017]空气弹簧的虚功可以进一步定义为：
[0018][0019]其中δW是充气前的虚功，δV是充气时的虚体积变化，δp是虚压强变化，V是气囊的实际体积，上式的微分为：
[0020][0021]其中pdδV是压力负载刚度，是空气体积相对于压力的变化率。由于空气弹簧内的每个单元都被认为处于相同的压力下，根据实际做功原理，虚拟功也可以用每个单元e的总虚拟功表示：
[0022][0023]优选地，所述S2中空气弹簧内部的空气密度可定义为：
[0024][0025]式中，θ
R
和p
R
分别为参考密度ρ
R
下的温度和压力，θ
A
是按热力学温度换算的环境温度，p
A
是环境压力。气囊内的空气体积可以计算为：
[0026][0027]其中m
e
是单元e的质量，m
t
是所有单元的总质量，然后空气体积相对于压力的变化率可以表示为：
[0028][0029]优选地，所述S3中空气单元的节点坐标可以定义为：
[0030][0031]其中N
i
是空气单元的形状函数，可以用等参坐标ξ和η表示，其中X
i
(x,y)为空气单元的参考点坐标，然后空气单元的雅可比行列式可以表示为：
[0032][0033]通过将空气单元的法向量Q与无穷小的面积dA相乘，可以得到以下方程：
[0034][0035]无穷小的体积可以表示为：
[0036][0037]其中X
R
是空气单元的参考坐标。积分之后，空气单元的体积可以定义为：
[0038][0039]设为相对节点坐标，则可以获得无穷小的体积为：
[0040][0041]对于本专利技术中的F3D4空气单元，空气单元体积的表达式可表述为：
[0042][0043]其中是F3D4空气单元的相对节点坐标，其边界无穷小的体积为:
[0044][0045]优选地，所述S4中气囊的上内壁与上板之间的接触以及气囊的下内壁与活塞之间的接触都定义为壳体单元与刚性表面之间的接触。设置气囊为主动体，设置上板和活塞为被动体，因此仅考虑接触表面上的法向压力，将接触表面视为光滑表面，并将数值模型视为无摩擦接触模型。
[0046]数值模型可分为四个部分：气囊模型、空气模型、上板模型和活塞模型。采用四节点壳单元(ABAQUS中的S4R)来定义气囊模型；关于空气模型，在假设空气弹簧内部的空气是理想气体之后，可以根据INP文件的编程用流体单元来定义。在空气模拟过程中，使用了两种类型的空气单元，它们是三维三节点(ABAQUS中的F3D3)和三维四节点(ABAQUS中的F3D4)；上板和活塞均由金属材料制成，这些材料被视为刚体，可以用三维三刚度单元(ABAQUS中的R3D3)定义。
[0047]本专利技术的有益效果在于：
[0048]1、本专利技术通过同时考虑体积变化、气固耦合以及气体和固体单元边界定义来建立膜式空气弹簧数值模型。
[0049]2、本专利技术对膜式空气弹簧的垂直特征进行了分析，该模型能够精确反映空气弹簧特性，能够用于模拟和分析膜式空气弹簧在不同配置和参数下的特性。
附图说明
[0050]为了更清楚的说明本专利技术的实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0051]图1为本专利技术设计流程图。
[0052]图2为本专利技术空气弹簧Rebar单元。
[0053]图3为本专利技术空气弹簧运动过程分析。
[0054]图4为本专利技术完成接触条件与材料类型定义的有限元模型。
[0055]图5为本专利技术膜式空气弹簧在0.3Mpa初始空气压力下，有限元模型(Finite Element Model)、分析模型(Analytical Model)以及实验结果(Experimental Result)的力
‑
位移曲线。
[0056]图6为本专利技术膜式空气弹簧在0.5Mpa初始空气压力下，有限元模型、分析模型以及实验结果的力
‑
位移曲线。
[0057]图7为本专利技术膜式空气弹簧在0.8Mpa初始空气压力下，有限元模型、分析模型以及实验结果的力
‑
位移曲线。
具体实施方式
[0058]需要说明的是，在不冲突的情况下，本专利技术中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于虚功理论的膜式空气弹簧建模方法，其特征在于，所述设计方法包括以下步骤：S1、分析膜式空气弹簧压缩与拉伸过程，定义虚位移并得到虚功表达式；S2、在S1的基础上，求解空气弹簧空气体积相对于气压的变化率；S3、在S2基础上，同时考虑空气单元和气囊外壳的相互作用，建立气固耦合模型；S4、在S3的基础上，定义膜式空气弹簧各部分间接触类型与摩擦类型；S5、在S4基础上，定义膜式空气弹簧材料类型，完成膜式空气弹簧建模。2.根据权利要求1所述的基于虚功理论的膜式空气弹簧建模方法，其特征在于：所述S1中，空气弹簧气囊内瞬时空气体积可定义为：其中p是空气弹簧的内部空气压力，θ是空气温度，假设粒子的数量为m，则虚功可以表示为：其中F
n
表示施加在第n个粒子上的力，d
n
表示沿第n个粒子最小作用，空气弹簧的虚功可以进一步定义为：其中δW是充气前的虚功，δV是充气时的虚体积变化，δp是虚压强变化，V是气囊的实际体积，上式的微分为：其中pdδV是压力负载刚度，是空气体积相对于压力的变化率，由于空气弹簧内的每个单元都被认为处于相同的压力下，根据实际做功原理，虚拟功也可以用每个单元e的总虚拟功表示：δW
*
＝∑[δW
e
‑
pδV
e
‑
δp(V
e
‑
V
e
)]。3.根据权利要求2所述的基于虚功理论的膜式空气弹簧建模方法，所述S2中，气弹簧内部的空气密度可定义为：式中，θ
R
和p
R
分别为参考密度ρ
R
下的温度和压力，θ
A
是按热力学温度换算的环境温度，p
A
是环境压力...

【专利技术属性】
技术研发人员：赵晶，安润兴，刘泰佑，刘锦灿，李国全，庞滔，李锦豪，
申请(专利权)人：广东溢康通空气弹簧有限公司东北大学，
类型：发明
国别省市：