发明名称 高速轨道车辆一系垂向悬架最优阻尼比的设计方法
摘要 本发明涉及高速轨道车辆一系垂向悬架最优阻尼比的设计方法,属于高速轨道车辆悬置技术领域。本发明通过建立高速轨道车辆1/4车体二自由度行驶垂向振动模型,建立了一系垂向悬架系统阻尼比的目标函数,进而计算得到一系垂向悬架系统的最佳阻尼比和最大许用阻尼比,并利用黄金分割原理,计算得到一系垂向悬架系统的最优阻尼比。通过设计实例及SIMPACK仿真验证可知,该方法可得到准确可靠的一系垂向悬架系统的最优阻尼比值,为高速轨道车辆一系垂向悬架阻尼比的设计提供了可靠的设计方法。利用该方法,不仅可提高高速轨道车辆悬置系统的设计水平及产品质量,提高车辆乘坐舒适性和安全性,还可降低产品设计及试验费用,增强我国轨道车辆的国际市场竞争力。
申请公布号 CN105117554A 申请公布日期 2015.12.02
申请号 CN201510559804.X 申请日期 2015.09.06
申请人 山东理工大学 发明人 周长城;于曰伟;赵雷雷
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 代理人
主权项 高速轨道车辆一系垂向悬架最优阻尼比的设计方法,其具体设计步骤如下:(1)建立高速轨道车辆1/4车体二自由度行驶垂向振动微分方程:根据1/4单节车体的满载质量m,悬架弹簧的等效刚度K<sub>1</sub>;待设计一系垂向悬架的阻尼比ξ,其中,悬架减振器的等效阻尼系数<img file="FDA0000796318740000011.GIF" wi="314" he="90" />悬架减振器的端部连接等效刚度K<sub>01</sub>;以减振器活塞杆的垂向位移z<sub>01</sub>及车体质心的垂向位移z<sub>1</sub>为坐标;以轨道高低不平顺随机输入z<sub>v</sub>为输入激励,建立高速轨道车辆1/4车体二自由度行驶垂向振动微分方程,即:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open = "{" close = ""><mtable><mtr><mtd><mi>m</mi><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mn>01</mn></msub><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>01</mn></msub><mo>)</mo><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mi>v</mi></msub><mo>)</mo><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mn>2</mn><mi>&xi;</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mo>(</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>01</mn></msub><mo>-</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mi>v</mi></msub><mo>)</mo><mo>-</mo><msub><mi>K</mi><mn>01</mn></msub><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>01</mn></msub><mo>)</mo><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000012.GIF" wi="782" he="173" /></maths>(2)确定车体垂向振动位移频率响应函数<img file="FDA0000796318740000018.GIF" wi="263" he="88" />根据步骤(1)中所建立的高速轨道车辆1/4车体二自由度行驶垂向振动微分方程,经傅里叶变换,确定车体垂向振动位移频率响应函数<img file="FDA0000796318740000019.GIF" wi="269" he="88" />即:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>H</mi><msub><mrow><mo>(</mo><mi>j</mi><mi>&omega;</mi><mo>)</mo></mrow><mrow><msub><mi>z</mi><mn>1</mn></msub><mo>~</mo><msub><mi>z</mi><mi>v</mi></msub></mrow></msub><mo>=</mo><mfrac><mrow><msub><mi>K</mi><mn>1</mn></msub><msub><mi>K</mi><mn>1</mn></msub><mo>+</mo><mn>2</mn><mi>&xi;</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mrow><mo>(</mo><msub><mi>K</mi><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mn>01</mn></msub><mo>)</mo></mrow><mi>j</mi><mi>&omega;</mi></mrow><mrow><mo>-</mo><mn>2</mn><mi>&xi;</mi><mi>m</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><msup><mi>j&omega;</mi><mn>3</mn></msup><mo>-</mo><msub><mi>K</mi><mn>01</mn></msub><msup><mi>m&omega;</mi><mn>2</mn></msup><mo>+</mo><mn>2</mn><mi>&xi;</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mrow><mo>(</mo><msub><mi>K</mi><mn>1</mn></msub><mo>+</mo><msub><mi>K</mi><mn>01</mn></msub><mo>)</mo></mrow><mi>j</mi><mi>&omega;</mi><mo>+</mo><msub><mi>K</mi><mn>1</mn></msub><msub><mi>K</mi><mn>01</mn></msub></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000013.GIF" wi="1468" he="171" /></maths>(3)建立一系垂向悬架系统阻尼比目标函数J(ξ)的解析表达式:I步骤:根据轨道高低不平顺大小幅值参数G,车辆行驶速度v,及步骤(2)中所确定的车体垂向振动位移频率响应函数<img file="FDA00007963187400000110.GIF" wi="272" he="98" />建立一系垂向悬架系统阻尼比的目标函数J(ξ),即:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>Gv&pi;</mi><mn>2</mn></msup><msubsup><mi>K</mi><mn>1</mn><mn>2</mn></msubsup><mo>&lsqb;</mo><mn>4</mn><mrow><mo>(</mo><msup><mi>&mu;</mi><mn>3</mn></msup><mo>+</mo><mn>3</mn><msup><mi>&mu;</mi><mn>2</mn></msup><mo>+</mo><mn>3</mn><mi>&mu;</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><msup><mi>&xi;</mi><mn>2</mn></msup><mo>+</mo><msup><mi>&mu;</mi><mn>2</mn></msup><mo>&rsqb;</mo></mrow><mrow><msup><mi>&mu;</mi><mn>2</mn></msup><mi>m</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mi>&xi;</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000014.GIF" wi="951" he="177" /></maths>式中,μ=K<sub>01</sub>/K<sub>1</sub>;II步骤:根据I步骤所建立的一系垂向悬架系统阻尼比的目标函数J(ξ),利用Matlab,求解得到方程<img file="FDA0000796318740000015.GIF" wi="204" he="141" />的正实数根,即μ=2;III步骤:根据I步骤所建立的一系垂向悬架系统阻尼比的目标函数J(ξ),及II步骤中求得的μ=2,建立一系垂向悬架系统阻尼比目标函数J(ξ)的解析表达式,即:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>J</mi><mrow><mo>(</mo><mi>&xi;</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msup><mi>Gv&pi;</mi><mn>2</mn></msup><msubsup><mi>K</mi><mn>1</mn><mn>2</mn></msubsup><mrow><mo>(</mo><mn>1</mn><mo>+</mo><mn>27</mn><msup><mi>&xi;</mi><mn>2</mn></msup><mo>)</mo></mrow></mrow><mrow><mi>m</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mi>&xi;</mi></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000016.GIF" wi="566" he="170" /></maths>(4)确定一系垂向悬架系统的最佳阻尼比ξ<sub>oc</sub>:根据步骤(3)中所建立的一系垂向悬架系统阻尼比目标函数J(ξ)的解析表达式,利用MATLAB,求解<img file="FDA0000796318740000017.GIF" wi="212" he="142" />的正实数根,便可得到一系垂向悬架系统的最佳阻尼比ξ<sub>oc</sub>;(5)确定一系垂向悬架系统的最大许用阻尼比ξ<sub>oa</sub>:①确定轨道高低不平顺大小幅值参数G:A步骤:根据步骤(1)中所建立的高速轨道车辆1/4车体二自由度行驶垂向振动微分方程,及步骤(3)中II步骤计算得到的μ=K<sub>01</sub>/K<sub>1</sub>=2,利用Matlab/Simulink仿真软件,构建高速轨道车辆1/4车体二自由度行驶垂向振动仿真模型;B步骤:根据步骤(4)中设计得到的最佳阻尼比ξ<sub>oc</sub>,<img file="FDA0000796318740000021.GIF" wi="349" he="100" />及A步骤中所建立的高速轨道车辆1/4车体二自由度行驶垂向振动仿真模型,以轨道高低不平顺随机输入为输入激励,仿真得到车体垂向运动的振动加速度均方根值<img file="FDA0000796318740000022.GIF" wi="161" he="72" />C步骤:根据步骤(4)中设计得到的最佳阻尼比ξ<sub>oc</sub>,<img file="FDA0000796318740000023.GIF" wi="349" he="93" />及A步骤中所建立的高速轨道车辆1/4车体二自由度行驶垂向振动仿真模型,以利用轨道高低不平顺随机输入简化函数<img file="FDA0000796318740000024.GIF" wi="285" he="132" />所合成的轨道高低不平顺随机信号为输入激励,对车体垂向运动的振动加速度均方根值进行仿真;D步骤:利用步骤B中所得到的车体垂向运动的振动加速度均方根值<img file="FDA0000796318740000025.GIF" wi="160" he="74" />及步骤C中仿真所得到的车体垂向运动的振动加速度均方根值<img file="FDA0000796318740000026.GIF" wi="155" he="70" />建立轨道高低不平顺大小幅值参数的优化设计目标函数J<sub>min</sub>,即:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>J</mi><mi>min</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><msub><mi>&sigma;</mi><mrow><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>s</mi></msub><mo>_</mo><mi>s</mi><mi>i</mi><mi>m</mi></mrow></msub><mo>-</mo><msub><mi>&sigma;</mi><mrow><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mi>s</mi></msub><mo>_</mo><mi>t</mi><mi>e</mi><mi>s</mi><mi>t</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000027.GIF" wi="529" he="116" /></maths>E步骤:根据A步骤中所建立的高速轨道车辆1/4车体二自由度行驶垂向振动仿真模型,以轨道高低不平顺大小幅值参数为设计变量,以轨道高低不平顺随机输入z<sub>v</sub>为输入激励,利用优化算法求D步骤中所建立轨道高低不平顺大小幅值参数的优化设计目标函数J<sub>min</sub>的最小值,所对应的优化变量即为轨道高低不平顺大小幅值参数G;②一系垂向悬架系统最大许用阻尼比ξ<sub>oa</sub>的设计:根据车辆参数,车体垂向振动加速度的最大设计要求值A,步骤(3)中所建立的一系垂向悬架系统最大许用阻尼比目标函数J(ξ)的解析表达式,及①步骤中求得的轨道高低不平顺大小幅值参数G,建立一系垂向悬架系统最大许用阻尼比ξ<sub>oa</sub>的设计数学模型,即:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mn>27</mn><msup><mi>Gv&pi;</mi><mn>2</mn></msup><msubsup><mi>K</mi><mn>1</mn><mn>2</mn></msubsup><msup><mi>&xi;</mi><mn>2</mn></msup><mo>-</mo><msup><mi>A</mi><mn>2</mn></msup><mi>m</mi><msqrt><mrow><msub><mi>K</mi><mn>1</mn></msub><mi>m</mi></mrow></msqrt><mi>&xi;</mi><mo>+</mo><msup><mi>Gv&pi;</mi><mn>2</mn></msup><msubsup><mi>K</mi><mn>1</mn><mn>2</mn></msubsup><mo>=</mo><mn>0</mn><mo>;</mo></mrow>]]></math><img file="FDA0000796318740000028.GIF" wi="879" he="101" /></maths>利用Matlab计算程序,求解上述关于ξ的方程的正实数根,便可得到一系垂向悬架系统最大许用阻尼比ξ<sub>oa</sub>的设计值;(6)一系垂向悬架最优阻尼比ξ<sub>o</sub>的设计:根据步骤(4)中设计得到的最佳阻尼比ξ<sub>oc</sub>,及步骤(5)中②步骤设计得到的最大许用阻尼比ξ<sub>oa</sub>,利用黄金分割原理,计算得到偏舒适性的一系垂向悬架系统的最优阻尼比ξ<sub>o</sub>,即:ξ<sub>o</sub>=ξ<sub>oc</sub>+(1‑0.618)(ξ<sub>oa</sub>‑ξ<sub>oc</sub>)。
地址 255086 山东省淄博市高新技术产业开发区高创园A座313室