Dynamic modeling for rotating composite Timoshenko beam and analysis on its bending-torsion coupled vibration

A formulation is presented for steady-state dynamic responses of rotating bending-torsion coupled composite Timoshenko beams (CTBs) subjected to distributed and/or concentrated harmonic loadings. The separation of cross section's mass center from its shear center and the introduced coupled rigidity of composite material lead to the bending-torsion coupled vibration of the beams. Considering those two coupling factors and based on Hamilton's principle, three partial differential non-homogeneous governing equations of vibration with arbitrary boundary conditions are formulated in terms of the flexural translation, torsional rotation and angle rotation of cross section of the beams. The parameters for the damping, axial load, shear deformation, rotation speed, hub radius and so forth are incorporated into those equations of motion. Subsequently, the Green's function element method (GFEM) is developed to solve these equations in matrix form, and the analytical Green's functions of the beams are given in terms of piecewise functions. Using the superposition principle, the explicit expressions of dynamic responses of the beams under various harmonic loadings are obtained. The present solving procedure for Timoshenko beams can be degenerated to deal with for Rayleigh and Euler beams by specifying the values of shear rigidity and rotational inertia. Cantilevers with bending-torsion coupled vibration are given as examples to verify the present theory and to illustrate the use of the present formulation. The influences of rotation speed, bending-torsion couplings and damping on the natural frequencies and/or shape functions of the beams are performed. The steady-state responses of the beam subjected to external harmonic excitation are given through numerical simulations. Remarkably, the symmetric property of the Green's functions is maintained for rotating bending-torsion coupled CTBs, but there will be a slight deviation in the numerical calculations.     

Interaction of surface water waves with a vertical elastic plate: a hypersingular integral equation approach

In this paper, we present an alternative method to investigate scattering of water waves by a submerged thin vertical elastic plate in the context of linear theory. The plate is submerged either in deep water or in the water of uniform finite depth. Using the condition on the plate, together with the end conditions, the derivative of the velocity potential in the direction of normal to the plate is expressed in terms of a Green’s function. This expression is compared with that obtained by employing Green’s integral theorem to the scattered velocity potential and the Green’s function for the fluid region. This produces a hypersingular integral equation of the first kind in the difference in potential across the plate. The reflection coefficients are computed using the solution of the hypersingular integral equation. We find good agreement when the results for these quantities are compared with those for a vertical elastic plate and submerged and partially immersed rigid plates. New results for the hydrodynamic force on the plate, the shear stress and the shear strain of the vertical elastic plate are also evaluated and represented graphically.     

Steady flow of incompressible fluid between two co-rotating disks

This article provides an analytical solution of the Navier–Stokes equations for the case of the steady flow of an incompressible fluid between two uniformly co-rotating disks. The solution is derived from the asymptotical evolution of unknown components of velocity and pressure in a radial direction – in contrast to the Briter–Pohlhausen analytical solution, which is supported by simplified Navier–Stokes equations. The obtained infinite system of ordinary differential equations forms recurrent relations from which unknown functions can be calculated successively. The first and second approximations of solution are solved analytically and the third and fourth approximations of solutions are solved numerically. The numerical example demonstrates agreements with results obtained by other authors using different methods.     

Numerical method in wave-body interactions

The application of Green's function in calculation of flow characteristics around submerged and floating bodies due to a regular wave is presented. It is assumed that the fluid is homogeneous, inviscid and incompressible, the flow is irrotational and all body motions are small. Two methods based on the boundary integral equation method (BIEM) are applied to solve associated problems. The first is a low order panel method with triangular flat patches and uniform distribution of velocity potential on each panel. The second method is a high order panel method in which the kernels of the integral equations are modified to make it nonsingular and amenable to solution by the Gaussian quadrature formula. The calculations are performed on a submerged sphere and some floating spheroids of different aspect ratios. The excellent level of agreement with the analytical solutions shows that the second method is more accurate and reliable.     

The rolling ship problem—revisited

Solution of dual integral equations involving trigonometric functions as kernel has been utilized here to reinvestigate the classical rolling ship problem which involves study of the wave motion due to small rolling oscillations of a thin vertical plate partially immersed in deep water. Well-known results are produced in a simple and straightforward manner. The analytical solution for the velocity potential is depicted graphically.     

Green's functions for forced vibration analysis of bending-torsion coupled Timoshenko beam

A method based on Green's functions is proposed for the analysis of the steady-state dynamic response of bending-torsion coupled Timoshenko beam subjected to distributed and/or concentrated loadings. Damping effects on the bending and torsional directions are taken into account in the vibration equations. The elastic boundary conditions with bending-torsion coupling and damping effects are derived and the classical boundary conditions can be obtained by setting the values of specific stiffness parameters of the artificial springs. The Laplace transform technology is employed to work out the Green's functions for the beam with arbitrary boundary conditions. The Green's functions are obtained for the beam subject to external lateral force and external torque, respectively. Coupling effects between bending and torsional vibrations of the beam can be studied conveniently through these analytical Green's functions. The direct expressions of the steady-state responses with various loadings are obtained by using the superposition principle. The present Green's functions for the Timoshenko beam can be reduced to those for Euler–Bernoulli beam by setting the values of shear rigidity and rotational inertia. In order to demonstrate the validity of the Green's functions proposed, results obtained for special cases are given for a comparison with those given in the literature and they agree with each other exactly. The influences of external loading frequency and eccentricity on Green's functions of bending-torsion coupled Timoshenko beam are investigated in terms of the numerical results for both simply supported and cantilever beams. Moreover, the symmetric property of the Green's functions and the damping effects on the amplitude of Green's functions of the beam are discussed particularly.     

部分潜入水中椭圆柱体振动分析*

本文在文献[1]的基础上,提出一种研究部分潜入水中椭圆柱体弯曲振动的解析法,并指出文献[2]的不足之处.作为特例,本文计算了考虑水体可压缩影响时的水中圆柱体的自振频率,给出了可压缩性影响范围.     

计及液面波动和液体可压缩性时部分潜入水中椭圆柱体的扭转振动

本文研究部分潜入水中椭圆柱体的扭转振动问题,同时计及了液面波动和液体可压缩性对椭圆柱体扭转振动的影响。利用马休(Mathieu)级数以及一组广义三角级数的正交完备性,导出了柱水耦联扭振振型函数和频率方程的精确解析解,可利用非线性代数方程的求根方法,数值解出各阶频率参数。     

Prediction of the dynamic response of a mini-cantilever beam partially submerged in viscous media using finite element method

In this work, the use of mini cantilever beams for characterization of rheological properties of viscous materials is demonstrated. The dynamic response of a mini cantilever beam partially submerged in air and water is measured experimentally using a duel channel PolyTec scanning vibrometer. The changes in dynamic response of the beam such as resonant frequency, and frequency amplitude are compared as functions of the rheological properties (density and viscosity) of fluid media. Next, finite element analysis (FEA) method is adopted to predict the dynamic response of the same cantilever beam. The numerical prediction is then compared with experimental results already performed to validate the FEA modeling scheme. Once the model is validated, further numerical analysis was conducted to investigate the variation in vibration response with changing fluid properties. Results obtained from this parametric study can be used to measure the rheological properties of any unknown viscous fluid.     

Computation of added mass and damping coefficients due to a heaving cylinder

We present the boundary value problem (BVP) for the heave motion due to a vertical circular cylinder in water of finite depth. The BVP is presented in terms of velocity potential function. The velocity potential is obtained by considering two regions, namely, interior region and exterior region. The solutions for these two regions are obtained by the method of separation of variables. The analytical expressions for the hydrodynamic coefficients are derived. Computational results are presented for various depth to radius and draft to radius ratios.     

Effect of Torsion and Warping on the Free Vibration of Sandwich Beams

The problem on free vibrations of wide sandwich beams is tackled in this paper. Torsional and warping effects in addition to flexure are included in the formulation of the dynamic problem. In order to show the effects of bending-torsion coupling and warping on the natural frequencies and the corresponding vibration modes, three cases are considered. First, the warping and torsion effects are neglected, second, the effect of warping on the coupling terms is neglected, and third, the effect of warping on the coupling terms is included. The viscoelastic core is modeled by elastic translational and rotational springs. The finite-difference method is used to solve the partial differential equations of motion with different boundary conditions for the top and bottom layers. Results for different materials, fiber orientations, depth-to-width ratios, and boundary conditions are found. The natural frequencies and the corresponding vibration modes obtained are in a good agreement with those cited in the literature. If the bending-torsion coupling is pronounced, the inclusion of warping affects the natural frequency considerably.__________Russian translation published in Mekhanika Kompozitnykh Materialov, Vol. 41, No. 2, pp. 163–176, April–May, 2005.     

积分形式非局部本构关系的界带分析方法

基于Hamilton体系研究了Eringen的非局部线弹性本构关系.Eringen的非局部线弹性理论存在积分型和微分型两类本构关系.由于方程的形式简单,目前多采用微分型本构;而积分型本构方程是典型的积分-微分方程,数值求解较为困难.在分析结构力学中提出的界带分析方法,成功求解了时间滞后问题的积分-微分方程.根据分析动力学与分析结构力学的模拟关系,将界带分析方法引入到非局部理论的积分型本构方程,可以实现积分-微分方程的数值求解.通过杆件的振动分析算例验证了该套理论算法的准确性和可行性,也指出了辛体系算法在非局部力学问题中的潜力.     

Nonlinear characteristics of vortex-induced vibration at low Reynolds number

Numerical simulations of vortex-induced vibration of a two-dimensional elastic circular cylinder under the uniform flow are calculated when Reynolds number is 200. In order to achieve the vortex-induced vibration, two-dimensional incompressible Navier–Stokes equations are solved with the space–time finite element method, the equations of the cylinder motion are solved with the new explicit integral method and the emeshing is achieved by the spring analogy technology. Considering vortex-induced vibration with the low reduced damping parameters, the variety trends of the lift coefficient, the drag coefficient, the displacement of cylinder are analyzed under different oscillating frequencies of cylinder. The nonlinear phenomena of locked-in, beat and phaseswith are captured successfully. The limit cycle and bifurcation of lift coefficient and displacement are analyzed. Besides, the Poincare sections of the lift coefficient are used for discussing the bifurcation of periodic solution. There are some differences in nonlinear characteristics between the results of the one degree of freedom cylinder model and those of the two degrees of freedom cylinder model. The streamwise vibration has a certain effect on the lateral vibration.     

Numerical solutions of two-dimensional nonlinear integral analysis

In this paper, the approximate solutions for two different type of two-dimensional nonlinear integral equations: two-dimensional nonlinear Volterra-Fredholm integral equations and the nonlinear mixed Volterra-Fredholm integral equations are obtained using the Laguerre wavelet method. To do this, these two-dimensional nonlinear integral equations are transformed into a system of nonlinear algebraic equations in matrix form. By solving these systems, unknown coefficients are obtained. Also, some theorems are proved for convergence analysis.Some numerical examples are presented and results are compared with the analytical solution to demonstrate the validity and applicability of the proposed method.     

An ultra-weak method for acoustic fluid–solid interaction

We introduce the ultra-weak variational formulation (UWVF) for fluid–solid vibration problems. In particular, we consider the scattering of time-harmonic acoustic pressure waves from solid, elastic objects. The problem is modeled using a coupled system of the Helmholtz and Navier equations. The transmission conditions on the fluid–solid interface are represented in an impedance-type form after which we can employ the well known ultra-weak formulations for the Helmholtz and Navier equations. The UWVF approximation for both equations is computed using a superposition of propagating plane waves. A condition number based criterion is used to define the plane wave basis dimension for each element. As a model problem we investigate the scattering of sound from an infinite elastic cylinder immersed in a fluid. A comparison of the UWVF approximation with the analytical solution shows that the method provides a means for solving wave problems on relatively coarse meshes. However, particular care is needed when the method is used for problems at frequencies near the resonance frequencies of the fluid–solid system.     

Stability analysis of the servohydraulic equations with linear state feedback and harmonic reference signal

The servohydraulic equations describe the dynamic behaviour of a hydrostatic drive featuring a servo‐valve and a hydraulic cylinder. Under the assumptions of constant supply pressures, a rigid support of the cylinder, and a rigid mass attached to the piston, the evolution of cylinder pressures as well as piston position and velocity is governed by a 4th order non‐linear ordinary differential equation. The servo valve is assumed to be much faster than the dynamics of the mass‐cylinder system. Thus, the valve dynamics is neglected. Adding a feedback control law makes up the servohydraulic equations. These equations are partially singularly perturbed. The perturbation parameter is associated with fluid compressibility and the reduced system corresponds to the case of an incompressible fluid. This paper deals with the stability analysis of periodic orbits arising in the case of harmonic reference position trajectories. Geometric singular perturbation theory is used to generate stability charts. Numerical computations verify the analytical results. For a given set of hydraulic system parameters, the stability of periodic orbits depends on controller gains as well as on the amplitude and frequency of the reference signal.     

上覆单相弹性层的饱和地基上刚性基础竖向振动的轴对称混合边值问题

运用作者提出的饱和土弹性波动方程,从理论上研究了上覆单相弹性土层的饱和地基上刚性基础的竖向振动轴对称问题,即采用Hankel积分变换技术并按混合边值条件建立了部分饱和地基上刚性基础竖向振动的对偶积分方程,并将其蜕化为完全饱和地基的情形;该对偶积分方程可化为易于数值计算的第二类Fredholm积分方程。文末的算例给出了地基表面动力柔度系数C_v随无量纲频率a₀的变化曲线。     

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam with time-dependent velocity

Non-linearly parametric resonances of an axially moving viscoelastic sandwich beam are investigated in this paper. The beam is moving with a time-dependent velocity, namely a harmonically varied velocity about the mean velocity. The partial differential equation is discretized into nonlinear ordinary differential equations via the method of Galerkin truncation and then the steady-state response is obtained using the method of multiple scales, an approximate analytical method. The tuning equations are obtained by eliminating secular terms and the amplitude of the vibration is derived from the tuning equations expressed in polar form, and two bifurcation points are obtained as well. Additionally, the stability conditions of trivial and nontrivial solutions are analyzed using the Routh–Hurwitz criterion. Eventually, the effects of various parameters such as the thickness of core layer, mean velocity, initial tension, and the amplitude of axially moving velocity on amplitude–frequency response curves and unstable regions are investigated.     

关于多裂纹圆柱体的扭转*

本文在文[1]基础上,导出了含有任意分布裂纹系的圆柱扭曲函数的解析表达式,从而把问题化为以未知位错密度函数表示的奇异积分方程组.文中利用奇异积分方程的数值方法[2,7],对带有多根裂纹的圆柱的抗扭刚度和应力强度因子作了若干数值计算.此外,本文还首次将裂纹切割法[5]推广用于求解矩形柱的扭转,数值结果表明方法是成功的.     

Stokes flow of a micropolar fluid exterior to several non-intersecting closed surfaces,but contained by an exterior contour

The problem of determining the Stokes flow of a micropolar fluid exterior to several closed surfaces but contained by an exterior contour that encloses all the interior surfaces, is formulated as a system of linear Fredholm integral equations of the second kind. These integral equations are obtained when the velocity and microrotation vector fields are represented by a double-layer potential with unknown density, and certain singular solutions of the Stokes' micropolar equations. This double-layer potential is defined over the union of all the surfaces involved including the exterior contour. The singularities, corresponding to a concentrated force and concentrated couple located within each interior surface, give rise to force and torque whose magnitudes are linearly dependent on the unknown density of the double layer. It is shown that the system possesses a unique continuous solution when the boundaries are Lyapunov surfaces and the boundary data is continuous.