Three-dimensional analysis of elastic stress distribution of indented ceramic surface by finite element method

The three-dimensional stress distributions in the area surrounding indentation pattern for three different materials, Al2O3, Si3N4 and SiC were analyzed by finite element method(FEM). Those theoretical results were also compared with the experimental ones by Rockwell hardness test. The effect of loading stress on the plastic deformation in specimens, surface was investigated on the assumption of shear strain energy theory by Huber-Mises when the materials were indented. The distributions of nomal stress, shear stress, and Mises stress were analysed with variations of loading conditions. It is clear that the analytical results for the stress distributions, the crack length and its density of probability are in good agreement with the experimental results.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

Finite element analysis of three-dimensional laser-induced transient thermal grating in diamond/ZnSe system

This paper uses finite element method to obtain the three-dimensional temperature field of laser-induced transient thermal grating （TTG） for two-layered structure of diamond film on ZnSe substrate. The numerical results indicate that unique two-times heating process is gradually experienced in the area between two adjacent grating stripes. However, there is a little change for the temperature field along the depth direction for the diamond film due to its great thermal conductivity. It further finds that the thickness of the diamond film has a significant influence on the temperature field in diamond/ZnSe system. The results are useful for the application of laser-induced TTG technique in film/substrate system.

Three-Dimensional Thermo-Elastic-Plastic Finite Element Method Modeling for Predicting Weld-Induced Residual Stresses and Distortions in Steel Stiffened-Plate Structures

The objective of the present paper is to develop nonlinear finite element method models for predicting the weld-induced initial deflection and residual stress of plating in steel stiffened-plate structures. For this purpose, three-dimensional thermo-elastic-plastic finite element method computations are performed with varying plate thickness and weld bead length (leg length) in welded plate panels, the latter being associated with weld heat input. The finite element models are verified by a comparison with experimental database which was obtained by the authors in separate studies with full scale measurements. It is concluded that the nonlinear finite element method models developed in the present paper are very accurate in terms of predicting the weld-induced initial imperfections of steel stiffened plate structures. Details of the numerical computations together with test database are documented.

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of one dimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.

Combining the complex variable reproducing kernel particle method and the finite element method for solving transient heat conduction problems

In this paper, the complex variable reproducing kernel particle （CVRKP） method and the finite element （FE） method are combined as the CVRKP-FE method to solve transient heat conduction problems. The CVRKP-FE method not only conveniently imposes the essential boundary conditions, but also exploits the advantages of the individual methods while avoiding their disadvantages, then the computational efficiency is higher. A hybrid approximation function is applied to combine the CVRKP method with the FE method, and the traditional difference method for two-point boundary value problems is selected as the time discretization scheme. The corresponding formulations of the CVRKP-FE method are presented in detail. Several selected numerical examples of the transient heat conduction problems are presented to illustrate the performance of the CVRKP-FE method.

Distributed Lagrange Multiplier/Fictitious Domain Finite Element Method for a Transient Stokes Interface Problem with Jump Coefficients

The distributed Lagrange multiplier/fictitious domain(DLM/FD)-mixed finite element method is developed and analyzed in this paper for a transient Stokes interface problem with jump coefficients.The semi-and fully discrete DLM/FD-mixed finite element scheme are developed for the first time for this problem with a moving interface,where the arbitrary Lagrangian-Eulerian(ALE)technique is employed to deal with the moving and immersed subdomain.Stability and optimal convergence properties are obtained for both schemes.Numerical experiments are carried out for different scenarios of jump coefficients,and all theoretical results are validated.

New nonconforming finite element method for solving transient Naiver-Stokes equations

For transient Naiver-Stokes problems, a stabilized nonconforming finite element method is presented, focusing on two pairs inf-sup unstable finite element spaces, i.e., pNC/pNC triangular and pNQ/pNQ quadrilateral finite element spaces. The semi- and full-discrete schemes of the stabilized method are studied based on the pressure projection and a variational multi-scale method. It has some attractive features： avoiding higher-order derivatives and edge-based data structures, adding a discrete velocity term only on the fine scale, being effective for high Reynolds number fluid flows, and avoiding increased computation cost. For the full-discrete scheme, it has second-order estimations of time and is unconditionally stable. The presented numerical results agree well with the theoretical results.

Parametric study on single shot peening by dimensional analysis method incorporated with finite element method

Shot peening is a widely used surface treatment method by generating compressive residual stress near the surface of metallic materials to increase fatigue life and re- sistance to corrosion fatigue, cracking, etc. Compressive re- sidual stress and dent profile are important factors to eval- uate the effectiveness of shot peening process. In this pa- per, the influence of dimensionless parameters on maximum compressive residual stress and maximum depth of the dent were investigated. Firstly, dimensionless relations of pro- cessing parameters that affect the maximum compressive residual stress and the maximum depth of the dent were de- duced by dimensional analysis method. Secondly, the in- fluence of each dimensionless parameter on dimensionless variables was investigated by the finite element method. Fur- thermore, related empirical formulas were given for each di- mensionless parameter based on the simulation results. Fi- nally, comparison was made and good agreement was found between the simulation results and the empirical formula, which shows that a useful approach is provided in this pa- per for analyzing the influence of each individual parameter.

Validation and application of three-dimensional discontinuous deformation analysis with tetrahedron finite element meshed block

In the last decade, three dimensional discontin- uous deformation analyses （3D DDA） has attracted more and more attention of researchers and geotechnical engineers worldwide. The original DDA formulation utilizes a linear displacement function to describe the block movement and deformation, which would cause block expansion under rigid body rotation and thus limit its capability to model block de- formation. In this paper, 3D DDA is coupled with tetrahe- dron finite elements to tackle these two problems. Tetrahe- dron is the simplest in the 3D domain and makes it easy to implement automatic discretization, even for complex topol- ogy shape. Furthermore, element faces will remain planar and element edges will remain straight after deformation for tetrahedron finite elements and polyhedral contact detection schemes can be used directly. The matrices of equilibrium equations for this coupled method are given in detail and an effective contact searching algorithm is suggested. Valida- tion is conducted by comparing the results of the proposed coupled method with that of physical model tests using one of the most common failure modes, i.e., wedge failure. Most of the failure modes predicted by the coupled method agree with the physical model results except for 4 cases out of the total 65 cases. Finally, a complex rockslide example demon- strates the robustness and versatility of the coupled method.

Finite Element Method for Transient Electric Field by Using Indirect Laplace Transform with High Accuracy

This paper focuses on the finite element method in the complex frequency domain(CFD-FEM)for the transient electric field.First,the initial value and boundary value problem of the transient electric field under the electroquasistatic field in the complex frequency domain is given.In addition,the finite element equation and the constrained electric field equation on the boundary are derived.Secondly,the indirect algorithm of the numerical inverse Laplace transform is introduced.Based on it,the calculation procedures of the CFD-FEM are illustrated in detail.Thirdly,the step response,zero-state response under the positive periodic square waveform(PPSW)voltage,and the zero-input response by the CFD-FEM with direct algorithm and indirect algorithm are compared.Finally,the reason for the numerical oscillations of the zero-state response under the PPSW voltage is analyzed,and the method to reduce oscillations is proposed.The results show that the numerical accuracy of the indirect algorithm of the CFD-FEM is more than an order of magnitude higher than that of the direct algorithm when calculating the step response of the transient electric field.The proposed method can significantly reduce the numerical oscillations of the zero-state response under the PPSW voltage.The proposed method is helpful for the calculation of the transient electric field,especially in the case of frequency-dependent parameters.

Numerical Investigation of Thermal Behavior of CNC Machine Tool and Its Effects on Dimensional Accuracy of Machined Parts

The dimensional accuracy of machined parts is strongly influenced by the thermal behavior of machine tools (MT). Minimizing this influence represents a key objective for any modern manufacturing industry. Thermally induced positioning error compensation remains the most effective and practical method in this context. However, the efficiency of the compensation process depends on the quality of the model used to predict the thermal errors. The model should consistently reflect the relationships between temperature distribution in the MT structure and thermally induced positioning errors. A judicious choice of the number and location of temperature sensitive points to represent heat distribution is a key factor for robust thermal error modeling. Therefore, in this paper, the temperature sensitive points are selected following a structured thermomechanical analysis carried out to evaluate the effects of various temperature gradients on MT structure deformation intensity. The MT thermal behavior is first modeled using finite element method and validated by various experimentally measured temperature fields using temperature sensors and thermal imaging. MT Thermal behavior validation shows a maximum error of less than 10% when comparing the numerical estimations with the experimental results even under changing operation conditions. The numerical model is used through several series of simulations carried out using varied working condition to explore possible relationships between temperature distribution and thermal deformation characteristics to select the most appropriate temperature sensitive points that will be considered for building an empirical prediction model for thermal errors as function of MT thermal state. Validation tests achieved using an artificial neural network based simplified model confirmed the efficiency of the proposed temperature sensitive points allowing the prediction of the thermally induced errors with an accuracy greater than 90%.

Analysis Process of Finite Element Method in Heat Transfer through Fabrics

According to heat transfer principle and the process of solving engineering problems by finite element method, examples were given to demonstrate how finite element analysis can be used to describe transient heat transfer through fabrics. Details were given to describe how conduction and convection affect temperature distribution and heat loss during heat transfer processes by taking advantage of the quick calculation of FEA software MSC.Marc. Experimental results show good agreement with the theoretical results.

An Efficient Multiple-Dimensional Finite Element Solution for Water Flow in Variably Saturated Soils

Multiple-dimensional water flow in variably saturated soils plays an important role in ecological systems such as irrigation and water uptake by plant roots; its quantitative description is usually based on the Richards＇ equation. Because of the nonlinearity of the Richards＇ equation and the complexity of natural soils, most practical simulations rely on numerical solutions with the nonlinearity solved by iterations. The commonly used iterations for solving the nonlinearity are Picard and Newton methods with the former converging at first-order rate and the later at second-order rate. A recent theoretical analysis by the authors, however, revealed that for solving the diffusive flow, the classical Picard method is actually a chord-Newton method, converging at a rate faster than first order; its linear convergence rate is due to the treatment of the gravity term. To improve computational efficiency, a similar chord-Newton method as for solving the diffusive term was proposed to solve the gravity term. Testing examples for one-dimensional flow showed significant improvement. The core of this method is to produce a diagonally dominant matrix in the linear system so as to improve the iteration-toiteration stability and hence the convergence. In this paper, we develop a similar method for multiple-dimensional flow and compare its performance with the classical Picard and Newton methods for water flow in soils characterised by a wide range of van Genuchten parameters.

Three dimensional tectonic stress field in North China

According to the latest data of geological structure, geophysics, in-situ stress measurement and focal mechanism,3-D tectonic stress field model in North China is built and 3-D tectonic stress field pattern of North China aresimulated by finite element method. Then the overall characteristics and regional specific feature of North Chinaare studied. Finally, the influences of the valid dynamic boundary conditions of North China Block, active faultsand the inhomogeneity of crustal medium on tectonic stress field of North China are investigated.

Forward simulation of coal fire model by transient electromagnetic method with electrical source

Uncontrolled coal fires are natural disasters that may cause mineral loss and environmental damage.The traditional loop source transient electromagnetic method can effectively detect the low-resistivity region of coal fires,but its detection efficiency is not so good for high-resistivity regions.In view of this limitation,a technique based on electrical source transient electromagnetics is proposed in this paper to detect high-resistivity regions in the spontaneous combustion process of coal.Considering the complex geometry of the coal fire area,an unstructured tetrahedral grid is used in this study to realize the spatial discretization of the model,and solve the electromagnetic field based on a vector finite element algorithm.Numerical analysis is used to investigate methods for detecting coal fires and the characteristics of effective anomalies are further examined to provide guidance for practical detection.

Three-dimensional stochastic seepage field for embankment engineering

Owing to the complexity of geo-engineering seepage problems influenced by different random factors, three-dimensional simulation and analysis of the stochastic seepage field plays an important role in engineering applications. A three-dimensional anisotropic heterogeneous steady random seepage model was developed on the basis of the finite element method. A statistical analysis of the distribution characteristics of soil parameters sampled from the main embankment of the Yangtze River in the Southern Jingzhou zone of China was conducted. The Kolomogorov-Smirnov test verified the statistical hypothesis that the permeability coefficient tensor has a Gaussian distribution. With the help of numerical analysis of the stochastic seepage field using the developed model, various statistical and random characteristics of the stochastic seepage field of the main embankment of the Yangtze River in the Southern Jingzhou zone of China were investigated. The model was also examined with statistical testing. Through the introduction of random variation of the upstream and downstream water levels into the model, the effects of the boundary randomness due to variation of the downstream and upstream water levels on the variation of simulated results presented with a vector series of the random seepage field were analyzed. Furthermore, the combined influence of the variation of the soil permeability coefficient and such seepage resistance measures as the cut-off wall and relief ditch on the hydraulic head distribution was analyzed and compared with the results obtained by determinate analysis. Meanwhile, sensitivities of the hydraulic gradient and downstream exit height to the variation of boundary water level were studied. The validity of the simulated results was verified by stochastic testing and measured data. The developed model provides more detail and a full stochastic algorithm to characterize and analyze three-dimensional stochastic seepage field problems.

Algebraic Calculation Method of One-Dimensional Steady Compressible Gas Flow

This paper describes a new method of calculation of one-dimensional steady compressible gas flows in channels with possible heat and mass exchange through perforated sidewalls. The channel is divided into small elements of a finite size for which mass, energy and momentum conservation laws are written in the integral form, assuming linear distribution of the parameters along the length. As a result, the calculation is reduced to finding the roots of a quadratic algebraic equation, thus providing an alternative to numerical methods based on differential equations. The advantage of this method is its high tolerance to coarse discretization of the calculation area as well as its good applicability for transonic flow calculations.

Finite Element Simulation of the Polymeric Liquid Flow in Coat-hanger Die

A three - dimenslonal finite element model is developed to deal with the polymeric liquid flow in coat - hanger die. This model is used to predict the flow behavior of the 2％ CMC/watsr solution in the coat - hanger die with linearly taper manifolds and its validity is experimentally verified quantitatively and qualitatively by using Laser Doppler Velocimetry and Particle Image Velocimetry respectively.