针对某电动轻便摩托车车架存在不满足垂直工况下静强度要求的问题,基于相对灵敏度分析,提出了电动轻便摩托车车架的轻量化方法。在建立电动轻便摩托车车架有限元模型的基础上,结合相对灵敏度分析方法确定了车架的优化设计变量。建立了电动轻便摩托车车架轻量化的优化模型。基于厚度尺寸的寻优策略,获得了电动轻便摩托车车架的最优结构参数。研究结果表明,优化后车架最大Von Mises应力比优化前的最大Von Mises应力值降低了51.15%,车身重量降低了1.11 Kg。优化后的摩托车车架还能满足电动轻便摩托车的时速、动力、续航能力的使用需求。 In response to the issue of insufficient static strength of the frame of a certain electric lightweight motorcycle under vertical conditions, a lightweight method based on relative sensitivity analysis for the frame of the electric lightweight motorcycle was proposed. After a finite element model estab-lished for the frame of the electric lightweight motorcycle, an optimization model was established and the optimal design variables of the frame were determined. Employing a thickness dimension optimization strategy, the optimal cross-sectional dimensions were sought. The research results show that the maximum Von Mises stress value of the optimized motorcycle frame is reduced by 51.15% compared with it before optimization, and the body weight is reduced by 1.11 kg. The op-timized motorcycle frame can also meet the needs of electric moped for speed, power, and endur-ance.
针对某电动轻便摩托车车架存在不满足垂直工况下静强度要求的问题,基于相对灵敏度分析,提出了电动轻便摩托车车架的轻量化方法。在建立电动轻便摩托车车架有限元模型的基础上,结合相对灵敏度分析方法确定了车架的优化设计变量。建立了电动轻便摩托车车架轻量化的优化模型。基于厚度尺寸的寻优策略,获得了电动轻便摩托车车架的最优结构参数。研究结果表明,优化后车架最大Von Mises应力比优化前的最大Von Mises应力值降低了51.15%,车身重量降低了1.11 Kg。优化后的摩托车车架还能满足电动轻便摩托车的时速、动力、续航能力的使用需求。
电动轻便摩托车,车架,垂直工况,轻量化,Von Mises应力
Daijian Xu*, Yelong Huang, Han Mei, Feng Cheng#
Liuzhou Wuling Automobile Industry Co., Ltd., Liuzhou Guangxi
Received: Nov. 10th, 2023; accepted: Dec. 15th, 2023; published: Dec. 21st, 2023
In response to the issue of insufficient static strength of the frame of a certain electric lightweight motorcycle under vertical conditions, a lightweight method based on relative sensitivity analysis for the frame of the electric lightweight motorcycle was proposed. After a finite element model established for the frame of the electric lightweight motorcycle, an optimization model was established and the optimal design variables of the frame were determined. Employing a thickness dimension optimization strategy, the optimal cross-sectional dimensions were sought. The research results show that the maximum Von Mises stress value of the optimized motorcycle frame is reduced by 51.15% compared with it before optimization, and the body weight is reduced by 1.11 kg. The optimized motorcycle frame can also meet the needs of electric moped for speed, power, and endurance.
Keywords:Electric Moped, Frame, Vertical Working Condition, Lightweight, Von Mises Stress
Copyright © 2023 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
2020年,中国政府提出实现碳达峰和碳中和的目标,使安全、绿色环保逐渐成为社会的主题,交通运输是碳排放的重要领域之一。随着我国经济的蓬勃发展和人们生活水平的提高,居民的出行需求逐渐增加,城市公共交通系统的线路限制和私家车的高成本导致居民的出行需求无法完全满足。在这种情况下,以电池为动力的电动摩托车就备受普通居民的青睐,逐渐发展成为道路交通中一项重要的交通工具,是加速行业向绿色低碳转型、推动交通运输高质量发展的重要抓手 [
在对车架的优化研究中主要采用有限元分析和优化设计相结合的方法,即根据分析结果对机械结构进行优化设计:清华大学的扶原放教授考虑了多种行驶工况的冲击载荷对车架的破坏作用,基于结构可靠性和有限元法对车架的结构参数进行了可靠性优化设计 [
已公开发表的参考文献表明,摩托车车架的轻量化研究中设计变量的选取存在盲目性,导致计算效率低的问题。本文使用基于相对灵敏度分析的车架轻量化方法,以车架强度为主要约束条件,可针对性地确定优化设计变量并能快速地将优化设计变量收敛到最优值附近,大幅度地提高优化效率。
依照新的国标规范,按摩托车设计的电轻摩车架结构的强度、质量会富裕,可能导致不必要的经济损失。因此,在满足市场需求的同时,也要达到减小排放和降低功耗的目标。将优化设计中优化目标函数定义为电轻摩车架的质量最小。
M g = min ∑ i = 1 N m ( x g i ) (1)
式(1)中, M g 为电轻摩车架的总质量, m ( x g i ) 为组成电轻摩车架的第i个零部件的质量,N为对应的电轻摩车架零部件设计变量的个数。
电轻摩车架由Q235材料制成,其屈服极限为235 Mpa,故选择的约束条件为:
231 Mpa ≤ σ ≤ [ σ y ] = 235 Mpa (2)
考虑电轻摩车架各部件的位置、结构和管材厚度等因素将优化问题中设计变量定为7个,如表1所示 [
序号 | 零部件名 | 初始值 单位(mm) |
---|---|---|
1 | 转换器支架 | 1.5 |
2 | 边管支架、中撑主管、边管、主梁管 | 2 |
3 | 中撑限位板 | 2.5 |
4 | 线夹 | 0.8 |
5 | 后平叉梁 | 1.8 |
6 | 车头立管、后轮泥瓦支架 | 3 |
7 | 侧面加强板 | 2.5 |
表1. 电轻摩车架轻量化分析的优化设计变量
进行有限元线性静态优化分析时,结构的性能参数 T j 对其设计参数 x i 的灵敏度定义如下 [
S e n ( T j x i ) = ∂ T j ∂ x i (3)
约束函数和目标函数都可以是通过求解静力平衡方程位移得到的响应,可表示为 T = T ( δ ) ,同时位移是设计变量的隐式函数,可用 δ = δ ( X ) 表示,则
T = T ( δ ( X ) ) (4)
静力平衡方程如式(5)所示:
[ K ] { δ } = [ F ] (5)
并对第i项设计变量 x i 进行偏微分处理,并移项得
[ K ] ∂ { δ } ∂ x i = ∂ [ F ] ∂ x i − ∂ [ K ] ∂ x i { δ } (6)
等式两边同乘 [ K ] − 1 得到质量对设计变量的灵敏度:
∂ { δ } ∂ x i = [ K ] − 1 [ ∂ ( F ) ∂ x i − ∂ [ K ] ∂ x i { δ } ] (7)
并将上式进行微分处理,可表示为:
Δ δ = ∑ i ∂ { δ } ∂ x i Δ x i = [ K ] − 1 ∑ i [ ∂ [ F ] ∂ m i − ∂ [ K ] ∂ m i { δ } ] Δ m i = [ K ] − 1 ( Δ F − Δ [ K ] { δ } ) (8)
由于静态分析中载荷向量F不随设计变量的变化而变化,即 Δ F = 0 , Δ δ = − [ K ] − 1 Δ [ K ] { δ } ,因此节点位移函数对应的性能参数(包括约束和目标函数)对于设计变量 x i 的灵敏度为
∂ T j ∂ x i ≈ Δ T j Δ x i = T j ( X 0 + Δ x i , δ 0 + Δ δ ) − T j ( X 0 , δ 0 ) Δ x i (9)
综上,本研究的优化设计数学模型为:
{ 求 x i ( i = 1 , 2 , ⋯ , n ) 使 得 M g = ∑ i = 1 k m ( x g i ) → min s .t . 231 Mpa ≤ σ ≤ [ σ y ] = 235 Mpa (10)
式中 x i 为车架零部件的壁厚。
研究对象是优化约束问题,其寻优策略为 [
x ( v + 1 ) = x ( v ) + α ( v ) s ( v ) (11)
式中,向量 x ( v ) , x ( v + 1 ) , s ( v ) ∈ E n , x ( v ) 为当前设计点, x ( v + 1 ) 为下一个设计点, α ( v ) 表示数量, s ( v ) 表示方向。
采用一维搜索方法,并将目标负梯度作为寻优方向,即 s = − ∇ f 。
代入到优化模型后得到
{ 求 ( x 1 , ⋯ , x n ) 使 f ( x 1 , ⋯ , x n ) → min s .t . g j ( x 1 , ⋯ , x n ) = 0 ( j = 1 , ⋯ , J ) h k ( x 1 , ⋯ , x n ) ≤ 0 ( k = 1 , ⋯ , K ) x l ≤ x ≤ x u (12)
得到:
{ 求 α 使 F ( α ) = f ( x 0 − α ∇ f ) → min s .t . G j ( α ) = g j ( x 0 − α ∇ f ) = 0 ( j = 1 , ⋯ , J ) H k ( α ) = h k ( x 0 − α ∇ f ) ≤ 0 ( k = 1 , ⋯ , K ) α l ≤ α ≤ α u (13)
其中,目标梯度采用差商代替
∇ f ( x ) = { ∂ f ∂ x 1 ⋮ ∂ f ∂ x m ⋮ ∂ f ∂ x n } ≅ { f ( x 1 + Δ x 1 , ⋯ , x m , ⋯ , x n ) Δ x 1 ⋮ f ( x 1 , ⋯ , x m + Δ x m , ⋯ , x n ) Δ x m ⋮ f ( x 1 , ⋯ , x m , ⋯ , x n + Δ x n ) Δ x n − f ( x 1 , ⋯ , x n ) Δ x 1 ⋮ − f ( x 1 , ⋯ , x n ) Δ x m ⋮ − f ( x 1 , ⋯ , x n ) Δ x n } (14)
通常电轻摩车架的结构复杂,其中相当一部分结构对仿真结果的影响微乎其微。为提高仿真速度、增强工作效率和实用性,在不影响计算结果准确性的前提下依据相应的准测对车架进行简化 [
图1. 电轻摩车架有限元模型
如图1所示,取车架行驶前进的方向为x轴,驾驶员侧向为y轴,竖直方向为z轴并建立坐标系。三种工况对应的边界条件如表2所示 [
工况 | 约束位置 | 约束 |
---|---|---|
水平工况 | 前轮轴 | y、z轴平动和x、z轴转动自由度 |
后轮轴 | 3个平动和x、z转动自由度 | |
座垫工况 | 前轮轴 | y、z轴平动和x、z轴转动自由度 |
后平叉轴 | 3个平动和x、z转动自由度 | |
垂直工况 | 前轮轴 | y、z轴平动和x、z轴转动自由度 |
后平叉轴 | 3个平动和x、z转动自由度 |
表2. 电轻摩车架三种工况下边界条件
基于静强度有限元模型,并考虑了设计变量、目标函数等多个优化信息得到的电轻摩车架优化模型如图1所示,主要设计变量参数如表3所示。
计算水平、座垫、垂直三种工况下电轻摩车架的静态应力分布,并得到相应的最大应力,如图2所示。在垂直工况下,车架的最大Von Mises应力超过了材料Q235的屈服极限,是进行轻量化设计的主要工况。
图2. 多工况下电轻摩车架静态最大Von Mises应力
序号 | 零部 件名 | 设计 变量 | 初始值 单位(mm) | 优化值 (mm) | 序号 | 零部 件名 | 设计 变量 | 初始值 单位(mm) | 优化值 (mm) |
---|---|---|---|---|---|---|---|---|---|
1 | 垂直锁支架 | x1 | 1.20 | 10 | 后平叉梁 | x10 | 1.80 | 1.72 | |
2 | 转换器支架 | x2 | 1.50 | 1.49 | 11 | 中撑主管 | x11 | 2.00 | 1.52 |
3 | 边管支架 | x3 | 2.00 | 1.28 | 12 | 边管 | x12 | 2.00 | 1.28 |
4 | 中撑限位板 | x4 | 2.50 | 1.89 | 13 | 主梁管 | x13 | 2.00 | 1.28 |
5 | 中轴管 | x5 | 3.50 | 14 | 车头立管 | x14 | 3.00 | 1.92 | |
6 | 边撑触地板 | x6 | 4.00 | 15 | 中轴支架 | x15 | 3.00 | ||
7 | 边撑U型板 | x7 | 6.00 | 16 | 后平叉管 | x16 | 3.00 | 5.75 | |
8 | 铭牌支架 | x8 | 1.00 | 17 | 后轮泥瓦支架 | x17 | 3.00 | 1.92 | |
9 | 线夹 | x9 | 0.80 | 0.97 | 18 | 侧面加强板 | x18 | 2.50 | 5.94 |
表3. 电轻摩车架主要设计参数的初始值和优化值
应用相对灵敏度分析方法,计算垂直工况下各项参数的灵敏度系数后,根据确定的参数信息进行后续的优化工作,其结果如表3所示。结果表明序号1、5、6、7、8和15零件质量的改变对电轻摩车架应力结果没有显著的影响;增加序号2、9、10、16、和18零件的质量会降低电轻摩车架的应力值,减少序号3、4、11、12、13、14和17零件的质量会降低电轻摩车架的应力值,是优化设计中着重关注的设计变量。
图3. 垂直工况下电轻摩车架应力对质量的相对灵敏度值图
由图2的三种工况可知,在垂直工况中,电轻摩车架最大Von Mises应力远超材料的屈服极限,为477 Mpa,须在不改变原车架结构的基础上且满足静强度的条件下进行车架的轻量化设计。在研究中,根据图3的计算结果确定优化设计变化量为x2、x3、x4、x9、x10、x11、x12、x13、x14、x16、x17、x18等12个设计变量。电轻摩车架厚度迭代终止于第15次。优化前后其设计变量值如表3所示。
图4的计算结果表明,优化过程中电轻摩车架质量不断地发生改变。车架的质量急剧地下降到21.89 Kg,随后缓慢逐渐地下降到21.5 Kg,经过八次迭代后,车架的质量上升到21.62 Kg,又在21.56 Kg左右起伏,升至21.71 Kg后在第14次迭代到21.61 Kg,经过第15次迭代,曲线趋于稳定,接近至优化目标。此时,在垂直工况中的电轻摩车架最大Von Mises降低至233 Mpa,小于材料的屈服极限。其质量由初始的22.72 Kg降低至21.61 Kg,质量减轻了约5%,达到了较为明显的优化效果。各设计变量在优化前后的值的变化情况如表3所示。
图4. 垂直工况下电轻摩车架质量随迭代次数变化曲线
在前面的优化设计后通过静强度分析对车架性能的提升进行验证。图5的垂直工况下仿真结果表明,优化后,电轻摩车架的Von Mises应力降至材料的屈服极限以下,相较于优化前下降了51.15%。
图5. 优化前后电轻摩车架的Von Mises应力云图。(a) 优化前电轻摩车架的应力云图;(b) 优化后电轻摩车架的应力云图
阶数 | 优化前 频率(Hz) | 优化后 频率(Hz) | 阶数 | 优化前 频率(Hz) | 优化后 频率(Hz) | |
---|---|---|---|---|---|---|
1 | 30.52 | 20.424 | 6 | 112.22 | 73.794 | |
2 | 36.259 | 26.498 | 7 | 121.84 | 94.183 | |
3 | 62.589 | 34.952 | 8 | 133.26 | 102.34 | |
4 | 83.562 | 67.391 | 9 | 136.49 | 107.15 | |
5 | 99.458 | 71.985 | 10 | 148.94 | 116.34 |
表4. 优化前后前10阶电轻摩车架的固有频率
进行电轻摩车架轻量化设计后,计算了优化前后前10阶电轻摩车架的固有频率,计算结果如表4所示。表4的数据表明,电轻摩车架质量的改变导致前10阶固有频率也发生相应的改变,相同阶次下自由模态频率均有降低。路面不平度是车辆振动的关键外部激励源之一,有必要进一步利用典型路面实测功率谱密度的频率分量,来确定优化后车架的固有频率成分是否落入在输入的电轻摩振动的时间频率范围内,验证优化结果是否可行。
当路面不平度的空间频率为Ω (次/米),电动轻便摩托的行驶速度为V (米/秒)时,输入的时间频率f (次/秒)是Ω和V的乘积,即
f = V × Ω (15)
若电轻摩振动的时间频率变化为与电轻摩车架的模态频率一致,电轻摩车架就会诱发共振,其共振车速为
V = 3.6 L W f ( km / h ) (16)
表5为在前两阶模态频率下,电轻摩在不同路面上行驶时的共振车速。
路面 | 未铺装路 | 碎石路 | 搓板路 | 平坦公路 | |
---|---|---|---|---|---|
道路不平度波长LW(m) | 0.77~2.5 | 0.32~6.4 | 0.74~5.6 | 1~6.3 | |
共振车速(km/h) | 20.424 Hz | 56.52~183.82 | 23.53~470.57 | 54.41~411.75 | 73.53~220.58 |
26.498 Hz | 73.5~238.48 | 30.53~610.51 | 70.59~534.2 | 95.39~600.98 |
表5. 电轻摩行驶在不同路面上的共振车速
GB/T24158-2018国标要求,电轻摩的行驶速度 ≤ 50 km,因此表5计算结果表明除碎石路外在我国典型路面上行驶,路面激励不会诱发车架共振,表明优化后电轻摩车架满足使用性能要求。
开展了基于相对灵敏度分析的电轻摩车架的轻量化研究,得到以下结论:
1) 相对灵敏度分析法能有针对性地确定优化设计变量,提出的电轻摩车架结构的优化模型能快速地将优化设计变量收敛到最优值附近,大幅度地提高了优化效率。
2) 优化后的电轻摩车架结构的质量减轻了1.11 Kg。轻量化的电轻摩车架最大Von Mises应力低于材料的屈服极限,相较于优化前的应力降低了51.15%,同时确保其满足行驶要求。
徐代俭,黄业龙,梅 寒,程 峰. 基于相对灵敏度分析的电动轻便摩托车车架的轻量化研究Research on the Lightweight of Electric Mo-ped Frame Based on Relative Sensitivity Analysis[J]. 机械工程与技术, 2023, 12(06): 528-538. https://doi.org/10.12677/MET.2023.126058