Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthqu...Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthquakes.However,no previous studies have assessed the mechanisms underlying seismic failure in rock slopes.In this study,large-scale shaking table tests and numerical simulations were conducted to delineate the seismic failure mechanism in terms of acceleration,displacement,and earth pressure responses combined with shattering failure phenomena.The results reveal that acceleration response mutations usually occur within weak interlayers owing to their inferior performance,and these mutations may transform into potential sliding surfaces,thereby intensifying the nonlinear seismic response characteristics.Cumulative permanent displacements at the internal corners of the berms can induce quasi-rigid displacements at the external corners,leading to greater permanent displacements at the internal corners.Therefore,the internal corners are identified as the most susceptible parts of the slope.In addition,the concept of baseline offset was utilized to explain the mechanism of earth pressure responses,and the result indicates that residual earth pressures at the internal corners play a dominant role in causing deformation or shattering damage.Four evolutionary deformation phases characterize the processes of seismic responses and shattering failure of the bedding parallel stepped rock slope,i.e.the formation of tensile cracks at the internal corners of the berm,expansion of tensile cracks and bedding surface dislocation,development of vertical tensile cracks at the rear edge,and rock mass slipping leading to slope instability.Overall,this study provides a scientific basis for the seismic design of engineering slopes and offers valuable insights for further studies on preventing seismic disasters in bedding parallel stepped rock slopes.展开更多
The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to ...The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.展开更多
During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configura...During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.展开更多
This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooli...This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooling systems.Through a combination of numerical simulation and experimental validation,the technical advantages and mechanisms of immersion cooling are systematically explored.Simulation results indicate that under a 3C fast-charging condition(inlet temperature 20℃,flow rate 36 L/min),the immersion cooling structure 3demonstrates a triple enhancement in thermal performance compared to the cold plate structure 1:a 13.06%reduction in peak temperature,a 31.67%decrease in overall maximum temperature difference,and a 47.62%decrease in single-cell temperature deviation,while also reducing flow resistance by 33.61%.Furthermore,based on the immersion cooling model,a small battery module comprising seven cylindrical cells was designed for thermal runaway testing via nail penetration.The results show that the peak temperature of the triggered cell was limited to 437.6℃,with a controllable temperature rise gradient of only 3.35℃/s and a rapid cooling rate of 0.6℃/s.The maximum temperature rise of adjacent cells was just 64.8℃,effectively inhibiting thermal propagation.Post-test disassembly revealed that the non-triggered cells retained>99.2%of their original voltage and>99%structural integrity,confirming the module’s ability to achieve“localized failure with global stability.”展开更多
A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocit...A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocity field during the filling process and the temperature field during the solidification process of the alloy melt under different process parameters were obtained.Based on the simulation results,a Zr-based amorphous alloy micro-gear was prepared via casting.The results indicate that increasing the suction casting temperature enhances the fluidity of alloy melt but induces unstable flow rate during filling,which is detrimental to complete filling.Zr-based amorphous micro-gears with a module of 0.6 mm,a tooth top diameter of 8 mm,and 10 teeth were prepared through the suction casting.X-ray diffraction and differential scanning calorimetry analyses confirm that the fabricated micro-gear exhibits characteristic amorphous structural features,demonstrating well-defined geometrical contours and satisfactory forming completeness.展开更多
Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing ...Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing a greater role in the performance of components. This paper investigates the effect of heat treatment on residual stresses induced in AISI 1025, manufactured using LENS. Finite element model was developed and simulated to analyze residual stress development. AISI 1025 samples suitable for tool and die applications in Fused Deposition Modelling (FDM) filament production, were fabricated using Laser Engineered Net Shaping (LENS) process, followed by heat treatment where annealing and quenching processes were done. The material’s microstructure, residual stress and hardness of heat-treated samples under investigation, were compared against the as-built samples. The results indicated that after annealing, tensile residual stresses were reduced by 93%, resulting in a reduced crack growth rate, compared to the as-built sample, although the hardness was reduced significantly by 25%. On the other hand, high tensile residual stresses of 425 ± 14 MPa were recorded after quenching process with an improvement of hardness by 21%.展开更多
In this study,we simulated the thermal behavior of the mud-brick walls of a Nubian vault.We used EnergyPlus software for the simulation.The results obtained showed that the indoor temperature varies from 25.5℃ to 26....In this study,we simulated the thermal behavior of the mud-brick walls of a Nubian vault.We used EnergyPlus software for the simulation.The results obtained showed that the indoor temperature varies from 25.5℃ to 26.5℃ for the period of January 2018.It varies from 33.2℃ to 33.6℃ with an average value of 33.1℃ for the month of April 2018.For the period of July 2018,it varies from 30.3℃ to 32.2℃ with an average value of 31.2℃..Relative humidity for the period of July ranged from 62.3%to 73.5%,with an average value of 67.9%.The simulation enabled us to compare simulated and measured temperature and humidity values.We found that the level of thermal comfort in the Nubian vault is acceptable in both cool and hot periods.In view of these results,we can say that the Nubian vault is an architecture suited to our climate.The technical concept of the Nubian vault is adapted to the climatic conditions and traditional know-how of the Sahel.We also found that the use of raw earth,a locally available material,and the Nubian vault architectural process,contribute to thermal comfort and a reappropriation of local and adapted know-how.展开更多
Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized...Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.展开更多
We proposed a new technique route of directional solidification for the manufacture of super slab.A 7-t laboratory-scale thick slab was casted and characterised for trial.To further understand the process,the evolutio...We proposed a new technique route of directional solidification for the manufacture of super slab.A 7-t laboratory-scale thick slab was casted and characterised for trial.To further understand the process,the evolution of the multiple physical fields during the directional solidification was simulated and verified.Similar to the convectional ingot casting,a negative segregated cone of equiaxed grains was formed at the bottom,and a seriously positive segregated region was formed beneath the top surface of the slab.Specific measures on the lateral walls,base plate,and free surface were strongly recommended to ensure that the slab is relatively directionally casted.A water-cooling copper base plate accelerates the solidification rate and the columnar growth along the vertical direction.It inhibits the sedimentation of equiaxed grains and enlarges the columnar zone.Based on the simulation analysis,it can be concluded that the directional solidification technique route is promising to manufacture super slab with lower segregation level,and less porosities and inclusions.展开更多
Inclusions in nickel-based superalloys significantly influence their mechanical properties which limit the application and development.A two-dimensional axisymmetric model coupling electromagnetic flow,heat transfer,a...Inclusions in nickel-based superalloys significantly influence their mechanical properties which limit the application and development.A two-dimensional axisymmetric model coupling electromagnetic flow,heat transfer,and inclusions transport was developed using the finite element method.The effects of current intensity and frequency on the transport behavior of inclusions and removal rate during vacuum induction melting were investigated using this model.To verify the accuracy of the mathematical model,experiments were conducted using the vacuum induction furnace model VIF200.A comparison of the experimental results with the simulation results reveals an excellent agreement.Four eddies exist in the central section of the molten pool,with any two eddies flowing in opposite directions.The drag force exerted by the eddies causes the inclusions converge toward the four corners of the cut surface.Due to buoyancy,an increase in the particle size of inclusions makes it easier for them to be adsorbed by the free surface,a phenomenon that is particularly noticeable for inclusions with a particle size of 400μm.An increase in current intensity changes the adsorption interface of certain 400-μm inclusions from the free surface to the sidewall.Increasing both current intensity and frequency enhances the crucible removal rate of inclusions,with current intensity exerting a more significant effect,increasing the removal rate by approximately 1% for every 150 A.Upon completion of melting,the ingot is typically obtained by pouring or bottom pouring,and inclusions adsorbed to the free surface are difficult to separate.Therefore,higher current intensity and frequency should be employed during melting to enhance the crucible removal rate.展开更多
In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components ...In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.展开更多
The growing complexity of electromagnetic(EM)interference has driven significant demand for next-generation absorbers that combine lightweight,flexibility,and good electromagnetic attenuation capability.The low-dimens...The growing complexity of electromagnetic(EM)interference has driven significant demand for next-generation absorbers that combine lightweight,flexibility,and good electromagnetic attenuation capability.The low-dimensional ternary Co_(3)ZnC/Co/CNT composites with hollow structures have been synthesized through in-situ polymerization and high-temperature carbonization.The unique integration of low-dimensional nanostructures and multicomponent heterointerfaces confers exceptional EM absorption properties,achieving a reflection loss of−70.0 dB and significantly reducing radar cross section(RCS)scattering signals.It is particularly meaningful that the numerical simulation of Co_(3)ZnC/Co/CNT metama-terial reveals ultrawideband absorption performance,achieving 10.7 GHz(7.3-18.0 GHz)at a thickness of 4.5 mm and extending to 15 GHz(3.0-18.0 GHz)with a 10.5 mm.Moreover,the Co_(3)ZnC/Co/CNT composites retain meritorious EM absorption properties after flexible film formation,broadening their usability and application scope.These investigations will provide seminal insights encompassing theoretical validation,experimental synthesis,and practical application for the next generation of absorbers.展开更多
The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbule...The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbulence models struggle to make accurate predictions for subsonic and supersonic flows in nozzles.In this study,we explored a novel model,the algebraic stress model k-kL-ARSM+J,to enhance the accuracy of turbulence numerical simulations.This new model was used to conduct numerical simulations of the design and off-design performance of a 3D supersonic asymmetric truncated nozzle designed in our laboratory,with the aim of providing a realistic pattern of changes.The research indicates that,compared to linear eddy viscosity turbulence models such as k-kL and shear stress transport(SST),the k-kL-ARSM+J algebraic stress model shows better accuracy in predicting the performance of supersonic nozzles.Its predictions were identical to the experimental values,enabling precise calculations of the nozzle.The performance trends of the nozzle are as follows:as the inlet Mach number increases,both thrust and pitching moment increase,but the rate of increase slows down.Lift peaks near the design Mach number and then rapidly decreases.With increasing inlet pressure,the nozzle thrust,lift,and pitching moment all show linear growth.As the flight altitude rises,the internal flow field within the nozzle remains relatively consistent due to the same supersonic nozzle inlet flow conditions.However,external to the nozzle,the change in external flow pressure results in the nozzle exit transitioning from over-expanded to under-expanded,leading to a shear layer behind the nozzle that initially converges towards the nozzle center and then diverges.展开更多
The local time-stepping(LTS)algorithm is an adaptive method that adjusts the time step by selecting suitable intervals for different regions based on the spatial scale of each cell and water depth and flow velocity be...The local time-stepping(LTS)algorithm is an adaptive method that adjusts the time step by selecting suitable intervals for different regions based on the spatial scale of each cell and water depth and flow velocity between cells.The method can be optimized by calculating the maximum power of two of the global time step increments in the domain,allowing the optimal time step to be approached throughout the grid.To verify the acceleration and accuracy of LTS in storm surge simulations,we developed a model to simulate astronomical storm surges along the southern coast of China.This model employs the shallow water equations as governing equations,numerical discretization using the finite volume method,and fluxes calculated by the Roe solver.By comparing the simulation results of the traditional global time-stepping algorithm with those of the LTS algorithm,we find that the latter fit the measured data better.Taking the calculation results of Typhoon Sally in 1996 as an example,we show that compared with the traditional global time-stepping algorithm,the LTS algorithm reduces computation time by 2.05 h and increases computation efficiency by 2.64 times while maintaining good accuracy.展开更多
As the integration of electronic components in high-performance servers increases,heat generation significantly impacts performance and raises failure rates.Therefore,heat dissipation has become a critical concern in ...As the integration of electronic components in high-performance servers increases,heat generation significantly impacts performance and raises failure rates.Therefore,heat dissipation has become a critical concern in electronic circuit design.This study uses numerical simulations to investigate the heat dissipation characteristics of electronic components in air-cooled servers.By adjusting airflow speed,heat sink configurations,and the arrangement of straight-fin heat sinks,we optimize heat dissipation performance and analyze the mechanisms at different airflow speeds.The results show that,at the same airflow speed,the temperature of the heat sink is lower than that of the electronic components,creating a temperature gradient that enhances heat transfer.Compared to a front-to-back arrangement of two straight-fin heat sinks,placing the heat sinks parallel to each other results in a lower maximum component temperature and better temperature uniformity.Heat sinks with fins significantly improve heat dissipation.The heat sink with semicylindrical fins on the rib surface provides the best cooling performance.Moreover,compared to natural convection,the maximum temperature of the electronic components decreases by 56.17%and 61%when the incoming flow velocity is 6 m/s with two parallel flat ribbed heat sinks and front-to-back arrangement,respectively.展开更多
As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-s...As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.展开更多
Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the ...Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the low-permeability characteristics with an average permeability of 5.5 mD and moderate heterogeneity,a 3-D geological model of heterogeneous Class I hydrate reservoirs with three-phase transition layers is established by Kriging interpolation and stochastic modeling method,and a numerical simulation model is used to describe the depressurization production performance of the reservoir.With the development of depressurization,a specific range of complete decomposition zones appear both in the hydrate and transition layers.The entire decomposition zone of the whole reservoir tends to outward and upward diffusion.There is apparent methane escape in the three-phase transition layer.Due to the improvement of local permeability caused by the phase transition of hydrate dissociation,some methane accumulation occurs at the bottom of the hydrate layer,forming a local methane enrichment zone.The methane migration trends in reservoirs are mainly characterized by movement toward production wells and hydrate layers under the influence of gravity.However,due to the permeability limitation of hydrate reservoirs,many fluids have not been effectively produced and remain in the reservoir.Therefore,to improve the effective pressure drop of the reservoir,the perforation method and pressure reduction method were optimized by analyzing the influencing factors based on the gas production rate.The comparative study demonstrates that perforating through the free gas layer combined with one-time depressurization can enhance the effective depressurization and improve production performance.The gas production rate from perforating through the free gas layer can be twice as high as that from perforating through the transition layer.This study can provide theoretical support for the utilization of marine energy.展开更多
The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the res...The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the results are reliable.Nozzle fluid is influenced by pipe structure,causing pressure and density to drop as speed increases.Differences in pressure and velocity between the jet and surrounding gas lead to jet velocity attenuation,flow expansion,deflection,and eddy currents.The optimal top blowing pressure is 0.6 MPa,and the center velocity and width of the jet are 345 m/s and 0.124 m,respectively,at 20De(De is the nozzle exit diameter).It achieves a maximum jet velocity of 456 m/s.The optimal nozzle Mach number is 1.75,with a maximum jet velocity of 451 m/s.At 20D_(e),the jet center velocity is 338 m/s,with a width of 0.12 m.展开更多
This paper presents a comprehensive experimental and numerical investigation of radiant floor heating(RFH)systems integrated with phase changematerial(PCM)-based thermal energy storage(TES).The study compares two unde...This paper presents a comprehensive experimental and numerical investigation of radiant floor heating(RFH)systems integrated with phase changematerial(PCM)-based thermal energy storage(TES).The study compares two underfloor pipe configurations:double serpentine and spiral.It also looks at how well a paraffin wax PCM system works with compact heat exchanger-type TES units during winter in Iraq.Key performance indicators including discharge temperature,heat transfer rate,liquid fraction evolution,and temperature uniformity were assessed through in situ experimental measurements and ANSYS fluent simulations.Results demonstrate that the spiral design provides slightly more uniform temperature distribution on the tile surface at an inlet water temperature of 55℃,with an average difference of approximately 0.5%,the serpentine layout exhibits higher slab temperature distribution by about 0.66%.Notably,the serpentine configuration shows superior thermal homogeneity and heat distribution,with a 15.05%increase in heat gain at a 55℃ inlet temperature compared to the spiral design.The performance gap between the two layouts narrows as the inlet temperature increases from 50℃ in 5℃ increments by approximately 4.1%,3.7%,and 1.7%,respectively.Higher inlet temperatures also improve PCM discharging and charging rates,improving energy storage utilization.The findings provide significant design guidelines for sustainable heating systems for cold climates.展开更多
The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope ...The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.展开更多
基金supported by the National Natural Science Foundation of China (Grant No.52108361)the Sichuan Science and Technology Program of China (Grant No.2023YFS0436)the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection Independent Research Project (Grant No.SKLGP2022Z015).
文摘Bedding parallel stepped rock slopes exist widely in nature and are used in slope engineering.They are characterized by complex topography and geological structure and are vulnerable to shattering under strong earthquakes.However,no previous studies have assessed the mechanisms underlying seismic failure in rock slopes.In this study,large-scale shaking table tests and numerical simulations were conducted to delineate the seismic failure mechanism in terms of acceleration,displacement,and earth pressure responses combined with shattering failure phenomena.The results reveal that acceleration response mutations usually occur within weak interlayers owing to their inferior performance,and these mutations may transform into potential sliding surfaces,thereby intensifying the nonlinear seismic response characteristics.Cumulative permanent displacements at the internal corners of the berms can induce quasi-rigid displacements at the external corners,leading to greater permanent displacements at the internal corners.Therefore,the internal corners are identified as the most susceptible parts of the slope.In addition,the concept of baseline offset was utilized to explain the mechanism of earth pressure responses,and the result indicates that residual earth pressures at the internal corners play a dominant role in causing deformation or shattering damage.Four evolutionary deformation phases characterize the processes of seismic responses and shattering failure of the bedding parallel stepped rock slope,i.e.the formation of tensile cracks at the internal corners of the berm,expansion of tensile cracks and bedding surface dislocation,development of vertical tensile cracks at the rear edge,and rock mass slipping leading to slope instability.Overall,this study provides a scientific basis for the seismic design of engineering slopes and offers valuable insights for further studies on preventing seismic disasters in bedding parallel stepped rock slopes.
基金supported by the Guangdong Major Project of Basic and Applied Basic Research(Grant No.2021B0301030001)the National Key Research and Development Program of China(Grant No.2021YFB3802300)the Foundation of National Key Laboratory of Shock Wave and Detonation Physics(Grant No.JCKYS2022212004)。
文摘The graded density impactor(GDI)dynamic loading technique is crucial for acquiring the dynamic physical property parameters of materials used in weapons.The accuracy and timeliness of GDI structural design are key to achieving controllable stress-strain rate loading.In this study,we have,for the first time,combined one-dimensional fluid computational software with machine learning methods.We first elucidated the mechanisms by which GDI structures control stress and strain rates.Subsequently,we constructed a machine learning model to create a structure-property response surface.The results show that altering the loading velocity and interlayer thickness has a pronounced regulatory effect on stress and strain rates.In contrast,the impedance distribution index and target thickness have less significant effects on stress regulation,although there is a matching relationship between target thickness and interlayer thickness.Compared with traditional design methods,the machine learning approach offers a10^(4)—10^(5)times increase in efficiency and the potential to achieve a global optimum,holding promise for guiding the design of GDI.
文摘During the propagation of high-power lasers within internal channels,the laser beam heats the propagation medium,causing the thermal blooming effect that degrades the beam quality at the output.The intricate configuration of the optical path within the internal channel necessitates complex and time-consuming efforts to assess the impact of thermal blooming effect on the optical path.To meet the engineering need for rapid evaluation of thermal blooming effect in optical paths,this study proposed a rapid simulation method for the thermal blooming effect in internal optical paths based on the finite element method.This method discretized the fluid region into infinitesimal elements and employed finite element method for flow field analysis.A simplified analytical model of the flow field region in complex internal channels was established,and regions with similar thermal blooming effect were divided within this model.Based on the calculated optical path differences within these regions,numerical simulations of phase distortion caused by thermal blooming were conducted.The calculated result were compared with those obtained using the existing methods.The findings reveal that for complex optical paths,the discrepancy between the two approaches is less than 3.6%,with similar phase distortion patterns observed.For L-type units,this method and the existing methods identify the same primary factors influencing aberrations and exhibit consistent trends in their variation.This method was used to analyze the impact of thermal blooming effect in a straight channel under different gravity directions.The results show that phase distortion varies with changes in the direction of gravity,and the magnitude of the phase difference is strongly correlated with the component of gravity perpendicular to the optical axis.Compared to the existing methods,this approach offers greater flexibility,obviates the need for complex custom analysis programming.The analytical results of this method enable a rapid assessment of the thermal blooming effect in optical paths within the internal channel.This is especially useful during the engineering design.These results also provide crucial references for developing strategies to suppress thermal blooming effect.
文摘This study focuses on the thermal management of 4680-type cylindrical lithium-ion battery packs utilizing NCM811 chemistry.It establishes coupled multi-physics models for both immersion and serpentine cold plate cooling systems.Through a combination of numerical simulation and experimental validation,the technical advantages and mechanisms of immersion cooling are systematically explored.Simulation results indicate that under a 3C fast-charging condition(inlet temperature 20℃,flow rate 36 L/min),the immersion cooling structure 3demonstrates a triple enhancement in thermal performance compared to the cold plate structure 1:a 13.06%reduction in peak temperature,a 31.67%decrease in overall maximum temperature difference,and a 47.62%decrease in single-cell temperature deviation,while also reducing flow resistance by 33.61%.Furthermore,based on the immersion cooling model,a small battery module comprising seven cylindrical cells was designed for thermal runaway testing via nail penetration.The results show that the peak temperature of the triggered cell was limited to 437.6℃,with a controllable temperature rise gradient of only 3.35℃/s and a rapid cooling rate of 0.6℃/s.The maximum temperature rise of adjacent cells was just 64.8℃,effectively inhibiting thermal propagation.Post-test disassembly revealed that the non-triggered cells retained>99.2%of their original voltage and>99%structural integrity,confirming the module’s ability to achieve“localized failure with global stability.”
基金National Natural Science Foundation of China(51971103)Key Research and Development Program in Gansu Province(20YF8GA052)。
文摘A suction casting experiment was conducted on Zr_(55)Cu_(30)Al_(10)Ni_(5)(at%)amorphous alloy.Using ProCAST software,numerical simulations were performed to analyze the filling and solidification processes.The velocity field during the filling process and the temperature field during the solidification process of the alloy melt under different process parameters were obtained.Based on the simulation results,a Zr-based amorphous alloy micro-gear was prepared via casting.The results indicate that increasing the suction casting temperature enhances the fluidity of alloy melt but induces unstable flow rate during filling,which is detrimental to complete filling.Zr-based amorphous micro-gears with a module of 0.6 mm,a tooth top diameter of 8 mm,and 10 teeth were prepared through the suction casting.X-ray diffraction and differential scanning calorimetry analyses confirm that the fabricated micro-gear exhibits characteristic amorphous structural features,demonstrating well-defined geometrical contours and satisfactory forming completeness.
文摘Heat treatment processes, such as annealing and quenching, are crucial in determining residual stress evolution, microstructural changes and mechanical properties of metallic materials, with residual stresses playing a greater role in the performance of components. This paper investigates the effect of heat treatment on residual stresses induced in AISI 1025, manufactured using LENS. Finite element model was developed and simulated to analyze residual stress development. AISI 1025 samples suitable for tool and die applications in Fused Deposition Modelling (FDM) filament production, were fabricated using Laser Engineered Net Shaping (LENS) process, followed by heat treatment where annealing and quenching processes were done. The material’s microstructure, residual stress and hardness of heat-treated samples under investigation, were compared against the as-built samples. The results indicated that after annealing, tensile residual stresses were reduced by 93%, resulting in a reduced crack growth rate, compared to the as-built sample, although the hardness was reduced significantly by 25%. On the other hand, high tensile residual stresses of 425 ± 14 MPa were recorded after quenching process with an improvement of hardness by 21%.
文摘In this study,we simulated the thermal behavior of the mud-brick walls of a Nubian vault.We used EnergyPlus software for the simulation.The results obtained showed that the indoor temperature varies from 25.5℃ to 26.5℃ for the period of January 2018.It varies from 33.2℃ to 33.6℃ with an average value of 33.1℃ for the month of April 2018.For the period of July 2018,it varies from 30.3℃ to 32.2℃ with an average value of 31.2℃..Relative humidity for the period of July ranged from 62.3%to 73.5%,with an average value of 67.9%.The simulation enabled us to compare simulated and measured temperature and humidity values.We found that the level of thermal comfort in the Nubian vault is acceptable in both cool and hot periods.In view of these results,we can say that the Nubian vault is an architecture suited to our climate.The technical concept of the Nubian vault is adapted to the climatic conditions and traditional know-how of the Sahel.We also found that the use of raw earth,a locally available material,and the Nubian vault architectural process,contribute to thermal comfort and a reappropriation of local and adapted know-how.
基金financially supported by the National Natural Science Foundation of China(No.51904026)the Fundamental Research Funds for the Central Universities(No.06500108)。
文摘Hydrogen displays the potential to partially replace pulverized coal injection(PCI)in the blast furnace,and it can reduce CO_(2)emissions.In this paper,a three-dimensional mathematical model of hydrogen and pulverized coal co-injection in blast furnace tuyere was established through numerical simulation,and the effect of hydrogen injection and oxygen enrichment interaction on pulverized coal combustion and raceway smelting was investigated.The simulation results indicate that when the coal injection rate decreased from 36 to 30t/h and the hydrogen injection increased from 0 to 3600 m^(3)/h,the CO_(2)emissions decreased from 1860 to 1551 kg/t,which represents a16.6%reduction,and the pulverized coal burnout decreased from 70.1%to 63.7%.The heat released from hydrogen combustion can not only promote the volatilization of pulverized coal but also affect the combustion reaction between volatilization and oxygen,which resulted in a decrease in the temperature at the end of the raceway.Co-injection of hydrogen with PCI increased the wall temperature near the upper half part of the raceway and at the outlet of the tuyere,which required a high cooling efficiency to extend the service life of the blast furnace.The increase in oxygen level compensated for the decreased average temperature in the raceway due to hydrogen injection.The increase in the oxygen content by 3%while maintaining constant hydrogen and PCI injection rates increased the burnout and average raceway temperature by 4.2%and 43 K,respectively.The mole fraction of CO and H_(2) production increased by 0.04 and 0.02,respectively.Burnout can be improved through optimization of the particle size distribution of pulverized coal.
基金the National Natural Science Foundation of China(No.52074182)Joint Funds of the National Natural Science Foundation of China(No.U23A20612).
文摘We proposed a new technique route of directional solidification for the manufacture of super slab.A 7-t laboratory-scale thick slab was casted and characterised for trial.To further understand the process,the evolution of the multiple physical fields during the directional solidification was simulated and verified.Similar to the convectional ingot casting,a negative segregated cone of equiaxed grains was formed at the bottom,and a seriously positive segregated region was formed beneath the top surface of the slab.Specific measures on the lateral walls,base plate,and free surface were strongly recommended to ensure that the slab is relatively directionally casted.A water-cooling copper base plate accelerates the solidification rate and the columnar growth along the vertical direction.It inhibits the sedimentation of equiaxed grains and enlarges the columnar zone.Based on the simulation analysis,it can be concluded that the directional solidification technique route is promising to manufacture super slab with lower segregation level,and less porosities and inclusions.
基金financially supported by the National Natural Science Foundation of China(Nos.52422408 and 52171031)Excellent Youth Fund of Liaoning Natural Science Foundation(No.2023JH3/10200001)Liaoning Xingliao Talents-Top-notch Young Talents Project(No.XLYC2203064).
文摘Inclusions in nickel-based superalloys significantly influence their mechanical properties which limit the application and development.A two-dimensional axisymmetric model coupling electromagnetic flow,heat transfer,and inclusions transport was developed using the finite element method.The effects of current intensity and frequency on the transport behavior of inclusions and removal rate during vacuum induction melting were investigated using this model.To verify the accuracy of the mathematical model,experiments were conducted using the vacuum induction furnace model VIF200.A comparison of the experimental results with the simulation results reveals an excellent agreement.Four eddies exist in the central section of the molten pool,with any two eddies flowing in opposite directions.The drag force exerted by the eddies causes the inclusions converge toward the four corners of the cut surface.Due to buoyancy,an increase in the particle size of inclusions makes it easier for them to be adsorbed by the free surface,a phenomenon that is particularly noticeable for inclusions with a particle size of 400μm.An increase in current intensity changes the adsorption interface of certain 400-μm inclusions from the free surface to the sidewall.Increasing both current intensity and frequency enhances the crucible removal rate of inclusions,with current intensity exerting a more significant effect,increasing the removal rate by approximately 1% for every 150 A.Upon completion of melting,the ingot is typically obtained by pouring or bottom pouring,and inclusions adsorbed to the free surface are difficult to separate.Therefore,higher current intensity and frequency should be employed during melting to enhance the crucible removal rate.
基金co-supported by the National Natural Science Foundation of China(Nos.52105411,52105400and 52305420)the China Postdoctoral Science Foundation(No.2023M742830)Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(No.CX2023008).
文摘In this study,a new linear friction welding(LFW)process,embedded LFW process,was put forward,which was mainly applied to combination manufacturing of long or overlong loadcarrying titanium alloy structural components in aircraft.The interfacial plastic flow behavior and bonding mechanism of this process were investigated by a developed coupling EulerianLagrangian numerical model using software ABAQUS and a novel thermo-physical simulation method with designed embedded hot compression specimen.In addition,the formation mechanism and control method of welding defects caused by uneven plastic flow were discussed.The results reveal that the plastic flow along oscillating direction of this process is even and sufficient.In the direction perpendicular to oscillation,thermo-plastic metals mainly flow downward along welding interface under coupling of shear stress and interfacial pressure,resulting in the interfacial plastic zone shown as an inverted“V”shape.The upward plastic flow in this direction is relatively weak,and only a small amount of flash is extruded from top of joint.Moreover,the wedge block and welding components at top of joint are always in un-steady friction stage,leading to nonuniform temperature field distribution and un-welded defects.According to the results of numerical simulation,high oscillating frequency combined with low pressure and small amplitude is considered as appropriate parameter selection scheme to improve the upward interfacial plastic flow at top of joint and suppress the un-welded defects.The results of thermo-physical simulation illustrate that continuous dynamic recrystallization(CDRX)induces the bonding of interface,accompanying by intense dislocation movement and creation of many low-angle grain boundaries.In the interfacial bonding area,grain orientation is random with relatively low texture density(5.0 mud)owing to CDRX.
基金supported by the National University of Defense Technology Independent Innovation Science Fund(Nos.24-ZZCXJDZ-42 and XJJC2024088)Hefei Municipal Natural Science Foundation(No.K130936103)The authors extend their gratitude to Ms.Zhang from Scientific Compass(www.shiyanjia.com)for providing invaluable assistance.
文摘The growing complexity of electromagnetic(EM)interference has driven significant demand for next-generation absorbers that combine lightweight,flexibility,and good electromagnetic attenuation capability.The low-dimensional ternary Co_(3)ZnC/Co/CNT composites with hollow structures have been synthesized through in-situ polymerization and high-temperature carbonization.The unique integration of low-dimensional nanostructures and multicomponent heterointerfaces confers exceptional EM absorption properties,achieving a reflection loss of−70.0 dB and significantly reducing radar cross section(RCS)scattering signals.It is particularly meaningful that the numerical simulation of Co_(3)ZnC/Co/CNT metama-terial reveals ultrawideband absorption performance,achieving 10.7 GHz(7.3-18.0 GHz)at a thickness of 4.5 mm and extending to 15 GHz(3.0-18.0 GHz)with a 10.5 mm.Moreover,the Co_(3)ZnC/Co/CNT composites retain meritorious EM absorption properties after flexible film formation,broadening their usability and application scope.These investigations will provide seminal insights encompassing theoretical validation,experimental synthesis,and practical application for the next generation of absorbers.
基金supported by the Zhejiang Provincial Key Research and Development Program of China(No.2020C01020).
文摘The nozzle is a critical component responsible for generating most of the net thrust in a scramjet engine.The quality of its design directly affects the performance of the entire propulsion system.However,most turbulence models struggle to make accurate predictions for subsonic and supersonic flows in nozzles.In this study,we explored a novel model,the algebraic stress model k-kL-ARSM+J,to enhance the accuracy of turbulence numerical simulations.This new model was used to conduct numerical simulations of the design and off-design performance of a 3D supersonic asymmetric truncated nozzle designed in our laboratory,with the aim of providing a realistic pattern of changes.The research indicates that,compared to linear eddy viscosity turbulence models such as k-kL and shear stress transport(SST),the k-kL-ARSM+J algebraic stress model shows better accuracy in predicting the performance of supersonic nozzles.Its predictions were identical to the experimental values,enabling precise calculations of the nozzle.The performance trends of the nozzle are as follows:as the inlet Mach number increases,both thrust and pitching moment increase,but the rate of increase slows down.Lift peaks near the design Mach number and then rapidly decreases.With increasing inlet pressure,the nozzle thrust,lift,and pitching moment all show linear growth.As the flight altitude rises,the internal flow field within the nozzle remains relatively consistent due to the same supersonic nozzle inlet flow conditions.However,external to the nozzle,the change in external flow pressure results in the nozzle exit transitioning from over-expanded to under-expanded,leading to a shear layer behind the nozzle that initially converges towards the nozzle center and then diverges.
基金National Natural Science Foundation of China(No.52071306)the Natural Science Foundation of Shandong Province(No.ZR2019MEE050)the Natural Science Foundation of Zhejiang Province(No.LZ22E090003).
文摘The local time-stepping(LTS)algorithm is an adaptive method that adjusts the time step by selecting suitable intervals for different regions based on the spatial scale of each cell and water depth and flow velocity between cells.The method can be optimized by calculating the maximum power of two of the global time step increments in the domain,allowing the optimal time step to be approached throughout the grid.To verify the acceleration and accuracy of LTS in storm surge simulations,we developed a model to simulate astronomical storm surges along the southern coast of China.This model employs the shallow water equations as governing equations,numerical discretization using the finite volume method,and fluxes calculated by the Roe solver.By comparing the simulation results of the traditional global time-stepping algorithm with those of the LTS algorithm,we find that the latter fit the measured data better.Taking the calculation results of Typhoon Sally in 1996 as an example,we show that compared with the traditional global time-stepping algorithm,the LTS algorithm reduces computation time by 2.05 h and increases computation efficiency by 2.64 times while maintaining good accuracy.
基金supported by the key technology project of China Southern Power Grid Corporation(GZKJXM20240009).
文摘As the integration of electronic components in high-performance servers increases,heat generation significantly impacts performance and raises failure rates.Therefore,heat dissipation has become a critical concern in electronic circuit design.This study uses numerical simulations to investigate the heat dissipation characteristics of electronic components in air-cooled servers.By adjusting airflow speed,heat sink configurations,and the arrangement of straight-fin heat sinks,we optimize heat dissipation performance and analyze the mechanisms at different airflow speeds.The results show that,at the same airflow speed,the temperature of the heat sink is lower than that of the electronic components,creating a temperature gradient that enhances heat transfer.Compared to a front-to-back arrangement of two straight-fin heat sinks,placing the heat sinks parallel to each other results in a lower maximum component temperature and better temperature uniformity.Heat sinks with fins significantly improve heat dissipation.The heat sink with semicylindrical fins on the rib surface provides the best cooling performance.Moreover,compared to natural convection,the maximum temperature of the electronic components decreases by 56.17%and 61%when the incoming flow velocity is 6 m/s with two parallel flat ribbed heat sinks and front-to-back arrangement,respectively.
基金supported by the National Key Research and Development Program of China(2022YFC2904400)Guangxi Science and Technology Major Project(Gui Ke AA23023033)。
文摘As a pyrometallurgical process,circulating fluidized bed(CFB) roasting has good potential for application in desulfurization of high-sulfur bauxite.The gas-solid distribution and reaction during CFB roasting of high-sulfur bauxite were simulated using the computational particle fluid dynamics(CPFD) method.The effect of primary air flow velocity on particle velocity,particle volume distribution,furnace temperature distribution and pressure distribution were investigated.Under the condition of the same total flow of natural gas,the impact of the number of inlets on the desulfurization efficiency,atmosphere mass fraction distribution and temperature distribution in the furnace was further investigated.
基金supported by the Sinopec Technology Research and Development Project(No.30000000-22-ZC0607-0235,No.33550000-22-ZC0607-0009)the National Natural Science Foundation of China(No.52334002).
文摘Natural gas hydrate widely exists in the South China Sea as clean energy.A three-phase transition layer widely exists in low permeability Class I hydrates in the Shenhu offshore area.Therefore,taking into account the low-permeability characteristics with an average permeability of 5.5 mD and moderate heterogeneity,a 3-D geological model of heterogeneous Class I hydrate reservoirs with three-phase transition layers is established by Kriging interpolation and stochastic modeling method,and a numerical simulation model is used to describe the depressurization production performance of the reservoir.With the development of depressurization,a specific range of complete decomposition zones appear both in the hydrate and transition layers.The entire decomposition zone of the whole reservoir tends to outward and upward diffusion.There is apparent methane escape in the three-phase transition layer.Due to the improvement of local permeability caused by the phase transition of hydrate dissociation,some methane accumulation occurs at the bottom of the hydrate layer,forming a local methane enrichment zone.The methane migration trends in reservoirs are mainly characterized by movement toward production wells and hydrate layers under the influence of gravity.However,due to the permeability limitation of hydrate reservoirs,many fluids have not been effectively produced and remain in the reservoir.Therefore,to improve the effective pressure drop of the reservoir,the perforation method and pressure reduction method were optimized by analyzing the influencing factors based on the gas production rate.The comparative study demonstrates that perforating through the free gas layer combined with one-time depressurization can enhance the effective depressurization and improve production performance.The gas production rate from perforating through the free gas layer can be twice as high as that from perforating through the transition layer.This study can provide theoretical support for the utilization of marine energy.
基金supported by the National Key Research and Development Project of China(No.2022YFC3902001)the National Natural Science Foundation of China(No.52004340)the Guangxi Innovation-driven Development Project,China(Nos.Gui 2021AA12006 and 2021AB26024)。
文摘The turbulent characteristics of the top-blown Laval nozzle and the influence of pressure and Mach number were studied through numerical simulation.With 2.72%error between the results and the empirical formula,the results are reliable.Nozzle fluid is influenced by pipe structure,causing pressure and density to drop as speed increases.Differences in pressure and velocity between the jet and surrounding gas lead to jet velocity attenuation,flow expansion,deflection,and eddy currents.The optimal top blowing pressure is 0.6 MPa,and the center velocity and width of the jet are 345 m/s and 0.124 m,respectively,at 20De(De is the nozzle exit diameter).It achieves a maximum jet velocity of 456 m/s.The optimal nozzle Mach number is 1.75,with a maximum jet velocity of 451 m/s.At 20D_(e),the jet center velocity is 338 m/s,with a width of 0.12 m.
文摘This paper presents a comprehensive experimental and numerical investigation of radiant floor heating(RFH)systems integrated with phase changematerial(PCM)-based thermal energy storage(TES).The study compares two underfloor pipe configurations:double serpentine and spiral.It also looks at how well a paraffin wax PCM system works with compact heat exchanger-type TES units during winter in Iraq.Key performance indicators including discharge temperature,heat transfer rate,liquid fraction evolution,and temperature uniformity were assessed through in situ experimental measurements and ANSYS fluent simulations.Results demonstrate that the spiral design provides slightly more uniform temperature distribution on the tile surface at an inlet water temperature of 55℃,with an average difference of approximately 0.5%,the serpentine layout exhibits higher slab temperature distribution by about 0.66%.Notably,the serpentine configuration shows superior thermal homogeneity and heat distribution,with a 15.05%increase in heat gain at a 55℃ inlet temperature compared to the spiral design.The performance gap between the two layouts narrows as the inlet temperature increases from 50℃ in 5℃ increments by approximately 4.1%,3.7%,and 1.7%,respectively.Higher inlet temperatures also improve PCM discharging and charging rates,improving energy storage utilization.The findings provide significant design guidelines for sustainable heating systems for cold climates.
基金Funded by the National Natural Science Foundation of China Academy of Engineering Physics and Jointly Setup"NSAF"Joint Fund(No.U1430119)。
文摘The multi-scale modeling combined with the cohesive zone model(CZM)and the molecular dynamics(MD)method were preformed to simulate the crack propagation in NiTi shape memory alloys(SMAs).The metallographic microscope and image processing technology were employed to achieve a quantitative grain size distribution of NiTi alloys so as to provide experimental data for molecular dynamics modeling at the atomic scale.Considering the size effect of molecular dynamics model on material properties,a reasonable modeling size was provided by taking into account three characteristic dimensions from the perspective of macro,meso,and micro scales according to the Buckinghamπtheorem.Then,the corresponding MD simulation on deformation and fracture behavior was investigated to derive a parameterized traction-separation(T-S)law,and then it was embedded into cohesive elements of finite element software.Thus,the crack propagation behavior in NiTi alloys was reproduced by the finite element method(FEM).The experimental results show that the predicted initiation fracture toughness is in good agreement with experimental data.In addition,it is found that the dynamics initiation fracture toughness increases with decreasing grain size and increasing loading velocity.
