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.     

A Three Dimensional Torso Model for the Numerical Simulation of ECG--Forward Problem of ECG(1)

In this paper,a numerical model of a three dimensional human torso has been presented.In this model,the torso part including anisotropic muscle layers has been discretized by means of the finite element method(FEM).The rest of torso have been divided into the surface elements by the boundary element method(BEM).This discretized model can be used in the field such as the forward problem and the inverse problem of the electrocardiogram(ECG).     

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.     

Numerical simulation of the internal dynamic excitation of a spiral bevel gear transmission

A 3-D geometry model and a static/dynamic finite element (FE) model of a spiral bevel gear transmission are established. Using ANSYS software, the mesh stiffness excitation of a tooth pair is calculated via the static contact FE method. Using LS     

Computation method of magnetic field inverse problem on grounding grids fault diagnosis

Grounding grids are the guarantee of substation running safety. The fault of grounding grids is a threat to people and equipment in substations. A method to diagnose the fault of grounding grids based on magnetic field inverse problem. The inverse problem of magnetic field is established by injecting and extracting a low frequency current. The morbid of inverse problem is solved through regularization and Newton iteration method to obtain the distribution of the mesh currents. The approximate distribution of the magnetic field is calculated to diagnose the fault of the grounding grid. A 4×4 substation grounding grids is taken as an example. By comparison of the forward problem result, the error is under 3%, so the prposed method of inverse problem is verified. The selection of measurement point is discussed, and 40% of the number of meshes is accounted to approximate the real distribution of the magnetic field on the grounding grids.     

Study on the Nonlinear Vibration System with Slow-changing Parameters

In this paper a posteriori error estimates for Galerkin approximation of general operator equations is firstly presented in the framework of Sobolev spaces. Then a practical posteriori error estimates formula for the adaptive boundary element method solving the acoustic scattering problem with a finite plane screen is obtained by triangulations. The mathematical model of this problem is the three dimensional Neumann boundary value problem of Helmholtz equation with finite plane boundary.     

Lateral Bearing Capacity and Simplified Equations of Dry Stack In filled Reinforcement Concrete Frame Structure

The mechanism and lateral bearing capacity of reinforcement concrete (RC) frame with dry stack in filled panel (DSIP) were investigated using quasi static experiments and finite element models. According to the parallel model and equivalent strut model, the lateral bearing capacity of RC frame and DSIP were researched separately. Results show that: 1) The plain stress element and interface element are applicable in finite element (FE) model analysis. According to the FE model, the failure of RC frame with DSIP is caused by the damage of frame; 2) Lateral bearing capacity of DSIP is mainly from the friction between bricks in the panel, which can be divided into 3 stages: constant stage, increasing stage and ultimate stage; 3) Equations for lateral bearing capacity of DSIP were proposed and verified by FE model results.     

基于超静定方程的橡胶材料本构模型参数识别

由于橡胶材料具有非线性和大变形特性,使得描述橡胶力学特性的本构模型参数的确定比较烦琐和困难。为了提高橡胶本构模型参数识别的准确性,基于超静定方程求解原理推导出一种新的识别方法。以某橡胶衬套为例,识别的参数应用于有限元分析,对比试验数据,结果表明该识别方法可以准确识别橡胶材料参数,并且精度相对于最小二乘法有了明显改善,显示了超静定识别法的有效性和可靠性。     

Finite element modeling of current stimulating response on the nerve axon

To investigate the characteristics of microphysiological electrical conduction on nerve cell axon, numerical simulation for the dynamic responses of neuronal axon to electrical stimulation is conducted on a three-dimensional (3D) finite element (FE) model. A 3D geometrical model for a segment of a hippocampal neuron axon is developed and assigned biophysical parameters, then the stimulation pulses of various amplitude and duration are imparted on the FE model of neuronal axon. The combination of Hodgkin-Huxley equations and Maxwell equations are performed to get action potential curve and 3D distribution of electric potential. The simulating results show that the resting potential of nerve axon is -65 mV and no action potential occurred under the stimuli of 2 ms duration and 0.01 A/m 2 current intensity, whereas the stimuli of (2 ms, 0.2 A/m 2), (20 ms, 0.01 A/m 2) and (20 ms, 0.2 A/m 2) induced action potentials, and the time arrived at the peak value are 0.012 s, 0.017 s and 0.012 s, respectively. The range of action potential is 100 mV and duration is 2 ms. The simulation results agree well with that of the experiment, which indicates that both the FE model of neuronal axon and the simulation analysis method are reliable and may contribute to further study on neural electrophysiology.     

Design optimization of magneto-rheological damper based on finite element parametric language

The mountain hazards like snow avalanches, landslides, rock falls, debris flows and so on all have strong power of destruction which seriously threaten human’s lives and belongings. Therefore, it is necessary to study more the development of these disasters in order to prevent them. Setting up obstacles is the primary measure to control the movement and deposition process of mountain hazards. For the study of the influence of this measure on the development of disasters, the numerical simulation calculation aiming at the flow and deposition process of the granulars flowing past different built obstacles was made by the theory of SH granular flow and the method of finite volume discretization based on Roe's Scheme. Futhermore, the influence of different settings of obstacles on the granular flow is discussed. The numerical results show that the settings of obstacles have a great effect on the process of granular flow. That is, the effective precaution against the disasters is no other than setting up the obstacles correctly. As the numerical simulation calculation can optimize the setting of obstacle, it provides a reasonable and economic reference scheme of disaster prevention and planning of mountain area.     

Elastic Buckling and Critical Load of an Pre- stressing Arch

According to the field theory of additional deformation on pre - stressed configuration , in the paper , the ordinary expression of the governing equation and variational equation of elastic buckling are brought forward . Under the theory system, through lowering dimensions, the governing equation and variational equation for the critical condition solution of elastic buckling of a plane arch are deduced, and the eigenvalues problem of the linear homogeneous differential equations corresponding to the equations are concluded While Abandoning the plane assumption and considering shearing deformation, the linear finite element method arithmetic of bending bar's cross section containing six degree of freedom is given. The process of derivation and calculational results show that, under this system info, the finite element equations of bending bar deduced are accurate and easy to be used to numeric calculations, and the conclusion achieved is more practical.     

Three Dimensional Eddy-current Field computation by Edge Finite Element Method of Mixed Elements

The edge finite element interpolation function of 1-forms for prism is derived.The electrical field intensity is used as the unknown variable.If the tangential component of electrical field intensity is an enforcing condition on the interface of eddy-current region and non-eddycurrent region,the tangential component of magnetic field intensity is a natural continuous condition.At last,the mixed elements of prism and hexahedron are used to compute a three dimensional eddy-current problem.     

A Study on the Response of Brick Masonry Structure under Dynamic and Static Shear-Compression Loading and Suggested Shearing Equation

Based on the studies and experiments of the responses of brick model under dynamic and static shear-compression loading, this paper brought forward an equation as a transition one for standard of masonry structure design in revision and the standard of aseismic building design. With its simple form, definite concept and similar calculation results in comparison with those obtained by equations from two standards, this equation has realized the correlation of the old and new ones. Besides it is rather save, so it has given a new advice for the revision of the standards.     

The Non-Linear Finite Element Analysis of Cable Structures Based on Four-Node Isoparametric Curved Element

In this paper,based on the geometric nonlinear behavior of cable structures,a finite element method with four-node isoparametric element was presented,the third order polynomial was used as displacement functions and initial curve of the element.By means of the principle of virtual work and the updated Lagrangian method,the authors derived the finite element equations,and solved them by Newton-Raphson method.The proposed model leads to high precision and can meet the engineering requirements.The model presented in this paper can be applied in the analysis of long-span tension structures,such as cable structures,cable domes and so forth.     

Modal Shape Analysis Methods of Machine Structure

The iterative Lanczos-reduce method for sensitivity analysis, the Guyan condensed method and the iterative synchronous vector methods in finite element modal are studied in order to solve the dynamic characteristic of complicated machine system. The primary reasons influenced on calculation precision are pointed out. The principia of master freedom for Guyan condensed method are put forward. The calculation precision of three methods is compared with the same standard finite element model. The iterative Lanczos-reduce method for sensitivity analysis is proved to be optimum analysis and the eigenvalue of arm in work device of hydraulic excavator in finite element model is solved with iterative Lanczos-reduce method for sensitivity analysis. The theory basis is provided for reasonable choose to solve the dynamic of complicated machine structure.     

Effects of Thorax Muscle Fiber on Body Surface Potentials

A numerical method calculating the forward problem of electrocardiogram (ECG) has been presented. By using this method, the torso part including anisotropic muscle layers has been discretized by means of the finite element method (FEM). The rest of torso have been divided into the surface elements by the boundary element method (BEM). The effects of the skelltel muscle layers over thorax wall, have been discussed in terms of a three-dimensional torso models.     

The Analysis of Abnormal State for Power Distribution Network Based on MATLAB

Along with the extensive use of power cable lines in power distribution system, the probability of fundamental frequency resonance has greatly increased. It may endanger the security operation and power supply reliability of the power system. On the basis of theoretical analysis and numerical calculations using MATLAB language combined with a real instance, it is indicated that the improper electrical parameter coordination caused by the abundance of power cable lines result in fundamental frequency current resonance, in addition, the countermeasures based on SVC to solve the problem are put forward. Further more, digital simulations have been carried out in detail based on the dynamic simulation model developed by the application of MATLAB/Simulink and its power system simulation software. Simulation results proves the effectiveness and practicability of the scheme presented.     

Simulation analysis on stimulation patterns of three-dimension electrical impedance tomography

To solve the lack of information and the error in forward and inverse problem caused by the approximation of actual three-dimensional field from two-dimensional field in electrical impedance tomography (EIT), this article presents a new EIT system. This system improves one-dimensional electrode array to two-dimensional array, and extends the solution domain to the entire three-dimensional field in order to increase the information quantity and eliminate the model error. Based on several current stimulation patterns in EIT which might yield the best sensitiveness for 3-D field, we calculate the forward problem by finite element method and the distinguishability and surface projection image method. The results show that the back electrode pattern has obvious advantages in detecting depth and precision, and has convenience in clinical application. The conclusion of this article can provide reference for the study of 3-D EIT.     

Coupled Thermo-Elastic Contact Finite Element Method and It''s Application in Gearing

A new finite element method is presented in this paper to analyze the thermo- elastic contact problem with heat input on meshing tooth faces. The computational model of meshing gear has been established and the frictional heat input on meshing tooth faces is estimated. In addition, 2-D finite element analysis program for this kind of problem is listed. The test and practical examples are given to show the validity and versatility of this method.     

Post-buckling and mode jumping analysis of angle-ply composite laminates

The governing partial differential equations (PDEs) are deduced from the asymptotically correct, geometrically nonlinear theory to research the bucking and mode jumping behavior of biaxially antisymmetric angle-ply composite laminates. The PDEs are transformed into a system of nonlinear ordinary differential equations (ODEs) by using the generalized Galerkin method. Then, the post-buckling behavior is analyzed by using the solution extension software. At last, the paths of buckling mode jumping and critical loads for the composite laminates under the in-plane boundary condition of straight edge are determined. An example of 4-layer composite laminates shows that the numerical results in the primary post-buckling region from the present method agree with the finite element analysis (FEA); while the FEA may lose its convergence when solution comes close to the secondary bifurcation point, the analytic method has the capability to explore deeply into the post-secondary bucking realm and capture the mode jumping phenomenon.