Neumann stochastic finite element method for calculating temperature field of frozen soil based on random field theory

To study the effect of uncertain factors on the temperature field of frozen soil, we propose a method to calculate the spatial average variance from just the point variance based on the local average theory of random fields. We model the heat transfer coefficient and specific heat capacity as spatially random fields instead of traditional random variables. An analysis for calculating the random temperature field of seasonal frozen soil is suggested by the Neumann stochastic finite element method, and here we provide the computational formulae of mathematical expectation, variance and variable coefficient. As shown in the calculation flow chart, the stochastic finite element calculation program for solving the random temperature field, as compiled by Matrix Laboratory （MATLAB） sottware, can directly output the statistical results of the temperature field of frozen soil. An example is presented to demonstrate the random effects from random field parameters, and the feasibility of the proposed approach is proven by compar- ing these results with the results derived when the random parameters are only modeled as random variables. The results show that the Neumann stochastic finite element method can efficiently solve the problem of random temperature fields of frozen soil based on random field theory, and it can reduce the variability of calculation results when the random parameters are modeled as spatial- ly random fields.

Application of Wavelet Finite Element Method to Simulation of the Temperature Field of Copier Paper

Simulation of the temperature field of copier paper in copier fusing is very important for improving the fusing property of reprography. The temperature field of copier paper varies with a high gradient when the copier paper is moving through the fusing rollers. By means of conventional shaft elements, the high gradient temperature variety causes the oscillation of the numerical solution. Based on the Daubechies scaling functions, a kind of wavelet based element is constructed for the above problem. The temperature field of the copier paper moving through the fusing rollers is simulated using the two methods. Comparison of the results shows the advantages of the wavelet finite element method, which provides a new method for improving the copier properties.

Finite element simulation of three-dimensional temperature field in underwater welding

Mathematical models of three-dimensional temperature fields in underwater welding with moving heat sources are built. Double ellipsoid Gauss model is proposed as heat sources models. Several factors which affect the temperature fields of underwater welding are analyzed. Water has little influence on thermal efftciency. Water convection coefftcient varies with the temperature difference between the water and the workpiece , and water convection makes molten pool freeze quickly. With the increase of water depth, the dimensions of heat sources model should be reduced as arc shrinks. Finite element technology is used to solve mathematical models. ANSYS software is used as finite element tool, and ANSYS Parametric Design Language is used to develop subprograms for loading the moving heat sources and the various convection coefftcients. Experiment results show that computational results by using double ellipsoid Gauss heat sources model accord well with the experimental results.

Numerical Simulation of Temperature Field and Thermal Stress Field of Work Roll During Hot Strip Rolling

Based on the thermal conduction equations, the three-dimensional （3D） temperature field of a work roll was investigated using finite element method （FEM）. The variations in the surface temperature of the work roll during hot strip rolling were described, and the thermal stress field of the work roll was also analyzed. The results showed that the highest roll surface temperature is 593 ℃, and the difference between the minimum and maximum values of thermal stress of the work roll surface is 145.7 MPa. Furthermore, the results of this analysis indicate that temperature and thermal stress are useful parameters for the investigation of roll thermal fatigue and also for improving the quality of strip during rolling.

THE COUPLED FEM ANALYSIS OF THE TRANSIENT TEMPERATURE FIELD DURING INERTIA FRICTION WELDING OF GH4169 ALLOY

The inertia friction welding process is a non-linear process because of the interaction between the temperature field and the material properties as well as the friction force. A thermo-mechanical coupled finite element model is established to simulate the temperature field of this process. The transient temperature distribution during the inertia friction welding process of two similar workpieces of GH4169 alloy is calculated. The region of the circular cross-section of the workpiece is divided into a number of four-nodded isoparametric elements. In this model, the temperature dependent thermal properties, time dependent heat inputs, contact condition of welding interface, and deformation of the flash were considered. At the same time, the convection and radiation heat losses at the surface of the workpieces were also considered. A temperature data acquisition system was developed. The temperature at some position near the welding interface was measured using this system. The calculated temperature agrees well with the experimental data. The deformation of the flash and the factor affecting the temperature distribution at the welding interface are also discussed.

Numerical simulation of temperature fields for T-joint during TIG welding of titanium alloy

Three-dimensional finite element model was established to simulate temperature fields of T-joint titanium sheets during TIG welding with finite element method （FEM） software. Temperature dependent material properties and the effect of latent heat were considered. A technique of element birth and death was used to simulate the process of welded metal filling. Dynamic variation process of temperature fields during T1G welding was achieved. The simulated results agreed well with the measured results.

FEM simulations and experimental studies of the temperature field in a large diamond crystal growth cell

We investigate the temperature field variation in the growth region of a diamond crystal in a sealed cell during the whole process of crystal growth by using the temperature gradient method （TGM） at high pressure and high temperature （HPHT）. We employ both the finite element method （FEM） and in situ experiments. Simulation results show that the temperature in the center area of the growth cell continues to decrease during the process of large diamond crystal growth. These results are in good agreement with our experimental data, which demonstrates that the finite element model can successfully predict the temperature field variations in the growth cell. The FEM simulation will be useful to grow larger high-quality diamond crystal by using the TGM. Furthermore, this method will be helpful in designing better cells and improving the growth process of gem-quality diamond crystal.

Effects of Process Parameters on the Temperature Field in Ti-6Al-4V Alloy Blade Precision Forging Process

Blade precision forging is a high temperature and large plastic deformation process. Process parameters have a great effect on temperature distribution in billet, so in this paper, by taking a Ti-6Al-4V alloy blade with a tenon as an object, the influence of process parameters on the temperature distribution in precision forging process was investigated using 3D coupled thermo-mechanical FEM （finite element method） code developed by the authors. The results obtained illustrate that： （1） the gradient of temperature distribution increases with increasing the deformation degree; （2） with increasing the initial temperature of the billet, the zones of high temperature become larger, and the gradient of temperature distribution hardly has any increase; （3） friction factors have little effect on the distribution of temperature field; （4） with increasing upper die velocity, temperature of the billet increases while the temperature gradient in billet decreases. The results are helpful to the design and optimization of the process parameters in precision forging process of Ti-alloy blade.

3D Eddy Current and Temperature Field Analysis of Large Hydro-generators in Leading Phase Operations

As a common practice,a large hydro-generator will operate in leading phase conditions to absorb the reactive power of the power grid.However,the accurate and precise prediction of the leading phase operation capacity of a large hydro-generator has always been a formidable challenge to engineers and academicians because it is extremely hard to compute the eddy currents and losses as well as the local overheating in the pressure plate and finger.To address this problem,a full three dimensional(3D)finite element model and method of the coupled eddy current and temperature fields in the end region of a large hydro-generator are developed.The equivalent medium parameters used in the computations are comprehensively discussed.Moreover,some numerically based solution methodologies for accurate computation of the field and armature currents under different leading phase conditions are proposed.Numerical results on the coupled eddy current and temperature fields in the end regions of a 250 MW hydro-generator confirm positively the feasibility of the present work.

Numerical Analysis of Silicon Micromachined Gas Pendulum Tilt Sense Organ Temperature Field

An analysis of the sensitive mechanism of silicon micromachined gas pendulum tilt sense organ is made.Adopting the method of FEA (finite element analysis),the temperature field at two points heat source,when the two-dimensional enclosure was inclined,was obtained by application of the program ANSYS-FLOTRAN CFD and a series of procedures,such as modeling,meshing,loading and equation solving.The numerical results show that in the level state,the temperatures at two points heat source are two points in the same isotherm;however,the temperatures are not the same when the enclosure is inclined.The difference of the temperatures will increase with the augment of the tilt angle,and contrarily it will decrease.That is the characteristic used to sense the transformation of obliquity.

Simulation of local effect of reinforcement on temperature field during solidification of aluminum metal matrix composite

The effect of reinforcement on the solidification of pure metal matrix composites (MMCs) was simulated using a two-dimensional solidification temperature field model by the finite element method. The concept of the character length was proposed to describe the size of reinforcement local heat influential zone in MMCs solidification according to the change of the morphologies of solid-liquid interface. The relationship between the character length and the geometrical conditions, the boundary condition and physical properties of the reinforcement were studied, respectively. The results show that the width of the unit and the cold boundary temperature have no effect on the character lengths but have effect on the distance between cold boundary and reinforcement (l) and the thermal parameters of the reinforcement. An experimental rule to predict the value of the character length was derived and applied.

Simulation of Temperature and Flow Field by Three Dimension Finite ElementMethod for Castex Process of AS Wire

Castex of AS wire is a new technology of near net shape. To study the variation of temperature and velocity of liquid (or semisolid) aluminum during dynamic solidification the numerical simulation was carried out with the theory of heat-transfer and hydrodynamics by means of 3-dimensional finite element method. From simulation results, it is found that the variation of temperature and velocityis mainly influenced by the casting temperature of aluminum, rotating speed of Castex wheel and flow of cooling water. Among theseinfluencing factors, the casting temperature distributes most to the length of liquid phase metal. Moreover, the faster the metal solidifies,the higher the metal there moves with the overall trend of descending from the bottom of the wheel to the shoe wall as well as from sidewalls to the center of wheel groove. In comparison with the practical value, the simulation is reliable.

An On-line Finite Element Temperature Field Model for Plate Ultra Fast Cooling Process

Taking the element specific-heat interpolation function into account, a one-dimensional （l-D） finite ele- ment temperature field model for the on-line control of the ultra fast cooling process was developed based on the heat transfer theory. This 1-D model was successfully implemented in one 4 300 mm plate production line. To improve the calculation accuracy of this model, the temperature-dependent material properties inside an element were considered during the modeling process. Furthermore, in order to satisfy the real-time requirements of the on-line model, the variable bandwidth storage method and the Cholesky decomposition method were used in the programming to storage the data and carry out the numerical solution. The on-line application of the proposed model indicated that the devia- tion between the calculated cooling stop temperature and the measured one was less than ± 15 ℃.

Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.

Plastic wrinkling model and characteristics of shear enforced Ti-alloy thin-walled tubes under combination die constraints and differential temperature fields

Plastic wrinkling predictions and shear enforced wrinkling characteristics of Ti-alloy thin-walled tubes under combination die constraints have become key problems urgently in need of solutions in order to improve forming quality in their shear bending processes under differential temperature fields. To address this, a wrinkling wave function was developed by considering their shear bend deformation characteristics. Based on this wave function and the thin shell theory, an energy prediction model for this type of wrinkling was established. This model enables consideration of the effects of shear deformation zone ranges, material parameters, loading modes, and friction coefficients between tube and dies on the minimum wrinkling energy. Tube wrinkling sensitive zones(WSZs) can be revealed by combining this wrinkling prediction model with a thermalmechanical coupled finite element model for simulating these bending processes. The reliability of this wrinkling prediction model was verified, and an investigation into the tube wrinkling characteristics was carried out based on the experimental conditions. This found that the WSZs are located on either a single side or both sides of the maximum shear stress zone. When the friction coefficients between the tube and the various dies coincide, the WSZs are located on both sides.The larger the value of the tube inner corner radius and/or the smaller the value of the outer cornerradius, the smaller the wrinkling probability. With an increase in the value of the moving die displacement, the wrinkling probability increases at first, and then decreases.

FINITE ELEMENT ANALYSIS OF AXIAL FEED BAR ROLLING

A flexible technique of hot working of bars by axial feed rolling was introduced. The process deformation, strain field, stress field, and temperature fie M of the parts are analyzed by finite element method （ FEM）-simulation software DEFORM-3D. The material flow rule and tool load have been investigated.

COUPLED NUMERICAL SIMULATION ON COLD ROLLER'S TEMPERATURE FIFLD-PHASE TRANSFORMATION - STRESS FIELD DURING ITS QUENCHING PSOCESS

Complicated changns occur inside the steel parts during quenching process. The abruptly changed boundary conditions make the temperature field,micro - structure and stress field change dramatically in very short time, and these variables take a contact interactions in the whole process. In this paper, a three dimensional non - linear mathematical model for queeching process has been founded and the numerical simulation on temperature field,microstructre and stress field has been realized.In the FEM analysis, the incremental iteration method is used to deal with such complicated nonlinear as boundary nonlinear, physical property nonlinear,transformation nonlinear etc.The effect of stress on transformation kinetics has been considered in the calculation of microstructure. In the stress field anal- ysis,a thermo- elasto - plastic model has been founded, which considers such factors as transforma- tion strain,transformation plastic strain, themal strain and the effect of temperature and transforma- tion on mechanical propertier etc. The transient temperature field, microstructure distribution and stress field of the roller on any time can be displayed vividly,and the cooling curve and the changes of stress on any position can also be given.

Numerical analysis for random temperature fields of embankment in cold regions

The stochastic finite element equations for random temperature are obtained using the first-order per-turbation technique taking into account the random thermal properties and boundary condition, based on heat transfer variational principle. The local average method for 2-D is used to discretize random fields. Then, the random temperature fields of embankment in cold regions are investigated on condi-tion that the thermal properties and boundary condition are taken as random fields, respectively, by using the program, which is written by the methods. The expected value of temperature field and the standard deviation of the temperature field of embankment in cold regions are obtained and analyzed.

Finite element analysis of the effect of micro-pore defect on linear friction welding of medium carbon steel

Micro-pore is a very common material defect. In the present paper, the temperature fields of medium carbon steel joints with and without micro-pore defect during linear friction welding （LFW） were investigated by using finite element method. The effect of micro-pore defect on the axial shortening of joints during LFW was examined. The x- and y-direction displacements of micro-pore during the LFW process were also studied. In addition, the shape of micro-pore after LFW was observed. The heat conducted from the weld inteace to the specimen interior. The fluctuation range of the temperature curves for the joint with micro-pore is larger than that without micro-pore. Position of micro-pore changes with the change of the friction time. The circular shape of micro-pore becomes oval after welding.