船舶有限元模型舱室识别方法.pdf

第45卷第2 0 期2023年10 月舰船科学技术SHIP SCIENCEANDTECHNOLOGYVol.45,No.20Oct.,2023船舶有限元模型舱室识别方法张志康，刘戈，陈乐昆，李闯，张思航，孙明宇!(1.中国船舶及海洋工程设计研究院，上海2 0 0 0 11；2.中国船舶集团有限公司系统工程研究院，北京10 0 0 0 0)摘要：在船舶与海洋结构物设计过程中，一般利用有限元法校核船体结构强度，其中对舱室进行定义的步骤费时费力，还易出错。针对此问题进行研究，提出一种高效的有限元模型舱室边界自动识别算法，称为“无自由边法”。基于FEMAP软件编程实现该算法，并以实船模型为例进行测试。研究表明：该算法高效实用，适用于任意形状复杂的舱室，所开发的程序能够快速、准确地识别出舱室边界的网格单元。关键词：舱室；空间识别；FEMAP；二次开发中图分类号：U662.2文章编号：16 7 2-7 6 49(2 0 2 3)2 0-0 0 40-0 4Research on method of ship compartment recognition with FE modelZHANG Zhi-kang,LIU Ge,CHEN Le-kun,LI Chuang,ZHANG Sihang,SUN Ming-yu(1.Marine Design and Research Institute of China,Shanghai 200011,China;2.CSSC Systems EngineeringAbstract:FEA method is often used for hull structure strength calculation in the ship and marine structure design.Itsnot only time-consuming for the definition of the ship compartment,but also easy to make mistakes.An algorithm helpful torecognize the boundary of the ship compartment efficiently called No Free Edge was proposed to solve this problem.A pro-gram was developed based on FEMAP software,and a practical engineer project was adopted as test example.It can befound that this algorithm was efficient and also suitable for compartments with very complex shapes.The developed pro-gram can recognize the boundary elements of the compartment quickly and accurately.Key words:ship compartment;space detect;FEMAP;secondary development0引言有限元分析是船舶结构设计中非常重要的校核方法，一般包括舱段有限元分析和全船有限元分析。在强度分析时，舱室内部水压力载荷的施加和构件腐蚀扣除的确定均与舱室相关，因此舱室定义是必不可少的环节。手动搜索舱室空间的边界单元，不仅费时费力，效率低下，而且还容易出错。尤其是近年来散货船和油船共同结构规范（HCSR）生效后，规范对散货船和油船的舱段校核范围要求由中间货舱扩展到全船各货舱!，不仅工作量加重，而且船首和船尾区域的线型变化剧烈，给舱室定义增加了难度。此外，对于集装箱船、客滚船和豪华邮轮等需进行全船有限元计算的船型，每条船需定义的各类大小舱室多达上百收稿日期：2 0 2 2-0 8-2 9作者简介：张志康（1991-），男，硕士，工程师，研究方向为船舶数字化设计。文献标识码：AResearch Institute,Beijing 100000,China)doi:10.3404/j.issn.1672-7649.2023.20.007个，非常考验设计人员的体力与耐心。因此，急需实用、高效和智能的舱室自动识别算法和程序，以期能够缩短操作时间，提升设计人员的效率。针对此类舱室空间边界自动识别的问题，国内外学者也进行了不少研究。单威俊等2 4 基于有限元模型，提出“切分拼接法”实现舱室边界单元的识别，并基于Patran开发自动化程序，该程序已经纳人CCS的HSCR规范校核软件体系中，但识别过程仍需用户人为干预，且对船首、船尾等曲面较复杂的区域，识别精度还有待提高。章志兵等5 提出一种基于CAD模型的舱室识别算法，并在NX平台上以C+语言编程实现，算法自动化与智能化程度高，有助于舱室快速建模。章志兵等6 通过设计网格拓扑关系和分析舱室的几何特征，基于半边半面结构，提出了半面第45卷封闭空间概念，并利用半面扩展算法，实现了CAE模型舱室的自动识别，但仅限于封闭舱室，且算法相对复杂。此外，法国船级社的商业软件VeristarHull可以对有限元网格进行舱室搜索，但仅限于水密舱室，内在算法也不得而知。本文基于拓扑学原理和船舶舱室特征，提出一种基于CAE的船舶舱室识别算法“无自由边法”，该算法适用于多种船型及任意复杂的舱室，且改进算法可适用于非封闭舱室。在此基础上开发出自动化程序，经过实船测试验证，能够大幅提高舱室定义效率，具有工程应用价值。1舱室识别船舶舱室最典型的特点为大多是水密舱室，因此本研究以该特点为突破点，从识别有限元三维封闭空间人手，提出一种“无自由边法”以实现船舶水密舱室的识别，并针对非水密舱室对算法进行改进与优化，以提高算法的适用性。1.1原理首先，根据某舱室的空间位置搜索该范围内的所有单元（Elem），统计这些单元的边（Edge）及其数量，其中数量等于1的一定是自由边（Freeedge）。根据该原则筛选出所有自由边，然后将其所属单元剔除出集合，接着重新统计边与数量，以此不断循环往复，直到没有自由边。最终，剩下的单元便形成封闭的三维空间，组成船舶舱室的搜索目标边界。现有研究的搜索方法大多属于“加法”，而本文提出的算法则采用的是“减法”思想。1.2实现过程舱室识别的流程如图1所示。步骤1根据船舶舱室边界的定义，获取该长方体空间内的所有单元Elem和及其Edge元素。步骤2 建立有限元模型的拓扑关系，并创建单元和边的数据结构储存拓扑关系。步骤3寻找唯一存在的Edge；若不存在，则结束。步骤4识别该Edge所属的单元，将其从表示封闭空间的单元数组中剔除，并更新边的数据储存Hash 表。步骤5将结果作为新的输人，回到步骤3，递归执行舱室识别过程。1.3建立拓扑关系船舶舱室是由一组单元组成的封闭空间区域，各个单元之间由单元边相互联系。为实现本文的算法，张志康，等：船舶有限元模型舱室识别方法No是否唯一EdgeYes剔除相应单元的Edges结束图1舱室识别流程Fig.1 Recognition process of ship compartment需建立的拓扑关系：单元与边的组合关系，边与边的重合关系。将单元边定义为从起点Startnode指向终点Endnode的有向边（Startnode，En d n o d e），其中起点的编号一定比终点的编号小。如图2 所示，单元Elem_1包含有向边(1,2)，(2,3)，(3,4)和(1,4)；单元Elem_2包含有向边(3,4)，（4,6)，(3,5)和(5,6)。它们存在共同边(3,4)，该边的数量为2，其余均为1。Node1-(1.4)(1,2)Elem_1Node2.Node3(2,3)图2 拓扑关系示意Fig.2Topological relation将Elem和其包含的有向边储存在二维数组中。另外以有向边的名称作为Key，数量作为Value，使用Hash表进行储存。1.4算法介绍“无自由边法”算法用于搜索三维模型中由单元构成的封闭空间。首先，确定搜索范围。舱室定义数据包括了舱室名称和x,和z方向的范围值。因此，对于单个舱室的搜索范围(Elem_i，它包含的单元坐标应满足式（1）式（3)。.41开始舱室定义搜索相应单元创建哈希表(Key:Edge,Value:Count)+Node4(3,4)/3.4)(3,5)(4,6)Elem2+Node5+Node6(5,6)42建立该范围内的拓扑关系，创建包含Elem和有向边的二维数组，及储存有向边和数量的Hash表。查询Hash表中Value等于1的Key，则该有向边即为自由边。通过二维数组获得该自由边所属的Elem，一方面将其从(Elem_i)内剔除，另一方面记录它包含的其他有向边，再次查询Hash表，各有向边对应的Value减一。把更新后的二维数组和Hash表作为新的输入，进行递归运算，直至无法找到自由边，则(Elem_i)内包含的单元即为该舱室的边界。运行过程如图3所示。(a)搜索范围模型(b）自由边及所属单元（循环1）(c)自由边及所属单元（循环2)图3“无自由边法”示意Fig.3No free edge algorithm2算法改进“无自由边法”本质是搜索有限元网格中的封闭空间。为了扩大算法的适用性，满足工程实际应用的需求，针对非封闭空间的识别问题，对本文的算法进行改进和优化。2.1舱室边界带有小孔在实际的舱室定义过程中，常常碰到有限元模型中水密舱室边界有小孔、不封闭的情况，此时类似的舱室搜索算法都会失效，而手动去寻找缺失的单元和修补也会非常耗时。导致此问题的原因有2 种：一种是由于建模人员在建模过程中的粗心，不小心遗漏了一些网格单元；另一种是考虑了水密舱室上的人孔。针对有限元模型中水密舱室边界上缺失多个且舰船科学技术Elem_i(x)=(x|x E(xXmin,Xmax),(1)Elem_i(y)=(yly E(Vmin,ymax),(2)Elem_i(z)=(zlz (zmin,Zmax)。(3)第45卷“独立”单元的情况，对本文提出的算法进行优化。大致的原理：不以有向边的数量来判断是否为自由边，而是以单元节点Node的数量作为依据，即搜索范围内Node数量2 的，其所属有向边即为自由边。基于该改进算法，便可准确识别此类带小孔的舱室。2.2货舱大开口对于散货船和集装箱船这类船型，货舱均带有大开口，对此现有的舱室搜索算法均不适用。文献5 中提出手动创建舱口盖进行封口的预处理方法，虽然能将此类大开口舱室转换为“封闭舱室”，但考虑到货舱的数量较多，因此处理起来还是过于繁琐。本文借鉴此预处理的思想，提出给此类货舱加一个虚拟“货舱盖”的改进算法。具体实现方法：在自由边的判断过程中，若有向边的2 个节点都在“货舱盖”的高度方向上，则该有向边必定不是自由边。该改进算法可由程序自动实现，避免了人工预处理的工作，有利于提高舱室识别效率。3应用实例FEMAP是一款先进的有限元分析软件，提供完全开放的二次开发接口，支持各种通用语言，非常适用于CAE的二次开发7 。基于该软件，采用VB.NET语言和WinForm客户端开发技术，实现本文提出的舱室识别算法。对某150 0 0 TEU集装箱船实船项目的全船有限元模型进行测试。全船有限元模型如图4所示，共计约62万个单元，需定义的舱室总计110 个。L图4测试集装箱船全船有限元模型Fig.4 Test global ship FEM of container vessel测试计算机的CPU为IntelXeonE5-26433.4GHz,内存容量为6 4GB，操作系统为Windows7（6 4位）。利用Napa软件输出总体专业的船舶舱室数据，并按照特定格式写人Excel文件，将其作为舱室必要数据的输人文件。导人成功后，一键运行程序即可快速、自动、准确地识别出绝大多数舱室，如首/尾尖舱、底/边压载水舱、LNG燃料舱等封闭舱室和集装箱货舱，如图5所示。第45卷(a)首尖舱(c)LNG燃料舱图5水密舱室边界Fig.5Watertight compartment boundary此外，对于机舱等个别非封闭舱室，需先采用手动填补开孔的方法进行预处理，然后再进行舱室边界自动搜索。机舱搜索结果如图6 所示。图6 机舱边界Fig.6Engine room boundary经多次测试，完成上述110 个舱室一键自动识别共需10 6 s，平均一个舱室耗时约1.0 s。经检查，结果无需修改的识别正确率约为9 0%，剩余的限于舱室特点和有限元模型质量，还需手动处理，耗时约半小时。总体而言，该程序识别率较高，所需时间较少，具有较强的工程应用价值。测试结果表明，对于具有不规则、复杂形状的舱室，本文的算法能很好地实现边界搜索。4结语结合实际设计和科研的需求，针对船舶有限元舱室自动识别问题进行研究，主要结论如下：1）提出的全新船舶有限元舱室识别算法“无自由边法”，可实现舱室空间的快速智能识别，编写程序对其有效性进行了验证；2）该算法本质是解决有限元三维封闭空间的识别问题，可拓展应用到类似的科研和开发项目中；张志康，等：船舶有限元模型舱室识别方法(b)底压载水舱(d)集装箱货舱433）通过实例验证，该算法实用可靠，除了水密舱室，对具有大开口或小开孔的非水密舱室也适用。程序操作简单，交互界面友好，可大幅提高设计人员的工作效率。参考文献：1高鹏，陈有芳.IACS协调共同结构规范的主要变化及影响.船海工程,2 0 14,43(6):2 530.GAO Peng,CHEN Youfang.Research on main differences andscantling impacts of IACS HCSR rulesJ.Ship and Ocean En-gineering,2014,43(6):25-30.2】单威俊,李锋,蒲海，等.船舶模型边界网格单元识别技术研究.计算机工程与设计，2 0 0 8(10):2 6 52-2 6 54.SHAN Weijun,LI Feng,PU Hai,et al.Research on boundaryelements of ship model identification techniquesJ.ComputerEngineering and Design,2008(10):2652-2654.3单威俊，许方，陈坤坤，等.海洋结构与船舶有限元模型空间识别方法.舰船科学技术,2 0 14,36(5):2 5-30+35.SHAN Weijun,XU Fang,CHEN Kunkun,et al.Research onspaces detection methods of offshore structure and ship FEMmodelJ.Ship Science and Technology,2014,36(5):25-30+35.4】单威俊,陈坤坤，许方，等.基于MSC.Patran 的舱室识别程序开发.船海工程，2 0 14,43(2)：7 7-8 0+8 4.SHAN Weijun,CHEN Kunkun,XU Fang,et al.Developmentof the cabin recognition program based on MSC.patranJJ.Shipand Ocean Engineering,2014,43(2):77-80+84.5章志兵,程登，孟凡冲,等.基于CAD模型的船舶舱室识别算法研究J.华中科技大学学报（自然科学版)，2 0 17，45(12)101-106.ZHANG Zhibing,CHENG Deng,MENG Fanchong,et al.Shipcompartment recognition algorithm based on CAD modelJ.Journal of Huazhong University of Science and Technology(Natural Science Edition),2017,45(12):101-106.6】章志兵，杨明润，王丽荣，等.基于CAE模型的船舶舱室自动识别高效算法.精密成形工程,2 0 2 0(4)：16 0-16 5.ZHANG Zhibing,YANG Mingrun,WANG Lirong,et al.An ef-ficient algorithm for automatic identification of ship cabinsbased on CAE modelJ.Journal of Netshape Forming Engi-neering,2020(4):160-165.7】张志康,司南,陈乐昆,等.FEMAP二次开发技术介绍与应用.自动化技术与应用,2 0 2 1,40(4)54-57.ZHANG Zhikang,SI Nan,CHEN Lekun,et al.Introduction andapplication of FEMAP secondary development technologyJ.Techniques of Automation and Applications,2021,40(4):