首页/文章/ 详情

HyperWorks在船舶甲板支撑结构设计中的应用(周维星)

5月前浏览283

周维星 徐峰 徐帅 朱眀罡

中国舰船研究设计中    湖北武汉   430064


摘要 :文章应用 HyperWorks 软件评估了液压折臂吊下甲板支撑构件的强度,并进一步应用OptiStruct 模块对液压折臂吊下的船舶甲板支撑结构进行了尺寸优化设计。优化过程中将甲板支撑构件尺寸参数作为变量,将中国船级社规范中规定的许用应力指标作为约束,将支撑结构总质量最小作为优化目标, 最终得到了满足规范要求的甲板支撑构件最优尺寸。

关键词:船舶结构, 尺寸优化, HyperWorks


1 引言

近年来,民船大型化,功能多样化已成为一种趋势。在船舶设计过程中,船体局部结构有限元强度计算任务较以往有所增加。设计部门需对船级社规范指定的结构,进行有限元强度计算,确保该结构应力符合规范要求,并编制好计算报告书送船级社审核。

在规范指定进行有限元强度校核的构件中,甲板设备支撑结构占了很大比重。其中,典型的甲板设备主要有锚机、起重机、吊杆、起重柱、系缆桩、导缆器和应急拖带装置等。按规范要求液压折臂吊属起重机吊杆一类,需进行支撑结构强度校核。本文以液压折臂吊为例,说明 HyperWorks 在船舶结构强度计算中的具体应用。

通常,在进行强度校核计算前,甲板支撑结构的构件尺寸已初步确定。构件尺寸的初始值是根据整条船的结构规范计算书得来的,这些尺寸主要是基于船舶种类、主尺度、骨架形式等全船性的参数根据规范计算出来的,没有考虑其上布置甲板设备带来的载荷。在以往的计算中,通常先校核构件初始尺寸的强度,如不满足规范满求,则逐步增大构件尺寸,直至满足规范要求。本文使用 HyperWorks 软件的优化功能,完成船舶甲板结构支撑构件的尺寸优化设计。该方法相对以往方法更加方便,在很大程度上提高了工作效率,通过将质量最小作为优化目标,得到的构件尺寸也更加合理。


2 有限元计算模型

在三维笛卡尔坐标系中建立液压折臂吊支撑结构舱段三维有限元计算模型。船长方向为X 轴,正方向由船尾指向船首; 船宽方向为 Y 轴,正方向由右舷指向左舷; 型深方向为 Z轴,正方向由基线指向甲板。规范要求计算模型采用局部立体结构模型,以液压折臂吊有效作用平面矩形的形心为中心, 向四周分别扩展至少一倍的该矩形相对应的长、宽距离,垂向应从液压折臂吊基座面扩展至甲板之下的第一个平台甲板或至少 D/4 处(D 为型深)。且应保证有限元计算模型范围延伸至结构主要支撑构件上。计算中将甲板纵桁、强横梁、强肋骨等强构件的腹板、主甲板板、舷侧板用壳单元离散;将甲板纵桁、强横梁、强肋骨等强构件的面板用梁单元离散;将甲板纵骨、舷侧纵骨及舱壁扶强材用梁单元离散。模型中位移单位为 mm,应力单位为 MPa。

船体材料主要采用屈服极限 235MPa  的 B  级船用结构钢,弹性模量 E=2 .06 × 105 N/mm2 ,泊松比 μ=0.3 。基于 HyperMesh 建立的液压折臂吊支撑结构舱段三维有限元计算模型见图 1。

image.png  image.png

 

                                         a)俯视图                               b)仰视图

                                 图 1 液压折臂吊支撑结构舱段有限元模型

3 初始尺寸强度校核

首先,针对液压折臂吊下甲板支撑结构初始尺寸进行强度校核。设计中将吊机布置在甲板强横梁和纵桁的交汇出,根据船体规范计算书算得的甲板强构件初始尺寸为:T 型材腹板9mm ,T 型材面板为 12x120 的扁钢。

3.1 载荷条件

根据规范要求,计算中主要考虑吊机的最大倾覆力矩、安全工作载荷和吊机自重三个方面的载荷。

1)吊机最大倾覆力矩:

“规范 ”中给出,倾覆力矩系指起重设备在安全载荷下作业,起重设备与船体结构连接处所计算得出的最大弯矩,本计算中折臂吊最大动弯矩为 308.2kN.m。

2)安全工作载荷:

安全工作载荷系指在任何规定的变幅长度范围内,起重机可起升的最大载荷。本吊机由设备商提供的安全工作载荷为 4t。根据“规范 ”,对于近海作业的起重设备,附加于起重设备自重的 150%安全工作载荷应进行校核,所以计算中安全工作载荷取6t。

3)液压折臂吊自重:

本船液压折臂吊自重为 2.38t。

3.2 边界条件

模型两端剖面各节点施加简支边界条件, 约束 x、y、z 方向线位移;模型下边界(包括主甲板以下 1.8m 处舷侧外板、肋骨以及支柱)施加简支边界条件,约束 x、y 、z 方向线位移;中纵剖面主甲板施加简支边界条件,约束 x、y、z 方向线位移。

3.3 分析结果

应用 HyperWorks 软件 RADIOSS 求解器进行计算,图2 给出了支撑结构初始尺寸下的应力云图。

image.png   image.png

                        a)Mises 应力云图                                    b)  剪切应力云图

image.png c)梁单元正应力云图

                                                图 2 初始尺寸强度校核结果

根据规范,构件的许用应力见表 1 。通过应力云图可见,吊机下甲板强横梁的腹板剪应力最大值为 95MPa,面板的正应力最大值为 224MPa,均不满足规范规定的许用应力值,需对吊机支撑结构进行加强。

image.png

其中:σS 为材料屈服应力。

4 优化分析


通过计算发现吊机载荷作用下,高应力区域集中出现在吊机下甲板交叉的强横梁和纵桁上,因此仅将 T 型材腹板厚度、T 型材面板宽度和面板厚度三个尺寸作为变量对支撑结构进行优化。考虑到工程实际,将各变量定义为离散型,如表 2。

image.png

在定义变量之后,指定响应作为优化约束条件和目标函数。针对本次研究,目标函数设置为甲板支撑结构质量最小; 约束条件为各应力分量不超过规范规定的许用值,表 3  给出了优化分析中约束的具体定义。

image.png

通过 OptiStruct  求解,最终优化结果为 T  型材腹板厚度 12mm ,T  型材面板厚度为17mm,T 型材面板宽度为 140mm。此时,液压折臂吊下甲板支撑结构应力云图如图 3。可见,优化尺寸后支撑结构的最大剪应力为 72.7MPa,最大正应力为 156.5Mpa 均满足了规范规定的许用值。此时,结构正应力值已经十分接近许用应力值,说明通过优化分析,我们得到了满足规范要求的最合理结构尺寸,使材料得到充分的利用, 节约了成本。

image.png  image.png

  

a)支撑结构剪切应力云图                         b)支撑结构正应力云图

                                  图 3 进行尺寸优化后结构的应力云图

5  结论

本文以液压折臂吊下甲板支撑结构为优化对象,采用 HyperWorks 软件进行强度评估并完成了优化分析,获得了满足船舶规范要求的结构最优尺寸。相对以往方法,该方法更加科学方便,获得的结构尺寸也更加合理。 除了船体结构强度方面,在船体振动、冲击响应等方向,还可以将该方法进行进一步推广应用,这也是本文作者下一步的研究内容。


HyperWorks Application in the Design of the Supporting Structure of the Ship Deck

Zhou Weixing       Xu Feng        Xu Shuai       Zhu Minggang

Abstract: This paper shows the application of HyperWorks software in evaluating the strength  of  the  supporting  structure  of  hydraulic  folding  arm,  and  the  application  of OptiStruct module in the optimization of it. Taking deck support structure size as design variable, allowed stress as constraint, minimize the total mass. Finally get an optimal size.

Keywords: Ship Structure       Size Optimization        HyperWorks


免责声明:


本页面/内容部分素材来源于互联网公 开 信 息,旨在传递更多信息,不代表本平台立场。

版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。

本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险。

船舶Star-CCM+
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2026-04-30
最近编辑:5月前
仿真支持爱好者
在仿真的路上越走越远
获赞 333粉丝 27文章 318课程 0
点赞
收藏
作者推荐

程序员专用 DSP 算法:数学原理与源码实现

程序员专用DSP算法:数学原理与源码实现((Dsp)AlgorithmsForProgrammers)摘要:《(Dsp)AlgorithmsForProgrammers》是一本面向程序员的DSP算法与数学原理实战手册,聚焦数字信号处理核心算法的数学原理与代码实现。全书系统讲解离散傅里叶变换、基2/基4/分裂基FFT、卷积、Z变换、哈特利变换、数论变换、沃尔什变换、哈尔变换等核心算法,给出递归与非递归实现、三角函数计算优化、内存访问改进及代码自动生成方法。同时涵盖位操作、排列、排序、组合算法、高精度运算等底层技术,提供可直接复用的伪代码与C++源码,兼顾理论推导与工程效率,适合DSP开发、算法实现与高性能计算参考。AlgorithmsforprogrammersideasandsourcecodeThisdocumentisworkinprogress:readthe”importantremarks”nearthebeginningJ¨orgArndtarndt@jjj.deThisdocument1wasLATEX’datSeptember26,2002Someimportantremarksaboutthisdocument.Thisdraftisintendedtoturnintoabookaboutselectedalgorithms.Theaudienceinmindarepro-grammerswhoareinterestedinthetreatedalgorithmsandactuallywanttohave/createworkingandreasonablyoptimizedcode.Theprintablefullversionwillalwaysstayonlineforfreedownload.ItisplannedtoalsomakepartsoftheTEXsources(plusthescriptsusedforautomation)available.RightnowafewfilesoftheTEXsourcesandallextractedpseudo-codesnippets1areonline.TheC++-sourcesareonlineaspartofFXTorhfloat(arithmeticalalgorithms).ThequalityandspeedofdevelopmentdoesdependonthefeedbackthatIreceivefromyou.Yourcriticismconcerninglanguage,style,correctness,omissions,technicalitiesandeventhegoalssethereisverywelcome.Thankstothose2whohelpedtoimprovethisdocumentsofar!Thanksalsotothepeoplewhosharetheirideas(orsourcecode)onthenet.Itrytogiveduereferencestooriginalsources/authorswhereverIcan.However,IaminnowayanexpertforhistoryofalgorithmsandIprettysurewillneverbeone.Soifyoufeelthatareferenceismissingsomewhere,letmeknow.Newchapters/sectionsappearassoonastheycontainanythinguseful,sometimesjustlistingsorremarksoutliningwhatistoappearthere.A”TBD:somethingtobedone”isaremindertomyselftofillinsomethingthatismissingorwouldbenicetohave.ThestylevariesfromchaptertochapterwhichIdonotconsiderbadperse:whilesometopics(e.g.FFTs)needaclearandexplicitintroductionothers(e.g.thebitwizardrychapter)seemtobebestpresentedbybasicallyshowingthecodewithjustafewcomments.Stillotherparts(e.g.sorting)arepresentedelsewhereextremelywellsoIwillintroducethebasicideasonlyveryshortlyandsupplysome(hopefully)usefulcode.Sprachewillpartlygoaway:using/includingtheactualcodefromFXTwillbebeneficialtoboththisdocumentandFXTitself.Thegoalistoautomaticallyincludethefunctionsreferenced.Clearly,thiswilldrasticallyreducethechanceoferrorsintheshowncode(andatthesametimedrasticallyreducetheworkloadforme).InitiallyIplannedtowriteaninterpreterforSprache,itjustneverhappened.AtthesametimeFXTwillbebetterdocumentedwhichitreallyneeds.AsaconsequenceSprachewillonlybeusedwhenthereisaclearadvantagetodoso,mainlywhenthecorrespondingC++doesnotappeartobeselfexplanatory.LargerpiecesofcodewillbepresentedinC++.AtinystarteraboutC++(somegoodreasonsinfavorofC++andsomeoftheverybasicsofclasses/overloading/templates)willbeincluded.Cprogrammersdonotneedtobeshockedbytheζ++’:onlyanratherminimalsetoftheC++featuresisused.Thetheorem-likeenvironmentforthecodesshallcompletelygoaway.Itleadstoduplicationofstate-ments,especiallywithnon-pseudocode(runningtext,descriptionintheenvironmentandcommentsatthebeginoftheactualcode).Enjoyreading!ListofimportantSymbols沦xrealpartofx$ximaginarypartofxx∗complexconjugateofxaasequence,e.g.{a0,a1,...,an−1},theindexalwaysstartswithzero.atransformed(e.g.Fouriertransformed)sequencememphasizethatthesequencestotheleftandrightarealloflengthm=F[a](=c)(discrete)Fouriertransform(FT)ofa,ck=Σx(n)0(1)axzxkwherez=e±2πi/nF−1[a]inverse(discrete)Fouriertransform(IFT)ofa,F−1[a]k=Σx(n)0(1)axz−xkSkaasequencecwithelementscx:=axe±k2πix/nH[a]discreteHartleytransform(HT)ofaasequencereversedaroundelementwithindexn/2aSthesymmetricpartofasequence:aS:=a+aaAtheantisymmetricpartofasequence:aA:=a−aZ[a]discretez-transform(ZT)ofaWv[a]discreteweightedtransformofa,weight(sequence)vW1[a]inversediscreteweightedtransformofa,weightva④bcyclic(orcircular)convolutionofsequenceawithsequenceba④acbacyclic(orlinear)convolutionofsequenceawithsequenceba④−bnegacyclic(orskewcircular)convolutionofsequenceawithsequenceba④{v}bweightedconvolutionofsequenceawithsequenceb,weightva④⊕bdyadicconvolutionofsequenceawithsequencebn\NndividesNn⊥mgcd(n,m)=1a(j%m)sequenceconsistingoftheelementsofawithindicesk:k≡jmodme.g.a(even),a(odd)a(0%2),a(1%2)a(j/m)sequenceconsistingoftheelementsofawithindicesk:j·n/m≤k<(j+1)·n/me.g.a(left),a(right)a(0/2),a(1/2)Chapter1TheFouriertransform1.1ThediscreteFouriertransformThediscreteFouriertransform(DFTorsimplyFT)ofacomplexsequenceaoflengthnisdefinedaszisann-throotofunity:zn=1.Backtransform(orinversediscreteFouriertransformIDFTorsimplyIFT)isthenToseethis,considerelementyoftheIFToftheFTofa:AsΣk(z①__y)k=nforx=yandzeroelse(becausezisann-throotofunity).ThereforethewholeexpressionisequaltowhereHerewewillcalltheFTwiththeplusintheexponenttheforwardtransform.Thechoiceisactuallyarbitrary1.CHAPTER1.THEFOURIERTRANSFORMThenormalizationfactorinfrontoftheFTsumsissometimesreplacedbyasingleinfrontoftheinverseFTsumwhichisoftenconvenientincomputation.Then,ofcourse,Parseval’sequationhastobemodifiedaccordingly.AstraightforwardimplementationofthediscreteFouriertransform,i.e.thecomputationofnsumseachoflengthnrequires~n2operations:[FXT:slowftinslow/slowft.cc]ismustbe+1(forwardtransform)or__1(backwardtransform),SinCos(x)returnsaComplex(cos(x),sin(x)).AfastFouriertransform(FFT)algorithmisanalgorithmthatimprovestheoperationcounttopropor-tionalnΣk__1),wheren=p1p2···pmisafactorizationofn.Incaseofapowern=pmthevaluecomputeston(p__1)logp(n).Inthespecialcasep=2evenn/2log2(n)(complex)multiplicationssuffice.ThereareseveraldifferentFFTalgorithmswithmanyvariants.1.2SymmetriesoftheFouriertransformAbitofnotationturnsouttobeuseful:Letabethesequencea(lengthn)reversedaroundelementwithindexn/2:LetaS,aAbethesymmetric,antisymmetricpartofthesequencea,respectively(Theelementswithindices0andn/2ofaAarezero).Nowleta∈R(meaningthateachelementofais∈R),thenTheFTofarealsymmetricsequenceisrealandsymmetricandtheFTofarealantisymmetricsequenceispurelyimaginaryandantisymmetric.TherebytheFTofageneralrealsequenceisthecomplexconjugateofitsreversed:Similarly,forapurelyimaginarysequenceb∈iR:TheFTofacomplexsymmetric/antisymmetricsequenceissymmetric/antisymmetric,respectively.1.3Radix2FFTalgorithms1.3.1AlittlebitofnotationAlwaysassumeaisalength-nsequence(napoweroftwo)inwhatfollows:Leta(even),a(odd)denotethe(length-n/2)subsequencesofthoseelementsofathathaveevenoroddindices,respectively.Leta(left)denotethesubsequenceofthoseelementsofathathaveindices0...n/2__1.Similarly,a(right)forindicesn/2n__1.LetSkadenotethesequencewithelementsaxe§k2πix/nwherenisthelengthofthesequenceaandthesignisthatofthetransform.ThesymbolSshallsuggestashiftoperator.InthenexttwosectionsonlyS1/2willappear.S0istheidentityoperator.1.3.2Decimationintime(DIT)FFTThefollowingobservationisthekeytothedecimationintime(DIT)FFT2algorithm:Forneventhek-thelementoftheFouriertransformiswherez=e§i2π/nandk∈{0,1,,n__1}.Thelastidentitytellsushowtocomputethek-thelementofthelength-nFouriertransformfromthelength-n/2Fouriertransformsoftheevenandoddindexedsubsequences.Toactuallyrewritethelength-nFTintermsoflength-n/2FTsonehastodistinguishthecases0≤k<n/2andn/2≤k<n,thereforewerewritek∈{0,1,2,...,n__1}ask=j+δwherej∈{0,1,...,n/2__1},δ∈{0,1}.Notingthatz2isjusttherootofunitythatappearsinalength-n/2FTonecanrewritethelasttwoequationsastheIdea1.1(FFTradix2DITstep)Radix2decimationintimestepfortheFFT:(Hereitissilentlyassumedthat’+’or’__’betweentwosequencesdenoteselementwiseadditionorsubtraction.)Thelength-ntransformhasbeenreplacedbytwotransformsoflengthn/2.Ifnisapowerof2thisschemecanbeappliedrecursivelyuntillength-onetransforms(identityoperation)arereached.Therebytheoperationcountisimprovedtoproportionaln·log2(n):Therearelog2(n)splittingsteps,theworkineachstepisproportionalton.Code1.1(recursiveradix2DITFFT)Pseudocodeforarecursiveprocedureofthe(radix2)DITFFTalgorithm,ismustbe+1(forwardtransform)or-1(backwardtransform):Thedatalengthnmustbeapowerof2.Theresultisinx[].Notethatnormalization(i.e.multiplicationofeachelementofx[]by1/√n)isnotincludedhere.[FXT:recursivedit2fftinslow/recfft2.cc]TheprocedureusesthesubroutineCode1.2(Fouriershift)Foreachelementinc[0..n-1]replacec[k]byc[k]timesev2πik/n.UsedwithU=§1/2fortheFouriertransform.TherecursiveFFT-procedureinvolvesnlog2(n)functioncalls,whichcanbeavoidedbyrewritingitinanon-recursiveway.Onecanevendoalloperationsinplace,notemporaryworkspaceisneededatall.Thepriceisthenecessityofanadditionaldatareordering:Theprocedurerevbin_permute(a[],n)rearrangesthearraya[]inawaythateachelementaxisswappedwitha˜x,where˜xisobtainedfromxbyreversingitsbinarydigits.Thisisdiscussedinsection8.1.Code1.3(radix2DITFFT,localized)Pseudocodeforanon-recursiveprocedureofthe(radix2)DITalgorithm,ismustbe-1or+1:[FXT:dit2-fft-localizedinfft/fftdit2.cc]Thisversionofanon-recursiveFFTprocedurealreadyavoidsthecallingoverheadanditworksinplace.Itworksasgiven,butisabitwasteful.The(expensive!)computatione:=exp(is*2*PI*I*j/m)isdonen/2·log2(n)times.Toreducethenumberoftrigonometriccomputations,onecansimplyswapthetwoinnerloops,leadingtothefirstζrealworld’FFTprocedurepresentedhere:Code1.4(radix2DITFFT)Pseudocodeforanon-recursiveprocedureofthe(radix2)DITalgo-rithm,ismustbe-1or+1:procedurefft_dit2(a[],ldn,is)//complexa[0..2**ldn-1]input,resultSwappingthetwoinnerloopsreducesthenumberoftrigonometric(exp())computationstonbutleadstoafeaturethatmanyFFTimplementationsshare:Memoryaccessishighlynonlocal.Foreachrecursionstage(valueofldm)thearrayistraversedmhtimeswithn/maccessesinstridesofmh.Asmhisapowerof2thiscan(oncomputersthatusememorycache)haveaverynegativeperformanceimpactforlargevaluesofn.OnacomputerwheretheCPUclock(366MHz,AMDK6/2)is5.5timesfasterthanthememoryclock(66MHz,EDO-RAM)IfoundthatindeedforsmallnthelocalizedFFTisslowerbyafactorofabout0.66,butforlargenthesameratioisinfavouroftheζnaive’procedure!Itisagoodideatoextracttheldm==1stageoftheoutermostloop,thisavoidscomplexmultiplicationswiththetrivialfactors1+0i:Replace更多内容见附件免责声明:本页面/内容部分素材来源于互联网公开信息,旨在传递更多信息,不代表本平台立场。版权归原作者或机构所有,如涉及侵权,请通过平台联系我们,我们将在核实后第一时间处理。本平台对转载内容的真实性、准确性不作任何保证,用户需自行判断并承担使用风险

有附件
未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习计划 福利任务
下载APP
联系我们
帮助与反馈