首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 0 毫秒
1.
To analyze angle-ply laminated composite and sandwich plates coupled bending and extension under thermo-mechanical loading, a refined global–local higher-order theory considering transverse normal strain is presented in this work. Hitherto, present theory for angle-ply laminates has never been reported in the literature, and this theory can satisfy continuity of transverse shear stresses at interfaces. In addition, the number of unknowns in present model is independent of layer numbers of the laminate. Based on this theory as well as methodology of the refined triangular discrete Kirchhoff plate element, a triangular laminated plate element satisfying the requirement of C1 continuity is presented. Numerical results show that the present refined theory can accurately analyze the bending problems of angle-ply composite and sandwich plates as well as thermal expansion problem of cross-ply plates, and the present refined theory is obviously superior to the existing global–local higher-order theory proposed by Li and Liu [Li XY, Liu D. Generalized laminate theories based on double superposition hypothesis. Int J Numer Meth Eng 1997;40:1197–212]. After ascertaining the accuracy of present model, the distributions of displacements and stresses for angle-ply laminated plates under temperature loads are also given in present work. These results can serve as a reference for future investigations.  相似文献   

2.
    
An efficient triangular element based on an inter‐laminar shear stress continuous plate theory is developed and applied to the analysis of composite and sandwich plates under different situations to study the performance of the element. The plate theory represents parabolic through thickness variation of transverse shear stresses where the continuity condition of these stresses are satisfied at the layer interfaces. It also satisfies transverse shear stress free condition at the top and bottom surfaces of the plate. The most attractive feature of the plate theory is that the basic unknowns are same as those used in first‐order shear deformation theory. The only problem lies with this elegant plate theory is found in its finite element implementation, as it requires C1 continuity of transverse displacement at the element interfaces. This is a well‐known problem of thin plate elements, which is also found in some other refined plate theories. Although there are some elements based on these refined plate theories but the number of such elements is very few and they possess certain drawbacks in general. Keeping these aspects in view, an attempt has been made in this study to develop a six‐noded triangular element having equal degrees of freedom at each node. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
基于精化锯齿理论,构造了六节点三角形协调板单元并推导了夹层板自由振动问题有限元列式。不同于已有锯齿理论,精化锯齿理论特点是面内位移不含有横向位移一阶导数,构造有限元时仅需要C0 插值函数。为验证单元性能,分析了软核夹层板自由振动问题。结果表明,该文构造的单元能准确计算软核夹层板固有频率,然而基于已有锯齿理论建立的不协调元计算结果精度较低。  相似文献   

4.
    
The existing global–local multiscale computational methods, using finite element discretization at both the macro‐scale and micro‐scale, are intensive both in terms of computational time and memory requirements and their parallelization using domain decomposition methods incur substantial communication overhead, limiting their application. We are interested in a class of explicit global–local multiscale methods whose architecture significantly reduces this communication overhead on massively parallel machines. However, a naïve task decomposition based on distributing individual macro‐scale integration points to a single group of processors is not optimal and leads to communication overheads and idling of processors. To overcome this problem, we have developed a novel coarse‐grained parallel algorithm in which groups of macro‐scale integration points are distributed to a layer of processors. Each processor in this layer communicates locally with a group of processors that are responsible for the micro‐scale computations. The overlapping groups of processors are shown to achieve optimal concurrency at significantly reduced communication overhead. Several example problems are presented to demonstrate the efficiency of the proposed algorithm. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
Based on the higher-order global–local theories, a finite element model is proposed to study the bending behavior of stiffened laminated plates. The proposed model treats the embedded stiffeners as the part of laminated plate, so that the compatibility of displacements and stresses between the plate and the stiffeners can be automatically satisfied. Distributions of the displacements and stresses through the thickness of laminates were also given for the first time, which can serve as references for future investigations as such information is lacking in the published literature. In addition, the impact of the stiffeners on the bending response of the stiffened laminated plates is also studied in terms of the quantity, the collocation and the geometry of stiffeners. Numerical results showed that the higher-order global–local theories are more suitable for predicting the bending response of thick and moderately thick stiffened laminated plates compared to the first order theory commonly used in engineering. By varying the quantity, the collocation and the geometry of stiffeners, the stiffness and the strength of stiffened laminated plates can be remarkably improved.  相似文献   

6.
    
The accuracy of stability analysis depends on the accuracy of both the element stiffness matrix and geometry stiffness matrix. Therefore, when carrying out the stability analysis of thin cylindrical shells using the finite element methods will require, firstly, a refined non‐conforming rectangular curved cylindrical shell element RCSR4 is proposed according to the refined non‐conforming FE method, in which both the C1 and C0 weak continuity conditions are satisfied and as a result, can ensure the convergence of computation. At the same time, a refined geometrical stiffness matrix is introduced to replace the standard consistent geometrical stiffness matrix. Simple expressions of the refined constant strain matrices with adjustable constants are introduced with respect to the weak continuity conditions. Numerical examples are presented to show that the present method can indeed improve the performance and the accuracy in stability analysis. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

7.
    
In this paper, we develop a block preconditioner for Jacobian‐free global–local multiscale methods, in which the explicit computation of the Jacobian may be circumvented at the macroscale by using a Newton–Krylov process. Effective preconditioning is necessary for the Krylov subspace iterations (e.g. GMRES) to enhance computational efficiency. This is, however, challenging since no explicit information regarding the Jacobian matrix is available. The block preconditioning technique developed in this paper circumvents this problem by effectively deflating the spectrum of the Jacobian matrix at the current Newton step using information about only the Krylov subspaces corresponding to the Jacobian matrices in the previous Newton steps and their representations on those subspaces. This approach is optimal and results in exponential convergence of the GMRES iterations within each Newton step, thus minimizing expensive microscale computations without requiring explicit Jacobian formation in any step. In terms of both computational cost and storage requirements, the action of a single block of the preconditioner per GMRES step scales linearly as the number of degrees of freedom of the macroscale problem as well as the dimension of the invariant subspace of the preconditioned Jacobian matrix. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

8.
    
A new method is proposed to place local meshes in a global mesh with the aid of the interface‐element method (IEM). The interface‐elements use moving least‐square (MLS)‐based shape functions to join partitioned finite‐element domains with non‐matching interfaces. The supports of nodes are defined to satisfy the continuity condition on the interfaces by introducing pseudonodes on the boundaries of interface regions. Particularly, the weight functions of nodes on the boundaries of interface regions span only neighbouring nodes, ensuring that the resulting shape functions are identical to those of adjoining finite‐elements. The completeness of the shape functions of the interface‐elements up to the order of basis provides a reasonable transfer of strain fields through the non‐matching interfaces between partitioned domains. Taking these great advantages of the IEM, local meshes can be easily inserted at arbitrary places in a global mesh. Several numerical examples show the effectiveness of this technique for modelling of local regions in a global domain. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

9.
A refined non-conforming quadrilateral thin plate bending element RPQ4 which can satisfy the requirement of convergence is established such that the non-conforming displacement function can be derived directly. A simple explicit expression of a refined constant strain matrix can be introduced into the formulation of the standard displacement element which results in the constraint condition of interelement continuity being satisfied in an average sense. Numerical examples are presented to show that the present model can pass the patch test and possesses high accuracy. © 1997 John Wiley & Sons, Ltd.  相似文献   

10.
    
Based on the refined non‐conforming element method, simple flat triangular elements with standard nodal displacement parameters are proposed for the analysis of shell structures. For ensuring the convergence of the elements a new coupled continuity condition at the inter‐element has been established in a weaker form. A common displacement for the inter‐element, an explicit expression of refined constant strain matrix, and an adjustable constant are introduced into the formulation, in which the coupled continuity requirement at the inter‐element is satisfied in the average sense. The non‐conforming displacement function of the well‐known triangular plate element BCIZ [1] and the membrane displacement of the constant strain triangular element CST [2] are employed to derive the refined flat shell elements RTS15, and the refined flat shell elements RTS18 is derived by using the element BCIZ and the Allman's triangular plane element [3] with the drilling degrees of freedom. A simple reduced higher‐order membrane strain matrix is proposed to avoid membrane locking of the element RTS18. An alternative new reduced higher‐order strain matrix method is developed to improve the accuracy of the elements RTS15 and RTS18. Numerical examples are given to show that the present methods have improved the accuracy of the shell analysis. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

11.
复合材料层合厚圆柱壳高阶理论的改进及其应用   总被引:3,自引:1,他引:2  
建立了一个改进的LCW型的精化高阶理论,以分析厚圆柱壳的振动。提出u,v为三次多项式、w为二次多项式的位移模式,并利用上、下自由表面横向剪应力为零的边界条件,对所假定的位移场作了化简,将三阶剪切变形理论的未知数缩减为7个,在此基础上建立了相应的有限元列式。通过一个典型算例,与Soldatos和Lam的高阶剪切变形理论的解析解作了比较,说明笔者的精化高阶理论是可行的,而且具有较高的精确性,比LCW高阶理论更具有实用性。还通过频率参数随长度半径比L/R的变化,说明由于考虑了法向应力和法向应变,本文方法更适用于长度半径比较小的结构。  相似文献   

12.
    
In his article a special form of gradient elasticity is presented that can be used to describe wave dispersion. This new format of gradient elasticity is an appropriate dynamic extension of the earlier static counterpart of the gradient elasticity theory advocated in the early 1990s by Aifantis and co‐workers. In order to capture dispersion of propagating waves, both higher‐order inertia and higher‐order stiffness contributions are included, a fact which implies (and is denoted as) dynamic consistency. The two higher‐order terms are accompanied by two associated length scales. To facilitate finite element implementations, the model is rewritten such that ??0‐continuity of the interpolation is sufficient. An auxiliary displacement field is introduced which allows the original fourth‐order equations to be split into two coupled sets of second‐order equations. Positive‐definiteness of the kinetic energy requires that the inertia length scale is larger than the stiffness length scale. The governing equations, boundary conditions and the discretized system of equations are presented. Finally, dispersive wave propagation in a one‐dimensional bar is considered in a numerical example. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

13.
    
This paper presents new developments on a weakly intrusive approach for the simplified implementation of space and time multiscale methods within an explicit dynamics software. The ‘substitution’ method proposed in previous works allows to take advantage of a global coarse model, typically used in an industrial context, running separate, refined in space and in time, local analyses only where needed. The proposed technique is iterative, but the explicit character of the method allows to perform the global computation only once per global time step, while a repeated solution is required for the small local problems only. Nevertheless, a desirable goal is to reach convergence with a reduced number of iterations. To this purpose, we propose here a new iterative algorithm based on an improved interface inertia operator. The new operator exploits a combined property of velocity Hermite time interpolation on the interface and of the central difference integration scheme, allowing the consistent upscaling of interface inertia contributions from the lower scale. This property is exploited to construct an improved mass matrix operator for the interface coupling, allowing to significantly enhance the convergence rate. The efficiency and robustness of the procedure are demonstrated through several examples of growing complexity. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

14.
C1连续曲面重构与光顺的有限元算法   总被引:3,自引:0,他引:3  
提出了一种基于离散的测量数据重建光顺自由曲面的有限元新方法。根据最佳逼近与能量光顺原理,建立正定的目标泛函,采用18自由度三角形板单元对泛函离散,进行极小化,求得最优解。根据有限元插值计算,重新构造出全场C^1连续的自由曲面。这种方法结合了能量光顺技术,有效地抑制了输入数据上误差噪声的影响,曲面重建的精度高、光顺性好,而且能给出合理的一阶导数。该方法计算简单、便于应用,所需的输入数据点少,并可用于处理曲线边界区域的问题。  相似文献   

15.
A displacement-based variable kinematic global–local finite element model is developed using hierarchical, multiple assumed displacement fields at two different levels: (1) at the element level, and (2) at the mesh level. The displacement field hierarchy contains both a conventional plate expansion (2-D) and a full layerwise (3-D) expansion. Depending on the accuracy desired, the variable kinematic element can use various terms from the composite displacement field, thus creating a hierarchy of different elements having a wide range of kinematic complexity and representing a number of different mathematical models. The VKFE is then combined with the mesh superposition technique to further increase the computational efficiency and robustness of the computational algorithm. These models are used to analyse a number of laminated composite plate problems that contain localized subregions where significant 3-D stress fields exist (e.g. free-edge effects).  相似文献   

16.
    
Three‐dimensional higher‐order eXtended finite element method (XFEM)‐computations still pose challenging computational geometry problems especially for moving interfaces. This paper provides a method for the localization of a higher‐order interface finite element (FE) mesh in an underlying three‐dimensional higher‐order FE mesh. Additionally, it demonstrates, how a subtetrahedralization of an intersected element can be obtained, which preserves the possibly curved interface and allows therefore exact numerical integration. The proposed interface algorithm collects initially a set of possibly intersecting elements by comparing their ‘eXtended axis‐aligned bounding boxes’. The intersection method is applied to a highly reduced number of intersection candidates. The resulting linearized interface is used as input for an elementwise constrained Delaunay tetrahedralization, which computes an appropriate subdivision for each intersected element. The curved interface is recovered from the linearized interface in the last step. The output comprises triangular integration cells representing the interface and tetrahedral integration cells for each intersected element. Application of the interface algorithm currently concentrates on fluid–structure interaction problems on low‐order and higher‐order FE meshes, which may be composed of any arbitrary element types such as hexahedra, tetrahedra, wedges, etc. Nevertheless, other XFEM‐problems with explicitly given interfaces or discontinuities may be tackled in addition. Multiple structures and interfaces per intersected element can be handled without any additional difficulties. Several parallelization strategies exist depending on the desired domain decomposition approach. Numerical test cases including various geometrical exceptions demonstrate the accuracy, robustness and efficiency of the interface handling. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

17.
    
A new formulation of gradient elasticity is derived and implemented. Padé approximations are used to introduce an implicit relation between non‐local strains and displacements. As a result, the finite element interpolation requires only ??0‐continuous (rather than ??1‐continuous) shape functions. The underlying energy functional is presented and it is found that the present formulation is of the mixed type, whereby the non‐local strains act as the primary unknowns and the displacements as the auxiliary unknowns. The implications on the interpolation are studied. Finally, the influence of the additional length scale parameter on the global response is assessed. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

18.
    
In gradient elasticity strain gradient terms appear in the expression of virtual work, leading to the need for C1 continuous interpolation in finite element discretizations of the displacement field only. Employing such interpolation is generally avoided in favour of the alternative methods that interpolate other quantities as well as displacement, due to the scarcity of C1 finite elements and their perceived computational cost. In this context, the lack of three‐dimensional C1 elements is of particular concern. In this paper we present a new C1 hexahedral element which, to the best of our knowledge, is the first three‐dimensional C1 element ever constructed. It is shown to pass the single element and patch tests, and to give excellent rates of convergence in benchmark boundary value problems of gradient elasticity. It is further shown that C1 elements are not necessarily more computationally expensive than alternative approaches, and it is argued that they may be more efficient in providing good‐quality solutions. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

19.
Abstract

The partial hybrid stress model is applied to the refined C 1 higher‐order plate theory in this paper. The displacement model is adopted in the flexural part and the hybrid stress model in the transverse shear part. The plate concept is introduced and the governing equations of plate are derived variationally from the modified Hellinger‐Reissner principle. This new plate element is demonstrated to be more accurate than displacement formulation in the analysis of orthotropic thick laminated plates. Moreover, the through thickness distribution of transverse shear stress is precisely predicted.  相似文献   

20.
A refined triangular discrete Kirchhoff thin plate bending element RDKT which can be used to improve the original triangular discrete Kirchhoff thin plate bending element DKT is presented. In order to improve the accuracy of the analysis a simple explicit expression of a refined constant strain matrix with an adjustable constant can be introduced into its formulation. The new element displacement function can be used to formulate a mass matrix called combined mass matrix for calculation of the natural frequency and in the same way a combined geometric stiffness matrix can be obtained for buckling analysis. Numerical examples are presented to show that the present methods indeed, can improve the accuracy of thin plate bending, vibration and buckling analysis. © 1998 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号