多工况下基于子模型的矿砂船舱口角隅形状尺寸耦合优化分析.docx

该【多工况下基于子模型的矿砂船舱口角隅形状尺寸耦合优化分析 】是由【科技星球】上传分享，beplayapp体育下载一共【24】页，该beplayapp体育下载可以免费在线阅读，需要了解更多关于【多工况下基于子模型的矿砂船舱口角隅形状尺寸耦合优化分析 】的内容，可以使用beplayapp体育下载的站内搜索功能，选择自己适合的beplayapp体育下载，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此beplayapp体育下载到您的设备，方便您编辑和打印。多工况下基于子模型的矿砂船舱口角隅形状尺寸耦合优化分析??陈雪,王德禹上海交通大学船舶海洋与建筑工程学院,上海2002400引言货运船甲板设置的大开口破坏了甲板的纵向连续性,且开口角隅处几何形状的突变致使该处产生了较高的应力集中。因此,改善舱口角隅的应力集中成为研究者关注的重点。目前的研究主要集中在角隅的形状优化方面,即通过改变角隅的几何形状来实现应力集中系数的降低。陈震等[1]对某船舶舱口角隅处结构强度进行了有限元细化分析,结果表明,细化分析后的结果远大于粗网格结构强度分析结果。俞铭华等[2]采用边界元法和带自适应移动界限的序列线性规划,对平面应力状态下的船舶甲板舱口角隅进行了形状优化,使舱口应力集中极小化,从而确定了舱口角隅的最佳形状。Okumoto等[3]针对具有椭圆形角隅的方形甲板开口,采用有限元方法计算了其应力集中系数,并据此对某集装箱船的开口角隅进行了优化,但处理的工况简单且单一。高上地等[4]对矩形耐压舱角隅结构进行形状优化,得到了弧形加强角隅形状,通过对形状优化后的角隅结构进行拓扑优化,有效降低了角隅结构的应力集中,并获得了创新性的角隅结构型式。张干锋等[5]结合参数化建模和改进Kriging近似模型,对某过渡肘板和舱口角隅边界进行了形状优化,在保证极小误差和缩小计算成本的情况下,得到了理想的应力分布和重量优化结果。可以看出,上述研究均是对舱口角隅进行形状优化,因而降低了船舶角隅结构的应力集中。通过分析发现,舱口角隅应力集中问题除了和角隅肘板的形状关系密切之外,还和角隅肘板周边的结构有关,例如甲板纵桁、甲板横梁、舱口围板、舱口围肘板等。因此,本文将主要关注某矿砂船的舱口角隅应力集中问题,在对角隅肘板进行形状优化的同时,考虑角隅肘板与周边关联结构的耦合效应;通过提出一种适于优化的子模型选取方法和多工况筛选方法,以及灵敏度分析方法,最终实现矿砂船舱口角隅的形状尺寸耦合优化。:第1种是直接将细化网格模型嵌入三舱段模型中进行局部强度分析,也即直接使用细化的三舱段模型进行局部强度分析;第2种是采用子模型技术创建局部模型来代替原模型进行计算分析。所谓子模型技术,是指从整体模型中截取一部分,并对局部关注区域网格予以细化,然后将整体模型分析的位移结果作为对应子模型的边界条件,并保留对应子模型上的外载荷,重新计算,从而得到更为精确的局部分析结果的一种技术。由于子模型可以在减小计算量的同时模拟结构的真实受力情况,所以在分析局部强度或者进行疲劳分析时,常使用子模型技术。对本文而言,第1种方法虽然计算精确,但计算量较大,并不适合需要进行大量重复计算的优化问题。而子模型技术只需对局部模型进行分析,可以极大地减少计算时间,适于后续的优化分析。所以,本文选用子模型方法进行角隅局部强度分析和后续优化计算。关于子模型技术,已有不少学者进行过相关研究。许允等[6]以散货船底边舱下折角的局部结构细化网格分析为例,分析了子模型范围、边界节点角位移约束、子模型新增节点的插值处理以及子模型边界节点所在位置对应力结果的影响。任慧龙等[7]通过有限元分析软件,建立了SWATH船舷台与支柱体连接处过渡结构关键部位的参数化子模型,并进行了结构优化设计,在满足细网格屈服衡准的情况下降低了结构质量。程远胜等[8]通过子模型技术,对船舶肘板节点结构进行了应力精细化分析,并在此基础上对肘板进行拓扑优化,提出了一种新型的肘板结构形式,有效降低了节点应力集中。但在许多采用子模型进行优化的研究中,都只验证了子模型对强度分析的适用性,而并未验证其对于优化问题的适用性。,因此,边界的选取通常需要进行多次尝试并予以验证。采用子模型方法时,子模型的最小范围为其边界对应于相邻主要支撑构件所在的位置[9]。本文中,角隅结构的形状和尺寸是随着优化迭代而不断变化的,如果切割边界较小,依据子模型原理,那么形状和尺寸的改变势必会对切割边界的位移产生较大影响,即原先加载在切割边界上的位移是不准确的,那么所得优化结果也是不可信的。因此,适于优化的子模型的切割边界应足够远离关注区域,以保证优化区域的形状尺寸改变后,子模型应力分析结果仍在需求的精度范围内。同时,适于优化的子模型的范围应尽可能小,以降低计算成本。因此,本文截取了不同范围的子模型,首先验证子模型方法的正确性以及子模型本身的计算精度,然后,通过对比形状尺寸变化后角隅的最大应力是否与完整模型误差较小,来验证子模型是否适于优化,并提出适于优化的子模型选取方法。。该船总长325m,型深25m,型宽57m,设计吃水18m,,宽17m。采用有限元软件Patran/Natran建模并进行强度分析。图1所示为该矿砂船的三舱段有限元模型。本文主要关注中间第5货舱舱口角隅的应力集中,需对该角隅处进行细化网格有限元分析。原舱段模型的舱口区域为无角隅肘板的粗网格模型,需首先建立椭圆型角隅肘板模型。规范规定,当强力甲板上机炉舱、货舱开口的角隅是抛物线形或椭圆形时,椭圆长轴与舱口宽度之比不得小于1/10,椭圆短轴与舱口宽度之比不得小于1/20[10],因此,本模型角隅的初始设计取为长轴1700mm,短轴850mm。然后,采用50mm×50mm的网格对角隅局部区域进行进行细化。最后,细网格区域内的所有板采用壳单元表示,并将纵骨一维梁单元替换为二维壳单元建模[9],以保证能够真实反映实际结构的力学性能。图2所示为细化网格前、后角隅局部模型。、:子模型a为最小范围的子模型,其边界横向对应于相邻的甲板横梁,纵向对应于相邻的甲板纵桁,其中多出支撑构件的1个单元格仅为方便子模型边界位移加载,如图3(a)所示;子模型b是在最小子模型区域的基础上向外延伸一个强构件间距,如图3(b)所示;子模型c为截取的三框架模型,即模型范围内包含完整的三个横框架,如图3(c)所示。-models本文选取满载轻货迎浪状态下的1_FULL_L_HSM1_SJ载荷工况作为验证工况。首先,在初始设计下,从原模型分析的计算结果中输出该工况下对应子模型边界位置处的六自由度节点位移信息,然后,通过Pcl命令语句将位移数据自动加载到子模型的边界节点上,同时保留子模型范围内相应的外载荷,以此作为子模型的边界条件和载荷。加载位移边界后的子模型b如图4所示。其中,模型初始设计参数为1700,850,42,27,12,15,12mm,该数据从左到右分别为椭圆肘板长轴长a、椭圆肘板短轴长b、肘板及甲板厚度t1、甲板纵桁厚度t2、甲板横梁厚度t3、舱口围肘板厚度t4和舱口围板厚度t5。-modelbafterloadingdisplacementboundaries将3个子模型角隅区域的最大中面应力与原模型的最大中面应力进行对比,结果显示,3个子模型角隅区域的最大中面应力与原模型的一致,均为523MPa,验证了子模型方法的正确性,同时也说明3个子模型均适用于局部强度分析。原舱段模型和子模型的计算结果应力云图如图5所示。,并对比其对子模型精度的影响,来筛选出用于后续优化的子模型。为有效验证,本文选取的2组实验的设计变量分别为第3节耦合优化设计变量的上、下限。-modelsatoriginaldesign为有效验证,本文首先采用拉丁超立方抽样方法,从设计变量的设计域中抽取了20个样本点,加上设计域的上、下限共计22个样本点。然后,始终保持子