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.展开更多
This study presents a new method to solve the difficult problem of precise machining a non-cylinder pinhole of a piston using embedded giant magnetostrictive material (GMM) in the component. We propose the finite elem...This study presents a new method to solve the difficult problem of precise machining a non-cylinder pinhole of a piston using embedded giant magnetostrictive material (GMM) in the component. We propose the finite element model of GMM smart component in electric, magnetic, and mechanical fields by step computation to optimize the design of GMM smart com-ponent. The proposed model is implemented by using COMSOL multi-physics V3.2a. The effects of the smart component on the deformation and the system resonance frequencies are studied. The results calculated by the model are in excellent agreement (relative errors are below 10%) with the experimental values.展开更多
Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulati...Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.展开更多
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.展开更多
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.展开更多
Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradu...Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradually. In fact, the dimension fluctuation exists in the forming process all the while. A mathematical model of Gauss distribution for processing parameters was employed and a finite element numerical model of dieless drawing process with non-steady processing parameters was established. Dieless drawing processing of Ni-Ti alloy wire was conducted for verifying the proposed model. The results indicated that the non-steady processing parameters model had higher simulation accuracy of the wire diameter than that given by the steady parameters model. Furthermore, the model could also be used to analyze the fluctuation characteristics in the whole dieless drawing process.展开更多
Capacitive micromechanical ultrasonic transducers(CMUTs)have been widely studied because they can be used as substitutes for piezoelectric ultrasonic transducers in imaging applications.However,it is unclear whether a...Capacitive micromechanical ultrasonic transducers(CMUTs)have been widely studied because they can be used as substitutes for piezoelectric ultrasonic transducers in imaging applications.However,it is unclear whether and how CMUTs can be developed for sensors incorporating other functions.For instance,researchers have proposed the use of CMUTs for pressure sensing,but fundamental and practical application issues remain unsolved.This study explored ways in which a pressure sensor can be properly developed based on a CMUT prototype using a simulation approach.A three-dimensional finite element model of CMUTs was designed using the COMSOL Multiphysics software by combining the working principle of CMUTs with pressure sensing characteristics in which the resonance frequency of the CMUT cell shifts accordingly when it is subjected to an external pressure.Simultaneously,when subjected to pressure,the CMUT membrane deforms,thus the pressure can be reflected by the change in the capacitance.展开更多
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.展开更多
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 multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between ...The multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.展开更多
基金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.
基金supported by the National Natural Science Foundation of China (No. 50575205)the Hi-Tech Research and Development (863) Program of China (Nos. 2006AA04Z233 and 2007AA04Z101)+1 种基金the Doctoral Foundation of Ministry of Education of China (No. 20070335204)the Zhejiang Provincial Natural Science Foundation of China (No. Z1080537)
文摘This study presents a new method to solve the difficult problem of precise machining a non-cylinder pinhole of a piston using embedded giant magnetostrictive material (GMM) in the component. We propose the finite element model of GMM smart component in electric, magnetic, and mechanical fields by step computation to optimize the design of GMM smart com-ponent. The proposed model is implemented by using COMSOL multi-physics V3.2a. The effects of the smart component on the deformation and the system resonance frequencies are studied. The results calculated by the model are in excellent agreement (relative errors are below 10%) with the experimental values.
基金financially supported by the Riset Pengembangan dan Penerapan(RPP),Diponegoro University 2023 research scheme with contract number 609-18/UN7.D2/PP/VIII/2023.
文摘Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.
基金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.
基金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.
基金National Basic Research Program of China (2006CB605200)National Natural Science Foundation of China (50634010, 50674008)Program for New Century Excellent Talents in University (NCET-06-0083)
文摘Since processing parameters have always been assumed to be stable in the current finite element numerical simulation of dieless drawing process, the simulation results for the product dimension tend to stabilize gradually. In fact, the dimension fluctuation exists in the forming process all the while. A mathematical model of Gauss distribution for processing parameters was employed and a finite element numerical model of dieless drawing process with non-steady processing parameters was established. Dieless drawing processing of Ni-Ti alloy wire was conducted for verifying the proposed model. The results indicated that the non-steady processing parameters model had higher simulation accuracy of the wire diameter than that given by the steady parameters model. Furthermore, the model could also be used to analyze the fluctuation characteristics in the whole dieless drawing process.
文摘Capacitive micromechanical ultrasonic transducers(CMUTs)have been widely studied because they can be used as substitutes for piezoelectric ultrasonic transducers in imaging applications.However,it is unclear whether and how CMUTs can be developed for sensors incorporating other functions.For instance,researchers have proposed the use of CMUTs for pressure sensing,but fundamental and practical application issues remain unsolved.This study explored ways in which a pressure sensor can be properly developed based on a CMUT prototype using a simulation approach.A three-dimensional finite element model of CMUTs was designed using the COMSOL Multiphysics software by combining the working principle of CMUTs with pressure sensing characteristics in which the resonance frequency of the CMUT cell shifts accordingly when it is subjected to an external pressure.Simultaneously,when subjected to pressure,the CMUT membrane deforms,thus the pressure can be reflected by the change in the capacitance.
文摘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.
文摘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 multi-physics simulation of coupled fluid-structure interaction problems, with disjoint fluid and solid domains, requires one to choose a method for enforcing the fluid-structure coupling at the interface between solid and fluid. While it is common knowledge that the choice of coupling technique can be very problem dependent, there exists no satisfactory coupling comparison methodology that allows for conclusions to be drawn with respect to the comparison of computational cost and solution accuracy for a given scenario. In this work, we develop a computational framework where all aspects of the computation can be held constant, save for the method in which the coupled nature of the fluid-structure equations is enforced. To enable a fair comparison of coupling methods, all simulations presented in this work are implemented within a single numerical framework within the deal.ii [1] finite element library. We have chosen the two-dimensional benchmark test problem of Turek and Hron [2] as an example to examine the relative accuracy of the coupling methods studied;however, the comparison technique is equally applicable to more complex problems. We show that for the specific case considered herein the monolithic approach outperforms partitioned and quasi-direct methods;however, this result is problem dependent and we discuss computational and modeling aspects which may affect other comparison studies.
