还剩11页未读,继续阅读
本资源只提供10页预览,全部文档请下载后查看!喜欢就下载吧,查找使用更方便
文本内容:
第18卷第10期船舶力学Vol.18No.102014年10月Journal ofShip MechanicsOct.2014文章编号1007-7294
(2014)10-1228-09激光焊接夹层甲板板格强度计算的子模型方法朱扬,程远胜,刘均(华中科技大学船舶与海洋工程学院,武汉430074)摘要提出了基于有限元软件ANSYS的激光焊接钢质夹层甲板板格结构强度计算的子模型方法分别对考虑激光焊接焊缝缺陷的I型夹层板格结构的壳单元计算模型和体单元子模型进行强度分析,并与全部体单元模型的芯层与上下面板连接处及面板中部处应力分布计算结果进行对比,验证壳单元计算模型和子模型方法用于计算夹层甲板板格强度的正确性计算结果表明,对强度特征关注区域,可建立多个体单元子模型,确定子模型边界影响区域范围,从而可较为准确地评估夹层甲板板格结构强度特性,包括焊缝处应力分布壳单元计算模型可获得较为精确的板格变形值,但无法考虑激光焊接焊缝缺陷,获得的焊缝处最大应力值明显偏小关键词激光焊接夹层板格;子模型;有限元中图分类号U
661.43文献标识码A doi
10.3969/j.issn.1007-
7294.
2014.
10.009Sub-model methodfor strengthcalculation of a laser-welded steelsandwich panelstructure ZHUYang,CHENG Yuan-sheng,LIU JunCollege ofNaval Architectureand OceanEngineering,Huazhong Universityof Scienceand Technology,Wuhan430074AbstractThis paperproposes sub-model methodfor calculating the strengthof alaser-welded steelsand-wich panelstructure byusing commercialcode ANSYS.The strengthcomputations of the sandwich panel us-ingthe all-shell modeland sub-model areconducted respectivelyand the stress distributionsnections between the core and upperbottom sheetand onesat centralpart bythose twocomputati models are comparedrespectively withthose obtained from theall-solid modelto verifythe theall-shell modeland the sub model.The resultsshow thatmultiple submodels arenecessary fordeter-mining the accurate zonesin whichthe stressesfrom thedifferent submodelsaresame orclosely similar.In thisway,theaccuratestresses ofinterested structureparts includingthe laser-welded defectscan beob-tained.The resultsalso verifythat thedeformation of the panelobtainedfrom theall-shell modelis accept-able,however,the stressesat the laser-welded defectsfrom theall-shell modelare muchles obtainedfromtheall-solid model.Key words:laser-welded steelsandwichpanel;sub-model;finite elementanalysis1引言金属夹层结构采用以激光焊接为基础的制造技术,能够显著提高生产效率和降低能耗,具有高比收稿日期2013-11-15基金项目国家自然科学基金
(51079058)作者简介朱扬(1987-),男,硕士研究生,E-mail:zhuyang334@
163.com;程远胜(1962-),男,教授,博士生导师,通讯作者,E-mail:yscheng@mail.hust.edu.cn;刘均(1981-),男,博士,讲师,E-mail:hustlj@mail.hust.edu.cn第10期朱扬等激光焊接夹层甲板板格强度计算的…1229刚度和高比强度的特点,同时,夹层结构应用在船舶甲板上能有效降低甲板结构自身高度,改善舱室空间的利用率
[1]此外,在耐撞击,抗冲击以及防火、隔热等方面也具有优良的性能,在船舶结构制造中有着广阔的应用前景
[2]但是在制备过程中,由于激光焊接焊缝缺陷的存在,导致在焊缝处产生明显的应力集中
[3],如图1所示目前,关于夹层板结构强度特性,已有多位学者展开研究刘均等
[4]在经典夹层板理论模型的基础上,运用能量原理推导出了方形蜂窝夹层板在面内载荷作用下的平衡微分方程,并得到了方形蜂窝夹层板结构总体屈曲的临界载荷图1激光焊接焊缝缺陷处应力分布示意图
[3]值何力和程远胜
[5]提出了基于整体位移场与局Fig.1The stressdistribution in thelaser-welded defects部位移场叠加的夹层板结构强度理论分析新方法,对波纹型、蜂窝型和X型芯层夹层板的强度、稳定性和自由振动特性进行分析计算以上几位学者在数值仿真研究中,均未考虑激光焊接焊缝缺陷对结构强度特性的影响Niklas
[6]采用平面应变模型,考虑激光焊接焊缝细节,分析了两种典型夹层板连接构件在拉压载荷作用下强度的变化规律,得出两种连接形式的最优设计方案,并进行了对比分析准确计算激光焊接焊缝缺陷处的应力,对正确评估夹层甲板强度特性有重要意义在数值仿真计算中采用子模型方法,则可以在有限的计算资源下有效模拟激光焊接焊缝缺陷,从而对夹层甲板板架结构强度进行精细分析Wolfgang和Paetzold
[7]在船舶结构疲劳强度实尺寸实验研究中,应用子模型方法对船舶结构的疲劳强度进行数值仿真计算,在子模型建模中考虑了焊接细节,从而获得与实验结果吻合较好的数值结果赵尚辉
[8]则利用子模型方法对救生船中央月池开口部位的应力集中进行数值仿真分析,计算结果与实验数据吻合较好王虎等
[9]则基于有限元数值仿真,分别采用壳体连接技术和子模型法,对I型金属夹层结构的两种典型连接构件进行了强度分析,提出了I型金属夹层结构连接构件强度数值计算方法本文提出了激光焊接夹层甲板板格结构强度数值仿真分析的有限元子模型方法,考虑了激光焊接焊缝缺陷对结构强度特性的影响,成功地克服了结构单元规模过大而无法计算的问题,提高了计算效率,为准确评估夹层甲板板格结构强度提供了数值仿真方法2“I”型夹层甲板板格结构强度计算的子模型法
2.1结构计算的子模型方法在夹层甲板板格强度数值仿真计算过程中,采用壳单元建模,无法模拟出激光焊接焊缝的缺陷,计算结果不足以作为强度校核的最后依据采用实体单元模拟焊缝缺陷,由于焊缝缺陷几何尺度较小,受限于计算机硬件及运算时间而无法计算故利用有限元高级分析技术子模型法,为关注的焊接缺陷结构区域建立子模型,作精细分析,不仅使计算量在计算机硬件能够承受的范围之内,而且能获得较准确的激光焊接焊缝缺陷处应力分布子模型方法基于圣维南原理提出
[10]应力集中只对其附近的区域产生影响,当子模型的边界离应力集中区较远时,就可通过子模型得到比较准确的结果实际计算中首先对结构对象采用壳单元建模,称之为主模型,并对其进行有限元计算分析再把关注的局部结构从主模型中取出来,作为一个独立的模型,采用壳单元或者实体单元重新建模,称之为子模型,单独对其精细划分网格并以主模型在切割边界上的计算位移值作为子模型的边界条件再次分析计算对于夹层甲板板格结构子模型,切割1230船舶力学第18卷第10期边界要完全避开应力集中区域很难做到,但是切割边界的应力集中仅在一定范围内影响子模型计算结果本文利用大型有限元软件ANSYS的子模型技术,对I型夹层板结构算例进行强度特性分析,并与全部体单元模型计算结果进行对比,讨论切割边界经过应力集中区域对子模型计算结果的影响,验证子模型方法用于计算夹层甲板板格结构强度的正确性
2.2I型夹层甲板板格结构强度分析本文分别在计算模型中部及x=0(如图2所示)边界处不同几何位置截取子模型,按照生成并分析较粗糙的主模型,生成子模型,提供切割边界插值,分析子模型,验证计算结果正确性的步骤,对I型夹层甲板板格结构算例强度特性进行分析
2.
2.1I型夹层甲板板格结构计算模型I型夹层甲板板格结构设计变量共5个参数上面板板厚tt、下面板板厚tb,夹层板格高度hc,芯层腹板板厚tc,芯层腹板间距b(芯层腹板数量n)如图2所示算例几何参数如表1所示图2I型夹层甲板板架结构参数表征Fig.2The parametercharacterization of the I-core sandwichstructure panel表1I型夹层甲板板格结构几何参数Tab.1The parameterof theI-core sandwichstructure plate长(mm)宽(mm)tc(mm)tt(mm)tb(mm)n hc(mm)
1000.
0720.
01.
002.
002.
00460.00材料弹性模量为
2.06E+11Pa,密度为7850kg/m3,泊松比为
0.3计算强度时,板格结构的边界条件处理为固支边界,在上面板上施加10kPa均布压力
2.
2.2板格强度计算的主模型取夹层结构相邻两腹板及之间面板作为胞元,计算模型沿胞元高度方向上划分4份单元,胞元上下面板各划分8份单元,整个计算模型共有壳单元6120个计算模型如图3,4所示图3主模型单元划分图4主模型强度计算加载情况Fig.3Element meshof the main model Fig.4The loadingapplied tothemain model第10期朱扬等激光焊接夹层甲板板格强度计算的…1231I型夹层甲板板格结构在10kPa载荷作用下的整体变形云图如图5所示,最大变形为
0.283mm;应力云图如图6所示,最大Von Mises应力
35.1MPa,位于腹板和上面板连接位置图5主模型在均布载荷作用下变形云图图6主模型在均布载荷作用下应力云图Fig.5Deformation distributionof themain model under Fig.6Stress distributionof themain modelunder the uniform out-of-plane loadthe uniform out-of-plane load
2.
2.3板格强度计算的子模型对于子模型,必须使用与主模型中一致的单元实常数和材料性能参数子模型和主模型采用完全一致的空间坐标系,并且子模型建模位置就是该局部区域在主模型中的空间位置,前后两者必须完全一致为了评估I型夹层板格结构强度特性,子模型应当能够包括x=0处边界中部区域,芯层腹板和面板连接区域,上下面板中部区域故建立子模型时选择I型夹层板格结构中部区域板条,保留中部两芯层腹板及腹板间面板结构,子模型位置如图7所示分别选取腹板间距1/4,1/2,3/4处为切割边界,在主模型相同几何位置,利用solid45单元分别建立子模型(下称板条子模型),在芯层腹板与上下面板连接处利用小圆孔模拟激光焊接缺陷,根据相关文献试验研究结果
[11],圆孔半径取芯层腹板厚度的1/10-1/5,本计算模型取小孔半径
0.1mm以腹板间距1/2处切割边界板条子模型为例进行分析,其余两个子模型可做类似处理为保证计算精度,子模型在截面上沿圆孔周向划分12份单元,如图8所示,沿胞元高度方向上划分60份单元,胞元上下面板各划分80份单元,上下面板及芯层厚度方向分别划分4份单元,并沿腹板长度方向掠扫,生成六面体网格,单元共计74万单元,如图9所示子距离模分型别位为腹板间距的1/4,1/2,3/4置图7面板中部应力计算子模型位置示意图图8子模型激光焊接缺陷处单元划分Fig.7Location of thesub-model formiddle partFig.8Element meshof thesub-model in thestresscomputation laser-welded defects在主模型固支边界相同位置,分别为子模型施加固支边界条件,并在子模型上表面施加10kPa均布压力根据主模型切割边界处节点的角位移和平动位移,为子模型提供切割边界插值,如图10所示1232船舶力学第18卷第10期图9子模型单元划分情况图图10子模型加载及边界条件示意Fig.9Element meshof thesub-model Fig.10Loading andboundary conditionsof thesub-model为研究子模型切割边界处应力集中对计算结果的影响,同时也为减小单元规模,提高计算效率,在研究夹层板格结构中部应力特性时,在板条子模型的基础上,选取其中部区域,在距固支边界分别1/3,1/4,1/5板格长度处为切割边界,经过面板和腹板连接处应力集中区域建立子模型(下称板条中部子模型),子模型位置如图11所示在研究夹层板格结构固支边界处应力特性时,选取板条左端区域,距离x=0边界处1/3,1/4,1/5板格长度处为切割边界,建立子模型(下称板条边界子模型),子模型位置如图12所示计算分析方法与板条子模型类似子应距力离模集分型中别位区为板格总长的1/3,1/4,1/5置域应子距力离模集分型中别位区为板格总长的1/3,1/4,1/5域置图11板条中部子模型位置示意图图12板条边界子模型位置示意图Fig.11Location of the centralpart sub-model inFig.12Location of the sidepart sub-model thecenter of the panelalong the core of the panel图13全体单元模型单元划分情况图14全体单元模型强度计算加载情况Fig.13Element meshof theall-solid modelFig.14Loading conditionsof theall-solid model
2.
2.4子模型方法计算结果验证及分析为了验证子模型方法的正确性,全部应用solid45单元建立计算模型,焊接缺陷处圆孔半径取为第10期朱扬等激光焊接夹层甲板板格强度计算的…
12330.1mm,与子模型一致材料参数,施加载荷及边界条件与主模型一致全部体单元网格密度与子模型保持一致,模型在截面上沿圆孔周向划分12份单元,并沿腹板长度方向掠扫,生成六面体网格完整模型共计149万单元,如图13,14所示对比全部体单元建立的整体模型和壳单元建立的主模型变形云图,如图15,16所示,全体模型和图15全体模型在均布载荷作用下变形云图图16主模型在均布载荷作用下应力云图板格轴向长度m板格轴向长度m Fig.15Deformation distributionof theall-solid modelFig.16Deformation distributionof themainmodelunder under the uniform out-of-plane loadthe uniform out-of-plane load合成应力值Pa合成应力值Pa图17板条子模型腹板与上面板连接处应力分布图18板条子模型腹板与下面板连接处应力分布板格轴向长度板格轴向长度m mFig.17Stress distributionin theconnected Fig.18Stress distributionin theconnected regionbe-tween the core and upper sheet at regionbe-tweenthe core and bottom sheet at the middle part of thepane the middle part of the panel图19板条子模型上面板中部处应力分布图20板条子模型下面板中部处应力分布Fig.19Stress distributionin the middle regionof theFig.20Stress distributionin the middle regionofthupper sheet at the middle part ofthe panel lowersheet atthemiddlepartofthe panel合成应力值Pa合成应力值Pa1234船舶力学第18卷第10期主模型变形云图吻合较好,主模型切割边界位置上的位移结果可以作为子模型的计算条件由于应力关注区域限于两芯层腹板及之间的面板区域,对主模型,全体模型和板条子模型,选取芯层腹板和上下面板连接处,以及上下面板中部,沿芯层长度方向每隔10mm选择测点,分别提取Von Mises应力,并提取上下面板中部位置位移,对比结果如图17-20所示由图21-24对比结果分析可知,板条子模型两腹板及之间面板区域应力分布和位移分布几乎不受切割边界处面板宽度的影响,各子模型与全体单元模型及主模型吻合较好,主模型应力偏小,在计算中建议选取1/2腹板间距板条子模型进行计算分析板格轴向长度m板格轴向长度m对板条中部子模型和板条边界子模型做进一步分析,同样选取芯层腹板和上下面板连接处,沿芯层长度方向每隔10mm选择测点,分别提取Von Mises应力,对比结果如图19-22所示图21板条中部子模型腹板与上面板连接处应力分布Fig.21Stress distributionin thecentral sub-model be-合成应力值Pa coreand uppersheet atthemiddle图22板条中部子模型腹板与下面板连接处应力分布Fig.22Stress distributionin thecentral sub-model be-板格轴向长度m tweenthecoreand bottomsheetatthemiddle板格轴向长度m partofthepanel partofthepanel图23板条边界子模型腹板与上面板连接处应力分布图24板条边界子模型腹板与下面板连接处应力分布Fig.23Stress distributionintheside sub-model betweenFig.24Stress distributionintheside sub-model betweenthecoreanduppersheetatthemiddlepartofthecoreandbottomsheetatthemiddlepartofthepanelthepanel由图21-24对比结果可知;由于切割边界直接经过腹板与面板连接应力集中位置,板条中部子模型及板条边界子模型在靠近切割边界处应力明显偏大,但是应力受影响区域极为有限,约为板格总长度的1/20在影响区域之外,与全体单元模型相比,应力计算结果最大误差
3.8%以1/4切割边界板条中部子模型为例,除去子模型两端板格总长度的1/20区域,与全部体单元模型相应区域对比位移云图及应力云图,如图25-28所示子模型与全部体单元云图较为一致经对比,三类模型最大位移值相当子模型与全部体单元模型应力分布较为一致,而主模型由于未考虑激光焊接焊缝缺陷,在腹板与两面板连接处应力值明显偏小合成应力值Pa合成应力值Pa第10期朱扬等激光焊接夹层甲板板格强度计算的…1235图25子模型在均布载荷作用下变形云图图26子模型在均布载荷作用下应力云图Fig.25Deformation distributionofthesub-modelFig.26Stress distributionofthesub-modelunderunder theuniformout-of-plane loadtheuniformout-of-plane load图27全体模型在均布载荷作用下变形云图图28全体模型在均布载荷作用下应力云图Fig.27Deformation distributionoftheall-solid modelFig.28Stress distributionoftheall-solid modelunder theuniformout-of-plane loadundertheuniformout-of-plane load在应用子模型方法计算夹层甲板板格强度时,由于全体单元模型单元规模过大,故可以根据关注区域选取板条子模型,板条中部子模型或者板条边界子模型进行计算,为了确定切割边界经过应力集中区域对计算结果的影响,可以在关注区域建立两个或多个子模型进行计算,对比其腹板和上下面板连接处应力分布情况,受影响区域应力较正常结果明显偏大,但是除去受影响区域,两个子模型相应位置应力分布较为一致,当关注区域处在应力一致区域范围中时,计算结果与全体单元模型吻合,可以较为准确地评估夹层甲板板格结构的强度特性对于芯层形式不同,但是芯层仅沿单一方向分布的夹层甲板板格结构,也可采用类似的子模型方法评估其强度特性4结论本文提出了夹层甲板板格结构强度计算的子模型方法,通过一算例,在不同位置截取体单元子模型,与全部体单元模型计算结果进行对比,研究了切割边界经过应力集中区域对计算结果的影响,验证子模型方法用于计算夹层甲板板格结构强度的正确性,获得如下结论
(1)切割边界位置对计算结果的影响区域极为有限,对强度特征关注区域建立多个子模型,确定影响区域范围,从而可较为准确地评估夹层甲板板格结构强度特性,同时提高计算效率
(2)夹层甲板板格结构强度壳单元计算模型可以获得较为精确的变形值,但是焊缝处最大应力值明显偏小,采用体单元子模型可以模拟激光焊接焊缝缺陷,并较为准确地计算焊缝处应力分布结果1236船舶力学第18卷第10期参考文献
[1]Roland F,Reinert T.Laser welded sandwich panelsfor theshipbuilding industry[C]//Lightweight Construction-L velopments.London,2000:1-
12.
[2]Denney P.Corrugated coreLASCOR structuresin shipbuildingand otherapplications[C]//Presentation Lightweightwich Seminar.Helsinki,Finland,
1998.
[3]Dariusz B,Szala J.Tests oflocal strainsin steellaser-weldedsandwichstructure[J].Polish MaritimeResear31-
36.
[4]Liu J,Cheng YS,Li RF,et al.A semi-analytical methodfor bending,buckling,and freevibration analysespanels withsquare-honeycomb cores[J].International Journalof StructuralStability andDynamics,2010,101:127-He L,Cheng YS,Liu J.Precise bendingstress analysis of corrugated-core,honeycomb-coreandX-core sandwichp els[J].Composite Structures,2012,945:1656-
1668.
[6]Kozak J.Forecasting offatigue lifeof laserwelded joints[J].Zagadnienia EksploatacjiMaszyn,2007,1491:8
[7]Wolfgang F,Paetzold H.Full-scale fatiguetests ofship structuresto validatethe S-N approachesfor fatiguesessment[J].Marine Structures,201023:115-
130.
[8]Zhao Shanghui,Ren Hong,Li Fujun,et al.Stress concentrationanalysisofmoon-pool inrescue ship[J].Journa Mechanics,2009,132:278-
284.in Chinese
[9]王虎,程远胜,刘均.I型金属夹层结构连接构件强度数值计算方法[J].中国舰船研究,2012,73:51-
56.
[10]刘涛,杨凤鹏,李贵敏,等.精通ANSYS[M].北京:清华大学出版社,
2002.
[11]Ryszard P.Strength assessmentofaversion ofjoint ofsandwich panels[J].Polish MaritimeResearch,2006s1。