超大型集装箱船中横隔板拓扑优化方法研究.pdf

第45卷第1 6 期2023年8 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.16Aug.,2023超大型集装箱船中横隔板拓扑优化方法研究李海洲，刘平1 2，李小灵，童晓旺1，李安琪2(1.江南造船（集团）有限责任公司，上海2 0 1 9 1 8；2.江苏科技大学土木工程与建筑学院，江苏镇江2 1 2 0 0 3)摘要：随着国际货运需求的增加，大型集装箱船的吨位越来越大，为节约建造成本，在设计阶段有必要对集装箱船结构进行优化。本文以超大型集装箱船的板架结构作为研究对象，探索拓扑优化在箱船中的应用。在指定应力强度的约束条件、以结构重量最轻为优化目标下，对中横隔板的几何形状进行优化计算。此外，为了使优化后的结构拓扑形式便于生产、加工，约束条件还增加了轴对称条件。迭代优化后的结果表明，横隔板下部材料可以去除，一方面节约船体空间，有利于管道、电缆的铺设；另一方面结构重量（相对于原横隔板）减少了1 8%，而最大应力只增加1.5%。在箱船横隔板设计中，可以采用相同的方法进行拓扑优化以获得更好的拓扑形状。关键词：箱船；中横隔梁；拓扑优化中图分类号：U674.131文章编号：1 6 7 2-7 6 49(2 0 2 3)1 6-0 0 2 7-0 5Typical topology optimization design of the transversal support structures of(1.Jiangnan Shipyard(Group)Co.,Ltd.,Shanghai 201918,China;2.School of Architecture and Civil Engineering,Abstract:With the study of typical main structure(shell and middle section structure),the application of structural to-pology optimization design method in the design is discussed.The shape of the transverse partition is optimized under theconstraints and the lightest structure weight,in addition,in order to make the optimized structure topology form to facilitateproduction and processing,the constraints add symmetry conditions.The results after iterative optimization show that thelower material of the structure can be removed,on the one hand,save the hull space and facilitate the laying of pipeline andcable;on the other hand,the structure weight(relative to the original transverse partition)is reduced by 28%.It is shown thatin the current transverse partition design,the same method for topological optimization design can be used to obtain excel-lent design forms.Key words:very large container ship;tranvers e support structure;shape optimization0引言优化设计思想很早就有学者提出，早在2 0 世纪60年代，优化设计方法就被一些学者应用于机械、飞机制造、船舶工业等领域。1 9 8 6 年，Moe等1-2 阐述了关于运输船甲板的优化计算方法，并第一次在船舶结构领域应用了数学规划理论。这项研究成果标志着船舶结构设计进入了优化思想阶段3-4。随着半个世纪的发展，针对舱口形状等局部结构的尺寸优化和拓扑收稿日期：2 0 2 2-0 4-0 7基金项目：国家自然科学基金资助项目（51 50 8 2 38）作者简介：李海洲（1 9 8 1-），男，研究员，研究方向为船舶工程结构研发。文献标识码：Avery large container shipLI Hai-zhou,LIU Ping2,LI Xiao-ling,TONG Xiao-wang,LI An-qi?Jiangsu University of Science and Technology,Zhenjiang 212003,China)doi:10.3404/j.issn.1672-7649.2023.16.006设计已经较为成熟5-。Michae7对船体的早期阶段设计优化进行了研究，表明在早期阶段考虑船体的优化可最高节约1 5%的建造成本。2 0 1 2 年，LingCaol8)提出了一种新型拓扑优化算法，并成功应用于船体的横梁剖面。几何优化的基本思想是以材料分布作为优化变量，基于优化对象设定约束条件，然后通过特定算法，找到最优的分布方案9-1 0 。目前，在连续假设基础上兴起的拓扑优化理论被广泛应用于航天工程、机28械制造、汽车工业等领域，也被少量应用于船舶行业-1 3。但是，关于船舶构件拓扑与尺寸优化分析的实际应用非常少。Masafumi等6 研究了发动机振动下的船甲板拓扑优化问题，采用多尺度优化策略分析结构的振动控制，结果表明通过对船体甲板的拓扑形状进行优化，噪声平均可以减少8%。SOGUT等1 41 对船体的拓扑设计与能耗进行研究，结果表明优良的拓扑形状在正常航行状态下可节约5.7%的能耗。孙志军1 5 提出了一种基于拓扑优化技术进行船舶机舱发电机的抗振加固研究，结果表明在减少1.2 t重量的情况下，振动效应减少约6%。钦伦洋1 6 应用变密度法拓扑优化理论，采用AnsysWorkbench对设计区域进行优化计算，并研究了船舶结构的轻量化问题，结果表明在不降低结构强度前提下，船体舱壁的结构用材料最高可减少到原来的50%。张会新1 7 利用拓扑优化方法对船体的振动特性、声学性能进行评估，结果表明在不改变声学约束的情况下，船体振动性能改善有限，但是优化后的重量可减少1 6%。吴贝尼1 8 以船体横剖面为研究对象，研究了多种插值函数的效能，并编程实现了拓扑优化。结果表明，G-BESO方法最为有效，效能最高。本文在大型集装箱船设计初始时假定中横剖面为矩形板壳结构，然后计算其优化结果。探讨结构形状以及拓扑优化方法在船舶领域中的应用前景，并根据实际工程需要及限制条件，给出了优化模型和优化对象的具体拓扑形状。1船体板架结构力学分析模型舰船科学技术距，这样优化结果基本上与边界效应无关。计算几何模型尺寸为：标准纵骨之间距离为0.7 5m；标准舱室间距为4.5m；纵骨腹板厚为1 2.5mm；外甲板厚度为1 7 mm而横舱壁厚为1 2.5m。在甲板层，纵骨面板厚度为36mm，大肘板厚度为1 2 mm而小肘板厚8 mm，扶强材厚1 0 mm，某超大型箱船横隔舱布置如图1 所示。010上第45卷203040m(a)船体轴视图某大型箱船结构构件主要包括外甲板、纵骨和横舱壁以及相关附属结构。大型箱船船体外围尺寸和主要设计尺寸如表1 所示。表1 箱船主要尺寸Tab.1Main dimension of container ship总长/m垂线间长/m计算船长/m型宽/m型深/m设计吃水/m结构吃水/m366355.4本计算采用的箱船坐标系统定义如下：船长为X轴，方向为船行进方向相反；型深方向为Z轴，正方向与重力方向相反；甲板船宽方向为Y轴，正方向由右手坐标系决定。计算模型选取则综合考虑了计算量与简化模型精度的平衡。在实际计算时，在X方向选取了6 个舱段长度，6 个舱段可以充分隔离中横隔剖面边界效应的影响。而在船宽方向选取5个纵骨间(b)中剖面图1 某超大型箱船横隔舱布置图Fig.1 The design of LNG ship计算箱船用材料为普通钢，其力学特性参数如表2所示。在计算中采用弹性方法进行计算，不考虑材料的塑性性能。350.65130.214.516.0表2 材料计算参数Tab.2Material calculation parameters材料密度/kgm3普通钢78502重量目标下船体板架尺寸和形状2.1变密度优化方法将优化对象的密度设定为设计变量，并将其值限弹性模/GPa210泊松比0.3第45卷定在0,1 之间。其中：0 表明材料为空；1 表明材料为真实保留。同时，将材料本身的物理属性与单元密度之间按某种函数关系相联系。这样，把真实结构的拓扑优化问题转化为材料密度分布的数学问题。具体公式如下：p=p1,P2,P3,pnT,Pminpi1.0,Mimimize(pisi),subject to:Q;0近似当作0 来处理，Pmin可取0.0 1。当pi=Pmin时，表明此时材料单元可删除；李海洲，等：超大型集装箱船中横隔板拓扑优化方法研究291）载荷工况板架结构所承受的外荷载主要包括结构自重、水体浮力以及弯矩。由于船舶上建部分在有限元模型中没有体现，因此以荷载的形式加载到横隔板上部；船轴方向的总纵弯曲取2 0 0 kNm，分别加载到横隔板上；水压力取水线在2 0 m位置，加载图3所示。B:Stutie structurntlydrostaticpreUnit:Paimme:1.sMaX:2000011/10/20211/57PM(1)Min:o200.000.077778.0155556.033333.011111.088888.966666.744444.422222.2Min:Max40m当pi=1.0时，表明此时单元必须保留；当0 pi1时，单元的存留需另行判断。由上述计算过程可知，在优化结束后，单元密度为0 或1 最好，为了实现这种目的，在变密度法引了罚函数，即在迭代计算中对单元密度乘以常罚函数，使单元密度趋向于0 或1。如图2 所示，惩罚因子为2 时，迭代次数与材料密度分布关系图。实际上，这种函数是在迭代时使得中间密度向0 与1 状态两侧靠近，迭代次数越多，处于中间值（0 pi1）的密度单元则越少，从而拓扑优化模型则越接近于0 1 状态。1.00.90.80.70.630.30.20.101.00.90.80.70.60.50.40.30.20.10密度分布图2 罚函数迭代次数-变密度分布图（惩罚因子取2）Fig.2 The curve of iteration number vs density(penalty约束条件方面，在弹性设计分析时，一般取在满足结构刚度的条件下，使所有的单元应力值满足;s，,为材料允许应力，即船体所有位置应力均不超过材料的设计强度值。2.2结构计算模型本计算的目标是确保满足应力约束条件下（结构最大应力不超过材料屈服强度）使得结构重量最小，从而实现最优化的经济效益。位移和载荷边界条件如下：.e-009N2）边界条件船体在水中漂浮时重力与水浮力平衡，实际上没有固定约束。为了在有限元中能够进行分析（为约束初始值一次送代二次选代三次送代N次送代number=2)图3水浮力及荷载加载图示Fig.3The load schematic刚体位移），采取了“弱弹簧”技术进行处理，弹簧刚度取1 E-6N/m。在实际计算中，读取弱弹簧的反作用力，取实际作用力的1/1 0 0 0 以下作为依据。船舶设计经济化是最自然的目标，即使船体结构质量最小为优化目标函数。本计算采用AnsysWork-bench软件进行拓扑优化计算。单元尺寸取0.3m，四边形网格为主，共有1 1 6 1 34个结点，1 2 2 42 6 个网格。网格质量分布如图4所示。可以看出，网格质量较好，有限元计算结果可信。优化分析采用如下约束条件进行计算：结构的最大等效应力不大于材料的4238940000Tri3-Quad4300002000010 00000.140.250.350.500.63 0.750.881.00Elementmetrics图4单元质量分布图Fig.4The metric of mesh30屈服强度，即结构应力不大于2 1 5MPa；中横隔板的几何形状约束，确保优化后横隔板几何形状对称，即中横隔板符合生产和施工工艺的要求。3优化过程在开始优化过程之前，首先对箱船结构进行静力分析，这样对船体结构的力学状态有初步了解，同时也可对优化结果有初步预判。3.1静力分析在优化计算前，首先对板架结构进行静力分析，其典型等效应力如图5所示。可以看出，结构大部分区域的等效应力值均小于2 0 3MPa。此外，应力最大值出现在中横剖面下部左下、与右下位置，中间部分应力较小。由此可见，中间部分结构所用材料并没有得到充分利用，有较大的减重潜力。可看出，约束条件为满足强度条件时，可有更好的空间材料分布。利用Ansys Workbench Topology模块对构件进行拓扑优化分析。202.9100181.3200159.7300138.1400116.540094.953073.362051.771030.18008.5894舰船科学技术本文提供的案例在2 0 个迭代之后就已经可以获得较为满意的结果，在一般情况下，结果是否可以接受取决于模型的复杂程度。3.3优化结果最终优化结果如图7 所示。可以看出，绝大部分密度是1.0，在密度小于0.4的材料单元，软件自动把它去除。根据拓扑优化结果，图中横隔板与原来面积相比减少约1 8%。另外，为了加工方面，两侧小块去除的部分一般仍然保留。Remove(0.0 to0.4)Marginal(0.4to0.6)Kcep(0.6to1.0)3.4优化后受力结果对优化后的构型进行静力分析，边界条件与荷载条件与优化之前相同。图8 为优化后区域的等效应力云图。可知，构型最大应力为2 0 6 MPa，与优化前的应力相比较，只增加了约1.5%，优化效果非常显著。第45卷图7 优化后的形状结果Fig.7The mesh after optimization05c+031e+041.5c+042e+04mm图5静力分析时的等效应力大小205.9600Fig.5 The equivalent stress of static analysis184.0900162.230140.3600118.50003.2选代过程96.634074.770052.906031.041.0结构经过优化后，变密度法过程最大密度一直在9.17721.0，但是最小密度则由1.0 趋向于0，迭代过程如图6所示。经过大约1 0 个迭代，最小密度约为0.1，然后随着迭代的深入，最小密度缓慢趋向于0。实际上，01.0000.8750.7500.6250.5000.3750.2507.5403E20Fig.6 The process curve of optimization图8 优化后形状应力结果Fig.8The stress contour after optimization4结语本文采用变密度优化算法，对某大型箱船横隔板1020图6 优化过程曲线304050进行了结构拓扑优化，得出最终的拓扑优化结果。研究结果表明，在指定应力强度的约束条件（最大应力不超过材料强度）以及结构质量最小的目标函数下，第45卷船舶横隔板板架拓扑结构在优化后质量可减小1 8%，而最大应力却只增加了1.5%。说明本文的变密度拓扑优化方法可以为船舶设计提供优化依据。本文的研究结果也可为类似结构件的设计提供方法借鉴。参考文献：【1 程远胜，刘甜甜，刘均.船舶肘板拓扑优化设计.中国舰船研究,2 0 1 5,1 0(5):53-58.CHEN Yuan-sheng,LIU Tian-tian,LIU Jun.Topology opti-mization of ship bracket structureJ.Chinese Journal of ShipReseraver,2015,10(5):53-58.【2 邓林义，林焰，陈明，等.采用拓扑遍历方法优化船舶设计计划.计算机集成制造系统，2 0 0 7,1 3(3)：431-436.DENG Lin-yi,LIN Yan,CHEN Ming,et al.Optimization onship design planning by topological ergodicityJJ.ComputerIntegrated Manufacturing System,2007,13(3):431-436.3张会新，杨德庆.典型船舶板架拓扑与形状优化设计.中国舰船研究,2 0 1 5,1 0(6):2 7-33.ZHANG Hui-xin,YANG De-qing.Typical shape and topolo-gy optimization design of the ship grillage structureJ.Chi-nese Journal of Ship Research,2015,10(6):27-33.【4 尚高峰.轻量化船舶结构极限强度研究D.北京：中国舰船研究院,2 0 1 1.【5王洁松.海上船舶物联网监测拓扑结构可靠性设计与优化分析1 .舰船科学技术,2 0 1 8,40(6):1 7 5-1 7 7.WANG Jiesong.Reliability design and optimization analysis ofmarine ship internet of things monitoring to pologyJ.ShipScience and Technology,2018:175-177.6 DAIFUKU M,NISHIZU T,et al.Design methodology usingtopology optimization for anti-vibration reinforcement of gen-erators in a ships engine roomJ.Proceedings of the Institu-李海洲，等：超大型集装箱船中横隔板拓扑优化方法研究滨：哈尔滨工程大学，2 0 2 0.10 RIGO P,CAPRACE J D.Optimization of Ship Structures.Marine Structures,2017,203(18):1193-1200.11】赵宏旺.多工况应力下船舶多载荷作用下连续体结构拓扑优化 .舰船科学技术,2 0 1 9,41(6)：1-3.ZHAO Hong-wang.Topology optimization of continumstructure under multi-load of ship under multi-conditionstress.Ship Science and Technology,2019,41(6):1-3.12 QU D,LI X,et al.The stress and stiffness analysis of di-aphragmJ.MATEC Web of Conferences,2017,95:1011.13 Zhou J,Wang C,et al.Multi-objective optimization of a springdiaphragm clutch on an automobile based on the non-dominat-ed sorting genetic algorithm(NSGA-II)J.Mathematical andComputational Applications,2016,21(4):47.14 SOGUT O S,ARSLAN O.Improving energy efficiency ofships through optimisation of ship operationsR.University,Istanbul Technical,2014.15】孙志军.基于拓扑优化的船舶机舱发电机抗振加固设计方法D.大连：大连理工大学，2 0 1 9.16 钅钦伦洋.基于拓扑优化的船舶结构轻量化研究D.大连：大连海事大学,2 0 1 6.17】张会新.船舶结构拓扑优化设计及声学评估方法研究D.上海：上海交通大学，2 0 1 6.18】吴贝尼.船舶中横剖面拓扑优化研究及程序系统开发D上海：上海交通大学，2 0 2 0.31tion of Mechanical Engineers,Part M:Journal of Engineeringfor the Maritime Environment,2016,230(1):216-226.7 PARSONS M G.Applications of optimization in early stageship designJ.Ship Science and Technology,2009,3(5):9-32.8 CAO L.Design and optimization of elastic diaphragm inminiature high-temperature pressure sensorJ.Advanced Ma-terials Research,2012,591-593:1260-1264.9】韩松言.基于多智能体的船舶电网拓扑结构优化D.哈尔