Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing addit...Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.展开更多
The coupling between heat and pressure is the kernel of inertia friction welding(IFW)and is still not fully understood.A novel 3D fully coupled finite element model based on a plastic friction pair was developed to si...The coupling between heat and pressure is the kernel of inertia friction welding(IFW)and is still not fully understood.A novel 3D fully coupled finite element model based on a plastic friction pair was developed to simulate the IFW process of a Ni-based superalloy and reveal the omnidirectional thermo-mechanical coupling mechanism of the friction interface.The numerical model successfully simulated the deceleration,deformation processes,and peak torsional moments in IFW and captured the evolution of temperature,contact pressure,and stress.The simulated results were validated through measured thermal history,optical macrography,and axial shortening.The results indicated that interfacial friction heat was the primary heat source,and plastic deformation energy only accounted for 4%of the total.The increase in initial rotational speed and friction pressure elevated the peak temperature,reaching a maximum of 1525.5K at an initial rotational speed of 2000 r/min and friction pressure of 400 MPa.The interface heat generation could form an axial temperature gradient exceeding 320K/mm.The radial inhomogeneities of heat generation and temperature were manifested in a concentric ring distribution with maximum heat flux and temperature ranging from 2/5 to 2/3 radius.The radial inhomogeneities were caused by increasing linear velocity along the radius and an opposite distribution of contact pressure,which could reach 1.7 times the set pressure at the center.The circumferential inhomogeneity of thermomechanical distribution during rotary friction welding was revealed for the first time,benefiting from the 3D model.The deflection and transformation of distribution in contact pressure and Mises stress were indicators of plastic deformation and transition of quasi-steady state welding.The critical Mises stress was 0.5 times the friction pressure in this study.The presented modeling provides a reliable insight into the thermo-mechanical coupling mechanism of IFW and lays a solid foundation for predicting the microstructures and mechanical properties of inertia friction welded joints.展开更多
This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element couplin...This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.展开更多
The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatu...The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatures,pressures and complex geological environments of deep strata frequently result in the coupling of multiple physical fields,including mechanical,thermal and hydraulic fields,during the fracturing of rocks.This review initially presents an overview of the coupling mechanisms of these physical fields,thereby elucidating the interaction processes ofmechanical,thermal,and hydraulic fields within rockmasses.Secondly,an in-depth analysis ofmulti-field coupling is conducted from both spatial and temporal perspectives,with the introduction of simulation methods for a range of scales.It emphasizes cross-scale coupling methodologies for the transfer of rock properties and physical field data,including homogenization techniques,nested coupling strategies and data-driven approaches.To address the discontinuous characteristics of the rock fracture process,the review provides a detailed explanation of continuousdiscontinuous couplingmethods,to elucidate the evolution of rock fracturing and deformationmore comprehensively.In conclusion,the review presents a summary of the principal points,challenges and future directions of multi-field coupling simulation research.It also puts forward the potential of integrating intelligent algorithms with multi-scale simulation techniques to enhance the accuracy and efficiency of multi-field coupling simulations.This offers novel insights into multi-field coupling simulation analysis in deep rock masses.展开更多
A novel modeling technique based on the coupled Eulerian-Lagrangian(CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mec...A novel modeling technique based on the coupled Eulerian-Lagrangian(CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mechanical properties during spudcan penetration in normally consolidated clay soil. In the CEL model, the normal method of assigning an increasing shear strength profile with depth(NA) is defective due to its Eulerian framework. In this paper, a new technique is proposed to update soil material properties by introducing thermo-mechanical coupled analysis(TMCA) to the CEL models. During establishment of the CEL models, the optimal penetration velocity and minimum mesh size are determined through parametric studies. Reasonability and accuracy are then verified through comparison of the preliminary results with the soil flow configuration and penetration resistance(Fv) of a centrifuge test, and the results of the proposed method are compared with those of the remeshing and interpolation technique with small strain(RITSS) method. To achieve a CEL model with satisfactory accuracy, the NA and TMCA methods implemented in the CEL models and the RITSS method are first adopted in weightless soil. Comparison of the findings with those obtained in previous studies shows that the TMCA method can update material properties and predict Fv. The TMCA method is then applied to soils with self-weight and different shear strength profiles. Results show that the proposed method is capable of accurately modeling the large deformation problem of spudcan penetration in non-homogeneous clay.展开更多
3D thermo-meehanical coupled simulation of whole rolling process for 60 kg/m heavy rail was accomplished by FEM method. The finite element model, physical parameters of U75V and parameter setting of simulation were in...3D thermo-meehanical coupled simulation of whole rolling process for 60 kg/m heavy rail was accomplished by FEM method. The finite element model, physical parameters of U75V and parameter setting of simulation were introduced in detail. The whole rolling process of 60 kg/m heavy rail was divided into 27 time cells to simulate respectively, and the model rebuilding and temperature inheritance method in intermediate pass were proceeded. Then, based on simulation results, the workpiece deformation result, metal flow, stress and strain of 60 kg/m heavy rail for typical passes were obtained. The temperature variation curves of whole rolling process for section key points of 60 kg/m heavy rail were plotted, and the temperature falling law of whole rolling process for 60 kg/m heavy rail was studied. In addition, temperature distribution of 60 kg/m heavy rail after whole rolling process was analyzed, and the results showed that temperature was highest at center of rail head and lowest at fringe of rail base. Moreover, the simulation results and measured results of rolling force for 60 kg/m heavy rail were compared, and the regularity was in good agreement.展开更多
A finite volume method is applied to simulate a closed die hot forging process of a cylinder billet. Since variation and distribution of temperature play very important role in hot forging, the code involves a methodo...A finite volume method is applied to simulate a closed die hot forging process of a cylinder billet. Since variation and distribution of temperature play very important role in hot forging, the code involves a methodology of a coupled system of mechanical and thermal equations. The simulated results are compared with the experimental ones. The distribution of temperature in the billet obtained from the simulation is also discussed.展开更多
In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,th...In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,the random distribution characteristics of particles,including the shape,size,orientation,spatial location,and volume fractions,are all considered.Firstly,the repre.sentation for the microscopic configuration of the statistically inhomogeneous materials is described.Secondly,the SSOTS formulation for the dynamic thermo-mechanical coupled problem is proposed in a constructive way,including the cell problems,effective thermal and mechanical parameters,homogenized problems,and the SSOTS formulas of the temperatures,displacements,heat flux densities and stresses.And then the algorithm procedure corresponding to the SSOTS method is brought forward.The numerical results obtained by using the SSOTS algorithm are compared with those by classical methods.In addition,the thermo-mechanical coupling effect is studied by comparing the results of coupled case with those of uncoupled case.It demonstrates that the coupling effect on the temperatures,heat flux densities,displacements,and stresses is very distinct.The results show that the SSOTS method is valid to predict the dynamic thermo-mechanical coupled performances of statistically inhomogeneous materials.展开更多
In this paper,the degradation of mechanical properties of engineering cementitious composites(ECCs)at elevated temperatures and the failure of fiber are considered.A failure model under coupled thermo-mechanical loads...In this paper,the degradation of mechanical properties of engineering cementitious composites(ECCs)at elevated temperatures and the failure of fiber are considered.A failure model under coupled thermo-mechanical loads for ECC is developed based on bond-based peridynamics.A semi-discrete model is constructed to describe fiber–matrix interactions and simulate thermal failure in ECC.The peridynamic differential operator(PDDO)is utilized for non-local modeling of thermal fluid flow and heat transfer.A multi-rate explicit time integration method is adopted to address thermo-mechanical coupling over different time scales.Model validation is achieved through simulating transient heat transfer in a homogeneous plate,with results aligning with analytical solutions.The damage behavior of a heated ECC plate in a borehole and under a fire scenario is analyzed,providing insights for enhancing fire resistance and high-temperature performance of ECC materials and structures.展开更多
The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo val...The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo valve would reach 120℃ and the valve core and valve sleeve deform in a short amount of time. So the control precision of servo valve significantly decreases and the clamping stagnation phenomenon of valve core appears. In order to solve the problem of degraded control accuracy and clamping stagnation of servo valve under large temperature difference circumstance, the numerical simulation of heat-fluid-solid coupling by using finite element method is done. The simulation result shows that zero position leakage of servo valve is basically impacted by oil temperature and change of fit clearance. The clamping stagnation is caused by warpage-deformation and fit clearance reduction of the valve core and valve sleeve. The distribution roles of the temperature and thermal-deformation of shell, valve core and valve sleeve and the pressure, velocity and temperature field of flow channel are also analyzed. Zero position leakage and electromagnet＇s current when valve core moves in full-stroke are tested using Electro-hydraulic Servo-valve Characteristic Test-bed of an aerospace sciences and technology corporation. The experimental results show that the change law of experimental current at different oil temperatures is roughly identical to simulation current. The current curve of the electromagnet is smooth when oil temperature is below 80℃, but the amplitude of current significantly increases and the hairy appears when oil temperature is above 80℃. The current becomes smooth again after the warped valve core and valve sleeve are reground. It indicates that clamping stagnation is caused by warpage-deformation and fit clearance reduction of valve core and valve sleeve. This paper simulates and tests the heat-fluid-solid coupling of double flapper-nozzle servo valve, and the obtained results provide the reference value for the design of double flapper-nozzle force feedback servo valve.展开更多
Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulati...Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.展开更多
The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite elemen...The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.展开更多
The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separat...The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separate energy equations for the gas and solid phases coupled with an experimental kinetic model. A lumped kinetic model containing four main species CH4, O2, COx （CO2, CO）, and C2 （C2H4 and C2H6） was used with a plug flow reactor model as well. The results from the model agreed with the experimental data. The model was used to analyze the influence of temperature and feed gas composition on the conversion and selectivity of the reactor performance. The analytical results indicate that the conversion decreases, whereas, C2 selectivity increases by increasing gas hourly space velocity （GHSV） and the methane conversion also decreases by increasing the methane to oxygen ratio.展开更多
The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges tha...The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges that affect its accuracy and efficiency.Firstly,the implicit handling of fluid-solid boundary conditions through parameter averaging strategy often results in low simulation accuracy.Secondly,surface topography can introduce staircase diffraction noise when grid spacing is large.To address these issues,this paper presents a novel approach.We derive an implicit expression for fluidsolid boundary conditions based on average medium theory,translating explicit boundary conditions into model parameter modification.This enables implicit handling of fluid-solid boundaries by modifying the parameters near the boundary.Furthermore,to mitigate staircase diffraction noise,we employ multiple interface discretization based on the superposition method.This effectively suppresses staircase diffraction noise without requiring grid refinement.The efficacy of our method in accurately modeling wave propagation phenomena in fluid-solid coupled media is demonstrated by numerical examples.Results align well with those obtained using the spectral element method(SEM),with significant reduction in staircase diffraction noise.展开更多
Urbanization and eco-environment coupling is a research hotspot.Dynamic simulation of urbanization and eco-environment coupling needs to be improved because the processes of coupling are complex and statistical method...Urbanization and eco-environment coupling is a research hotspot.Dynamic simulation of urbanization and eco-environment coupling needs to be improved because the processes of coupling are complex and statistical methods are limited.Systems science and cross-scale coupling allow us to define the coupled urbanization and eco-environment system as an open complex giant system with multiple feedback loops.We review the current state of dynamic simulation of urbanization and eco-environment coupling and find that:(1)The use of dynamic simulation is an increasing trend,the relevant theory is being developed,and modeling processes are being improved;(2)Dynamic simulation technology has become diversified,refined,intelligent and integrated;(3)Simulation is mainly performed for three aspects of the coupling,multiple regions and multiple elements,local coupling and telecoupling,and regional synergy.However,we also found some shortcomings:(1)Basic theories are inadequately developed and insufficiently integrated;(2)The methods of unifying systems and sharing data are behind the times;(3)Coupling relations and the dynamic characteristics of the main driving elements are not fully understood or completely identified.Additionally,simulation of telecoupling does not quantify parameters and is not systemically unified,and therefore cannot be used to represent spatial synergy.In the future,we must promote communication between research networks,technology integration and data sharing to identify the processes governing change in coupled relations and in the main driving elements in urban agglomerations.Finally,we must build decision support systems to plan and ensure regional sustainable urbanization.展开更多
The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanica...The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.展开更多
Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated...Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.展开更多
In this research,a detailed multi-physics study has been carried out by numerically simulating a solid fractured gun barrel for 20 thermo-mechanical cycles.The numerical model is based on thermal effects,mechanical st...In this research,a detailed multi-physics study has been carried out by numerically simulating a solid fractured gun barrel for 20 thermo-mechanical cycles.The numerical model is based on thermal effects,mechanical stress fields and fatigue crack mechanics.Elastic-plastic material data of modified AISI 4340 at temperatures ranging from 25 to 1200℃and at strain rates of 4,16,32 and 48 s^(-1) was acquired from high-temperature compression tests.This was used as material property data in the simulation model.The boundary conditions applied are kept similar to the working gun barrel during continuous firing.A methodology has been provided to define thermo-mechanically active surface-to-surface type interface between the crack faces for a better approximation of stresses at the crack tip.Comparison of results from non-autofrettaged and autofrettaged simulation models provide useful information about the evolution of strains and stresses in the barrel at different points under combined thermo-mechanical loading cycles in both cases.The effect of thermal fatigue under already induced compressive yield due to autofrettage and the progressive degradation of the accumulated stresses due to thermo-mechanical cyclic loads on the internal surface of the gun barrel(mimicking the continuous firing scenario)has been analyzed.Comparison between energy release rate at tips of varying crack lengths due to cyclic thermo-mechanical loading in the non-autofrettaged and autofrettaged gun has been carried out.展开更多
Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal...Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal around the roadway was divided into the seepage open zone, seepage orientation zone, seepage decay zone and original seepage zone along the radial direction of the roadway. The loaded gassy coal was treated as a viscoelastic and plastic softened medium, and the mechanical behaviors of the viscoelastic zone, plastic softened zone and broken zone around the roadway were analyzed with the consideration of the loading creep, softening and expansion effect of the gassy coal. According to the law of conservation of mass and the Darcy law, the flow-solid coupled model for the gas transportation of the coal around the roadway was established considering the dynamic evolution of the adsorption characteristics, porosity and permeability of the coal, and the simulation software COMSOL was utilized to numerically simulate the stress state and gas flow regularity around the coal, which provided meaningful reference for investigating the stability of the coal and rock around the roadway.展开更多
In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing appr...In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.展开更多
基金National Key Research and Development Program of China(2022YFB4600902)Shandong Provincial Science Foundation for Outstanding Young Scholars(ZR2024YQ020)。
文摘Wire arc additive manufacturing(WAAM)has emerged as a promising approach for fabricating large-scale components.However,conventional WAAM still faces challenges in optimizing microstructural evolution,minimizing additive-induced defects,and alleviating residual stress and deformation,all of which are critical for enhancing the mechanical performance of the manufactured parts.Integrating interlayer friction stir processing(FSP)into WAAM significantly enhances the quality of deposited materials.However,numerical simulation research focusing on elucidating the associated thermomechanical coupling mechanisms remains insufficient.A comprehensive numerical model was developed to simulate the thermomechanical coupling behavior in friction stir-assisted WAAM.The influence of post-deposition FSP on the coupled thermomechanical response of the WAAM process was analyzed quantitatively.Moreover,the residual stress distribution and deformation behavior under both single-layer and multilayer deposition conditions were investigated.Thermal analysis of different deposition layers in WAAM and friction stir-assisted WAAM was conducted.Results show that subsequent layer deposition induces partial remelting of the previously solidified layer,whereas FSP does not cause such remelting.Furthermore,thermal stress and deformation analysis confirm that interlayer FSP effectively mitigates residual stresses and distortion in WAAM components,thereby improving their structural integrity and mechanical properties.
基金supported by the National Key Research and Development Program of China(Grant No.2022YFB3404904)。
文摘The coupling between heat and pressure is the kernel of inertia friction welding(IFW)and is still not fully understood.A novel 3D fully coupled finite element model based on a plastic friction pair was developed to simulate the IFW process of a Ni-based superalloy and reveal the omnidirectional thermo-mechanical coupling mechanism of the friction interface.The numerical model successfully simulated the deceleration,deformation processes,and peak torsional moments in IFW and captured the evolution of temperature,contact pressure,and stress.The simulated results were validated through measured thermal history,optical macrography,and axial shortening.The results indicated that interfacial friction heat was the primary heat source,and plastic deformation energy only accounted for 4%of the total.The increase in initial rotational speed and friction pressure elevated the peak temperature,reaching a maximum of 1525.5K at an initial rotational speed of 2000 r/min and friction pressure of 400 MPa.The interface heat generation could form an axial temperature gradient exceeding 320K/mm.The radial inhomogeneities of heat generation and temperature were manifested in a concentric ring distribution with maximum heat flux and temperature ranging from 2/5 to 2/3 radius.The radial inhomogeneities were caused by increasing linear velocity along the radius and an opposite distribution of contact pressure,which could reach 1.7 times the set pressure at the center.The circumferential inhomogeneity of thermomechanical distribution during rotary friction welding was revealed for the first time,benefiting from the 3D model.The deflection and transformation of distribution in contact pressure and Mises stress were indicators of plastic deformation and transition of quasi-steady state welding.The critical Mises stress was 0.5 times the friction pressure in this study.The presented modeling provides a reliable insight into the thermo-mechanical coupling mechanism of IFW and lays a solid foundation for predicting the microstructures and mechanical properties of inertia friction welded joints.
基金supported by grants from the National Natural Science Foundation of China (51478130)the Guangzhou Municipal Education Bureau’s Scientific Research Project, China (2024312217)+1 种基金the China Scholarship Council (201808440070)the 111 Project of China (D21021).
文摘This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.
基金supported by the National Natural Science Foundation of China(Grant Nos.42477185,41602308)the Zhejiang Provincial Natural Science Foundation of China(Grant No.LY20E080005)the Postgraduate Course Construction Project of Zhejiang University of Science and Technology(Grant No.2021yjskj05).
文摘The utilization of multi-field coupling simulation methods has become a pivotal approach for the investigation of intricate fracture behavior and interaction mechanisms of rock masses in deep strata.The high temperatures,pressures and complex geological environments of deep strata frequently result in the coupling of multiple physical fields,including mechanical,thermal and hydraulic fields,during the fracturing of rocks.This review initially presents an overview of the coupling mechanisms of these physical fields,thereby elucidating the interaction processes ofmechanical,thermal,and hydraulic fields within rockmasses.Secondly,an in-depth analysis ofmulti-field coupling is conducted from both spatial and temporal perspectives,with the introduction of simulation methods for a range of scales.It emphasizes cross-scale coupling methodologies for the transfer of rock properties and physical field data,including homogenization techniques,nested coupling strategies and data-driven approaches.To address the discontinuous characteristics of the rock fracture process,the review provides a detailed explanation of continuousdiscontinuous couplingmethods,to elucidate the evolution of rock fracturing and deformationmore comprehensively.In conclusion,the review presents a summary of the principal points,challenges and future directions of multi-field coupling simulation research.It also puts forward the potential of integrating intelligent algorithms with multi-scale simulation techniques to enhance the accuracy and efficiency of multi-field coupling simulations.This offers novel insights into multi-field coupling simulation analysis in deep rock masses.
基金supported by the National Natural Science Foundation of China (No. 51779236)the NSFC-Shandong Joint Fund Project (No. U1706226)funded by the China Scholarship Council (No. 201606330049)
文摘A novel modeling technique based on the coupled Eulerian-Lagrangian(CEL) method is provided to solve the geotechnical problems with large deformations. The technique is intended to solve the update problem of soil mechanical properties during spudcan penetration in normally consolidated clay soil. In the CEL model, the normal method of assigning an increasing shear strength profile with depth(NA) is defective due to its Eulerian framework. In this paper, a new technique is proposed to update soil material properties by introducing thermo-mechanical coupled analysis(TMCA) to the CEL models. During establishment of the CEL models, the optimal penetration velocity and minimum mesh size are determined through parametric studies. Reasonability and accuracy are then verified through comparison of the preliminary results with the soil flow configuration and penetration resistance(Fv) of a centrifuge test, and the results of the proposed method are compared with those of the remeshing and interpolation technique with small strain(RITSS) method. To achieve a CEL model with satisfactory accuracy, the NA and TMCA methods implemented in the CEL models and the RITSS method are first adopted in weightless soil. Comparison of the findings with those obtained in previous studies shows that the TMCA method can update material properties and predict Fv. The TMCA method is then applied to soils with self-weight and different shear strength profiles. Results show that the proposed method is capable of accurately modeling the large deformation problem of spudcan penetration in non-homogeneous clay.
文摘3D thermo-meehanical coupled simulation of whole rolling process for 60 kg/m heavy rail was accomplished by FEM method. The finite element model, physical parameters of U75V and parameter setting of simulation were introduced in detail. The whole rolling process of 60 kg/m heavy rail was divided into 27 time cells to simulate respectively, and the model rebuilding and temperature inheritance method in intermediate pass were proceeded. Then, based on simulation results, the workpiece deformation result, metal flow, stress and strain of 60 kg/m heavy rail for typical passes were obtained. The temperature variation curves of whole rolling process for section key points of 60 kg/m heavy rail were plotted, and the temperature falling law of whole rolling process for 60 kg/m heavy rail was studied. In addition, temperature distribution of 60 kg/m heavy rail after whole rolling process was analyzed, and the results showed that temperature was highest at center of rail head and lowest at fringe of rail base. Moreover, the simulation results and measured results of rolling force for 60 kg/m heavy rail were compared, and the regularity was in good agreement.
文摘A finite volume method is applied to simulate a closed die hot forging process of a cylinder billet. Since variation and distribution of temperature play very important role in hot forging, the code involves a methodology of a coupled system of mechanical and thermal equations. The simulated results are compared with the experimental ones. The distribution of temperature in the billet obtained from the simulation is also discussed.
基金supported by the National Natural Science Foundation of China(Grants 11471262,11202032)the Basic Research Project of National Defense(Grant B 1520132013)supported by the State Key Laboratory of Science and Engineering Computing and Center for high performance computing of Northwestem Polytechnical University
文摘In this paper,a statistical second-order twoscale（SSOTS） method is developed to simulate the dynamic thcrmo-mechanical performances of the statistically inhomogeneous materials.For this kind of composite material,the random distribution characteristics of particles,including the shape,size,orientation,spatial location,and volume fractions,are all considered.Firstly,the repre.sentation for the microscopic configuration of the statistically inhomogeneous materials is described.Secondly,the SSOTS formulation for the dynamic thermo-mechanical coupled problem is proposed in a constructive way,including the cell problems,effective thermal and mechanical parameters,homogenized problems,and the SSOTS formulas of the temperatures,displacements,heat flux densities and stresses.And then the algorithm procedure corresponding to the SSOTS method is brought forward.The numerical results obtained by using the SSOTS algorithm are compared with those by classical methods.In addition,the thermo-mechanical coupling effect is studied by comparing the results of coupled case with those of uncoupled case.It demonstrates that the coupling effect on the temperatures,heat flux densities,displacements,and stresses is very distinct.The results show that the SSOTS method is valid to predict the dynamic thermo-mechanical coupled performances of statistically inhomogeneous materials.
基金supported by National Natural Science Foundation of China(Grant Numbers 11872339,11472248).
文摘In this paper,the degradation of mechanical properties of engineering cementitious composites(ECCs)at elevated temperatures and the failure of fiber are considered.A failure model under coupled thermo-mechanical loads for ECC is developed based on bond-based peridynamics.A semi-discrete model is constructed to describe fiber–matrix interactions and simulate thermal failure in ECC.The peridynamic differential operator(PDDO)is utilized for non-local modeling of thermal fluid flow and heat transfer.A multi-rate explicit time integration method is adopted to address thermo-mechanical coupling over different time scales.Model validation is achieved through simulating transient heat transfer in a homogeneous plate,with results aligning with analytical solutions.The damage behavior of a heated ECC plate in a borehole and under a fire scenario is analyzed,providing insights for enhancing fire resistance and high-temperature performance of ECC materials and structures.
基金Supposed by National Natural Science Foundation of China(Grant No.51075348)Hebei Provincial Natural Science Foundation of China(Grant No.E2011203151)Research Fund for Doctoral Program of Higher Education of China(Grant No.20101333110002)
文摘The double flapper-nozzle servo valve is widely used to launch and guide the equipment. Due to the large instantaneous flow rate of servo valve working under specific operating conditions, the temperature of servo valve would reach 120℃ and the valve core and valve sleeve deform in a short amount of time. So the control precision of servo valve significantly decreases and the clamping stagnation phenomenon of valve core appears. In order to solve the problem of degraded control accuracy and clamping stagnation of servo valve under large temperature difference circumstance, the numerical simulation of heat-fluid-solid coupling by using finite element method is done. The simulation result shows that zero position leakage of servo valve is basically impacted by oil temperature and change of fit clearance. The clamping stagnation is caused by warpage-deformation and fit clearance reduction of the valve core and valve sleeve. The distribution roles of the temperature and thermal-deformation of shell, valve core and valve sleeve and the pressure, velocity and temperature field of flow channel are also analyzed. Zero position leakage and electromagnet＇s current when valve core moves in full-stroke are tested using Electro-hydraulic Servo-valve Characteristic Test-bed of an aerospace sciences and technology corporation. The experimental results show that the change law of experimental current at different oil temperatures is roughly identical to simulation current. The current curve of the electromagnet is smooth when oil temperature is below 80℃, but the amplitude of current significantly increases and the hairy appears when oil temperature is above 80℃. The current becomes smooth again after the warped valve core and valve sleeve are reground. It indicates that clamping stagnation is caused by warpage-deformation and fit clearance reduction of valve core and valve sleeve. This paper simulates and tests the heat-fluid-solid coupling of double flapper-nozzle servo valve, and the obtained results provide the reference value for the design of double flapper-nozzle force feedback servo valve.
基金financially supported by the Riset Pengembangan dan Penerapan(RPP),Diponegoro University 2023 research scheme with contract number 609-18/UN7.D2/PP/VIII/2023.
文摘Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.
基金supported by the Hebei Provincial Natural Science Foundation of China(No.E2007000591).
文摘The aim of this study was to simulate the solidification process of beam blank continuous casting, and then find the reasons for the typical defects of the beam blank. A two-dimensional transient coupled finite element model has been developed to compute the temperature and stress profile in beam blank continuous casting. The enthalpy method was used in the heat conduction equation. The thermo-mechanical property in the mushy zone was taken into consideration in this calculation. It is shown that at the mold exit the thickness of the shell had its maximum value at the flange tip and its minimum value at the fillet. The temperature had a great fluctuation on the surface of the beam blank in the secondary cooling zone. At the unbending point, the surface temperature of the web was in the brittleness temperature range under the present condition. To ensure the quality, it is necessary to weaken the intensity of secondary cooling. At the mold exit the equivalent stress and strain have higher values at the flange tip and at the web. From the spray 1 to the unbending point, the maximum values of stress and strain gradually moved to the internal section of the flange tip and the web. However, whenever, there were bigger stress and strain values near the flange tip and the web than in the other parts, it must be very easy to generate cracks at those positions. Now, online verification of this simulation has been developed, which has proved to be very useful and efficient to instruct the practical production of beam blank continuous casting.
文摘The oxidative coupling of methane （OCM） to ethylene over a perovskite titanate catalyst in a fixed bed reactor was studied experimentally and numerically. The two-dimensional steady state model accounted for separate energy equations for the gas and solid phases coupled with an experimental kinetic model. A lumped kinetic model containing four main species CH4, O2, COx （CO2, CO）, and C2 （C2H4 and C2H6） was used with a plug flow reactor model as well. The results from the model agreed with the experimental data. The model was used to analyze the influence of temperature and feed gas composition on the conversion and selectivity of the reactor performance. The analytical results indicate that the conversion decreases, whereas, C2 selectivity increases by increasing gas hourly space velocity （GHSV） and the methane conversion also decreases by increasing the methane to oxygen ratio.
基金supported by the National Natural Science Foundation of China(Nos.42404134,U24B2031,42174160)the China Postdoctoral Science Foundation(No.2024M753204)the National Key R&D Program of China(Nos.2021YFA0716901,2022YFB3904601)。
文摘The finite-difference method(FDM)is an essential tool in exploration geophysics,particularly for simulating wave propagation in fluid-solid coupled media.Despite its widespread use,FDM faces significant challenges that affect its accuracy and efficiency.Firstly,the implicit handling of fluid-solid boundary conditions through parameter averaging strategy often results in low simulation accuracy.Secondly,surface topography can introduce staircase diffraction noise when grid spacing is large.To address these issues,this paper presents a novel approach.We derive an implicit expression for fluidsolid boundary conditions based on average medium theory,translating explicit boundary conditions into model parameter modification.This enables implicit handling of fluid-solid boundaries by modifying the parameters near the boundary.Furthermore,to mitigate staircase diffraction noise,we employ multiple interface discretization based on the superposition method.This effectively suppresses staircase diffraction noise without requiring grid refinement.The efficacy of our method in accurately modeling wave propagation phenomena in fluid-solid coupled media is demonstrated by numerical examples.Results align well with those obtained using the spectral element method(SEM),with significant reduction in staircase diffraction noise.
基金Major Program of National Natural Science Foundation of China,No.41590840,No.41590842。
文摘Urbanization and eco-environment coupling is a research hotspot.Dynamic simulation of urbanization and eco-environment coupling needs to be improved because the processes of coupling are complex and statistical methods are limited.Systems science and cross-scale coupling allow us to define the coupled urbanization and eco-environment system as an open complex giant system with multiple feedback loops.We review the current state of dynamic simulation of urbanization and eco-environment coupling and find that:(1)The use of dynamic simulation is an increasing trend,the relevant theory is being developed,and modeling processes are being improved;(2)Dynamic simulation technology has become diversified,refined,intelligent and integrated;(3)Simulation is mainly performed for three aspects of the coupling,multiple regions and multiple elements,local coupling and telecoupling,and regional synergy.However,we also found some shortcomings:(1)Basic theories are inadequately developed and insufficiently integrated;(2)The methods of unifying systems and sharing data are behind the times;(3)Coupling relations and the dynamic characteristics of the main driving elements are not fully understood or completely identified.Additionally,simulation of telecoupling does not quantify parameters and is not systemically unified,and therefore cannot be used to represent spatial synergy.In the future,we must promote communication between research networks,technology integration and data sharing to identify the processes governing change in coupled relations and in the main driving elements in urban agglomerations.Finally,we must build decision support systems to plan and ensure regional sustainable urbanization.
基金Project(51605234)supported by the National Natural Science Foundation of ChinaProjects(2019JJ50510,2019JJ70077)supported by the Natural Science Foundation of Hunan Province,ChinaProjects(18B285,18B552)supported by Scientific Research Fund of Hunan Provincial Education Department,China。
文摘The interfacial heat transfer coefficient between hot profile surface and cooling water was determined by using inverse heat conduction model combined with end quenching experiment. Then, a Deform-3 D thermo-mechanical coupling model for simulating the on-line water quenching of extruded profile with unequal and large thicknesses was developed. The temperature field, residual stress field and distortion of profile during quenching were investigated systematically. The results show that heat transfer coefficient increases as water flow rate increases. The peak heat transfer coefficient with higher water flow rates appears at lower interface temperatures. The temperature distribution across the cross-section of profile during quenching is severe nonuniform and the maximum temperature difference is 300 ℃ at quenching time of 3.49 s. The temperature difference through the thickness of different parts of profile first increases sharply to a maximum value, and then gradually decreases. The temperature gradient increases obviously with the increase of thickness of parts. After quenching, there exist large residual stresses on the inner side of joints of profile and the two ends of part with thickness of 10 mm. The profile presents a twisting-type distortion across the cross-section under non-uniform cooling and the maximum twisting angle during quenching is 2.78°.
基金Project(20080431380) supported by China Postdoctoral Science Foundation
文摘Abstract： A joint solution model of variabk：-mass flow in two-phase region and fluid-solid coupling heat transfer, concerned about the charge process of variable-mass thermodynamic system, is built up and calculated by the finite element method （FEM）. The results are basically consistent with relative experimental data. The calculated average heat transfer coefficient reaches 1.7~105 W/（m2. K）. When the equal percentage valve is used, the system needs the minimum requirements of valve control, but brings the highest construction cost. With the： decrease of initial steam pressure, the heat transfer intensity also weakens but the steam flow increases. With the initial water filling coefficient increasing or the temperature of steam supply decreasing, the amount of accumulative steam flow increases with the growth of steam pressure. When the pressure of steam supply drops, the steam flow gradient increases during the maximum opening period of control valve, and causes the maximum steam flow to increase.
基金the DAAD Faculty Development for Ph.D.Candidates(Balochistan)2016(57245990)-HRDI-UESTP’s/UET’s funding scheme in cooperation with the Higher Education Commission of Pakistan(HEC)for sponsoring the stay at IMF TU Freiberg,Germany.
文摘In this research,a detailed multi-physics study has been carried out by numerically simulating a solid fractured gun barrel for 20 thermo-mechanical cycles.The numerical model is based on thermal effects,mechanical stress fields and fatigue crack mechanics.Elastic-plastic material data of modified AISI 4340 at temperatures ranging from 25 to 1200℃and at strain rates of 4,16,32 and 48 s^(-1) was acquired from high-temperature compression tests.This was used as material property data in the simulation model.The boundary conditions applied are kept similar to the working gun barrel during continuous firing.A methodology has been provided to define thermo-mechanically active surface-to-surface type interface between the crack faces for a better approximation of stresses at the crack tip.Comparison of results from non-autofrettaged and autofrettaged simulation models provide useful information about the evolution of strains and stresses in the barrel at different points under combined thermo-mechanical loading cycles in both cases.The effect of thermal fatigue under already induced compressive yield due to autofrettage and the progressive degradation of the accumulated stresses due to thermo-mechanical cyclic loads on the internal surface of the gun barrel(mimicking the continuous firing scenario)has been analyzed.Comparison between energy release rate at tips of varying crack lengths due to cyclic thermo-mechanical loading in the non-autofrettaged and autofrettaged gun has been carried out.
基金the financial support from the National Natural Science Foundation for Young Scientists of China (Nos.51604116 and 51604096)Natural Science Foundation ofHenbei Province (No.E2016508036)+1 种基金Hebei State Key Laboratory of Mine Disaster Prevention (No.KJZH2017K08)Basic and Frontier Technology Research Project of Henan Province in 2016 (No.162300410031)
文摘Due to the change of initial stress state caused by roadway excavation, the permeability of the coal body may be changed during the excavation process. In this paper, according to the different stress states, the coal around the roadway was divided into the seepage open zone, seepage orientation zone, seepage decay zone and original seepage zone along the radial direction of the roadway. The loaded gassy coal was treated as a viscoelastic and plastic softened medium, and the mechanical behaviors of the viscoelastic zone, plastic softened zone and broken zone around the roadway were analyzed with the consideration of the loading creep, softening and expansion effect of the gassy coal. According to the law of conservation of mass and the Darcy law, the flow-solid coupled model for the gas transportation of the coal around the roadway was established considering the dynamic evolution of the adsorption characteristics, porosity and permeability of the coal, and the simulation software COMSOL was utilized to numerically simulate the stress state and gas flow regularity around the coal, which provided meaningful reference for investigating the stability of the coal and rock around the roadway.
文摘In this paper, the iterative coupling approach is proposed for applications to solving multiphase flow equation systems in reservoir simulation, as it provides a more flexible time-stepping strategy than existing approaches. The iterative method decouples the whole equation systems into pressure and saturation/concentration equations, and then solves them in sequence, implicitly and semi-implicitly. At each time step, a series of iterations are computed, which involve solving linearized equations using specific tolerances that are iteration dependent. Following convergence of subproblems, material balance is checked. Convergence of time steps is based on material balance errors. Key components of the iterative method include phase scaling for deriving a pressure equation and use of several advanced numerical techniques. The iterative model is implemented for parallel computing platforms and shows high parallel efficiency and scalability.
