THE CONVERGENT CONDITION AND UNITED FORMULA OF STEP REDUCTION METHOD

The step reduction method was first suggested by Prof. Yeh Kai-yuan. This method has more advantages than other numerical methods. By this method, the analytic expression of solution can he obtained for solving nonuniform elastic mechanics. At the same time. its ealculuting time is very short and convergent speed very fast. In this paper, the convergent condition and nited formula of step reduction method are given by mathemutical method. It is proved that the solution of displacement and stress resultants obtained by this method can eonverge to exact solution uniformly, when the convergent condition is sutisfied. By united formula, the analytic solution solution can be expressed as matrix form, and therefore the former complicated expression can be avoMed. Two numerical examples are given at the end of this paper which indicate that. by the theory in this paper, a right model can be obtained for step reduction method.     

New matrix method for analyzing vibration and damping effect of sandwich circular cylindrical shell with viscoelastic core

Based on the linear theories of thin cylindrical shells and viscoelastic materials, a governing equation describing vibration of a sandwich circular cylindrical shell with a viscoelastic core under harmonic excitation is derived. The equation can be written as a matrix differential equation of the first order, and is obtained by considering the energy dissipation due to the shear deformation of the viscoelastic core layer and the interaction between all layers. A new matrix method for solving the governing equation is then presented With an extended homogeneous capacity precision integration approach. Having obtained these, vibration characteristics and damping effect of the sandwich cylindrical shell can be studied. The method differs from a recently published work as the state vector in the governing equation is composed of displacements and internal forces of the sandwich shell rather than displacements and their derivatives. So the present method can be applied to solve dynamic problems of the kind of sandwich shells with various boundary conditions and partially constrained layer damping. Numerical examples show that the proposed approach is effective and reliable compared with the existing methods.     

General analytic solution of dynamic response of beams with nonhomogeneity and variable cross-section

In this paper, a new method, the step-reduction method, is proposed to investigate the dynamic response of the Bernoulli-Euler beams with arbitrary nonhomogeneity and arbitrary variable cross-section under arbitrary loads. Both free vibration and forced vibration of such beams are studied. The new method requires to discretize the space domain into a number of elements. Each element can be treated as a homogeneous one with uniform thickness. Therefore, the general analytical solution of homogeneous beams with uniform cross-section can be used in each element. Then, the general analytic solution of the whole beam in terms of initial parameters can be obtained by satisfying the physical and geometric continuity conditions at the adjacent elements. In the case of free vibration, the frequency equation in analytic form can be obtained, and in the case of forced vibration, a final solution in analytical form can also be obtained which is involved in solving a set of simultaneous algebraic equations with only     

NONSTATIONARY PANDOM VIBRATION ANALYSIS OF LINEAR ELASTIC STRUCTURES WITH FINITE ELEMENT METHOD

At present,the finite element method is an efficient method for an-alyzing structural dynamic problems.When the physical quantitiessuch as displacements and stresses are resolved in the spectra andthe dynamic matrices are obtained in spectral resolving form,the re-lative equations cannot be solved by the vibration mode resolvingmethod as usual.For solving such problems,a general method is putforward in this paper.The excitations considered with respect tononstationary processes are as follows:P(t) ={P_i(t)}, P_i(t)=a_i(t)P_i~o(t),a_i(t) is a time function already known.We make Fourier transforma-tion for the discretized equations obtained by finite element method,and by utilizing the behaviour of orthcgonal increment of spectralquantities in random process,some formulas of relations about thespectra of excitation and response are derived.The cross power spec-tral denisty matrices of responses can be found by these formulas,then the structrual safety analysis can be made.When a_(i)(t)=1(i=1,2,…,n),the.m     

BOUNDARY INTEGRAL FORMULAS FOR ELASTIC PLANE PROBLEM OF EXTERIOR CIRCULAR DOMAIN

After the stress function and the normal derivative on the boundary for the plane problem of exterior circular domain are expanded into Laurent series, comparing them with the Laurent series of the complex stress function and making use of some formulas in Fourier series and the convolutions, the boundary integral formula of the stress function is derived further. Then the stress function can be obtained directly by the integration of the stress function and its normal derivative on the boundary. Some examples are given. It shows that the boundary integral formula of the stress function is convenient to be used for solving the elastic plane problem of exterior circular domain.     

Analytical and numerical methods of symplectic system for Stokes flow in two-dimensional rectangular domain

In this paper, a new analytical method of symplectic system, Hamiltonian system, is introduced for solving the problem of the Stokes flow in a two-dimensional rectangular domain. In the system, the fundamental problem is reduced to an eigenvalue and eigensolution problem. The solution and boundary conditions can be expanded by eigensolutions using adjoint relationships of the symplectic ortho-normalization between the eigensolutions. A closed method of the symplectic eigensolution is presented based on completeness of the symplectic eigensolution space. The results show that fundamental flows can be described by zero eigenvalue eigensolutions, and local effects by nonzero eigenvalue eigensolutions. Numerical examples give various flows in a rectangular domain and show effectiveness of the method for solving a variety of problems. Meanwhile, the method can be used in solving other problems.     

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations(ODEs) are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.     

THE PERTURBATION FINITE ELEMENT METHOD FOR SOLVING THE PLANE PROBLEM IN CONSIDERATION OF LINEAR CREEP

In this paper, amethod (PFMC) for solving plane problem of linear creep is presented by using perturbation finite element. It can be used in plane problem in consideration of creep, such as reinforced concrete beam, prestressed concrete beam, reinforced concrete cylinder and reinforced concrete tunnel in elastic or visco-elastic medium, as well as underground building and so on. In the presented method, the assumption made in the general increment method that variables remain constant in a divtded time interval is not taken. The accuracy is improved and the length of time step becomes larger. The computer storage can be reduced and the calculating efficiency can be increased. Perturbation finite element formulae for four-node quadrilateral isoparametric element including reinforcement are established and five numerical examples are given. As contrasted with the analytical solution, the accuracy is satisfactory.     

A SEMI-ANALYTICAL METHOD FOR THE VIBRATION OF AND SOUND RADIATION FROM A TWO-DIMENTIONAL BEAM-STIFFENED PLATE

A semi-analytical method based on space harmonics to investigate the vibration of and sound radiation from an infinite,fluid-loaded plate is presented.The plate is reinforced with two sets of orthogonally and equally spaced beam stiffeners,which are assumed to be line forces.The response of the stiffened plate to a convected harmonic pressure in the wave-number space is obtained by adopting the Green’s function and Fourier transform methods.Using the boundary conditions and space harmonic method,we establish the relationship between the stiffener forces and the vibration displacement of the plate.In this paper,the stiffener forces are expressed in terms of harmonic amplitudes of the plate displacement,which are calculated by using a numerical reduction technique.Finally,the Fourier inverse transform is employed to find expressions of the vibration and sound radiation in physical space.Agreements with existing results prove the validity of this approach and more numerical results are presented to show that this method converges rapidly.     

WAVE SUPERPOSITION METHOD BASED ON VIRTUAL SOURCE BOUNDARY WITH COMPLEX RADIUS VECTOR FOR SOLVING ACOUSTIC RADIATION PROBLEMS

By virtue of the comparability between the wave superposition method and the dynamic analysis of structures, a general format for overcoming the non-uniqueness of solution induced by the wave superposition method at the eigenfrequencies of the corresponding interior problems is proposed. By adding appropriate damp to the virtual source system of the wave superposition method, the unique solutions for all wave numbers can be ensured. Based on this thought, a novel method-wave superposition method with complex radius vector is constructed. Not only is the computational time of this method approximately equal to that of the standard wave superposition method, but also the accuracy is much higher compared with other correlative methods. Finally, by taking the pulsating sphere and oscillating sphere as examples, the results of calculation show that the present method can effectively overcome the non-uniqueness problem.     

求解二维结构-声耦合问题的一种半数值半解析方法

基于传递矩阵法和虚拟源强模拟技术提出了一种求解在谐激励作用下二维结构-声相互作用问题的半数值半解析法．在足够小的积分步长内，文中对任意形状弹性环沿周向曲线坐标的非齐次状态微分方程组，建立了一种齐次扩容方法．对于外声场，采用多圆形虚拟源强配置方案。并在每一条圆形配置曲线上将源强密度函数用Fourier级数展开，同时结合快速Fourier变换法，提出了一种高精度、高效率求解任意形状二维孔穴Helmholtz外问题的快速算法．在耦合方程的求解方面，根据叠加原理，将外激励和虚拟源强的Fourier级数展开项作为广义力分别作用在弹性环上，借助齐次扩容方法和精细积分法求得弹性环的状态向量，再利用流固交接条件和最小二乘法直接建立了耦合系统的求解方程．文中给出了二个典型弹性环在集中谐激励力作用下声辐射算例，计算结果表明该文方法较通常采用的混合FE-BE法更为有效．     

埋入水中旋转薄壳谐耦振分析的传递矩阵法

基于一般线弹性薄壳和势流理论,导出了旋转壳状态向量的一阶常微分矩阵方程和水动压力表达式,再借助齐次扩容技术和精细积分法,应用推广的传递矩阵法对埋入水中旋转薄壳的流固耦振进行了数值求解,并研究了一些因素对精度的影响.算例表明传递矩阵法和有限元方法相比不仅精度良好,而且有较高的计算效率.     

求解三维裂纹前缘SIF分布时间历程的通用权函数法和有限变分法

将三维热权函数法扩展为适用于表面力、体积力和温度载荷的通用权函数法(UWF).推导出以变分型积分方程表达的UWF法基本方程,从变分的角度,将求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE)改造为可以适用于一般的变分型积分方程的一类新型数值方法--有限变分法(FVM).在FVM中可以引入无穷多种线性无关的局部变分模式,可以根据计算要求在求解域中插入任意多个计算节点,单一型裂纹问题FVM所得到的最终方程组的系数矩阵总是一个对称的窄带矩阵,而且对角元总是大数,具有良好的数值计算性能.FVM对于SIF沿裂纹前缘急剧变化的复杂情况具有较好的数值模拟能力和较高的计算精度,利用自身一致性,可以求得三维裂纹前缘SIF的高精度解.     

具有约束层阻尼贮液圆柱容器的耦振与减振特性分析

该文在对空圆柱层合壳(不充液)振动分析研究的基础上,借助线性势流理论,进一步考虑液固耦合效应,导出了敷有CLD贮液圆柱层合壳谐耦振的一阶常微分矩阵方程,该方程右边多了液动压力项.由于它不能预先给定,导致方程中出现未知项.为了克服这一困难,文中研究了液动压力的解式,并采用新型齐次扩容精细积分技术,提出了一种高效率和高精度的半解析半数值解法.进而还求解了CLD贮液圆柱容器在地面谐运动激励下的响应;通过大量数值计算分析了CLD的厚度、长度、敷设位置以及粘弹芯的复剪切模量对减振效果的影响.计算还表明了有液体的存在,这些影响和空容器时是不一样的.文中方法为CLD贮液容器的控制优化提供了有力手段.     

复数矢径虚拟边界谱方法在二维空穴声辐射和散射问题中的应用

基于复数矢径虚拟边界积分法,通过将虚拟积分曲线上的未知源强密度函数用Fourier级数展开,同时借助快速数值Fourier变换计算程序,提出了一种求解二维任意形状空穴声辐射和散射问题的复数矢径虚拟边界谱方法．该方法具有以下特点：(1)不存在奇异积分处理;(2)采用复数矢径虚拟边界积分方法,不仅保证了解的唯一性,而且由于虚拟源强密度函数采用Fourier级数展开,克服了用单元离散方法不能用于较高频率范围的缺点;(3)采用快速数值Fourier变换技术使计算效率大幅度提高．文中给出的计算结果表明：在求解任意形状二维空穴声辐射和散射问题上较通常采用的FEM、BEM和VBEM更为有效．     

New matrix method for response analysis of circumferentially stiffened non-circular cylindrical shells under harmonic pressure

Based on the governing equation of vibration of a kind of cylindrical shells written in a matrix differential equation of the first order,a new matrix method is pre- sented for steady-state vibration analysis of a noncircular cylindrical shell simply sup- ported at two ends and circumferentially stiffened by rings under harmonic pressure.Its difference from the existing works by Yamada and Irie is that the matrix differential equation is solved by using the extended homogeneous capacity precision integration ap- proach other than the Runge-Kutta-Gill integration method.The transfer matrix can easily be determined by a high precision integration scheme.In addition,besides the normal interacting forces,which were commonly adopted by researchers earlier,the tan- gential interacting forces between the cylindrical shell and the rings are considered at the same time by means of the Dirac-δfunction.The effects of the exciting frequencies on displacements and stresses responses have been investigated.Numerical results show that the proposed method is more efficient than the aforementioned method.     

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation effectiveness and precision, is presented for solving the acoustic radiation from a submerged infinite non-circular cylindrical shell stiffened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work,besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The effects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efficient than the existing mixed FE-BE method.     

求解水下纵向加肋无限长非圆柱壳声辐射问题的一种新的半解析方法

基于齐次扩容精细积分法和复数矢径虚拟边界谱方法，利用Fourier积分变换和稳相法，提出了一种具有较高效率和精度的新的求解水下纵向加肋无限长非圆柱壳声辐射问题的半解析方法．考虑了非圆柱壳和肋骨之间同时存在多种相互作用力和力偶矩，较已往很多学者仅计及法向相互作用力更加符合实际．不仅比较了该文方法和精确解计算纵向加肋圆柱壳在集中点力激励下的声辐射计算结果，同时还研究了肋骨数量、大小以及椭圆柱壳横截面椭圆度对声辐射特性的影响．数值计算结果表明该文方法较已有的混合FE-BE法更为有效．     

求解变分型积分方程的一种新型数值方法——有限变分法

将作者提出的多虚拟裂纹扩展法(MVCE法)拓展为求解变分型积分方程问题的一种新型数值方法——有限变分法(FVM)。它的基本思想是,给定有限个(N个)局部变分模式,将所求解的未知量用适当的方法离散化,针对这N个局部变分模式列出N个方程,求解N个未知系数,从而求得未知量。单一未知变量FVM的最终方程组的系数矩阵通常是一个对称的窄带矩阵,对角元是大数,有很好的数值计算性能。用FVM求解了三维I型裂纹前缘的应力强度因子(SIF)分布。利用基于FVM的通用权函数法计算程序,可以高精度和高效率地求解表面力、体积力和温度载荷共同作用情况下三维裂纹前缘SIF的分布及其时间历程。FVM可以被推广到更广泛的领域,是一个求解变分型积分方程问题的普遍适用的新型数值方法。