首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
    
Whereas the higher-order versions of the finite element method (p- and hp-versions) are fairly well established as highly efficient methods for monitoring and controlling the discretization error in linear problems, little has been done to exploit their benefits in elasto-plastic structural analysis. In this paper, we discuss which aspects of incremental elasto-plastic finite element analysis are particularly amenable to improvements by the p-version. These theoretical considerations are supported by several numerical experiments. First, we study an example for which an analytical solution is available. It is demonstrated that the p-version performs very well even in cycles of elasto-plastic loading and unloading, not only as compared with the traditional h-version but also with respect to the exact solution. Finally, an example of considerable practical importance—the analysis of a cold-working lug—is presented which demonstrates how the modelling tools offered by higher-order finite element techniques can contribute to an improved approximation of practical problems.  相似文献   

2.
    
  相似文献   

3.
A finite element analysis model for material and geometrical non-linearities due to large plastic deformations of ductile materials is presented using the continuum damage mechanics approach. To overcome limitations of the conventional plastic analysis, a fourth-order tensor damage, defined in Part I of this paper to represent the stiffness degradation in the finite strain regime, is incorporated. General forms of an updated Lagrangian (U.L.) finite element procedure are formulated to solve the governing equations of the coupled elastic–plastic-damage analysis, and a computer program is developed for two-dimensional plane stress/strain problems. A numerical algorithm to treat the anisotropic damage is proposed in addition to the non-linear incremental solution algorithm of the U.L. formulation. Selected examples, compared with published results, show the validity of the presented finite element approach. Finally, the necking phenomenon of a plate with a hole is studied to explore plastic damage in large strain deformations. This revised version was published online in August 2006 with corrections to the Cover Date.  相似文献   

4.
    
  相似文献   

5.
    
An Erratum has been published for this article in International Journal for Numerical Methods in Engineering 2005, 63(8): 1228. We present a new formulation and a numerical procedure for the quasi‐static analysis of three‐dimensional crack propagation in brittle and quasi‐brittle solids. The extended finite element method (XFEM) is combined with linear tetrahedral elements. A viscosity‐regularized continuum damage constitutive model is used and coupled with the XFEM formulation resulting in a regularized ‘crack‐band’ version of XFEM. The evolving discontinuity surface is discretized through a C0 surface formed by the union of the triangles and quadrilaterals that separate each cracked element in two. The element's properties allow a closed form integration and a particularly efficient implementation allowing large‐scale 3D problems to be studied. Several examples of crack propagation are shown, illustrating the good results that can be achieved. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
  相似文献   

7.
    
We propose a new computational framework for the treatment of acousto‐magneto‐mechanical coupling that arises in low‐frequency electro‐magneto‐mechanical systems such as magnetic resonance imaging scanners. Our transient Newton–Raphson strategy involves the solution of a monolithic system obtained from the linearisation of the coupled system of equations. Moreover, this framework, in the case of excitation from static and harmonic current sources, allows us to propose a simple linearised system and rigorously motivate a single‐step strategy for understanding the response of systems under different frequencies of excitation. Motivated by the need to solve industrial problems rapidly, we restrict ourselves to solving problems consisting of axisymmetric geometries and current sources. Our treatment also discusses in detail the computational requirements for the solution of these coupled problems on unbounded domains and the accurate discretisation of the fields using hp–finite elements. We include a set of academic and industrially relevant examples to benchmark and illustrate our approach. Copyright © 2017 The Authors. International Journal for Numerical Methods in Engineering Published by John Wiley & Sons, Ltd.  相似文献   

8.
    
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.  相似文献   

9.
    
This paper presents improvements to three‐dimensional crack propagation simulation capabilities of the generalized finite element method. In particular, it presents new update algorithms suitable for explicit crack surface representations and simulations in which the initial crack surfaces grow significantly in size (one order of magnitude or more). These simulations pose problems in regard to robust crack surface/front representation throughout the propagation analysis. The proposed techniques are appropriate for propagation of highly non‐convex crack fronts and simulations involving significantly different crack front speeds. Furthermore, the algorithms are able to handle computational difficulties arising from the coalescence of non‐planar crack surfaces and their interactions with domain boundaries. An approach based on moving least squares approximations is developed to handle highly non‐convex crack fronts after crack surface coalescence. Several numerical examples are provided, which illustrate the robustness and capabilities of the proposed approaches and some of its potential engineering applications. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

10.
    
In this paper the formulation of an electric–mechanical beam‐to‐beam contact element is presented. Beams with circular cross‐sections are assumed to get in contact in a point‐wise manner and with clean metallic surfaces. The voltage distribution is influenced by the contact mechanics, since the current flow is constricted to small contacting spots. Therefore, the solution is governed by the contacting areas and hence by the contact forces. As a consequence the problem is semi‐coupled with the mechanical field influencing the electric one. The electric–mechanical contact constraints are enforced with the penalty method within the finite element technique. The virtual work equations for the mechanical and electric fields are written and consistently linearized to achieve a good level of computational efficiency with the finite element method. The set of equations is solved with a monolithic approach. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

11.
    
A mixed membrane eight-node quadrilateral finite element for the analysis of masonry walls is presented. Assuming that a nonlinear and history-dependent 2D stress-strain constitutive law is used to model masonry material, the element derivation is based on a Hu-Washizu variational statement, involving displacement, strain, and stress fields as primary variables. As the behavior of masonry structures is often characterized by strain localization phenomena, due to strain softening at material level, a discontinuous, piecewise constant interpolation of the strain field is considered at element level, to capture highly nonlinear strain spatial distributions also within finite elements. Newton's method of solution is adopted for the element state determination problem. For avoiding pathological sensitivity to the finite element mesh, a novel algorithm is proposed to perform an integral-type nonlocal regularization of the constitutive equations in the present mixed formulation. By the comparison with competing serendipity displacement-based formulation, numerical simulations prove high performances of the proposed finite element, especially when coarse meshes are adopted.  相似文献   

12.
扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。  相似文献   

13.
    
In this paper the numerical solution for ideal irrotational incompressible flow around axisymmetric 3D obstacles is discussed with the artificial boundary method. By introducing an artificial boundary, we divide the exterior unbounded domain into a bounded part and an unbounded part. After setting up a proper artificial boundary condition, we get an approximate reduced problem defined on the bounded part. Both non‐local and local artificial boundary conditions are designed. Numerical experiment is also presented, and its result demonstrates the effectiveness of these artificial boundary conditions. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

14.
    
A computational procedure for remapping material state information from one finite element mesh to another is described. The procedure is useful in connection with evolving meshes for inelastic problems, as for example occur in the context of fracture simulation and adaptive mesh refinement. The proposed method is based on weak enforcement of equality between corresponding fields on the two meshes, where piecewise‐constant fields on both meshes are generalized from the quadrature‐point values. The essential algorithmic problem is that of calculating the volume partition of an arbitrary convex region with respect to a covering set of disjoint convex regions. Instead of geometrically resolving the associated intersections, the problem is herein approximated by a constrained optimization problem, which may be readily and efficiently solved computationally. This formulation is a main contribution of the paper. Computational examples are given that illustrate the effectiveness of the proposed procedure. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

15.
    
This paper presents a coupling method between a discrete element code CeaMka3D and a finite element code Sem. The coupling is based on a least‐squares method, which adds terms of forces to finite element code and imposes the velocity at coupling particles. For each coupling face, a small linear system with a constant matrix is solved. This method remains conservative in energy and shows good results in applications. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

16.
    
In this paper we introduce and numerically solve a mathematical model for numerical simulation of electro‐magnetic field in a three‐phase electric reduction furnace. The model allows us to compute the current distribution on a cross‐section of the three electrodes. A combined boundary element/finite element method is used. Numerical results for real industrial furnaces are shown. As a by‐product we compute the torque on the electrodes due to the Lorentz electromagnetic force. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

17.
    
A general local/nonlocal implicit coupling technique called the dual-based approach is proposed to couple peridynamics (PD) with classical continuum mechanics. In the present method, physical information is transmitted mutually from local to nonlocal regions through the coupling elements; no transition region is introduced. For different mesh discretizations, two coupling methods are achieved with simplicity and effectivity. To obtain the stiffness matrix of the coupled model, without loss of generality, the implicit dual-horizon ordinary state-based peridynamic model is proposed, in which the linearization of dual-horizon ordinary state-based PD is derived and the dual assembly algorithm of the peridynamic stiffness matrix is developed. It will be seen that the implicit dual-based coupling approach provides a new implicit coupling method that is easy to implement and makes full use of the internal connection between PD and classical continuum mechanics. Several numerical examples involving static crack propagation are investigated, and the satisfactory results show both quantitative and qualitative agreement with either the analytic solution or the available experiment.  相似文献   

18.
    
This paper extends the double porosity deformable model previously developed by one of the authors to include temperature effects. Unlike the double temperature approach, the model assumes a single representative thermodynamics continuum which includes both conduction and convection heat transfer in a fully saturated porous media. Examples are presented to evaluate the suitability of the model. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

19.
         下载免费PDF全文
In order to understand the underlying mechanisms of inelastic material behavior and nonlinear surface interactions, which can be observed on macroscale as damping, softening, fracture, delamination, frictional contact etc., it is necessary to examine the molecular scale. Force fields can be applied to simulate the rearrangement of chemical and physical bonds. However, a simulation of the atomic interactions is very costly so that classical molecular dynamics (MD) is restricted to structures containing a low number of atoms such as carbon nanotubes. The objective of this paper is to show how MD simulations can be integrated into the finite element method (FEM) which is used to simulate engineering structures such as an aircraft panel or a vehicle chassis. A new type of finite element is required for force fields that include multi-body potentials. These elements take into account not only bond stretch but also bending, torsion and inversion without using rotational degrees of freedom. Since natural lengths and angles are implemented as intrinsic material parameters, the developed molecular dynamic finite element method (MDFEM) starts with a conformational analysis. By means of carbon nanotubes and elastomeric material it is demonstrated that this pre-step is needed to find an equilibrium configuration before the structure can be deformed in a succeeding loading step.  相似文献   

20.
    
A computational model for simulation of surface pitting of mechanical elements subjected to rolling and sliding contact conditions is presented. The two-dimensional computational model is restricted to modelling of high-precision mechanical components with fine surface finishing and good lubrication, where the cracks leading to pitting are initiated in the area of largest contact stresses at certain depth under the contacting surface. Hertz contact conditions with addition of friction forces are assumed and the position and magnitude of the maximum equivalent stress is determined by the finite element method. When the maximum equivalent stress exceeds the local material strength, it is assumed that the initial crack develops along the slip line in a single-crystal grain. The Virtual Crack Extension method in the framework of finite element analysis is then used for two-dimensional simulation of the fatigue crack propagation under contact loading from the initial crack up to the formation of the surface pit. The pit shapes and relationships between the stress intensity factor and crack length are determined for various combinations of contacting surface curvatures and loadings. The model is applied to simulation of surface pitting of two meshing gear teeth. Numerically predicted pit shapes in the face of gear teeth show a good agreement with the experimental observations. © 1998 John Wiley & Sons, Ltd.  相似文献   

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

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