一种潜水器圆柱壳无模自由胀形过程的性能预测方法

摘要

本发明公开了一种潜水器圆柱壳无模自由胀形过程的性能预测方法，设计计算了轴向曲率半径、周向曲率半径、表面积、中点厚度、中点等效应变、内部体积、胀形压力和中点等效应力在无模自由胀形过程中的变化规律，为潜水器圆柱壳无模自由胀形工艺参数设计和优化提供理论数据支撑。

说明书

技术领域

本发明属于板壳塑性加工技术领域，特别涉及了潜水器圆柱壳无模自由胀形工艺。

背景技术

潜水器是大洋勘查开发的重要装备，耐压壳是潜水器关键部件和浮力单元，在潜水器下潜过程中耐压壳可以确保内部设备不被破坏以及工作人员的安全，同时提供正浮力。圆柱壳是耐压壳的基础承压单元，具有空间利用率高、水动力学特性好、设计计算方便等优点。

然而，圆柱壳极限承载能力对初始几何缺陷非常敏感，导致安全性下降。通过对圆柱壳施加内压，进行圆柱壳无模自由胀形，可以消除过大初始几何缺陷、增加材料屈服强度、平衡壳体各处力学性能，进而提高圆柱壳的极限承载能力，提高壳体安全性，同时也增大了内部空间。但是，业界仍然缺乏有关圆柱壳无模自由胀形工艺的性能预测方法。

发明内容

为了解决上述背景技术提到的技术问题，本发明提出了一种潜水器圆柱壳无模自由胀形过程的性能预测方法。

为了实现上述技术目的，本发明的技术方案为：

一种潜水器圆柱壳无模自由胀形过程的性能预测方法，包括以下步骤：

(1)设计圆柱壳轴向曲率半径和周向曲率半径的计算式，并将胀形量代入计算式中，得到圆柱壳轴向曲率半径与周向曲率半径随着胀形量变化的历程；

(2)设计圆柱壳表面积的计算式，并将胀形量代入计算式中，得到圆柱壳表面积随着胀形量变化的历程；

(3)设计圆柱壳中点厚度的计算式，并将胀形量代入计算式中，得到圆柱壳中点厚度随着胀形量变化的历程；

(4)设计圆柱壳内部体积的计算式，并将胀形量代入计算式中，得到圆柱壳内部体积随着胀形量变化的历程；

(5)设计圆柱壳中点等效应变的计算式，并将胀形量代入计算式中，得到圆柱壳中点等效应变随着胀形量变化的历程；

(6)设计圆柱壳胀形压力的计算式，并将胀形量代入计算式中，得到圆柱壳胀形压力随着胀形量变化的历程；

(7)设计圆柱壳中点等效应力的计算式，并将胀形量代入计算式中，得到圆柱壳中点等效应力随着胀形量变化的历程。

进一步地，在步骤(1)中，圆柱壳轴向曲率半径和周向曲率半径的计算式如下：

上式中，R

进一步地，在步骤(2)中，圆柱壳表面积的计算式如下：

上式中，y为胀形曲线，x为圆柱壳中点垂直方向位移，S为圆柱壳胀形表面积，y′

进一步地，在步骤(3)中，圆柱壳中点厚度的计算式如下：

上式中，t为圆柱壳中点厚度，t

进一步地，厚度系数k的表达式如下：

k＝2.79154E-06h

上式中，E为弹性模量。

进一步地，在步骤(4)中，圆柱壳内部体积的计算式如下：

上式中，V

进一步地，在步骤(5)中，圆柱壳中点等效应变的计算式如下：

上式中，ε

进一步地，在步骤(6)中，圆柱壳胀形压力的计算式如下：

上式中，σ

进一步地，在步骤(7)中，圆柱壳中点等效应力的计算式如下：

上式中，

进一步地，在步骤(1)-(7)中，胀形量的取值范围为胀形半径的百分比。

采用上述技术方案带来的有益效果：

本发明设计计算轴向曲率半径、周向曲率半径、表面积、中点厚度、中点等效应变、内部体积、胀形压力和中点等效应力在无模自由胀形过程中的变化规律，为潜水器圆柱壳无模自由胀形工艺参数设计和优化提供理论数据支撑。本发明的圆柱壳无模自由胀形性能预测公式可推广到其他材料或领域的圆柱壳无模自由胀形预测。

附图说明

图1为本发明的整体方法流程图；

图2为圆柱壳无模自由胀形结构示意图；

图3为胀形量h与轴向曲率半径R

图4为胀形量h与圆柱壳胀形表面积S的曲线图；

图5为胀形量h与圆柱壳胀形中点厚度t的曲线图；

图6为胀形量h与圆柱壳胀形体积V′的曲线图；

图7为胀形量h与圆柱壳胀形中点等效应变ε

图8为胀形量h与圆柱壳胀形压力p的曲线图；

图9为胀形量h与圆柱壳胀形中点等效应力σ

图10为圆柱壳无模自由胀形加压原理图。

具体实施方式

以下将结合附图，对本发明的技术方案进行详细说明。

根据图1所示的圆柱壳无模自由胀形性能预测方法流程图，对本发明作进一步的详细阐述。实施例的圆柱壳的尺寸参数如表1所示，壳体胀形结构如图2所示，壳体材料为不锈钢，其材料参数表2所示。具体实施过程采用解析解的方法，进行操作实现，并通过三个实物柱壳的胀形试验，验证解析方法的正确性。

表1.圆柱壳公称尺寸

表2.标准不锈钢材料试样的弹塑性参数

σ

步骤一(S1)：计算圆柱壳轴向曲率半径和周向曲率半径，公式如下：

其中：R

具体操作如下：

a.将半柱高l＝50mm代入轴向曲率半径R

b.将圆柱壳初始外圆半径r

c.将胀形量h(1-30mm)分别代入轴向曲率半径R

步骤二(S2)：计算圆柱壳表面积，公式如下：

其中：y为圆柱壳水平方向位移，x为圆柱壳中点垂直方向位移，S为圆柱壳胀形表面积。

具体操作如下：

a.基于上述步骤一求出的轴向曲率半径R

b.再将胀形半径y公式代入圆柱壳胀形表面积S计算公式，得到圆柱壳胀形表面积随着胀形量变化的解析方程，即：

c.将胀形量h(1-30mm)代入圆柱壳胀形表面积S解析方程求解并在商业软件Origin中，输出圆柱壳表面积随着胀形量变化的历程，即胀形表面积-胀形量曲线，如图4所示。

步骤三(S3)：计算圆柱壳中点厚度，公式如下：

其中，k为厚度系数，表达式为k＝2.79154E-06h

具体操作如下：

a.基于上述步骤二求出的圆柱壳胀形半径r，厚度系数k，圆柱壳初始外圆半径r

b.将厚度系数k表达式，胀形量h(1-30mm)代入圆柱壳中点厚度t解析方程求解并在商业软件Origin中，输出圆柱壳中点厚度t随着胀形量变化的历程，即中点厚度-胀形量曲线，如图5所示。

步骤四(S4)：计算圆柱壳内部体积，公式如下：

其中，V

具体操作如下：

a.在初始状态下，胀形量h为零，将圆周率π，半柱高l＝50mm，圆柱壳初始外圆半径r

b.基于上述步骤一求出的轴向曲率半径R

c.基于上述步骤一求出的轴向曲率半径R

步骤五(S5)：计算圆柱壳中点等效应变，公式如下：

其中，ε

具体操作如下：

a.在整个胀形过程中，圆柱壳是各向同性的。同一点的法向应变，周向应变，纬向应变三者之和为零，因此可用法向应变ε

b.将周向应变ε

c.基于上述步骤二求出的胀形半径r，步骤三求出的中点厚度t随胀形量h(1-30mm)的解代入圆柱壳中点等效应变ε

步骤六(S6)：计算圆柱壳胀形压力，公式如下：

其中，σ

具体操作如下：

a.基于上述步骤一求出的轴向曲率半径R

b.将胀形量h(1-30mm)，圆柱壳屈服强度σ

c.在商业软件Origin中，输出圆柱壳胀形压力p随着胀形量变化的历程，即胀形压力-胀形量曲线，如图8所示。

步骤七(S7)：计算圆柱壳中点等效应力，公式如下：

其中，

具体操作如下：

a.将圆柱壳周向应力

b.基于上述步骤一所求的圆柱壳轴向曲率半径R

c.在商业软件Origin中，输出圆柱壳中点等效应力σ

为了验证上述数值计算方法的正确性，根据表1参数，加工了三个不锈钢圆柱壳，如图10所示，分别施加8.5MPa、10.5MPa、12MPa的胀形压力。采用三维扫描仪测得三种胀形压力下的外表面积和内部体积数据，分别如图4和6所示；采用无损测厚仪测得三种胀形压力下的壁厚数据，如图5所示。由图4、5、6、8可知，解析结果和试验结果具有良好的一致性，证明了本发明方法的正确性。

实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案基础上所做的任何改动，均落入本发明保护范围之内。