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.展开更多
The structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress field...The structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress fields and basement tectonics are two fundamental factors to control the cover tectonic framework. This paper uses the finite-element method with a elastic-plastic pIan problem model to simulate the three periods of stress fields resulting from field geological study’ Based on these works, the formation and evolution of tectonic framework of Lu’ an mining area have been discussed.展开更多
A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The resul...A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.展开更多
The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-E...The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.展开更多
Spontaneous combustion of coal is a problem that affects the mining operation and generates environ-mental,economic,social and geotechnical impacts.This phenomenon has been divided into two pro-cesses:ignition and pro...Spontaneous combustion of coal is a problem that affects the mining operation and generates environ-mental,economic,social and geotechnical impacts.This phenomenon has been divided into two pro-cesses:ignition and propagation.Fire propagation develops in coal seams because of a set of factors such as direction and wind speed,fracturing and temperature.In this work,heat transfer and chemical kinetics are studied from conservation equations of energy and species,respectively,using the software COMSOL Multiphysics to simulate the propagation of fires in coal seams.Two possible scenarios were analyzed that usually occur in the walls of the coal seams,such as fire focus and fire complete screens.It was found that the propagation kinetics of the fire changes depending on the temperature,the fractur-ing of rock mass and the area of fire influence.For temperature values lower than 300℃,there is con-sumption around 250 cm^3/h,values around 700℃,the consumption is 1500 cm^3/h,and for fires of 1200℃ have values of 3000 cm^3/h.Depending on the speed of propagation can vary from 4 to 17cm/day,considering on the level and fracturing of the final wall of the open pit.展开更多
Soft robots have become important members of the robot community with many potential applications owing to their unique flexibility and security embedded at the material level.An increasing number of researchers are i...Soft robots have become important members of the robot community with many potential applications owing to their unique flexibility and security embedded at the material level.An increasing number of researchers are interested in their designing,manufacturing,modeling,and control.However,the dynamic simulation of soft robots is difficult owing to their infinite degrees of freedom and nonlinear characteristics that are associated with soft materials and flexible geometric structures.In this study,a novel multi-flexible body dynamic modeling and simulation technique is introduced for soft robots.Various actuators for soft robots are modeled in a virtual environment,including soft cable-driven,spring actuation,and pneumatic driving.A pneumatic driving simulation was demonstrated by the bending modules with different materials.A cable-driven soft robot arm prototype and a cylindrical soft module actuated by shape memory alley springs inspired by an octopus were manufactured and used to validate the simulation model,and the experimental results demonstrated adequate accuracy.The proposed technique can be widely applied for the modeling and dynamic simulation of other soft robots,including hybrid actuated robots and rigid-flexible coupling robots.This study also provides a fundamental framework for simulating soft mobile robots and soft manipulators in contact with the environment.展开更多
Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that...Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.展开更多
As hydraulic fracturing was forbidden in some countries due to possible negative environmental impacts and enhanced coalbed methane(ECBM)was restricted by in-situ conditions,microwave heating was proposed to enhance c...As hydraulic fracturing was forbidden in some countries due to possible negative environmental impacts and enhanced coalbed methane(ECBM)was restricted by in-situ conditions,microwave heating was proposed to enhance coalbed permeability.One of the mechanisms of improving coal permeability with microwave irradiation is that thermal expansion caused by microwave heating.To study the influence of microwave’s heating effect of coal samples,the simulations were conducted using a coupled electromagnetic,thermal and mechanical model in this paper.The temperature,Von-Mises stress and strain distribution of coal sample are recorded every 10 s.The permeability distribution is also obtained based on the relationship between strain and permeability from articles.It was found that volume average temperature,stress,strain and permeability increase almost linearly with time.The average permeability increased from 1.65 10×16 m^2 to 3.63 10×16 m^2 under 2.45 GHz and 500Wmicrowave radiation after 300 s.The significant increase proved microwave to be effective in coal seam permeability enhancement.展开更多
Central shrinkage crack is a common defect encountered in steel ingot casting. It is necessary to limit the degree of crack in case of further propagation in forging. A 234-t steel ingot was dissected to check the int...Central shrinkage crack is a common defect encountered in steel ingot casting. It is necessary to limit the degree of crack in case of further propagation in forging. A 234-t steel ingot was dissected to check the internal quality, and a central shrinkage crack band of 1,400 mm in height and 120 mm in width, was found at a distance of 450 mm under the riser bottom line. Then, thermo-mechanical simulation using an elasto-viscoplastic finite-element model was conducted to analyze the stress-strain evolution during ingot solidification. A new criterion considering mush mechanical property in the brittle temperature range as well as shrinkage porosity was used to identify the shrinkage crack potential, where the degree of shrinkage porosity is regarded as a probability factor using a modified sigmoid function. Different casting processes, such as pouring speed, mould preheating and riser insulation, were optimized with the simulation model. The results show that fast pouring, proper mould preheating and good riser insulation can alleviate shrinkage crack potential in the ingot center.展开更多
The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carr...The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.展开更多
Based on the Finite Element Theory of Rigid Plastic,relevant problems during plas-tic simulation on sheet metal forming were technologically studied and simplified;a simplified model of the blank holder during the dra...Based on the Finite Element Theory of Rigid Plastic,relevant problems during plas-tic simulation on sheet metal forming were technologically studied and simplified;a simplified model of the blank holder during the drawing process was established and the effects of related parameters on the forming processes were also studied.At the same time,a finite--element numerical simulation program SPID was developed.The distribution of strain and relationship of stress--stroke simulated were compared with experimented,the results are well coincided with each other.It is verified that the analytical program is reliable.展开更多
The corrosion behavior of B10 copper-nickel alloy welded joints in seawater pipeline system was analyzed under local turbulence induced by weld residual height.The corrosion behavior was evaluated by array electrode t...The corrosion behavior of B10 copper-nickel alloy welded joints in seawater pipeline system was analyzed under local turbulence induced by weld residual height.The corrosion behavior was evaluated by array electrode technology,mor-phology and elemental characterization,and COMSOL Multiphysics simulation.The results provide a theoretical basis for the corrosion and leakage of B10 alloy in seawater pipeline under the action of turbulence.The results show that residual height-induced turbulence exhibits a significant effect on the corrosion behavior in different areas of welded joints in B10 alloy.Turbulence can damage some surfaces,causing polarity deflection followed by acceleration of corrosion,or it is easier to form a protective film to slow down corrosion.COMSOL Multiphysics results show that the shear rate and turbulent kinetic energy increase linearly with the increase in residual height and velocity.The corrosion behavior of alloy surface is influenced by controlling the mass transfer rate and surface state.展开更多
The capacitively coupled radio frequency(CCRF)plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the ...The capacitively coupled radio frequency(CCRF)plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the behavior in plasma.In this paper,a glass vacuum chamber and a pair of plate electrodes were designed and fabricated,using 13.56 MHz radio frequency(RF)discharge technology to ionize the working gas of Ar.This discharge was mathematically described with equivalent circuit model.The discharge voltage and current of the plasma were measured atdifferent pressures and different powers.Based on the capacitively coupled homogeneous discharge model,the equivalent circuit and the analytical formula were established.The plasma density and temperature were calculated by using the equivalent impedance principle and energy balance equation.The experimental results show that when RF discharge power is 50–300 W and pressure is 25–250 Pa,the average electron temperature is about 1.7–2.1 e V and the average electron density is about 0.5?×10^17–3.6?×10^17m^-3.Agreement was found when the results were compared to those given by optical emission spectroscopy and COMSOL simulation.展开更多
Kraft lignin has the potential to replace traditional fossil resources for the preparation of high-value chemicals because it is rich in aromatic rings and active functional groups.An effective method for the pyrolysi...Kraft lignin has the potential to replace traditional fossil resources for the preparation of high-value chemicals because it is rich in aromatic rings and active functional groups.An effective method for the pyrolysis of kraft lignin into chemicals/fuels is microwave-assisted depolymerization.A simulation model is urgently needed to illustrate the coupling effect and mechanism of lignin conversion during the depolymerization process.In this study,COMSOL Multiphysics was used to simulate the microwave-assisted depolymerization process.The results showed that microwave power had a significant effect on the electric field and temperature distribution in the microwave cavity,while the reaction time had little effect on the electric field.The effect of the nitrogen flow rate on the electric field and temperature was negligible.The intensity of the electric field,heating rate of lignin,and final temperature of lignin depolymerization increased with increasing microwave power.展开更多
We investigated the deformation behaviors of Zr_65Cu_17.5Ni_10Al_7.5 in superplastic forming in silicon mould via numerical modeling and experiments. The data needed for the constitutive formulation were obtained from...We investigated the deformation behaviors of Zr_65Cu_17.5Ni_10Al_7.5 in superplastic forming in silicon mould via numerical modeling and experiments. The data needed for the constitutive formulation were obtained from compressive tests to establish a material library for finite-element simulation using a DEFORM 3D software. A constant speed forming process of a micro gear was modeled where the loading force, feature size and amount of deformation in the micro gear in silicon mould were analyzed in detail for the optimal requirements of micro gear forming and the protection of silicon mould. Guided by the modeling parameters, an amorphous metal micro gear was successfully obtained by our home-made superplastic forming system with the optimized parameters (temperature of 683 K, top speed of 0.003 mm/s until the load force reaching limiting value at 1960 N, and a gradually decelerating process for holding the force to the end). Our work gives a good example for optimization of superplastic forming and fabrication of BMGs in microparts.展开更多
A numerical model of 2.5D non-isothermal resin transfer molding simulation is developed for thin part based on the control volume/finite element method. The non-uniform temperature distribution and the heat generation...A numerical model of 2.5D non-isothermal resin transfer molding simulation is developed for thin part based on the control volume/finite element method. The non-uniform temperature distribution and the heat generation during the filling stage are modeled with the lumped temperature system and the species balance. Numerical algorithm of the simulation are studied. The molding simulation for a part is performed to show the effectiveness of simulating filling time, temperature distribution and curing degree.展开更多
Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues ...Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues for precise treatment within intricate regions of the human body.展开更多
This study investigates the breakdown voltage characteristics in sulfur hexafluoride(SF6)circuit breakers,employing a novel approach that integrates both experimental investigations and finite element simulations.Util...This study investigates the breakdown voltage characteristics in sulfur hexafluoride(SF6)circuit breakers,employing a novel approach that integrates both experimental investigations and finite element simulations.Utilizing a sphere-sphere electrode configuration,we meticulously measured the relationship between breakdown voltage and electrode gap distances ranging from 1 cm to 4.5 cm.Subsequent simulations,conducted using COMSOL Multiphysics,mirrored the experimental setup to validate the model’s accuracy through a comparison of the breakdown voltage-electrode gap distance curves.The simulation results not only aligned closely with the experimental data but also allowed the extraction of detailed electric field strength,electric potential contours,and electric current flow curves at the breakdown voltage for gap distances extending from 1 to 4.5 cm.Extending the analysis,the study explored the electric field and potential distribution at a constant voltage of 72.5 kV for gap distances between 1 to 10 cm,identifying the maximum electric field strength.A comprehensive comparison of five different electrode configurations(sphere-sphere,sphere-rod,sphere-plane,rod-plane,rod-rod)at 72.5 kV and a gap distance of 1.84 cm underscored the significant influence of electrode geometry on the breakdown process.Moreover,the research contrasts the breakdown voltage in SF6 with that in air,emphasizing SF6’s superior insulating properties.This investigation not only elucidates the intricate dynamics of electrical breakdown in SF6 circuit breakers but also contributes valuable insights into the optimal electrode configurations and the potential for alternative insulating gases,steering future advancements in high-voltage circuit breaker technology.展开更多
基金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 structural de formation of Lu’ an mining area is characterized by a remarkable feature of zoning along E-W direction, in the east limb of Qinshui basin, Shanxi Province, China. The re gional tectonic stress fields and basement tectonics are two fundamental factors to control the cover tectonic framework. This paper uses the finite-element method with a elastic-plastic pIan problem model to simulate the three periods of stress fields resulting from field geological study’ Based on these works, the formation and evolution of tectonic framework of Lu’ an mining area have been discussed.
文摘A 3-D finite-element numerical simulation model of temperature field for CIESC casting solidification process was developed with the aid of ANSYS software and a series of corresponding experiments were made. The results showed that the good agreement was obtained between the numerical simulation and the experiments. Based on the numerical simulation results, the characteristics of temperature distribution in the castings during CIESC solidification process were analyzed and summarized. According to the G/R-1/2 method and numerical simulation results, there is no any shrinkage defect in the CIESC casting and structure or casting is fine and compact.
基金financially supported by the National Natural Science Foundation of China(No.41330638)
文摘The technology used to enhance coalbed methane(CBM) recovery by injecting CO_2(CO_2-ECBM) with heat, combining heat injection with CO_2 injection, is still in its infancy; therefore, theoretical studies of this CO_2-ECBM technology should be perused. First, the coupling equations of the di usion–adsorption–seepage–heat transfer fields of gas are established. The displacement processes under di erent pressures and temperatures are simulated by COMSOL. Finally, the displacement effects, a comparison of the CO_2 storage capacity with the CH_4 output and the e ective influencing radius of CO_2 injection are analyzed and discussed. The results show that(1) the displacement pressure and temperature are two key factors influencing the CH_4 output and the CO_2 storage capacity, and the increase in the CO_2 storage capacity is more sensitive to temperature and pressure than the CH_4 output.(2) The gas flow direction is from the injection hole to the discharge hole during the displacement process, and the regions with high velocity are concentrated at the injection hole and the discharge hole.(3) A reduction in the CH_4 concentration and an increase in the CO_2 concentration are obvious during the displacement process.(4) The e ective influencing radius of injecting CO_2 with heat increases with the increase in time and pressure. The relationship between the e ective influencing radius and the injection time of CO_2 has a power exponential function, and there is a linear relationship between the functional coe cient and the injection pressure of CO_2. This numerical simulation study on enhancing CBM recovery by injecting CO_2 with heat can further promote the implementation of CO_2-ECBM project in deep coal seams.
基金CIMEX Mineral Institute of the National University of Colombia for the trust and support provided for the development of this researchthe Cerrejón company for financing this project
文摘Spontaneous combustion of coal is a problem that affects the mining operation and generates environ-mental,economic,social and geotechnical impacts.This phenomenon has been divided into two pro-cesses:ignition and propagation.Fire propagation develops in coal seams because of a set of factors such as direction and wind speed,fracturing and temperature.In this work,heat transfer and chemical kinetics are studied from conservation equations of energy and species,respectively,using the software COMSOL Multiphysics to simulate the propagation of fires in coal seams.Two possible scenarios were analyzed that usually occur in the walls of the coal seams,such as fire focus and fire complete screens.It was found that the propagation kinetics of the fire changes depending on the temperature,the fractur-ing of rock mass and the area of fire influence.For temperature values lower than 300℃,there is con-sumption around 250 cm^3/h,values around 700℃,the consumption is 1500 cm^3/h,and for fires of 1200℃ have values of 3000 cm^3/h.Depending on the speed of propagation can vary from 4 to 17cm/day,considering on the level and fracturing of the final wall of the open pit.
基金Supported by the National Natural Science Foundation of China(Grant Nos.51822502 and 91948202)the National Key Research and Development Program of China(No.2019YFB1309500)the“111 Project”(Grant No.B07018).
文摘Soft robots have become important members of the robot community with many potential applications owing to their unique flexibility and security embedded at the material level.An increasing number of researchers are interested in their designing,manufacturing,modeling,and control.However,the dynamic simulation of soft robots is difficult owing to their infinite degrees of freedom and nonlinear characteristics that are associated with soft materials and flexible geometric structures.In this study,a novel multi-flexible body dynamic modeling and simulation technique is introduced for soft robots.Various actuators for soft robots are modeled in a virtual environment,including soft cable-driven,spring actuation,and pneumatic driving.A pneumatic driving simulation was demonstrated by the bending modules with different materials.A cable-driven soft robot arm prototype and a cylindrical soft module actuated by shape memory alley springs inspired by an octopus were manufactured and used to validate the simulation model,and the experimental results demonstrated adequate accuracy.The proposed technique can be widely applied for the modeling and dynamic simulation of other soft robots,including hybrid actuated robots and rigid-flexible coupling robots.This study also provides a fundamental framework for simulating soft mobile robots and soft manipulators in contact with the environment.
基金Project(2017ZX05013002-002)supported by Major National Science and Technology Projects of ChinaProject(RIPED-2016-JS-276)supported by Petro-China Research Institute of Petroleum Exploration and Development
文摘Crustal stresses play an important role in both exploration and development in the oil and gas industry.However,it is difficult to simulate crustal stress distributions accurately,because of the incompatibilities that exist among different software.Here,a series of algorithms is developed and integrated in the Petrel2ANSYS to carry out two-way conversions between the 3D attribute models that employ corner-point grids used in Petrel and the 3D finite-element grids used in ANSYS.Furthermore,a modified method of simulating stress characteristics and analyzing stress fields using the finite-element method and multiple finely resolved 3D models is proposed.Compared to the traditional finite-element simulation-based approach,which involves describing the heterogeneous within a rock body or sedimentary facies in detail and simulating the stress distribution,the single grid cell-based approach focuses on a greater degree on combining the rock mechanics described by 3D corner-point grid models with the finely resolved material characteristics of 3D finite-element models.Different models that use structured and unstructured grids are verified in Petrel2ANSYS to assess the feasibility.In addition,with minor modifications,platforms based on the present algorithms can be extended to other models to convert corner-point grids to the finite-element grids constructed by other software.
文摘As hydraulic fracturing was forbidden in some countries due to possible negative environmental impacts and enhanced coalbed methane(ECBM)was restricted by in-situ conditions,microwave heating was proposed to enhance coalbed permeability.One of the mechanisms of improving coal permeability with microwave irradiation is that thermal expansion caused by microwave heating.To study the influence of microwave’s heating effect of coal samples,the simulations were conducted using a coupled electromagnetic,thermal and mechanical model in this paper.The temperature,Von-Mises stress and strain distribution of coal sample are recorded every 10 s.The permeability distribution is also obtained based on the relationship between strain and permeability from articles.It was found that volume average temperature,stress,strain and permeability increase almost linearly with time.The average permeability increased from 1.65 10×16 m^2 to 3.63 10×16 m^2 under 2.45 GHz and 500Wmicrowave radiation after 300 s.The significant increase proved microwave to be effective in coal seam permeability enhancement.
基金supported by the NSFC-Liaoning Joint Fund(U1508215)the project to strengthen industrial development at the grass-roots level of MIIT China(TC160A310/21)
文摘Central shrinkage crack is a common defect encountered in steel ingot casting. It is necessary to limit the degree of crack in case of further propagation in forging. A 234-t steel ingot was dissected to check the internal quality, and a central shrinkage crack band of 1,400 mm in height and 120 mm in width, was found at a distance of 450 mm under the riser bottom line. Then, thermo-mechanical simulation using an elasto-viscoplastic finite-element model was conducted to analyze the stress-strain evolution during ingot solidification. A new criterion considering mush mechanical property in the brittle temperature range as well as shrinkage porosity was used to identify the shrinkage crack potential, where the degree of shrinkage porosity is regarded as a probability factor using a modified sigmoid function. Different casting processes, such as pouring speed, mould preheating and riser insulation, were optimized with the simulation model. The results show that fast pouring, proper mould preheating and good riser insulation can alleviate shrinkage crack potential in the ingot center.
基金financially supported by the National Key Research and Development Plan(No.2018YFB0605601)the National Natural Science Foundation of China(No.41972168)。
文摘The three-dimensional(3 D) structures of pores directly affect the CH4 flow.Therefore,it is very important to analyze the3 D spatial structure of pores and to simulate the CH4 flow with the connected pores as the carrier.The result shows that the equivalent radius of pores and throats are 1-16 μm and 1.03-8.9 μm,respectively,and the throat length is 3.28-231.25 μm.The coordination number of pores concentrates around three,and the intersection point between the connectivity function and the X-axis is 3-4 μm,which indicate the macro-pores have good connectivity.During the single-channel flow,the pressure decreases along the direction of CH4 flow,and the flow velocity of CH4 decreases from the pore center to the wall.Under the dual-channel and the multi-channel flows,the pressure also decreases along the CH4 flow direction,while the velocity increases.The mean flow pressure gradually decreases with the increase of the distance from the inlet slice.The change of mean flow pressure is relatively stable in the direction horizontal to the bedding plane,while it is relatively large in the direction perpendicular to the bedding plane.The mean flow velocity in the direction horizontal to the bedding plane(Y-axis) is the largest,followed by that in the direction horizontal to the bedding plane(X-axis),and the mean flow velocity in the direction perpendicular to the bedding plane is the smallest.
文摘Based on the Finite Element Theory of Rigid Plastic,relevant problems during plas-tic simulation on sheet metal forming were technologically studied and simplified;a simplified model of the blank holder during the drawing process was established and the effects of related parameters on the forming processes were also studied.At the same time,a finite--element numerical simulation program SPID was developed.The distribution of strain and relationship of stress--stroke simulated were compared with experimented,the results are well coincided with each other.It is verified that the analytical program is reliable.
基金support from the National Natural Science Foundation of China(Grant No.42176209)the Natural Science Foundation of Shandong Province(Grant No.ZR2021MD064)the Fundamental Research Funds for the Central Universities(Grant No.19CX05001A).
文摘The corrosion behavior of B10 copper-nickel alloy welded joints in seawater pipeline system was analyzed under local turbulence induced by weld residual height.The corrosion behavior was evaluated by array electrode technology,mor-phology and elemental characterization,and COMSOL Multiphysics simulation.The results provide a theoretical basis for the corrosion and leakage of B10 alloy in seawater pipeline under the action of turbulence.The results show that residual height-induced turbulence exhibits a significant effect on the corrosion behavior in different areas of welded joints in B10 alloy.Turbulence can damage some surfaces,causing polarity deflection followed by acceleration of corrosion,or it is easier to form a protective film to slow down corrosion.COMSOL Multiphysics results show that the shear rate and turbulent kinetic energy increase linearly with the increase in residual height and velocity.The corrosion behavior of alloy surface is influenced by controlling the mass transfer rate and surface state.
基金supported by National Natural Science Foundation of China(Grant No.61378037)the Fundamental Research Funds for the Central Universities(Nos.2013B33614,2017B15214)+1 种基金the Research Funds of Innovation and Entrepreneurship Education Reform for Chinese Universities(No.16CCJG01Z004)the Changzhou Science and Technology Program(No.CJ20160027)
文摘The capacitively coupled radio frequency(CCRF)plasma has been widely used in various fields.In some cases,it requires us to estimate the range of key plasma parameters simpler and quicker in order to understand the behavior in plasma.In this paper,a glass vacuum chamber and a pair of plate electrodes were designed and fabricated,using 13.56 MHz radio frequency(RF)discharge technology to ionize the working gas of Ar.This discharge was mathematically described with equivalent circuit model.The discharge voltage and current of the plasma were measured atdifferent pressures and different powers.Based on the capacitively coupled homogeneous discharge model,the equivalent circuit and the analytical formula were established.The plasma density and temperature were calculated by using the equivalent impedance principle and energy balance equation.The experimental results show that when RF discharge power is 50–300 W and pressure is 25–250 Pa,the average electron temperature is about 1.7–2.1 e V and the average electron density is about 0.5?×10^17–3.6?×10^17m^-3.Agreement was found when the results were compared to those given by optical emission spectroscopy and COMSOL simulation.
基金This work was supported by the Foundation of Key Laboratory of Pulp and Paper Science and Technology of the Ministry of Education of China(No.KF201917)the National Natural Science Foundation of China(31800497).
文摘Kraft lignin has the potential to replace traditional fossil resources for the preparation of high-value chemicals because it is rich in aromatic rings and active functional groups.An effective method for the pyrolysis of kraft lignin into chemicals/fuels is microwave-assisted depolymerization.A simulation model is urgently needed to illustrate the coupling effect and mechanism of lignin conversion during the depolymerization process.In this study,COMSOL Multiphysics was used to simulate the microwave-assisted depolymerization process.The results showed that microwave power had a significant effect on the electric field and temperature distribution in the microwave cavity,while the reaction time had little effect on the electric field.The effect of the nitrogen flow rate on the electric field and temperature was negligible.The intensity of the electric field,heating rate of lignin,and final temperature of lignin depolymerization increased with increasing microwave power.
基金Funded by the National Natural Science Foundation of China(Nos.51222508,51175211)
文摘We investigated the deformation behaviors of Zr_65Cu_17.5Ni_10Al_7.5 in superplastic forming in silicon mould via numerical modeling and experiments. The data needed for the constitutive formulation were obtained from compressive tests to establish a material library for finite-element simulation using a DEFORM 3D software. A constant speed forming process of a micro gear was modeled where the loading force, feature size and amount of deformation in the micro gear in silicon mould were analyzed in detail for the optimal requirements of micro gear forming and the protection of silicon mould. Guided by the modeling parameters, an amorphous metal micro gear was successfully obtained by our home-made superplastic forming system with the optimized parameters (temperature of 683 K, top speed of 0.003 mm/s until the load force reaching limiting value at 1960 N, and a gradually decelerating process for holding the force to the end). Our work gives a good example for optimization of superplastic forming and fabrication of BMGs in microparts.
文摘A numerical model of 2.5D non-isothermal resin transfer molding simulation is developed for thin part based on the control volume/finite element method. The non-uniform temperature distribution and the heat generation during the filling stage are modeled with the lumped temperature system and the species balance. Numerical algorithm of the simulation are studied. The molding simulation for a part is performed to show the effectiveness of simulating filling time, temperature distribution and curing degree.
基金supported by NSFC(62273019,52072015,12332019,U20A20390)the 111 Project(B13003)。
文摘Dear Editor,This letter presents a biocompatible cross-shaped magnetic soft robot and investigates its deformation mode control strategy through COMSOL modeling and simulation.Magnetic soft robots offer novel avenues for precise treatment within intricate regions of the human body.
基金Ningbo Science and Technology Plan Project(Grant No.2023Z043)。
文摘This study investigates the breakdown voltage characteristics in sulfur hexafluoride(SF6)circuit breakers,employing a novel approach that integrates both experimental investigations and finite element simulations.Utilizing a sphere-sphere electrode configuration,we meticulously measured the relationship between breakdown voltage and electrode gap distances ranging from 1 cm to 4.5 cm.Subsequent simulations,conducted using COMSOL Multiphysics,mirrored the experimental setup to validate the model’s accuracy through a comparison of the breakdown voltage-electrode gap distance curves.The simulation results not only aligned closely with the experimental data but also allowed the extraction of detailed electric field strength,electric potential contours,and electric current flow curves at the breakdown voltage for gap distances extending from 1 to 4.5 cm.Extending the analysis,the study explored the electric field and potential distribution at a constant voltage of 72.5 kV for gap distances between 1 to 10 cm,identifying the maximum electric field strength.A comprehensive comparison of five different electrode configurations(sphere-sphere,sphere-rod,sphere-plane,rod-plane,rod-rod)at 72.5 kV and a gap distance of 1.84 cm underscored the significant influence of electrode geometry on the breakdown process.Moreover,the research contrasts the breakdown voltage in SF6 with that in air,emphasizing SF6’s superior insulating properties.This investigation not only elucidates the intricate dynamics of electrical breakdown in SF6 circuit breakers but also contributes valuable insights into the optimal electrode configurations and the potential for alternative insulating gases,steering future advancements in high-voltage circuit breaker technology.