基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析.pdf

第4 0 卷第4 期2023年8 月文章编号：10 0 5-0 52 3(2 0 2 3)0 4-0 0 8 6-0 7基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析华东交通大学学报Journal of East China Jiaotong UniversityVol.40 No.4Aug.,2023杜雨恒1,贡哲蓉2,罗秋雨3,刘启洲1,岳琪玻1,肖兵1,彭必友1,刘俊4,杨杰4,陈刚1（1.西华大学材料科学与工程学院，四川成都6 10 0 39;2.阿里云计算有限公司，浙江杭州310 0 30；3.西华大学西华学院，四川成都6 10 0 39；4.四川攀钢嘉德精工科技有限公司，四川成都6 10 0 0 0摘要：零件轻量化是机械制造领域的重要研究方向。以中国大学生方程式汽车大赛BTR-X的悬架立柱作为例，在分析其实际受力情况的基础上，利用Altair inspire form软件以最大化刚度为目标对其进行了减重设计和拓扑结构优化。得到了带刹入弯工况条件下两种不同设计方案的应力、应变值及安全系数，并对比了其优缺点。结果表明，在最佳优化方案中,BTR-X悬架立柱最大有效应力为557.4 MPa，减重18.4%，提升了车辆的操稳性，实现了其轻量化的目标，对进一步提升FSCC赛车的性能有一定的参考价值。最后，通过选择性激光烧结成功实现了优化悬架立柱结构的3D打印成型，为FSCC赛车的轻量化设计及加工提供了新思路。关键词：FSCC赛车；悬架立柱；AltairInspireForm；拓扑优化CAE分析中图分类号：U214;TF142本文引用格式：杜雨恒，贡哲蓉，罗秋雨，等.基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析.华东交通大学学报，2023,40(4):86-92.文献标志码：AOptimal Design and CAE Analysis of FSCC Racing SuspensionColumn Based on Altair Inspire Form SoftwareDu Yuheng,Gong Zherong,Luo Qiuyu,Liu Qizhou,Yue Qibo,Xiao Bing,Peng Biyou,Liu Junt,Yang Jiet,Chen Gang(1.School of Materials Science and Engineering,Xihua University,Chengdu 610039,China;2.Alibaba Cloud Computing Co.,Ltd.,Hangzhou 310030,China;3.Xihua College,Xihua University,Chengdu 610039,China;4.Panzhihua Iron and Stel ResearchInstitute,Chengdu 610000,China)Abstract:The practice of lightweight parts is an important research direction in the field of mechanical manu-facturing.In this study,the suspension column of the Formula Student Combustion China(FSCC)BTR-X wastaken as an example.Based on the analysis of its actual force,the Altair inspire form software was employed toperform weight reduction design and topology optimization with the goal of maximizing stiffness.The stress,strainvalues and safety factors of two different design schemes were obtained under the condition of entering the bendwith brakes,and then their advantages and disadvantages were compared.In the best optimization scheme ob-tained,the maximum effective stress of the BTR-X suspension column was 557.4 MPa,and the weight was re-duced by 18.4%,thus improving the stability of the vehicle.Accordingly,the lightweight goal was achieved,which has certain reference value for further improving the performance of FSCC racing cars.Finally,the 3D收稿日期：2 0 2 2-12-12基金项目：四川省科技厅重点研发项目（2 0 2 2 YFG0366）；四川省粉末冶金工程技术研究中心项目（SC-FMYJ2021-06）第4期printing molding process of the optimized suspension column structure has been successfully realized by selectivelaser sintering(SLS),providing a new route for the lightweight design and processing of FSCC racing cars.Key words:FSCC racing;suspension column;Altair Inspire Form;topology optimization;CAE analysisCitation format:DU Y H,GONG Z R,LUO Q Y,et al.Optimal design and CAE analysis of FSCC racing sus-pension column based on altair inspire form softwareJ.Journal of East China Jiaotong University,2023,40(4):86-92.杜雨恒，等：基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析87中国大学生方程式汽车大赛（formula studentcombustionchina,FSCC）是中国汽车行业水平最高的大学生赛事，旨在培养具有优异创新设计能力、模拟仿真能力及实践动手能力的综合性人才，促进中国制造向中国创造转型1-2 。悬架立柱用于连接车架与车轮，其作用是把作用在车轮上的各种力（如支撑力、驱动力、制动力等）及其力矩传递到车架上3。因此,悬架立柱的结构对赛车的驾驶操稳性有十分重要的影响。基于此，目前已有研究者从悬架结构形式选取、参数计算以及结构优化角度人手，提升悬架的结构和性能 48 。拓扑结构优化的原则是在不影响制件结构和刚度的前提下，重新设计结构的外形和几何形状，以获得最小成本或最小重量。如在汽车领域，研究人员使用拓扑结构优化方法对汽车底板结构进行重新设计，实现了减轻汽车总重量的目标。不仅提高了燃油效率,还有效降低了生产和维护成本。在航空领域，使用拓扑结构优化可获得更轻量化的飞机结构。如波音公司使用拓扑结构优化的方法重新设计了7 7 7 X飞机的机翼结构，有效改善了机翼的空气动力学性能。实现零部件的轻量化设计对提升整个部件的贡献巨大9。拓扑结构优化不仅在航空和汽车等领域已经得到了广泛的应用和研究,取得了重要的成果。同时，在FSC赛车领域，拓扑结构优化可以被用来重新设计底盘结构、车身外壳等部件，以实现更轻量化的设计和更好的性能。如Jiang等利用有限元法对FSC赛车车架的刚度和振动模态进行了分析，优化后的车架减重7.18%l0。关亮亮等首先利用CA-TIA对悬架立柱进行了三维建模，之后利用ANSYS对悬架立柱进行了静力学分析，发现前后立柱的最大应力均小于材料的屈服极限45 5 MPa，安全系数为3.3,满足使用需求。Shi等利用MATLAB设计了FSC赛车的转向梯形机制，减少了轮胎的磨损，保证了良好的转向能力和抓地力。Yang等计算了后轮轮毂在不同工况下的力学参数，结果表明优化设计后的后轮轮毂可减重约2 0%13 。Ye等利用Hypermesh分析了汽车转弯时所需的承载构件扭转刚度，并优化了零部件的结构，发现模拟横向碰撞的最大变形为14.5 7 mm，最大应力为13 6 MPa，可保证驾驶员的安全 14。AltairInspireForm软件可用于子优化设计，并模拟制造工艺过程，引起了越来越多轻量化设计研究者的关注 15-18 。本文以第十届FSCC赛车BTR-X的悬架立柱作为研究对象，利用CATIA建立起三维模型，之后对赛车参数和悬架立柱受力进行了分析，得到符合实际的工况条件。在此基础上，采用Altair InspireForm软件，以满足强度、刚度要求为设计目标，完成了悬架立柱的应力分析及优化重构，并实现其结构轻量化。1赛车主要参数本文所用的FSCC赛车主要参数如表1所示。表中参数除轴距、轮距在设计与制造装配过程中能够较精确保证和控制之外，其余参数均为目标参数。表1BTR-X主要参数Tab.1 Main parameters of BTR-XParameterValueTotal mass245 kgCurb weight180 kgWheelbase1 560 mmFront track1 160 mmRear track width1 120 mmCentroid heightMinimum turningLength2.700 mmWidth1 350 mmParameterHeightFront overhangRear overhangFront and rearaxle load ratioradiusValue1 080 mm740 mm370 mm45:55280 mm3.3 m882悬架立柱受力计算在FSCC赛事中,为提高每圈的赛道成绩，赛道行驶设置很多弯，但在进弯前又不可把速度降到允许的最大弯速度附近，所以参赛选手会普遍用到带刹人弯 19 。带刹人弯也称为“循迹刹车”,带着刹车人弯，使赛车重心前移，前轮抓地力提升，最大限度地逼近最大转弯速度，到了急转弯处刚好减速到位，缩短过弯的时间。可见，带刹人弯是大学生方程式赛车中较为常见的一种车辆行驶工况，也是较为恶劣的一种工况。因此，本文选择的受力工况为带刹人弯，在此基础上分析前悬立柱的受力情况,如表2所示。表2 前悬零件受力及参数表Tab.2Forces and parameters of front suspension partsParameterFFLLL3TBTBDaFFL4L682.1制动卡钳安装座受力计算FFBuy工FBuyFBUL3FBLyFBIz图1卡钳安装座受力简图Fig.1 Force diagram of caliper mounting seat华东交通大学学报图1为卡钳安装座受力简图，在冲击因数为1的情况下，由受力平衡可得FBuz+FBL:=FBF L3-FBuyL2=0FBuz=FBlz=Fr/2=2 053.2 NFruy=Fely=(Fg Ls)/L2=1 231.93 N2.2轴承座受力计算LIOFH区FBFBlzFFValueParameter1 455.3 NF1787.9 N39 mm70.8 mm21 mm79.74 mm47.25 mm1 874.9 N4.52.516.9 N11.6 mm70.6 mm202023年Fa区区FBlzValueF2524.9 N2Lo50.8mmL270 mmR225mmTHa50 mmF2 461.1 NF4 106.4 NB3F1 453.3 NLs60.5 mm5.43.6FBLyFu图2 Z轴方向上轴承座受力简图Fig.2 Force diagram of the bearing housing along theZ-axisFFLL区FBox图3 X轴方向上轴承座受力简图Fig.3 Force diagram of the bearing housing along theX-axis图2 和图3 为轴承的受力简图，根据受力平衡条件计算可得-F,L-FBo Lp+F,r+Fg-Lo=OFro,=Fboz=(F,r+FgLo-F/L1)/L,=16 143.3 N-F,Li+Feox:Li=0Feox=FBo=(F L,)/L,=2 624 NF+FBoz-FB-FBL,=0FBL=FBL=F,+FBoz-Fg=13 912 NF,+FeLx-FBo=0Fex=FeL=Feox-F,=1 186.7 NF,-Fi=0Fo区FBk区第4期3悬架立柱优化AltairInspire内部搭载了功能强大的Altair杜雨恒，等：基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析Fiy=Fuy=F,=2 461.1 N89上OptiStruct求解器,是一个经过工业验证的线性和非线性静力学及振动力学求解器，已广泛应用于工业结构设计及优化设计领域如基于应力、疲劳的优化等。一般来说优化设计主要包括优化目标、工况条件及约束条件3 个要素，优化目标不同对应的数学模型也不同。利用AltairInspire软件进行仿真驱动设计分为以下几个步骤：草绘或输人模型；简化零件；设置相应的材料和载荷；产生理想的形状；确认性能；将概念设计输出为CAD几何模型。拓扑结构优化的基本原理是基于有限元分析的最小势能原理，即对于一个特定的载荷和边界条件，通过在结构中去除不必要的材料，最大限度地减少其势能，并使结构保持足够的刚度和强度。本研究以悬架立柱的静力学刚度最大化为目标，以设计空间体积分数为约束条件，其优化原理的数学语言如下 2 0 FindX=(xi,x2,x3,x;/TeminC(X)=-UKU2NS.:Zv(a.)Vi=1F=KU,0 xminx;xmx1,(i=1,2,.,N)式中：X为设计变量，2 为设计变量的集合，C为柔度,U为结构位移,K为刚度矩阵,F为载荷矩阵,V（x）为体积约束函数,V*为体积限制分数值，i为单元变量，xmin为设计变量下限值，xmx为设计变量上限值。基于Altair Inspire的优化机理，本研究的悬架立柱优化原则和思路如下。降低主销偏移距：由于选用的赛车车轮的偏心距较小，因此在设计前立柱时，需尽量使上、下臂与立柱的铰接点在不与其他零件发生干涉的前提下尽量靠近车轮中心线，以减小转向主销偏移距和转向所需要的力矩；采用轻质材料降低质量：选用钛合金，以减小立柱的整体质量；由于几何空间干涉原因，在保证转向推杆安装位置的前提下，没有多余的空间安装支耳连接立柱和转向推杆，因此只对悬架立柱支耳和立柱体进行结构优化。(a)Load constraintdistribution(c)First principal(d)First principal(e)Safety factorstress图4未优化前悬架立柱的有限元分析Fig.4Finite element analysis of unoptimized frontsuspension pllar首先对未优化前的悬架立柱进行了CAE分析,如图4所示。按照前述受力计算结果，分别在相应位置施加载荷,如图4(a)所示。由于悬架立柱上下两部分的受力情况与形状特征不同，所以将其分解成上下两部分，并采用不对称的方式，如图4(b)所示。设置其分析网格尺寸大小为1.5 mm，并选用TC4钛合金材料，其杨氏模量为116.5 GPa，屈服强度10 2 9 MPa。利用Inspire软件得到有限元模拟(1)结果,如图4(c)图4(e)所示。前立柱的最大应力值为2 10 MPa,最大应变值为0.0 0 2,但考虑到进一步提升整个FSC赛车的性能，需对悬架立柱进行拓扑结构优化,以期达到减重的目的。(2)在Inspire的优化模块中选择拓扑优化，以支架的刚度最大化为设计目标，以设计空间的总体积为设计变量，以无频率约束、最小厚度为4.5 mm、滑动接触为约束条件，质量目标为设计空间的2 0%。在Inspire软件中,厚度约束指的是优化后结构的最小厚度，其值越小优化精度越高，但所需时间也越久。此外，厚度约束值与后台的网格数量密切相关，厚度约束值越小则后台生成的网格数量也越多。Inspire软件默认后台网格尺寸大小为厚度约束的1/3。在得到优化结果后,对该结果进行了分析验证，如图5 所示。由结果可知，轻量化后的悬架立柱具有典型的拓扑结构，其主应力最大为147 7 MPa,最小安全系数为0.5 7 2 3。对方案一的优化结构进行了包覆处理，采用Inspire自带的PolyNURBS模块进行建模，结合多边形建模的简便性与NURBS的精度和灵活性，对方案一的结构进行重构。根据方案一的仿真优化参(b)Design spacedecompositionstrain90(a)Topology optimization(c)First principal strain图5 方案一优化结果及其CAE分析Fig.5Optimization results of scheme 1 and its CAEanalysis数，进行局部优化得到了方案二，并将两种方案的位移值、安全系数、第一主应力进行对比，结果如图6 所示。方案一最大位移量为0.5 7 mm，而方案二最大仅为0.40 mm，有效抑制了2 9.8%的变形。由主应变云图得方案一的最大主应变为0.0 0 48 1，方案二为0.0 0 48 7,最大主应变差别不大。由应力值变云图得方案一的最大应力为6 3 7.4MPa，方案二为557.4MPa,表明经过局部结构调整，有效降低了悬华东交通大学学报架立柱的应力值，提高了结构的稳定性。悬架立柱选用的材料为TC4钛合金，其屈服极限为10 2 9 MPa，经过结构优化后的等效应力值小于材料的屈服强度，保证了悬架立柱的使用安全性。对悬架立柱的两种结构优化方案进行了静力学强度分析校核，安全系数n用如下公式进行表示(b)Safetyfactorn=Omax式中：，为材料屈服强度，mx为最大应力值，n，为安全系数。一般来说对于钛合金,其最小安全系数为2.8。对比两种方案，大部分区域的安全系数均大于3，表明两种优化方案基本安全。但部分区域的安(d)First principal stress全系数较小，需在使用过程中监测其疲劳状态，如方案一和方案二的最小安全系数分别为1.2 9,1.7 8。4悬架立柱的3 D打印加工采用选择性激光烧结3 D打印对优化后的悬架立柱进行了加工成型,如图7 所示。图7(a)为悬架立柱的三维模型,其相应的STL模型如图7(b)所示。图7(c)为打印制件的切片过程,将原始的悬架立柱的三维模型离散成多个二维平面，其切片参数如下：烧结温度为15 9，预热温度140,扫描间距0.3 mm，扫描速度为2.7 m/s，打印层厚为0.1mm。选用华曙高科的HS403P设备进行3 D打印加工，材料选用尼龙2023年=ns(3)(a)Smoothingeffect diagram(b)Displacementvalue cloud diagramof scheme 1(c)Displacementvalue diagram ofscheme 2(a)The 3D model of thesuspension column(b)The STL model of thesuspension column(d)Safety factorcloud diagram ofscheme 1(e)Safetyfactor of scheme2 Cloud diagram(f)Scheme 1strain valuecloud diagram(c)The slicing processof the part(d)The digital photo ofthe 3D printing process ofthe suspension column10134567(g)Scheme 2(h)Scheme 1 strainvaluestress valuecloud diagramcloud diagram图6 方案一和方案二优化结果分析Fig.6 Analysis of optimization results of scheme 1 andscheme 2(i)Scheme 2stress valuecloud diagram(e)Front view of theprinted part图7 优化悬架立柱的3 D打印加工Fig.7 Optimizing the 3D printing process of the suspensioncolumn(f)Top view of thesuspension column第4期12,其烧结过程图片如图7(d)所示，可以看到整个打印过程稳定，成型缸中无翘曲现象，这归因于尼龙12良好的烧结性能。经过3 D打印的低维构建、层层叠加,最终得到了优化后的悬架立柱实物，如图7(e)和图7(f),可以看到打印制件形状保持较好且外观规整，可为其最终的成型提供模型参考。需注意的是在本文中所使用的尼龙12 屈服强度仅为45 MPa,难以满足悬架立柱的实际使用工况，因此若考虑实际使用还需使用金属3 D打印技术如激光选取熔化等。5结论本文以FSCC赛车BTR-X的悬架立柱作为研究对象，利用Altair Inspire Form对原模型结构进行拓扑优化设计及数值模拟分析，对比了两种优化设计方案，并根据最优方案对零件模型进行重构，结论如下。1）优化后的零件相较于悬架立柱传统切削加工方案，零件的安全系数、最大应力值、最大应变值均得到了改善。重构后的零件总重为47 3 g，比传统切削加工方案降低了10 6.7 g，减重比达18.4%，实现了悬架立柱轻量化的目标。此外，在各项设计参数均达标的前提下，减轻了簧下质量，提升了车辆的操稳性。2）本文选用打印精度较高的选择性激光烧结技术,以尼龙12 为打印材料，成功地打印出了优化后的悬架立柱，为赛车的轻量化设计及加工提供了参考。参考文献：1王淑旺,夏麒翔,罗建辉.大学生方程式赛车电子差速控制策略研究 J.合肥工业大学学报（自然科学版）,2 0 2 2，45(3):289-293.WANG S W,XIA Q X,LUO J H.Research on electronicdifferential control strategy of Formula Student racing carJ.Journal of Hefei University of Technology(Natural ScienceEdition),2022,45(3):289-293.2张瑞乾.关于大学生方程式赛车车队的组织与管理的探讨 J.教育教学论坛,2 0 19(2 2):7-8.ZHANG R Q.Discussion on the organization and managementof Formula Student racing teamJJ.Education and TeachingForum,2019(22):7-8.3王世权，任豪放.大学生方程式赛车整车优化设计 J科技与创新,2 0 2 0(16):44-46.WANG S Q,REN H F.Vehicle optimization design of Formula Student racing carJ.Technology and Innovation,杜雨恒，等：基于Inspire的FSCC赛车悬架立柱优化设计及CAE分析nology,2021,48(2):87-94.5张璇，韩忠浩，曾飞云，等.FSC赛车空气动力学套件的设计与分析 J.辽宁工业大学学报：自然科学版，2 0 15,3 5(5):308-312.ZHANG X,HAN Z H,ZENG F Y,et al.Design and analy-sis of FSC racing aerodynamic kitJ.Journal of LiaoningUniversity of Technology:Natural Science Edition,2015,35(5):308-312.6阎力，史青录,连晋毅.FSC赛车车架的静态结构与模态分析 J.太原科技大学学报，2 0 17,3 8(2)：9 8-10 3.YAN L,SHI Q L,LIAN J Y.Static structure and modalanalysis of FSC racing car frameJ.Journal of Taiyuan Uni-versity of Science and Technology,2017,38(2):98-103.7马芳武,陈实现,赵花垒，等.强度,刚度与模态约束下的FSC赛车车架轻量化 J湖南大学学报（自然科学版），2018,45(4):18-25.MAF W,CHEN S X,ZHAO H L,et al.Lightweight FSCracing frame under strength,stiffness and modal constraintsJ.Journal of Hunan University(Natural Science Edition),2018,45(4):18-25.8唐亦鹏,张嘉杰,刘新田,等.FSC 赛车踏板总成设计与改进 J.上海工程技术大学学报，2 0 19,3 3(3)：2 5 6-2 6 1.TANG Y P,ZHANG J J,LIU X T,et al.Design and improve-ment of FSC racing pedal assemblyJ.Journal of Shanghai U-niversity of Engineering Technology,2019,33(3):256-261.9胡国良，吴礼繁.磁流变制动器结构设计研究现状分析J华东交通大学学报,2 0 2 0,3 7(5)：1-8.HU G L,WU L F.Analysis on the research status of magnetorheological brake structure designJ.Journal of EastChina Jiaotong University,2020,37(5):1-8.10 JIANG L,WANG G,GONG G,et al.Proceedings of theFISITA 2012 World Automotive Congress,November 27-29,2012Cj.Beijing:Springer,2013.11关亮亮，王本善,王思朋，等.小型赛车悬架立柱有限元分析 J.汽车实用技术,2 0 19（17）：47-48.GUAN L L,WANG B S,WANG S M,et al.Finite elementanalysis of suspension column of small racing carJJ.Auto-motive Practical Technology,2019(17):47-48.12 SHI C Z,SHI Z X,WANG T T.Design of steering trape-zoidal mechanism for FSC racing base on matlab J.Ad-vanced Materials Research,2013(647):885-890.13 YANG Y,YANG S,ZHAI Y.Proceedings of the 2022 5thInternational Conference on Intelligent Robotics and Con-trol Engineering(IRCE),September 23-25,2022CJ.Tian-jin:IEEE,2022.912020(16):44-46.4刘明,朱瑞峰,乔晓亮.基于ANSYS的FSAE赛车仿真分析 J.应用科技,2 0 2 1,48(2):8 7-9 4.LIU M,ZHU R F,QIAO X L.Simulation analysis of FSAEracing car based on ANSYSJ.Applied Science and Tech-9214 YE P,ZHOU W,ZHAO C.Design and lightweight basedon FSC racing monocoqueJ.Journal of Physics ConferenceSeries,2021,1885(5):052031.15 SRINIVAS G L,JAVED A.A novel method to synthesize asingle topology for dynamically loaded membersJ.Journalof Mechanical Science and Technology,2021,35(4):1549-1557.16罗勇，杜平，朱丽君，等.基于Inspire软件的拓扑优化设计案例分析 J.制造技术与机床,2 0 2 1(11):3 1-3 4.LUO Y,DU P,ZHU L J,et al.Case analysis of topologyoptimization design based on inspire softwareJ.Manufac-turing Technology and Machine Tools,2021(11):31-34.17唐锋，熊建武，胡智清，等.摩托车金属连接件ABS轻量化设计与FDM增材制造 .工程塑料应用，2 0 2 2,5 0(3)：86-90.TANG F,XIONG J W,HU Z Q,et al.ABS lightweight de-sign and FDM additive manufacturing of motorcycle metalconnectorsJ.Engineering Plastics Application,2022,50(3):86-90.18赵阳,陈敏超,王震.面向增材制造的索杆结构节点拓扑优化设计 .建筑结构学报，2 0 19,40(3)：5 8-6 8.ZHAO Y,CHEN M C,WANG Z.Topology optimizationdesign of cable-strut structure nodes for additive manufac-turingJ.Journal of Building Structures,2019,40(3):58-68.19冯红晶.基于Cruise的电动车整车性能参数匹配及仿真分析 J.西华大学学报（自然科学版）,2 0 17,3 6（3）：16-22.FENG H J.Performance parameter matching and simula-tion analysis of electric vehicle based on CruiseJ.Journalof Xihua University(Natural Science Edition),2017,36(3):16-22.华东交通大学学报20周冉，沈嵘枫，谢诗妍，等.以Altairinspire对育苗成型机纵封部支架的拓扑优化 J.河南工程学院学报（自然科学版）,2 0 2 0,3 2(2):7 1-7 3.ZHOU R,SHEN R F,XIE S Y,et al.Topology optimiza-tion of vertical sealing part support of seedling formingmachine with Altair inspireJ.Journal of Henan Instituteof Engineering(Natural Science Edition),2020,32(2):71-73.第一作者：杜雨恒（19 9 5 一）,男，硕士研究生，研究方向为增材制造及智能装备。E-mail：49 7 0 7 3 3 2 7 q q.c o m。通信作者：陈刚（19 9 1一），男，副教授，博士，硕士生导师,研究方向为增材制造及CAE仿真模拟。E-mail:。（责任编辑：熊玲玲）2023年