Nonlinear thermo-elastic buckling characteristics of cross-ply laminated joined conical–cylindrical shells

Here, the nonlinear thermo-elastic buckling/post-buckling characteristics of laminated circular conical–cylindrical/conical–cylindrical–conical joined shells subjected to uniform temperature rise are studied employing semi-analytical finite element approach. The nonlinear governing equations, considering geometric nonlinearity based on von Karman’s assumption for moderately large deformation, are solved using Newton–Raphson iteration procedure coupled with displacement control method to trace the pre-buckling/post-buckling equilibrium path. The presence of asymmetric perturbation in the form of small magnitude load spatially proportional to the linear buckling mode shape is assumed to initiate the bifurcation of the shell deformation. The study is carried out to highlight the influences of semi-cone angle, material properties and number of circumferential waves on the nonlinear thermo-elastic response of the different joined shell systems.     

Free vibration analysis of an elliptical plate with eccentric elliptical cut-outs

The elaborated collocation multipole method is employed to obtain a semi-analytical solution, involving proper products of angular and radial Mathieu functions, for the free flexural vibrations of a fully clamped thin elastic plate of elliptical planform containing multiple elliptical cutouts of arbitrary size, location, and orientation. The problem boundary conditions are satisfied by uniformly collocating points on the boundaries, and exactly calculating the normal derivative of plate displacement at the collocation points through use of appropriate directional derivative in each coordinate system. The multipole expansion is truncated to yield a coupled algebraic linear system of equations that is then solved for the nontrivial eigensolutions. Extensive numerical simulations present the first three calculated natural frequencies and the associated deformed mode shapes of an elliptical plate with elliptical/circular cutouts, for a wide range of plate/cutout aspect ratios, and cutout location/orientation parameters. The accuracy of solutions is checked through appropriate convergence studies, and the validity of results is established with the aid of a commercial finite element package as well as by comparison with the data available in the existing literature.     

Finite element methods for one-dimensional combustion problems

Three adaptive finite element methods based on equidistribution, elliptic grid generation and hybrid techniques are used to study a system of reaction–diffusion equations. It is shown that these techniques must employ sub-equidistributing meshes in order to avoid ill-conditioned matrices and ensure the convergence of the Newton method. It is also shown that elliptic grid generation methods require much longer computer times than hybrid and static rezoning procedures. The paper also includes characteristic, Petrov–Galerkin and flux-corrected transport algorithms which are used to study a linear convection–reaction–diffusion equation that has an analytical solution. The flux-corrected transport technique yields monotonic solutions in good agreement with the analytical solution, whereas the Petrov–Galerkin method with quadratic upstream-weighted functions results in very diffused temperature profiles. The characteristic finite element method which uses a Lagrangian–Eulerian formulation overpredicts the flame front location and exhibits overshoots and undershoots near the temperature discontinuity. These overshoots and undershoots are due to the interpolation of the results of the Lagrangian operator onto the fixed Eulerian grid used to solve the reaction–diffusion operator, and indicate that characteristic finite element methods are not able to eliminate numerical diffusion entirely.     

复合材料正交加筋板在横向载荷下挠度计算的半解析法

本文采用一种新的半解析法,即独特利用Heaviside函数建立与加筋板等效的变刚度模型来开展复合材料双向正交加筋板在横向载荷下的弯曲挠度分析．此模型可以准确地描述筋条在板面上的分布,以及由于筋条的存在而导致的板面刚度不均匀分布．使用Galerkin加权残值法求解该模型的控制方程,得到不同边界条件和载荷情况下的级数解．对于双向正交加筋板,将此半解析法的结果与传统均匀化方法和使用商业有限元软件ABAQUS建立的有限元模型所得到的弯曲挠度结果比较,验证了此方法的准确性和优越性．不同于传统均匀化方法,本双向正交加筋板的弯曲挠度半解析法可精确、有效地获取加筋间的局部弯曲挠度,可以促进复合材料结构的设计分析与优化的研究进展．     

二阶非自伴两点边值问题Garlerkin有限元的超收敛算法

提出了基于改进位移模式的二阶非自伴两点边值问题Garlerkin有限元的超收敛算法. 用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,基于Garlerkin 方法,采用积分形式推导了单元平衡方程. 对于线性单元,本文给出了有代表性的算例,结点和单元的位移、导数都达到了h⁴阶的超收敛精度.     

无单元Galerkin法和边界元耦合法

无网格方法与有限元法或边界元法耦合是无网格方法处理边界条件的方法之一,在无网格方法中研究无网格方法与有限元法或边界元法耦合的研究显得非常重要.本文在无单元Galerkin法和边界元法的基础上,基于无单元Galerkin法子域和边界元法子域的界面上位移连续和面力平衡条件,提出了一种新的无单元Galerkin法和边界元法的直接耦合方法,对弹性力学问题详细推导了在整个求解域上的耦合公式.与以往的耦合法相比,这种方法简单直观,不需要增加新的耦合区域,也不需要建立新的逼近函数来保证界面位移的连续性.算例结果表明,该方法具有较好的计算精度.     

Three-dimensional magneto-thermo-elastic analysis of functionally graded cylindrical shell

The present paper presents the three-dimensional magneto-thermo-elastic analysis of the functionally graded cylindrical shell immersed in applied thermal and magnetic fields under non-uniform internal pressure. The inhomogeneity of the shell is assumed to vary along the radial direction according to a power law function, whereas Poisson's ratio is supposed to be constant through the thickness. The existing equations in terms of the displacement components, temperature, and magnetic parameters are derived, and then the effective differential quadrature method(DQM) is used to acquire the analytical solution. Based on the DQM, the governing heat differential equations and edge boundary conditions are transformed into algebraic equations, and discretized in the series form. The effects of the gradient index and rapid temperature on the displacement,stress components, temperature, and induced magnetic field are graphically illustrated.The fast convergence of the method is demonstrated and compared with the results obtained by the finite element method(FEM).     

基于改进位移模式的一维有限元超收敛算法

将多尺度方法的思想与超收敛计算的解析公式结合起来,提出了改进有限元位移模式的算法。利用超收敛计算的解析公式,将高阶有限元解的位移模式用常规有限元解的位移模式表示。用常规有限元解的位移模式与高阶有限元解的位移模式之和构造新的位移模式,采用积分形式推导了单元刚度矩阵。该算法在前处理和后处理两个阶段都使用超收敛计算公式,在常规试函数的基础上,增加了高阶试函数,使得单元内平衡方程的残差减少,从而达到提高精度的目标。对于线性单元,本文结点和单元的位移、导数都达到了h4阶的超收敛精度。     

PERTURBATION FINITE ELEMENT METHOD FOR SOLVING GEOMETRICALLY NONLINEAR PROBLEMS OF AXISYMMETRICAL SHELL

In analysing the geometrically nonlinear problem of an axisymmetrical thin-walled shell, the paper combines the perturbation method with the finite element method by introducing the former into the variational equation to obtain a series of linear equations of different orders and then solving the equations with the latter. It is well-known that the finite element method can be used to deal with difficult problems as in the case of structures with complicated shapes or boundary conditions, and the perturbation method can change the nonlinear problems into linear ones. Evidently the combination of the two methods will give an efficient solution to many difficult nonlinear problems and clear away some obstacles resulted from using any of the two methods solely. The paper derives all the formulas concerning an axisym-metric shell of large deformation by means of the perturbation finite element method and gives two numerical examples,the results of which show good convergence characteristics.     

结构优化半解析灵敏度分析的改进算法

提出了结构半解析灵敏度分析的改进算法,该算法实现简便,对于设计变量摄动步长具有极佳的数值稳定特性。首先,从总体角度推导静力问题半解析法灵敏度分析新算法,提出了位移与应力灵敏度列式,并给出了算法实施途径;然后,将此思路推广于自振频率、屈曲临界荷载和瞬态响应等多种问题,提出了相应的计算步骤。以梁单元与壳单元等典型结构为例,开展了多个算例测试。测试表明,改进算法计算精度和效率均有提升,特别是设计变量步长有更大的数值稳定区域,为复杂工程结构形状优化的灵敏度分析提供了新途径。     

A finite wavelet domain method for wave propagation analysis in thick laminated composite and sandwich plates

An efficient numerical method is developed for the simulation of three dimensional transient dynamic response in thick laminated composite and sandwich plate structures involving very high frequencies and wave numbers. The proposed method incorporates Daubechies wavelet scaling functions for the interpolation of the in-plane displacements with a Galerkin formulation. It further explores the orthonormality and compact support of wavelet scaling functions to produce near diagonal consistent mass matrices and banded stiffness matrices. Hence, an uncoupled equivalent discrete spatial dynamic system is formulated, synthesized and rapidly solved in the wavelet domain using an explicit time integration scheme. The in-plane wavelet interpolation is further combined with an efficient high order layerwise laminate plate theory, that implements Hermite cubic splines for the through-the-thickness approximation of displacement fields. Numerical results are presented on the prediction of guided waves in laminated and thick sandwich composite plates and compared with respective solutions obtained by analytical, semi-analytical and time domain spectral element models. The method yielded higher convergence rates and substantial reductions in computational effort compared to respective time domain spectral finite elements.     

饱和土一维动力响应分析的弱式微分求积元法

基于伽辽金加权残值法，本文首先建立一维饱和土动力学控制微分方程的弱形式，而后分别采用微分求积法和有限元法将其空间坐标离散，得到以土体骨架位移、流体-土骨架相对位移和孔隙流体压力为自由度的单元离散方程，从而采用Crank-Nicolson 法求解．数值算例一方面通过与解析解的对比，验证了离散方程和数值程序的正确性．另一方面，通过地表位移和基底孔隙压力的收敛性分析，检验了求积元和有限元法的收敛效率．数值结果表明：所建立的弱式微分求积法在饱和土动力分析中不仅具有显著优于常规有限元法的收敛效率，而且还具有可变阶的收敛性能，为今后高效率分析提供了一种可能．     

有限元—康托洛维奇杂交法及其应用

本文将有限元法与康托洛维奇方法进行适当组合,吸收二者的主要优点,提出在单元的位移插值函数中附加康托洛维奇项。康托洛维奇方法是一种半解析法,其在解析方向具有较高的精度,在一定程度上弥补了原有限元方法中插值函数选取的盲目性,能够较好地反映微分方程的固有性质,提高其适应性。在单元的位移插值函数中附加内部无节点的位移项,无需增加新的单元与节点,使用较少的单元即可获得较高的精度。并且,这些附加项满足单元边界条件为零,故其在单元与单元的交界面上是保证协调的。本文通过算例充分说明了此方法的特点和优越性。     

变角度纤维复合材料圆柱壳的轴压屈曲

变角度纤维复合材料的纤维方向角可沿铺层面内连续变化,因此相应结构的性能具有更高的设计灵活性和更大的优化空间．本文假设纤维方向角沿圆柱壳的轴向呈正弦函数变化,对变角度纤维复合材料圆柱壳在两端简支边界条件下的轴压屈曲问题进行研究．基于Donnell经典壳体理论,推导变角度纤维复合材料圆柱壳的前屈曲控制方程并运用伽辽金法进行求解,然后采用瑞利里兹法求解屈曲问题．通过和现有文献及有限元数值结果的对比,验证了本文模型的收敛性和正确性,通过数值算例分析了纤维起始角和终止角的变化对圆柱壳的屈曲临界荷载的影响．本文的研究结果可为变角度纤维复合材料圆柱壳的分析和设计提供一定的参考．     

Shaft design: A semi-analytical finite element approach

This work deals with the practical use of semi-analytical finite elements in the machine design. The case of mechanical shafts is considered. The most usual loading condition characterized by the presence of axial, torsional, bending, and shear loads can be modeled by over imposing an axi-symmetric, an axi-antisymmetric and a harmonic load, corresponding to the first three terms of the Fourier series expansion, if semi-analytical plane finite element is used. A practical case is presented and the advantages, with respect to the three-dimensional approach in terms of computational time and accuracy for stress and displacement evaluation, is put in evidence.     

基于等几何分析的比例边界有限元方法

提出了一种具有比例边界有限元的半解析特性和等几何分析的几何特性的新方法。该新方法是在比例边界有限元框架中用NURBS曲线或曲面精确描述域边界几何形状,同时域边界位移场采用描述几何形状的NURBS形函数等参构造。这种新方法具有比例边界有限元固有的径向解析特性和NURBS的高阶连续性的优点。数值算例显示,与传统的比例边界有限元相比,基于等几何分析的比例边界有限元方法提高了域边界单元和域内应力场的连续性,减少了计算自由度。应用此方法可以用较少的计算自由度获得更高连续阶和更高精度的位移、应力和应变场。     

H 1 space-time discontinuous finite element method for convection-diffusion equations

An H~1 space-time discontinuous Galerkin (STDG) scheme for convectiondiffusion equations in one spatial dimension is constructed and analyzed. This method is formulated by combining the H~1 Galerkin method and the space-time discontinuous finite element method that is discontinuous in time and continuous in space. The existence and the uniqueness of the approximate solution are proved. The convergence of the scheme is analyzed by using the techniques in the finite difference and finite element methods. An optimal a-priori error estimate in the L~∞ (H~1 ) norm is derived. The numerical exper- iments are presented to verify the theoretical results.     

结构优化半解析灵敏度及误差修正改进算法

提出结构半解析灵敏度分析及其针对刚体位移的误差修正方法的改进算法, 构建灵敏度分析与误差修正项可分离形式. 该方法实现简便, 数值精度不受摄动步长与单元数目的影响. 首先从总体角度推得静力问题的误差修正半解析灵敏度分析方法, 提出了位移误差修正灵敏度列式, 并给出算法实施途径; 然后将此思路推广于自振频率、屈曲临界载荷问题, 提出了相应的计算步骤. 随后, 给出梁单元与壳单元误差修正项的具体推导方法, 并分别使用两种单元构建有限元模型进行算例测试. 结果表明, 该方法适用于多种分析类型, 数值精度不受单元数目与摄动步长的影响. 由于灵敏度分析与误差修正项可以分开计算, 该方法支持将误差修正项直接叠加于灵敏度求解结果进行误差修正, 使已有灵敏度分析程序得到充分利用. 尤其对于复杂工程结构的优化设计, 特别是形状优化设计以及尺寸、形状混合优化设计, 相比于原误差修正方法, 实现更为简便, 效率有所提升, 能为半解析灵敏度分析方法及其程序实现提供新的思路.