二维复合材料结构模态分析的边界元法

将边界元方法用于分析二维复合材料结构的自由振动模态，利用特解处理体积力（惯性力）仅需静态基本解就可求解问题，对一各向同性悬臂梁，用该法得到的结果与用有限元或各向性边界元法得到的结果符合得很好，但该法可解各向异性问题，对层状复合材料简支梁，用该法得到了数值结果与用一维层状复合材料梁的理论解的比较表明，当结构的长厚比大于２０时，二者符合得很好，当结构的长厚比小于２０时，一维层状材料梁的理论将产生很大的     

Boundary element method for arbitrary elastic thin shell

In this paper, a boundary element scheme for arbitrary elastic thin shells is elaborated, Based on BEM of 3D linear elasticity and Kirchhoff's hypothesis, boundary integral equations for shells are deduced. As a result, only Kelvin's solution is used, the difficulty in finding a fundamental solution of arbitrary shells is successfully avoided. The project supported by National Natural Science Foundation of China     

BEM for modal analysis of 3-D anisotropic structures

The boundary element method for the modal anlysis of free vibration for 3-D anisotropic structures using particular solutions has been developed. The complete polynomials of order two are used to construct the particular solutions for general anisotropic materials. The numerical results for 3-D free vibration analysis of an isotropic cantilever beam by the method presented is in good agreement with the results using the Ritz technique. For anisotropic materials, the numerical results calculated from the proposed method are in good agreement with the results from MSC.NASTRAN. Project supported by the National Natural Science Foundation of China (No. 19872011).     

Stochastic boundary element method for reliability analysis of plate and beams composite structures

I.IntroductionTheengineeringstructuresareoftells,'1>:'rectcdtotheactionofthestochasticloadingthatvarieswiththetime,forexample,theengineeringstructuresactedonbytheearthquake,theoceanstructuresactedonbydynamicpressure,andthevehiclesofthetransportationinflue…     

复合材料周期结构数学均匀化方法的一种新型单胞边界条件

数学均匀化方法是计算周期复合材料结构的有效方法之一,单胞边界条件施加的合理性直接决定了影响函数控制方程的计算效率和精度,进而影响均匀化弹性参数和摄动位移的计算精度.本文首先将单胞影响函数作为虚拟位移处理,给出了单胞在结构中真实的边界条件,结果表明,四边固支适合作为二维结构单胞边界条件;其次,针对二维结构提出了超单胞周期边界条件,有效提高了影响函数的计算精度,并使用与虚拟位移相对应的虚拟势能泛函验证超单胞周期边界条件的有效性;最后,利用数值分析验证多尺度渐进展开方法的计算精度,强调了二阶摄动的必要性.     

基于局部移动Kriging无网格方法的层合板自由振动分析

利用基于局部移动Kriging插值无网格法对层合板自由振动进行了数值分析,基于一阶剪切层合理论导出了层合板振动的控制方程和边界条件,进一步得到了自由振动的离散化特征方程。由于Kriging插值函数具有Kronecker delta函数性质,可以直接施加本质边界条件。通过本文给出的方法,对不同边界条件、不同跨厚比、不同材料参数和铺设角度的层合板的振动频率进行了计算,均得到满意结果。最后用该方法对层合板的铺设角度进行优化设计,得到了与已有文献完全一致的优化结果。数值结果充分表明了无网格Kriging方法分析层合板自由振动问题的有效性和高精确度。     

复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法

基于一阶剪切变形理论,提出了复合材料层合板自由振动分析的无网格自然邻接点Petrov-Galerkin法。计算时在复合材料层合板中面上仅需要布置一系列的离散节点,并利用这些节点构建插值函数。在板中面上的局部多边形子域上,采用加权余量法建立复合材料层合板自由振动分析的离散化控制方程,并且这些子域可由Delaunay三角形方便创建。自然邻接点插值形函数具有Kronecker delta函数性质,因而无需经过特别处理就能准确地施加本质边界条件。对不同边界条件、不同跨厚比、不同材料参数和不同铺设角度的复合材料层合板,由本文提出的无网格自然邻接点Petrov-Galerkin法进行自由振动分析时均可得到满意的结果。数值算例结果表明,本文方法求解复合材料层合板的自由振动问题是行之有效的。     

Free vibration analyses by boundary element method using component mode synthesis

The boundary element method for dynamic analysis of two dimensional structures was developed by Nardini and Brebbia in 1983. Free vibration analysis of structures by BEM results in a generalized eigenproblem with coefficient matrices nonsymmetric and full. So, it costs a lot of CPU time to solve the eigenproblem. In order to reduce the cost the improved free interface coupling method has been adopted. By aid of lower eigenvectors and residual modes the degrees of freedom of substructure equaion of dynamic equilibrium can be reduced. By assembling the reduced substructure equations, the system of equations for the whole structure will be obtained. By using the traction and displacement compatibility conditions on the interface, the second transformation matrix has been developed for further reduction of system equations. The technique has been applied to some specific examples which give the evidence that the method is satisfactory and efficient.     

Global composite element iteration for analysis of seepage free surface

ＩｎｔｒｏｄｕｃｔｉｏｎＳｅｅｐａｇｅｆｒｅｅｓｕｒｆａｃｅｉｓｏｆｔｅｎｅｎｃｏｕｎｔｅｒｅｄｉｎｔｈｅｃｏｕｒｓｅｏｆｓｅｅｐａｇｅｐｒｏｂｌｅｍｔｈｒｏｕｇｈｍｅｄｉａｏｆｓｌｏｐｅ,ｅａｒｔｈｄａｍａｎｄｔｕｎｎｅｌｅｎｇｉｎｅｅｒｉｎｇ．Ｔｈｅｐｈｅｎｏｍｅｎａｏｆｓｅｅｐａｇｅｆｒｅｅｓｕｒｆａｃｅｉｓｏｎｅｏｆｔｈｅｉｍｐｏｒｔａｎｔａｎｄｄｉｆｆｉｃｕｌｔｔｏｐｉｃｓｏｆｒｏｃｋｍａｓｓｈｙｄｒｏｍｅｃｈａｎｉｃｓ［１～３］,ａｎｄｑｕｉｔｅａｎｉｍｐｏｒｔａｎｔｒｅｓｅａｒｃｈ…     

求解不连续中厚板自由振动的微分容积单元法

基于区域叠加原理和微分容积法,发展了一种新型的数值方法——微分容积单元法,用以分析具有不连续几何特征的中厚板的自由振动。根据板的不连续情况将其划分为若干单元,在每个单元内用微分容积法将控制微分方程离散成为一组线性代数方程．在相邻的单元连接处应用位移连续条件和平衡条件,引入边界约束条件后得到一套关于各配点位移的齐次线性代数方程,由此可导出求解系统固有频率的特征方程。本文用子空间迭代法求解特征方程,并以开孔板、混合边界条件板和突变厚度板为例研究了方法的收敛性和计算精度。     

Boundary element analysis of contact problems using artificial boundary node approach

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.     

压电介质边界元法及奇异性处理

本文从压电材料的基本方程出发,利用功的互等原理推导了边界积分方程,并详细地讨论了边界元的计算步骤,利用等参变换,着重研究了在边界元计算中基本解的奇异性问题,对各种情况讨论了系数矩阵H和G的算法,并给出院具体的表达式,作为算例,选取了均匀薄板和开孔薄板PZY-4压电材料,计算结果表明,本文提出的边界元的计算格式和奇异性的处理方法相当有效。     

The 3-D boundary element method of roller bearing by plate element analogue

The mill roller bearing is made up of an internal ring, middle rolls and an external ring, the analysis of which is a multi-bodies contact problem. In this paper, based on the three-dimensional elastic contact BEM without friction, and using the structural characteristics of roller bearings, middle rolls are described by elastic plate units of different shapes, which is placed on the internal ring. The discontinuous traction problem of the contact element can be dealt with by the traction sub-element method. Therefore the contact problem can be changed from one of multi-bodies to one of two-bodies. As an elastic deformation of the plate element, it can be calculated by an elastic contact formula and is regarded as an equivalent gap to be put into the total matrix equality. The face-to-face contact condition of displacement and traction on the element will serve as the convergence judging condition. So convergence solution can be obtained before any geometry inconsistent on the contact element may take place. Therefore, the three-dimensional distribution of load of the mill roller bearing can be achieved. It has been proved that the calculating model of this method is visualized and simple, has fast convergence solution and higher accuracy, making it an effective numerical method of designing and analyzing load characteristics of the mill roller bearing. Project supported by the National Natural Science Foundation of China (50075075).     

BOUNDARY ELEMENT METHOD FOR SOLVING DYNAMICAL RESPONSE OF VISCOELASTIC THIN PLATE(Ⅱ)──THEORETICAL ANALYSIS

In this paper, the necessary theoretical analysis for the approximation boundary element method to solve dynamical response of viscoelastic thin plate presented in [1] is discussed. The theorem of existence and uniqueness of the approximate solution and the error estimation are also obtained. Based on these conclusions, the principle for choosing the mesh size and the number of truncated terms in the fundamental solution are given. It is shown that the theoretical analysis in this paper are consistent with the numerical results in [1]. Project supported by the National Natural Science Fundation of China     

NUMERICAL SIMULATION OF 2D FIBER-REINFORCED COMPOSITES USING BOUNDARY ELEMENT METHOD

Introduction Withthedevelopmentofmodernindustry,compositesareincreasinglybeingappliedto agreatnumberofimportantstructures.Todeterminethemacroscopicaleffective characteristicsofcompositesisanessentialprobleminmanyengineeringapplications.The macroscopicalef…     

SOLVING THE FREE BOUNDARY PROBLEM IN CONTINUOUS CASTING BY USING BOUNDARY ELEMENT METHOD

ＳＯＬＶＩＮＧＴＨＥＦＲＥＥＢＯＵＮＤＡＲＹＰＲＯＢＬＥＭＩＮＣＯＮＴＩＮＵＯＵＳＣＡＳＴＩＮＧＢＹＵＳＩＮＧＢＯＵＮＤＡＲＹＥＬＥＭＥＮＴＭＥＴＨＯＤＬｉＹａｏｙｏｎｇ（李耀勇）；ＺｈａｎｇＺｈｉｌｉ（张自立）（ＲｅｃｅｉｖｅｄＪｕｎｅ，１８，１９...     

Boundary element method for moving and rolling contact of 2D elastic bodies with defects

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes. The project supported by the National Natural Science Foundation of China (19772025)     

Boundary element method to solve torsion problem of cracked cylinder

Separating the discontinuous solution by use of the single crack solution, together with the regular solution of harmonic function, the torsion problem of a cracked cylinder is reduced to solving a set of mixed-type integral equations and its numerical technique is then proposed by combining the numerical method of singular integral equation with the boundary element method. Several numerical examples are calculated which will be useful to engineering practice. The method proposed is characterized by its fine accuracy and convenience for using, which can be extended to the cases of multiple crack.The project supported by National Natural Science Foundation of China.     

螺栓松动边界下纤维增强复合薄板的振动测试方法

提出了基于预埋压力传感器的量化测试方法,研究了螺栓松动边界对纤维增强复合薄板振动特性的影响。首先,自主设计并开发了带有预埋压力传感器的螺栓松动边界下复合薄板的振动测试系统,并详细介绍了系统各个部件的组成和功能;然后,归纳出一套合理、规范的松动边界下复合薄板的振动测试流程,并对HF10碳纤维/树脂复合薄板进行了实际测试。结果表明:随着螺栓松动程度的不断增加,复合薄板的固有频率逐渐降低,模态振型的节线位置也发生了不同程度的变化,但其阻尼结果呈现先增大后减小的趋势;而共振和非共振响应呈现先减小后增大的趋势。     

Transfer matrix approach of vibration isolation analysis of periodic composite structure

The transmission properties of elastic waves propagating in a three-dimensional composite structure embedded periodically with spherical inclusions are analyzed by the transfer matrix method in this paper. Firstly, the periodic composite structures are divided into many layers, the transfer matrix of monolayer structure is deduced by the wave equations, and the transfer matrix of the entire structure is obtained in the case of boundary conditions of displacement and stress continuity between layers. Then, the effective impedance of the structure is analyzed to calculate its reflectivity and transmissivity of vibration isolation. Finally, numerical simulation is carried out; the experiment results validate the accuracy and feasibility of the method adopted in the paper and some useful conclusions are obtained. Project (No. 50075029) supported by the National Natural Science Foundation of China.