CN106599522B - 非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 - Google Patents
非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 Download PDFInfo
- Publication number
- CN106599522B CN106599522B CN201710002608.1A CN201710002608A CN106599522B CN 106599522 B CN106599522 B CN 106599522B CN 201710002608 A CN201710002608 A CN 201710002608A CN 106599522 B CN106599522 B CN 106599522B
- Authority
- CN
- China
- Prior art keywords
- spring
- simulation calculation
- order gradient
- main
- main spring
- 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
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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F1/00—Springs
- F16F1/02—Springs made of steel or other material having low internal friction; Wound, torsion, leaf, cup, ring or the like springs, the material of the spring not being relevant
- F16F1/18—Leaf springs
- F16F1/185—Leaf springs characterised by shape or design of individual leaves
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F1/00—Springs
- F16F1/02—Springs made of steel or other material having low internal friction; Wound, torsion, leaf, cup, ring or the like springs, the material of the spring not being relevant
- F16F1/18—Leaf springs
- F16F1/26—Attachments or mountings
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
- F16F—SPRINGS; SHOCK-ABSORBERS; MEANS FOR DAMPING VIBRATION
- F16F3/00—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic
- F16F3/02—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic with springs made of steel or of other material having low internal friction
- F16F3/023—Spring units consisting of several springs, e.g. for obtaining a desired spring characteristic with springs made of steel or of other material having low internal friction composed only of leaf springs
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Mechanical Engineering (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Vehicle Body Suspensions (AREA)
Abstract
本发明涉及非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,属于悬架钢板弹簧技术领域。本发明可根据各片主簧和副簧的结构参数、弹性模量、主簧和副簧的初始切线弧高设计值,在接触载荷仿真计算的基础上,对非等偏频一级渐变刚度板簧悬架的刚度特性进行仿真计算。通过实例仿真计算和样机刚度特性试验可知,表明本发明所提供的非等偏频一级渐变刚度板簧悬架的挠度特性的仿真计算法是正确的,确保刚度特性满足设计要求,为非等偏频一级渐变刚度板簧悬架设计和特性仿真验证及CAD软件开发奠定了的技术基础。利用该方法可提产品的设计水平、质量和性能,提高车辆行驶平顺性和安全性;同时,还可以降低设计和试验费用,加快产品开发速度。
Description
技术领域
本发明涉及车辆悬架板簧,特别是非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法。
背景技术
为了满足一级渐变刚度板簧的主簧强度的要求,通常使副簧尽早起作用承担载荷而降低主簧应力,即采用非等偏频一级渐变刚度板簧悬架,其中,悬架在不同载荷下的刚度特性不仅影响悬架偏频和车辆行驶平顺性,而且还影响悬架静挠度、动挠度、最大限位挠度及车辆行驶安全性;对于给定设计结构板簧悬架的刚度特性的仿真计算,可确保其刚度特性满足悬架系统的设计要求,而且为主簧和副簧的初始切线弧高及主副簧间隙的仿真验证奠定可靠的计算基础。然而,由于受渐变接触过程中的挠度和渐变刚度计算及接触载荷仿真计算的制约,先前一直未能给出非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,不能满足车辆行业快速发展及现代化CAD软件开发的要求。随着车辆行驶速度及其对平顺性要求的不断提高,对非等偏频一级渐变刚度板簧悬架提出了更高要求,因此,必须建立一种精确、可靠的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,为非等偏频一级渐变刚度板簧设计、特性仿真计算及CAD软件开发奠定可靠的技术基础,满足车辆行业快速发展、车辆行驶平顺性及对非等偏频一级渐变刚度板簧设计的要求,进一步提高非等偏频一级渐变刚度板簧的设计水平、产品质量和可靠性,及车辆行驶平顺性和安全性;同时,降低产品设计及试验费用,加快产品开发速度。
发明内容
针对上述现有技术中存在的缺陷,本发明所要解决的技术问题是提供一种简便、可靠的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,其仿真计算流程图,如图1所示。非等偏频一级渐变刚度板簧悬架的一半对称结构如图2所示,是由主簧1和副簧2所组成的,一级渐变刚度板簧的一半总跨度,即为首片主簧的一半作用长度为L1t,骑马螺栓夹紧距的一半为L0,板簧的宽度为b,弹性模量为E。主簧1的片数为n,各片主簧的厚度为hi,一半作用长度为Lit,一半夹紧长度Li=Lit-L0/2,i=1,2,…n。副簧2的片数为m,各片副簧的厚度为hAj,一半作用长度为LAjt,一半夹紧长度LAj=LAjt-L0/2,j=1,2,…m。通过主簧和副簧初始切线弧高,确保副簧首片端部上表面与主簧末片端部下表面之间设置有一定的主副簧间隙δMA,以满足渐变刚度板簧开始接触载荷和完全接触载荷、主簧应力强度和悬架渐变刚度的设计要求,并且还应该满足板簧安装及在额定载荷下剩余切线弧高的设计要求。非等偏频一级渐变刚度板簧的空载载荷P0,开始接触载荷为Pk,完全接触载荷为Pw;为了满足主簧应力强度的要求,悬架开始接触载荷偏频f0k与完全接触载荷偏频f0w不相等,即板簧悬架设计为非等偏频型。根据各片主簧和副簧的结构参数、弹性模量、主簧和副簧的初始切线弧高,在接触载荷仿真计算的基础上,对非等偏频一级渐变刚度板簧悬架的刚度特性进行仿真计算。
为解决上述技术问题,本发明所提供的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,其特征在于采用以下仿真计算步骤:
(1)非等偏频一级渐变刚度板簧悬架的主簧夹紧刚度及主副簧复合夹紧刚度的仿真计算:
I步骤:各不同片数重叠段的等效厚度计算
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;副簧片数m,各片副簧的厚度hAj,j=1,2,…,m;主副簧的总片数N=n+m,对非等偏频一级渐变刚度板簧各不同片数k重叠段的等效厚度hke,k=1,2,…,N,分别进行计算,即
其中,主簧根部重叠部分的等效厚度主副簧根部重叠部分的等效厚度
II步骤:主簧夹紧刚度KM的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;主簧片数n,各片主簧的一半夹紧长度Li,i=1,2,…,n,及I步骤中仿真计算得到的hke,k=i=1,2,…,n,对主簧夹紧刚度KM进行计算,即
III步骤:主副簧复合夹紧刚度KMA的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;主簧片数n,各片主簧的一半夹紧长度Li,i=1,2,…,n;副簧片数m,各片副簧的一半夹紧长度分别为LAj=Ln+j,j=1,2,…,m;主副簧的总片数N=n+m,及I步骤中仿真计算得到的hke,k=1,2,…,N,对主副簧的复合夹紧刚度KMA进行仿真计算,即
(2)非等偏频一级渐变刚度板簧悬架的开始接触载荷Pk的仿真计算:
A步骤:末片主簧下表面初始曲率半径RM0b的确定
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;首片主簧的一半夹紧长度L1,主簧初始切线弧高HgM0,确定末片主簧下表面初始曲率半径RM0b,即
B步骤:首片副簧上表面的初始曲率半径RA0a的确定
根据首片副簧的一半夹紧长度LA1,副簧初始切线弧高HgA0,确定首片副簧上表面的初始曲率半径RA0a,即
C步骤:开始接触载荷Pk的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;首片主簧的一半夹紧长度L1,步骤(1)的I步骤中计算得到的hMe;步骤(2)的A步骤中所确定的RM0b,B步骤中所确定的RA0a,对开始接触载荷Pk进行仿真计算,即
(3)非等偏频一级渐变刚度板簧悬架的完全接触载荷Pw的仿真计算:
a步骤:主副簧完全接触时的主簧切线弧高表达式HgMw的建立:
根据主簧初始切线弧高HgM0,步骤(1)中仿真计算得到的KM和KMA,步骤(2)中仿真计算得到的Pk,以完全接触载荷Pw为参变量,建立完全接触时的主簧切线弧高表达式HgMw,即
式中,A、B和C为所定义的渐变挠度计算的中间参数,B=-CPk,
b步骤:主副簧完全接触时的末片主簧下表面曲率半径RMwb表达式的建立:
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;首片主簧的一半夹紧长度L1,a步骤中所建立的HgMw,以完全接触载荷Pw为参变量,建立完全接触时的末片主簧下表面曲率半径表达式RMwb,即
c步骤:完全接触载荷Pw的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;首片主簧的一半夹紧长度L1,步骤(1)的I步骤中计算得到的hMe;步骤(2)中所确定的RM0b和Pk,及b步骤中所建立的RMwb,以完全接触载荷Pw为参变量,建立完全接触载荷Pw的仿真计算数学模型,即
求解上述数学模型,便可得到非等偏频一级渐变板簧的完全接触载荷Pw;
(4)非等偏频一级渐变刚度板簧悬架不同载荷下的夹紧刚度特性的仿真计算:
根据步骤(1)中仿真计算得到的KM和KMA,步骤(2)中仿真计算得到的Pk,步骤(3)中仿真计算得到的Pw,对非等偏频一级渐变刚度板簧悬架在不同载荷P下的夹紧刚度特性进行计算,即
本发明比现有技术具有的优点
由于受渐变接触过程中的挠度和渐变刚度计算及接触载荷仿真计算的制约,先前一直未能给出非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,不能满足车辆行业快速发展及现代化CAD软件开发的要求。悬架在不同载荷下的刚度特性不仅影响悬架偏频和车辆行驶平顺性,而且还影响悬架静挠度、动挠度、最大限位挠度及车辆行驶安全性;对于给定设计结构板簧悬架的刚度特性的仿真计算,可确保其刚度特性满足悬架系统的设计要求。本发明可根据各片主簧和副簧的结构参数、弹性模量、主簧和副簧的初始切线弧高、额定载荷,在接触载荷仿真计算的基础上,对非等偏频一级渐变刚度板簧悬架的刚度特性进行仿真计算。通过实例仿真验证和样机加载挠度及刚度特性试验测试可知,本发明所提供的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法是正确的,为非等偏频一级渐变刚度板簧的主副簧初始切线弧高和主副簧间隙仿真验证提供了可靠的技术基础。利用该方法可可提产品的设计水平、质量和性能,进一步提高车辆的行驶平顺性和安全性;同时,还可以降低设计和试验测试费用,加快产品开发速度。
附图说明
为了更好地理解本发明,下面结合附图做进一步的说明。
图1是非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算流程图;
图2是非等偏频一级渐变刚度板簧的一半对称结构示意图;
图3是实施例的仿真得到的主簧ANSYS变形仿真云图;
图4是实施例的仿真得到的主副簧ANSYS变形仿真云图;
图5是实施例的仿真计算得到的非等偏频一级渐变刚度板簧悬架的夹紧刚度K随载荷P变化的特性曲线。
具体实施方案
下面通过实施例对本发明作进一步详细说明。
实施例:某非等偏频一级渐变刚度板簧悬架的宽度b=63mm,弹性模量E=200GPa;跨度的一半即一半作用长度L1t=525mm,骑马螺栓夹紧距的一半L0=50mm。主簧片数n=3片,副簧片数m=2片,主副簧的总片数N=n+m=5。其中,各片主簧的厚度h1=h2=h3=8mm,各片主簧的一半作用长度分别为L1t=525mm,L2t=450mm,L3t=350mm;一半夹紧长度分别为L1=L1t-L0/2=500mm,L2=L2t-L0/2=425mm,L3=L3t-L0/2=325mm。各片副簧的厚度hA1=hA2=13mm,各片副簧的一半作用长度分别为LA1t=250mm,LA2t=150mm,一半夹紧长度分别为LA1=L4=LA1t-L0/2=225mm,LA2=L5=LA2t-L0/2=125mm。主簧初始切线弧高HgM0=102mm,副簧的初始切线弧高HgA0=12mm,额定载荷PN=7227N。根据各片主簧和副簧的结构参数、弹性模量、主簧和副簧的初始切线弧高、额定载荷,对该非等偏频一级渐变刚度板簧悬架的刚度特性进行仿真计算。
本发明实例所提供的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,其仿真计算流程如图1所示,具体仿真计算步骤如下:
(1)非等偏频一级渐变刚度板簧悬架的主簧夹紧刚度及主副簧复合夹紧刚度的仿真计算:
I步骤:各不同片数重叠段的等效厚度计算
根据主簧片数n=3,各片主簧的厚度h1=h2=h3=8mm;副簧片数m=2,各片副簧的厚度hA1=hA2=13mm,主副簧的总片数N=n+m=5,对非等偏频一级渐变刚度板簧各不同片数k重叠段的等效厚度hke,k=1,2,…,N,分别进行计算,即
h1e=h1=8.0mm
其中,主簧根部重叠部分等效厚度hMe=hne=h3e=11.5mm;主副簧根部重叠部分等效厚度hMAe=hNe=h5e=18.1mm;
II步骤:主簧夹紧刚度KM的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b=63mm,弹性模量E=200GPa;主簧片数n=3,各片主簧的一半夹紧长度L1=500mm,L2=425mm,L3=325mm,及I步骤中仿真计算得到的h1e=8.0mm,h2e=10.1mm,h3e=11.5mm,k=i=1,2,...,n,对主簧夹紧刚度KM进行计算,即
根据主簧各片的厚度和一半夹紧长度,建立主簧的一半对称夹紧结构的ANSYS仿真模型,在端部施加一集中力F=1330N,进行ANSYS变形仿真和刚度验证,仿真得到的主簧ANSYS变形仿真云图,如图3所示,其中,端部最大挠度fMmax=34.984mm,因此,主簧夹紧刚度ANSYS仿真验证值KM=2F/fMmax=76.034N/mm,与计算值KM=75.4N/mm的相对偏差仅为0.84%,结果表明主簧夹紧刚度KM的仿真计算值是准确可靠的。
III步骤:主副簧复合夹紧刚度KMA的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b=63mm,弹性模量E=200GPa;主簧片数n=3,各片主簧的一半夹紧长度L1=500mm,L2=425mm,L3=325mm;副簧片数m=2,各片副簧的一半夹紧长度分别为LA1=L4=225mm,LA2=L5=125mm,主副簧的总片数N=5,及I步骤中仿真计算得到的h1e=8.0mm,h2e=10.1mm,h3e=11.5mm,h4e=15.5mm,h5e=18.1mm,k=1,2,...,N,对主副簧的复合夹紧刚度KMA进行仿真计算,即
根据主簧和副簧的各片的厚度和一半夹紧长度,弹性模量E,建立主副簧的一半对称夹紧结构的ANSYS仿真模型,在端部施加一集中力F=4000N,进行ANSYS变形仿真和刚度验证,仿真得到的ANSYS变形仿真云图,如图4所示,其中,端部最大挠度fMA2max=45.44mm,因此,主副簧的复合夹紧刚度ANSYS仿真验证值KMA=2F/fMA2max=176.05N/mm,与仿真计算值KMA1=172.9N/mm的相对偏差仅为1.82%,结果表明主副簧复合夹紧刚度KMA的仿真计算值是准确可靠的。
(2)非等偏频一级渐变刚度板簧悬架的开始接触载荷Pk的仿真计算:
A步骤:末片主簧下表面初始曲率半径RM0b的确定
根据主簧片数n=3,各片主簧的厚度h1=h2=h3=8mm,首片主簧的一半夹紧长度L1=500mm,初始切线弧高HgM0=102mm,确定末片主簧下表面初始曲率半径RM0b,即
B步骤:首片副簧上表面初始曲率半径RA0a的确定
根据首片副簧的一半夹紧长度LA1=225mm,副簧初始切线弧高HgA0=12mm,确定首片副簧上表面初始曲率半径RA0a,即
C步骤:开始接触载荷Pk的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b=63mm,弹性模量E=200GPa;首片主簧的一半夹紧长度L1=500mm,步骤(1)的I步骤中计算得到的hMe=11.5mm;步骤(2)的A步骤中所确定的RM0b=1300.5mm,B步骤中所确定的RA0a=2115.4mm,对开始接触载荷Pk进行仿真计算,即
(3)非等偏频一级渐变刚度板簧悬架的完全接触载荷Pw的仿真计算:
a步骤:主副簧完全接触时的主簧切线弧高表达式HgMw的建立:
根据主簧初始切线弧高HgM0=102mm,步骤(1)中仿真计算得到的KM=75.4N/mm和KMA=172.9N/mm,步骤(2)中仿真计算得到的Pk=1911N,以完全接触载荷Pw为参变量,建立完全接触时的主簧切线弧高表达式HgMw,即
式中,B=-1911C,
b步骤:主副簧完全接触时的末片主簧下表面曲率半径RMwb表达式的建立:
根据主簧片数n=3,各片主簧的厚度hi=8mm,i=1,2,…n;首片主簧的一半夹紧长度L1=500mm,a步骤中所建立的HgMw,以完全接触载荷Pw为参变量,建立完全接触时的末片主簧下表面曲率半径表达式RMwb,即
c步骤:完全接触载荷Pw的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b=63mm,弹性模量E=200GPa;首片主簧的一半夹紧长度L1=500mm,步骤(1)中计算得到的hMe=11.5mm,步骤(2)中所确定的RM0b=1300.5mm和Pk=1911N,及b步骤中所建立的RMwb,以完全接触载荷Pw为参变量,建立完全接触载荷仿真Pw的计算数学模型,即
利用Matlab计算程序,求解上述数学模型,便可得到该非等偏频一级渐变板簧的完全接触载荷的仿真计算值Pw=3843N。
(4)非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算:
根据额定载荷PN=7227N,步骤(1)中仿真计算得到的KM=75.4N/mm和KMA=172.9N/mm,步骤(2)中仿真计算得到的Pk=1911N,步骤(3)中仿真计算得到的Pw=3843N,对该非等偏频一级渐变刚度板簧悬架在不同载荷下的夹紧刚度特性进行仿真计算,即
利用Matlab计算程序,仿真计算得到的该非等偏频一级渐变刚度板簧悬架的夹紧刚度K随载荷P变化的特性曲线,如图5所示,其中,在开始接触载荷Pk=1911N下的夹紧刚度分别为Kk=KM=75.4N/mm,在完全接触载荷Pw=3843N和额定载荷PN=7227N下的夹紧刚度Kw=KN=KMA=172.9N/mm,在区间[Pk,Pw]时,悬架板簧的夹紧刚度K随载荷P而变化。
通过样机加载挠度及刚度特性试验测试可知,表明本发明所提供的非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法是正确的,为非等偏频一级渐变刚度板簧悬架的夹紧刚度特性仿真奠定了可靠的技术方法。利用该方法可得到准确可靠的在不同载荷下的夹紧刚度仿真计算值,可提高产品的设计水平、质量和性能及车辆的行驶平顺性;同时,降低设计及试验费用,加快产品开发速度。
Claims (1)
1.非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法,其中,各片板簧为以中心穿装孔对称的结构,安装夹紧距的一半为骑马螺栓夹紧距的一半;通过主簧和副簧的初始切线弧高及渐变间隙,确保满足悬架偏频特性和主簧应力强度设计要求,即非等偏频一级渐变刚度板簧悬架;根据各片主簧和副簧的结构参数、弹性模量、主簧和副簧的初始切线弧高,在接触载荷仿真计算的基础上,对非等偏频一级渐变刚度板簧悬架的刚度特性进行仿真计算,其特征在于采用以下具体仿真计算步骤:
(1)非等偏频一级渐变刚度板簧悬架的主簧夹紧刚度及主副簧复合夹紧刚度的仿真计算:
I步骤:各不同片数重叠段的等效厚度计算
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;副簧片数m,各片副簧的厚度hAj,j=1,2,…,m;主副簧的总片数N=n+m,对非等偏频一级渐变刚度板簧各不同片数k重叠段的等效厚度hke,k=1,2,…,N,分别进行计算,即
其中,主簧根部重叠部分的等效厚度主副簧根部重叠部分的等效厚度
II步骤:主簧夹紧刚度KM的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;主簧片数n,各片主簧的一半夹紧长度Li,i=1,2,…,n,及I步骤中仿真计算得到的hke,k=i=1,2,…,n,对主簧夹紧刚度KM进行计算,即
III步骤:主副簧复合夹紧刚度KMA的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;主簧片数n,各片主簧的一半夹紧长度Li,i=1,2,…,n;副簧片数m,各片副簧的一半夹紧长度分别为LAj=Ln+j,j=1,2,…,m;主副簧的总片数N=n+m,及I步骤中仿真计算得到的hke,k=1,2,…,N,对主副簧的复合夹紧刚度KMA进行仿真计算,即
(2)非等偏频一级渐变刚度板簧悬架的开始接触载荷Pk的仿真计算:
A步骤:末片主簧下表面初始曲率半径RM0b的确定
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;首片主簧的一半夹紧长度L1,主簧初始切线弧高HgM0,确定末片主簧下表面初始曲率半径RM0b,即
B步骤:首片副簧上表面的初始曲率半径RA0a的确定
根据首片副簧的一半夹紧长度LA1,副簧初始切线弧高HgA0,确定首片副簧上表面的初始曲率半径RA0a,即
C步骤:开始接触载荷Pk的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;首片主簧的一半夹紧长度L1,步骤(1)的I步骤中计算得到的hMe;步骤(2)的A步骤中所确定的RM0b,B步骤中所确定的RA0a,对开始接触载荷Pk进行仿真计算,即
(3)非等偏频一级渐变刚度板簧悬架的完全接触载荷Pw的仿真计算:
a步骤:主副簧完全接触时的主簧切线弧高表达式HgMw的建立:
根据主簧初始切线弧高HgM0,步骤(1)中仿真计算得到的KM和KMA,步骤(2)中仿真计算得到的Pk,以完全接触载荷Pw为参变量,建立完全接触时的主簧切线弧高表达式HgMw,即
式中,A、B和C为所定义的渐变挠度计算的中间参数,B=-CPk,
b步骤:主副簧完全接触时的末片主簧下表面曲率半径RMwb表达式的建立:
根据主簧片数n,各片主簧的厚度hi,i=1,2,…,n;首片主簧的一半夹紧长度L1,a步骤中所建立的HgMw,以完全接触载荷Pw为参变量,建立完全接触时的末片主簧下表面曲率半径表达式RMwb,即
c步骤:完全接触载荷Pw的仿真计算
根据非等偏频一级渐变刚度板簧的宽度b,弹性模量E;首片主簧的一半夹紧长度L1,步骤(1)的I步骤中计算得到的hMe;步骤(2)中所确定的RM0b和Pk,及b步骤中所建立的RMwb,以完全接触载荷Pw为参变量,建立完全接触载荷Pw的仿真计算数学模型,即
求解上述数学模型,便可得到非等偏频一级渐变板簧的完全接触载荷Pw;
(4)非等偏频一级渐变刚度板簧悬架不同载荷下的夹紧刚度特性的仿真计算:
根据步骤(1)中仿真计算得到的KM和KMA,步骤(2)中仿真计算得到的Pk,步骤(3)中仿真计算得到的Pw,对非等偏频一级渐变刚度板簧悬架在不同载荷P下的夹紧刚度特性进行计算,即
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710002608.1A CN106599522B (zh) | 2017-01-03 | 2017-01-03 | 非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710002608.1A CN106599522B (zh) | 2017-01-03 | 2017-01-03 | 非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106599522A CN106599522A (zh) | 2017-04-26 |
CN106599522B true CN106599522B (zh) | 2019-07-26 |
Family
ID=58582236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710002608.1A Expired - Fee Related CN106599522B (zh) | 2017-01-03 | 2017-01-03 | 非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106599522B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102817951A (zh) * | 2012-07-17 | 2012-12-12 | 山东理工大学 | 汽车钢板弹簧厚度的拆分设计方法 |
US8347505B2 (en) * | 2008-10-13 | 2013-01-08 | Baker Hughes Incorporated | Method for fabricating a cylindrical spring by compressive force |
CN105550483A (zh) * | 2016-02-23 | 2016-05-04 | 山东恒日悬架弹簧有限公司 | 端部非等构少片斜线型变截面钢板弹簧的设计方法 |
CN105840702A (zh) * | 2016-06-12 | 2016-08-10 | 山东理工大学 | 非端部接触式少片抛物线型变截面副簧长度的设计方法 |
CN106033484A (zh) * | 2016-05-13 | 2016-10-19 | 山东理工大学 | 端部接触式少片抛物线型变截面副簧长度的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030137087A1 (en) * | 2001-12-19 | 2003-07-24 | Hadi Rod G. | Hydraulic bushing with springs in series |
-
2017
- 2017-01-03 CN CN201710002608.1A patent/CN106599522B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8347505B2 (en) * | 2008-10-13 | 2013-01-08 | Baker Hughes Incorporated | Method for fabricating a cylindrical spring by compressive force |
CN102817951A (zh) * | 2012-07-17 | 2012-12-12 | 山东理工大学 | 汽车钢板弹簧厚度的拆分设计方法 |
CN105550483A (zh) * | 2016-02-23 | 2016-05-04 | 山东恒日悬架弹簧有限公司 | 端部非等构少片斜线型变截面钢板弹簧的设计方法 |
CN106033484A (zh) * | 2016-05-13 | 2016-10-19 | 山东理工大学 | 端部接触式少片抛物线型变截面副簧长度的设计方法 |
CN105840702A (zh) * | 2016-06-12 | 2016-08-10 | 山东理工大学 | 非端部接触式少片抛物线型变截面副簧长度的设计方法 |
Non-Patent Citations (2)
Title |
---|
三段式渐变刚度钢板弹簧的参数化设计与研究匹配研究;陈明媚;《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》;20150815(第08期);第2-4章 |
钢板弹簧重叠部分应力及许用厚度计算理论的研究;宋群等;《山东理工大学学报(自然科学版)》;20160731;第30卷(第4期);第11-14页 |
Also Published As
Publication number | Publication date |
---|---|
CN106599522A (zh) | 2017-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106599522B (zh) | 非等偏频一级渐变刚度板簧悬架的刚度特性的仿真计算法 | |
CN106650168A (zh) | 基于限位挠度的非等偏频一级渐变板簧应力强度的校核法 | |
CN106682337B (zh) | 基于偏频仿真的一级渐变刚度板簧接触载荷的调整设计法 | |
CN106844925A (zh) | 基于偏频仿真的两级副簧式渐变刚度板簧接触载荷的调整设计法 | |
CN106802996A (zh) | 两级副簧式非等偏频型渐变刚度板簧接触载荷的验算方法 | |
CN107013616B (zh) | 高强度一级渐变刚度板簧夹紧刚度特性的仿真计算方法 | |
CN106777793B (zh) | 两级副簧式非等偏频型渐变刚度板簧刚度特性的计算方法 | |
CN106763378B (zh) | 非等偏频一级渐变刚度板簧最大限位挠度的仿真验算法 | |
CN106812846B (zh) | 基于偏频仿真的两级主簧式渐变刚度板簧接触载荷调整设计法 | |
CN106777789B (zh) | 非等偏频型渐变刚度板簧接触载荷的仿真计算方法 | |
CN106777804B (zh) | 基于偏频仿真的三级渐变刚度板簧接触载荷的调整设计法 | |
CN106678223B (zh) | 非等偏频一级渐变刚度板簧的接触载荷的仿真验算法 | |
CN106682357B (zh) | 高强度三级渐变刚度板簧悬架系统偏频特性的仿真计算方法 | |
CN106650166A (zh) | 非等偏频一级渐变刚度板簧悬架的偏频特性的仿真计算法 | |
CN106763386B (zh) | 高强度两级渐变刚度板簧悬架系统偏频特性的仿真计算法 | |
CN106777798B (zh) | 高强度两级渐变刚度板簧的刚度特性的计算方法 | |
CN106545609B (zh) | 两级副簧式非等偏频渐变刚度板簧刚度特性的仿真计算法 | |
CN106802994B (zh) | 两级副簧式非等偏频型渐变刚度板簧根部最大应力的仿真计算法 | |
CN106777800B (zh) | 高强度两级渐变刚度板簧的刚度特性的仿真计算方法 | |
CN106855907B (zh) | 两级副簧式非等偏频型渐变刚度板簧悬架偏频特性的仿真计算方法 | |
CN106704428B (zh) | 基于刚度和应力的一级渐变刚度板簧各片主簧的设计方法 | |
CN106682360B (zh) | 高强度两级渐变刚度主副簧的最大应力特性的仿真计算法 | |
CN106682359B (zh) | 两级副簧式非等偏频型渐变刚度板簧主簧挠度的计算方法 | |
CN107061585B (zh) | 高强度两级渐变刚度板簧的主簧初始切线弧高的设计方法 | |
CN106777794B (zh) | 高强度两级渐变刚度板簧的主簧挠度的计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for 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 |
Granted publication date: 20190726 Termination date: 20220103 |
|
CF01 | Termination of patent right due to non-payment of annual fee |