随机Bagley-Torvik方程的非平稳解析解
Non-stationary analytic solution of the stochastic Bagley‑Torvik equation
-
摘要: Bagley‑Torvik(B‑T)方程是一种带分数(3/2)阶导数项的运动微分方程,被应用于描述刚性板在牛顿流体中的振动状态。发展非齐次项为随机过程的B‑T方程非平稳解析解,将B‑T方程转化为矩阵形式的半阶状态空间方程并进行特征分析,得到复特征值和特征向量;引入广义坐标变换将方程解耦为独立的1/2阶微分方程组,并利用Laplace变换求解得到广义坐标下的解;将广义坐标解转换为自然坐标解,得到脉冲或阶跃响应函数的解析解。方程非齐次项为随机过程时,可利用Laplace变换求解时变频响函数,并基于激励与响应功率谱密度之间的关系得到非平稳随机响应解析解。以Spanos‑Solomos完全非平稳随机激励为例,通过数值案例验证方法的正确性。Abstract: The Bagley-Torvik(B-T) equation is a differential equation of motion with fractional (3/2)-order derivative terms that is applied to describe the motion of a rigid plate in Newtonian, viscous fluid. In this paper, we develop non-stationary analytic solutions of the B-T equation whose inhomogeneous term is a stochastic process. The B-T equation is transformed into a half-order state-space equation in matrix form and eigen-analysis is performed to obtain complex eigenvalues and eigenvectors. Subsequently, the generalized coordinate transformation is introduced to decouple the equation into a system of independent 1/2-order differential equations which are solved by Laplace transform to obtain the solution in generalized coordinates; The generalized coordinate solution is converted into a natural coordinate solution to obtain the impulse or step response function. When the inhomogeneous term of the equation is a stochastic process, the Laplace transform can be used to derive the time-varying frequency response function from which the analytical solution of the non-stationary stochastic response can be obtained by relying on the relationship between the excitation and the response power spectral density. The correctness of the method is verified by numerical cases using the Spanos-Solomos fully non-statoionary stochastic excitation as an example.