充液圆柱壳固有频率预报方法技术

本发明专利技术的目的在于提供充液圆柱壳固有频率预报方法，包括以下步骤：采用Kennard薄壳理论建立圆柱壳的自由振动方程；通过求解圆柱坐标系下的Helmholtz波动方程，将流体作为径向载荷作用于圆柱壳上，建立充液圆柱壳动力学模型；选取改进状态向量；由Kennard薄壳理论得到传统状态向量的元素位移、转角、力等与中面位移的关系，从而建立改进状态向量与传统状态向量之间的转换矩阵。引入充液圆柱壳两端边界条件预报充液圆柱壳固有频率。本发明专利技术降低了一阶导数方程和传递矩阵推导工作量及推导错误出现的可能性，适用范围广，实现过程简便，非常有利于编程计算；预报精度高，整个求解过程不会导致变量自由度数的增加，保证了较高的预报效率。

Prediction method of natural frequency of liquid filled cylindrical shell

【技术实现步骤摘要】
充液圆柱壳固有频率预报方法
本专利技术涉及的是一种固有频率预测方法。
技术介绍
充液圆柱壳结构在建筑及船舶管路系统中普遍存在，对其固有频率的预报也存在较多的方法。传递矩阵方法因其具有涉及变量自由度少、求解简便快捷等优点，成为求解链式结构动力学问题的常用方法，在求解梁或管结构动力学、管内声传播及管路系统流固耦合等问题时有较多应用。若将结构控制微分方程写成一阶微分方程组的形式，则起点与终点的状态向量间即可用传递矩阵建立简单关系，利用边界条件即可进行充液圆柱壳固有频率预报。因此，传递矩阵法计算精度主要取决于状态向量的选取与矩阵方程的积分。动力方程积分方法研究较多，如精细时程积分法、龙格库塔法、齐次扩容积分法。随着计算机技术的发展，积分计算精度不断提高，已基本能满足要求。而随着传递矩阵法在多自由度系统中的应用，状态向量一阶方程的推导却越加困难。传统状态向量一般由位移、速度、力、力矩等组成，而状态向量并不能从振动方程中直接获得，造成状态向量一阶微分方程推导过程非常复杂，不仅工作量增加，且常有错误出现。罗文(罗文，输流管路流固耦合计算研究.哈尔滨工程大学硕士学位论文.2014.)应用传递矩阵法求解了充液圆柱壳的无量纲固有频率，但其状态向量一阶微分方程推导过程较为复杂，且其传递矩阵为周向传递矩阵，与一般传递矩阵沿链式结构长度方向不同。万浩川等人(万浩川，李以农，郑玲，改进的结构振动传递矩阵法，振动与冲击，2013,32(1):173–177.)提出了一种使用改进传递矩阵法计算圆柱壳振动固有频率的方法，他们提出的解法较传统解法推导与求解过程得到简化，但并未考虑内部充液对圆柱壳改进传递矩阵和固有频率的影响。
技术实现思路
本专利技术的目的在于提供对传统传递矩阵法进行改进的充液圆柱壳固有频率预报方法。本专利技术的目的是这样实现的：本专利技术充液圆柱壳固有频率预报方法，其特征是：(1)充液圆柱壳动力学建模：首先采用Kennard薄壳理论建立圆柱壳的自由振动方程同时将流体作为径向载荷作用于圆柱壳上；之后通过求解圆柱坐标系下的Helmholtz波动方程，根据耦合面上流体、固体的边界条件最终推导得到充液圆柱壳动力学模型：考虑壳体惯性力的影响并且有流体载荷作用，Kennard薄壳理论振动微分方程形式如下：其中为与圆柱壳体材料相同的平板中伸缩波的传播相速度；p为管内流体作用在圆柱壳内壁上的水动压力；对轴向半波数为m、周向波数为n的振动模态，ωmn为频率，上式的解写为：圆柱壳内部流体方程满足：假设内部流体声压的解写为：为径向波数，满足下式：其中：kf＝ω/cf，cf为管内流体自由声场波速，km表示轴向波数，与系统的边界条件有关。在r＝R处，连续性条件为：将式w(x，θ，t)和式p(θ，r，x，t)代入连续性条件得：(2)改进状态向量选取：将式(2)u(x，θ，t)、v(x，θ，t)、w(x，θ，t)和式p(x，θ，t)|r＝R代入Kennard薄壳理论振动微分方程得：直接由上式选择改进状态向量:利用改进状态向量改写：C为8阶方阵：C(1,2)＝C(3,4)＝C(5,6)＝C(6,7)＝C(7,8)＝1；(3)引入转换矩阵：壳模型的状态向量为：Kennard薄壳理论给出了传统状态向量的元素位移、转角、力与中面位移的关系，从而得到改进状态向量与传统状态向量满足：Φ＝DξD为转换矩阵，其非0元素满足：D(1，1)＝D(2，3)＝D(3，5)＝D(4，6)＝1；(4)固有频率预报：由式及转换矩阵可以求出充液圆柱壳传递矩阵T：Φ(x)＝Dξ(χ)＝DeCsξ(0)＝DeCsD-1Dξ(0)＝DeCsD-1Φ(0)＝TΦ(0)引入充液圆柱壳两端边界条件对充液圆柱壳固有频率进行预报；其中，x，r,θ为柱坐标方向；μ为泊松比；ω为频率变量、R为圆柱壳中面半径、E为杨氏模量、n为周向波数、km为轴向波数、Jn为贝塞尔函数、Nx为中面等效薄膜力、Mx为中面等效弯矩、u，v，w为沿柱坐标方向的位移、ρ为密度、h为圆柱壳厚度、m为轴向半波数、为径向波数、t为时间变量、cf为流体自由声场波速、J′n为贝塞尔函数的一阶导数、Fx为横向合剪力、Sx为薄膜合剪力。本专利技术还可以包括：1、所述的固有频率预报，两端简支边界为：有：令:有非零解条件为T′＝0，而T′为频率函数由此可得各阶模态固有频率值；对不同的边界约束，可得不同的T′，但均有T′＝0。本专利技术的优势在于：本专利技术降低了一阶导数方程和传递矩阵推导工作量及推导错误出现的可能性，适用范围广，实现过程简便，非常有利于编程计算；预报精度高，整个求解过程不会导致变量自由度数的增加，保证了较高的预报效率。附图说明图1为本专利技术的流程图；图2为基于薄壳模型的充液圆柱壳结构示意图。具体实施方式下面结合附图举例对本专利技术做更详细地描述：结合图1-2，本专利技术充液圆柱壳固有频率预报的改进传递矩阵法，是基于传统传递矩阵方法和改进状态向量实现的。本专利技术的预报方法含有以下步骤：(1)充液圆柱壳动力学建模。首先采用Kennard薄壳理论建立圆柱壳的自由振动方程同时将流体作为径向载荷作用于圆柱壳上；之后通过求解圆柱坐标系下的Helmholtz波动方程，根据耦合面上流体、固体的边界条件最终推导得到充液圆柱壳动力学模型。考虑壳体惯性力的影响并且有流体载荷作用时，Kennard薄壳理论振动微分方程形式如下(其中径向运动以远离壳体轴线为正)：其中β2＝h2/12R2，为与圆柱壳体材料相同的平板中伸缩波的传播相速度；p为管内流体作用在圆柱壳内壁上的水动压力，反映了管内流体的影响。对轴向半波数为m、周向波数为n的振动模态，ωmn为频率，式(1)的解可写为：圆柱壳内部流体方程满足：假设内部流体声压的解可写为：为径向波数，满足下式：其中：kf＝ω/cf，cf为管内流体自由声场波速，km表示轴向波数，与系统的边界条件有关。连续性条件为(在r＝R处)：将式(2c)和式(4)代入连续性条件得：(2)改进状态向量选取将式(2)和式(7)代入式(1)得：直接由式(8)选择改进状态向量:...

【技术保护点】
1.充液圆柱壳固有频率预报方法，其特征是：/n(1)充液圆柱壳动力学建模：/n首先采用Kennard薄壳理论建立圆柱壳的自由振动方程同时将流体作为径向载荷作用于圆柱壳上；之后通过求解圆柱坐标系下的Helmholtz波动方程，根据耦合面上流体、固体的边界条件最终推导得到充液圆柱壳动力学模型：/n考虑壳体惯性力的影响并且有流体载荷作用，Kennard薄壳理论振动微分方程形式如下：/n

【技术特征摘要】
1.充液圆柱壳固有频率预报方法，其特征是：
(1)充液圆柱壳动力学建模：
首先采用Kennard薄壳理论建立圆柱壳的自由振动方程同时将流体作为径向载荷作用于圆柱壳上；之后通过求解圆柱坐标系下的Helmholtz波动方程，根据耦合面上流体、固体的边界条件最终推导得到充液圆柱壳动力学模型：
考虑壳体惯性力的影响并且有流体载荷作用，Kennard薄壳理论振动微分方程形式如下：









其中为与圆柱壳体材料相同的平板中伸缩波的传播相速度；p为管内流体作用在圆柱壳内壁上的水动压力；
对轴向半波数为m、周向波数为n的振动模态，ωmn为频率，上式的解写为：









圆柱壳内部流体方程满足：



假设内部流体声压的解写为：




为径向波数，满足下式：



其中：kf＝ω/cf，cf为管内流体自由声场波速，km表示轴向波数，与系统的边界条件有关。
在r＝R处，连续性条件为：



将式w(x，θ，t)和式p(θ，r，x，t)代入连续性条件得：



(2)改进状态向量选取：
将式(2)u(x，θ，t)、v(x，θ，t)、w(x，θ，t)和式p(x，θ，t)|r＝R代入Kennard薄壳理论振动微分方程得：









直接由上式选择改进状态向量:



利用改进状态向量改写：



C为8阶方阵：
C(1，2)＝C(3，4)＝C(5，...

【专利技术属性】
技术研发人员：柳贡民，曹银行，张文平，张新玉，明平剑，曹贻鹏，国杰，赵晓臣，
申请(专利权)人：哈尔滨工程大学，
类型：发明
国别省市：黑龙;23