【技术实现步骤摘要】
本专利技术涉及流固耦合，涉及高压水下环境固体结构的模拟和动态响应的计算，具体地，涉及一种基于虚拟流体方法的深海耐压结构内爆计算方法及系统。
技术介绍
1、深海耐压结构液-固-气三相流固耦合的水下内爆，是流固耦合问题中的一种，常见于水下压力容器的测试、深海潜水器和潜艇的研发。对深海耐压结构液-固-气三相流固耦合的水下内爆问题进行数值模拟，可以预测深海耐压结构的动态响应和变形以有材料的各种，从而进一步为深海载人设备的设计和研发提供参考。
2、用于数值模拟深海耐压结构液-固-气三相流固耦合的现有技术有存在以下缺陷：流固边界处需要加密网格，对于流固界面的物理量，仅考虑使用当前时刻的值，不考虑使用下一时刻的值。
技术实现思路
1、针对现有技术中的缺陷，本专利技术提供一种基于虚拟流体方法的深海耐压结构内爆计算方法及系统。
2、根据本专利技术提供的一种基于虚拟流体方法的深海耐压结构内爆计算方法及系统，所述方案如下：
3、第一方面，提供了一种基于虚拟流体方法的深海耐压结构内爆计算方法，所述方法包括：
4、步骤s1：设定流场初始条件，得到深海耐压结构的计算环境；
5、步骤s2：设定固体初始条件，得到深海耐压结构的几何模型；
6、步骤s3：用步骤s1和s2中的流场初始条件以及固体初始条件，或者用上一时间步的已知结果，结合虚拟流体黎曼解算器，迭代求解流固边界，模拟出深海环境水和结构体的相互作用；
7、步骤s4：用虚拟流体
8、步骤s5：用步骤s1中的流场初始条件，或者用上一时间步的已知结果，结合tvd格式，对流场求解，刻画波在流场中的传播过程；
9、步骤s6：用步骤s2中固体初始条件，或者使用上一时间步的已知结果，结合有限元法，对固体网格结点求解，得到结构体的动态响应和变形过程。
10、优选地，所述步骤s1包括：
11、深海耐压结构，所处的流场初始条件为：
12、pw＝ρwgh+p0 1)
13、其中，ρw为海水密度，g为重力加速度，h；深海耐压结构所处的水深，p0为海平面大气压。
14、优选地，所述步骤s2包括：
15、设定固体初始条件：密度ρs，杨氏模量ys,泊松比γs，屈服应力sy0,剪切模量gs；
16、给出得到深海耐压结构的几何模型，耐压结构内部的空气物理状态参数满足：
17、pa＝(γa-1)ρaea 2)
18、其中，pa为气体压强，γa为气体绝热指数，ρa为气体密度，ea为气体内能。
19、优选地，所述步骤s3包括：
20、步骤s3.1：在流固界面处，固体获得流体的压强，ps＝pl，ps为固体压强，pl为流体压强；代入固体求解器式3)～式16)求解新的固体物理状态；
21、采用八节点六面体有限单元，插值函数选择三维等参插值函数：
22、
23、其中，a，b和c分别为待定系数，不同的物理量对应不同的系数；对于任一物理量q(x,y,z)，采用八个单元节点的插值形式：
24、
25、其中，qi为第i个节点物理量；加速度ak、位移xk及其导数的变分表示为：
26、
27、
28、
29、其中，δ为变分符号，下标i表示爱因斯坦求和约定，下标k和上标k均表示第k个单元，m表示沿x轴、y轴或z轴的分量。加速度ak、位移xk及其导数添加下标i后,表示第k个单元在第i个节点的加速度；表示第k个单元在第i个节点的位移；表示第k个单元在m轴上速度；将整个区域划分为ne个单元，
30、
31、其中，f为待求物理量，vn为单元体积，sn为单元表面积，[a]是由插值函数组成的矩阵，其形式为，
32、
33、引入偏导数矩阵[b]，
34、
35、和应力向量[t]，
36、[t]＝(σ11,σ22,σ33,σ23,σ31,σ12)t11)
37、其中，σ11,σ22,σ33,σ23,σ31,σ12分别表示不同方向的应力，下标1对应x轴，下标2对应y轴，下标3对应z轴；对于时间的推进，采用显示时间差分，
38、vn+1＝vn+an(δtn+δtn+1)/2 12)
39、(δr)n+1＝(δr)n+vn+1·δtn+1 13)
40、其中，δtn和δtn+1分别为n和n+1时刻的时间步长，δr和v分别为位移和速度。
41、用沙漏矢量抑制沙漏的产生，
42、
43、
44、
45、其中，下标k表示网格节点，i表示x轴、y轴或z轴的方向，α表示沙漏的各个方向，a＝qρl2c/4，q为常数，c是固体声速；l是单元特征长度；γαk是沙漏基矢量，γαk为人工粘性系数，为人工粘性项，nk,i为形函数，giα为沙漏质量，表示速度。
46、步骤s3.2：将新的固体加速度as代入虚拟流体黎曼解算器式14)～16)，求解新的流体压强pl；
47、沿着右特征线
48、
49、的相容关系为：
50、
51、将界面两侧的ρc视为常数，
52、
53、其中，ui、pi、ρil分别为界面流体速度、界面流体压力、界面流体密度，ul为流体相对于流固界面的速度；cil表示界面流体声速；pl表示界面附近流体压强；ρl表示表示界面附近流体密度；cl表示表示界面附近流体声速。
54、步骤s3.3：回到步骤s3.1，直到误差小于设定值。
55、优选地，所述步骤s4包括：
56、定义为符号距离函数，
57、
58、其中，ω+表示虚拟区域，ω-表示真实区域，表示物质界面，dis表示距离，表示位置矢量。
59、将步骤s3得到的流固界面压强和流固界面速度，用下述方程外推：
60、
61、其中，q为需要外推的函数；为流固界面法向单位向量，从流体指向固体；t为人工时间；把i从ω-外推到ω+时，方程18)取+；把i从ω+外推到ω-时，方程18)取-。
62、优选地，所述步骤s5包括：
63、使用二阶tvd格式，设待求网格节点为(i0,j,k)，流场三个方向的通量f，g，h都需要离散，
64、f＝(ρu,ρu2+p,ρuv,ρuw,u(e+p))t 22)
65、g＝(ρv,ρvu,ρv2+p,ρvw,v(e+p))t 本文档来自技高网...
【技术保护点】
1.一种基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，包括：
2.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤S1包括：
3.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤S2包括：
4.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤S3包括：
5.根据权利要求4所述的基于虚拟流体方法的深海耐压结构内爆计算方法特征在于，所述步骤S4包括：
6.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法特征在于，所述步骤S5包括：
7.根据权利要求4所述的基于虚拟流体方法的深海耐压结构内爆计算方法特征在于，所述步骤S6中使用公式3)至16)求解结构体动态响应曲线和结构体的变形。
8.一种基于虚拟流体方法的深海耐压结构内爆计算系统，其特征在于，包括：
9.根据权利要求8所述的基于虚拟流体方法的深海耐压结构内爆计算系统，其特征在于，所述模块M1包括：
...
【技术特征摘要】
1.一种基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，包括：
2.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤s1包括：
3.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤s2包括：
4.根据权利要求1所述的基于虚拟流体方法的深海耐压结构内爆计算方法，其特征在于，所述步骤s3包括：
5.根据权利要求4所述的基于虚拟流体方法的深海耐压结构内爆计算方法特征在于，所述步骤s4包括：
6.根据权利要求1所...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。