A combined finite element-stiffness equation transfer method for steady state vibration response analysis of structures

An extended finite element transfer matrix method, in combination with stiffness equation transfer, is applied to dynamic response analysis of the structures under periodic excitations. In the present method, the transfer of state vectors from left to right in a combined finite element-transfer matrix (FE-TM) method is changed into the transfer of general stiffness equations of every section from left to right. This method has the advantages of reducing the order of standard transfer equation systems, and minimizing the propagation of round-off errors occurring in recursive multiplication of transfer and point matrices. Furthermore, the drawback that in the ordinary FE-TM method, the number of degrees of freedom on the left boundary be the same on the right boundary, is now avoided. A FESET program based on this method using microcomputers is developed. Finally, numerical examples are presented to demonstrate the accuracy as well as the potential of the proposed method for steady state vibration response analysis of structures.     

Vibration analysis of liquid-filled pipelines with elastic constraints

In this paper, a transfer matrix method (TMM) in frequency domain considering fluid-structure interaction of liquid-filled pipelines with elastic constraints is proposed. The time-domain equations considering fluid-structure interaction, are transformed into frequency domain by Laplace transformation, and then twelve fourth-order ordinary differential equations and two second-order ordinary differential equations are deduced from the frequency-domain equations. The results of the fourteen frequency-domain equations are assembled into a transfer matrix, which represents the motion of a single pipe section. Combined with point matrices that describe specified boundary conditions, an overall transfer matrix for liquid-filled pipeline system can be assembled. Using the method, all the pipeline with no and rigid constraints can be easily calculated by simply setting the stiffness of the restraining springs from zero to a large number. Taking into account the longitudinal vibration, transverse vibration and torsional vibration, the proposed method can be used to analyze the pipelines with bends. Several numerical examples with different constraints are presented here to illustrate the application of the proposed method. The results are validated by measured and simulation data. Through the numerical examples, it is shown that the proposed method is efficient.     

Coupled flexural-torsional vibration band gap in periodic beam including warping effect

The propagation of coupled flexural-torsional vibration in the periodic beam including warping effect is investigated with the transfer matrix theory. The band structures of the periodic beam, both including warping effect and ignoring warping effect, are obtained. The frequency response function of the finite periodic beams is simulated with finite element method, which shows large vibration attenuation in the frequency range of the gap as expected. The effect of warping stiffness on the band structure is studied and it is concluded that substantial error can be produced in high frequency range if the effect is ignored. The result including warping effect agrees quite well with the simulated result.     

Position optimization of particular solution sources for distributed source boundary point method by volume velocity-matching

Choosing particular solution source and its position have great influence on accuracy of sound field prediction in distributed source boundary point method.An optimization method for determining the position of particular solution sources is proposed to get high accuracy prediction result.In this method,tripole is chosen as the particular solution.The upper limit frequency of calculation is predicted by setting 1%volume velocity relative error limit using vibration velocity of structure surface.Then,the optimal position of particular solution sources,in which the relative error of volume velocity is minimum,is determined within the range of upper limit frequency by searching algorithm using volume velocity matching.The transfer matrix between pressure and surface volume velocity is constructed in the optimal position.After that,the sound radiation of structure is calculated by the matrix.The results of numerical simulation show that the calculation error is significantly reduced by the proposed method.When there are vibration velocity measurement errors,the calculation errors can be controlled within 5% by the method.     

An extended transfer matrix-finite element method for free vibration of plates

A combination of extended transfer matrix and finite element methods is proposed for obtaining vibration frequencies of structures. This method yields the value of the frequency once a trial value is assumed. By using this technique, the number of nodes required in the regular finite element method is reduced and therefore a smaller computer can be used. Besides, no plotting of the values of the determinants corresponding to each assumed frequency is necessary. A worked example is given for the case of vibration of a cantilever plate. The results show fast convergence from the assumed value to the true natural frequency.     

Stable reformulation of transfer matrix method for wave propagation in layered anisotropic media.

The numerical instability problem in the standard transfer matrix method has been resolved by introducing the layer stiffness matrix and using an efficient recursive algorithm to calculate the global stiffness matrix for an arbitrary anisotropic layered structure. For general anisotropy the computational algorithm is formulated in matrix form. In the plane of symmetry of an orthotropic layer the layer stiffness matrix is represented analytically. It is shown that the elements of the stiffness matrix are as simple as those of the transfer matrix and only six of them are independent. Reflection and transmission coefficients for layered media bounded by liquid or solid semi-spaces are formulated as functions of the total stiffness matrix elements. It has been demonstrated that this algorithm is unconditionally stable and more efficient than the standard transfer matrix method. The stiffness matrix formulation is convenient in satisfying boundary conditions for different layered media cases and in obtaining modal solutions. Based on this method characteristic equations for Lamb and surface waves in multilayered orthotropic media have been obtained. Due to the stability of the stiffness matrix method, the solutions of the characteristic equations are numerically stable and efficient. Numerical examples are given.     

An optimization of frame structures with exact dynamic constraints based on Timoshenko beam theory

Sensitivity analyses of eigensolutions and eigenfunctions of 3-D frame structures using the exact frequency equation from the transfer dynamic stiffness matrix that was derived on Timoshenko beam theory were developed in this paper. Based on the sensitivity data of frame structures, the minimum weight design with an exact frequency constraint can be carried out efficiently. Three examples that demonstrated the results obtained by the proposed method, are in good agreement with those computed by ANSYS.     

Coupled bending and torsional vibration of axially loaded Bernoulli-Euler beams including warping effects

The dynamic transfer matrix method for determining natural frequencies and mode shapes of the bending-torsion coupled vibration of axially loaded thin-walled beams with monosymmetrical cross sections is developed by using a general solution of the governing differential equations of motion based on Bernoulli-Euler beam theory. This method takes into account the effect of warping stiffness and gives allowance to the presence of axial force. The dynamic transfer matrix is derived in detail. Two illustrative examples on the application of the present theory are given for bending-torsion coupled beams with thin-walled open cross sections. The effects of axial load and warping stiffness on coupled bending-torsional frequencies are discussed. Compared with those available in the literature, numerical results demonstrate the accuracy and effectiveness of the proposed method.     

Free vibration of joined conical-cylindrical shells

An analysis is presented for the free vibration of joined conical-cylindrical shells. The governing equations of vibration of a conical shell, including a cylindrical shell as a special case, are written as a coupled set of first order differential equations by using the transfer matrix of the shell. Once the matrix has been determined, the entire structure matrix is obtained by the product of the transfer matrices of the shells and the point matrix at the joint, and the frequency equation is derived with terms of the elements of the structure matrix under the boundary conditions. The method has been applied to a joined truncated conical-cylindrical shell and an annular plate-cylindrical shell system, and the natural frequencies and the mode shapes of vibration calculated numerically. The results are presented.     

Dynamic stiffness method for inplane free vibration of rotating beams including Coriolis effects

The paper addresses the in-plane free vibration analysis of rotating beams using an exact dynamic stiffness method. The analysis includes the Coriolis effects in the free vibratory motion as well as the effects of an arbitrary hub radius and an outboard force. The investigation focuses on the formulation of the frequency dependent dynamic stiffness matrix to perform exact modal analysis of rotating beams or beam assemblies. The governing differential equations of motion, derived from Hamilton's principle, are solved using the Frobenius method. Natural boundary conditions resulting from the Hamiltonian formulation enable expressions for nodal forces to be obtained in terms of arbitrary constants. The dynamic stiffness matrix is developed by relating the amplitudes of the nodal forces to those of the corresponding responses, thereby eliminating the arbitrary constants. Then the natural frequencies and mode shapes follow from the application of the Wittrick–Williams algorithm. Numerical results for an individual rotating beam for cantilever boundary condition are given and some results are validated. The influences of Coriolis effects, rotational speed and hub radius on the natural frequencies and mode shapes are illustrated.     

体积速度匹配分布源边界点法的特解源位置优化方法研究

在分布源边界点法中,选择特解源及特解源位置对声场计算的精度具有很大的影响。为减小声场计算的误差,提出了一种确定特解源位置的优化方法。在该方法中,采用三极子作为特解源,以1%体积速度相对误差对结构声辐射计算的上限频率进行预测,然后在计算上限频率范围内通过体积速度匹配搜索,得到体积速度相对误差最小意义下的特解源优化位置,在该优化位置处构造声压与表面体积速度之间的传递矩阵来计算声场,以减小计算误差。数值仿真结果表明:在求取的计算上限频率范围内,特解源处于优化位置时,计算误差得以显著减小。当振速测量存在一定误差时,采用该方法仍可将计算误差控制在5%以内。      

Power flow boundary element analysis for multi-domain problems in vibrational built-up structures

A power flow boundary element methodology for the treatment of multi-domain problems is developed to predict the vibrational responses of coupled structures in medium-to-high frequency ranges. In the proposed methodology, a matrix equation is formulated by utilizing the power flow coupling relationships at the interface of the domains. This analysis method is successfully applied to the treatment of simply supported coupled beams and coupled plates. The vibration analysis of these structures is then considered for all wave components. The vibrational energy density and intensity of two coupled structures obtained through numerical simulations are compared with those of the original power flow analysis, and these results show good agreement.     

The design of resonators with specified reactances

The transfer matrix method is used to analyse the vibration of a thin straight rod made up of a number of linearly tapered cones. This analysis is used to design a bar with a specified end reactance in a certain range of frequency. Results are given for three specified reactances and it is shown how they may be used in the design of a contact impedance meter and a stiffness measuring device.     

Stop-pass behavior of acoustic waves in a 1D fractured system

This study examines the dispersion and the stop-pass band behavior of acoustic waves propagating across periodically spaced and non periodically spaced parallel fractures. Laboratory ultrasonic wave measurements performed on a stack of synthetic fractures (identical steel plates with roughened interfaces) and numerical propagator matrix simulations show spectra with distinct stop-pass band structures that develop with decreasing fracture stiffness. To understand the physics behind these observations, an exact dispersion equation for wave propagation through an infinite series of equally spaced fractures is derived using displacement-discontinuity boundary conditions to model the constitutive behavior of the fractures and Floquet's (Bloch's) theory for the periodic boundary conditions. Both the measured and numerically simulated stop-pass band structures show good agreement with the theoretical predictions. Furthermore, the theory reveals that the left boundary of the stop-bands contains information about the fracture stiffness, suggesting the possibility of determining the stiffness of the parallel fractures from seismic waves. This paper also discusses the effects of fractured systems with random distributions of fracture spacings and stiffnesses on the stop-pass band structures of seismic waves in fractured rock.     

THE DYNAMIC STIFFNESS MATRIX METHOD IN FORCED VIBRATION ANALYSIS OF MULTIPLE-CRACKED BEAM

The dynamic behaviour of a beam with numerous transverse cracks is studied. Based on the equivalent rotational spring model of crack and the transfer matrix for beam, the dynamic stiffness matrix method has been developed for spectral analysis of forced vibration of a multiple cracked beam. As a particular case, when the excitation frequency is close to zero, the solution for static response of beam with an arbitrary number of cracks has been obtained exactly in an analytical form. In general case, the effect of crack number and depth on the dynamic response of beam was analyzed numerically.     

全新的电导率特征矩阵方法及其在石墨烯THz频率光学特性上的应用

基于麦克斯韦方程组所要求的电磁场边界条件首次从理论上严格推导得到超薄导电体及其复合多层介质结构光学特性的一般计算方法及其特征矩阵公式, 其优点在于只要借助于导电体的电导率而无需知道其介电常数和磁导率即可计算得到反射、透射和吸收等光学特性, 克服了传统的传输矩阵方法必需知道组成材料的介电常数和磁导率才能获得其光学性质的问题, 并利用此方法获得了石墨烯及其复合多层结构在THz频率范围内反射、透射和吸收等光学行为.     

Dynamic stiffness method for exact inplane free vibration analysis of plates and plate assemblies

The inplane free vibration behaviour of plates is investigated using the dynamic stiffness method. Some distinctive modes which went unnoticed in earlier investigations using the dynamic stiffness method have been addressed by revisiting the problem and focusing on the special set of missing solutions. Results are validated extensively both by published results as well as by numerical studies using NASTRAN and ABAQUS. The accuracy of the finite element method for inplane free vibration analysis is assessed and critically examined through the provision of benchmark solutions. Some representative modes that are missed by well-established dynamic-stiffness-based computer programs are presented. The inplane dynamic stiffness matrix presented is of great importance when combined with the out of plane matrix in order to obtain the closed-form solution for free vibration analysis of structures with complex geometries.     

PREDICTION AND MEASUREMENT OF SOME DYNAMIC PROPERTIES OF SANDWICH STRUCTURES WITH HONEYCOMB AND FOAM CORES

Some dynamical properties of sandwich beams and plates are discussed. The types of elements investigated are three-layered structures with lightweight honeycomb or foam cores with thin laminates bonded to each side of the core. A six order differential equation governing the apparent bending of sandwich beams is derived using Hamilton's principle. Bending, shear and rotation are considered. Boundary conditions for free, clamped and simply supported beams are formulated. The apparent bending stiffness of sandwich beams is found to depend on the frequency and the boundary conditions for the structure. Simple measurements on sandwich beams are used to determine the bending stiffness of the entire structure and at the same time the bending stiffness of the laminates as well as the shear stiffness of the core. A method for the prediction of eigenfrequencies and modes of vibration are presented. Eigenfrequencies for rectangular and orthotropic sandwich plates are calculated using the Rayleigh-Ritz technique assuming frequency dependent material parameters. Predicted and measured results are compared.     

Dynamic stiffness matrix of thin-walled composite I-beam with symmetric and arbitrary laminations

For the spatially coupled free vibration analysis of thin-walled composite I-beam with symmetric and arbitrary laminations, the exact dynamic stiffness matrix based on the solution of the simultaneous ordinary differential equations is presented. For this, a general theory for the vibration analysis of composite beam with arbitrary lamination including the restrained warping torsion is developed by introducing Vlasov's assumption. Next, the equations of motion and force–displacement relationships are derived from the energy principle and the first order of transformed simultaneous differential equations are constructed by using the displacement state vector consisting of 14 displacement parameters. Then explicit expressions for displacement parameters are derived and the exact dynamic stiffness matrix is determined using force–displacement relationships. In addition, the finite-element (FE) procedure based on Hermitian interpolation polynomials is developed. To verify the validity and the accuracy of this study, the numerical solutions are presented and compared with analytical solutions, the results from available references and the FE analysis using the thin-walled Hermitian beam elements. Particular emphasis is given in showing the phenomenon of vibrational mode change, the effects of increase of the modulus and the bending–twisting coupling stiffness for beams with various boundary conditions.     

A distributed-mass approach for dynamic analysis of Bernoulli-Euler plane frames

The exact dynamic analysis of plane frames should consider the effect of mass distribution in beam elements, which can be achieved by using the dynamic stiffness method. Solving for the natural frequencies and mode shapes from the dynamic stiffness matrix is a nonlinear eigenproblem. The Wittrick-Williams algorithm is a reliable tool to identify the natural frequencies. A deflated matrix method to determine the mode shapes is presented. The dynamic stiffness matrix may create some null modes in which the joints of beam elements have null deformation. Adding an interior node at the middle of beam elements can eliminate the null modes of flexural vibration, but does not eliminate the null modes of axial vibration. A force equilibrium approach to solve for the null modes of axial vibration is presented. Orthogonal conditions of vibration modes in the Bernoulli-Euler plane frames, which are required in solving the transient response, are theoretically derived. The decoupling process for the vibration modes of the same natural frequency is also presented.