Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the...Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.展开更多
Projection welding is a variation of electric resistance welding with the dynamic changes of the flow paths for heat and electrical properties with changing temperature caused by the large plastic deformation collapse...Projection welding is a variation of electric resistance welding with the dynamic changes of the flow paths for heat and electrical properties with changing temperature caused by the large plastic deformation collapse of projection. As the joint type between the auto door hinge and the inner plate, projection welding may bring welding distortions and would affect the assembly quality of auto body. A comprehensive electric-thermal-mechanical numerical simulation was performed to quantitatively simulate the processes of projection welding by using a coupled finite element method. The mechanism of projection collapse and the formation process of nugget were discussed and good conclusions have been achieved comparing with the test results.展开更多
Purpose Charged particle cascade simulations,in conjunction with finite element thermal and mechanical calculations,are essential for addressing engineering challenges associated with the design of accelerator beam du...Purpose Charged particle cascade simulations,in conjunction with finite element thermal and mechanical calculations,are essential for addressing engineering challenges associated with the design of accelerator beam dumps.This study aims to investigate the integration of the FLUKA Monte Carlo program and the COMSOL finite element program as indispensable tools in the development of advanced models for spatial energy deposition distributions,temperature assessments,and stress analysis within the DALS beam dump.Methods The paper delves into aspects of model development,data transfer,and practical applications,focusing on the successful coupling of FLUKA Monte Carlo simulations with COMSOL finite element analyses.Tens of millions of elements were created and utilized to estimate spatial energy deposition distributions,assess temperatures,and analyze stresses within the critical absorber of the DALS beam dump.The study comprehensively analyzes the process from the impact of high-energy electron beam on the beam dump to the spatial distribution of energy deposition,providing input for subsequent thermal and structural analyses.Results The successful coupling of FLUKA and COMSOL enabled the calculation of spatial temperature distributions and structural analyses of the absorber within the DALS beam dump.The utilization of tens of millions of defined bins ensured seamless data transfer from the particle cascade simulation to finite element analysis,guaranteeing high resolution and accuracy in the calculations.The results provide valuable insights into the thermal and mechanical behavior of the beam dump absorber,a critical safety component in accelerator systems.Conclusions The study demonstrates that the integration of FLUKA Monte Carlo simulations with COMSOL finite element analyses is a dependable and efficient tool for addressing real-world engineering challenges,particularly those related to the design of beam dumps in accelerator systems with charged particle beams.The advanced analytical approach provides crucial information for the optimal design and safety assessment of accelerator components.展开更多
Finite element simulations were conducted to study the mechanism of spark plasma sintering. The spark plasma sintering of SiC ceramics was simulated by the Marc software based on the load current curve and temperature...Finite element simulations were conducted to study the mechanism of spark plasma sintering. The spark plasma sintering of SiC ceramics was simulated by the Marc software based on the load current curve and temperature-time curve deserved by SPS experiment. The concept of equivalent radiation coefficient was presented and applied during the simulation. The temperature distribution regularity of SiC ceramics sintered by SPS technology was got by thermal-electrical coupled finite element simulation. The experimental results show that by thermal-electrical coupled finite element analysis, the temperature rising and distribution regularity of nonconductive material can be preferable forecasted in the sintering process of SPS. In the initial stage of the heat preservation, the temperature of the central part of the sample has achieved sintering temperature, but now, the temperature of the sample is not uniform. The temperature for each part of the die is also quite different and the sample temperature in the center is higher than that in the edge. In the end of heat preservation, the central temperature of the sample is 50 ℃higher than the required sintering temperature, and the temperature gap for each part of the die decreases gradually.展开更多
A coupled water and heat transport mode is established based on the Richards equation to study water flow and heat transport in soil during freezing process. Both the finite difference and finite element method are us...A coupled water and heat transport mode is established based on the Richards equation to study water flow and heat transport in soil during freezing process. Both the finite difference and finite element method are used in the discretization, respectively. Two different computer programs are written and used to simulate an indoor unidirectional frozen test. The freezing depth, freezing rate and temperature variation are compared among lab tests, finite difference calculation simulation and finite element calculation simulation. Result shows that: the finite difference method has a better performance in freezing depth simulation while the finite element method has a better performance in numerical stability in one-dimensional freezing simulation.展开更多
The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field...The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.展开更多
For establishing the refined numerical simulation model for coupled vibration between vehicle and bridge, the refined three-dimensional vehicle model is setup by multi-body system dynamics method, and finite element m...For establishing the refined numerical simulation model for coupled vibration between vehicle and bridge, the refined three-dimensional vehicle model is setup by multi-body system dynamics method, and finite element method of dynamic model is adopted to model the bridge. Taking Yujiang River Bridge on Nanning-Guangzhou railway line in China as study background, the?refined numerical simulation model of whole vehicle and whole bridge system for coupled vibration analysis is set up. The dynamic analysis model of the cable-stayed bridge is established by finite element method, and the natural vibration properties of the bridge are analyzed. The German ICE Electric Multiple Unit (EMU) train refined three-dimensional space vehicle model is set up by multi-system dynamics software SIMPACK, and the multiple non-linear properties are considered. The space vibration responses are calculated by co-simulation based on multi-body system dynamics and finite element method when the ICE EMU train passes the long span cable-stayed bridge at different speeds. In order to test if the bridge has the sufficient lateral or vertical rigidity and the operation stability is fine. The calculation results show: The operation safety can be guaranteed, and comfort?index is “excellent”. The bridge has sufficient rigidity, and vibration is in good condition.展开更多
The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-el...The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-element modeling, we evaluated the relationship between plate kinematics and present-day deformation of Luzon Island and adjacent sea areas. The concept of coupling rate was applied to define the boundary velocities along the subduction zones. The distribution of velocity fields calculated in our models was compared with the velocity field revealed by recent geodetic (GPS) observations. The best model was obtained that accounts for the observed velocity field within the limits of acceptable mechanical parameters and reasonable boundary conditions. Sensitivity of the selection of parameters and boundary conditions were evaluated. The model is sensitive to the direction of convergence between the South China Sea and the Philippine Sea plates, and to different coupling rates in the Manila trench, Philippine trench and eastern Luzon trough. We suggest that a change of ±15° of the direction of motion of the Philippine Sea plate can induce important changes in the distribution of the computed displacement trajectories, and the movement of the Philippine Sea plate toward azimuth 330° best explains the velocity pattern observed in Luzon Island. In addition, through sensitivity analysis we conclude that the coupling rate in the Manila trench is much smaller compared with the rates in the eastern Luzon trough and the Philippine trench. This indicates that a significant part of momentum of the Philippine Sea plate motion has been absorbed by the Manila trench; whereas, a part of the momentum has been transmitted into Luzon Island through the eastern Luzon trough and the Philippine trench.展开更多
In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations ...In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations which are derived from the Biot’s consolidation theory.Fluid pore pressure,solid displacement,and temperature are chosen as initial variables in these equations,and the finite element method in combination with the interface element is used for spatial discretization of continuous and discontinuities(fault)parts of the medium to solve the equations.The main purpose of this study is providing precise formulations,applicability,and ability of the triple-node zero-thickness interface element in THM modeling of faults.It should be noted that the system of equations is solved using a computer code written in Matlab program.In order to verify the developed method,simulations of index problems such as Mandel’s problem,and coupled modeling of a faulted porous medium and a faulted aquifer are presented.The modeling results obtained from the developed method show a very good agreement with those by other modeling methods,which indicates its accuracy.展开更多
基金the National High Technical Reasearch and Development Programme of China (No. 2003AA327140) the National Natural Science Foundation of China (No. 50374081).
文摘Two full 3D steady mathematical models are developed by finite element method (FEM) to calcalate coupled physics fields. the electro-magnetic model is built and solved first and so is the fluid motion model with the acquired electromagnetic force as source body forces in Navier-Stokes equations. Effects caused by the ferromagnetic shell, busbar system around, and open boundary problem as well as inside induced current were considered in terms of the magnetic field. Furthermore, a new modeling method is found to set up solid models and then mesh them entirely with so-called structuralized grids, namely hex-mesh. Examples of 75kA prebaked cell with two kinds of busbar arrangements are presented. Results agree with those disclosed in the literature and confirm that the coupled simulation is valid. It is also concluded that the usage of these models facilitates the consistent analysis of the electric field to magnetic field and then flow motion to the greater extent, local distributions of current density and magnetic flux density are very much dependent on the cell structure, the steel shell is a shield to reduce the magnetic field and flow pattern is two dimensional in the main body of the metal pad.
基金The National Natural Science Foundation of China(No.50575140)
文摘Projection welding is a variation of electric resistance welding with the dynamic changes of the flow paths for heat and electrical properties with changing temperature caused by the large plastic deformation collapse of projection. As the joint type between the auto door hinge and the inner plate, projection welding may bring welding distortions and would affect the assembly quality of auto body. A comprehensive electric-thermal-mechanical numerical simulation was performed to quantitatively simulate the processes of projection welding by using a coupled finite element method. The mechanism of projection collapse and the formation process of nugget were discussed and good conclusions have been achieved comparing with the test results.
文摘Purpose Charged particle cascade simulations,in conjunction with finite element thermal and mechanical calculations,are essential for addressing engineering challenges associated with the design of accelerator beam dumps.This study aims to investigate the integration of the FLUKA Monte Carlo program and the COMSOL finite element program as indispensable tools in the development of advanced models for spatial energy deposition distributions,temperature assessments,and stress analysis within the DALS beam dump.Methods The paper delves into aspects of model development,data transfer,and practical applications,focusing on the successful coupling of FLUKA Monte Carlo simulations with COMSOL finite element analyses.Tens of millions of elements were created and utilized to estimate spatial energy deposition distributions,assess temperatures,and analyze stresses within the critical absorber of the DALS beam dump.The study comprehensively analyzes the process from the impact of high-energy electron beam on the beam dump to the spatial distribution of energy deposition,providing input for subsequent thermal and structural analyses.Results The successful coupling of FLUKA and COMSOL enabled the calculation of spatial temperature distributions and structural analyses of the absorber within the DALS beam dump.The utilization of tens of millions of defined bins ensured seamless data transfer from the particle cascade simulation to finite element analysis,guaranteeing high resolution and accuracy in the calculations.The results provide valuable insights into the thermal and mechanical behavior of the beam dump absorber,a critical safety component in accelerator systems.Conclusions The study demonstrates that the integration of FLUKA Monte Carlo simulations with COMSOL finite element analyses is a dependable and efficient tool for addressing real-world engineering challenges,particularly those related to the design of beam dumps in accelerator systems with charged particle beams.The advanced analytical approach provides crucial information for the optimal design and safety assessment of accelerator components.
基金Funded by the Natural Science Foundation of Hebei Province, China (No.E2012203086)
文摘Finite element simulations were conducted to study the mechanism of spark plasma sintering. The spark plasma sintering of SiC ceramics was simulated by the Marc software based on the load current curve and temperature-time curve deserved by SPS experiment. The concept of equivalent radiation coefficient was presented and applied during the simulation. The temperature distribution regularity of SiC ceramics sintered by SPS technology was got by thermal-electrical coupled finite element simulation. The experimental results show that by thermal-electrical coupled finite element analysis, the temperature rising and distribution regularity of nonconductive material can be preferable forecasted in the sintering process of SPS. In the initial stage of the heat preservation, the temperature of the central part of the sample has achieved sintering temperature, but now, the temperature of the sample is not uniform. The temperature for each part of the die is also quite different and the sample temperature in the center is higher than that in the edge. In the end of heat preservation, the central temperature of the sample is 50 ℃higher than the required sintering temperature, and the temperature gap for each part of the die decreases gradually.
基金the support and motivation provided by National 973 Project of China (No. 2012CB026104)National Natural Science Foundation of China (No. 41171064)
文摘A coupled water and heat transport mode is established based on the Richards equation to study water flow and heat transport in soil during freezing process. Both the finite difference and finite element method are used in the discretization, respectively. Two different computer programs are written and used to simulate an indoor unidirectional frozen test. The freezing depth, freezing rate and temperature variation are compared among lab tests, finite difference calculation simulation and finite element calculation simulation. Result shows that: the finite difference method has a better performance in freezing depth simulation while the finite element method has a better performance in numerical stability in one-dimensional freezing simulation.
文摘The electric field intensity (EFI) is important characteristic quantity for evaluating the internal insulation state of cable joints. Based on finite element method, this paper proposes two EFI research methods, field-circuit coupling method and equivalent circuit method. The average EFI of the inner surface of the outer semi-conducting shield can be calculated from the current in the measuring circuit. The relative error between these two methods is about 15%, which roughly proves the consistency of the two methods. Further practical application research enables online monitoring of cable joints.
文摘For establishing the refined numerical simulation model for coupled vibration between vehicle and bridge, the refined three-dimensional vehicle model is setup by multi-body system dynamics method, and finite element method of dynamic model is adopted to model the bridge. Taking Yujiang River Bridge on Nanning-Guangzhou railway line in China as study background, the?refined numerical simulation model of whole vehicle and whole bridge system for coupled vibration analysis is set up. The dynamic analysis model of the cable-stayed bridge is established by finite element method, and the natural vibration properties of the bridge are analyzed. The German ICE Electric Multiple Unit (EMU) train refined three-dimensional space vehicle model is set up by multi-system dynamics software SIMPACK, and the multiple non-linear properties are considered. The space vibration responses are calculated by co-simulation based on multi-body system dynamics and finite element method when the ICE EMU train passes the long span cable-stayed bridge at different speeds. In order to test if the bridge has the sufficient lateral or vertical rigidity and the operation stability is fine. The calculation results show: The operation safety can be guaranteed, and comfort?index is “excellent”. The bridge has sufficient rigidity, and vibration is in good condition.
基金supported by the National Scientific and Tech-nological Support Project of China (No. 2006BAB19B02)the Guangdong Natural Science Foundation (No. 07004206)the National Natural Science Foundation of China (No. 40476026)
文摘The Luzon Island is a volcanic arc sandwiched by the eastward subducting South China Sea and the northwestward subducting Philippine Sea plate. Through experiments of plane-stress, elastic, and 2-dimensional finite-element modeling, we evaluated the relationship between plate kinematics and present-day deformation of Luzon Island and adjacent sea areas. The concept of coupling rate was applied to define the boundary velocities along the subduction zones. The distribution of velocity fields calculated in our models was compared with the velocity field revealed by recent geodetic (GPS) observations. The best model was obtained that accounts for the observed velocity field within the limits of acceptable mechanical parameters and reasonable boundary conditions. Sensitivity of the selection of parameters and boundary conditions were evaluated. The model is sensitive to the direction of convergence between the South China Sea and the Philippine Sea plates, and to different coupling rates in the Manila trench, Philippine trench and eastern Luzon trough. We suggest that a change of ±15° of the direction of motion of the Philippine Sea plate can induce important changes in the distribution of the computed displacement trajectories, and the movement of the Philippine Sea plate toward azimuth 330° best explains the velocity pattern observed in Luzon Island. In addition, through sensitivity analysis we conclude that the coupling rate in the Manila trench is much smaller compared with the rates in the eastern Luzon trough and the Philippine trench. This indicates that a significant part of momentum of the Philippine Sea plate motion has been absorbed by the Manila trench; whereas, a part of the momentum has been transmitted into Luzon Island through the eastern Luzon trough and the Philippine trench.
文摘In this paper,a coupled thermo-hydro-mechanical(THM)simulation in a faulted deformable porous medium is presented.This model involves solving the mass conservation,linear momentum balance,and energy balance equations which are derived from the Biot’s consolidation theory.Fluid pore pressure,solid displacement,and temperature are chosen as initial variables in these equations,and the finite element method in combination with the interface element is used for spatial discretization of continuous and discontinuities(fault)parts of the medium to solve the equations.The main purpose of this study is providing precise formulations,applicability,and ability of the triple-node zero-thickness interface element in THM modeling of faults.It should be noted that the system of equations is solved using a computer code written in Matlab program.In order to verify the developed method,simulations of index problems such as Mandel’s problem,and coupled modeling of a faulted porous medium and a faulted aquifer are presented.The modeling results obtained from the developed method show a very good agreement with those by other modeling methods,which indicates its accuracy.