Research of the Higher-order Finite Element Arithmetic for Radiation Exchange

With the development of aeronautic and astronautic techniques, radiation becomes much more significant while the structure is exposed to the higher and higher temperature. Most of the current finite element software packages treat it using the net-radiation method or absorbed radiation method based on the assumption of isothermal surface with uniform radiation heat flux, which brings the conflict between the precision and the quantity of grids. Using integral method to compute the variable radiation heat flux in higher-order finite element, the precision can be improved greatly while using the same quantity of grids, because it is more consistent with the distribution of real temperature. In this paper, the integral is only processed on the same integral points as those used for solving the finite element equations, so it may be of high efficiency. In an academic testing model, the result is contrast to which get in ANSYS, proving the high precision of the method. Then an actual sandwich panel used in the thermal protection system is analyzed with the method, and the error is comparatively low to the analytical answer while the computation being of high efficiency.

Experimental research and finite element analysis of bridge piers failed in flexure-shear modes

In recent earthquakes, a large number of reinforced concrete （RC） bridges were severely damaged due to mixed flexure-shear failure modes of the bridge piers. An integrated experimental and finite element （FE） analysis study is described in this paper to study the seismic performance of the bridge piers that failed in flexure-shear modes. In the first part, a nonlinear cyclic loading test on six RC bridge piers with circular cross sections is carried out experimentally. The damage states, ductility and energy dissipation parameters, stiffness degradation and shear strength of the piers are studied and compared with each other. The experimental results suggest that all the piers exhibit stable flexural response at displacement ductilities up to four before exhibiting brittle shear failure. The ultimate performance of the piers is dominated by shear capacity due to significant shear cracking, and in some cases, rupturing of spiral bars. In the second part, modeling approaches describing the hysteretic behavior of the piers are investigated by using ANSYS software. A set of models with different parameters is selected and evaluated through comparison with experimental results. The influences of the shear retention coefficients between concrete cracks, the Bauschinger effect in longitudinal reinforcement, the bond-slip relationship between the longitudinal reinforcement and the concrete and the concrete failure surface on the simulated hysteretic curves are discussed. Then, a modified analysis model is presented and its accuracy is verified by comparing the simulated results with experimental ones. This research uses models available in commercial FE codes and is intended for researchers and engineers interested in using ANSYS software to predict the hysteretic behavior of reinforced concrete structures.

A Method Combining Numerical Analysis and Limit Equilibrium Theory to Determine Potential Slip Surfaces in Soil Slopes

This paper describes a precise method combining numerical analysis and limit equilibrium theory to determine potential slip surfaces in soil slopes. In this method, the direction of the critical slip surface at any point in a slope is determined using the Coulomb’s strength principle and the extremum principle based on the ratio of the shear strength to the shear stress at that point. The ratio, which is considered as an analysis index, can be computed once the stress field of the soil slope is obtained. The critical slip direction at any point in the slope must be the tangential direction of a potential slip surface passing through the point. Therefore, starting from a point on the top of the slope surface or on the horizontal segment outside the slope toe, the increment with a small distance into the slope is used to choose another point and the corresponding slip direction at the point is computed. Connecting all the points used in the computation forms a potential slip surface exiting at the starting point. Then the factor of safety for any potential slip surface can be computed using limit equilibrium method like Spencer method. After factors of safety for all the potential slip surfaces are obtained, the minimum one is the factor of safety for the slope and the corresponding potential slip surface is the critical slip surface of the slope. The proposed method does not need to pre-assume the shape of potential slip surfaces. Thus it is suitable for any shape of slip surfaces. Moreover the method is very simple to be applied. Examples are presented in this paper to illustrate the feasibility of the proposed method programmed in ANSYS software by macro commands.

FINITE ELEMENT NUMERICAL SIMULATION OF TEMPERATURE FIELD IN METAL PATTERN CASTING SYSTEM AND "REVERSE METHOD" OF DEFINING THE THERMAL PHYSICAL COEFFICIENT

With the technology support of virtual reality and ANSYS software, an example on the simulation of temperature distribution of casting system during the solidification process was provided, which took the latent heat of phase change, the conditions for convection, and the interface heat transfer coefficient into consideration. The result of ANSYS was found to agree well with the test data. This research offers an unorthodox way or ＂reverse method＂ of defining the relevant thermal physical coefficient.

Two-dimensional Numerical Modeling Research on Continent Subduction Dynamics

Continent subduction is one of the hot research problems in geoscience. New models presented here have been set up and two-dimensional numerical modeling research on the possibility of continental subduction has been made with the finite element software, ANSYS, based on documentary evidence and reasonable assumptions that the subduction of oceanic crust has occurred, the subduction of continental crust can take place and the process can be simplified to a discontinuous plane strain theory model. The modeling results show that it is completely possible for continental crust to be subducted to a depth of 120 km under certain circumstances and conditions. At the same time, the simulations of continental subduction under a single dynamical factor have also been made, including the pull force of the subducted oceanic lithosphere, the drag force connected with mantle convection and the push force of the mid-ocean ridge. These experiments show that the drag force connected with mantle convection is critical for continent subduction.

Influences of Randomly Distributed Wall Thickness of Beverage Can on Its Strength and Stiffness

This paper describes the research undertaken on the strength and stiffness of fluctuation on the wall thickness of steel beverage cans using the Monte Carlo stochastic finite element method. Sample distributions were firstly assumed and then proven using the data observations of the wall thickness, the APDL language was then applied, and the stresses and displacements of the can were calculated by using the ANSYS software. It is concluded that the structural reliability of a steel making beverage can be estimated accurately.

Seismic Response Analysis of Portal Water Injection Sheet Pile

To further the study on the newly developed portal water injection sheet pile under static loads, in this paper, by adopting the nonlinear calculation module of FEM software ANSYS, a model for the interaction between the soil and the sheet piles is set up, and the seismic response analysis for this type of space-retaining structure is performed. The effects of the embedded depth and the distance between the front pile and the back pile on the dynamic characteristics of the portal water injection sheet pile are studied.

Optimization and Simulation of Inclined Channel Area Lining for CDQ Ovens

The mechanical force borne by the lintel and“brackets”in the inclined channel area of CDQ and the thermal stress generated by the top-down temperature gradient were simulated by ANSYS software.The following conclusions are drawn.(1)The mechanical force has very mild effect on the damage of the“brackets”.(2)The temperature gradient caused by the uneven temperature distribution of CDQ oven is the key factor affecting the“bracket”damage.The different expansion of refractory materials in various parts results in the thermal stress concentration which tears the“brackets”.(3)The structure adjustment of the“brackets”has very little effect on the thermal stress distribution.(4)The multi-phase nitride bonded silicon carbide for CDQ can effectively alleviate the stress concentration at the“brackets”,realizing the long-life and stable operation.

3-D Finite-Element Numerical Simulation of Centrifugal Induction Electrosalg Casting Solidification Process

A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.

Analysis of Demagnetizer Operation Using the Finite Element Method

Permanent magnets with high energy products are widely used in a variety of electromagnetic devices. Such devices can be found in marine, aerospace, and robotic applications which require the minimization of weight and volume of the electromagnetic device. During the magnetizing process, the magnet may not be perfectly magnetized. Therefore, it needs to be demagnetized. Because of high coercivity of some permanent magnets, the demagnetization process requires the intense magnetic fields in close proximity with the magnetic material. The fields must be produced for a short period of time （millisecond range） and they also must be bidirectional in order to overcome the coercivity of magnetic material. Different parameters have been known to affect the demagnetizer operation such as the core shape, core material, turn number, cross section of air gap, magnet type and so on. The amplitude and waveform of the fixture current is obtained from PSPICE simulations and also from experimental measurement. A 2D finite element analysis is developed to simulate the magnetic fields and the mechanical forces. In this paper, the main parameters affecting the optimal design of the demagnetizer are discussed. Simulation results show that the core structure, air gap width, and turn numbers are the most important parameters when designing such a device.