• Journal of Magnetics
• /
• v.11 no.4
• /
• pp.151-155
• /
• 2006

The electromagnetic engine valve actuator is a key technology to achieve variable valve timing in internal combustion engine and the steel core and clapper of the electromagnetic engine valve actuator are laminated to reduce the eddy current loss. To design and characterize the performance of the electromagnetic engine valve actuator, FE (finite element) analysis is the most effective way, but FE (finite element) 3-D modeling of real lamination needs very fine meshes resulting in countless meshes for modeling and numerous computations. In this paper, the equivalent FE 2-D model of electromagnetic engine valve actuator is introduced and FE analysis is performed using the equivalent FE 2-D model.

• Steel and Composite Structures
• /
• v.49 no.5
• /
• pp.587-602
• /
• 2023

This article presents the numerical modelling of transient heat transfer in highly heterogeneous composite materials where the thermal conductivity, specific heat and density are assumed to be directional-dependent. This article uses a coupled finite element-finite difference scheme to perform the transient heat transfer analysis of unidirectional (1D) and multidirectional (2D/3D) functionally graded composite panels. Here, 1D/2D/3D functionally graded structures are subjected to nonuniform heat source and inhomogeneous boundary conditions. Here, the multidirectional functionally graded materials are modelled by varying material properties in individual or in-combination of spatial directions. Here, fully spatial-dependent material properties are evaluated using Voigt's micromechanics scheme via multivariable power-law functions. The weak form is obtained through the Galerkin method and solved further via the element-space and time-step discretisation through the 2D-isoparametric finite element and the implicit backward finite difference schemes, respectively. The present model is verified by comparing it with the previously reported results and the commercially available finite element tool. The numerous illustrations confirm the significance of boundary conditions and material heterogeneity on the transient temperature responses of 1D/2D/3D functionally graded panels.

• Interaction and multiscale mechanics
• /
• v.6 no.2
• /
• pp.237-255
• /
• 2013

In this paper, a meshfree-enriched finite element method (ME-FEM) is introduced for the large deformation analysis of nonlinear path-dependent problems involving contact. In linear ME-FEM, the element formulation is established by introducing a meshfree convex approximation into the linear triangular element in 2D and linear tetrahedron element in 3D along with an enriched meshfree node. In nonlinear formulation, the area-weighted smoothing scheme for deformation gradient is then developed in conjunction with the meshfree-enriched element interpolation functions to yield a discrete divergence-free property at the integration points, which is essential to enhance the stress calculation in the stage of plastic deformation. A modified variational formulation using the smoothed deformation gradient is developed for path-dependent material analysis. In the industrial metal forming problems, the mortar contact algorithm is implemented in the explicit formulation. Since the meshfree-enriched element shape functions are constructed using the meshfree convex approximation, they pose the desired Kronecker-delta property at the element edge thus requires no special treatments in the enforcement of essential boundary condition as well as the contact conditions. As a result, this approach can be easily incorporated into a conventional displacement-based finite element code. Two elasto-plastic problems are studied and the numerical results indicated that ME-FEM is capable of delivering a volumetric locking-free and pressure oscillation-free solutions for the large deformation problems in metal forming analysis.

• Transactions of Materials Processing
• /
• v.13 no.6 s.70
• /
• pp.484-489
• /
• 2004

This paper employed a systematic analysis using a 2-D hybrid special finite element containing an edge crack in order to describe the fracture behavior of spot-welds in automotive structures. The 2-D hybrid special finite element is derived form a mixed formulation with a complex potential function with the description of the singularity of a stress field. The hybrid special finite element containing an edge crack can give a better description of its singularity with only one hybrid element surrounding one crack. The advantage of this special element is that it can greatly simplify the numerical modeling of the spot welds. Some numerical examples demonstrate the validity and versatility of the present analysis method. The lap-shear, lap-tension and angle-clip specimens are analyzed and some useful fracture parameters such as the stress intensity factor and the initial direction of crack growth are obtained simultaneously.

• Transactions of Materials Processing
• /
• v.12 no.3
• /
• pp.194-201
• /
• 2003

An automatic quadrilateral mesh generator for large deformation finite element analysis such as metal forming simulation was developed. The NURBS interpolation method is used for modeling arbitrary 2-D free surface. This mesh generation technique is the modified paving algorithm, which is an advancing front technique with element-by-element resolving method for paving boundary intersection problem. The mesh density for higher analysis accuracy and less analysis time can be easily controlled with high-density points, maximum and minimum element size. A couple of application to large deformation finite element analysis is given as an example, which shows versatility and applicability of the proposed approach and the developed mesh generator for large deformation finite element analysis.

• Journal of the Korean Geotechnical Society
• /
• v.15 no.4
• /
• pp.221-234
• /
• 1999

This paper deals with the application of joint element in the finite element modeling of discontinuities encountered during rock tunneling. A nodal displacement joint element was implemented in a two dimensional finite element program GEOFE2D. The applicability of the joint element for modeling of discontinuities and the numerical stability of the implemented algorithm were examined by comparing the results of reduced small scale model tests as well as commercially available FEM program. The GEOFE2D was then used to analyze a tunnel crossed by a major discontinuity for the purpose of understanding the effect of discontinuity on the tunnel behavior. In addition, a modeling technique for the junction of discontinuity and shotcrete lining was presented. The results of analysis indicated that the stress-strain field around the tunnel is significantly altered by the presence of discontinuity, and that the stresses in the shotcrete lining considerably increase at the junction of the shotcrete lining and the discontinuity. It is therefore concluded that the major discontinuities must be carefully modeled in the finite element analysis of a tunneling problem in order to obtain more reliable results close to actual tunnel behavior.

• Structural Engineering and Mechanics
• /
• v.40 no.2
• /
• pp.257-277
• /
• 2011

In this paper the geometrically nonlinear continuum plate finite element model, hitherto not reported in the literature, is developed using the total Lagrange formulation. With the layerwise displacement field of Reddy, nonlinear Green-Lagrange small strain large displacements relations (in the von Karman sense) and linear elastic orthotropic material properties for each lamina, the 3D elasticity equations are reduced to 2D problem and the nonlinear equilibrium integral form is obtained. By performing the linearization on nonlinear integral form and then the discretization on linearized integral form, tangent stiffness matrix is obtained with less manipulation and in more consistent form, compared to the one obtained using laminated element approach. Symmetric tangent stiffness matrixes, together with internal force vector are then utilized in Newton Raphson's method for the numerical solution of nonlinear incremental finite element equilibrium equations. Despite of its complex layer dependent numerical nature, the present model has no shear locking problems, compared to ESL (Equivalent Single Layer) models, or aspect ratio problems, as the 3D finite element may have when analyzing thin plate behavior. The originally coded MATLAB computer program for the finite element solution is used to verify the accuracy of the numerical model, by calculating nonlinear response of plates with different mechanical properties, which are isotropic, orthotropic and anisotropic (cross ply and angle ply), different plate thickness, different boundary conditions and different load direction (unloading/loading). The obtained results are compared with available results from the literature and the linear solutions from the author's previous papers.

• Geomechanics and Engineering
• /
• v.17 no.4
• /
• pp.317-322
• /
• 2019

This paper presents the results of parametric analyses to compute the axial capacity of a suction pile using 2D and 3D finite element approaches. The study is intended to simplify the process of analyzing suction piles from 3D to 2D model. The research focuses on obtaining the coefficient to be applied into the 2D model in order to obtain results that are as close as possible to the 3D model. Two 2D models were used in the analysis, namely the plane strain and axisymmetric models. The analyses were performed using two actual offshore soil data of the North and West Java Indonesia. The study reveals that the simplification of model through 2D Finite Element is achievable by applying the appropriate coefficient to the stiffness parameters. The results show that the simplified model of the 2D FEA provides more conservative results (with the difference between 2% to 7%) than the 3D FEA.

• Earthquakes and Structures
• /
• v.9 no.1
• /
• pp.173-194
• /
• 2015

Site response analysis is an important topic in earthquake engineering. A time-domain numerical method called as one-dimensional (1D) finite element artificial boundary method is proposed to simulate the homogeneous plane elastic wave propagation in a layered half space subjected to the obliquely incident plane body wave. In this method, an exact artificial boundary condition combining the absorbing boundary condition with the inputting boundary condition is developed to model the wave absorption and input effects of the truncated half space under layer system. The spatially two-dimensional (2D) problem consisting of the layer system with the artificial boundary condition is transformed equivalently into a 1D one along the vertical direction according to Snell's law. The resulting 1D problem is solved by the finite element method with a new explicit time integration algorithm. The 1D finite element artificial boundary method is verified by analyzing two engineering sites in time domain and by comparing with the frequency-domain transfer matrix method with fast Fourier transform.

• Journal of Mechanical Science and Technology
• /
• v.20 no.3
• /
• pp.319-328
• /
• 2006

A simple approximation method for the stress intensity factor at the tip of the axial semielliptical cracks on the cylindrical vessel is developed. The approximation methods, incorporated in VINTIN (Vessel INTegrity analysis-INner flaws), utilizes the influence coefficients to calculate the stress intensity factor at the crack tip. This method has been compared with other solution methods including 3-D finite element analysis for internal pressure, cooldown, and pressurized thermal shock loading conditions. For these, 3-D finite-element analyses are performed to obtain the stress intensity factors for various surface cracks with t/R=0.1. The approximation solutions are within $\pm2.5%$ of the those of finite element analysis using symmetric model of one-forth of a vessel under pressure loading, and 1-3% higher under pressurized thermal shock condition. The analysis results confirm that the approximation method provides sufficiently accurate stress intensity factor values for the axial semi-elliptical flaws on the surface of the reactor pressure vessel.