Thermal buckling analysis of functionally graded cylindrical shells

Thermal buckling behavior of cylindrical shell made of functionally graded material(FGM) is studied. The material constituents are composed of ceramic and metal.The material properties across the shell thickness are assumed to be graded according to a simple power law distribution in terms of the volume fraction rule of mixtures. Based on the Donnell shell theory, a system of dimensionless partial differential equations of buckling in terms of displacement components is derived. The method of separation of variables is used to transform the governing equations to ordinary differential equations(ODEs). A shooting method is used to search for the numerical solutions of the differential equations under two types of boundary conditions. Effects of the power law index, the dimensionless geometrical parameters, and the temperature ratio on the critical buckling temperature are discussed in detail.     

The displacement solution of conical shell and its application

In this paper,the displacement solution method of the conical shell is presented.Fromthe differential equations in displacement form of conical shell and by introducing adisplacement function,U(s,θ),the differential equations are changed into an eight-ordersoluble partial differential equation about the displacement function U(s,θ)in which thecoefficients are variable.At the same time,the expressions of the displacement and internalforce components of the shell are also given by the displacement function.As special casesof this paper,the displacement function introduced by V.Z.Vlasov in circular cylindricalshell,the basic equation of the cylindrical shell and that of the circular plate are directlyderived.Under the arbitrary loads and boundary conditions,the general bending problem of theconical shell is reduced to finding the displacement function U(s,θ),and the generalsolution of the governing equation is obtained in generalized hypergeometric function,Forthe axisymmetric bending deformation of the     

The exact solution for the general bending problems of conical shells on the elastic foundation

The general bending problem of conical shells on the elastic foundation (Winkler Medium) is not solved. In this paper, the displacement solution method for this problem is presented. From the governing differential equations in displacement form of conical shell and by introducing a displacement function U(s,θ), the differential equations are changed into an eight-order soluble partial differential equation about the displacement function U(s,θ) in which the coefficients are variable. At the same time, the expressions of the displacement and internal force components of the shell are also given by the displacement function U(s θ). As special cases of this paper, the displacement function introduced by V.S. Vlasov in circular cylindrical shell^[5], the basic equation of the cylindrical shell on the elastic foundation and that of the circular plates on the elastic foundation are directly derived.Under the arbitrary loads and boundary conditions, the general bending problem of the conical shell on the elastic foundation is reduced to find the displacement function U(s,θ).The general solution of the eight-order differential equation is obtained in series form. For the symmetric bending deformation of the conical shell on the elastic foundation, which has been widely usedinpractice,the detailed numerical results and boundary influence coefficients for edge loads have been obtained. These results have important meaning in analysis of conical shell combination construction on the elastic foundation,and provide a valuable judgement for the numerical solution accuracy of some of the same type of the existing problem.     

角度非均匀材料平面V形切口应力奇性分析

提出了一种确定角度非均匀材料平面V形切口尖端应力奇性指数的有效方法。首先,在弹性力学基本方程中引入V形切口尖端位移场的级数渐近展开,建立以位移为特征函数的变系数和非线性微分方程组。然后,采用微分求积法(DQM)求解微分方程组,可得到多阶应力奇性指数及其相对应的特征函数,该法具有公式简单、编程方便、计算量少和精度高等优点,可处理任意开口角度和任意材料组合的V形切口。典型算例验证了微分求积法的有效性和精确性。     

Three-dimensional free vibration analysis of rotating laminated conical shells: layerwise differential quadrature (LW-DQ) method

This paper focuses on the free vibration analysis of thick, rotating laminated composite conical shells with different boundary conditions based on the three-dimensional theory, using the layerwise differential quadrature method (LW-DQM). The equations of motion are derived applying the Hamilton’s principle. In order to accurately account for the thickness effects, the layerwise theory is used to discretize the equations of motion and the related boundary conditions through the thickness of the shells. Then, the equations of motion as well as the boundary condition equations are transformed into a set of algebraic equation applying the DQM in the meridional direction. This study demonstrates the applicability, accuracy, stability and the fast rate of convergence of the present method, for free vibration analyses of rotating thick laminated conical shells. The presented results are compared with those of other shell theories obtained using conventional methods and a special case where the angle of the conical shell approaches zero, that is, a cylindrical shell and excellent agreements are achieved.     

轴向均布载荷下压杆稳定问题的DQ解

叙述了微分求积法(differential quadrature method)的一般方法，研究用微分求积法求解在均布轴向载荷下细长杆的稳定问题．通过Newton-Raphson法求解非线性方程组，以及对问题进行线性假设后求解广义特征值方程，得到了精度很高的后屈曲挠度数值和临界载荷数值．与解析解和其他近似解相比，微分求积法具有较高的精度和简便性．     

Solution to transient Navier–Stokes equations by the coupling of differential quadrature time integration scheme with dual reciprocity boundary element method

The two‐dimensional time‐dependent Navier–Stokes equations in terms of the vorticity and the stream function are solved numerically by using the coupling of the dual reciprocity boundary element method (DRBEM) in space with the differential quadrature method (DQM) in time. In DRBEM application, the convective and the time derivative terms in the vorticity transport equation are considered as the nonhomogeneity in the equation and are approximated by radial basis functions. The solution to the Poisson equation, which links stream function and vorticity with an initial vorticity guess, produces velocity components in turn for the solution to vorticity transport equation. The DRBEM formulation of the vorticity transport equation results in an initial value problem represented by a system of first‐order ordinary differential equations in time. When the DQM discretizes this system in time direction, we obtain a system of linear algebraic equations, which gives the solution vector for vorticity at any required time level. The procedure outlined here is also applied to solve the problem of two‐dimensional natural convection in a cavity by utilizing an iteration among the stream function, the vorticity transport and the energy equations as well. The test problems include two‐dimensional flow in a cavity when a force is present, the lid‐driven cavity and the natural convection in a square cavity. The numerical results are visualized in terms of stream function, vorticity and temperature contours for several values of Reynolds (Re) and Rayleigh (Ra) numbers. Copyright © 2008 John Wiley & Sons, Ltd.     

基于边界值方法的微分动力系统数值计算方法

对高维非线性初值问题,微分求积法在每一步的积分过程中需要求解一个更高维的非线性方程组,因而计算量巨大。基于微分求积法与边界值方法两者之间的关系,可以将广义向后差分方法和扩展的隐式梯形积分方法看作是经典微分求积法的稀疏表达形式。将广义向后差分方法以及扩展的隐式梯形积分方法这两类边界值方法应用于微分动力系统的数值计算,提出了一类新的数值计算方法。理论分析及算例结果表明,对高维非线性微分初值问题的数值计算,本文方法相对于经典的微分求积法具有更高的计算效率。     

Free vibration and critical speed of moderately thick rotating laminated composite conical shell using generalized differential quadrature method

The generalized differential quadrature method (GDQM) is employed to consider the free vibration and critical speed of moderately thick rotating laminated composite conical shells with different boundary conditions developed from the first-order shear deformation theory (FSDT). The equations of motion are obtained applying Hamilton’s concept, which contain the influence of the centrifugal force, the Coriolis acceleration, and the preliminary hoop stress. In addition, the axial load is applied to the conical shell as a ratio of the global critical buckling load. The governing partial differential equations are given in the expressions of five components of displacement related to the points lying on the reference surface of the shell. Afterward, the governing differential equations are converted into a group of algebraic equations by using the GDQM. The outcomes are achieved considering the effects of stacking sequences, thickness of the shell, rotating velocities, half-vertex cone angle, and boundary conditions. Furthermore, the outcomes indicate that the rate of the convergence of frequencies is swift, and the numerical technique is superior stable. Three comparisons between the selected outcomes and those of other research are accomplished, and excellent agreement is achieved.     

板锥网壳结构的连续化分析

板锥网壳结构是一种受力性能合理，技术经济效益良好的新型空间结构形式。本文将板锥网壳结构连续化为能共同作用的特殊形式的三层薄壳，按薄壳理论建立其位移法和混合法的基本微分方程。通过对微分方程的求解，计算其整体位移及结构内力。该法具有一定的精度，可宏观地了解结构的力学性能，并可用于板锥网壳结构的初步设计。     

A mixed modal-differential quadrature method for free and forced vibration of beams in contact with fluid

A simple and accurate mixed modal-differential quadrature formulation is proposed to study the dynamic behavior of beams in contact with fluid. Both free and forced vibration problems are considered. The proposed mixed methodology uses the modal technique for the structural domain while it applies the differential quadrature method (DQM) to the fluid domain. Thus, the governing partial differential equations of the beam and fluid are reduced to a set of ordinary differential equations in time. In the case of forced vibration, the Newmark time integration scheme is employed to solve the resulting system of ordinary differential equations. The proposed formulation, in general, combines the simplicity of the modal method and high accuracy and efficiency of the DQM. Its application is shown by solving some beam-fluid interaction problems. Comparisons with analytical solutions show that the present method is very accurate and reliable. To demonstrate its efficiency, the test problems are also solved using the finite element method (FEM). It is found that the proposed method can produce better accuracy than the FEM using less computational time. The technique presented in this investigation is general and can be used to solve various fluid-structure interaction problems.     

Boundary element solution of unsteady magnetohydrodynamic duct flow with differential quadrature time integration scheme

A numerical scheme which is a combination of the dual reciprocity boundary element method (DRBEM) and the differential quadrature method (DQM), is proposed for the solution of unsteady magnetohydrodynamic (MHD) flow problem in a rectangular duct with insulating walls. The coupled MHD equations in velocity and induced magnetic field are transformed first into the decoupled time‐dependent convection–diffusion‐type equations. These equations are solved by using DRBEM which treats the time and the space derivatives as nonhomogeneity and then by using DQM for the resulting system of initial value problems. The resulting linear system of equations is overdetermined due to the imposition of both boundary and initial conditions. Employing the least square method to this system the solution is obtained directly at any time level without the need of step‐by‐step computation with respect to time. Computations have been carried out for moderate values of Hartmann number (M?50) at transient and the steady‐state levels. As M increases boundary layers are formed for both the velocity and the induced magnetic field and the velocity becomes uniform at the centre of the duct. Also, the higher the value of M is the smaller the value of time for reaching steady‐state solution. Copyright © 2005 John Wiley & Sons, Ltd.     

Effect of temperature dependent viscosity on revolving axi-symmetric ferrofluid flow with heat transfer

The prime objective of the present study is to examine the effect of temperature dependent viscosity μ(T) on the revolving axi-symmetric laminar boundary layer flow of an incompressible, electrically non-conducting ferrofluid in the presence of a stationary plate subjected to a magnetic field and maintained at a uniform temperature. To serve this purpose, the non-linear coupled partial differential equations are firstly converted into the ordinary differential equations using well-known similarity transformations. The popular finite difference method is employed to discretize the non-linear coupled differential equations. These discretized equations are then solved using the Newton method in MATLAB, for which an initial guess is made with the help of the Flex PDE Solver. Along with the velocity profiles, the effects of temperature dependent viscosity are also examined on the skin friction, the heat transfer, and the boundary layer displacement thickness. The obtained results are presented numerically as well as graphically.     

含孔隙变厚度FG圆板的湿热力学响应

功能梯度材料(functionally graded materials, FGM)是组份含量按特定方向连续变化的非均匀复合材料,可有效解决传统复合材料组份之间结合能力弱和不同组份性能难以协调等问题,达到诸如缓和应力集中和优化应力分布等效果,使整体材料在保持细观结构完整性的同时充分发挥各组份材料的性能优势.由于制备技术等原因或出于特殊功能的需要,微孔或孔隙是各类型FGM中的常见缺陷.从细观结构上看,多孔FGM中的孔隙包含了单一组份内的材料孔隙和组份微粒间的结构间隙,这些孔隙将对FGM的力学性能,尤其是在湿热环境下的力学行为产生影响.本文考虑FGM中的两类细观孔隙(材料孔隙和结构孔隙),提出了令各类孔隙依赖于各自组份变化,再线性叠加得到的整体孔隙计算式.考虑组份材料和孔隙填充物(液相水和水蒸气)性质的温度相关性,建立了湿热相关FGM材料模型.针对厚度沿径向变化的旋转圆板结构,应用该FGM材料模型,推导了圆板的非线性稳态湿热控制方程及考虑湿热弹性本构的位移控制方程,采用微分求积法(differential quadrature method,DQM),获得了圆板的湿热场、位移场和应力分布.在数值算例中,利用退化模型的解析解对本文的数值计算方法进行了验证,继而通过改变各关键参数,讨论了两类孔隙率、梯度指数和圆板厚度变化对含孔隙FGM变厚度旋转圆板湿热力学响应的影响规律.      

导电圆柱薄壳的热磁弹性效应分析

研究了磁场环境中受机械载荷作用的导电圆柱薄壳的热磁弹性问题．首先,根据电动力学方程和广义Ohm定律,得到了导电薄壳电流密度的分布,考虑到Joule热效应及热平衡方程,得到了导电薄壳的温度分布．其次,通过几何方程、物理方程、运动方程和电动力学方程导出了导电薄壳在机械场、电磁场以及温度场作用下的基本方程．最后,采用差分法及准线性化方法,得到了可以应用离散正交法求解的准线性微分方程组．对于导电圆柱薄壳,得到了Lorentz力表达式,并且推导了温度场积分特征值．讨论了导电圆柱薄壳应力、温度及变形随外加电磁参量的变化规律,并通过实例证实了可以通过改变电、磁、力场的参数来实现对薄壳的应力、应变、温度的控制．     

ANALYSES ON NONLINEAR COUPLING OF MAGNETO-THERMO-ELASTICITY OF FERROMAGNETIC THIN SHELL I： GENERALIZED VARIATIONAL THEORETICAL MODELING

Based on the generalized variational principle of magneto-thermo-elasticity of the ferromagnetic elastic medium, a nonlinear coupling theoretical modeling for a ferromagnetic thin shell is developed. All governing equations and boundary conditions for the ferromagnetic shell are obtained from the variational manipulations on the magnetic scalar potential, temperature and the elastic displacement related to the total energy functional. The multi-field couplings and geometrical nonlinearity of the ferromagnetic thin shell are taken into account in the modeling. The general modeling can be further deduced to existing models of the magneto-elasticity and the thermo-elasticity of a ferromagnetic shell and magneto-thermo-elasticity of a ferromagnetic plate, which are coincident with the ones in literature.     

铁磁薄板在磁场中变形和应力的理论分析与实验研究

在二维板壳磁弹性理论基本方程基础上,以二维铁磁薄板磁弹性问题为例,建立了含有10个基本未知量的偏微分方程组,再利用Newmark有限等差式,得到了可以应用DOM方法求解的标准型方程组.求解得到了载流铁磁薄板的位移及应力与各电磁量之间的关系,计算了载流铁板在磁场中的变形和应力.同时进行了载流铁磁性薄板在电磁场中变形和应力的实验研究,介绍了电磁场中铁磁性薄板的实验装置和实验方法,给出了实验数据,并将实验结果与理论计算结果进行了分析对比.     

ANALYSES ON NONLINEAR COUPLING OF MAGNETO-THERMO-ELASTICITY OF FERROMAGNETIC THIN SHELL-Ⅰ: GENERALIZED VARIATIONAL THEORETICAL MODELING

Based on the generalized variational principle of magneto-thermo-elasticity of the ferromagnetic elastic medium, a nonlinear coupling theoretical modeling for a ferromagnetic thin shell is developed. All governing equations and boundary conditions for the ferromagnetic shell are obtained from the variational manipulations on the magnetic scalar potential, temperature and the elastic displacement related to the total energy functional. The multi-field couplings and geometrical nonlinearity of the ferromagnetic thin shell are taken into account in the modeling. The general modeling can be further deduced to existing models of the magneto-elasticity and the thermo-elasticity of a ferromagnetic shell and magneto-thermo-elasticity of a ferromagnetic plate, which axe coincident with the ones in literature.     

Magneto-thermoelastic problem in non-homogeneous isotropic cylinder

The present paper concerns the investigation of the stress, temperature and magnetic fields in an isotropic elastic cylinder in a primary magenetic field when the curved surface of the cylinder subject to certain boundary conditions.The system of fundamental equations is solved by means of a finite difference method and the numerical calculations are carried out for the temperature, the components of displacement and the components of stresses with time and through the thickness of the cylinder. The results indicate that the effects of inhomogeneity and magnetic field are very pronounced.     

Winkler-Pasternak弹性地基梁自由振动的二维弹性分析

基于二维线弹性理论,应用Hamilton原理,获得Winkler-Pasternak弹性地基梁自由振动的控制微分方程,应用微分求积法(DQM)数值研究了梁自由振动的无量纲频率特性。计算结果与已有的结果(Bernoulli-Euler梁和Timoshenko梁)比较表明,本文的分析方法对弹性地基长梁和短梁自由振动的研究都有效。最后考虑了几何参数对梁频率的影响,以及不同边界条件下地基系数对频率的影响和收敛性。