2. 哈尔滨工程大学 船舶工程学院,黑龙江 哈尔滨 150001;
3. 哈尔滨工程大学 青岛船舶科技有限公司,山东 青岛 266400
2. College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China;
3. Harbin Engineering University, Qingdao Ship Science and Technology Co., Ltd., Qingdao 266400, China
我国部分水面舰船由于服役海区海况条件较好,实际航行时间也低于设计要求,同时加上维护保养得当,舰船平台状态良好,有很大的潜力超期服役[1]。此类超役舰船一部分会被改装为实体靶船,用以部队训练、试验等科研活动。极限强度是舰船结构强度的重要方面,也是实体靶船生命力的重要保障,关系到靶船航行安全及试训活动组织。目前,国军标等相关水面舰艇规范及标准缺少实体靶船极限强度评估的明确方法[2]。
船体的极限承载能力一般指总纵弯曲作用下船体梁的极限弯矩,常用的计算方法有解析法、简化逐步破坏法和非线性有限元法等[3]。其中,非线性有限元法通过建立结构有限元模型,模拟结构的初始缺陷,并考虑材料非线性和几何非线性,计算精度高,被国际船舶与海洋工程结构大会(ISSC)作为极限承载能力预报的衡准方法。如果考虑材料机械性能、结构尺寸和腐蚀等参数的不确定性,极限承载能力为符合一定概率分布的随机变量[4]。响应面法结合蒙特卡罗法随机抽样是一种求解随机变量统计特征值的较好方法,该方法先通过有限元数值计算并采用最小二乘法等方法拟合一个响应面来替代未知的、真实的极限状态方程,再根据随机变量的概率分布运用蒙特卡罗法抽样获得目标变量的概率分布,二者的结合使结构分析的工作量大大减少,克服了传统蒙特卡罗法庞大的有限元计算工作量和传统响应面因为线性化带来误差的缺点[5]。
为了合理地评估实体靶船的承载能力,考虑材料机械性能和结构腐蚀等参数的不确定性,本文将采用响应面-蒙特卡罗法结合非线性有限元法研究建立一种实体靶船极限强度的计算方法,分析响应面设计点的选取方案,给出实体靶船极限承载能力的概率分布。
1 基本原理与方法 1.1 可靠性分析的响应面-蒙特卡罗法工程结构的失效概率可以表示为[6]:
${P_f} = P\left\{ {G(X) < 0} \right\} = \int_{{D_f}} {f\left( X \right){\rm{d}}X}\text{。} $ | (1) |
式中:
用蒙特卡罗法计算结构可靠性时,式(1)可表达为:
${\hat P_f} = \frac{1}{N}\sum\limits_{i = 1}^N {I\left[ {G\left( {{{{\hat X}}_i}} \right)} \right]} \text{。}$ | (2) |
式中:
$\varepsilon = \frac{{\left| {{{\hat P}_f} - {P_f}} \right|}}{{{P_f}}} < 2\sqrt {\frac{{\left( {1 - {{\hat P}_f}} \right)}}{{N{{\hat P}_f}}}} \text{,}$ | (3) |
当给定
响应面法通过拟合一个响应面
$\hat Z = \hat G(x,y) \approx \sum\limits_{{\begin{array}{l} {i + j = 0}\\ {i + j \leqslant 3} \end{array}}}^3 {{p_{ij}}{x^i}{y^j}} \text{。}$ | (4) |
式中:
为使响应面函数更好地逼近真实状态,设计点选取非常重要。一方面尽量要少,以降低计算工作量,另一方面还要尽量包含较多极限状态信息,以保证响应面的拟合精度。
响应面-蒙特卡罗法[7]首先根据随机变量的概率分布选取确定功能函数值的设计点。设计点一般以随机变量的均值点为中心,根据
基于非线性有限元法的船体极限强度计算精度主要取决于模型构建的准确性和非线性求解器的选择与参数设置[8]。
为了避免边界条件的影响,模型范围一般选择舱段模型。板和主要纵向构件应采用4节点单元,且尽量保证网格形状接近正方形。为了准确模拟极限状态结构的塑性大变形和失效模式演化,骨材腹板一般要求至少布置3~4个单元,并以此为基准划分整个舱段模型的网格往往可以取得较好的效果。同时,极限强度计算需要考虑结构的初始挠度[9],包括:
板架的整体变形
${w_{S0}} = {B_0}\sin \frac{{\text{π} x}}{a}\sin \frac{{\text{π} y}}{S}\text{;}$ |
板的局部变形
${w_{P0}} = {A_0}\sin \frac{{m\text{π} x}}{a}\sin \frac{{\text{π} y}}{b}\text{;}$ |
骨材的侧倾变形
${w_{T0}} = {C_0}\frac{y}{{{h_w}}}\sin \frac{{\text{π} x}}{a}\text{。}$ |
式中:
由于极限强度计算需要考虑材料的非线性,在缺少船用钢材料拉伸曲线的情况下,一般采用理想弹塑性或理想塑性强化模型以考虑材料的非线性效应。
目前,极限强度分析的非线性求解方法主要有Newton-Raphson迭代法、Riks法和显式动态算法等。其中,Newton-Raphson法以载荷增量为控制量,根据节点力平衡判断收敛,但结构大变形引起的节点力不平衡容易造成不收敛。Riks法采用弧长作为控制变量,改进了Newton-Raphson法的迭代策略,计算效率和收敛性主要取决于弧长相关参数的设置。显式动态算法采用中心差分法,用上一步和当前步的结果计算下一步的计算结果,无需迭代运算,收敛性好,计算效率较高,目前被广泛地应用于结构静态和准静态的非线性响应分析。
2 实体靶船极限强度可靠性分析 2.1 计算模型以某型超役水面舰艇为例,该船已服役近30年,将改装成为实体靶船。考虑材料屈服限和结构腐蚀的随机变化,采用响应面-蒙特卡罗法结合非线性有限元法,开展了极限强度可靠性分析。
选取船中舱段为分析对象,采用Abaqus通用有限元软件建立有限元模型,网格尺寸整体上取50 mm×50 mm,同时纵骨腹板划分3个单元,单元类型选用4节点薄壳单元S4R,如图1所示。该舰材料屈服限为235 MPa,杨氏模量为206 GPa,泊松比0.3
由于是实体靶船,因此结构腐蚀是重点考虑的一个随机变量,同时由于材料屈服限对于舰船极限强度影响较大,因此主要考虑腐蚀量和材料屈服限的随机变化。材料屈服限服从正态分布,变异系数取0.05,即
国内外学者通过对结构腐蚀统计分析和试验研究,提出了一些腐蚀模型[10]。秦圣平等[11]在现有腐蚀模型比较分析的基础上,提出了一种适用于船体结构时变可靠性分析的非线性腐蚀模型,该模型可以较好模拟钢结构在腐蚀环境下的腐蚀损伤过程。表达式为:
$d(t) = \left\{ \begin{array}{l} 0,0 \leqslant t \leqslant {T_{st}}, \\ {d_\infty }\left\{ {\left. {1 - \exp \left[ - {{\left(\frac{{t - {T_{st}}}}{\eta }\right)}^\beta }\right]} \right\}} \right. ,{\rm{ }}{T_{st}} \leqslant t \leqslant {T_L} \text{。} \\ \end{array} \right.$ | (5) |
式中:
将表中相关参数代入腐蚀模型,采用蒙特卡罗抽样,抽样结果如图2所示,获得了服役30年的腐蚀量服从正态分布,
采用显式动态算法分析该舰中垂状态的极限强度,图3为舱段的弯矩-转角曲线,图4为极限状态结构应力云图。可以发现,该舰的中垂极限弯矩为36.21 MN·m,极限状态该舰强横梁间的甲板板架发生了屈曲破坏。
根据材料屈服限和结构腐蚀量的概率分布及其统计特征值,采用
从舱段极限强度的计算结果可以发现:随着材料屈服限的提高,极限强度逐步增加;随着结构腐蚀量的增加,极限强度不断降低。总体而言,材料屈服限对极限强度的影响较腐蚀量的影响要大。
在此基础上,采用三次多项式表达式对响应面进行拟合,设计了10种不同设计点方案,分别包括49,25,21,21,17,17,17,13,13和13个设计点,如图5所示。表3为相关设计点方案的响应面拟合结果,图6为部分设计点方案拟合得到的响应面。
为了评价响应面拟合效果,引入了局部拟合优度
$R_{Loc}^2 = 1 - \frac{{\displaystyle\sum\limits_{i = 1}^k {{{({z_i} - {{\hat z}_i})}^2}} }}{{\displaystyle\sum\limits_{i = 1}^k {{{({z_i} - \bar z)}^2}} }}\text{,}$ |
$R_G^2 = 1 - \frac{{\displaystyle\sum\limits_{i = 1}^{49} {{{({z_i} - {{\hat z}_i})}^2}} }}{{\displaystyle\sum\limits_{i = 1}^{49} {{{({z_i} - \bar z)}^2}} }}\text{。}$ |
式中:
从响应面的拟合结果可以发现:前9个方案的局部拟合效果均大于0.99;同时,前8个方案的整体拟合优度均大于0.99,且方案1、方案4、方案6和方案8达到了0.9998,对应的设计点数量分别为49,21,17和13,说明响应面的整体拟合优度与设计点的选取息息相关,并不是设计点数量越多越好。方案10的设计点选取没有考虑材料屈服限和腐蚀量的相关性,拟合得到的响应面已失真,如图6(d)所示,局部拟合优度仅0.482。方案8虽然只选取了13个设计点,但其整体拟合优度与49个设计点一致,达到了0.9998,局部拟合优度甚至达到了1.0,拟合优度比方案2(25个设计点)和方案3(21个设计点)都要好。
因此,响应面设计点的合理选取,不仅可以有效降低计算工作量,同时可以保证响应面拟合精度。
2.4 极限强度的可靠性分析获得极限强度的响应面后,结合相关参数的概率分布,采用蒙特卡罗法随机抽样,进一步分析实体靶船极限强度的概率分布。
首先,开展了蒙特卡罗法的收敛性分析,抽样次数分别取1000,10000,100000和1000000次,样本的统计特征值见表5。可以发现,当样本数量达到10000次,极限强度的均值与标准差已趋于稳定;当样本数量增加到100000次或更高,极限强度的统计特征值变化不大。
在正态分布的基础上,进一步采用对数正态分布和Weibull分布拟合极限强度的概率分布,拟合结果如图7所示。其中,正态分布和对数正态分布的分布函数分别为
可以发现:极限强度的随机分布与Weibull分布偏差较大。因此,实体靶船极限强度的概率分布更趋近于正态分布或对数正态分布。
3 结 语本文采用蒙特卡罗法与响应面法相结合的可靠性分析方法,开展考虑材料屈服限和腐蚀影响的舱段极限强度研究,研究表明:在保证拟合优度的前提下,响应面设计点的合理选取,可以有效降低计算工作量;考虑材料屈服限和腐蚀联合影响的实体靶船极限强度概率分布更趋近于正态分布或对数正态分布。本文研究建立的基于响应面-蒙特卡罗法的可靠性分析方法,对于实体靶船极限强度评估具有指导意义。
[1] |
张克辉, 滑林. 老龄舰船总纵强度状态评估可行性研究[J]. 海军工程大学学报, 2015, 27(5): 94-98. |
[2] |
中国人民解放军总装备部. 舰船通用规范 1组 船体结构: GJB4000-2000[S]. 上海: 海军论证中心标准规范所舰船通用办公室, 2000.
|
[3] |
李陈峰, 高超, 马开开, 等. 考虑非对称剖面中和轴偏转的改进Smith法研究[J]. 船舶力学, 2019, 23(7): 810-819. DOI:10.3969/j.issn.1007-7294.2019.07.007 |
[4] |
任慧龙, 李陈峰, 李辉, 等. 破损舰船剩余强度的可靠性评估方法研究[J]. 船舶力学, 2010, 14(7): 757-764. DOI:10.3969/j.issn.1007-7294.2010.07.008 |
[5] |
李远瑛, 张德生. 基于响应面和蒙特卡罗法结构位移可靠度[J]. 辽宁工程技术大学学报(自然科学版), 2011, 30(3): 392-395. |
[6] |
赵国藩, 金伟良, 贡金鑫. 结构可靠度理论[M]. 北京: 中国建筑工业出版社, 2000.
|
[7] |
苏成, 李鹏飞, 韩大建. 基于Neumann展开响应面技术的重要抽样蒙特卡罗法[J]. 工程力学, 2009, 26(12): 1-5+11. |
[8] |
李陈峰, 任慧龙, 冯国庆, 等. 外载荷作用下的船体梁极限弯矩[J]. 哈尔滨工程大学学报, 2010, 31(1): 26-31. DOI:10.3969/j.issn.1006-7043.2010.01.005 |
[9] |
中国船级社. 船体梁极限强度非线性有限元方法评估指南[S]. 中国船级社, 2014.
|
[10] |
GUEDES SOARES C, GARBATOV Y. Reliability of maintained, corrosion protected plates subjected to non-linear corrosion and compressive loads[J]. Marine Structures, 1999, 12(06): 425-445. DOI:10.1016/S0951-8339(99)00028-3 |
[11] |
秦圣平, 崔维成, 沈凯. 船舶结构时变可靠性分析的一种非线性腐蚀模型[J]. 船舶力学, 2003, 7(1): 94-103. DOI:10.3969/j.issn.1007-7294.2003.01.012 |