含初应力复合材料柱壳结构的双稳态特性

为分析初应力对复合材料圆柱壳结构双稳态特性的影响,采用经典板壳理论建立复合材料圆柱壳力学模型,基于层合结构本构关系推导用双参数表达的系统应变能公式;根据最小势能原理得到双稳态产生的条件和稳态时的曲率表达式。利用Abaqus软件构建圆柱壳的有限元模型,通过附加边界弯矩对柱壳稳态跃迁过程进行模拟。理论计算结果与有限元结果的对比验证理论模型的正确性。分析结果表明：当初应力满足一定条件时,复合材料柱壳结构在其变形过程中有2个稳定平衡位置,并且在稳定平衡位置结构都不产生扭转变形;2个稳定平衡位置的曲率方向可以相同或相反,这与无初应力时反对称复合材料柱壳双稳态曲率方向只能相同的情况有区别。     

A unified approach to nonlinear buckling optimization of composite structures

A unified approach to nonlinear buckling fiber angle optimization of laminated composite shell structures is presented. The method includes loss of stability due to bifurcation and limiting behaviour. The optimization formulation is formulated as a mathematical programming problem and solved using gradient-based techniques. Buckling of a well-known cylindrical shell benchmark problem is studied and the solutions found in literature are proved to be incorrect. The nonlinear buckling optimization formulation is benchmarked against the traditional linear buckling optimization formulation through several numerical optimization cases of a composite cylindrical shell panel which clearly illustrates the advantage and potential of the presented approach.     

Failure analysis of laminated composite cylindrical/spherical shell panels subjected to low-velocity impact

A 4-noded, 48 d.o.f. doubly curved quadrilateral shell finite element based on Kirchhoff–Love shell theory, is used in the nonlinear finite element analysis to predict the damage of laminated composite cylindrical/spherical shell panels subjected to low-velocity impact. The large displacement stiffness matrix is formed using Green's strain tensor based on total Lagrangian approach. An incremental/iterative scheme is used for solving resulting nonlinear algebraic equations by Newton–Raphson method. The damage analysis is performed by applying Tsai–Wu quadratic failure criterion at all Gauss points and the mode of failure is identified using maximum stress criteria. The modes of failure considered are fiber breakage and matrix cracking. The progressive failure analysis is carried out by degrading the stiffness of the material suitably at all failed Gauss points. The load due to low-velocity impact is treated as an equivalent quasi-static load and Hertzian law of contact is used for finding the maximum contact force. After evaluating the nonlinear finite element analysis thoroughly for typical problems, damage analysis was carried out for cross-ply and quasi-isotropic cylindrical/spherical shell panels.     

Instability of short stiffened composite cylindrical shells under bending with prebuckling displacements

This paper presents results for buckling of a stiffened cylindrical shell with cutouts and both isotropic and composite shells without cutouts acting under end bending moments. The STAGS-C program has been used in the analysis.     

A Ritz procedure for optimisation of cylindrical shells, formed by a nearly symmetric and balanced angle-ply composite laminate, with fixed minimum frequency

The paper deals with the problem of optimisation of a cylindrical shell profile under a frequency constraint. The minimum value of the thickness has been established a priori. The structure considered is typical of aerospace craft vessels.The same value of the lowest vibration frequency of the reference cylindrical shell with uniform thickness, has been imposed. That is the minimization procedure of the structure weight must not affect its lowest vibration frequency.Instead of the currently applied finite element method (FEM), Ritz series expansions have been utilized in the analytical developments both for the dynamic variables and for the thickness axial distribution over the shell surface. Lagrange multipliers, together with governing equations and objective function, have been utilized to form the Lagrangian functional, as in the classical Euler–Lagrange method. Imposing the stationary conditions with respect to the Lagrangian degrees of freedom gives a non-linear algebraic equations system, whose solution can be found with an appropriate algorithm.A series of repeated optimisation operations have been performed to arrive at the minimized weight profile, but with the pre-established minimum value of the shell thickness.A simplified nearly symmetric and balanced multilayer composite angle-ply laminate of the shell structure is supposed, as in the case of the uniform thickness reference shell, previously considered for the dynamic analysis. Significant results of some computation application cases can be helpful to evaluate the efficiency of the proposed optimisation procedure applied to cylindrical structures.     

Postbuckling of multilayer cylindrical and spherical shell panels reinforced with graphene platelet by isogeometric analysis

The present work fills a gap on the postbuckling behavior of multilayer functionally graded graphene platelet reinforced composite (FG-GPLRC) cylindrical and spherical shell panels resting on elastic foundations subjected to central pinching forces and pressure loadings. Based on a higher-order shear deformation theory and the von Kármán’s nonlinear strain–displacement relations, the governing equations of the FG-GPLRC cylindrical and spherical shell panels are established by the principle of virtual work. The non-uniform rational B-spline (NURBS) based isogeometric analysis (IGA), the modified arc-length method and the Newton’s iteration method are employed synthetically to obtain nonlinear load–deflection curves for the panels numerically. Several comparative examples are performed to test reliability and accuracy of IGA and arc-length method in present formulation and programming implementation. Parametric investigations are carried out to illustrate the effects of dispersion type of the graphene platelet (GPL), weight fraction of the GPL, thickness of the panel, radius of the panel and parameters of elastic foundation on the load–deflection curves of the FG-GPLRC shell panels. Some complex load–deflection curves of the FG-GPLRC cylindrical and spherical shell panels resting on elastic foundations may be useful for future references.

     

Effect of boundary conditions on natural frequencies for rotating composite cylindrical shells with orthogonal stiffeners

The effect of the boundary conditions on the natural frequencies for rotating composite cylindrical shells with the orthogonal stiffeners is investigated using Love’s shell theory and the discrete stiffener theory. The frequency equation is derived using the Rayleigh–Ritz procedure based on the energy method. The considered boundary conditions are four sets, namely: (1) clamped–clamped; (2) clamped–simply supported; (3) clamped–sliding; and (4) clamped–free. The beam modal function is used for the axial vibration mode and the trigonometric functions are used for the circumferential vibration mode. The composite shells are stiffened with uniform intervals and the stiffeners have the same material. By comparison with the previously published analytical results for the rotating composite shell without stiffeners and the orthogonally stiffened isotropic cylindrical shells, it is shown that natural frequencies can be determined with adequate precision.     

Linear programming and limit analysis of shell roofs

Limit analysis for cylindrical shell roofs has been formulated as a linear programming problem based on lower bound theorem. The differential equations of equilibrium for a circular cylindrical shell element are transformed into algebraic equations by finite differences. The equilibrium equations and the linearized non-linear yield conditions at various points of the shell are linear functions of the stress resultants. These form the linear constraints of the problem. The load parameter is taken as an objective function and it is maximized using revised simplex method. For a shell of given geometry, stress resultants at various points are obtained to give the optimum collapse load. Thus the versatile technique avoids various trial solutions to achieve best lower bound for complicated shell problems.     

Multi-objective optimization of laminated composite shells for minimum mass/cost and maximum buckling pressure with failure criteria under external hydrostatic pressure

In the present study, Multi-objective optimization of composite cylindrical shell under external hydrostatic pressure was investigated. Parameters of mass, cost and buckling pressure as fitness functions and failure criteria as optimization criterion were considered. The objective function of buckling has been used by performing the analytical energy equations and Tsai-Wu and Hashin failure criteria have been considered. Multi-objective optimization was performed by improving the evolutionary algorithm of NSGA-II. Also the kind of material, quantity of layers and fiber orientations have been considered as design variables. After optimizing, Pareto front and corresponding points to Pareto front are presented. Trade of points which have optimized mass and cost were selected by determining the specified pressure as design criteria. Finally, an optimized model of composite cylindrical shell with the optimum pattern of fiber orientations having appropriate cost and mass is presented which can tolerate the maximum external hydrostatic pressure.

     

Finite element analysis of bimodulus composite stiffened thin shells of revolution

This paper is a sequel to the work published by the first and third authors[l] on stiffened laminated shells of revolution made of unimodular materials (materials having identical properties in tension and compression). A finite element analysis of laminated bimodulus composite thin shells of revolution, reinforced by laminated bimodulus composite stiffeners is reported herein. A 48 dot doubly curved quadrilateral laminated anisotropic shell of revolution finite element and it's two compatible 16 dof stiffener finite elements namely: (i) a laminated anisotropic parallel circle stiffener element (PCSE) and (ii) a laminated anisotropic meridional stiffener element (MSE) have been used iteratively.The constitutive relationship of each layer is assumed to depend on whether the fiberdirection strain is tensile or compressive. The true state of strain or stress is realized when the locations of the neutral surfaces in the shell and the stiffeners remain unaltered (to a specified accuracy) between two successive iterations. The solutions for static loading of a stiffened plate, a stiffened cylindrical shell. and a stiffened spherical shell, all made of bimodulus composite materials, have been presented.     

A conical shell finite element

In the analysis of rocket and missiles structures one frequently encounters cylindrical and cornica' shells. A simple finite element which fits the above configuration is obviously a conical shell finite element. In this paper stiffness matrix for a conical shell finite element is derived using Novozhilov's strain-displacement relations for a conical shell. Numerical integration is carried out to ge. the stiffness matrix. The element has 28 degrees-of-freedom and is nonconforming. An eigenvalue analysis of the stiffness matrix showed that it contains all the rigid body modes (six in this case) adequately, which is one of the convergence criteria. An advantage of this element is that a cylindrical shell, an annular segment flat plate, a rectangular flat plate elements can easily be obtained as degenerate cases. The effectiveness of this element is shown through a variety of numerical examples pertaining to annular plate, cylindrical shell and conical shell problems. Comparison of the present solution is made with the existing ones wherever possible. The comparison shows that the present element is superior in some respects to the existing elements     

A geometrically nonlinear shell finite element for tire vibration analysis

An axisymmetric finite element is developed which includes such features as orthotropic material properties, doubly curved geometry, and both the first and second order nonlinear stiffness terms. This element can be used to predict the equilibrium state of an axisymmetric shell structure with geometrically nonlinear large displacements. Small amplitude vibration analysis can then be performed based on this equilibrium state. The nonlinear path is predicted by using the self-correcting incremental procedure and any point on the path can be checked by using the Newton-Raphson iterative scheme. The present formulation and solution procedure are evaluated by analyzing a series of examples with results compared with alternative known solutions. Examples include: free vibration of an isotropic cylindrical shell, a conical frustum, and an orthotropic cylindrical shell; buckling of a cylindrical shell; large deflection of a clamped disk, a spherical cap, and a steel belted radial tire. The final example is a free vibration analysis of the inflated tire and the natural frequencies obtained compared well with published experimental data.     

An exact solution for free vibration of cross-ply laminated composite cylindrical shells with elastic restraint ends

In this paper, a new exact solution based on the classical shell theory (CST) for free vibration of cross-ply laminated composite cylindrical shells with elastic restraint ends is proposed. The present exact solution can be summed up in the following steps: Firstly, the displacement functions are constructed by the governing differential equations with the exact closed form solutions; Then, the artificial spring technology is introduced to simulate the general boundary conditions of the two end edges of shell; Thirdly, the equation for natural frequencies is obtained by means of the method of reverberation-ray matrix (MRRM); Lastly, the vibration results are presented by the modified golden section search (MGSS) algorithm. By comparing the present method with published papers, the accuracy of present method is verified. On the basis of that, some new exact nature frequencies and mode shapes of the cross-ply laminated composite cylindrical shells with various classical boundary conditions and elastic restraints are performed and they can be served as the benchmark data for the future.     

对接圆柱壳结构模态特性分析

对接圆柱壳结构在航空航天、船舶、土木和机械等工程领域得到广泛应用,对其模态特性的分析是研究其动力学特性的重要方向.本文简要介绍了模态分析技术和最小二乘复频域法(PolyMAX)的基本原理,并对对接圆柱壳结构进行了计算模态分析和实验模态分析.实验模态分析过程中,针对自由边界的实验实现、分析结果的正确性进行了讨论,并将实验结果与计算模态分析结果进行对比,对比结果表明：对接圆柱壳结构具有圆柱壳结构一般振动特性的同时,由于对接形式的存在出现了以法兰面为分界的非对称振动.     

Optimum design of composite shells subject to natural frequency constraints

A technique for the optimum (minimum weight) design of a composite shell subject to constraints on its natural frequencies is presented. The optimization problem is posed as a general mathematical programming problem in which one or more of the inequality constraints involves the shell natural frequencies, which must be evaluated numerically during the optimization. For this reason, a method for numerically evaluating the natural frequencies of composite shells is also presented. The method is based upon the finite element method of structural analysis and Rayleigh's principle. Because the element used is applicable to anisotropic shells of arbitrary shape, the method is very general. By using Rayleigh's principle, the necessity of assembling overall mass and stiffness matrices for the shell is eliminated. The optimization is performed by nondimensionalizing the mathematical programming problem and using the penalty function method of Fiacco and McCormick to transform the problem to a sequence of unconstrained minimizations having solutions which converge to the solution of the original (constrained) problem. The unconstrained minimizations are performed using the variable metric method of Fletcher and Powell. Derivatives of the nondimensional frequency constraints are evaluated numerically using difference equations. The frequency calculation method is demonstrated by calculating the fundamental frequency for the transverse vibration mode of a multilayered cylindrical shell with fixed overall geometry and variable composite geometry. Results indicate that the frequency increases with increasing fiber orientation angle, fiber volume fraction, or lamina thickness. The optimization technique is demonstrated by minimizing the weight of the shell discussed above subject to a constraint on its fundamental transverse frequency. The design variables are the fiber orientation angle, the fiber volume fraction, and the lamina thickness. Results are presented and explained in terms of the physical aspects of the problem.     

Optimization of inelastic cylindrical shells with internal supports

A non-linear programming method is developed for optimization of inelastic cylindrical shells with internal ring supports. The shells under consideration are subjected to internal pressure loading and axial tension. The material of shells is a composite which is considered as an anisotropic inelastic material obeying the yield condition suggested by Lance and Robinson. Taking geometrical non/linearity of the structure into account optimal locations of internal ring supports are determined so that the cost function attains its minimum value. A particular problem of minimization of the mean deflection of the shell with weakened singular cross sections is treated in a greater detail.     

Free vibration analysis of multi-symmetric stiffened shells

The application of structural symmetry techniques to the free vibration analysis of cylindrical and conical shells for the prediction of natural frequencies and mode shapes is described. Appropriate boundary conditions have been developed for the analysis of only a part of the shell and have been shown to yield results comparable to the full shell analysis. Half and quarter models of the shell have been developed and analysed using semi-loof and facet shell finite elements. Unstiffened and stiffened circular cylindrical shells and stiffened conical shells have been considered.     

A comparison of analytical–numerical models for nonlinear vibrations of cylindrical shells

The nonlinear flexural vibration behaviour of cylindrical shells has received considerable attention to date. It is pointed out that, although in a well-known reference case there seems to be a reasonable agreement, there are unresolved discrepancies between the results obtained by different authors. In the present paper, the problem is studied using various analytical–numerical models with different levels of accuracy and complexity. The frequency–amplitude curves from the different analysis models developed are compared both for isotropic shells and for an orthotropic composite shell. Secondary modes can play an important role. In more complicated cases modal interactions may significantly influence the nonlinear vibration behaviour, and the results obtained strongly depend on the analysis model chosen.     

基于有限差分法的端部约束同轴圆柱壳的流致失稳分析

本文基于经典壳体和理想势流理论建立圆柱壳流固耦合系统的运动方程,并引入有限差分法(FDM)对运动方程进行离散.为将壳体表面的扰动压力离散到差分网格节点,基于差分离散的原理,提出了以分段函数作为基函数的展开方法.本文对结构与压力控制方程均采用FDM方法进行求解,发展了一种基于FDM的同轴圆柱壳流固耦合求解策略.首先,利用本文方法计算了同轴圆柱壳在静态流体环境中的振动频率,并与有限元软件(ANSYS)的计算结果相比较,验证了本文方法的正确性；然后,探究了在静流体环境中同轴圆柱壳的结构参数变化对其振动频率的影响规律；最后,研究了同轴圆柱壳系统在运动流体环境中的流弹失稳问题.     

On post-buckling analysis and experimental correlation of cylindrical composite shells with Reissner-Mindlin-Von Kármán type facet model

In this paper post-buckling analysis of carbon fibre reinforced plastic cylindrical shells under axial compression is considered. Reissner-Mindlin-Von Kármán type shell facet model is used in the computations. The effect of geometric imperfection shape and amplitude on nonlinear analysis results is discussed. Numerical-experimental correlation is performed using the results of experimental buckling tests found in the literature. Results show that bringing the diamond shape geometric imperfection in the model significantly improves the correlation and gives good accuracy in simulating cylindrical shell post-buckling behaviour.