薄板小波有限元理论及其应用

利用样条小波尺度函数构造了常用的三角形和矩形薄板单元的位移函数,得到了利用小波函数表示的形函数。采用合理的局部坐标,对单元进行压缩,使单元在局部坐标区间上有其值,成功地推导出了分域的三角形和矩形薄板小波有限元列式。在此基础上,提出了弹性地基薄板的小波有限元求解方法。通过两个算例对薄板的挠度和弯矩进行了计算,数值结果表明,求解结果具有收敛快、精度高的特点。     

Effect of the order of the finite element and selection of the constrained modes in deformable body dynamics

Finite elements with different orders can be used in the analysis of constrained deformable bodies that undergo large rigid body displacements. The constrained mode shapes resulting from the use of finite elements with different orders differ in the way the stiffness of the body bending and extension are defined. The constrained modes also depend on the selection of the boundary conditions. Using the same type of finite element, different sets of boundary conditions lead to different sets of constrained modes. In this investigation, the effect of the order of the element as well as the selection of the constrained mode shapes is examined numerically. To this end, the constant strain three node triangular element and the quadratic six node triangular element are used. The results obtained using the three node triangular element are compared with the higher order six node triangular element. The equations of motion for the three and six node triangular elements are formulated from assumed linear and quadratic displacement fields, respectively. Both assumed displacement fields can describe large rigid body translational and rotational displacements. Consequently, the dynamic formulation presented in this investigation can also be used in the large deformation analysis. Using the finite element displacement field, the mass, stiffness, and inertia invariants of the three and six-node triangular elements are formulated. Standard finite element assembly techniques are used to formulate the differential equations of motion for mechanical systems consisting of interconnected deformable bodies. Using a multibody four bar mechanism, numerical results of the different elements and their respective performance are presented. These results indicate that the three node triangular element does not perform well in bending modes of deformation.     

On the formulation of the planar ANCF triangular finite elements

In this paper, new planar isoparametric triangular finite elements (FE) based on the absolute nodal coordinate formulation (ANCF) are developed. The proposed ANCF elements have six coordinates per node: two position coordinates that define the absolute position vector of the node and four gradient coordinates that define vectors tangent to coordinate lines (parameters) at the same node. To shed light on the importance of the element geometry and to facilitate the development of some of the new elements presented in this paper, two different parametric definitions of the gradient vectors are used. The first parametrization, called area parameterization, is based on coordinate lines along the sides of the element in the reference configuration, while the second parameterization, called Cartesian parameterization, employs coordinate lines defined along the axes of the structure (body) coordinate system. The fundamental differences between the ANCF parameterizations used in this investigation and the parametrizations used for conventional finite elements are highlighted. The Cartesian parameterization serves as a unique standard for the triangular FE assembly. To this end, a transformation matrix that defines the relationship between the area and the Cartesian parameterizations is introduced for each element in order to allow for the use of standard FE assembly procedure and define the structure (body) inertia and elastic forces. Using Bezier geometry and a linear mapping, cubic displacement fields of the new ANCF triangular elements are systematically developed. Specifically, two new ANCF triangular finite elements are developed in this investigation, namely four-node mixed-coordinate and three-node ANCF triangles. The performance of the proposed new ANCF elements is evaluated by comparison with the conventional linear and quadratic triangular elements as well as previously developed ANCF rectangular and triangular elements. The results obtained in this investigation show that in the case of small and large deformations as well as finite rotations, all the elements considered can produce correct results, which are in a good agreement if appropriate mesh sizes are used.     

NONLINEAR DYNAMIC ANALYSIS OF A LAMINATED HYBRID COMPOSITE PLATE SUBJECTED TO TIME-DEPENDENT EXTERNAL PULSES

Nonlinear dynamic responses of a laminated hybrid composite plate subjected to time-dependent pulses are investigated. Dynamic equations of the plate are derived by the use of the virtual work principle. The geometric nonlinearity effects are taken into account with the von Kármán large deflection theory of thin plates. Approximate solutions for a clamped plate are assumed for the space domain. The single term approximation functions are selected by considering the nonlinear static deformation of plate obtained using the finite element method. The Galerkin Method is used to obtain the nonlinear differential equations in the time domain and a MATLAB software code is written to solve nonlinear coupled equations by using the Newmark Method. The results of approximate-numerical analysis are obtained and compared with the finite element results. Transient loading conditions considered include blast, sine, rectangular, and triangular pulses. A parametric study is conducted considering the effects of peak pressure, aspect ratio, fiber orientation and thicknesses.     

Bond stresses for a composite material reinforced with various arrays of fibers

The problem treated here is that of an isotropic body having a doubly periodic rectangular or triangular array of perfectly bonded circular elastic inclusions. The body is in tension or compression. This simulates a composite material wherein a relatively weak matrix is reinforced by stronger (and more rigid) fibers. Bond stresses for both rectangular and triangular arrays have been calculated using either boundary point matching or boundary point least squares techniques. Numerical results based on a plane strain analysis are given in graphical form.     

Maximum frequency design of laminated plates with mixed boundary conditions

A layerwise optimization (LO) approach is extended to accommodate the finite element analysis for optimizing the free vibration behavior of laminated composite plates with discontinuities along the boundaries. The classical non-conforming finite element is modified to fit into the LO procedure and is used to calculate natural frequencies of symmetrically laminated plates. This combined LO-FEM approach is applied to laminated rectangular plates with combinations of free, simply supported and clamped edges along parts of the plate boundary. The fundamental frequency, an object function for the present study, is maximized by optimizing design variables that are a set of fiber orientation angles in the layers. For illustrative purpose nine examples of square and rectangular plates with various types of mixed boundary conditions are considered, and a comprehensive set of results are presented for the optimum fiber orientation angles and the maximum fundamental frequencies of the 8-layer and 24-layer plates.     

One-dimensional finite element heat transfer solution of a fin with triangular perforations of bases parallel and towered its base

Heat transfer dissipation from a horizontal rectangular fin embedded with equilateral triangular perforations is computed numerically using one-dimensional finite element technique. The bases of the triangles are parallel and toward the fin base. The body of the fin is discretized into a number of subdivisions (finite elements). The number of these elements can be altered as required according to the automatic mesh generation. The heat dissipation of the perforated fin is computed and compared with that of the solid one of the same dimensions and same thermal properties. The comparison refers to acceptable results and heat dissipation enhancement due to certain perforation.     

Numerical solution of the Navier-Stokes equations using a finite element method

A finite element algorithm for solving the Navier-Stokes equations is presented for the analysis of high-speed viscous flows. The algorithm uses triangular elements. The unsteady equations are integrated to steady state with a Runge-Kutta time-marching scheme. A postprocessing artificial dissipation term is introduced to stabilize the computations and to dampen dissipation errors. Numerical results are compared with the calculation of uniform flow on a rectangular region which encounters an embedded oblique shock. A shock/turbulent boundary layer problem is also solved and results are compared with experimental data. It is shown that the postprocessing smoothing term and boundary conditions similar to the finite difference method work well in the present numerical studies.     

Flow/damage surfaces for fiber-reinforced metals having different periodic microstructures

Flow/damage surfaces can be defined in terms of stress, inelastic strain rate, and internal variables using a thermodynamics framework. A macroscale definition relevant to thermodynamics and usable in an experimental program is employed to map out surfaces of constant inelastic power in various stress planes. The inelastic flow of a model silicon carbide/titanium composite system having rectangular, hexagonal, and square diagonal fiber packing arrays subjected to biaxial stresses is quantified by flow/damage surfaces that are determined numerically from micromechanics, using both finite element analysis and the generalized method of cells. Residual stresses from processing are explicitly included and damage in the form of fiber–matrix debonding under transverse tensile and/or shear loading is represented by a simple interface model. The influence of microstructural architecture is largest whenever fiber-matrix debonding is not an issue; for example in the presence of transverse compressive stresses. Additionally, as the fiber volume fraction increases, so does the effect of microstructural architecture. These results indicate that microstructural architecture needs to be accounted for in an accurate continuum model, thus complicating development of such a model. With regard to the micromechanics analysis, the overall inelastic flow predicted by the generalized method of cells is in excellent agreement with that predicted using a large number of displacement-based finite elements.     

基于高阶位移模式的压电层合板动力有限元模型

建立了含压电片层合板的有限元动力学模型。以位于压电层上下表面处的电场强度和层间电压为未知量,给出了三次函数的电势分布模式,采用Reddy的高阶剪切理论描述板的位移场,假设板厚度方向的正应力为零给出了减缩的本构方程,采用有限元方法,基于Hamilton原理导出结构的动力学方程,然后用静态缩聚的方法压缩掉电场自由度和次要的位移自由度。最后用四边形矩形单元求解了一对称铺层和非对称铺层悬臂板的固有频率,并与ANSYS结果对比验证了本文模型的精确性。     

基于高阶剪切变形理论的四边形求积元板单元及其应用

近年来由各类新型复合材料或功能梯度材料构成的板结构在工程领域得到了广泛应用,其显著特点是材料性能沿板厚变化.为合理考虑横向剪切应变,许多学者基于Reddy高阶剪切变形理论,构建了不同的有限元单元对该类板结构进行分析,但其中满足$C^{1}$连续条件的单元相对较少.本文基于Reddy高阶剪切变形理论,采用求积元方法,建立了$C^{1}$连续的四边形板单元.利用该单元对均质材料、复合材料、功能梯度材料构成的等厚度矩形板、变厚度矩形板及等厚度斜板的线弹性弯曲和自由振动问题进行了计算分析,并与现有文献中的相应计算结果进行了对比.研究表明:基于高阶剪切变形理论的四边形求积元板单元具有较高的计算效率和良好的适应性,文中各类材料构成的等变厚度矩形板及等厚度斜板均只需1个单元即可得到理想的计算结果.对于等/变厚度矩形板,可仅使用9$\times$9个积分点,而对于等厚度斜板,随着斜角的增大,所需积分点的数目逐渐增多至15$\times $15.该四边形求积元板单元可进一步用于新型复合材料板的非线性分析.      

Measurement of inhomogeneous strain fields by fiber optic sensors embedded in a polymer composite material

Experimental results of strain field measurement in polymer composite specimens by Bragg grating fiber optic strain sensors embedded in the material are considered. A rectangular plate and a rectangular plate with “butterfly” shaped cuts are used as specimens. The results of uniaxial strain experiments with rectangular plates show that fiber optic strain sensors can be used to measure the strains, and these results can be used to calculate the calibration coefficients for fiber optic strain sensors. A gradient strain field is attained in a plate with cuts, and the possibility of measuring this field by fiber optic strain sensors is the main goal of this paper. The results of measurements of gradient strain fields in the plate with cuts are compared with the results obtained by using the three-dimensional digital optic system Vix-3D and with the results of numerical computations based on finite element methods. It is shown that the difference between the strain values obtained by these three methods does not exceed 5%.     

多边形有限单元形函数的比较研究

多边形有限单元形函数有wachspress插值、Laplace插值和平均值插值三种类型.本文对三种多边形有限单元形函数的性质作了比较研究,给出了三种形函数各自的优点和局限性.Waclaspress和Laplace形函数是有理函数形式,而平均值形函数是无理函数形式.三种形函数均满足单位分解性、线性完备性,且在单元边界上呈线性.在三角形单元上,它们都等价于三角形面积坐标插值.在矩形单元上,Wachspress和Laplace形函数等价于双线性多项式插值形函数.Wachspress和平均值形函数适用于任意凸多边形单元,Laplace形函数更适用于圆内接多边形单元.Wachspress形函数不能推广到含有边节点的单元,平均值形函数可以直接推广到含有边节点的单元.数值试验,验证了本文理论分析的结论.     

Finite element implementation and numerical issues of strain gradient plasticity with application to metal matrix composites

A framework of finite element equations for strain gradient plasticity is presented. The theoretical framework requires plastic strain degrees of freedom in addition to displacements and a plane strain version is implemented into a commercial finite element code. A couple of different elements of quadrilateral type are examined and a few numerical issues are addressed related to these elements as well as to strain gradient plasticity theories in general. Numerical results are presented for an idealized cell model of a metal matrix composite under shear loading. It is shown that strengthening due to fiber size is captured but strengthening due to fiber shape is not. A few modelling aspects of this problem are discussed as well. An analytic solution is also presented which illustrates similarities to other theories.     

ANALYSIS OF COMPOSITE LAMINATE BEAMS USING COUPLING CROSS-SECTION FINITE ELEMENT METHOD

Beams and plates manufactured from laminates of composite materials have distinct advantages in a significant number of applications. However, the anisotropy arising from these materials adds a significant degree of complexity, and thus time, to the stress and deformation analyses of such components, even using numerical approaches such as finite elements. The analysis of composite laminate beams subjected to uniform extension, bending, and/or twisting loads was performed by a novel implementation of the usual finite element method. Due to the symmetric features of the deformations, only a thin slice of the beam to be analysed needs to be modelled. Conventional three-dimensional solid finite elements were used for the structural discretization. The accurate deformation relationships were formulated and implemented through the coupling of nodal translational degrees of freedom in the numerical analysis. A sample solution for a rectangular composite laminate beam is presented to show the validity and accuracy of the proposed method.     

纤维增强复合材料弹性性能预测的域分解方法及应用

提出了新的有限元建模方法,即域分解方法,用于预测纤维增强复合材料单向带T300/BSL914C(环氧树脂)和AS4/3501-6(环氧树脂)的弹性性能。域分解方法基于区域叠合技术,分别建立单胞的整体域与纤维域模型用于代替传统有限元建模方法中单胞的基体域与纤维域模型。整体域是真实基体体积与纤维体积的叠加,两区域网格独立划分,互不影响。采用MSC.Nastran中的多节点约束Explicit单元,在整体域与纤维域节点之间建立位移连接属性模拟单胞基体域与纤维域之间的位移约束关系,从而实现两区域的耦合计算。计算结果表明:域分解方法单胞模型纤维增强方向弹性模量Ez预测值与试验值误差在7%以内,其余弹性常数也都与试验值吻合较好。域分解方法不仅可以大大简化纤维增强复合材料的细观力学建模,而且可以准确地预测纤维增强复合材料的弹性性能。     

Measurement and modeling of stress concentration around a circular notch

When designing composite materials, the presence of stress concentration at locations such as circular nothces is unavoidable. Such locations in structural elements arise from joints required to form a structure. The stress concentration, observed around the notch, is quantified by the stress concentration factorK. This quantity is normally calculated analytically and/or numerically and is an important design parameter. In this work, the experimental technique of remote laser Raman microscopy is used for the in situ measurement ofK in Kevlar 49 fiber/epoxy composite plates containing a circular hole. The results obtained by this technique are compared with those calculated analytically and by finite element analysis. Both analytical and numerical methods underestimate the experimental results for maximumK by approximately 10 percent, which is considered reasonable within experimental error. In addition, very good agreement between analytical and experimental data is obtained for the decay of the stress concentration factor as a function of distance from the edge of the hole. The numerical results, however, overestimate the decayK with distance from the notch boundary and only converge at relatively large distances.B. P. Arjyal has obtained his PhD from the Materials Department, Queen Mary & Westfield College, University of London, London, E1 4NS, United Kingdom.     

含内埋矩形脱层复合材料圆柱壳屈曲分析的混合变量条形传递函数方法

提出了一种分析含内埋矩形脱层正交各向异性圆柱壳稳定性问题的混合变量条形传递函数方法。首先基于Mindlin一阶剪切壳理论,通过定义圆柱壳的广义力变量和混合变量,建立了壳的改进混合变量能量泛函;然后,为了便于脱层壳的分区求解,通过引入条形单元,创建了基于混合变量条形传递函数解的含脱层和不合脱层两种超级壳单元;在此基础上,将含内埋矩形脱层的复合材料层合壳划分成两种超级壳单元的组合体,通过各超级壳单元相互之间连接结点处的位移连续和力平衡条件得到脱层壳的屈曲方程;最后由屈曲方程计算含内埋矩形脱层壳的屈曲载荷和屈曲模态。算例分析的结果验证了本方法的正确性,并给出了几种因素对屈曲载荷和屈曲模态的影响。     

The criterion for determining the direction of crack propagation in a random pattern composites

Presented paper contains results of fracture analysis of brittle composite materials with a random distribution of grains. The composite structure has been modelled as an isotropic matrix that surrounds circular grains with random diameters and space position. Analyses were preformed for the rectangular “numerical sample” by finite element method. FE mesh for the examples were generated using the authors’ computer program RandomGrain. Fracture analyses were accomplished with the authors’ computer program CrackPath3 executing the “fine mesh window” technique. Calculations were preformed in 2D space assuming the plane stress state. Current efforts focus on brittle materials such as rocks or concrete.     

用拟协调元直接构造平面任意四边形单元——进入有限元的禁区

利用拟协调元方法,在直角坐标系下直接构造了一族平面任意四边形单元,对其收敛性进行了分析,并与平面等参元进行了对比研究.结果证明平面任意四边形单元可采用多项式基函数直接列式,并可以保障单元的收敛性;拟协调元列式可以使平面问题的有限元方法得到统一.与平面等参元相比,单元列式简单,性能稳定,具有显式的刚度阵,计算量小,这说明对于有限元平面问题拟协调元是一个更正确、有效的做法.