共查询到20条相似文献，搜索用时 0 毫秒
1.
The Institute of Civil Engineers, U.K., have announced the winner of the 2008 O. C. Zienkiewicz Award. The Zienkiewicz Numerical Methods in Engineering Prize was introduced in 1998 by John Wiley & Sons Ltd. publishers of the International Journal of Numerical Methods in Engineering. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
2.
Xufeng Xiao Xinlong Feng Zhilin Li 《International journal for numerical methods in engineering》2019,120(7):901-917
In this paper, we present an adaptive mesh refinement method for solving convection-diffusion-reaction equations on surfaces, which is a fundamental subproblem in many models for simulating the transport of substances on biological films and solid surfaces. The method considered is a combination of well-known techniques: the surface finite element method, streamline diffusion stabilization, and the gradient recovery–based Zienkiewicz-Zhu error estimator. The streamline diffusion method overcomes the instability issue of the finite element method for the dominance of the convection. The gradient recovery–based adaptive mesh refinement strategy enables the method to provide high-resolution numerical solutions by relatively fewer degrees of freedom. Moreover, the implementation detail of a surface mesh refinement technique is presented. Various numerical examples, including the convection-dominated diffusion problems with large variations of solutions, nearly singular solutions, discontinuous sources, and internal layers on surfaces, are presented to demonstrate the efficacy and accuracy of the proposed method. 相似文献
3.
4.
Vasco Varduhn Ming‐Chen Hsu Martin Ruess Dominik Schillinger 《International journal for numerical methods in engineering》2016,107(12):1054-1079
The finite cell method (FCM) is an immersed domain finite element method that combines higher‐order non‐boundary‐fitted meshes, weak enforcement of Dirichlet boundary conditions, and adaptive quadrature based on recursive subdivision. Because of its ability to improve the geometric resolution of intersected elements, it can be characterized as an immersogeometric method. In this paper, we extend the FCM, so far only used with Cartesian hexahedral elements, to higher‐order non‐boundary‐fitted tetrahedral meshes, based on a reformulation of the octree‐based subdivision algorithm for tetrahedral elements. We show that the resulting TetFCM scheme is fully accurate in an immersogeometric sense, that is, the solution fields achieve optimal and exponential rates of convergence for h‐refinement and p‐refinement, if the immersed geometry is resolved with sufficient accuracy. TetFCM can leverage the natural ability of tetrahedral elements for local mesh refinement in three dimensions. Its suitability for problems with sharp gradients and highly localized features is illustrated by the immersogeometric phase‐field fracture analysis of a human femur bone. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
5.
Carsten Carstensen Kerstin Weinberg 《International journal for numerical methods in engineering》2003,56(15):2313-2330
Adaptive algorithms are important tools for efficient finite‐element mesh design. In this paper, an error controlled adaptive mesh‐refining algorithm is proposed for a non‐conforming low‐order finite‐element method for the Reissner–Mindlin plate model. The algorithm is controlled by a reliable and efficient residual‐based a posteriori error estimate, which is robust with respect to the plate's thickness. Numerical evidence for this and the efficiency of the new algorithm is provided in the sense that non‐optimal convergence rates are optimally improved in our numerical experiments. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
6.
Mofdi El‐Amrani Mohammed Seaïd 《International journal for numerical methods in engineering》2012,90(12):1485-1507
High‐order accurate methods for convection‐dominated problems have the potential to reduce the computational effort required for a given order of solution accuracy. The state of the art in this field is more advanced for Eulerian methods than for semi‐Lagrangian (SLAG) methods. In this paper, we introduce a new SLAG method that is based on combining the modified method of characteristics with a high‐order interpolating procedure. The method employs the finite element method on triangular meshes for the spatial discretization. An L2 interpolation procedure is developed by tracking the feet of the characteristic lines from the integration nodes. Numerical results are illustrated for a linear advection–diffusion equation with known analytical solution and for the viscous Burgers’ equation. The computed results support our expectations for a robust and highly accurate finite element SLAG method. Copyright © 2012 John Wiley & Sons, Ltd. 相似文献
7.
本文给出了极坐标系下弹性力学平面问题的Hamilton正则方程,并提出一种求解该方程的状态空间有限元法。文中通过对Hellinger-Reissner混和变分原理的修正,导出了Hamilton正则方程及其对应能量泛函,然后采用分离变量法对其场变量进行分离变量,这样就可在θ方向采用通常的有限元插值,而沿半径r方向采用状态空间法给出解析解答,从而实现了有限元法与控论制中状态空间的结合。通过计算表明,本文方法精度高。 相似文献
8.
Javier de Frutos Julia Novo 《International journal for numerical methods in engineering》2005,63(8):1146-1173
We show that the issue of a posteriori estimate the errors in the numerical simulation of non‐linear parabolic equations can be reduced to a posteriori estimate the errors in the approximation of an elliptic problem with the right‐hand side depending on known data of the problem and the computed numerical solution. A procedure to obtain local error estimates for the p version of the finite element method by solving small discrete elliptic problems with right‐hand side the residual of the p‐FEM solution is introduced. The boundary conditions are inherited by those of the space of hierarchical bases to which the error estimator belongs. We prove that the error in the numerical solution can be reduced by adding the estimators that behave as a locally defined correction to the computed approximation. When the error being estimated is that of a elliptic problem constant free local lower bounds are obtained. The local error estimation procedure is applied to non‐linear parabolic differential equations in several space dimensions. Some numerical experiments for both the elliptic and the non‐linear parabolic cases are provided. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
9.
10.
G. Ferté P. Massin N. Moës 《International journal for numerical methods in engineering》2014,100(11):834-870
The aim of this paper is to propose a procedure to accurately compute curved interfaces problems within the extended finite element method and with quadratic elements. It is dedicated to gradient discontinuous problems, which cover the case of bimaterials as the main application. We focus on the use of Lagrange multipliers to enforce adherence at the interface, which makes this strategy applicable to cohesive laws or unilateral contact. Convergence then occurs under the condition that a discrete inf‐sup condition is passed. A dedicated P1 multiplier space intended for use with P2 displacements is introduced. Analytical proof that it passes the inf‐sup condition is presented in the two‐dimensional case. Under the assumption that this inf‐sup condition holds, a priori error estimates are derived for linear or quadratic elements as functions of the curved interface resolution and of the interpolation properties of the discrete Lagrange multipliers space. The estimates are successfully checked against several numerical experiments: disparities, when they occur, are explained in the literature. Besides, the new multiplier space is able to produce quadratic convergence from P2 displacements and quadratic geometry resolution. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
11.
J. Kim A. Simone C. A. Duarte 《International journal for numerical methods in engineering》2017,109(2):235-258
A robust and efficient strategy is proposed to simulate mechanical problems involving cohesive fractures. This class of problems is characterized by a global structural behavior that is strongly affected by localized nonlinearities at relatively small‐sized critical regions. The proposed approach is based on the division of a simulation into a suitable number of sub‐simulations where adaptive mesh refinement is performed only once based on refinement window(s) around crack front process zone(s). The initialization of Newton‐Raphson nonlinear iterations at the start of each sub‐simulation is accomplished by solving a linear problem based on a secant stiffness, rather than a volume mapping of nonlinear solutions between meshes. The secant stiffness is evaluated using material state information stored/read on crack surface facets which are employed to explicitly represent the geometry of the discontinuity surface independently of the volume mesh within the generalized finite element method framework. Moreover, a simplified version of the algorithm is proposed for its straightforward implementation into existing commercial software. Data transfer between sub‐simulations is not required in the simplified strategy. The computational efficiency, accuracy, and robustness of the proposed strategies are demonstrated by an application to cohesive fracture simulations in 3‐D. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
12.
Mostafa E. Mobasher Haim Waisman 《International journal for numerical methods in engineering》2016,105(8):599-619
We propose a coupled boundary element method (BEM) and a finite element method (FEM) for modelling localized damage growth in structures. BEM offers the flexibility of modelling large domains efficiently, while the non‐linear damage growth is accurately accounted by a local FEM mesh. An integral‐type nonlocal continuum damage mechanics with adapting FEM mesh is used to model multiple damage zones and follow their propagation in the structure. Strong form coupling, BEM hosted, is achieved using Lagrange multipliers. Because the non‐linearity is isolated in the FEM part of the system of equations, the system size is reduced using Schur complement approach, then the solution is obtained by a monolithic Newton method that is used to solve both domains simultaneously. The coupled BEM/FEM approach is verified by a set of convergence studies, where the reference solution is obtained by a fine FEM. In addition, the method is applied to multiple fractures growth benchmark problems and shows good agreement with the literature. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
13.
14.
F. K. F. Radtke A. Simone L. J. Sluys 《International journal for numerical methods in engineering》2010,84(6):708-732
The numerical analysis of large numbers of arbitrarily distributed discrete thin fibres embedded in a continuum is a computationally demanding process. In this contribution, we propose an approach based on the partition of unity property of finite element shape functions that can handle discrete thin fibres in a continuum matrix without meshing them. This is made possible by a special enrichment function that represents the action of each individual fibre on the matrix. Our approach allows to model fibre‐reinforced materials considering matrix, fibres and interfaces between matrix and fibres individually, each with its own elastic constitutive law. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
15.
Vitaly E. Bulgakov Marina V. Bulgakova 《International journal for numerical methods in engineering》1998,43(3):533-548
A finite element constructed on the basis of boundary integral equations is proposed. This element has a flexible shape and arbitrary number of nodes. It also has good approximation properties. A procedure of constructing an element stiffness matrix is demonstrated first for one-dimensional case and then for two-dimensional steady-state heat conduction problem. Numerical examples demonstrate applicability and advantages of the method. © 1998 John Wiley & Sons, Ltd. 相似文献
16.
17.
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. 相似文献
18.
Thu Hang Vu Andrew J. Deeks 《International journal for numerical methods in engineering》2006,65(10):1714-1733
The scaled boundary finite element method is a novel semi‐analytical technique, whose versatility, accuracy and efficiency are not only equal to, but potentially better than the finite element method and the boundary element method for certain problems. This paper investigates the possibility of using higher‐order polynomial functions for the shape functions. Two techniques for generating the higher‐order shape functions are investigated. In the first, the spectral element approach is used with Lagrange interpolation functions. In the second, hierarchical polynomial shape functions are employed to add new degrees of freedom into the domain without changing the existing ones, as in the p‐version of the finite element method. To check the accuracy of the proposed procedures, a plane strain problem for which an exact solution is available is employed. A more complex example involving three scaled boundary subdomains is also addressed. The rates of convergence of these examples under p‐refinement are compared with the corresponding rates of convergence achieved when uniform h‐refinement is used, allowing direct comparison of the computational cost of the two approaches. The results show that it is advantageous to use higher‐order elements, and that higher rates of convergence can be obtained using p‐refinement instead of h‐refinement. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
19.
Jian‐Gang Han Wei‐Xin Ren Yih Huang 《International journal for numerical methods in engineering》2006,66(1):166-190
The wavelet‐based methods are powerful to analyse the field problems with changes in gradients and singularities due to the excellent multi‐resolution properties of wavelet functions. Wavelet‐based finite elements are often constructed in the wavelet space where field displacements are expressed as a product of wavelet functions and wavelet coefficients. When a complex structural problem is analysed, the interface between different elements and boundary conditions cannot be easily treated as in the case of conventional finite‐element methods (FEMs). A new wavelet‐based FEM in structural mechanics is proposed in the paper by using the spline wavelets, in which the formulation is developed in a similar way of conventional displacement‐based FEM. The spline wavelet functions are used as the element displacement interpolation functions and the shape functions are expressed by wavelets. The detailed formulations of typical spline wavelet elements such as plane beam element, in‐plane triangular element, in‐plane rectangular element, tetrahedral solid element, and hexahedral solid element are derived. The numerical examples have illustrated that the proposed spline wavelet finite‐element formulation achieves a high numerical accuracy and fast convergence rate. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
20.
Q. Tong H. L. Wang R. Xu B. Liu D.N. Fang 《International journal for numerical methods in engineering》2014,98(6):445-468
Finite element method (FEM) with fixed representative volume element (RVE) encounters some difficulties in simulating the periodical postbuckling behaviors of infinite long beam or infinite large film on soft substrate under compression, because the wavelength and pattern of buckling are not known before simulation and will change with the increase of compression strain. In this paper, an adaptive periodical RVE is constructed in a mapping space to avoid remeshing in the real space, and the mapping coefficients, that is, the dimension and shape of RVE in real space, are treated as variables in average energy density minimization to obtain correct postbuckling configurations. The validness and efficiency of the proposed algorithm have been demonstrated by our numerical examples. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献