发明名称 基于有限元分析的船体结构焊脚高度的确定方法
摘要 基于有限元分析的船体结构焊脚高度的确定方法,包括如下五个步骤:步骤1,构建船舶结构角焊缝应力的计算模型,提出一种从舱段有限元模型中提取焊缝应力的方法;步骤2,确定舱段有限元模型的板单元剪应力与焊缝平均剪应力之间的关系;步骤3,提出船体结构焊缝强度的应力标准;步骤4,计算和统计焊接系数利用因子,确定其分布函数;步骤5,由焊接系数利用因子保证率和设计极值量,确定焊接系数f<sub>weld</sub>,计算焊脚高度。
申请公布号 CN104050357B 申请公布日期 2017.02.22
申请号 CN201410211731.0 申请日期 2014.05.19
申请人 浙江工业大学 发明人 吴剑国;余东方
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 杭州天正专利事务所有限公司 33201 代理人 王兵;黄美娟
主权项 基于有限元分析的船体结构焊脚高度的确定方法,包括如下五个步骤:步骤1,构建船舶结构角焊缝应力的计算模型,提出一种从舱段有限元模型中提取焊缝应力的方法;针对船舶结构角焊缝的受力特点,提出以下七个计算模型:①扶强材与带板的焊缝计算模型;②扶强材端部的焊缝计算模型;③主要支撑构件与外板的焊缝计算模型;④主要支撑构件之间的焊缝计算模型;⑤平面板材或平面舱壁周界的焊缝计算模型;⑥船体梁计算模型;⑦支柱等以拉应力为主的焊缝计算模型;计算模型①、③、⑥匹配的校核准则为下式(1),计算模型⑤、⑦匹配的校核准则为下式(2),计算模型②、④匹配的校核准则为下式(3);<maths num="0001"><math><![CDATA[<mrow><msub><mi>&tau;</mi><mi>f</mi></msub><mo>&le;</mo><msubsup><mi>f</mi><mi>f</mi><mi>w</mi></msubsup><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000011.GIF" wi="958" he="63" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>&sigma;</mi><mi>f</mi></msub><mo>&le;</mo><msub><mi>&beta;</mi><mi>f</mi></msub><msubsup><mi>f</mi><mi>f</mi><mi>w</mi></msubsup><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000012.GIF" wi="974" he="63" /></maths><maths num="0003"><math><![CDATA[<mrow><msqrt><mrow><msup><mrow><mo>(</mo><mfrac><msub><mi>&sigma;</mi><mi>f</mi></msub><msub><mi>&beta;</mi><mi>f</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msubsup><mi>&tau;</mi><mi>f</mi><mn>2</mn></msubsup></mrow></msqrt><mo>&le;</mo><msubsup><mi>f</mi><mi>f</mi><mi>w</mi></msubsup><mo>.</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000013.GIF" wi="1006" he="139" /></maths>式(1)、(2)、(3)中:τ<sub>f</sub>为平均剪应力,<img file="FDA0001122646820000014.GIF" wi="58" he="56" />为焊缝强度设计值,σ<sub>f</sub>为有效截面上平均应力,β<sub>f</sub>为正面角焊缝强度设计值增大系数;根据焊缝的部位选取合适的计算模型,不同的计算模型有不同校核准则,故提取焊缝应力的类型也不同;计算模型①、③、⑥中需提取立板的剪应力;计算模型⑤、⑦需提取立板的正应力;计算模型②、④中需提取立板的剪应力与面板的正应力;步骤2,确定舱段有限元模型的板单元剪应力与焊缝平均剪应力之间的关系;建立了若干尺寸梁的有限元模型,比较800×800mm单元的有限元计算结果与焊缝强度计算公式结果之间的比例关系,约为1.1;步骤3,提出船体结构焊缝强度的应力标准;通常的焊缝强度许用值是根据建筑结构规范确定的,不能直接应用于船舶结构的有限元计算;二者存在两个差别:一是载荷与强度的体系不同;二是计算模型也不同,一个是基于近似的焊缝平均应力公式,一个是基于结构的有限元计算;通过类比,使得船舶结构中焊缝的强度标准随同构件的强度标准变化,得出船舶中各个构件的焊缝许用应力;步骤4,计算和统计焊接系数利用因子,确定其分布函数;4.1用步骤1所提供的计算模型,根据焊缝所处的部位提取出实船舱段有限元模型中所需的板应力;依据假设焊接系数和步骤2中的比例关系将板应力换算为焊缝应力,通过步骤3得到焊缝许用应力;4.2计算焊接系数利用因子;将焊缝应力、焊缝许用应力代入下式(4)即可计算焊接系数利用因子的数值;<img file="FDA0001122646820000021.GIF" wi="1205" he="104" />4.3利用Microsoft Excel中的插入图表功能形成贡贝尔概率纸,对焊接系数的利用因子样本进行贡贝尔Ⅰ型分布检验;4.4根据已形成的极值概率纸,确定贡贝尔Ⅰ型尔分布函数的参数B、C,代入下式(5),即可得贡贝尔Ⅰ型分布函数;p(x)=exp(‑exp(‑B(x‑C))) ‑∞<x<+∞  (5)若对焊接系数利用因子样本乘以系数k,即可得到新的分布函数,如下式(6);<maths num="0004"><math><![CDATA[<mrow><mtable><mtr><mtd><mrow><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>exp</mi><mo>(</mo><mrow><mo>-</mo><mfrac><mi>B</mi><mi>k</mi></mfrac><mrow><mo>(</mo><mrow><mi>x</mi><mo>-</mo><mi>k</mi><mi>C</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow></mrow></mtd><mtd><mrow><mo>-</mo><mi>&infin;</mi><mo>&lt;</mo><mi>x</mi><mo>&lt;</mo><mo>+</mo><mi>&infin;</mi></mrow></mtd></mtr></mtable><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000022.GIF" wi="1342" he="119" /></maths>步骤5,由焊接系数利用因子保证率和设计极值量,确定焊接系数f<sub>weld</sub>,计算焊脚高度;5.1设定保证率p(x<sub>α</sub>),令设计极值x<sub>α</sub>=1,代入上式(6)可计算出系数k;<maths num="0005"><math><![CDATA[<mrow><mi>p</mi><mrow><mo>(</mo><msub><mi>x</mi><mi>&alpha;</mi></msub><mo>)</mo></mrow><mo>=</mo><mi>p</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><mi>exp</mi><mo>(</mo><mrow><mo>-</mo><mfrac><mi>B</mi><mi>k</mi></mfrac><mrow><mo>(</mo><mrow><mn>1</mn><mo>-</mo><mi>k</mi><mi>C</mi></mrow><mo>)</mo></mrow></mrow><mo>)</mo><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000023.GIF" wi="837" he="119" /></maths>5.2根据利用因子的改变系数k,代入下式(7)计算焊接系数f<sub>weld</sub>;<maths num="0006"><math><![CDATA[<mrow><msub><mi>f</mi><mrow><mi>w</mi><mi>e</mi><mi>l</mi><mi>d</mi></mrow></msub><mo>=</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><msub><mi>f</mi><mrow><mi>w</mi><mi>e</mi><mi>l</mi><mi>d</mi><mn>1</mn></mrow></msub><mo>+</mo><mrow><mo>(</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mfrac><mn>2</mn><mrow><mn>0.8</mn><msub><mi>t</mi><mrow><mi>p</mi><mo>-</mo><mi>g</mi><mi>r</mi><mi>s</mi></mrow></msub></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001122646820000024.GIF" wi="1141" he="135" /></maths>式(7)中:f<sub>weld1</sub>为假设焊接系数,t<sub>p‑grs</sub>为立板建造厚度;5.3通过焊接系数f<sub>weld</sub>,代入下式(8)计算焊脚高度l<sub>leg</sub>;l<sub>leg</sub>=f<sub>yd</sub>f<sub>weld</sub>t<sub>p‑grs</sub>+2  (8)式(8)中:f<sub>yd</sub>为考虑焊缝熔敷金属屈服强度的修正因数。
地址 310014 浙江省杭州市下城区潮王路18号