新型空间薄壁梁单元

基于Timoshenko梁理论和Vlasov薄壁杆件约束扭转理论,建立了具有内部结点的新型空间薄壁截面梁单元．通过对弯曲转角和翘曲角采取独立插值的方法,考虑了横向剪切变形,扭转剪切变形及其耦合作用,弯曲变形和扭转变形的耦合以及二次剪应力等因素影响,由Hellinger-Reissner广义变分原理,推得单元刚度矩阵．算例表明所建模型具有良好的精度,可用于空间薄壁杆系结构的有限元分析．     

基于曲率插值的大变形梁单元

线性梁单元的形函数在单元大转动时会引起虚假应变,不适用于几何非线性分析．传统的几何非线性梁单元由于位移插值和转角插值的相干性,常常引起剪切闭锁等问题．该文 提出了一种平面大变形梁单元,通过单元域内的曲率插值以及曲率与节点位移之间的函数关系,将单元节点力和节点位移表示为节点曲率的函数．由于曲率插值本质上是对梁的应变进行插值,保证了单元任意刚体运动不会产生虚假的节点力;且将梁的截面形心位移表示为曲率的函数,避免了传统单元中的剪切闭锁问题．因而所提方法特别适用于梁的几何非线性分析．数值算例说明了所提方法的正确性和有效性.     

翘曲空间曲线梁自然坐标精确解

基于现有空间曲线梁理论,考虑与扭转有关的翘曲变形和横向剪切变形的影响,建立了自然标架下空间曲线梁的内力和变形的解析解答．将该解答应用于受均布扭矩和竖向分布荷载的平面曲线梁的分析,将所得结果与Heins解答进行比较,证明了理论的正确．并应用该理论分析了解析式中翘曲和横向剪切变形项的影响．     

纤维体积分数可变的复合材料梁计及转动惯量和剪切变形时的固有频率

研究纤维体积分数沿着厚度可变的对称复合材料梁的振动.分析中考虑了一阶剪切变形和转动惯量.该解法可适应任意边界条件.纤维体积分数沿着梁的厚度方向以坐标的m幂次多项式形式连续渐变.可变的纤维体积分数,在对称复合材料梁中形成功能梯度材料(FGM),会引起梁的某些振动特性的改变.结果显示,剪切变形、纤维体积分数和边界条件,对复合材料梁的固有频率和振型的影响.     

一种新的计算Timoshenko梁截面剪切系数的方法

Timoshenko梁理论中考虑了截面剪切变形的影响,推导了一种新的计算剪切系数的方法,首先采用悬臂梁纯弯曲变形条件下截面剪应力分布的精确解,基于能量原理得到了各种梁截面剪切系数新的表达式,然后推导了弯扭耦合变形条件下截面剪应力分布的精确解,进一步获得了该条件下截面的剪切系数.结果表明,悬臂梁端面作用力偏离截面的弯曲中心将使剪切系数变小,通过与Cowper计算结果的对比发现结果偏小,其原因是Cowper没有考虑与外力垂直的剪应力的影响,因此新的计算结果更优越.     

一般杆系结构的非线性数值分析

本文在total-Lagrange坐标系下,对Kirchhoff梁给出了考虑几何非线性的两种梁单元刚度的显式表达式.一种是一般的非线性梁元,它既考虑了应变增量和位移增量间的二次项,又计及了刚体位移的影响,另一种是简化的非线性梁元,它只在线性梁的平衡方程中直接加入了轴力对弯曲的影响.非线性方程采用混合法求解,文中通过一些算例的数值计算,对两种单元作了比较详细的分析和评估.     

一种新的计算Timoshenko梁截面剪切系数的方法

Timoshenko梁理论中考虑了截面剪切变形的影响,推导了一种新的计算剪切系数的方法.首先采用悬臂梁纯弯曲变形条件下截面剪应力分布的精确解,基于能量原理得到了各种梁截面剪切系数新的表达式,然后推导了弯扭耦合变形条件下截面剪应力分布的精确解,进一步获得了该条件下截面的剪切系数.结果表明,悬臂梁端面作用力偏离截面的弯曲中心将使剪切系数变小,通过与Cowper计算结果的对比发现结果偏小,其原因是Cowper没有考虑与外力垂直的剪应力的影响,因此新的计算结果更优越.     

考虑纤维弯曲刚度的橡胶-帘线复合材料各向异性超弹性本构模型

针对纤维增强复合材料的有限变形,基于Spencer的连续介质力学不变量理论,提出了一种考虑纤维弯曲刚度的非线性超弹性本构模型.通过引入变形后纤维方向向量的梯度项,把单位体积的自由应变能分解为便于参数识别的体积变形、等容变形、各向异性变形和弯曲刚度4部分.理论和实验分析均表明传统的基于连续介质力学的纤维增强复合材料有限变形理论不适用于弯曲变形,必须考虑纤维弯曲刚度的影响.数值仿真结果也验证了在应变能函数中增加弯曲刚度项是必要的.     

基于纵横弯曲有限元法的套管扶正器安放计算

套管扶正器的安放问题一直以来都是石油作业的主要问题之一.基于纵横弯曲理论,考虑轴向载荷对梁弯曲变形的影响,将套管离散成BEAM188梁单元,进行几何非线性有限元分析.通过水平段、斜直段、曲率段算例与解析解对比,验证了有限元法的准确性.对某一水平井进行了有限元计算,结果表明:共安放60个扶正器,套管弯曲变形小于许可偏心距,满足工程要求.     

单广义位移的深梁理论和中厚板理论

经典的梁板弯曲理论由于未考虑横向剪切变形的影响而只能适用于细长梁和薄板,传统的多广义位移的深梁理论和中厚板理论由于忽视了转角与挠度之间的内在关系而只能适用于短粗梁和中厚板。这两种理论存在着转角的独立性与不独立性之间的矛盾,因而不相兼容。鉴于此从基本假设出发,既考虑了横向剪切影响,又确定了转角与挠度的关系,导出了单广义位移的深梁理论和中厚板理论,给出了几种简单梁的解析解,并用数值算例验证了这一理论的适用性。     

Nonlinear finite element analysis of functionally graded plates integrated with patches of piezoelectric fiber reinforced composite

In this paper, a nonlinear static finite element analysis of simply supported smart functionally graded (FG) plates in the presence/absence of the thermal environment has been presented. The substrate FG plate is integrated with the patches of piezoelectric fiber reinforced composite (PFRC) material which act as the distributed actuators of the plate. The material properties of the FG substrate plate are assumed to be temperature dependent and graded along the thickness direction according to a simple power-law distribution in terms of the volume fractions of the constituents. The derivation of this nonlinear thermo-electro-mechanical coupled finite element model is based on the first order shear deformation theory and the Von Karman type geometric nonlinearity. The numerical solutions of the nonlinear equations of the finite element model are obtained by employing the direct iteration method. The numerical illustrations suggest the potential use of the distributed actuator made of the PFRC material for active control of nonlinear deformations of smart FG structures. The effects of volume fraction index of the FG material of the substrate plates and the locations of the PFRC patches on the control authority of the patches are investigated. Emphasis has also been placed on investigating the effect of variation of piezoelectric fiber orientation angle in the PFRC patches on their actuation capability for counteracting the large deflections of FG plates.     

Large-amplitude free vibrations of functionally graded beams by means of a finite element formulation

The large-amplitude free vibration analysis of functionally graded beams is investigated by means of a finite element formulation. The Von-Karman type nonlinear strain–displacement relationships are employed where the ends of the beam are constrained to move axially. The effects of the transverse shear deformation and rotary inertia are included based upon the Timoshenko beam theory. The material properties are assumed to be graded in the thickness direction according to the power-law distribution. A statically exact beam element which devoid the shear locking effect with displacement fields based on the first order shear deformation theory is used to study the geometric nonlinear effects on the vibrational characteristics of functionally graded beams. The finite element method is employed to discretize the nonlinear governing equations, which are then solved by the direct numerical integration technique in order to obtain the nonlinear vibration frequencies of functionally graded beams with different boundary conditions. The influences of power-law exponent, vibration amplitude, beam geometrical parameters and end supports on the free vibration frequencies are studied. The present numerical results compare very well with the results available from the literature where possible. Some new results for the nonlinear natural frequencies are presented in both tabular and graphical forms which can be used for future references.     

复合材料层合剪切圆柱曲板在侧压作用下的后屈曲

基于Reddy高阶剪切变形理论的Kármám-Donnell型非线性壳体方程,给出复合材料层合剪切圆柱曲板在侧压作用下的后屈曲分析。将壳体屈曲的边界层理论推广到复合材料层合剪切圆柱曲板受侧压作用的情况。相应的奇异摄动法,用于确定圆柱曲板的屈曲荷载和后屈曲平衡路径。分析中同时考虑非线性前屈曲变形和初始几何缺陷的影响。数值算例给出完善和非完善,中等厚度正交铺设层合圆柱曲板的后屈曲荷载-挠度曲线。讨论了横向剪切变形,曲板几何参数,铺层数,铺展方式和初始几何缺陷等各种参数变化的影响。     

A geometrically nonlinear (3,2) reduced degree-of-freedom unified zigzag laminated beam element for large deformation analysis

A geometrically nonlinear (3,2) unified zigzag beam element is developed with a reduced number of degree-of-freedom for the large deformation analysis. The main merit of the beam element model is the Kirchhoff and Cauchy shear stress solution for large deformation and large strain analysis is more accurate. The geometrically nonlinearity is considered in the calculation of the zigzag coefficients. Thus, the results of shear Cauchy stress are matching well with solid element analysis in case of the beam with aspect ratio greater than 20 under large deformation. The zigzag coefficients are derived explicitly. The Green strain and the second Piola Kirchhoff stress are used. The second Piola Kirchhoff shear stress is continuous at the interface between adjacent layers priori. The bottom surface second Piola Kirchhoff shear stress condition is used to determine the zigzag coefficient and the top surface second Piola Kirchhoff shear stress condition is used to reduce one degree-of-freedom. The nonlinear finite element equations are derived. In the numerical tests, several benchmark problems with large deformation are solved to verify the accuracy. It is observed that the proposed beam has accurate solution for beam with aspect ratio greater than 20. The second Piola Kirchhoff and Cauchy shear stress accuracy is also good. A convergence study is also presented.     

Construction of beam elements considering von Kármán nonlinear strains using B-spline wavelet on the interval

The current paper proposes the formulation of beam elements using B-spline wavelet on the interval based wavelet finite element method by incorporating von Kármán nonlinear strains. Formulation is proposed for both Euler–Bernoulli beam theory and Timoshenko beam theory. A background cell based Gauss quadrature is proposed for numerical integration. Numerical examples are solved for transverse deflections and stresses in axial direction, and are compared with the existing converged results from finite element method. The issues of membrane and shear locking for the proposed elements are examined and solution techniques are suggested to overcome the issues.     

Dynamic analysis of functionally graded nanocomposite beams reinforced by randomly oriented carbon nanotube under the action of moving load

This work deals with a study of the vibrational properties of functionally graded nanocomposite beams reinforced by randomly oriented straight single-walled carbon nanotubes (SWCNTs) under the actions of moving load. Timoshenko and Euler-Bernoulli beam theories are used to evaluate dynamic characteristics of the beam. The Eshelby-Mori-Tanaka approach based on an equivalent fiber is used to investigate the material properties of the beam. An embedded carbon nanotube in a polymer matrix and its surrounding inter-phase is replaced with an equivalent fiber for predicting the mechanical properties of the carbon nanotube/polymer composite. The primary contribution of the present work deals with the global elastic properties of nano-structured composite beams. The system of equations of motion is derived by using Hamilton’s principle under the assumptions of the Timoshenko beam theory. The finite element method is employed to discretize the model and obtain a numerical approximation of the motion equation. In order to evaluate time response of the system, Newmark method is also used. Numerical results are presented in both tabular and graphical forms to figure out the effects of various material distributions, carbon nanotube orientations, velocity of the moving load, shear deformation, slenderness ratios and boundary conditions on the dynamic characteristics of the beam. The results show that the above mentioned effects play very important role on the dynamic behavior of the beam and it is believed that new results are presented for dynamics of FG nano-structure beams under moving loads which are of interest to the scientific and engineering community in the area of FGM nano-structures.     

GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS GLOBAL FINITE ELEMENT NONLINEAR GALERKIN METHOD FOR THE PENALIZED NAVIER-STOKES EQUATIONS

1. IntroductionIn the numerical simulation of the Navier-Stokes equations one encounters three seriousdifficulties in the case of large Reynolds numbers f the treatment of the incomPressibility con-dition divu = 0, the treatment of the noIilinear terms and the large time integration. For thetreatment of the incoInPressibility condition, one use the penalty method in the case of finiteelemellts [1--2l and for the treatmen of the noulinar terms and the large tfor integration, oneuse the nonlin…     

Nonlinear flexural response of laminated composite plates under hygro-thermo-mechanical loading

The paper deals with Chebyshev series based analytical solution for the nonlinear flexural response of the elastically supported moderately thick laminated composite rectangular plates subjected to hygro-thermo-mechanical loading. The mathematical formulation is based on higher order shear deformation theory (HSDT) and von-Karman nonlinear kinematics. The elastic foundation is modeled as shear deformable with cubic nonlinearity. The elastic and hygrothermal properties of the fiber reinforced composite material are considered to be dependent on temperature and moisture concentration and have been evaluated utilizing micromechanics model. The quadratic extrapolation technique is used for linearization and fast converging finite double Chebyshev series is used for spatial discretization of the governing nonlinear equations of equilibrium. The effects of Winkler and Pasternak foundation parameters, temperature and moisture concentration on nonlinear flexural response of the laminated composite rectangular plate with different lamination scheme and boundary conditions are presented.     

Linear and Nonlinear Finite Element Analysis of Active Composite Laminates

The paper presents a finite element concept for analysis of thin-walled active structures featuring fiber reinforced composite laminate as a passive structural material. The structure is rendered active by embedding piezoelectric material as a multifunctional material. A 9-node degenerated shell element based on the first order shear deformation theory is developed as a modelling tool capable of predicting the general behavior of the structure for controlling purposes. The von-Kármán type nonlinearities are considered. The solution strategy of the geometrical nonlinear analysis is based on the incremental approach using the updated Lagrangian formulation. Some numerical results are given to demonstrate the behavior of the element. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)