发明名称 一种轮毂电驱动汽车乘坐舒适性分析方法及电机选型应用
摘要 本发明公开了一种轮毂电驱动汽车乘坐舒适性分析方法及电机选型应用,通过建立轮毂电驱动汽车的车身与车轮两自由度四分之一车辆振动系统模型,推导出车身加速度相对路面输入速度的幅频特性曲线的不变点,通过建模仿真分析指出轮毂电驱动汽车乘坐舒适性。同时公开了一种轮毂电驱动汽车乘坐舒适性分析方法在电机选型上的应用,分析得出由电机引起的非簧载质量增加是引起轮毂电驱动汽车乘坐舒适性变差的原因,并由不变点频率表达式给出轮毂驱动电机选型对电机质量的要求,改善高频段的乘坐舒适性。
申请公布号 CN106407507A 申请公布日期 2017.02.15
申请号 CN201610740307.4 申请日期 2016.08.26
申请人 吉林大学 发明人 郭洪艳;刘风;申忱;陈虹;郑琦;吴港
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 长春吉大专利代理有限责任公司 22201 代理人 杜森垚
主权项 一种轮毂电驱动汽车乘坐舒适性分析方法,其特征在于,包括以下步骤:步骤一、建立路面、车身与车轮两自由度四分之一车辆振动系统模型:1.1)建立时域滤波白噪声路面输入系统数学模型:<maths num="0001"><math><![CDATA[<mrow><mover><mi>q</mi><mo>&CenterDot;</mo></mover><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mn>2</mn><msub><mi>&pi;n</mi><mn>00</mn></msub><mi>u</mi><mi>q</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mn>2</mn><msub><mi>&pi;n</mi><mn>0</mn></msub><msqrt><mrow><msub><mi>G</mi><mi>q</mi></msub><mrow><mo>(</mo><msub><mi>n</mi><mn>0</mn></msub><mo>)</mo></mrow><mi>u</mi></mrow></msqrt><mi>W</mi><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001094096660000011.GIF" wi="910" he="94" /></maths>其中:n<sub>00</sub>为下截止空间频率,n<sub>00</sub>=0.011m<sup>‑1</sup>;q(t)为路面随机高度位移,m;W(t)为均值为零的高斯白噪声信号;1.2)建立车身与车轮两自由度四分之一车辆振动系统模型,并列出坐标原点选在各自平衡位置时的运动方程,所述运动方程如下:<maths num="0002"><math><![CDATA[<mrow><msub><mi>m</mi><mn>2</mn></msub><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mn>2</mn></msub><mo>+</mo><mi>c</mi><mrow><mo>(</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>-</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001094096660000012.GIF" wi="708" he="63" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>m</mi><mn>1</mn></msub><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mn>1</mn></msub><mo>+</mo><mi>c</mi><mrow><mo>(</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>1</mn></msub><mo>-</mo><msub><mover><mi>z</mi><mo>&CenterDot;</mo></mover><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><mi>k</mi><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><msub><mi>z</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mi>t</mi></msub><mrow><mo>(</mo><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><mi>q</mi><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mrow>]]></math><img file="FDA0001094096660000013.GIF" wi="948" he="62" /></maths>其中:m<sub>1</sub>和m<sub>2</sub>分别是车轮质量和车身质量,kg;z<sub>1</sub>、<img file="FDA0001094096660000014.GIF" wi="38" he="54" />和<img file="FDA0001094096660000015.GIF" wi="40" he="53" />分别是车轮位移、车轮速度和车轮加速度,m、m/s、m/s<sup>2</sup>;z<sub>2</sub>、<img file="FDA0001094096660000016.GIF" wi="46" he="55" />和<img file="FDA0001094096660000017.GIF" wi="43" he="55" />分别是车轮位移、车轮速度和车轮加速度,m、m/s、m/s<sup>2</sup>;C为减震器阻尼系数,N/(m/s);k和k<sub>t</sub>分别为悬架刚度和车轮刚度,N/m;q为路面不平度,m;步骤二、计算分析不变点大小、位置:将所述步骤1.2)列出的运动方程变形后进行拉氏变换:m<sub>2</sub>s<sup>2</sup>z<sub>2</sub>(s)+m<sub>1</sub>s<sup>2</sup>z<sub>1</sub>(s)+k<sub>t</sub>[z<sub>1</sub>(s)‑q(s)]=0;定义三个传递函数:<maths num="0004"><math><![CDATA[<mrow><msub><mi>H</mi><mi>A</mi></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><msub><mover><mi>z</mi><mo>&CenterDot;&CenterDot;</mo></mover><mn>2</mn></msub><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mfrac><mo>=</mo><mfrac><mrow><msub><mi>sz</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>q</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001094096660000018.GIF" wi="454" he="146" /></maths><maths num="0005"><math><![CDATA[<mrow><msub><mi>H</mi><mrow><mi>S</mi><mi>D</mi></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mn>2</mn></msub><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub></mrow><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mfrac><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mn>2</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>s</mi><mi>q</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA0001094096660000019.GIF" wi="701" he="142" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>H</mi><mrow><mi>T</mi><mi>D</mi></mrow></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mn>1</mn></msub><mo>-</mo><mi>q</mi></mrow><mover><mi>q</mi><mo>&CenterDot;</mo></mover></mfrac><mo>=</mo><mfrac><mrow><msub><mi>z</mi><mn>1</mn></msub><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow><mo>-</mo><mi>q</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow><mrow><mi>s</mi><mi>q</mi><mrow><mo>(</mo><mi>s</mi><mo>)</mo></mrow></mrow></mfrac></mrow>]]></math><img file="FDA00010940966600000110.GIF" wi="662" he="139" /></maths>其中:H<sub>A</sub>、H<sub>SD</sub>和H<sub>TD</sub>分别是车身加速度传递函数、悬架弹簧动挠度传递函数和轮胎弹性动挠度传递函数;s=jω;综合上述公式后得到关于|H<sub>A</sub>(jω)|的不变点频率表达式:<maths num="0007"><math><![CDATA[<mrow><msub><mi>f</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>&pi;</mi></mrow></mfrac><msqrt><mfrac><msub><mi>k</mi><mi>t</mi></msub><msub><mi>m</mi><mn>1</mn></msub></mfrac></msqrt><mo>,</mo></mrow>]]></math><img file="FDA0001094096660000021.GIF" wi="316" he="173" /></maths>其中:k<sub>t</sub>为轮胎刚度,m<sub>1</sub>为非簧载质量;步骤三、路面、车辆仿真建模:3.1)在simulink中搭建滤波白噪声形式的时域路面输入模型模块,并做成子模块用于不同路面或不同行驶速度的行驶工况;3.2)根据所述步骤一中建立的车身与车轮两自由度四分之一车辆振动系统模型列出的运动方程在simulink模块中搭建模型,引入所述步骤3.1)中搭建的滤波白噪声形式的时域路面输入模型模块,并给出车辆模型参数,并向workspace中输出目标参数;根据搭建的仿真模型进行仿真,将仿真结果绘制曲线与计算结果对比得出乘坐舒适性分析结果。
地址 130012 吉林省长春市前进大街2699号