电动汽车车架的多目标拓扑优化-尺寸优化-精细化设计.pdf

以某电动旅游观光车桁架式车架为例进行多目标优化，先后采用了基于梁单元的多次局部选代拓扑优化方法、基于梁单元的横截面尺寸优化方法和基于壳单元车架强度分析的精细化设计方法。具体步骤为，首先建立基于梁单元的车架有限元模型，以车架的刚度性能为优化目标建立多目标优化模型，通过多次局部选代拓扑优化确定梁结构的布置位置，然后以梁的横截面长、宽、壁厚作为优化变量进行尺寸优化设计，确定车架梁的横截面尺寸；在此基础上，建立基于壳单元的车架有限元模型进行强度分析，基于强度分析结果对车架的高应力区域进行精细化设计，最终优化后的车架质量降低了1 0.5%，性能至少提升了9.3%，并且直接可以用于生产制造。该优化设计方法综合运用了梁单元模型和壳单元模型分析的优势，分别应用于概念设计阶段和工程设计阶段，为桁架式车架的优化提供参考。关键词：车架；梁单元；壳单元；拓扑优化；尺寸优化；精细化设计中图分类号：TH16;U463.32Multi-Objective Topology Optimization,Size Optimization andDetailed Design of Frame for an Electric CarFANG Chao,ZHAO Shang-yi,YAN Cel,QIN Guo-feng*?(1.Liuzhou Wuling Automobile Industry Co.,Ltd.,Guangxi Liuzhou 545027,China;2.Post-Doctoral Research Center,Guangxi Automobile Group Co.,Ltd.,Guangxi Liuzhou 545027,China;3.Teachers College for Vocational and Technical Education,Guangxi Normal University,Guangxi Guilin 541004,China)Abstract:A truss frame of the electric tourism vehicle was taken as an example to carry out multi-objective optimization.The to-pology optimization method based on multiple local iterations of beam element,size optimization method of cross section sizebased on beam element and the detailed design method based on the strength analysis of shell element were used successively.Theconcrete steps were as follows:firstly,the finite element model of the frame was established based on the beam element,and amulti-objective optimization model was established to optimize the frame stiffness performance.The layout position of the beamstructure was determined by several local iterative topology optimization,and then the size of the beam cross section was deter-mined by the size optimization design by taking the length,width and wall thickness of the beam as optimization variables.Onthis basis,the finite element model of the frame based on the shell element was established for strength analysis,which was usedfor the detailed design of the frame.The mass of the final optimized frame reduced by 10.5%,while the performance improved byat least 9.3%,and the frame was able to be directly used for manufacturing.The optimal design method proposed comprehensive-ly utilizes the advantages of beam element model for conceptual design stage and shell element model for detailed design stage,which provides reference for the optimization of truss frame.Key Words:Frame;Beam Element;Shell Element;Topology Optimization;Size Optimization;Detailed Design文献标识码：A文章编号：1 0 0 1-39 9 7(2 0 2 3)0 8-0 0 1 6-0 71 引言车架是汽车关键部件之一，是其它系统安装的基础，在汽车行驶过程中承受各种方向的载荷。客车、旅游观光车等桁架式车架通常由矩形钢管、圆钢管、槽钢、工字钢等工业标准型材焊接或铆接而成，具有着承载能力强、工艺简单、成本低等优点，但车架自来稿日期:2 0 2 2-0 6-1 1基金项目：国家自然科学基金项目（5 2 1 0 2 4 7 3)；广西科技基地和人才专项（桂科AD20159010)作者简介：秦国锋，（1 9 9 0-）,男，四川广安人，博士研究生，副教授，主要研究方向：汽车结构设计理论与轻量化技术研究身质量大，占汽车总质量的(2 0 4 0)%l1-3，根据有关资料显示，当燃油汽车整车质量减少1 0%，燃油消耗将相应减少（6 8)%；电动汽车整车质量减少1 0 0 kg,续航里程提升(6 1 1)%2 ,因此，车架轻量化设计对于节能减排具有重要意义。车架轻量化方法主要有采用轻质高强材料、先进加工制造技术、车架结构优化设计三第8 期种方法4 。优化设计方法主要有拓扑优化、形状优化、形貌优化和尺寸优化四种5 ,其中拓扑优化和尺寸优化分别为概念设计和详细设计，在车架结构优化设计方面应用最为广泛。目前，在对客车等桁架式车架进行优化时，通常采用梁单元或者壳单元车架有限元模型。在使用梁单元方面，文献7 提出了一种改进遗传算法能够对桁架式梁单元车架进行拓扑优化。文献8 针对多材料客车车架，提出一种基于梁单元的拓扑优化和截面尺寸优化方法。文献9 基于NSGA-算法采用梁单元对客车底架进行了拓扑优化。文献1 0 采用梁单元对车架进行了拓扑优化，然后分别对横截面半径和壁厚进行了尺寸优化。在使用壳单元方面，文献!基于折衷规划法和SMP密度函数插值模型，采用壳单元对客车车架进行了多工况拓扑优化。文献1 2-1 4 在灵敏度分析的基础上，以壳单元客车车架的壁厚作为设计变量进行了多目标尺寸优化设计。文献1 5 采用壳单元基于局部拓扑优化、灵敏度优化和尺寸优化对电动客车车架进行了多目标优化。梁单元建立的车架有限元模型易修改、计算速度快，尺寸优化时，结构梁截面主要参数作为设计变量，如矩形钢管的宽、高、厚三个截面参数，也可以指定若干截面参数作为变量，其具有较高的灵活性，变量优化空间大，但简化处理过多，无法反映连接处位置的应力大小及分布情况，适合用于概念设计；壳单元具有模拟精度高的优点，能真实反映车架应力分布情况，但在尺寸优化时，壳单元模型的只能将壁厚作为设计变量，变量单一，优化效果有限,因此壳单元适合用于结构合理性验证1 6 。现有的研究以梁单元或者壳单元进行优化都存在一定的局限性，并不完全能发挥出车架优化的潜能。综合运用基于梁单元和壳单元的仿真分析模型，以某电动旅游观光车为研究对象进行多目标优化设计，建立梁单元车架模型进行多次局部迭代拓扑优化，然后以梁单元车架的横截面长、宽、壁厚为设计变量进行尺寸优化，在此基础上建立壳单元车架模型，提出基于强度分析的精细化设计，最终完成车架的轻量化设计。2基于梁单元的原车架性能分析2.1模型建立梁单元模型是在几何线框的基础上创建，首先需要对车架三维模型进行简化，如图1 所示。图1 原车架简化模型Fig.1 Simplified Model of Original Frame采用hypermesh软件建立基于梁单元的车架有限元模型，赋予单元材料和属性，车架结构材料皆为Q235A钢材，其材料属性，如表1 所示。网格大小平均为5 mm，共节点连接，基于梁单元的车架结构3D显示，如图2 所示。方朝等：电动汽车车架的多目标拓扑优化-尺寸优化-精细化设计弯曲刚度分析位置左前悬架右前悬架约束Z载荷施加位置载荷大小位置约束载荷施加位置载荷大小原车架的刚度分析与模态分析结果,如表3所示。其中，模态分析只选取对车架整体影响大的一阶扭转和一阶弯曲，弯曲刚度、扭转刚度分析和模态分析的位移云图，如图3 图6 所示。表3刚度、模态分析结果Tab.3 Stiffness and Modal Analysis Results弯曲位移0.668mmContour PlotDisplacement(Mag)Analysis system6.811E-016.054E:015.297E-014.541E-013.784E-013.027E-012.270E:011.514E-017.568E-021.327E-19NOReSUMax-6.811E-01Gds31304m-1327E-19Fig.3 Analysis Results of Bending Stiffness of Original Frame17表1 材料属性参数Tab.1Material Property Parameters名称数值弹性模量/MPa2.1 104泊松比0.3密度/(kgm=3)7.85 10-6屈服强度/(MPa)235图2 基于梁单元的车架结构3D显示Fig.2 3D Display of Frame Structure Based on Beam Element2.2原车架性能分析首先采用梁单元模型进行车架的刚度和模态分析，确定车架的初始性能，为后期的多目标拓扑优化和尺寸优化提供目标性能。车架刚度分析的载荷及约束，如表2 所示。模态分析采用自由模态。表2 原车架刚度分析的载荷及约束Tab.2 Loads and Constraints of OriginalFrame Stiffness Analysis左后悬架YZX、Z节点编号31 2 9 2节点编号30 9 33-1500-1500扭转刚度分析车架前端横梁左后悬架X.Y.ZX.Y.Z节点编号2 4 8 1 5节点编号2 4 8 1 4-1710-1710扭转位移一阶弯曲2.322mm16.21Hz31292Global2-0.67030933GlobalZ-0.666图3原车架弯曲刚度分析结果右后悬架X、Y、Z右后悬架X、Z一阶扭转22.91HzNo.818机械设计与制造Contour PlotDisplacement(Mag)Analysissystem8.314E+007.390E+006.466E+005.543E+004.619E+003.695E+002.771E+001.848E+009.238E-014.882E-21NOH.DUMax-8.314E+00Grigs 28297Kp-4.882E-21Fig.4 Analysis Results of Torsional Stiffness of Original FrameContour PlotEigen Mode(Mag)Analysissystem4.276E+003.801E+003.326E+002.851E+002.376E+001.901E+001.426E+009.507E-014.757E-01NoResul6.249E-04Max-4.276E+00Grids25590in6.249E-04Grds109Fig.5First Order Torsion Model of Original FrameContourPlotEigen Mode(Mag)Analysis system4.276E+003.801E+003.326E+002.851E+002.376E+001.901E+001.426E+009.507E-014.757E-016.249E:04NoResultMax=4.276E+00Min-6.249E-04Grids25590Fig.6 First Order Bending Model of Original Frame3基于迭代思想的多目标局部拓扑优化3.1 基本理论拓扑优化可以在均匀材料的优化空间内找到最优的材料分布方案，该方案在拓扑优化中体现为“最大刚度 设计1 7 。通常把刚度最大问题等效为柔度最小问题来研究，柔度则用应变能来定义。体积约束下求最小柔度的多工况拓扑优化问题可表示为1 8 :find:p=pi,P2,pnsin:m(p)=W,c(p),V=JV。subject to:0pminpii=1,2,N式中：p一第i个有限单元的相对密度；pmn一设计变量的取值下限；N一有限单元的数量；Cmin一总柔度；C,和W一结构在第j工况下的柔度和权重因子；V和V。一结构体积和初始结构体积；F一约束因子。3.2选代思想迭代思想是指重复反馈过程的活动，其目的通常是为了逼近所需目标或结果。每一次对过程的重复称为一次“迭代”,而每一次迭代得到的结果会作为下一次迭代的初始值。桁架结构的车架由多种不同截面尺寸和形状的梁组成，在对Aug.2023车架进行拓扑优化时，由于拓扑基础模型过于复杂，在给定的物理条件下，设计区域里无法拓扑出有效的结构，通过添加制造工艺约束(成员尺寸、对称约束、挤出约束)也无法解决，如图7(a)所示。因此，本次研究将采用迭代的思想解决该问题，根据拓扑优化Global结果,先对材料密度接近1 且结构形状呈现明显的设计域添加结24814GlobaZ-2.326图4 原车架扭转刚度分析结果图5 原车架一阶扭转模态图6 原车架一阶弯曲模态j=1构，根据工程经验来确定结构截面形状及截面尺寸，改变模型的传力路径，然后以此作为下一次拓扑优化的基础模型,再次进行拓扑优化,直至拓扑优化结果中没有明显的结构形状呈现,因此，便可在设计区域拓扑出具有参考意义的结构，如图7(b)所示。Sub(a)Fig.7 Topology Optimization Diagram3.3拓扑优化分析在满足总布置设计和造型的前提下,其余梁所在区域则可以作为设计区域，首先去掉设计区域的梁单元，然后在设计区域建立壳单元，如图8 中车架的蓝色所示。图8 车架结构拓扑优化模型Fig.8 Topological Optimization Model of Frame Structure拓扑优化定义如下：优化目标：车架弯曲刚度和扭转刚度最大；优化约束：体积分数上限为0.3；设计变量：指定设计域的壳单元的材料密度制造工艺约束：避免拓扑优化结果中出现一些细小的传力路径，添加了成员尺寸控制，最小成员尺寸定义为1 5，最大成员尺寸设置为8 0,并添加对称约束，保证拓扑优化后车架左右材料的分布对称。本次研究车架的拓扑优化结果中，较多设计区域存在不确定性、无参考价值的结构，无法通过一次拓扑优化获取理想状态下的最优结构布置方案，所以采用迭代的思想完成车架的拓扑优(b)图7 拓扑优化示意图No.8Aug.2023化。首先将拓扑优化结果材料密度设置为0.9 5，选定优化结果中,结构呈现最明显的设计区域,在模型里去掉该设计域的材料，添加结构，如图9 所示，完成一次迭代，然后此模型作为下一次拓扑优化的基模型，剩余可设计区域作为设计变量进行拓扑优化。机械设计与制造尺寸优化的数学模型为：find:X=X,i=1,2,.,Nmin mize:f(x)=Msubjectto:GCj,j=,2,N,XIX,XIX,0式中：M一车架总质量；G一模型第j工况的性能参数，如位移、频率应力等；G一模型优化时第j工况的性能参数的约束值；X一设计变量，x一第x,变量的梁截面尺寸，上角标带L为变量下限，上角标带U为变量上限。本次研究的车架为承载式车架，主要由矩形截面梁和少部分异型截面梁组成，本次优化只针对矩形截面梁，其截面形状，如图1 1 所示。将横截面的长、宽、壁厚作为尺寸优化的变量。19厚长图1 1 矩形梁截面Fig.11 Rectangular Section of Beam研究车架的梁截面厚度，主要集中在1.5 mm，考虑焊接工艺、使用寿命等问题，优化后厚度不应小于1 mm。本次优化共有1 7 6 个设计，在优化前设置设计变量的离散值，避免了优化结果无参考价值。尺寸优化定义及优化结果，如表4 所示。由于篇幅原因，梁截面尺寸变化只展示部分。表4 尺寸优化定义及优化结果Tab.4 Dimension Optimization图9 送代优化过程示例Fig.9 Example of Iterative Optimization Process车架经过7 次的迭代，拓扑优化结果中已没有明显的结构形状呈现，也没有材料密度高、具有参考价值的设计区域，车架最终结构布置，如图1 0 所示。图1 0 车架拓扑优化最终结构Fig.10 Final Structure After Topology Optimization4基于梁单元的多目标尺寸优化完成拓扑优化得到的车架结构在一定程度上只反应了材料分布的趋势，不能作为最终设计的结果，需要基于拓扑优化车架上进行尺寸优化设计。尺寸优化通过优化一维梁单元的截面尺寸、二维板单元的厚度等，寻求目标结构件的最佳截面尺寸并能满足相应的性能要求,尺寸优化整个过程中车架的形状拓扑结构没有发生变化1 9 。Definition and Optimization Results名称代号目标质量kg弯曲位移mm扭转位移mm约束第7 阶频率Hz第8 阶频率HzID1-宽ID1-长ID1-厚设计ID2-宽变量ID2-长ID2-厚ID3-厚经过优化，车架从最初的2 4 7.2 kg下降到了1 9 9.1 kg，共减少48.1kg，占总质量的1 9.4%，说明该车架拥有较大的优化空间。5基于壳单元模型强度分析的精细化设计桁架的车架受到载荷时，结构的连接处最为容易出现应力过大,超出连接部位材料的屈服极限值，产生塑性变形、断裂等状下限初值上限结果247.2199.10.670.652.312.3116.2122.022.9123.02030205011.52040204011.513100100460604420301602011No.820机械设计与制造况，而梁单元车架则无法精确分析车架连接处的应力情况，因此需要参考梁单元车架模型的尺寸优化结果，建立壳单元车架模型完成车架的强度分析。对于应力超出材料屈服极限处，若采用尺寸优化，约束连接处应力，通过改变整根梁的截面尺寸来降低连Aug.202310119接处的局部应力,则大大削弱了轻量化效果,因此采用添加角钢、斜撑梁等小结构件增加连接处强度，解决应力集中问题，称为精细化设计。经过精细优化后的车架质量、刚度等将发生变化，而梁单元模型和壳单元模型存在较大的差别,无法做比较，所以在建立拓扑优化和尺寸优化后车架壳单元模型时，也需要建立原车架的壳单元模型进行初始性能分析，用于对比分析。5.1基于壳单元建立拓扑优化-尺寸优化后的车架模型由于车架整体改动较大，所以先使用UG完成拓扑优化-尺寸优化后的车架结构三维建模。因该车架属于量产车型，因此对该车架进行优化时，需要考虑工艺、成本、改动影响范围等条件，在参考拓扑优化-尺寸优化后梁单元车架优化结果建立车架结构三维模型时，对于一些不符合实际要求的结构截面尺寸进行适当修改。对于修改后的三维模型车架建立基于壳单元的优化车架，简化车架构件较小的圆角、孔等特征，然后划分网格，几何形状简单的构件采用四边形网格，复杂的则采用三角形和四边形混合网格，网格平均大小为5 mm,结构连接采用RBE2单元。5.2基于壳单元对拓扑优化-尺寸优化后的车架强度分析在汽车行驶的过程中，典型的工况有扭转工况、制动工况、转向工况，垂向工况，观光车的道路行驶条件情况良好，所以只进行了满载制动工况和满载转向工况的分析，车架上的载荷主要为电池的重量、乘客的重量以及安装在车架上各附件的重量，载荷的重量采用集中载荷的方式添加，电池、附件、乘客都设计有固定的位置，在其质心位置赋予重量，然后使用Rbe3对质心和车架进行连接，以优化车架为例，各质心赋予的重量及位置具体情况，如表5、图1 2 所示表5 汽车各系统质量Tab.5Massof Automotive Systems质心代号12345678910112Fig.12 Centroid Position of Each Load5.2.1 制动工况制动工况是车辆进行紧急制动时，分析车架各结构所产生的应力分布情况，通过施加一个纵向的加速度和重力加速度的来进行模拟2 0 ,加速度和约束条件，如表6 所示。由图1 3知，车架制动工况产生的最大应力值为34 2.5 MPa，超过了材料的屈服应力。表6 制动工况加速度和约束条件Tab.6 Accelerations and Constraints of Braking Condition约束位置左前悬架右前悬架左后悬架右后悬架优化车架制动工况分析结果，如图1 3所示。ContourPlotElement Stresses(2D&3D)ivonMises,Max)AnalysissystemSimpleAverage3.425E+023.045E022.664E+022.284E021.903E021.522E+021.142E+027.612E+013.806E+013.737E-13Max-3.425E+02GRds.636940Xdh-3.37E-13ords-414216Fig.13 Stress Cloud Map of Braking Condition After Optimization5.2.2转向工况电动车在转弯时，整个车身在离心力的作用下有向转弯侧倾名称质量kg仪表台附件重量52副驾驶位乘客重量70前车架电池重量144中部电池+第一排乘客重量426中部电池+第二排乘客重量426后车架电池重量72电机控制器10第四排乘客重量210第三排乘客重量210车顶重量42副驾驶位乘客重量703图1 2 各载荷质心位置约束自由度方向X向ZY向Xy/zZ向StaticMax.Value342.547图1 3优化车架制动工况应力云图斜的趋势。通过施加一个横向的加速度和重力的加速度来模拟转向工况2 0 的应力分布情况，加速度和约束条件，如表7 所示。表7 转向工况加速度和约束条件Tab.7 Accelerations and Constraints of Steering Condition约束位置左前悬架右前悬架左后悬架右后悬架优化车架转向工况分析结果，如图1 4 所示。由图1 4 知，优化车架转向工况产生的最大应力值为4 2 5.3MPa，超过了材料的屈服应力。845约束自由度67方向X向Y向Z向加速度大小0.8g-1g加速度大小0.8g-1gNo.8Aug.2023ContourPlotElement Stresses(2D&3DKvonMises.Max)AnalysissystemSimpleAverage4.253E+023.308E+023.780E+022.835E+022.363E021.890E+021.418E+029.451E+014.725E014.856E-13NoReulMax=4.253E+02Gids43713GMdsm1348n-4.856E-13机械设计与制造ContourPlotDisplacement(Mag)Analysissystem1.094E+019.720E+008.505E+007.290E+006.075E+004.860E+003.645E+002.430E+001.215E+001.014E-16Max=1.094E+01425.284Grds273603X-1.014E-16Grids85754721图1 7 优化后车架扭转刚度分析位移云图图1 4 优化车架转向工况应力云图Fig.14 Stress Cloud Map of SteeringCondition After Optimization5.3基于强度分析的精细化设计在优化车架中，对应力值超过原车架应力的结构位置处添加斜撑梁、角钢等构件来分散应力集中，此时新车架制动工况和转向的最大应力极限2 2 2.2 MPa和2 1 4.0 MPa，达到了车架优化时，其性能不降低的要求。添加斜撑梁前、添加斜撑梁后应力变化示意图,如图1 5 所示。(b)添加斜撑梁后应力变化图1 5 添加斜撑梁前、添加斜撑梁后应力变化示意图Fig.15 Schematic Diagram of Stress ChangesBefore and After Adding Inclined Beam6优化后车架性能分析对优化车架的弯曲刚度和扭转刚度进行了分析，得出平均位移分别为0.9 4 9 mm和3.5 7 7 mm，如图1 6、图1 7 所示。Contour PlotDisplacement(Mag)Analysis system9.850E-018.756E-017.661E-015.472E-016.567E-014.378E-013.283E-012.189E-011.094E-011.520E-15NoResiltMax 9.850E-01Grigds635028M-1.520E-15Gfds.072035图1 6 优化后车架弯曲刚度分析位移云图Fig.16 Displacement Cloud Map of BendingStiffness Analysis After OptimizationFig.17 Displacement Cloud Map of TorsionalStiffness Analysis AfterOptimization对优化后的车架进行自由模态分析，优化车架一阶扭转频率为1 7.36，一阶弯曲频率为1 8.31,如图1 8、图1 9 所示。本次研究的电动汽车的行驶速度最高为5 0 km/h，平时一般行驶为30 km/h，即8.3m/s，输人时间频率得1 6.6 Hz，优化车架的一阶扭转和一阶弯曲频率皆高于1 6.6 Hz，避开了路面激励引起车架的共振。contourPlotEigen Mode(Mag)Analysissystem4.702E+004.184E+003.666E+003.148E+002.631E+002.113E+001.595E+001.077E+005.591E-014.121E-02Max-4.702E+00Grads85760463690s20576-4:121E-02StaticMax.Value=342.547(a)添加斜撑梁前1145163StaticMax.Value=222.201GlobalZ=-0.951Subca图1 8 优化车架一阶扭转模态Fig.18 First Order Torsion Model After OptimizationContourPlotEigenMode(Mag)Analysissystem3.748E+003.332E+002.917E+002.501E+002.086E+001.670E+001.255E+008.396E-014.241E-01NORSu8.705E-03Max-3.748E+00Grigls 400467Gnds.84757Xig-8.205E-03Fig.19 First Order Bending Mode After Optimization由表8 可知，车架在刚度性能不下降情况下，应力最大值在材料屈服极限值内，以及模态分析一阶扭转和一阶弯曲频率避开了路面激励，降低2 5 kg质量，约1 0.5%，达到了轻量化的目的。表8 车架优化前后性能对比Tab.8 Comparison of Frame PerformanceBefore and After Optimization优化车架质量kg237弯曲位移mm1.047扭转位移mm4.985制动工况应力MPa255.4转向工况应力MPa382.2一阶扭转Hz10.50一阶弯曲Hz16.227结论采用拓扑优化、尺寸优化和精细化设计对某电动观光车车架进行优化设计，最终实现质量降低了1 0.5%，性能至少提升了9.3%，达到了较好的轻量化效果。Subcase3(图1 9 优化车架一阶弯曲振型原车架2120.9493.577222.2214.017.3618.31变化率-10.5%+9.3%+28.2%-12.9%-44.0%+65.3%+12.8%No.822机械设计与制造(1)采用基于梁单元的多次局部迭代拓扑优化方法，优先确chassis frame based on NSGA-IlJ.Journal of Tongji University:Natu-定材料密度高、路径分布清晰的结构位置，经过多次迭代，逐步ralscience,2017,45(11):1664-1669.)确定设计区域的结构分布,具有较为理想的拓扑优化效果。10高超.拓扑优化驱动的某纯电动汽车车身骨架设计研究D.秦皇岛：燕山大学，2 0 1 4.(2)在拓扑优化的基础上，采用基于梁单元的车架尺寸优化(Gao Chao.A study of electric vehicles body skeleton structure design方法，能够对梁横截面的长、宽、高尺寸进行优化，能够大幅度降driven by topology optimizationDJ.Qinhuangdao:Yanshan University,低车架质量。2014.)（3)对拓扑优化和尺寸优化后的梁单元车架，考虑实际生产11范文杰，范子杰，桂良进，等.多工况下客车车架结构多刚度拓扑优化制造过程中的工艺、成本等因素，建立壳单元车架进行强度分析，设计研究J.汽车工程，2 0 0 8(6)：5 31-5 33.能够很好地反映车架接头处的应力分布，在此基础上进行精细化(Fan Wen-jie,Fan Zi-jie,Gui Liang-jin,et al.Multi-stiffiness topol-设计，能够降低布局应力。ogy optimiration of bus frame with multiple loading conditions JJ.Auto-参考文献motive Engineering,2008(6):531-533.)12 Fu C L,Bai Y C,Lin C,et al.Design optimization of a newly developed1陈瑞.复杂形状断面的薄壁梁车身骨架结构优化及程序设计D.长aluminum-steel multi-material electric bus body structure J.Struc-春：吉林大学，2 0 2 0.tural and Multidisciplinary Optimization,2019,60(5):2177-2187.(Chen Rui.Structural optimization and code design of thin-walled auto-13】龙江启，李毅，俞平，等.基于结构优化的客车骨架轻量化设计研究motive body frame with complex cross sectionD.Changchun:Jilin Uni-J.机械科学与技术,2 0 1 3,32(4):5 5 8-5 6 3.versity,2020.)(Long Jiang-qi,Li Yi,Yu Ping,et al.Study on the lightweight design of2赵笠程.纯电动汽车白车身多目标轻量化优化设计研究D.长沙：湖coach skeleton based on the structura optimization methodJJ.Mechani-南大学，2 0 1 9.cal Science and Technology for Aerospace Engineering,2013,32(4):(Zhao Li-cheng.Research on multi-objective lightweight optimization558563.)design of pure electric vehicle body-in-white D.Changsha:Hunan14】吴胜军，叶欣钰，江兴洋.混合动力客车车身骨架轻量化设计J.机University,2019.)械设计与制造,2 0 2 1(7)：2 2 1-2 2 4.3石朝亮，史建鹏.基于HyperWorks白车身灵敏度分析及结构优化J.(Wu Sheng-jun,Ye Xin-yu,Jiang Xing-yang.Lightweight design of汽车科技,2 0 1 2(6):6 1-6 4.hybrid bus body frameJ.Machinery Design&Manufacture,2021(7):(Shi Chao-liang,Shi Jian-peng.Sensitivity analysis and structure opti-221-224.)mization of the body-in-white by hyperworks J.Automotive Science15毛爱华.纯电动大客车骨架结构轻量化多目标优化设计D.长春：吉and Technology,2012(6):61-64.)林大学，2 0 1 5.4李光雾,刘新玲.汽车轻量化技术的研究现状综述J.材料科学与工(Mao Ai-hua.Multi-objective optimization design for electric bus body艺,2 0 2 0,2 8(5):4 7-6 1.skeleton structure lightweightD.Changchun:Jilin University,2015.(Li Guang-ji,Liu Xin-ling.Literature review on research and develop-16王松.某商用客车车架有限元分析与结构优化D.武汉：武汉科技大ment of automo