弹性——粘弹性复合结构模态理论

本文研究弹性-粘弹性复合结构动力学基本问题.复合结构动力学方程是一组微分积分方程,引入增广状态变量.将其变换为常规的状态方程.研究了状态方程特征解的性质.提出了"振荡模态"和"蠕变模态"慨念.给出了脉冲响应矩阵和传递函数矩阵,讨论了它们的特性.复合结构模态理论为其动特性和动响应分析提供理论依据.     

水平成层介质中粘弹性波的计算

本文提出了一种水平成层介质中弹性波粘弹性波计算的新方法，分别研究了线源和点源作用情况。该方法适用于各种粘弹性模型，数值计算简单，可模拟任意震源及所产生的各种体波、面波，数值结果表明具有很高的计算精度和计算效率。     

EIGEN THEORY OF VISCOELASTIC MECHANICS FOR ANISOTROPIC SOLIDS

Anisotropic viscoelastic mechanics is studied under anisotropic subspace. It is proved that there also exist the eigen properties for viscoelastic medium. The modal Maxwell's equation, modal dynamical equation (or modal equilibrium equation) and modal compatibility equation are obtained. Based on them, a new theory of anisotropic viscoelastic mechanics is presented. The advantages of the theory are as follows: 1) the equations are all scalar, and independent of each other. The number of equations is equal to that of anisotropic subspaces, 2) no matter how complicated the anisotropy of solids may be, the form of the definite equation and the boundary condition are in common and explicit, 3) there is no distinction between the force method and the displacement method for statics, that is, the equilibrium equation and the compatibility equation are indistinguishable under the mechanical space, 4) each model equation has a definite physical meaning, for example, the modal equations of order one and order two express the volume change and shear deformation respectively for isotropic solids, 5) there also exist the potential functions which are similar to the stress functions of elastic mechanics for viscoelastic mechanics, but they are not man-made, 6) the final solution of stress or strain is given in the form of modal superimposition, which is suitable to the proximate calculation in engineering.     

Geometrically nonlinear,steady state vibration of viscoelastic beams

The problem of geometrically non-linear steady state vibrations of beams excited by harmonic forces is considered in this paper. The beams are made of a viscoelastic material defined by the classic Zener rheological model - the simplest model that takes into account all the basic properties of real viscoelastic materials. The constitutive stress-strain relationship for this type of material is given as a differential equation containing derivatives of both stress and strain. This significantly complicates the solution to the problem. The von Karman theory is applied to describe the effects of geometric nonlinearities of beam deformations. The equations of motions are derived using the finite element methodology. A polynomial approximation of bending moments is used. The order of basis functions is set so as to obtain a coherent approximation of moments and displacements. In the steady-state solution of equations of motion, only one harmonic is taken into account. The matrix equations of amplitudes are derived using the harmonic balance method and the continuation method is applied for solving them. The tangent matrix of equations of amplitudes is determined in an explicit form. The stability of steady-state solution is also examined. The resonance curves for beams supported in a different way are shown and the results of calculation are briefly discussed.     

INVESTIGATION OF THE THICKNESS VARIABILITY AND MATERIAL HETEROGENEITY EFFECTS ON FREE VIBRATION OF THE VISCOELASTIC CIRCULAR PLATES

The present study deals with free vibration analysis of variable thickness viscoelastic circular plates made of heterogeneous materials and resting on two-parameter elastic foundations in addition to their edge conditions.It is assumed that the viscoelastic material properties vary in the transverse and radial directions simultaneously.The complex modulus approach is employed in conjunction with the elastic-viscoelastic correspondence principle to obtain the solution.The governing equations are solved by means of a power series solution.Finally,a sensitivity analysis including evaluation of effects of various edge conditions,thickness variations,coefficients of the elastic foundation,and material loss factor and heterogeneity on the natural frequencies and modal loss factors is accomplished.     

The residual formulation of the method of manufactured solutions for computationally efficient solution verification

The method of manufactured solutions (MMS) is a solution verification methodology for determining whether the implementation of a discretization method is achieving its theoretical order of accuracy. This methodology combines the advantages of grid refinement studies and comparison with exact solution, by modifying the governing equations solved within a code by adding a source term to drive the solution towards a predetermined analytic function. By solving the modified equations on a sequence of grids and comparing the differences between the converged solution and manufactured solution, the order of accuracy of the implementation can be investigated. However, in its current form, converged solutions on a sequence of grids are required which can be quite costly and difficult to obtain. In this paper, by comparing the MMS to the method for determining the theoretical order of accuracy of a discretization method, the residual formulation of the MMS is developed. This new formulation only requires that the residual of the discretized governing equations to be calculated and not the solution to the discretized equations, thus avoiding the computational cost and difficulties inherent in obtaining converged solutions. Furthermore, since only the residuals are interrogated, individual components of the flow solver can be tested, in isolation, allowing the MMS to be used more effectively in locating errors within the code. This new approach is demonstrated to yield the same order of accuracy as the original MMS using three different cases—one-dimensional porous media equation, one-dimensional St Venant equations and two-dimensional unstructured Navier–Stokes simulations.     

Identification of inertia properties from the results of output-only modal analysis

Accurate knowledge of the rigid body properties of a structure including the mass, the location of mass center and the moments of inertia is important in machine design, vibration analysis, optimization and modeling of mechanical systems. However, estimation of these properties through theoretical methods is difficult when the structure has a complicated shape. In practice, the inertia properties can be estimated using the conventional modal testing methods by extracting the rigid body modes when the structure is tested in free-free boundary condition. However, all the rigid body modes are not always detectable, due to this fact that the structure is not excited at all degrees of freedom. In order to obtain all of the rigid body modes, many activities have been conducted for selecting the type and location of excitation without much success. In operational modal analysis (OMA), the structure can be excited at any arbitrary point and in different directions. In this paper, a new approach is introduced for estimation of the inertia properties from OMA. The data from OMA are adequate to extract all the rigid body modes of structure. A modal method is used for estimating the inertia properties from the rigid body modes extracted from OMA. The suggested approach is applied to a numerical model of a two-dimensional steel beam as well as a numerical model of a 3D frame and the accuracy of results is evaluated. It is shown that OMA can provide enough data to extract the inertia properties. A real beam is also tested in order to evaluate the performance of the method in practice, needless of a complicated procedure as for conventional methods.     

Dynamics of an exponentially damped solid rod: Analytic solution and finite element formulations

This paper presents the analytic solution together with two finite element formulations to obtain the dynamic response of an exponentially damped solid rod. The analysis is carried out in time domain, the material being modeled by the standard linear viscoelastic solid model. On the one hand, the analytic solution is reached by means of modal superposition by two different ways: with and without internal variables. On the other hand, two different finite element formulations are proposed: the first one is conceived for direct integration methods, and the other one is developed to apply modal superposition. The interest of modal superposition method lies on investigating the influence of non-viscous modes.     

主动约束层阻尼梁的数值模型

为对主动约束层阻尼结构建立精确完善的数学模型，采用有限元建模，并考虑到压电材料的机电耦合效应和粘弹性材料的本构关系随温度、频率的变化而变化的特点，将有限元方法与粘弹性材料的GHA模型相结合，从而避免因粘弹性材料导致的非线性微分方程，能直接求解模态频率、模态阻尼及结构响应。为进一步设计控制器，先在物理空间进行动力缩聚，将系统降至适当的维数，然后在状态空间用鲁棒防阶的方法进一步降阶。这样既能大大降低系统维数，又能保证降阶后系统稳定、可控、可观。这对于重量轻、柔度大、低频密集的大型空间柔性结构尤其重要。     

The elastic–viscoelastic correspondence principle for non-homogeneous materials with time translation non-invariant properties

This paper revisits the elastic–viscoelastic correspondence principle for non-homogeneous materials. Several recent publications discussed this principle for functionally graded materials (FGMs) with time translation invariant viscoelastic properties. It was demonstrated that the correspondence principle is valid only for the FGMs with separable relaxation moduli (moduli in separable form in space and time). This paper reconsiders this issue. It shows that the correspondence principle is valid even for non-homogeneous materials with separable relaxation moduli even if the time-dependences of the relaxation moduli in shear and dilatation are not necessarily time translation invariant. The property of similarity of Volterra operators is used to obtain the corresponding elastic solution. The correspondence is established between the elastic solution and the operator-transformed viscoelastic solution. The transformation operators are combinations of the Laplace transform operator and additional integral operators.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Modal acoustic impedance force on a spherical source near a rigid interface

The modal acoustic radiation load on a spherical surface undergoing angularly periodic axisymmetric harmonic vibrations while immersed in an acoustic halfspace with a rigid (infinite impedance) planar boundary is analyzed in an exact fashion using the classical technique of separation of variables. The formulation utilizes the appropriate wave field expansions, the classical method of images and the appropriate translational addition theorem to simulate the relevant boundary conditions for the given configuration. The associated acoustic field quantities such as the modal impedance matrix and the modal acoustic radiation force acting on the spherical surface are determined. The analytical results are illustrated with a numerical example in which the spherical surface, excited in vibrational modes of various orders, is immersed near an impervious rigid wall. The presented solution could eventually be used to validate those obtained by numerical approximation techniques.     

Response of a flexible tube to transverse impact I. Theory

The large non-linear response of a flexible viscoelastic tube under central transverse impact was studied by means of theoretical models, computer simulations and experiments. The present paper is concerned with the analytical solution of the system represented either as an elastic or viscoelastic (standard linear solid) string or Timoshenko beam with the local pinching deformation at the impact point expressed by either a linear elastic spring or a viscoelastic standard linear solid. These four models were solved numerically by the method of characteristics. A few comparisons of the results for the various models are presented here; correlations with corresponding finite element solutions and experimental data are described in a companion paper.     

一般粘弹结构的模态分析

分析了一般粘弹结构特征值问题的特点 ,建立了一般粘弹结构的模态分析方法。与粘弹结构已有的模态分析方法相比 ,该方法通用于更一般的粘弹结构 ,在形式上不涉及粘弹本构关系项 ,并只涉及一种模态向量     

Time–transient response for ultrasonic guided waves propagating in damped cylinders

Herein, an enhanced spectral finite element (SFE) formulation to calculate the time–transient response in cylindrical waveguides is proposed. The original aspect over SFE-based formulations consists in the possibility to account for the effect of material absorption, i.e. guided waves attenuation, on the calculation of the time–transient response.First, the damped steady-state response is constructed by a weighted superposition of the waveguide modal properties obtained from the spectral decomposition of the governing wave equation. To this purpose an enhanced spectrally formulated finite element is developed, in which material damping is included allowing for complex stress–strain viscoelastic constitutive relations in force of the correspondence principle. Dispersive modal properties for the damped waveguide (phase velocity, energy velocity, attenuation and wavestructures) follow straightforwardly by simple formulae. Next, the frequency response of the problem is calculated by weighting the modal data and the spectrum of the applied time-dependent force via Cauchy residue theorem. Finally, the inverse Fourier transform of the frequency response leads to the time–transient response for propagative damped guided waves.The approach is not restricted to any anisotropy degree, holds for any linear viscoelastic constitutive relation that can be characterized and formulated in the frequency domain and it can be applied to SFE formulations for arbitrary cross-section waveguides. A study on guided waves propagating in a scheduled 4.in-40 ANSI steel pipe is presented, where the steel is considered first as perfectly elastic and then as an hysteretic viscoelastic medium, in order to show the effect of material absorption on the time–transient response.     

离散型设置粘弹阻尼结构的复模态综合法

本文利用粘弹材料的本构方程,对粘弹阻尼结构的模态分析表达式进行了理论推导,进而将Craig的复模态综合法推广到粘弹阻尼系统.本文还发展了一种处理刚体模态的方法,从而克服了Craig法的不足,解决了含粘性阻尼、粘弹阻尼的复杂结构的隔振设计问题.     

Complex point sources in probe-corrected cylindrical near-field scanning

A probe-corrected theory based on complex point sources is presented for computing the acoustic field of an arbitrary finite source from measurements of the near field on a cylindrical scanning surface. The complex point sources are used both as basis functions for the expansion of the field outside the scan cylinder and for the representation of the probe. The resulting probe-corrected formulas are considerably simpler than the standard probe-corrected formulas based on cylindrical waves. The new formulation makes simulations of near-field scanning systems much less computationally intensive than simulations based on standard theory. The complex point-source theory is validated through numerical examples involving a baffled circular piston transducer probe.     

Experimental evaluation of modal damping in automotive components with different constraint conditions

The selection of reliable damping values for mechanical components of complex shape is very difficult from a theoretical point of view, above all in the case of complex constraint conditions. This work presents an example of experimental procedure for estimating the modal damping in a family of mechanical components. An Experimental Modal Analysis (EMA) program has been carried out on different types of brackets used on diesel engines, and in different constraint conditions with the aim of collecting the damping values in a database. In particular, the EMA has been carried out on brackets in the freely-supported condition and in the clamped condition actually used in the engine, highlighting the variations in modal damping with the different constraint conditions. The results of the experimental measurements have been processed with different modal analysis algorithms in order to increase the robustness of the solution. The resulting damping plots represent very useful data in the field of numerical analysis.     

GENERAL SECOND ORDER FLUID FLOW IN A PIPE

ＧＥＮＥＲＡＬＳＥＣＯＮＤＯＲＤＥＲＦＬＵＩＤＦＬＯＷＩＮＡＰＩＰＥＨｅＧｕａｎｇｙｕ（何光渝）（ＤｅｐａｒｔｍｅｎｔｏｆＰｅｔｒｏｌｅｕｍＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ＇ａｎＰｅｔｒｏｌｅｕｍＩｎｓｔｉｔｕｔｅ，Ｘｉ＇ａｎ７１００６１，Ｐ．Ｒ．Ｃｈｉ...