An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) ...An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.展开更多
Focusing on the vibration of the roadbed and ground induced by high-speed train load, a three dimensional finite element model which includes the roadbed and horizontal layered site is established to study how the sit...Focusing on the vibration of the roadbed and ground induced by high-speed train load, a three dimensional finite element model which includes the roadbed and horizontal layered site is established to study how the site conditions, the load moving speed and the depth of the soil element influence the soil element stress response. Based on a track-subsoil analytical model in which the rail is simulated as an Euler-Bernoulli beam resting on Winkler foundation in the vertical plane, the reaction force between the sleeper and roadbed excited by a single axle is presented, and then that is exerted on relevant elements to simulate the moving load. The dynamic response in the roadbed and subsoil excited by a single axle moving load is computed based on the parallel computing platform of the ABAQUS finite element software, and the stress time-history, stress path and curves of the principal stress axes rotation of the soil element under the track are presented. The results show that： the soil element stress path is an apple-shaped curve in the horizontal shear stress τd versus the stress difference （σsh - σch ）/2 coordinate system; the principal stress axes rotate 180° for the soil element under the load moving line during the load running, and the stress state changes from the pure shear to triaxial shear and then back to the pure shear again. The element dynamic stress increases as the moving load speed increases, which increases sharply when the load speed approaches the Rayleigh wave velocity of soil layer; the site conditions and the soil element depth affect the soil element stress path significantly.展开更多
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 an...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.展开更多
This paper studies the damage-viscoelastic behavior of composite solid propellants of solid rocket motors(SRM).Based on viscoelastic theories and strain equivalent hypothesis in damage mechanics,a three-dimensional(3-...This paper studies the damage-viscoelastic behavior of composite solid propellants of solid rocket motors(SRM).Based on viscoelastic theories and strain equivalent hypothesis in damage mechanics,a three-dimensional(3-D)nonlinear viscoelastic constitutive model incorporating with damage is developed.The resulting viscoelastic constitutive equations are numerically discretized by integration algorithm,and a stress-updating method is presented by solving nonlinear equations according to the Newton-Raphson method.A material subroutine of stress-updating is made up and embedded into commercial code of Abaqus.The material subroutine is validated through typical examples.Our results indicate that the finite element results are in good agreement with the analytical ones and have high accuracy,and the suggested method and designed subroutine are efficient and can be further applied to damage-coupling structural analysis of practical SRM grain.展开更多
In this paper the variation of interface stress intensity factors is considered under arbitrary material combinations for several problems such as central internal interface cracks, periodic interface cracks, and edge...In this paper the variation of interface stress intensity factors is considered under arbitrary material combinations for several problems such as central internal interface cracks, periodic interface cracks, and edge interface cracks. The stress intensity factors for all these problems have been determined in terms of Dundurs＇ parameters a, b by the application of FEM in the authors＇ previous papers. In this paper, the variations of the stress intensity factors for all these problems are discussed under arbitrary a, b, then the maximum and minimum values of dimensionless stress intensity factors FI,max, FII,max, FI,min, FII,max are indicated on the map of a, b. As an example, it is found that FI,max always appears when （a, b）= （0.2,0.3）, and FI,min always appears when （a, b）= （1,0）, for the central internal interface cracks independent of crack length.展开更多
基金Project(2007CB714202) supported by the National Key Basic Research Program of ChinaProject(SLDRCE10-B-07) supported by theMinistry of Science and Technology of China
文摘An element coupling model (ECM) method was proposed to simulate the global behavior and local damage of a structure.In order to reflect the local damage and improve the computational efficiency,three-dimensional (3D) solid elements and one-dimensional (1D) beam element were coupled by the multi-point constraint equations.A reduced scale 1?8 model test was simulated by the ECM and a full three dimensional model (3DM) contrastively.The results show that the global behavior and local damages of ECM agree well with the test and 3DM.It is indicated that the proposed method can be used in the structural nonlinear analysis accurately and efficiently.
文摘Focusing on the vibration of the roadbed and ground induced by high-speed train load, a three dimensional finite element model which includes the roadbed and horizontal layered site is established to study how the site conditions, the load moving speed and the depth of the soil element influence the soil element stress response. Based on a track-subsoil analytical model in which the rail is simulated as an Euler-Bernoulli beam resting on Winkler foundation in the vertical plane, the reaction force between the sleeper and roadbed excited by a single axle is presented, and then that is exerted on relevant elements to simulate the moving load. The dynamic response in the roadbed and subsoil excited by a single axle moving load is computed based on the parallel computing platform of the ABAQUS finite element software, and the stress time-history, stress path and curves of the principal stress axes rotation of the soil element under the track are presented. The results show that： the soil element stress path is an apple-shaped curve in the horizontal shear stress τd versus the stress difference （σsh - σch ）/2 coordinate system; the principal stress axes rotate 180° for the soil element under the load moving line during the load running, and the stress state changes from the pure shear to triaxial shear and then back to the pure shear again. The element dynamic stress increases as the moving load speed increases, which increases sharply when the load speed approaches the Rayleigh wave velocity of soil layer; the site conditions and the soil element depth affect the soil element stress path significantly.
文摘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.
基金supported by the National Natural Science Foundation of China(Grant No.11132012)
文摘This paper studies the damage-viscoelastic behavior of composite solid propellants of solid rocket motors(SRM).Based on viscoelastic theories and strain equivalent hypothesis in damage mechanics,a three-dimensional(3-D)nonlinear viscoelastic constitutive model incorporating with damage is developed.The resulting viscoelastic constitutive equations are numerically discretized by integration algorithm,and a stress-updating method is presented by solving nonlinear equations according to the Newton-Raphson method.A material subroutine of stress-updating is made up and embedded into commercial code of Abaqus.The material subroutine is validated through typical examples.Our results indicate that the finite element results are in good agreement with the analytical ones and have high accuracy,and the suggested method and designed subroutine are efficient and can be further applied to damage-coupling structural analysis of practical SRM grain.
基金supported by the Japanese Government (Monbukagukusho) Scholarshipthe National Basic Research Program of China (Grant No.2011CB013702)
文摘In this paper the variation of interface stress intensity factors is considered under arbitrary material combinations for several problems such as central internal interface cracks, periodic interface cracks, and edge interface cracks. The stress intensity factors for all these problems have been determined in terms of Dundurs＇ parameters a, b by the application of FEM in the authors＇ previous papers. In this paper, the variations of the stress intensity factors for all these problems are discussed under arbitrary a, b, then the maximum and minimum values of dimensionless stress intensity factors FI,max, FII,max, FI,min, FII,max are indicated on the map of a, b. As an example, it is found that FI,max always appears when （a, b）= （0.2,0.3）, and FI,min always appears when （a, b）= （1,0）, for the central internal interface cracks independent of crack length.
