流体介质中复杂组合壳结构振动响应定量计算方法及系统技术方案

本发明专利技术公开了一种流体介质中复杂组合壳结构振动响应定量计算方法及系统，其中方法包括以下步骤：获取流体介质中组合壳系统的结构参数和激励参数；根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；根据组合壳系统两端边界条件和结构连续条件，建立流体介质-组合壳系统的动力学模型；通过Moore-Penrose广义逆矩阵求解得到流体介质-组合壳系统的振动响应。本发明专利技术基于精细积分和传递矩阵法，可以定量计算流体介质中复杂组合壳受外激励产生的振动响应，为水下航行器的振动特性分析提供参考和依据。

【技术实现步骤摘要】
流体介质中复杂组合壳结构振动响应定量计算方法及系统
本专利技术涉及船舶结构振动响应分析
，具体设计一种基于精细算法的复杂组合壳结构振动响应定量计算方法。
技术介绍
船舶特别是水下航行器振动噪声一直是国内外研究的热点，一方面其会降低船舶的舒适性、另一方面也影响船舶的安全性、还损害仪器设备的使用性能，当然还危害舰船的声隐身性。随着船舶的大功率化、高速化和轻型化，振动噪声问题也日益突出、另外随着船员和乘客舒适性的不断提高，振动噪声问题成为亟待解决的难题。因此，在船舶结构声学设计时，需要进行振动响应计算。船舶结构振动计算分为自由和强迫振动两个方面内容，前者计算船舶结构的固有频率和固有阵型，后者计算船舶结构在激振力作用下的振动变形。船舶结构振动响应计算通常采用的计算方法包括经验公式、解析法、有限元法和传递矩阵法。现有的船舶总体自由振动计算方法(CB/T3472-92)规定了船舶总体2～5节点垂向和水平振动固有频率的近视和数值计算方法，并没有给出船舶强迫振动的计算方法。近年来，国内外相关学者给出了复杂组合壳结构振动响应的计算方法，比如文献(两端圆板封闭圆柱壳自由振动的半解析解.2012，16(11)：1306-1313.)和(基于区域分解的环肋圆柱壳-圆锥壳组合结构振动分析2013，30(1)：166-172.)。但是，这些方法仅适用于较简单结构的组合，同时也没有结构周围流体介质的流固耦合的影响，因此上述方法不能很好的模拟船舶特别是水下航行器振动响应。关于包含球壳及舱壁等结构和流体介质在内复杂组合壳系统的计算和处理方式尚未见报道。
技术实现思路
针对上述缺陷或者不足，本专利技术提供了一种基于精细算法的复杂组合壳结构振动响应定量计算方法，能够定量计算由圆柱壳、圆锥壳、球壳、加强筋及舱壁等组成的复杂组合壳的振动响应，为船舶特别是水下航行器的振动综合评估提供依据。本专利技术解决其技术问题所采用的技术方案是：提供一种基于精细算法的流体介质中复杂组合壳结构振动响应定量计算方法，包括以下步骤：(a)获取流体介质中组合壳系统的结构参数和激励参数，所述结构参数包括组合壳系统中各个结构的几何尺寸、材料特性参数；激励参数包括激励力的幅值及分布位置；(b)根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；(c)根据组合壳系统两端边界条件和结构连续条件，建立流体介质-组合壳系统的动力学模型；所述流体介质-组合壳系统动力学模型为：式中，Ti和Pi(i＝2...n)为流体介质-组合壳系统中各配点的场传递矩阵和外激励，Zi(i＝1...n)为流体介质-组合壳系统中各配点处的状态向量，I为单位矩阵；(d)通过Moore-Penrose广义逆矩阵求解得到流体介质-组合壳系统的振动响应。本专利技术所述的方法中，所述步骤(b)中，在计算组合壳系统各配点间的场传递矩阵和点传递矩阵前，先判断各相邻配点间是否有加强筋或舱壁结构存在，若是，则计算各配点间的场传递矩阵时需考虑点传递矩阵，根据精细积分办法计算获得各配点间的传递矩阵。本专利技术所述的方法中，所述组合壳系统中各个结构包括圆柱壳、锥壳、球壳、加强筋及舱壁。本专利技术所述的方法中，所述步骤(c)中，利用圆柱壳、锥壳、球壳分别进行组合壳的建模；将环肋简化为加强筋，采用Timeshenko梁考虑三个方向平衡进行建模；将舱壁简化为圆形板进行建模。本专利技术还提供一种基于精细算法的流体介质中复杂组合壳结构振动响应定量计算系统，包括：原始数据获取模块，用于获取流体介质中组合壳系统的结构参数和激励参数，所述结构参数包括组合壳系统中各个结构的几何尺寸、材料特性参数；激励参数包括激励力的幅值及分布位置；传递矩阵及外力向量计算模块，用于根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；模型建立模块，用于根据组合壳系统两端边界条件和结构连续条件，建立流体介质-组合壳系统的动力学模型；所述流体介质-组合壳系统动力学模型为：式中，Ti和Pi(i＝2...n)为流体介质-组合壳系统中各配点的场传递矩阵和外激励，Zi(i＝1...n)为流体介质-组合壳系统中各配点处的状态向量，I为单位矩阵；振动响应计算模块，用于通过Moore-Penrose广义逆矩阵求解得到流体介质-组合壳系统的振动响应。本专利技术所述的系统中，还包括：判断模块，用于在计算组合壳系统各配点间的场传递矩阵和点传递矩阵前，先判断各相邻配点间是否有加强筋或舱壁结构存在，若是，则通过传递矩阵及外力向量计算模块计算各配点间的场传递矩阵时需考虑点传递矩阵，根据精细积分办法计算获得各配点间的传递矩阵。本专利技术所述的系统中，所述组合壳系统中各个结构包括圆柱壳、锥壳、球壳、加强筋及舱壁。本专利技术所述的系统中，所述模型建立模块建模时，利用圆柱壳、锥壳、球壳分别进行组合壳的建模；将环肋简化为加强筋，采用Timeshenko梁考虑三个方向平衡进行建模；将舱壁简化为圆形板进行建模。本专利技术产生的有益效果是：本专利技术基于精细算法的复杂组合壳结构振动响应定量计算方法运算简易，易于实现，动力学模型的建立是基于解析法进行的，求解是基于数值进行求解的，这种半数值半解析方法能够有效提高计算效率，不受计算频带所限，并避免了采用数值进行复杂组合壳结构振动响应计算时繁琐的建模、需获取满足计算需求庞大规模的质量阵、刚度阵等。，由于本专利技术中采用精细算法进行求解结构的传递关系，使得计算结果与现有结果相比更加准确。附图说明下面将结合附图及实施例对本专利技术作进一步说明，附图中：图1是本专利技术实施例基于精细算法的流体介质中复杂组合壳结构振动响应定量计算方法的流程图；图2是本专利技术实施例的计算流程图；图3是本专利技术一个实施例的复杂组合壳结构系统模型图；图4是本专利技术实施例计算组合壳结构时，不同的轴向、周向波数时的组合壳结构的振动响应对比图；图5是本专利技术实施例计算组合壳结构时，相同结构参数时，不同组合壳结构之间的振动响应对比图；图6是本专利技术实施例基于精细算法的流体介质中复杂组合壳结构振动响应定量计算系统示意图；图7是本专利技术另一实施例基于精细算法的流体介质中复杂组合壳结构振动响应定量计算系统示意图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。本专利技术实施例基于精细算法的流体介质中复杂组合壳结构振动响应定量计算方法如图1所示，包括以下步骤：步骤S1、获取流体介质中组合壳系统的结构参数和激励参数，所述结构参数包括组合壳系统中各个结构的几何尺寸、材料特性参数；激励参数包括激励力的幅值及分布位置；步骤S2、根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；步骤S3、根据系统两端边界条件和结构连续条件，考虑流固耦合的影响，对场传递矩阵进行组装，建立组合壳系统的动力学模型；所述流体介质-组合壳系统动力学模型为：式中，Ti和Pi(i＝2...n)为流体介质-组合本文档来自技高网...

【技术保护点】
一种基于精细算法的流体介质中复杂组合壳结构振动响应定量计算方法，其特征在于，包括以下步骤：(a)获取流体介质中组合壳系统的结构参数和激励参数，所述结构参数包括组合壳系统中各个结构的几何尺寸、材料特性参数；激励参数包括激励力的幅值及分布位置；(b)根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；(c)根据组合壳系统两端边界条件和结构连续条件，建立流体介质‑组合壳系统的动力学模型；所述流体介质‑组合壳系统动力学模型为：式中，Ti和Pi(i＝2...n)为流体介质‑组合壳系统中各配点的场传递矩阵和外激励，Zi(i＝1...n)为流体介质‑组合壳系统中各配点处的状态向量，I为单位矩阵；(d)通过Moore‑Penrose广义逆矩阵求解得到流体介质‑组合壳系统的振动响应。

【技术特征摘要】
1.一种基于精细算法的流体介质中复杂组合壳结构振动响应定量计算方法，其特征在于，包括以下步骤：(a)获取流体介质中组合壳系统的结构参数和激励参数，所述结构参数包括组合壳系统中各个结构的几何尺寸、材料特性参数；激励参数包括激励力的幅值及分布位置；(b)根据获取的结构参数和激励参数，在组合壳表面上沿母线方向进行配点，采用精细积分方法计算组合壳系统各配点间的场传递矩阵和点传递矩阵，并计算外力向量；(c)根据组合壳系统两端边界条件和结构连续条件，建立流体介质-组合壳系统的动力学模型；所述流体介质-组合壳系统动力学模型为：式中，Ti和Pi为流体介质-组合壳系统中各配点的场传递矩阵和外激励，其中，i＝2…n；Z(ξj)为流体介质-组合壳系统中各配点处的状态向量，其中j＝1…n，I为单位矩阵；(d)通过Moore-Penrose广义逆矩阵求解得到流体介质-组合壳系统的振动响应；其中，步骤(b)中，在计算模型组合壳外表面进行配点，沿母线方向人为离散为n个分段，在范围内进行离散，其中离散点沿母线方向的坐标位置，依次表示为：ξ0ξ1ξ2…ξn-1ξn；对该分段两端状态向量的传递关系如下：式中，状态向量Z(ξi)为结构的相应位移和内力分量，U(ξ)为与ξ有关的系数矩阵，f(τ)为组合壳结构受到的外界激励力向量。2.根据权利要求1所述的方法，其特征在于，所述步骤(b)中，在计算组合壳系统各配点间的场传递矩阵和点传递矩阵前，先判断各相邻配点间是否有加强筋或舱壁结构存在，若是，则计算各配点间的场传递矩阵时需考虑点传递矩阵，根据精细积分办法计算获得各配点间的传递矩阵。3.根据权利要求1所述的方法，其特征在于，所述组合壳系统中各个结构包括圆柱壳、锥壳、球壳、加强筋及舱壁。4.根据权利要求3所述的方法，其特征在于，所述步骤(c)中，利用圆柱壳、锥壳、球壳分别进行组合壳的建模；将环肋简化为加强筋，采用Timeshenko梁考虑三个方向平衡进行建模；将舱壁简化为圆形板进行建模。5.一种基于精细算法的流体介质中复杂组合壳结构振动响应定量计算...

【专利技术属性】
技术研发人员：王献忠，吴卫国，漆琼芳，许瑞阳，周雍，马丽，
申请(专利权)人：武汉理工大学，
类型：发明
国别省市：湖北;42