三维声弹性广义水动力系数的高效并行算法制造技术

本发明专利技术涉及三维声弹性广义水动力系数的高效并行算法，包括如下步骤：基于有限元软件得到船舶干模态的振型位移矩阵；将频率总数和总进程数分别分至各个频率分组中；对各频率分组内的各个频率得到Green函数矩阵及其偏导数矩阵；针对单个频率构建湿面元源强矩阵方程，对该方程中Green函数偏导数矩阵和振型位移法向分量矩阵进行二维并行划分，计算得湿面元源强矩阵；计算广义水动力系数中干模态附连水质量和附连水阻尼。本发明专利技术针对船舶三维声弹性分析方法，将频率总数分至各频率分组中并行计算，并对Green函数偏导数矩阵和振型位移法向分量矩阵进行二维并行划分，有效解决了计算量大的瓶颈问题，大大提升了船舶三维声弹性分析方法的大规模工程计算适用性。

Efficient parallel algorithm for generalized hydrodynamic coefficients of three-dimensional acoustoelasticity

【技术实现步骤摘要】
三维声弹性广义水动力系数的高效并行算法
本专利技术涉及船舶水下辐射噪声计算与声学边界元计算
，尤其是一种三维声弹性广义水动力系数的高效并行算法。
技术介绍
船舶三维声弹性频域分析方法继承了经典的浮体水弹性力学理论(Wu,Hydroelasticityoffloatingbodies,1984)的流固耦合处理方式，选用易于求解且具有正交完备性的干模态(结构在真空中的模态)作为广义基函数，通过计算干模态的广义水动力系数(包括附连水质量和附连水阻尼)，形成以干模态主坐标为待求量的频域内的声弹耦合广义动力学方程；其中广义水动力系数，即附连水质量和附连水阻尼的计算是形成频域内声弹耦合广义动力学方程的关键。现有技术中的船舶三维声弹性分析方法中，采用边界积分方法中的简单源数值方法计算广义水动力系数。简单源数值方法是将船舶湿表面离散成N块四边形的湿面元，每个湿面元的作用等效到其几何中心点上。每个湿面元相当于是一个单极子，通过求解一个N×N维度的矩阵方程进而得到每个湿面元的源强。获取源强之后，再进行干模态广义水动力系数的计算。当干模态数量较多时，高斯消元法作为一种可以采用的计算量和计算稳定性综合较优的算法被用于矩阵求解。当船舶湿面元数目N较大时(例如N>100000)，上述湿面元源强的计算量将非常大，耗时很长，严重影响到了船舶三维声弹性分析方法的工程实际应用时的计算速度和计算时长，成为一个瓶颈问题。
技术实现思路
本申请人针对上述现有生产技术中的缺点，提供一种结构合理的三维声弹性广义水动力系数的高效并行算法，从而通过高效的并行计算方法，有效提高船舶三维声弹性广义水动力系数的计算速度，工程适用性强。本专利技术所采用的技术方案如下：一种三维声弹性广义水动力系数的高效并行算法，包括如下步骤：第一步：基于有限元软件计算船舶结构的干模态，得到干模态的振型位移矩阵；第二步：设计算的频率总数为K，频率分组数为ngro，基于负载均衡，将频率总数K分至各频率分组内，设每组频率分组内分到的频率数为Kn，则ngro个频率分组中频率数Kn的总和为频率总数K，其中n＝1,...,ngro；完成频率总数进行并行计算的分组；设总进程数为npro，则每组频率分组分到的进程数为ngp＝npro/ngro；第三步：针对第二步中的每组频率分组，基于简单源数值方法，将船舶湿表面划分成N个湿面元，以表示场点坐标、表示源点坐标，计算得到第k个频率ωk上的Green函数矩阵以及Green函数偏导数矩阵，每个矩阵的维度为N×N；Green函数矩阵为：其中，k＝Kn0,…,Kn0+Kn-1，Kn0为第n个频率分组中的起始频率号，Kn为该第n个频率分组中的频率数；i,j＝1,2,…,N；Green函数偏导数矩阵为：实现每组频率分组中各个频率的Green函数和Green函数偏导数矩阵的一维并行处理；第四步：基于第三步中的Green函数偏导数矩阵，构建第k个频率ωk上的湿面元源强的矩阵方程，如下：其中，为虚数符号，({D1},…,{Dm})为干模态的振型位移法向分量矩阵，({σ1},…,{σm})为与各个频率相对应的待求解的湿面元源强矩阵，m为截断的干模态数；对Green函数偏导数矩阵和干模态的振型位移矩阵的法向分量矩阵进行二维并行划分计算，从而获得各组频率分组中各个频率对应的湿面元源强矩阵({σ1},…,{σm})；第五步：根据第四步中获得的与各个频率对应的湿面元源强矩阵，进行各个频率对应的广义水动力系数的计算，包括干模态附连水质量和附连水阻尼的计算，如下：其中，ρ0为流体密度，[A(ωk)]为附连水质量矩阵，[B(ωk)]为附连水阻尼矩阵，Re{}表示取复数的实部，Im{}表示取复数的虚部，上标T表示取矩阵的转置，为湿面元面积的对角矩阵；通过计算式的求解，从而获得干模态附连水质量和附连水阻尼系数。作为上述技术方案的进一步改进：第四步中，对Green函数偏导数矩阵进行二维并行划分的方法包括如下步骤：步骤一：设单个频率分组内的行进程数为Nrow，该组频率分组内的列进程数为Ncol，则该频率分组中的进程数ngp满足以下关系：ngp＝Nrow×Ncol；步骤二：基于行进程数Nrow和列进程数Ncol，对该频率分组中第k个频率ωk上的Green函数偏导数矩阵中的各个元素进行行向和列向的间隔取数，形成与进程数ngp对应的各个分块，用myid对各个分块进行标号，则myid＝0,....,ngp-1；从而实现第k个频率ωk上的Green函数偏导数矩阵的二维并行划分；步骤三：同步骤二对干模态的振型位移矩阵的法向分量矩阵({D1},…,{Dm})进行二维并行划分；进而求得该第k个频率ωk上的湿面元源强矩阵({σ1},…,{σm})；步骤四：对步骤一中频率分组中的Kn个频率，按照步骤二至步骤三依次重复，从而计算得到该频率分组内与各个频率相对应的湿面元源强矩阵({σ1},…,{σm})。第二步中，将K个频率分配至ngro个频率分组中，设ωk为第k个频率，则第n个频率分组中的频率ωk满足k＝Kn0,…,Kn0+Kn-1，其中Kn0为该第n个频率分组中的起始频率号，Kn为该第n个频率分组中的频率数。第二步中总进程数npro被频率分组数ngro整除。第二步中将各个频率分至各频率分组，各组频率分组并行各自计算；各组频率分组各自经第三步、第四步和第五步，获得各自频率分组内各个频率对应的湿面元源强矩阵，并最终获得各个频率对应的干模态附连水质量和附连水阻尼系数。第四步中干模态的振型位移法向分量矩阵为第一步中干模态的振型位移矩阵的一个分量，振型位移法向分量矩阵在振型位移矩阵基础上获得。第五步中，对Green函数矩阵采用行并行划分，进行一维并行计算处理。本专利技术的有益效果如下：本专利技术结构紧凑、合理，操作方便，通过在船舶三维声弹性分析中提出二层并行的高效算法，有效解决了计算量大的瓶颈问题，大大提升了广义水动力系数的计算速度，适于大规模的工程计算应用；第一层并行为频率级的并行，将频率总数分至各频率分组中并行计算，以频率分组为基并行计算，整体上实现了各频率的高效一维并行计算；第二层并行为：在进行湿面元源强矩阵求解时，通过间隔取数的方式分别对Green函数偏导数矩阵与干模态振型位移矩阵进行分块，实现了二维并行划分，有效助力于实现高斯消元法求解时的负载均衡以及通信和计算隐藏。本专利技术还包括如下优点：在第五步中求解干模态附连水质量和附连水阻尼系数时，可对Green函数矩阵采用行并行划分的一维并行计算，以进一步提升整体的计算速度。附图说明图1为本专利技术三维声弹性广义水动力系数高效并行算法的流程图。图2为本专利技术基于频率级并行的二级并行处理示意图。图3为本专利技术Green函数和Green函数本文档来自技高网...

【技术保护点】
1.一种三维声弹性广义水动力系数的高效并行算法，其特征在于：包括如下步骤：/n第一步：基于有限元软件计算船舶结构的干模态，得到干模态的振型位移矩阵；/n第二步：设计算的频率总数为K，频率分组数为n

【技术特征摘要】
1.一种三维声弹性广义水动力系数的高效并行算法，其特征在于：包括如下步骤：
第一步：基于有限元软件计算船舶结构的干模态，得到干模态的振型位移矩阵；
第二步：设计算的频率总数为K，频率分组数为ngro，基于负载均衡，将频率总数K分至各频率分组内，设每组频率分组内分到的频率数为Kn，则ngro个频率分组中频率数Kn的总和为频率总数K，其中n＝1,...,ngro；完成频率总数进行并行计算的分组；
设总进程数为npro，则每组频率分组分到的进程数为ngp＝npro/ngro；
第三步：针对第二步中的每组频率分组，基于简单源数值方法，将船舶湿表面划分成N个湿面元，以表示场点坐标、表示源点坐标，计算得到第k个频率ωk上的Green函数矩阵以及Green函数偏导数矩阵，每个矩阵的维度为N×N；
Green函数矩阵为：



其中，k＝Kn0,…,Kn0+Kn-1，Kn0为第n个频率分组中的起始频率号，Kn为该第n个频率分组中的频率数；i,j＝1,2,…,N；
Green函数偏导数矩阵为：



实现每组频率分组中各个频率的Green函数和Green函数偏导数矩阵的一维并行处理；
第四步：基于第三步中的Green函数偏导数矩阵，构建第k个频率ωk上的湿面元源强的矩阵方程，如下：



其中，为虚数符号，({D1},…,{Dm})为干模态的振型位移法向分量矩阵，({σ1},…,{σm})为与各个频率相对应的待求解的湿面元源强矩阵，m为截断的干模态数；
对Green函数偏导数矩阵和干模态的振型位移矩阵的法向分量矩阵进行二维并行划分计算，从而获得各组频率分组中各个频率对应的湿面元源强矩阵({σ1},…,{σm})；
第五步：根据第四步中获得的与各个频率对应的湿面元源强矩阵，进行各个频率对应的广义水动力系数的计算，包括干模态附连水质量和附连水阻尼的计算，如下：






其中，ρ0为流体密度，[A(ωk)]为附连水质量矩阵，[B(ωk)]为附连水阻尼矩阵，Re{}表示取复数的实部，Im{}表示取复数的虚部，上标T表示取矩阵的转置，为湿面元面积的对角矩阵；
通过计算式的求解，从而获得干模态附连水质量和附连水阻尼系数。


2.如权利要求1所述的三维声弹性广义水动力系...

【专利技术属性】
技术研发人员：邹明松，
申请(专利权)人：中国船舶科学研究中心中国船舶重工集团公司第七零二研究所，
类型：发明
国别省市：江苏;32