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

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

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

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

辐射换热下瞬态热一结构分析的一种空间薄壁杆单元

发展了一种用于辐射换热条件下瞬态热一结构分析的空间薄壁杆单元，其截面形式可以是任意形状的闭口截面和单支开口截面。该单元温度场分解为平均温度和多谐摄动温度，沿杆轴方向采用两结点线性插值，沿杆截面周向用三角函数展开，每结点含多个解耦的自由度，其中结点平均温度方程同传统一维温度有限元方程为非线性，各谐摄动温度方程为线性，然后利用Wilson-θ法求解结构的瞬态温度场。本文选择了两节点Bernoulli直梁单元得到准静态热弹性有限元方程并求解，针对非对称开口截面考虑了翘曲变形及弯扭耦合的影响。温度场引起的等效热载荷不仅包括常规的热轴力，还包括热弯矩以及热双力矩。本文针对不同截面形式的粱单元给出了瞬态温度场以及热变形的验证算例，并通过与商业程序中二维壳元计算结果的比较说明了本文所提出方法的正确性和高效性。     

半无限大体广义电磁热弹耦合的二维问题

基于Lord和Shulman广义热弹性理论,研究了热、电可导的半无限大体电磁热弹耦合的二维问题。半无限大体受热和外加恒定磁场的作用,文中建立了电磁热弹性耦合的控制方程,利用正则模态法求解得到了所考虑物理量的解析解,并用图形反映了各物理量的分布规律,从分布图上可以看出,介质中出现了电磁热弹耦合效应,各物理量的非零值仅在一个有限的区域内。     

层状弹性半空间非轴对称动力问题的奇异解

在柱坐标系下，利用关于方位角的Ｆｏｕｒｉｅｒ变换及关于径向的Ｈａｎｋｅｌ变换，将弹性力学基本方程组转化为非齐次的一阶常微分方程组的标准形式．采用求解微分方程组的矩阵法，建立了介质层的传递矩阵．由层间完全接触条件，导出了在任意埋藏源作用下层状弹性半空间频域奇异解，时域奇异解可通过关于频率的Ｆｏｕｒｉｅｒ积分得到．该方法可应用到固体、流体层的情况     

横向磁场中轴向运动铁磁梁的磁弹性主共振

研究在外激励力与磁场作用下轴向运动铁磁梁的磁弹性非线性主共振问题．基于弹性理论和电磁理论，给出梁的动能和弹性势能表达式，根据哈密顿原理，推导出磁场中轴向运动铁磁梁的磁弹性双向耦合非线性振动方程．通过伽辽金积分法进行离散，得出两端简支边界条件下铁磁梁磁弹性非线性强迫振动方程．应用多尺度法对方程进行求解，得出幅频响应方程．最后通过算例，给出铁磁梁的幅频特性曲线、振幅-磁感应强度和振幅-外激励力曲线并进行分析．结果显示，在幅频响应曲线中铁磁梁的轴向运动速度、外激励力、轴向拉力越大，共振振幅越大；而磁感应强度越大，振幅越小．     

半空间弹性波双参数强非线性逆散射

半空间弹性波强散射源逆问题是一个实用而又难度很大的强非线性问题，而经典的Born迭代方法只适用于弱非线性问题．基于DBI方法（Distorted-Born Iterative方法），提出了一种求解弹性波强非线性逆散射问题的迭代方法．在数值模拟运算时利用矩量法进行离散处理，并采用正则化原则避免求解病态矩阵方程。计算结果表明该方法具有较快的收敛速度及较高的精度。     

辐射换热下瞬态热-结构分析的一种空间薄壁杆单元

发展了一种用于辐射换热条件下瞬态热-结构分析的空间薄壁杆单元,其截面形式可以是任意形状的闭口截面和单支开口截面。该单元温度场分解为平均温度和多谐摄动温度,沿杆轴方向采用两结点线性插值,沿杆截面周向用三角函数展开,每结点含多个解耦的自由度,其中结点平均温度方程同传统一维温度有限元方程为非线性,各谐摄动温度方程为线性,然后利用Wilson-θ法求解结构的瞬态温度场。本文选择了两节点Ber-noulli直梁单元得到准静态热弹性有限元方程并求解,针对非对称开口截面考虑了翘曲变形及弯扭耦合的影响。温度场引起的等效热载荷不仅包括常规的热轴力,还包括热弯矩以及热双力矩。本文针对不同截面形式的梁单元给出了瞬态温度场以及热变形的验证算例,并通过与商业程序中二维壳元计算结果的比较说明了本文所提出方法的正确性和高效性。     

多层横观各向同性圆柱壳的三维弹性理论解

本文从三维弹性理论出发,用特征函数法研究多层横观各向同性圆柱壳的轴对称问题.把位移和应力分量的齐次解表达成特征函数展开式,并把特解部分用Fourier级数表示.以多层圆柱壳的内、外柱面作为齐次边界,同时考虑层间的连续条件,推导出问题的特征方程并用Muller法求解.文中运用传递矩阵技术处理多层问题,并用边界型最小二乘配点法处理端部边界条件.作为实例,对双层圆柱壳作了数值计算.     

多层弹性半空间问题解的精确刚度矩阵法

利用二维Fourier积分变换和刚度矩阵法,得到了直角坐标系下多层弹性半空间问题的精确解.在推导的过程中,直接从空间弹性力学的基本方程出发,利用积分变换等数学手段,首先推导出了单层空间弹性问题的刚度矩阵,然后按有限元法组成总体刚度矩阵.通过求解由总体刚度矩阵所构成的代数方程和积分逆变换,就可解出在任意静荷载作用下多层弹性半空间问题的精确解.由于刚度矩阵中不含有正指数项,计算时不会出现溢出现象,从而克服了传递矩阵法的缺点.由于在推导过程中摒弃了应力函数的选择,使得问题的求解更加合理化.最后还给出了计算实例来证明推导结果的准确性.     

A SEMI－ANALYTICAL AND SEMI－NUMERICAL METHOD FOR SOLVING 2－D SOUND－STRUCTURE INTERACTION PROBLEMS

Based on the transfer matrix method and the virtual source simulation technique,this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound-structure interaction problems under a harmonic excitation. Within any integration segment,as long as its length is small enough, along the circumferential curvilinear coordinate, the nonhomogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper. For the exterior fluid domain, the multi-circular virtual source simulation technique is adopted. The source density distributed on each virtual circular curve may be ex-panded as the Fourier‘‘s series. Combining with the inverse fast Fourier transformation, a higher accuracy and efficiency method for solving 2-D exterior Helmholtz‘s problems is presented in this paper. In the aspect of solution to the coupling equations, the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper. According to the superposition principle and compatibility conditions at the interface between the elastic ring and fluid, the algebraic equation of system can be directly constructed by using the least square approximation. Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented. Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.     

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

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

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.     

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

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

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.     

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

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

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.     

TRANSFER MATRIX METHOD FOR ANALYZING VIBRATION AND DAMPING CHARACTERISTICS OF ROTATIONAL SHELL WITH PASSIVE CONSTRAINED LAYER DAMPING TREATMENT

The first order differential matrix equations of the host shell and constrained layer for a sandwich rotational shell are derived based on the thin shell theory.Employing the layer wise principle and first order shear deformation theory, only considering the shearing deformation of the viscoelastic layer, the integrated first order differential matrix equation of a passive constrained layer damping rotational shell is established by combining with the normal equilibrium equation of the viscoelastic layer.A highly precise transfer matrix method is developed by extended homogeneous capacity precision integration technology.The numerical results show that present method is accurate and effective.