Finite Element Analysis on Electromechanical Behaviours in Bone Tissues

According to the importance of electromechancial behavior of bone on clinical medical treatment, the physiological structure of bone tissue is studied. Then the biphasic porous media theory associated with electro kinetics is used to model the electromechanical effects in bone tissue. The finite element governing equations are built up through Galerkin method and the penalty method. The distortion of skeleton, interstitial fluid flow and the streaming potential in bone tissue are numerically computed and simulated.     

Mixed Finite Element Method for Dynamic Problems of Fluid Saturated Biphasic Porous Media

Based on the fluid saturated biphasic porous media model deduced from mixture theory,a finite element formulation with u S u F p variables for dynamic response analysis is given out.This method overcomes the difficulty of choosing suitable penalty parameter in penalty finite element method,and the accuracy of pressure distribution obtained with the mixed finite element method is higher than the penalty finile element method.An iterative solution method is suggested to solve the system of equations whose coefficient matrices are non positive definite.It is concluded from patch test that the order of interpolation function for pressure variables must be higher than that of displacement variables of solid phase.     

Numerical Simulation of Mechanical Behaviors of Articular Cartilages

Articular cartilage, which is composed with a solid phase of collagen fibers and proteoglycan and a pore fluid phase, can be depicted with two-phase porous medium model based on mixture theory. This paper analyzed the creep and stress relaxation responses of confined compression problems of articular cartilage with the developed finite element method. The obtained responses of velocities of both solid and fluid phases and the effective stresses of solid phase as well as the pore pressures during the creep and stress relaxation processes of articular cartilage tissues under confined compression can be used as a theoretical reference for man-made cartilage materials, and behave important significance for the research in bio-mechanical engineering area.     

Hyperelastic Constitutive Relations and Finite Element Formulation for Porous Silicone Rubber

The strain energy density function which was decoupled into isochoric parts and volumetric parts is presented for describing mechanical behaviors of compressible porous silicone rubber material,the constitutive model is then attained.The finite element formulation for isotropic hyperelastic porous silicone rubber is developed by use of total Lagrangean method.     

THE APPLICATION OF FINITE - ELEMENT METHODS FOR THE ANALYSIS OF GALLBLADDER BILE FLOW IN VIVO

The transient flow fields of gallgladder bile in vivo are calculated by finite element methods. Computer tomograph, ultrasound B and contrast radiography are applied to determine the moving boundary conditions, The penalty formulation and moving triangle elements have been successfully adopted and be flow fields for axisymmetric Newtonian bile flow have been obtained. The results show that there exists obvions reverse flow at the cross section (Perpendicular to Center axis) and thus forms axisymmetric vartexes inside a contracting gallbladder in vivo,Which well agree with the previous prediction.     

A Finite Element Method for Quasi Static Problems of Two phase Porous Media

For the quasi static problems of two phase porous media,the velocity variable of fluid phase can be eliminated as the viscosity of fluid is neglected.The corresponding governing equations,with displacement of solid phase and pore pressure as basic variables,as well as the boundary and initial conditions are given out.Finally a mixed finite element formula with u S-p variables is obtained with weighted residual method.This method has lesser nodal variables and higher accuracy of pressure distribution.     

Investigation on flow and heat transfer characteristics in capillary porous wick based on alkali metal thermal-to-electric converter (AMTEC)

An axial symmetric invariable temperature phase change interface model of capillary porous wick based on alkali metal thermal-to-electric converter (AMTEC) is established to simulate flow and heat transfer characteristics in capillary porous wick and liquid channel by solving the mass and heat transport controlling equations. The effects of the working fluid flow rate, inlet temperature, porous wick thickness, porosity and effective pore diameter on the distributions of pressure, velocity and temperature are analyzed. In addition, the available effective pore size of the porous wick and the measures to improve the performance of the porous wick are given by analyzing the relationship between the maximum available capillary pressure and the loop pressure drop. The results show that the model presented in gives a good agreement with the three dimensional two-phase evaporation model. The flow and heat transfer characteristics of the capillary porous wick with the alkali metal liquid as working fluid are different from that with the traditional working fluid in some ways.     

Corrected extended finite element method for two phase flow

An extended finite element method based on SUPG/PSPG is proposed to simulate the two phase flow problems. A corrected XFEM is introduced to ensure the blending element to satisfy the Partition of Unity in processing the discontinuity of the interface. Level set method is adopted to track the kinetic phase interface as the fluid flow. Free oscillation is numerically simulated with the presented method. The obtained numerical results are in consistent with the analytical and experimental results. Additionally, breaking dam problem is considered, the numerical solutions agree well with experimental results which illustrate the correctness and efficiency of the proposed method. No re-meshing is needed during the simulation of the two phase fluid flow and the moving interface can be accurately tracked by means of the presented method.     

Computational Simulation on Temperature Field of Slab Continuous Casting Mold

On the core concepts that the distributing of temperature and solidified slab shell in mold are a stable process, a new method to calculate the temperature field in mold has been promoted. Consequently, a coupling model on fluid filed and temperature field in stable state has also been established. The concept of effective heat capacity is introduced to deal with the latent heat source term from phase transformation. During the programming for computation, the dynamic updating method is also applied to simulate the influence of solidification on physical properties. The resolution results of this model are analyzed and compared with the tested thickness of solidified shell from casting broke out slab, which demonstrates the correctness and applicability of this model to the practical continuous casting process. Additionally, the thickness distribution of solidified shell in slab mold are also studied in detail.     

Finite Element Analyses of Spherical Node of Coal Pipeline System of Blast Furnace

The mechanical model of spherical node of coal pipelines of blast furnace is built. Its boundary conditions and loads are presented. Its finite element analyses are used for this structure. Its stress distributing rules are obtained and the security of this structure is evaluated.     

A FINITE ELEMENT RESEARCH ON.A LINKAGE OF A COMPOSITE STRUCTURE OF SHELLS AND PLATES

In this paper,the penalty element stiffness matrix of linkage of plates and shells is deduced accrding to the principle of the penalty elementmethod.This stiffness matrix can be also used for a linkage between shells andsolids.It is of some practical value in the finite element calculations of engineeringstructure.     

Mixed Finite Element Formulation of the Linear Biphasic Theory for Soft Tissues

The behavior of soft tissues of human musculoskeletal system can be describedwith the biphasic model based on a continuum theory of mktures. This paper. using Galerkin weight-ed residual method.obtains a mixed finite element formulation for the linear biphasic model of smalldeformation. and. in turn. gives out the iterative scheme solving the system equations. The results ofnumerical analysis for the constrained compression problem are consistent with those obtained bytheory. which illustrates the correctnas and feasibility of the derived mixed finite element formula-tion. concludingly. this formulation provides an effective means of numerical analysis for the motfonmechanism of human articulating joints.     

THE THREE-DIMENSIONAL CALCULATING MODEL OF INVERSE PROBLEM IN ECG AND A SPECIAL BOUNDRYVALUE PROBTEM

In this paper, the mathematical model and the three-dimensional finite clement formulation of inverse problem in electrocardiography (ECG) under a general inhomogcneous and anisotropic configuration of the torso conductor is focused on. The following points are especially considered: 1. The mathamatical model is described by a differential equation instead of the Fredholm integral equation. 2. A kind of special boundary value problem is defined as ill-posed boundary value problem. 3. A group of finite element formulation, in which the method that treats anisotropic medium, is built up. 4. The calculating method of ill-posed boundary value problem is discussed.     

A Treatment to the Isoline for Finite Element Solution to Ferrite-Martensite Microstructure

The numerical solution of sample point stress field in finite element network of ferrite-martensite microstructure under finite deformation has a high-gradient distributive characteristic whose stress gradient between ferrite and martensite crystalline element is very high. If the interpolation method is used in the whole region,the isogram will be seriously smooth so that the original high-gradient sample point stress field will become untrue. This paper presents a new method for drawing isogram of ferrite-martensite microstructiire stress field by interpolating sample points in an element one by one. This method can make the post-treatment of high-gradient discrete stress field in crystalline network approach to a truer condition and solve the problem of arbitrary curve boundary of the deformed two phase steel.     

Coupled Vibration of Saturated Fractional Derivative Type Viscoelastic Soil of a Circular Tunnel with Partially Sealed Shell Lining

Considering the compression of medium and fluid, coupled vibration of saturated fractional derivative type viscoelastic soil and a circular tunnel with partially sealed shell lining in the frequency domain is investigated according to theories of Biot and elastic shell. The stress and displacement constitutive behavior of the soil skeleton is described by fractional derivative model which reflects the rheological properties of the medium while regarding soil as a liquid-solid saturated porous medium. The expressions of displacement, stress and pore water pressure are obtained while the inner boundary of circular tunnel is subjected to axially symmetric radial traction and axially symmetric fluid pressure respectively by introducing a partially permeable boundary condition. With the parameter analysis, it is revealed that the order of fractional derivative model on the responses for the system subjected to the symmetric radial traction is much greater than that of the system under the axially symmetric fluid pressure. And resonance phenomenon occurs obviously. Nevertheless the system responses do not have remarkable resonance phenomenon under axially symmetric fluid pressure.     

Numerical Simulation of Flow field and Temperature field in Combustion Chamber Widened Suddenly

The flowfield and temperature field are numerically Simulated by the mathematical model of the fluid flow and heat transfer in circular combustion chamber widened suddenly. The distributions and influence factors of the flow field and temperature field are studied in it. The results have the guiding significance for the design of combustion chamber. The mathematical models can serve as the quantitative analysis methods for design.     

Deformation of Soft Ground with Consideration of the Concomitant Damage

On the basis of continuum damage mechanics theory, the effective stress theory of damaged porous media was presented and its modified self contained seepage differential equation was deduced. The traditional Drucker Prager yield criterion subsequently was modified, and a coupled seepage and damage finite element computer program compiled that was combined with the modified Biot consolidation equation. A simple roadbed foundation model then was analyzed and calculated using the computer program we compiled, and the typical node settlement, pore pressure, and damage variables were discussed and compared with those from the analytic solution. The results show that the damage finite element program we compiled is reasonable and effective.     

The Finite Element Simulation of Pipe Bend Forming

A dynamic explicit finite element method is used to pipe bend forming processes on bender in this paper.The degenerated shell element and co-rotational coordinate system are employed.The material of pipe is assumed to satisfy Hill's anisotropic elasto-plastic yield condition.The bend die and press die are assumed as rigid surfaces and the contact force is computed by penalty method.Several numerical examples are given.They show that this method is effective.     

Study on the Reduction of Thermocapillary Convection of the Single-crystal Growth of Silicon under Microgravity

Float-zone method for growth of semi-conductor silicon is considering very promising under microgravity,but the effect of thermocapillary convection induced by surface tension gradient on the crystal quality is quite highlight.In the present work,a nwe method referred to as surface-cut method is suggested to reduce the melt flow.A mathmatical model to describe the thermocapillary convection in a half float-zone with surface-cut is proposed and numerical simulation for the temperature and the flow field under different surface-cut way is conducted using the finite element method.The results show that the effective reduction can be achieved and the flow will be weaked by 70% by proper increase of the surface-cut time.     

Simulation on parallel conditions for small single phase synchronous generators

The controller of paralleling small single phase synchronous generators should restrict the dash currents in generators below a defined value when it switches on. A field-circuit coupling time-stepping finite element model of the parallel system is built. The computation program is written in APDL languages of ANSYS and the transient process of paralleling two single phase synchronous generators in various parallel parameters are studied. The value and variation of the dash currents in two generators when switch on at various differences of three parallel parameters are analyzed and the accurate parallel conditions are provided to the parallel controller.