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 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°.展开更多
This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,...This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,and the material coefficients are quasi-periodic,i.e.,they depend not only on the microscale information but also on the macro location.Also,a mutual interaction between displacement and temperature fields is considered in the problem,which is our particular interest in this study.The multiscale asymptotic expansions of the temperature and displacement fields are constructed and associated error estimation in nearly pointwise sense is presented.Then,a finite element-difference algorithm based on the multiscale analysis method is brought forward in detail.Finally,some numerical examples are given.And the numerical results show that the multiscale method presented in this paper is effective and reliable to study the nonlinear thermo-mechanical coupling problem of composite structures with quasiperiodic properties.展开更多
Based on the theory of elastic-plastic finite element method, the high-speed hot continuous rolling process of a billet is simulated and analyzed in vertical and horizontal passes. The billet is dragged into the passe...Based on the theory of elastic-plastic finite element method, the high-speed hot continuous rolling process of a billet is simulated and analyzed in vertical and horizontal passes. The billet is dragged into the passes by contact friction force between the billet and rollers. The rollers and billet are represented by respectively rigid and deformable bodies, and three-dimensional models are developed for the billet and rollers. The distribution of deformation field, effective strain, rolling force and temperature field are accurately calculated for the whole rolling process （including unstable and stable stages）. In addition, the rolling pressure on the width symmetry center is compared with that in the in-situ experimental measurements. It is revealed that various heat exchange phenomena among the billet, rollers and surroundings can result in unbalanced temperature distribution on the cross section. Rolling force and strain can change significantly when the billet is moved towards or away from the roller gap, and keep almost invariable in the stable stage. It is expected that the simulation results would be useful for practical manufacture and provide the theoretical foundation for improvement of process planning and optimization of process parameters.展开更多
Semiconductor-based electronic devices usually work under multiphysics fields rendering complex electromagnetic-thermo-mechanical coupling.In this work,we develop a penalty function method based on a finite element an...Semiconductor-based electronic devices usually work under multiphysics fields rendering complex electromagnetic-thermo-mechanical coupling.In this work,we develop a penalty function method based on a finite element analysis to tackle this coupling behavior in a metal/semiconductor bilayer plate-the representative unit of semiconductor antenna,which receives strong and pulsed electromagnetic signals.Under these pulses,eddy current is generated,of which the magnitude varies remarkably from one plate to another due to the difference in electrical conductivity.In the concerned system,the metal layer generates much larger current,resulting in the large temperature rise and the nonnegligible Lorentz force,which could lead to delamination and failure of the semiconductor-based electronic device.This study provides theoretical guidance for the design and protection of semiconductor-based electronic devices in complex environments.展开更多
Permafrost is widely distributed in China and around the world.In permafrost regions,soil frost heave and thawing are severe and frequent,and can destabilize pile foundations.To this end,a finite element model of a si...Permafrost is widely distributed in China and around the world.In permafrost regions,soil frost heave and thawing are severe and frequent,and can destabilize pile foundations.To this end,a finite element model of a single pile in frozen soil is established to investigate the frost heave and frost jacking response to ensure its safety in the Qinghai-Tibet Plateau.Firstly,a hydro-thermal coupling model of a single pile in frozen soil is established based on coupling parameters and initial and boundary conditions.Then the temperature and moisture distributions are analyzed through the established coupling model.A hydro-thermo-mechanical coupling model is developed by importing the ice content and temperature results.Simulation results indicate that the amount of frost heave is greater at locations closer to the ground surface,and the displacement is smaller for frozen soil that is closer to the side of the pile.The results on frost jacking behavior of piles from this study can serve as a reference for the design,construction and maintenance of foundations.展开更多
The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear...The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.展开更多
Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon di...Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.展开更多
The transient finite element technique is applied, and a transient heat conduction model of wet brake friction disk is established. For obtaining the accurate heat flow density mathematic model and avoiding possibly i...The transient finite element technique is applied, and a transient heat conduction model of wet brake friction disk is established. For obtaining the accurate heat flow density mathematic model and avoiding possibly instable thermoelastic stress produced by the non uniform contact pressure of friction pair, a test method is applied to collect accurate contact pressure between the dual sheet steel and friction disk in the combining process. And then the heat-flow density and transient ther mo mechanical coupling simulation are analyzed. At the same time all possible boundary conditions are considered, such as the heat generation, heat conduction problem, relation between friction and contact, variation in load and heat change problem etc. The simulation results show that the me chanical model of thermo mechanical coupling can express well the dynamic characteristics of fric tion disk, and gives perfect reference for more study on thermoelastic distortion of brake friction pairs.展开更多
The structure optimization design under thermo-mechanical coupling is a difficult problem in the topology optimization field.An adaptive growth algorithm has become a more effective approach for structural topology op...The structure optimization design under thermo-mechanical coupling is a difficult problem in the topology optimization field.An adaptive growth algorithm has become a more effective approach for structural topology optimization.This paper proposed a topology optimization method by an adaptive growth algorithm for the stiffener layout design of box type load-bearing components under thermo-mechanical coupling.Based on the stiffness diffusion theory,both the load stiffness matrix and the heat conduction stiffness matrix of the stiffener are spread at the same time to make sure the stiffener grows freely and obtain an optimal stiffener layout design.Meanwhile,the objectives of optimization are the minimization of strain energy and thermal compliance of the whole structure,and thermo-mechanical coupling is considered.Numerical studies for square shells clearly show the effectiveness of the proposed method for stiffener layout optimization under thermo-mechanical coupling.Finally,the method is applied to optimize the stiffener layout of box type load-bearing component of themachining center.The optimization results show that both the structural deformation and temperature of the load-bearing component with the growth stiffener layout,which are optimized by the adaptive growth algorithm,are less than the stiffener layout of shape‘#’stiffener layout.It provides a new solution approach for stiffener layout optimization design of box type load-bearing components under thermo-mechanical coupling.展开更多
The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidanc...The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidance system.This paper presents a new method that uses the Finite Difference Time Domain(FDTD)method to calculate the electrical performance of radomes under Thermo-Mechanical-Electrical(TME)coupling.This method can accurately characterize the effects of material dielectric temperature drift and structural deformation on the electrical performance of the radome under flight conditions,enabling high-precision full-wave calculations of the broadband electrical performance of the radome.The method initiates by utilizing a Finite Element Grid Model(FE-GM)of the radome to sequentially acquire the radome's response temperature field and structural deformation field through thermal and mechanical simulations.Subsequently,spatial mapping techniques are developed to accurately incorporate the dielectric temperature drift and structural deformation of the radome into its Yee grid Electromagnetic(EM)simulation model.A verification case was designed to test the proposed method,and the results confirmed its high computational accuracy.Additionally,the effectiveness and necessity of the method were further demonstrated by analyzing the electrical performance of a fused silica ceramic radome used on a high-speed aircraft.展开更多
During high-speed flight,both thermal and mechani-cal loads can degrade the electrical performance of the antenna-radome system,which can subsequently affect the performance of the guidance system.This paper presents ...During high-speed flight,both thermal and mechani-cal loads can degrade the electrical performance of the antenna-radome system,which can subsequently affect the performance of the guidance system.This paper presents a method for evalu-ating the electrical performance of the radome when subjected to thermo-mechanical-electrical(TME)coupling.The method involves establishing a TME coupling model(TME-CM)based on the TME sharing mesh model(TME-SMM)generated by the tetrahedral mesh partitioning of the radome structure.The effects of dielectric temperature drift and structural deformation on the radome’s electrical performance are also considered.Firstly,the temperature field of the radome is obtained by tran-sient thermal analysis while the deformation field of the radome is obtained by static analysis.Subsequently,the dielectric varia-tion and structural deformation of the radome are accurately incorporated into the electrical simulation model based on the TME-SMM.The three-dimensional(3D)ray tracing method with the aperture integration technique is used to calculate the radome’s electrical performance.A representative example is provided to illustrate the superiority and necessity of the pro-posed method.This is achieved by calculating and analyzing the changes in the radome’s electrical performance over time dur-ing high-speed flight.展开更多
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.展开更多
The wide application of carbon fiber reinforced plastic(CFRP)components in modern aerospace manufacturing field puts high demands on the manufacturing process.Especially,the temperature increase during continuous mill...The wide application of carbon fiber reinforced plastic(CFRP)components in modern aerospace manufacturing field puts high demands on the manufacturing process.Especially,the temperature increase during continuous milling process becomes a key factor affecting the performance of composites,and the high milling temperature induces a variety of processing defects.This paper obtained the temperature variation data during the end milling process of CFRP laminates through experiments.After data fitting,the data were transformed into a function of heat flux density varying with time.In the finite element analysis,a double-ellipsoid moving heat source model was introduced,and a moving heat source subrou-tine was written based on the time-varying function of heat flux density to more accurately describe the thermal effects dur-ing the milling process and simulate the changes in the temperature field during milling.The Hashin failure criterion is a-dopted as the basis of fiber and matrix failure,and the simulation results of the temperature field are input into the thermal-force coupling simulation model as the predefined field conditions for solving and analyzing by means of sequential thermal-force coupling,so as to establish a thermal-force coupling simulation and analysis model for milling processing of CFRP end faces.The model simulation results can provide a basis for exploring the damage evolution law of CFRP material under the influence of temperature.展开更多
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.展开更多
Defects-rich heterointerfaces integrated with adjustable crystalline phases and atom vacancies,as well as veiled dielectric-responsive character,are instrumental in electromagnetic dissipation.Conventional methods,how...Defects-rich heterointerfaces integrated with adjustable crystalline phases and atom vacancies,as well as veiled dielectric-responsive character,are instrumental in electromagnetic dissipation.Conventional methods,however,constrain their delicate constructions.Herein,an innovative alternative is proposed:carrageenan-assistant cations-regulated(CACR)strategy,which induces a series of sulfides nanoparticles rooted in situ on the surface of carbon matrix.This unique configuration originates from strategic vacancy formation energy of sulfides and strong sulfides-carbon support interaction,benefiting the delicate construction of defects-rich heterostructures in M_(x)S_(y)/carbon composites(M-CAs).Impressively,these generated sulfur vacancies are firstly found to strengthen electron accumulation/consumption ability at heterointerfaces and,simultaneously,induct local asymmetry of electronic structure to evoke large dipole moment,ultimately leading to polarization coupling,i.e.,defect-type interfacial polarization.Such“Janus effect”(Janus effect means versatility,as in the Greek two-headed Janus)of interfacial sulfur vacancies is intuitively confirmed by both theoretical and experimental investigations for the first time.Consequently,the sulfur vacancies-rich heterostructured Co/Ni-CAs displays broad absorption bandwidth of 6.76 GHz at only 1.8 mm,compared to sulfur vacancies-free CAs without any dielectric response.Harnessing defects-rich heterostructures,this one-pot CACR strategy may steer the design and development of advanced nanomaterials,boosting functionality across diverse application domains beyond electromagnetic response.展开更多
Grinding is a generally utilized method,removing excess materials through effective abrasives.The grinding abrasives with multiple shape-position characteristics play a dominant role in determining the thermo-mechanic...Grinding is a generally utilized method,removing excess materials through effective abrasives.The grinding abrasives with multiple shape-position characteristics play a dominant role in determining the thermo-mechanical coupling,which may influence the surface quality directly.To investigate this correlated influence mechanism,this paper focuses on the abrasive shape and position characteristic on the grinding thermo-mechanical process with the analytic single abrasive interaction force model by considering the abrasive shape and its distribution information.It can be found that the mapped dynamic grinding temperature is actually discretized on the workpiece surface,which is on account of the diversity of the abrasive shape and its distribution.Moreover,higher spherical and conical abrasive particles,as well as lower pyramid shaped abrasive particle ratios,can generate greater specific grinding energy with a discretized temperature distribution,when compared with a higher proportion of pyramid shaped abrasive.The study can be utilized to provide valuable theoretical foundation for engineering practice by preparing structural wheel and its grinding property.展开更多
The occurrence of top-down(TD)cracking has gradually become a prevalent issue in semi-rigid base asphalt pavements after prolonged service.A coupled simulation model integrating the finite difference method(FDM)and di...The occurrence of top-down(TD)cracking has gradually become a prevalent issue in semi-rigid base asphalt pavements after prolonged service.A coupled simulation model integrating the finite difference method(FDM)and discrete element method(DEM)was employed to investigate the mechanical behavior of asphalt pavement containing a pre-existing TD crack.The mesoscopic parameters of the model were calibrated based on the mixture modulus and the static mechanical response on the MLS66 test road.Finally,an analysis was performed to assess how variations in TD crack depth and longitudinal length affect the distribution patterns of transverse tensile stress,vertical shear stress,and vertical compressive stress.The results indicate that the vertical propagation of TD crack significantly increases both the tensile stress value and range on the middle surface,while the longitudinal development of TD crack has minimal impact.This phenomenon may result in more severe fatigue failure on the middle surface.With the vertical and longitudinal development of TD crack,the vertical shear stress and compressive stress show obvious"two-stage"characteristics.When the crack's vertical length reaches 40 mm,there is a sharp increase in stress on the upper surface.As the crack continues to propagate vertically,the growth of stress on the upper surface becomes negligible,while the stress in the middle and lower layers increased significantly.Conversely,for longitudinal development of TD crack,any changes in stress are insignificant when their length is less than 180 mm;however,as they continue to develop longitudinally beyond this threshold,there is a sharp increase in stress levels.These findings hold great significance for understanding pavement structure deterioration and maintenance behavior associated with TD crack.展开更多
Electrochemical conversion of lignin for the production of high-value heterocyclic aromatic compounds has great potential.We demonstrate the targeted synthesis and cation modulation of NiCo_(2)O_(4)spinel nanoboxes,sy...Electrochemical conversion of lignin for the production of high-value heterocyclic aromatic compounds has great potential.We demonstrate the targeted synthesis and cation modulation of NiCo_(2)O_(4)spinel nanoboxes,synthesized via cation exchange and calcination oxidation.These catalysts exhibit excellent efficacy in the electrocatalytic conversion of lignin model compounds,specifically 2-phenoxy-1-phenylethanol,into nitrogen-containing aromatics,achieving high conversion rates and selectivities.These catalysts were synthesized via a cation exchange and calcination oxidation process,using Prussian blue nanocubes as precursors.The porous architecture and polymetallic composition of the NiCo_(2)O_(4)spinel demonstrated superior performance in electrocatalytic oxidative coupling,achieving a 99.2 wt%conversion rate of the 2-phenoxy-1-phenylethanol with selectivities of 37.5 wt%for quinoline derivatives and 31.5 wt%for phenol.Key innovations include the development of a sustainable one-pot synthesis method for quinoline derivatives,the elucidation of a multistage reaction pathway involving CAO bond cleavage,hydroxyaldol condensation,and CAN bond formation,and a deeper mechanistic understanding derived from DFT simulations.This work establishes a new strategy for lignin valorization,offering a sustainable route to produce high-value nitrogen-containing aromatics from renewable biomass under mild conditions,without the need for additional reagents.展开更多
Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression an...Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.展开更多
基金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.
基金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°.
基金financially supported by the National Natural Science Foundation of China(11501449)the Fundamental Research Funds for the Central Universities(3102017zy043)+2 种基金the China Postdoctoral Science Foundation(2016T91019)the fund of the State Key Laboratory of Solidification Processing in NWPU(SKLSP201628)the Scientific Research Program Funded by Shaanxi Provincial Education Department(14JK1353).
文摘This paper reports a multiscale analysis method to predict the thermomechanical coupling performance of composite structures with quasi-periodic properties.In these material structures,the configurations are periodic,and the material coefficients are quasi-periodic,i.e.,they depend not only on the microscale information but also on the macro location.Also,a mutual interaction between displacement and temperature fields is considered in the problem,which is our particular interest in this study.The multiscale asymptotic expansions of the temperature and displacement fields are constructed and associated error estimation in nearly pointwise sense is presented.Then,a finite element-difference algorithm based on the multiscale analysis method is brought forward in detail.Finally,some numerical examples are given.And the numerical results show that the multiscale method presented in this paper is effective and reliable to study the nonlinear thermo-mechanical coupling problem of composite structures with quasiperiodic properties.
基金Funded by National Natural Science Foundation of China (No. 51004047)Scientific Research Fund of Hunan Provincial Education Department (No. 10B020)Provincial Natural Science Foundation of Hunan (No. 09jj4024)
文摘Based on the theory of elastic-plastic finite element method, the high-speed hot continuous rolling process of a billet is simulated and analyzed in vertical and horizontal passes. The billet is dragged into the passes by contact friction force between the billet and rollers. The rollers and billet are represented by respectively rigid and deformable bodies, and three-dimensional models are developed for the billet and rollers. The distribution of deformation field, effective strain, rolling force and temperature field are accurately calculated for the whole rolling process （including unstable and stable stages）. In addition, the rolling pressure on the width symmetry center is compared with that in the in-situ experimental measurements. It is revealed that various heat exchange phenomena among the billet, rollers and surroundings can result in unbalanced temperature distribution on the cross section. Rolling force and strain can change significantly when the billet is moved towards or away from the roller gap, and keep almost invariable in the stable stage. It is expected that the simulation results would be useful for practical manufacture and provide the theoretical foundation for improvement of process planning and optimization of process parameters.
基金the National Natural Science Foundation of China(Grant nos.11772294,11621062)the Fundamental Research Funds for the Central Universities(Grant no.2017QNA4031).
文摘Semiconductor-based electronic devices usually work under multiphysics fields rendering complex electromagnetic-thermo-mechanical coupling.In this work,we develop a penalty function method based on a finite element analysis to tackle this coupling behavior in a metal/semiconductor bilayer plate-the representative unit of semiconductor antenna,which receives strong and pulsed electromagnetic signals.Under these pulses,eddy current is generated,of which the magnitude varies remarkably from one plate to another due to the difference in electrical conductivity.In the concerned system,the metal layer generates much larger current,resulting in the large temperature rise and the nonnegligible Lorentz force,which could lead to delamination and failure of the semiconductor-based electronic device.This study provides theoretical guidance for the design and protection of semiconductor-based electronic devices in complex environments.
基金supported by the National Natural Science Foundation of China(Nos.42071078,41731281 and 41701068)the Natural Science Foundation of Qinghai Province,China(No.2021-ZJ-908).
文摘Permafrost is widely distributed in China and around the world.In permafrost regions,soil frost heave and thawing are severe and frequent,and can destabilize pile foundations.To this end,a finite element model of a single pile in frozen soil is established to investigate the frost heave and frost jacking response to ensure its safety in the Qinghai-Tibet Plateau.Firstly,a hydro-thermal coupling model of a single pile in frozen soil is established based on coupling parameters and initial and boundary conditions.Then the temperature and moisture distributions are analyzed through the established coupling model.A hydro-thermo-mechanical coupling model is developed by importing the ice content and temperature results.Simulation results indicate that the amount of frost heave is greater at locations closer to the ground surface,and the displacement is smaller for frozen soil that is closer to the side of the pile.The results on frost jacking behavior of piles from this study can serve as a reference for the design,construction and maintenance of foundations.
基金supported by the National Natural Science Foundation of China(Grant No.51974173)the Natural Science Foundation of Shandong Province,China(Grant No.ZR2020QD122).
文摘The shear failure of intact rock under thermo-mechanical(TM)coupling conditions is common,such as in enhanced geothermal mining and deep mine construction.Under the effect of a continuous engineering disturbance,shear-formed fractures are prone to secondary instability,posing a severe threat to deep engineering.Although numerous studies regarding three-dimensional(3D)morphologies of fracture surfaces have been conducted,the understanding of shear-formed fractures under TM coupling conditions is limited.In this study,direct shear tests of intact granite under various TM coupling conditions were conducted,followed by 3D laser scanning tests of shear-formed fractures.Test results demonstrated that the peak shear strength of intact granite is positively correlated with the normal stress,whereas it is negatively correlated with the temperature.The internal friction angle and cohesion of intact granite significantly decrease with an increase in the temperature.The anisotropy,roughness value,and height of the asperities on the fracture surfaces are reduced as the normal stress increases,whereas their variation trends are the opposite as the temperature increases.The macroscopic failure mode of intact granite under TM coupling conditions is dominated by mixed tensileeshear and shear failures.As the normal stress increases,intragranular fractures are developed ranging from a local to a global distribution,and the macroscopic failure mode of intact granite changes from mixed tensileeshear to shear failure.Finally,3D morphological characteristics of the asperities on the shear-formed fracture surfaces were analyzed,and a quadrangular pyramid conceptual model representing these asperities was proposed and sufficiently verified.
基金the National Key R&D Program of China(No.2019YFB1504102).
文摘Geothermal energy is a kind of renewable,sustainable and clean energy resource.Geothermal energy is abundant in carbonate reservoirs.However,low matrix permeability limits its exploitation.The super-critical carbon dioxide(SC-CO_(2))jet fracturing is expected to efficiently stimulate the carbonate geothermal reservoirs and achieve the storage of CO_(2) simultaneously.In this paper,we established a transient seepage and fluid-thermo-mechanical coupled model to analyze the impact performance of sc-CO_(2) jet fracturing.The mesh-based parallel code coupling interface was employed to couple the fluid and solid domains by exchanging the data through the mesh interface.The physical properties change of sC-CO_(2) with temperature were considered in the numerical model.Results showed that SC-CO_(2) jet frac-turing is superior to water-jet fracturing with respect to jetting velocity,particle trajectory and pene-trability.Besides,stress distribution on the carbonate rock showed that the tensile and shear failure would more easily occur by SC-CO_(2) jet than that by water jet.Moreover,pressure and temperature control the jet field and seepage field of sC-CO_(2) simultaneously.Increasing the jet temperature can effectively enhance the impingement effect and seepage process by decreasing the viscosity and density of SC-CO_(2).The key findings are expected to provide a theoretical basis and design reference for applying SC-CO_(2) jet fracturing in carbonate geothermal reservoirs.
基金Supported by the National Basic Research Program of China("973"Program)(613002)
文摘The transient finite element technique is applied, and a transient heat conduction model of wet brake friction disk is established. For obtaining the accurate heat flow density mathematic model and avoiding possibly instable thermoelastic stress produced by the non uniform contact pressure of friction pair, a test method is applied to collect accurate contact pressure between the dual sheet steel and friction disk in the combining process. And then the heat-flow density and transient ther mo mechanical coupling simulation are analyzed. At the same time all possible boundary conditions are considered, such as the heat generation, heat conduction problem, relation between friction and contact, variation in load and heat change problem etc. The simulation results show that the me chanical model of thermo mechanical coupling can express well the dynamic characteristics of fric tion disk, and gives perfect reference for more study on thermoelastic distortion of brake friction pairs.
基金supported by National Natural Science Foundation of China (No.52075445)Science,Technology and Innovation Commission of Shenzhen Municipality (No.JCYJ20190806151013025).
文摘The structure optimization design under thermo-mechanical coupling is a difficult problem in the topology optimization field.An adaptive growth algorithm has become a more effective approach for structural topology optimization.This paper proposed a topology optimization method by an adaptive growth algorithm for the stiffener layout design of box type load-bearing components under thermo-mechanical coupling.Based on the stiffness diffusion theory,both the load stiffness matrix and the heat conduction stiffness matrix of the stiffener are spread at the same time to make sure the stiffener grows freely and obtain an optimal stiffener layout design.Meanwhile,the objectives of optimization are the minimization of strain energy and thermal compliance of the whole structure,and thermo-mechanical coupling is considered.Numerical studies for square shells clearly show the effectiveness of the proposed method for stiffener layout optimization under thermo-mechanical coupling.Finally,the method is applied to optimize the stiffener layout of box type load-bearing component of themachining center.The optimization results show that both the structural deformation and temperature of the load-bearing component with the growth stiffener layout,which are optimized by the adaptive growth algorithm,are less than the stiffener layout of shape‘#’stiffener layout.It provides a new solution approach for stiffener layout optimization design of box type load-bearing components under thermo-mechanical coupling.
文摘The electrical performance of radomes on high-speed aircraft can be influenced by the thermal and mechanical loads produced during high-speed flight,which can affect the detection dis-tance and accuracy of the guidance system.This paper presents a new method that uses the Finite Difference Time Domain(FDTD)method to calculate the electrical performance of radomes under Thermo-Mechanical-Electrical(TME)coupling.This method can accurately characterize the effects of material dielectric temperature drift and structural deformation on the electrical performance of the radome under flight conditions,enabling high-precision full-wave calculations of the broadband electrical performance of the radome.The method initiates by utilizing a Finite Element Grid Model(FE-GM)of the radome to sequentially acquire the radome's response temperature field and structural deformation field through thermal and mechanical simulations.Subsequently,spatial mapping techniques are developed to accurately incorporate the dielectric temperature drift and structural deformation of the radome into its Yee grid Electromagnetic(EM)simulation model.A verification case was designed to test the proposed method,and the results confirmed its high computational accuracy.Additionally,the effectiveness and necessity of the method were further demonstrated by analyzing the electrical performance of a fused silica ceramic radome used on a high-speed aircraft.
文摘During high-speed flight,both thermal and mechani-cal loads can degrade the electrical performance of the antenna-radome system,which can subsequently affect the performance of the guidance system.This paper presents a method for evalu-ating the electrical performance of the radome when subjected to thermo-mechanical-electrical(TME)coupling.The method involves establishing a TME coupling model(TME-CM)based on the TME sharing mesh model(TME-SMM)generated by the tetrahedral mesh partitioning of the radome structure.The effects of dielectric temperature drift and structural deformation on the radome’s electrical performance are also considered.Firstly,the temperature field of the radome is obtained by tran-sient thermal analysis while the deformation field of the radome is obtained by static analysis.Subsequently,the dielectric varia-tion and structural deformation of the radome are accurately incorporated into the electrical simulation model based on the TME-SMM.The three-dimensional(3D)ray tracing method with the aperture integration technique is used to calculate the radome’s electrical performance.A representative example is provided to illustrate the superiority and necessity of the pro-posed method.This is achieved by calculating and analyzing the changes in the radome’s electrical performance over time dur-ing high-speed flight.
基金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.
文摘The wide application of carbon fiber reinforced plastic(CFRP)components in modern aerospace manufacturing field puts high demands on the manufacturing process.Especially,the temperature increase during continuous milling process becomes a key factor affecting the performance of composites,and the high milling temperature induces a variety of processing defects.This paper obtained the temperature variation data during the end milling process of CFRP laminates through experiments.After data fitting,the data were transformed into a function of heat flux density varying with time.In the finite element analysis,a double-ellipsoid moving heat source model was introduced,and a moving heat source subrou-tine was written based on the time-varying function of heat flux density to more accurately describe the thermal effects dur-ing the milling process and simulate the changes in the temperature field during milling.The Hashin failure criterion is a-dopted as the basis of fiber and matrix failure,and the simulation results of the temperature field are input into the thermal-force coupling simulation model as the predefined field conditions for solving and analyzing by means of sequential thermal-force coupling,so as to establish a thermal-force coupling simulation and analysis model for milling processing of CFRP end faces.The model simulation results can provide a basis for exploring the damage evolution law of CFRP material under the influence of temperature.
基金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.
基金financially supported by the National Natural Science Foundation of China(Grants nos.62201411,62371378,22205168,52302150 and 62304171)the China Postdoctoral Science Foundation(2022M722500)+1 种基金the Fundamental Research Funds for the Central Universities(Grants nos.ZYTS2308 and 20103237929)Startup Foundation of Xidian University(10251220001).
文摘Defects-rich heterointerfaces integrated with adjustable crystalline phases and atom vacancies,as well as veiled dielectric-responsive character,are instrumental in electromagnetic dissipation.Conventional methods,however,constrain their delicate constructions.Herein,an innovative alternative is proposed:carrageenan-assistant cations-regulated(CACR)strategy,which induces a series of sulfides nanoparticles rooted in situ on the surface of carbon matrix.This unique configuration originates from strategic vacancy formation energy of sulfides and strong sulfides-carbon support interaction,benefiting the delicate construction of defects-rich heterostructures in M_(x)S_(y)/carbon composites(M-CAs).Impressively,these generated sulfur vacancies are firstly found to strengthen electron accumulation/consumption ability at heterointerfaces and,simultaneously,induct local asymmetry of electronic structure to evoke large dipole moment,ultimately leading to polarization coupling,i.e.,defect-type interfacial polarization.Such“Janus effect”(Janus effect means versatility,as in the Greek two-headed Janus)of interfacial sulfur vacancies is intuitively confirmed by both theoretical and experimental investigations for the first time.Consequently,the sulfur vacancies-rich heterostructured Co/Ni-CAs displays broad absorption bandwidth of 6.76 GHz at only 1.8 mm,compared to sulfur vacancies-free CAs without any dielectric response.Harnessing defects-rich heterostructures,this one-pot CACR strategy may steer the design and development of advanced nanomaterials,boosting functionality across diverse application domains beyond electromagnetic response.
基金supported by the National Natural Science Foundation of China(Grant No.52105433)the National Natural Science Foundation of China(Grant No.52175383)the Henan Key Laboratory of Superhard Abrasives and Grinding Equipment,Henan University of Technology,China(JDKFJJ2022006).
文摘Grinding is a generally utilized method,removing excess materials through effective abrasives.The grinding abrasives with multiple shape-position characteristics play a dominant role in determining the thermo-mechanical coupling,which may influence the surface quality directly.To investigate this correlated influence mechanism,this paper focuses on the abrasive shape and position characteristic on the grinding thermo-mechanical process with the analytic single abrasive interaction force model by considering the abrasive shape and its distribution information.It can be found that the mapped dynamic grinding temperature is actually discretized on the workpiece surface,which is on account of the diversity of the abrasive shape and its distribution.Moreover,higher spherical and conical abrasive particles,as well as lower pyramid shaped abrasive particle ratios,can generate greater specific grinding energy with a discretized temperature distribution,when compared with a higher proportion of pyramid shaped abrasive.The study can be utilized to provide valuable theoretical foundation for engineering practice by preparing structural wheel and its grinding property.
基金supported by National Key R&D Program of China(Grant No.2021YFB2601200)Open Fund of National Engineering Research Center of Highway Maintenance Technology(Changsha University of Science&Technology)(No.kfj230207).
文摘The occurrence of top-down(TD)cracking has gradually become a prevalent issue in semi-rigid base asphalt pavements after prolonged service.A coupled simulation model integrating the finite difference method(FDM)and discrete element method(DEM)was employed to investigate the mechanical behavior of asphalt pavement containing a pre-existing TD crack.The mesoscopic parameters of the model were calibrated based on the mixture modulus and the static mechanical response on the MLS66 test road.Finally,an analysis was performed to assess how variations in TD crack depth and longitudinal length affect the distribution patterns of transverse tensile stress,vertical shear stress,and vertical compressive stress.The results indicate that the vertical propagation of TD crack significantly increases both the tensile stress value and range on the middle surface,while the longitudinal development of TD crack has minimal impact.This phenomenon may result in more severe fatigue failure on the middle surface.With the vertical and longitudinal development of TD crack,the vertical shear stress and compressive stress show obvious"two-stage"characteristics.When the crack's vertical length reaches 40 mm,there is a sharp increase in stress on the upper surface.As the crack continues to propagate vertically,the growth of stress on the upper surface becomes negligible,while the stress in the middle and lower layers increased significantly.Conversely,for longitudinal development of TD crack,any changes in stress are insignificant when their length is less than 180 mm;however,as they continue to develop longitudinally beyond this threshold,there is a sharp increase in stress levels.These findings hold great significance for understanding pavement structure deterioration and maintenance behavior associated with TD crack.
基金National Natural Science Foundation of China (U23A6005 and 22078069)Project funded by China Postdoctoral Science Foundation (GZB20230172 and 2023M740748)。
文摘Electrochemical conversion of lignin for the production of high-value heterocyclic aromatic compounds has great potential.We demonstrate the targeted synthesis and cation modulation of NiCo_(2)O_(4)spinel nanoboxes,synthesized via cation exchange and calcination oxidation.These catalysts exhibit excellent efficacy in the electrocatalytic conversion of lignin model compounds,specifically 2-phenoxy-1-phenylethanol,into nitrogen-containing aromatics,achieving high conversion rates and selectivities.These catalysts were synthesized via a cation exchange and calcination oxidation process,using Prussian blue nanocubes as precursors.The porous architecture and polymetallic composition of the NiCo_(2)O_(4)spinel demonstrated superior performance in electrocatalytic oxidative coupling,achieving a 99.2 wt%conversion rate of the 2-phenoxy-1-phenylethanol with selectivities of 37.5 wt%for quinoline derivatives and 31.5 wt%for phenol.Key innovations include the development of a sustainable one-pot synthesis method for quinoline derivatives,the elucidation of a multistage reaction pathway involving CAO bond cleavage,hydroxyaldol condensation,and CAN bond formation,and a deeper mechanistic understanding derived from DFT simulations.This work establishes a new strategy for lignin valorization,offering a sustainable route to produce high-value nitrogen-containing aromatics from renewable biomass under mild conditions,without the need for additional reagents.
基金financially supported by the National Natural Science Foundation of China(Nos.42272153 and 42472195)the Research Fund of PetroChina Tarim Oilfield Company(No.671023060003)the Research Fund of China National Petroleum Corporation Limited(No.2023ZZ16YJ04).
文摘Deep tight reservoirs exhibit complex stress and seepage fields due to varying pore structures,thus the seepage characteristics are significant for enhancing oil production.This study conducted triaxial compression and permeability tests to investigate the mechanical and seepage properties of tight sandstone.A digital core of tight sandstone was built using Computed Tomography(CT)scanning,which was divided into matrix and pore phases by a pore equivalent diameter threshold.A fluid-solid coupling model was established to investigate the seepage characteristics at micro-scale.The results showed that increasing the confining pressure decreased porosity,permeability,and flow velocity,with the pore phase becoming the dominant seepage channel.Cracks and large pores closed first under increasing pressure,resulted in a steep drop in permeability.However,permeability slightly decreased under high confining pressure,which followed a first-order exponential function.Flow velocity increased with seepage pressure.And the damage mainly occurred in stress-concentration regions under low seepage pressure.Seepage behavior followed linear Darcy flow,the damage emerged at seepage entrances under high pressure,which decreased rock elastic modulus and significantly increased permeability.
