Semi-analytical solution of magneto-thermo-elastic stresses for functionally graded variable thickness rotating disks

In this paper, a semi-analytical solution for magneto-thermo-elastic problem in functionally graded (FG) hollow rotating disks with variable thickness placed in uniform magnetic and thermal fields is presented. Stresses and perturbation of magnetic field vector in FG rotating disks are determined using infinitesimal theory of magneto-thermo-elasticity under plane stress conditions. The material properties except Poisson’s ratio are modeled as power-law distribution of volume fraction. The profile of disk thickness is assumed to be a parabolic function of radius. The non-dimensional distribution of temperature, displacement, stresses and perturbation of magnetic field vector throughout radius are shown. Effects of material grading index, geometry of the disk and magnetic field on the stress and displacement fields are investigated. The results of stresses and radial displacements for two different boundary conditions with and without the effect of magnetic field are compared for a FG rotating disk with concave thickness profile. It has been found that imposing a magnetic field significantly decreases tensile circumferential stresses. Therefore the fatigue life of the disk will be significantly improved by applying the magnetic field. Results of this investigation could be applied for optimum design of FG hollow rotating disks with variable thickness.     

Vibration analysis of piezoelectric FGM sensors using an accurate method

In this study, free vibration analysis of moderately thick smart FG annular/circular plates with different boundary conditions is presented on the basis of the Mindlin plate theory. This structure comprised a host FG plate and two bonded piezoelectric layers. Piezoelectric layers are open circuit therefore this plate can be used as a sensor. According to power-law distribution of the volume fraction of the constituents, material properties vary continuously through the thickness of host plate while Poisson's ratio is set to be constant. Using Hamilton's principle and Maxwell electrostatic equation yields six complex coupled equations which are solved via an exact closed-form method. The accuracy of the frequencies is verified by the available literature, finite element method (FEM) and the Kirchhoff theory. The effects of plate parameters like boundary condition and gradient index are investigated and significance of coupling between in-plane and transverse displacements on the resonant frequency is proved.     

斜直井中钻柱非线性屈曲分析的有限元增量加权迭代法

建立了直井中有重钻柱非线性屈曲平衡方程及相应的泛函表达式，构造了一种能避免发散并加速收敛的有限元增量加权迭代法求解斜直井中钻柱屈曲问题，首次采用直接求解非线性特征值的方法来分析钻柱屈曲问题。计算结果表明，井斜角和屈曲位移幅值中任意一项增大，钻柱的非线性屈曲临界载荷也随之增大，且呈非线性；重力线密度对屈曲载荷的影响与屈曲状态有关；井壁约束力和特征向量的幅值随屈曲位移幅值增大而增大；扭矩对屈曲的影响很小，可以略去。研究结果为实际钻井作业中屈曲载荷的确定提供了参考。     

纤维增强FGM梁的自由振动和临界屈曲载荷分析

基于经典梁理论（CBT）研究轴向力作用下纤维增强功能梯度材料（FGM）梁的横向自由振动和临界屈曲载荷问题。首先考虑由混合律模型来表征纤维增强FGM梁的材料属性,其次利用Hamilton原理推导轴向力作用下纤维增强FGM梁横向自由振动和临界屈曲载荷的控制微分方程,并应用微分变换法（DTM）对控制微分方程及边界条件进行变换,计算了纤维增强FGM梁在固定-固定（C-C）、固定-简支（C-S）和简支-简支（S-S）3种边界条件下横向自由振动的无量纲固有频率和无量纲临界屈曲载荷。退化为各向同性梁和FGM梁,并与已有文献结果进行对比,验证了本文方法的有效性。最后讨论在不同边界条件下纤维增强FGM梁的刚度比、纤维体积分数和无量纲压载荷对无量纲固有频率的影响以及各参数对无量纲临界屈曲载荷的影响。     

受扭圆轴弹性屈曲分析时边界条件的不一致性

受扭圆轴广泛存在于机械结构中。对于细长圆轴,当其两端受扭矩作用时将发生弹性屈曲。在两端简支条件下,文献中已有的解析解和基于能量法的数值近似解结果不吻合,原因是人们没有注意到边界条件的不一致性。笔者用有限元法进行了分析,给出了单元的线性刚度矩阵和增量刚度矩阵。分析发现:由能量变分方程所得到的应自动满足的自然边界条件———力边界条件和得到解析解的二阶平衡微分方程所应满足的力边界条件两者在简支情况下是不一致的。考虑到与解析解的力边界条件的等效性后用有限元数值分析方法得到的结果与解析解极为吻合。有限元解与解析解间的差异,有时并非单元性能的原因或计算误差造成的。     

Buckling analysis of FGM plates under uniform,linear and non-linear in-plane loading

The paper investigates the buckling responses of functionally graded material (FGM) plate subjected to uniform, linear, and non-linear in-plane loads. New nonlinear in-plane load models are proposed based on trigonometric and exponential function. Non-dimensional critical buckling loads are evaluated using non-polynomial based higher order shear deformation theory. Navier’s method, which assures minimum numerical error, is employed to get an accurate explicit solution. The equilibrium conditions are determined utilizing the principle of virtual displacements and material property are graded in the thickness direction using simple Voigt model or exponential law. The present formulation is accurate and efficient in analyzing the behavior of thin, thick and moderately thick FGM plate for buckling analysis. It is found that with the help of displacement-buckling load curve, critical buckling load can be derived and maximum displacement due to the instability of inplane load can be obtained. Also, the randomness in the values of transverse displacement due to inplane load increases as the extent of uniformity of the load on the plate is disturbed. Furthermore, the parametric varying studies are performed to analyse the effect of span-to-thickness ratio, volume fraction exponent, aspect ratio, the shape parameter for non-uniform inplane load, and non-dimensional load parameter on the non-dimensional deflections, stresses, and critical buckling load for FGM plates.

     

Exact static analysis of partially composite beams and beam-columns

The ordinary differential equations and general solutions for the deflection and internal actions and, especially, the pertaining consistent boundary conditions for partially composite Euler–Bernoulli beams and beam-columns are presented. Static loading conditions, including transverse and axial loading and first- and second-order analyses are considered. The theoretical procedure is applicable to general loading and boundary conditions for uniform composite beams and beam-columns with interlayer slip. Further, the exact closed form characteristic equations and their associated exact buckling length coefficients for composite columns with interlayer slip are derived for the four Euler boundary conditions. It is shown that these coefficients are the same as those for ordinary fully composite (solid) columns, except for the Euler clamped-pinned case. For the clamped-pinned case, the difference between the exact buckling length coefficient and the corresponding value for solid columns is less than 1.8% depending on the so-called composite action parameter and relative bending stiffness parameter. Correspondingly, the maximum deviation between the exact and approximate buckling load is at most 2.5%. These small differences can in most practical cases be neglected. Also, the maximum theoretical range for the relative bending stiffness for partially composite beams and beam-columns is derived. An effective bending stiffness, valuable in the determination of the critical buckling load for partially composite members, is derived. This effective bending stiffness is also suitable for analysing approximate deflections and internal actions or stresses in composite beams with flexible shear connection. The beam-column analysis is applied to a specific case. The difference in the approaches to the first- and second-order analysis is illustrated and the results clearly show the magnification in the actions and displacements due to the second-order effect. The magnification of the internal axial forces is different from magnifications obtained for the other internal actions, since only that portion of an internal axial force that is induced by bending is magnified by the second-order effect.     

Buckling of rectangular plates with various boundary conditions loaded by non-uniform inplane loads

In the present paper, buckling loads of rectangular composite plates having nine sets of different boundary conditions and subjected to non-uniform inplane loading are presented considering higher order shear deformation theory (HSDT). As the applied inplane load is non-uniform, the buckling load is evaluated in two steps. In the first step the plane elasticity problem is solved to evaluate the stress distribution within the prebuckling range. Using the above stress distribution the plate buckling equations are derived from the principle of minimum total potential energy. Adopting Galerkin's approximation, the governing partial differential equations are converted into a set of homogeneous linear algebraic equations. The critical buckling load is obtained from the solution of the associated linear eigenvalue problem. The present buckling loads are compared with the published results wherever available. The buckling loads obtained from the present method for plate with various boundary conditions and subjected to non-uniform inplane loading are found to be in excellent agreement with those obtained from commercial software ANSYS. Buckling mode shapes of plate for different boundary conditions with non-uniform inplane loadings are also presented.     

Three-dimensional analysis of thermal stresses in four-parameter continuous grading fiber reinforced cylindrical panels

In this paper, three-dimensional analysis of thermal stresses in four-parameter continuous grading fiber reinforced cylindrical panel subjected to thermal load is studied. The cylindrical panel is assumed to be made of an orthotropic material. The continuous grading fiber reinforced panel has a smooth variation in matrix volume fraction in the radial direction. A generalization of the power-law distribution presented in literature is proposed. Symmetric and asymmetric volume fraction profiles are presented in this paper. Suitable temperature and displacement functions that identically satisfy the boundary conditions at the simply supported edges are used to reduce the thermoelastic equilibrium equations to a set of coupled ordinary differential equations with variable coefficients, which can be solved by generalized differential quadrature method. The fast rate of convergence of the method is demonstrated and to validate the results, comparisons are made with the available solutions for orthotropic shells. The main contribution of this work is to illustrate the influence of the power-law exponent, of the power-law distribution choice and of the choice of the four parameters on the thermal behaviour of continuous grading fiber reinforced cylindrical panels.     

Exact piezothermoelastic solution of simply-supported orthotropic flat panel in cylindrical bending

Exact piezothermoelastic solution is presented for infinitely long, simply-supported, orthotropic, piezoelectric, flat panel in cylindrical bending under pressure, thermal and electrostatic excitation. The displacements, electric potential and temperature are expanded in Fourier series to satisfy the boundary conditions at the longitudinal edges. The governing equations reduce to ordinary differential equations with constant coefficients. The six constants in their general solution for each Fourier component are obtained from the boundary conditions at the lateral surfaces. The solution of the inverse problem of inferring the applied temperature from the measured electrical potential difference between the lateral surfaces of the panel is also presented. Results are presented for sinusoidal thermal and electrostatic loadings and the effect of width to thickness ratio has been investigated.     

Buckling analysis of laminated plates using the extended Kantorovich method and a system of first-order differential equations

The extended Kantorovich method using multi-term displacement functions is applied to the buckling problem of laminated plates with various boundary conditions. The out-of-plane displacement of the buckled plate is written as a series of products of functions of parameter x and functions of parameter y. With known functions in parameter x or parameter y, a set of governing equations and a set of boundary conditions are obtained after applying the variational principle to the total potential energy of the system. The higher order differential equations are then transformed into a set of first-order differential equations and solved for the buckling load and mode. Since the governing equations are first-order differential equations, solutions can be obtained analytically with the out-of-plane displacement written in the form of an exponential function. The solutions from the proposed technique are verified with solutions from the literature and FEM solutions. The bucking loads correspond very well to other available solutions in most of the comparisons. The buckling modes also compare very well with the finite element solutions. The proposed solution technique transforms higher-order differential equations to first-order differential equations, and they are analytically solved for out-of-plane displacement in the form of an exponential function. Therefore, the proposed solution technique yields a solution which can be considered as an analytical solution.     

热-机荷载作用的P-FGM扁球壳跳跃屈曲分析

基于经典壳体理论和Sanders非线性应变-位移关系,导出了幂律型功能梯度材料(P-FGM)扁球壳在热-机械荷载作用下的几何非线性常微分控制方程。推导过程考虑了沿厚度存在一维热传导温度场和法向均布荷载作用。采用打靶法求解了由控制方程和固定夹紧边界条件构成的两点边值问题。得到了FGM扁球壳的一些典型的屈曲平衡路径和双稳态构形。对热-机械荷载作用的FGM扁球壳的跳跃屈曲行为进行了参数影响分析。结果表明：温度上升时,球壳上临界荷载显著增加、下临界荷载变化不明显。梯度指数增加时,球壳上、下临界荷载均显著减小。组分材料模量增加时,球壳上、下临界荷载均显著增加。当底圆半径和厚度给定时,随壳体中面曲率半径增加,球壳上、下临界荷载均显著增加。当中面曲率半径和厚度给定时,随底圆半径增加,球壳下临界荷载显著减小,上临界荷载几乎不变。     

Investigating the thermal environment effects on geometrically nonlinear vibration of smart functionally graded plates

An analytical solution for a sandwich circular FGM plate coupled with piezoelectric layers under one-dimensional heat conduction is presented. All materials of the device may be of any functional gradients in the direction of thickness. The solution exactly satisfies all the equilibrium conditions and continuity conditions for the stress, displacement and electric displacement as well as electric potential on the interfaces between adjacency layers. A nonlinear static problem is solved first to determine the initial stress state and pre-vibration deformations of the FG plate that is subjected to in-plane forces and applied actuator voltage in thermal environment in the case of simply supported boundary conditions. By adding an incremental dynamic state to the pre-vibration state, the differential equations that govern the nonlinear vibration behavior of pre-stressed piezoelectric coupled FGM plates are derived. The role of thermal environment as well as control effects on nonlinear static deflections and natural frequencies imposed by the piezoelectric actuators using high input voltages are investigated. Numerical examples are provided and simulation results are discussed. Numerical results for FGM plates with a mixture of metal and ceramic are presented in dimensionless forms. The good agreement between the results of this paper and those of the finite element (FE) analyses validated the presented approach. In a parametric study the emphasis is placed on investigating the effect of varying the applied actuator voltage and thermal environment as well as gradient index of FG plate on the dynamics and control characteristics of the structure.     

Stochastic hybrid numerical method for transient analysis of stress field in functionally graded thick hollow cylinders subjected to shock loading

This investigation aims to study the random stresses in a functionally graded (FG) thick hollow cylinder with uncertain material properties subjected to mechanical shock loading using a hybrid numerical method. The mechanical properties are considered to vary across thickness of FG cylinder as a nonlinear power function of radius. The stresses are obtained by solving Navier equation and using Galerkin finite element and Newmark finite difference methods. The Monte Carlo simulation is used to generate the random mechanical properties for the problem. The failure probabilities and time history analysis of stresses are determined for various coefficient of variation considering various grading patterns of mechanical properties. The presented hybrid numerical method is effective, with high capability for stochastic analysis of dynamic and transient analysis of FG structures with various boundary conditions.     

水平井中钻柱屈曲的非线性有限元分析

建立了直井中钻柱屈曲的平衡方程及对应的泛函表达式,使用有限元法对不同约束下水平井中钻柱从稳定到非线性屈曲的整个过程进行了分析。力学模型中考虑了重力和扭矩对屈曲的影响。分析结果表明:钻柱的屈曲是一个从局部屈曲到总体屈曲的过程,屈曲段的井壁约束力、钻柱弯矩和屈曲位移呈周期性变化;重力在水平井中对屈曲有较强的抑制作用;边界约束对屈曲的影响不可忽略;扭矩对屈曲的影响很小,可以忽略。     

Pseudospectral analysis of buckling and free vibration of rectangular Mindlin plates

A study of buckling and free vibration of rectangular Mindlin plates is presented. The analysis is based on the pseudospectral method, which uses basis functions that satisfy the boundary conditions. The equations of motion are collocated to yield a set of algebraic equations that are solved for the critical buckling load and for the natural frequencies in the presence of the in-plane loads. Numerical examples of rectangular plates with SS-C-SS-C boundary conditions are provided for various aspect ratios and thickness ratios, which show good agreement with those of the classical plate theory when the thickness ratio is very small. This paper was recommended for publication in revised form by Associate Editor Eung-Soo Shin Jinhee Lee received B.S. and M.S. degrees from Seoul National University and KAIST in 1982 and 1984, respectively. He received his Ph.D. degree from the University of Michigan, Ann Arbor in 1992 and joined the Dept. of Mechanical and Design Engineering of Hongik University in Choongnam, Korea. His research interests include inverse problems, pseudospectral method, vibration and dynamic systems.     

Studies on buckling behavior of honeycomb sandwich panel

In the present work, an avant-garde study on Buckling of honeycomb sandwich panel of hexagonal cell has been carried out by considering simply supported boundary constraints explicitly. A novel theoretical formulation has been developed for the buckling of honeycomb panel as per Mil Standard (MIL-HDBK-23A). The laminate element type has been considered for linear static analysis using Lanczos Method (finite element method (FEM)). The buckling analysis carried out by FEM has been compared and validated with theoretical and experimental investigation for different panel aspect ratio under static load.     

基于偶应力理论的压杆屈曲载荷的尺寸效应

许多微尺度实验证实,一些材料在微尺度下的力学性能存在尺寸效应.基于修正偶应力理论研究细长压杆屈曲载荷的尺寸效应,利用变分原理得到细长压杆屈曲的控制方程和边界条件,求解相应的屈曲边界值问题.结果表明,当细长压杆的特征尺寸与其材料本征尺度参数相当时,其屈曲载荷显著增大,表现出明显的尺寸效应;而当细长压杆的特征尺寸远大于其材料本征尺度参数时,其屈曲载荷与传统理论值基本相等,即其尺寸效应消失.     

Thermal post-buckling analysis of slender columns using the concept of coupled displacement field

Thermal post-buckling analysis of columns with an axially immovable ends is studied using the Rayleigh-Ritz (R-R) method, where the admissible displacement functions are chosen based on the concept of coupled displacement field (CDF) criteria. Geometric non-linearity is considered using the von-Karman strain displacement relations of the beam. Furthermore, the displacement fields derived from CDF criteria are used in an intuitive formulation, where the thermal post-buckling behavior can be predicted by using two parameters namely tension developed in the column and linear buckling load. An exhaustive set of column boundary conditions are considered namely classical such as hinged-hinged, clamped-clamped, clamped-hinged and non-classical such as clamped-guided and hinged-guided. Post-buckling analysis results are presented in the form of closed form expressions, where the ratio of post-buckling load to linear buckling load parameter is expressed as a function of central amplitude of the column for all the boundary conditions considered. The amount of non-linearity predicted using the present formulations (R-R method and intuitive method) based on the concept of coupled displacement field (CDF) criteria shows an excellent agreement with the available literature results for both classical and non-classical boundary conditions.     

The analytical solution of the buckling of composite truncated conical shells under combined external pressure and axial compression?

The objective of this research is determining the buckling load of composite truncated conical shells under external loading by theoretical and numerical methods. The boundary conditions are assumed to be clamped. At first, basic equations and stability relations of conical shells were derived. The analysis is carried out using Donnel-type stability equations for thin cross-ply conical shells. By applying Galerkin??s method, these equations are converted to a system of ordinary time dependent differential equations. Ritz method is employed for finding the dynamic stability load. Finally, the critical static and dynamic buckling loads and the corresponding wave numbers have been found analytically. Then comparison of results is considered. Results of analytical calculations are compared with numerical results and with other researchers?? analytical results. The effects of geometric parameters, the cone semi-vertex angle, number of layers and material of fibers on buckling loads are discussed.