空间杆系与叠层板组合结构的非线性稳定分析

本文用非线性有限元法研究一种空间杆系与复合材料叠层板组合结构的非线性稳定问题，将杆系离散成空间梁柱单元，将复合材料叠层板离散成复合材料叠层板单元，并考虑了几何非线性和板的成层正交各向异性性质。为了获得较高的计算效率并顺利通过稳定临界点，求解非线性平衡方程采用了变增量步长和柱面弧长法迭代策略。在考证了程序的正确性之后，具体分析了几个工程模型，得出了于工程设计有益的结论。     

钢筋混凝土框架—剪力墙结构非线性抗震分析的一种空间力学模型

将钢筋混凝土框架－剪力墙结构离散为能模拟梁、柱、墙抗震性能的单元，采用杆系－层间模型进行结构的非线性抗震分析。本文的空间力学模型可考虑局部楼板变形对结构地震反应的影响，并可沿任意角度输入相互垂直的两个地震动水平分量，适用于框架一剪力墙复杂结构的非线性抗震分析。     

基于有限单元柔度法的材料与几何双重非线性空间梁柱单元

从有限单元柔度法的基本思想出发,基于完全拉格朗日格式(TL格式),建立了能考虑材料与几何双重非线性的一般化空间梁柱单元,适用于满足Euler-Bernoulli梁柱二阶分析理论假定的空间杆系结构非线性分析。另外,在梁柱单元截面分析中引入纤维模型,使其适用于解决钢筋混凝土这类复合材料结构的非线性分析问题。对钢筋混凝土双向偏心受压柱试验结果的模拟分析表明,本文所提方法是正确、可靠的,能有效地分析钢筋混凝土框架柱的材料与几何双重非线性问题。     

圆柱坐标系中弹性力学的哈密顿正则方程及其哈密顿单元方法

给出圆柱坐标系中弹性力学的哈密顿形式体系，用有限元法得到了哈密顿单元的半解析解，可用于求解复合材料叠层圆柱曲板问题。     

圆柱坐标系中弹性力学的哈密顿正刚方程及其哈密顿单元…

给出圆柱坐标系中弹性力学的哈密顿形式体系，用有限元法得到了哈密顿单元的半解析解，可用于求解复合材料叠层圆柱曲板问题。     

空间杆系结构实用几何非线性分析

从简单实用的角度论述了空间杆系结构的几何非线性分析方法。文中分析了非线性有限元方法的求解过程,特别强调决定几何非线性收敛结果的关键问题,即由节点位移增量计算单元的内力增量。通过引入转随转坐标系,论述了平面和空间梁单元小应变时单元内力增量的计算问题。针对杆系结构的大应变问题,从有限应变理论出发进行分析,提出了对该问题的有效处理方法,并且用实例进行了验证。计算结果表明,该实用几何非线性分析方法是可靠和有效的。     

复合材料叠层板单向拉压强度分析

1. 前言复合材料叠层板是一种非均匀的各向异性材料,其铺层方式及顺序可按外载的需要进行设计,因此每种叠层板的强度特性是不同的。进行复合材料结构的强度分析时,需要知道叠层板元件的强度特性。这些特性可以从试验获得,但如能在给定单向层板特性的基础上,用计算的方法较准确的算出各种叠层元件的强度特性,这将给结构设计及分析带来很大的方便。叠层板是由很多正交各向异性的单层(lamina)按各种不同的方法叠压而成,因此,叠层板本身就是一个复杂的“结构”,加之它的破坏是个逐级破坏的过程,破坏机理又比较     

刚接与铰接混合连接杆系结构的几何非线性分析

本文提出用子结构原理解决具有刚接与铰接混合连接空间杆系结构的几何非线性分析,实现其非线性稳定性分析的载荷-位移全过程跟踪。该法无须单独推导刚接、铰接以及一端刚接一端铰接单元的弹性刚度矩阵和几何刚度矩阵,而可以直接由空间梁单元退化得到,而且可以将平面问题与空间问题、刚接与铰接混合连接体系进行统一处理,算例表明,本文方法对于杆系结构的统一和整体分析是有效的。     

基于Layerwise离散层理论的复合材料夹层板屈曲分析

复合材料夹层结构由于面板和芯层力学特性差异较大,屈曲分析时要分层考虑各层的剪切变形。基于Reddy的Layerwise离散层理论,假设每一层变形服从一阶剪切变形理论,在统一的位移场描述下,推导建立了一种用于复合材料夹层结构屈曲分析的四节点四边形板单元,并采用混合插值方法对单元的剪切锁定进行了修正。分别对三种典型的夹层板结构进行线性屈曲有限元分析,并将计算结果与文献中已有结果进行了对比。结果表明：本文的分析方法能离散考虑各层的力学特性,将结构离散为多层时,计算结果与三维弹性理论或高阶板理论吻合;将结构等效为单层时,计算结果与基于一阶剪切变形理论的文献结构吻合,验证了单元的有效性。     

超级有限元法在桁架组合结构分析中的应用

本文针对复杂杆系组合结构静动力问题讨论一种简捷有效的数值计算方法——超级有限元法,首先介绍其一般原理,继而针对一大类复杂桁架空间杆系组合结构准三维化处理(用一维变量加以表征),并将其离散成一系列单元,考虑弯曲、剪切、挤压、拉伸压缩、扭转等多种非经典变形效应,通过采用构件端部自由度向超级元自由度的转换而把复杂多构件问题的求解变为少量一维结点变量问题的求解,既达到了简化的目的,又保证了精度,而且可方便地与通常的有限元法结合使用。文中还给出了有关杆系组合结构分析的数值算例。     

Nonlinear vibration analysis of isotropic cantilever plate with viscoelastic laminate

The nonlinear vibration of an isotropic cantilever plate with viscoelastic laminate is investigated in this article. Based on the Von Karman’s nonlinear geometry and using the methods of multiple scales and finite difference, the dimensionless nonlinear equations of motion are analyzed and solved. The solvability condition of nonlinear equations is obtained by eliminating secular terms and, finally, nonlinear natural frequencies and mode-shapes are obtained. Knowing that the linear vibration of this type of plate does not have exact solution, Ritz method is employed to obtain semi-analytical nonlinear mode-shapes of transverse vibration of this plate. Airy stress function and Galerkin method are employed to reduce nonlinear PDEs into an ODE of duffing type. Stability of plate and chaotic behavior are investigated by Runge–Kutta method. Poincare section diagrams are in good agreement with results of Lyapunov criteria.     

A coupled interface-multilayer approach for mixed mode delamination and contact analysis in laminated composites

An accurate laminate model developed by using multilayer shear deformable plate modeling and interface elements, based on fracture mechanics and contact mechanics, is proposed to analyze mixed mode delamination in composite laminates. Perfect adhesion along the undelaminated portion of the delamination plane is simulated by treating interface stiffnesses as penalty parameters, whereas to enforce interface displacement continuity between plate elements constituting each sub-laminate above or below the delamination plane, the Lagrange multiplier method is used. The governing differential equations are derived through a variational procedure by using a modified total potential energy functional. Results are obtained by numerical integration of the non-linear three-point boundary value problem modeling mixed-mode delamination of the laminate plate subjected to end loading, which accounts also for the frictionless contact condition.The coupling of a penalty procedure with the Lagrange multiplier method, results in an accurate and direct energy release rate evaluation. Comparisons with results available from the literature obtained with a local continuum approach, show that mode partition may be performed to the desired accuracy by refining multilayer plate models for each sub-laminate. In addition, original analytical formulas for mode partition are obtained by coupling the interface approach and fracture mechanics concepts, evidencing the effectiveness of the proposed approach and gaining a better insight into the influence of shear effects on mode decomposition. Numerical computations for practical problems, evidence both the relative simplicity and the efficiency of the proposed model to represent mixed mode interlaminar fracturing as well as crack–face interaction.     

Time-varying nonlinear dynamics of a deploying piezoelectric laminated composite plate under aerodynamic force

Using Reddy’s high-order shear theory for laminated plates and Hamilton’s principle, a nonlinear partial differential equation for the dynamics of a deploying cantilevered piezoelectric laminated composite plate, under the combined action of aerodynamic load and piezoelectric excitation, is introduced. Two-degree of freedom (DOF) nonlinear dynamic models for the time-varying coefficients describing the transverse vibration of the deploying laminate under the combined actions of a first-order aerodynamic force and piezoelectric excitation were obtained by selecting a suitable time-dependent modal function satisfying the displacement boundary conditions and applying second-order discretization using the Galerkin method. Using a numerical method, the time history curves of the deploying laminate were obtained, and its nonlinear dynamic characteristics, including extension speed and different piezoelectric excitations, were studied. The results suggest that the piezoelectric excitation has a clear effect on the change of the nonlinear dynamic characteristics of such piezoelectric laminated composite plates. The nonlinear vibration of the deploying cantilevered laminate can be effectively suppressed by choosing a suitable voltage and polarity.     

Corotational nonlinear analyses of laminated shell structures using a 4-node quadrilateral flat shell element with drilling stiffness

A new 4-node quadrilateral flat shell element is developed for geometrically nonlinear analyses of thin and moderately thick laminated shell structures. The fiat shell element is constructed by combining a quadrilateral area co- ordinate method （QAC） based membrane element AGQ6- II, and a Timoshenko beam function （TBF） method based shear deformable plate bending element ARS-Q12. In order to model folded plates and connect with beam elements, the drilling stiffness is added to the element stiffness matrix based on the mixed variational principle. The transverse shear rigidity matrix, based on the first-order shear deformation theory （FSDT）, for the laminated composite plate is evaluated using the transverse equilibrium conditions, while the shear correction factors are not needed. The conventional TBF methods are also modified to efficiently calculate the element stiffness for laminate. The new shell element is extended to large deflection and post-buckling analyses of isotropic and laminated composite shells based on the element independent corotational formulation. Numerical re- sults show that the present shell element has an excellent numerical performance for the test examples, and is applicable to stiffened plates.     

Coupled buckling and postbuckling analysis of active laminated piezoelectric composite plates

A theoretical framework for analyzing the pre- and postbuckling response of composite laminates and plates with piezoactuators and sensors is presented. The mechanics include nonlinear effects due to large rotations and stress stiffening, and are incorporated into a coupled mixed-field piezoelectric laminate theory. Using the previous mechanics, a nonlinear finite element method and an incremental-iterative solution are formulated for the analysis of nonlinear adaptive plate structures subject to in-plane electromechanical loading. A novel eight-node nonlinear plate finite element is also developed. Evaluation cases predict the buckling and postbuckling response of adaptive composite beams and plates with piezoelectric actuators and sensors. The case of piezoelectric buckling and postbuckling induced by the actuators is addressed and quantified. Finally, the possibility to actively mitigate the mechanical buckling and postbuckling response of adaptive piezocomposite plates is illustrated.     

埋入压电元件的自适应层合板形状控制的数值分析

把压电元件埋入复合材料层合板中，层合板就变成压电自适应层合板，它除了具有承载能力外，还具有检测、动作、变形等功能。本文用有限元法对压电自适应层合板进行了分析，压电元件的驱动效应被等效为力学载荷。计算过程因此得到了简化，所得算例１的计算结果得到了实验结果的验证。在此基础上，利用非线性规划方法，对结构形状的最优控制进行了分析，并给出了算例２的分析结果。     

破片冲击作用下舰船复合材料结构损伤的近场动力学模拟

基于近场动力学方法，综合分析了破片的速度、层合板的铺层方式、加筋板的筋条尺寸和破片相对筋条的冲击位置对结构损伤模式和破片剩余速度的影响。结果显示:高速破片冲击作用下，层合板会发生侵彻和穿透现象，层合板的损伤模式以基体损伤为主，且随着破片冲击速度的增大，板上下表面的损伤区域呈现出一种先增大后减小的趋势；高速破片冲击作用下，层合的板损伤扩展方向和纤维铺设方向有关，对于纤维铺层方向相同的层合板，其上下表面的损伤扩展方向一般与纤维方向相同；加筋板通过增加少量质量可以获得比层合板更好的抗破片冲击性能，且加筋板的筋条尺寸和破片相对筋条的冲击位置对加筋板的损伤具有明显影响。     

亚音速气流中复合材料悬臂板的非线性振动响应研究

随着材料科学的发展,越来越多的新型材料应用到了工程实践中.在气流激励的作用下,对于以航空航天工程为背景、采用复合材料的板壳结构的非线性动力学问题仍是动力学领域的研究热点.本文研究了复合材料悬臂板在亚音速气流条件下的非线性振动和响应.根据理想不可压缩流体的流动条件和 Kutta--Joukowski升力定理,基于升力面理论,利用涡格法计算了三维有限长平板机翼上的亚音速气动升力.将亚音速气动力施加到复合材料悬臂板上,利用Hamilton原理,考虑Reddy三阶剪切变形理论并引入冯$\cdot$卡门非线性应变位移关系,建立了有限长平板的非线性动力学微分方程.利用有限元方法考察了不同几何参数下层合板悬臂板的固有特性,通过比较不同材料和几何参数的线性系统的固有频率,得到不同比例的内共振关系.利用Galerkin方法将偏微分方程截断为两自由度非线性常微分方程,在这里考虑了1:2的内部共振关系并利用多尺度法进行了摄动分析.对应多个选取参数,得到了频率响应曲线.结果展示了硬化弹簧型行为和跳跃现象.      

Nonlinear buckling optimization of composite structures considering “worst” shape imperfections

Nonlinear buckling optimization is introduced as a method for doing laminate optimization on generalized composite shell structures exhibiting nonlinear behaviour where the objective is to maximize the buckling load. The method is based on geometrically nonlinear analyses and uses gradient information of the nonlinear buckling load in combination with mathematical programming to solve the problem. Thin-walled optimal laminated structures may have risk of a relatively high sensitivity to geometric imperfections. This is investigated by the concepts of “worst” imperfections and an optimization method to determine the “worst” shape imperfections is presented where the objective is to minimize the buckling load subject to imperfection amplitude constraints. The ability of the nonlinear buckling optimization formulation to solve the laminate problem and determine the “worst” shape imperfections is illustrated by several numerical examples of composite laminated structures and the application of both formulations gives useful insight into the interaction between laminate design and geometric imperfections.     

Studies of nonlinear deformation and stability of noncircular cylindrical shells in transverse bending

The variational finite element method in displacements is used to solve the problem of geometrically nonlinear deformation and stability of cylindrical shells with a noncircular contour of the cross-section. Quadrangle finite elements of shells of natural curvature are used. In the approximations of element displacements, the displacements of elements as solids are explicitly separated. The variational Lagrange principle is used to obtain a nonlinear system of algebraic equations for the unknown nodal finite elements. The system is solved by the method of successive loadings and by the Newton-Kantorovich linearization method. The linear system is solved by the Crout method. The critical loads are determined in the process of solving the nonlinear problem by using the Sylvester stability criterion. An algorithm and a computer program are developed to study the problem numerically. The nonlinear deformation and stability of shells with oval and elliptic cross-sections are investigated in a broad range of variation of the elongation and ellipticity parameters. The shell critical loads and buckling modes are determined. The influence of the deformation nonlinearity, elongation, and ellipticity of the shell on the critical loads is examined.