共查询到20条相似文献,搜索用时 0 毫秒
1.
It is well known that the finite element method (FEM) suffers severely from the volumetric locking problem for incompressible materials in topology optimization owing to its numerical ‘overly stiff’ property. In this article, two typical smoothed FEMs with a certain softened effect, namely the node-based smoothed finite element method (NS-FEM) and the cell-based smoothed finite element method, are formulated to model the compressible and incompressible materials for topology optimization. Numerical examples have demonstrated that the NS-FEM with an ‘overly soft’ property is fairly effective in tackling the volumetric locking problem in topology optimization when both compressible and incompressible materials are involved. 相似文献
2.
3.
4.
G. Castellazzi P. Krysl 《International journal for numerical methods in engineering》2012,90(13):1618-1635
A finite element model for linear‐elastic small deformation problems is presented. The formulation is based on a weighted residual that requires a priori the satisfaction of the kinematic equation. In this approach, an averaged strain‐displacement matrix is constructed for each node of the mesh by defining an appropriate patch of elements, yielding a smooth representation of strain and stress fields. Connections with traditional and similar procedure are explored. Linear quadrilateral four‐node and linear hexahedral eight‐node elements are derived. Various numerical tests show the accuracy and convergence properties of the proposed elements in comparison with extant finite elements and analytic solutions. Specific examples are also included to illustrate the ability to resist numerical locking in the incompressible limit and insensitive response in the presence of shape distortion. Furthermore, the numerical inf‐sup test is applied to a selection of problems to show the stability of the present formulation. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
5.
6.
C.T. Wu W. Hu J.S. Chen 《International journal for numerical methods in engineering》2012,90(7):882-914
In this paper, a two‐dimensional displacement‐based meshfree‐enriched FEM (ME‐FEM) is presented for the linear analysis of compressible and near‐incompressible planar elasticity. The ME‐FEM element is established by injecting a first‐order convex meshfree approximation into a low‐order finite element with an additional node. The convex meshfree approximation is constructed using the generalized meshfree approximation method and it possesses the Kronecker‐delta property on the element boundaries. The gradient matrix of ME‐FEM element satisfies the integration constraint for nodal integration and the resultant ME‐FEM formulation is shown to pass the constant stress test for the compressible media. The ME‐FEM interpolation is an element‐wise meshfree interpolation and is proven to be discrete divergence‐free in the incompressible limit. To prevent possible pressure oscillation in the near‐incompressible problems, an area‐weighted strain smoothing scheme incorporated with the divergence‐free ME‐FEM interpolation is introduced to provide the smoothing on strains and pressure. With this smoothed strain field, the discrete equations are derived based on a modified Hu–Washizu variational principle. Several numerical examples are presented to demonstrate the effectiveness of the proposed method for the compressible and near‐incompressible problems. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
7.
An assumed‐strain finite element technique is presented for linear, elastic small‐deformation models. Weighted residual method (reminiscent of the strain–displacement functional) is used to weakly enforce the balance equation with the natural boundary condition and the kinematic equation (the strain–displacement relationship). A priori satisfaction of the kinematic weighted residual serves as a condition from which strain–displacement operators are derived via nodal integration. A variety of element shapes is treated: linear triangles, quadrilaterals, tetrahedra, hexahedra, and quadratic (six‐node) triangles and (27‐node) hexahedra. The degrees of freedom are only the primitive variables (displacements at the nodes). The formulation allows for general anisotropic materials. A straightforward constraint count can partially explain the insensitivity of the resulting finite element models to locking in the incompressible limit. Furthermore, the numerical inf–sup test is applied in select problems and several variants of the proposed formulations (linear triangles, quadrilaterals, tetrahedra, hexahedra, and 27‐node hexahedra) pass the test. Examples are used to illustrate the performance with respect to sensitivity to shape distortion and the ability to resist volumetric locking. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
8.
J.‐H. Hommel G. Meschke 《International journal for numerical methods in engineering》2008,75(12):1416-1434
The deficiency of volumetric locking phenomena in finite elements using higher‐order shell element formulations based on Lagrangean polynomials and a linear finite shell kinematics cannot be avoided by the existent enhanced assumed strain (EAS) concept established for low‐order elements. In this paper a consistent modification of the EAS concept is proposed to extend its applicability to higher‐order shell elements. This modification, affecting the transversal normal strain for polynomial orders p>1, eliminates pathological modes caused by volumetric locking. The efficiency of the proposed extended EAS method is demonstrated by means of eigenvalue analyses and two representative numerical examples. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
9.
A novel finite element (FE) formulation with adaptive mesh rezoning for large deformation problems is proposed. The proposed method takes the advantage of the selective smoothed FE method (S‐FEM), which has been recently developed as a locking‐free FE formulation with strain smoothing technique. We adopt the selective face‐based smoothed/node‐based smoothed FEM (FS/NS‐FEM‐T4) and edge‐based smoothed/node‐based smoothed FEM (ES/NS‐FEM‐T3) basically but modify them partly so that our method can handle any kind of material constitutive models other than elastic models. We also present an adaptive mesh rezoning method specialized for our S‐FEM formulation with material constitutive models in total form. Because of the modification of the selective S‐FEMs and specialization of adaptive mesh rezoning, our method is locking‐free for severely large deformation problems even with the use of tetrahedral and triangular meshes. The formulation details for static implicit analysis and several examples of analysis of the proposed method are presented in this paper to demonstrate its efficiency. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
10.
鉴于双头电火花加工机床部件结构复杂,很难建立力学模型采分析,随着计算机技术的不断发展,有限元方法已经广泛应用于工程计算,而Y轴横梁是三轴双头电火花加工机床的关键部件,因此应用有限元方法对其进行静力学分析,得到Y轴横梁在不同位置时导轨滑块处所受的载重以及主轴头端面的变形量,为直线导轨的选择和安装提供了理论依据. 相似文献
11.
12.
三维动态八结点有限单元的建立及理论推导 总被引:5,自引:0,他引:5
本文为扩大动态有限单元法的应用范围,在 J.S.Przemieniecki 和 K.K.Gupta 等人的研究基础上建立了一种新的三维动态八结点单元,并对其进行了理论推导。通过算例验证了本文理论结果的正确性。 相似文献
13.
目的对影响强生口腔护理产品Rembrandt软塑包装热封强度的因素进行研究。方法基于Ansys有限元方法分析内部气压大小对软塑管热封性能的影响。结果仿真结果表明,当管尾长度在44~54 mm之间,管身尺寸在65~70 mm之间,管壁厚度在0.4~0.5 mm之间时,软塑管的应变满足初始条件的要求;软塑包装的热封失效与材料的密度无关,随着弹性模量的增大,材料的最大应变值逐渐减小。结论要想获得理想的热封(良好的强度和外观)并满足生产效率,需要材料、结构、温度、压力、时间等多个方面共同作用。 相似文献
14.
针对SHDM结构(由依赖应变历史材料制成的结构)有限元非线性静力平衡方程组动力松弛法(DRM)迭代求解时的积分点应力更新步骤,提出非线性弹性增量算法,即在一个静力增量步内固定材料的加卸载路径,使之在该增量步内成为非线性弹性材料。该应力更新算法能使包括收敛解在内的迭代序列中不含虚假应变历史。此外,该算法还可避免静力解答精度依赖于静力增量步长的局限性。通过三个SHDM结构的数值试验对该算法进行了验证。该算法可望对SHDM结构非线性有限元静力问题DRM分析技术的发展起促进作用。 相似文献
15.
Antonio Huerta Sonia Fernndez‐Mndez 《International journal for numerical methods in engineering》2001,51(11):1361-1383
Volumetric locking (locking in the incompressible limit) for linear elastic isotropic materials is studied in the context of the element‐free Galerkin method. The modal analysis developed here shows that the number of non‐physical locking modes is independent of the dilation parameter (support of the interpolation functions). Thus increasing the dilation parameter does not suppress locking. Nevertheless, an increase in the dilation parameter does reduce the energy associated with the non‐physical locking modes; thus, in part, it alleviates the locking phenomena. This is shown for linear and quadratic orders of consistency. Moreover, the biquadratic order of consistency, as in finite elements, improves the locking behaviour. Although more locking modes are present in the element‐free Galerkin method with quadratic consistency than with standard biquadratic finite elements. Finally, numerical examples are shown to validate the modal analysis. In particular, the conclusions of the modal analysis are also confirmed in an elastoplastic example. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
16.
R. Hauptmann K. Schweizerhof S. Doll 《International journal for numerical methods in engineering》2000,49(9):1121-1141
In the present contribution we extend a previously proposed so‐called solid–shell concept which incorporates only displacement degrees of freedom to the simulation of large elastic and large elastoplastic deformations of shells. Therefore, the modifications necessary for hyper‐elastic or elastoplastic material laws are discussed. These modifications concern the right Cauchy–Green tensor for large elastic deformations, respectively, the deformation gradient for elastoplasticity which then are consistent to the modified Green–Lagrange strains that are necessary for transverse shear and membrane locking free solid–shell element formulations. However, in addition to the locking mentioned above especially in the range of plasticity incompressibility locking becomes important. Thus, the second major aspect of this contribution is the discussion of several ways to avoid incompressibility locking also including the investigation of eigenmodes. Finally, a selective reduced integration scheme with reduced integration for the volumetric term is employed and described in detail, although it is limited to material laws which allow the decomposition into a volumetric and a deviatoric part. Some numerical examples show the range of application for the proposed elements. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献
17.
M. Broccardo M. Micheloni P. Krysl 《International journal for numerical methods in engineering》2009,78(9):1113-1134
An assumed‐strain finite element technique for non‐linear finite deformation is presented. The weighted‐residual method enforces weakly the balance equation with the natural boundary condition and also the kinematic equation that links the elementwise and the assumed‐deformation gradient. Assumed gradient operators are derived via nodal integration from the kinematic‐weighted residual. A variety of finite element shapes fits the derived framework: four‐node tetrahedra, eight‐, 27‐, and 64‐node hexahedra are presented here. Since the assumed‐deformation gradients are expressed entirely in terms of the nodal displacements, the degrees of freedom are only the primitive variables (displacements at the nodes). The formulation allows for general anisotropic materials and no volumetric/deviatoric split is required. The consistent tangent operator is inexpensive and symmetric. Furthermore, the material update and the tangent moduli computation are carried out exactly as for classical displacement‐based models; the only deviation is the consistent use of the assumed‐deformation gradient in place of the displacement‐derived deformation gradient. Examples illustrate the performance with respect to the ability of the present technique to resist volumetric locking. A constraint count can partially explain the insensitivity of the resulting finite element models to locking in the incompressible limit. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
18.
Y. Onishi R. Iida K. Amaya 《International journal for numerical methods in engineering》2017,109(11):1582-1606
A new smoothed finite element method (S‐FEM) with tetrahedral elements for finite strain analysis of nearly incompressible solids is proposed. The proposed method is basically a combination of the F‐bar method and edge‐based S‐FEM with tetrahedral elements (ES‐FEM‐T4) and is named ‘F‐barES‐FEM‐T4’. F‐barES‐FEM‐T4 inherits the accuracy and shear locking‐free property of ES‐FEM‐T4. At the same time, it also inherits the volumetric locking‐free property of the F‐bar method. The isovolumetric part of the deformation gradient ( F iso) is derived from the F of ES‐FEM‐T4, whereas the volumetric part ( F vol) is derived from the cyclic smoothing of J(=det( F )) between elements and nodes. Some demonstration analyses confirm that F‐barES‐FEM‐T4 with a sufficient number of cyclic smoothings suppresses the pressure oscillation in nearly incompressible materials successfully with no increase in DOF. Moreover, they reveal that our method is capable of relaxing the corner locking issue arising at the corner in the cylinder barreling analysis. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
19.
M. A. Wheel 《International journal for numerical methods in engineering》1999,44(12):1843-1861
A finite volume formulation for determining small strain deformations in incompressible materials is presented in detail. The formulation includes displacement and hydrostatic pressure variables. The displacement field varies linearly along and across each cell face. The hydrostatic pressure field associated with each face is uniform. The cells that discretize the structure are geometrically unrestricted, each cell can have an arbitrary number of faces. The formulation is tested on a number of linear elastic plane strain benchmark problems. This testing reveals that when meshes of multifaceted cells are employed to represent the structure then locking behaviour is exhibited, but when triangular cells are used then accurate predictions of the displacement and stress fields are produced. Copyright © 1999 John Wiley & Sons, Ltd. 相似文献
20.
O. J. B. Almeida Pereira J. P. Moitinho de Almeida 《International journal for numerical methods in engineering》2001,50(4):885-898
In the assumed displacement, or primal, hybrid finite element method, the requirements of continuity of displacements across the sides are regarded as constraints, imposed using Lagrange multipliers. In this paper, such a formulation for linear elasticity, in which the polynomial approximation functions are not associated with nodes, is presented. Elements with any number of sides may be easily used to create meshes with irregular vertices, when performing a non‐uniform h‐refinement. Meshes of non‐uniform degree may be easily created, when performing an hp‐refinement. The occurrence of spurious static modes in meshes of triangular elements, when compatibility is strongly enforced, is discussed. An algorithm for the automatic selection, based on the topology of a mesh of triangular elements, of the sides in which to decrease the degree of the approximation functions, in order to eliminate all these spurious modes and preserve compatibility, is presented. A similar discussion is presented for the occurrence of spurious static modes in meshes of tetrahedral elements. An algorithm, based on heuristic criteria, that succeeded in eliminating these spurious modes and preserving compatibility in all the meshes of tetrahedral elements of uniform degree that were tested, is also presented. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献