Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is sti...Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.展开更多
In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible plat...In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.展开更多
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.展开更多
Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie conditio...Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.展开更多
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 Distributed Propulsion Wing(DPW)presents prominent advantages in terms of energy conservation during flight,but the intense integration of propulsive internal flow with aerodynamic external flow brings significant...The Distributed Propulsion Wing(DPW)presents prominent advantages in terms of energy conservation during flight,but the intense integration of propulsive internal flow with aerodynamic external flow brings significant design challenges.To tackle this issue,this paper undertakes a comprehensive investigation of the aero-propulsive coupling performance of the DPW under both hovering and cruising conditions,and subsequently proposes a multi-level collaboration optimization design method based on the decomposition principle.Specifically,the complex 3D surfaces of DPW are systematically dissociated into simple 2D curves with inherent relationships for design.The decomposition is achieved based on the analysis results of the aero-propulsive coupling characteristics.And a DPW design case is conducted and subsequently analyzed in order to further validate the effectiveness and feasibility of the proposed design method.It is shown that a 115.75%drag reduction of DPW can be achieved at cruise under a specified thrust level.Furthermore,the DPW exhibits inherent characteristics of consistent lift-to-drag ratio with the thrust-drag balance constraint,regardless of variations in incoming flow velocity or total thrust.展开更多
The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computatio...The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.展开更多
A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a...A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.展开更多
Magneto-electro-elastic (MEE) materials, a new type of composite intelligent materials, exhibit excellent multifield coupling effects. Due to the heterogeneity of the materials, it is challenging to use the traditiona...Magneto-electro-elastic (MEE) materials, a new type of composite intelligent materials, exhibit excellent multifield coupling effects. Due to the heterogeneity of the materials, it is challenging to use the traditional finite element method (FEM) for mechanical analysis. Additionally, the MEE materials are often in a complex service environment, especially under the influence of the thermal field with thermoelectric and thermomagnetic effects, which affect its mechanical properties. Therefore, this paper proposes the efficient multiscale computational method for the multifield coupling problem of heterogeneous MEE structures under the thermal environment. The method constructs a multi-physics field with numerical base functions (the displacement, electric potential, and magnetic potential multiscale base functions). It equates a single cell of heterogeneous MEE materials to a macroscopic unit and supplements the macroscopic model with a microscopic model. This allows the problem to be solved directly on a macroscopic scale. Finally, the numerical simulation results demonstrate that compared with the traditional FEM, the multiscale finite element method (MsFEM) can achieve the purpose of ensuring accuracy and reducing the degree of freedom, and significantly improving the calculation efficiency.展开更多
A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the key...A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.展开更多
The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. ...The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.展开更多
A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed cou...A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.展开更多
An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted F...An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.展开更多
One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary condit...One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpola- tion can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin（EFG）-radial point interpola- tion method （RPIM） is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.展开更多
A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(F...A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.展开更多
The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surfa...The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surface conditions to regional climate models( RCMs). There are two methods of downscaling: offline coupling and online coupling. The two kinds of coupling methods are described in detail by coupling the Weather Research and Forecasting model( WRF) with the Institute of Atmospheric Physics of Chinese Academy of Sciences Atmospheric General Circulation Model Version 4. 0( IAP AGCM4. 0) in the study. And the extreme precipitation event over Beijing on July 212012 is simulated by using the two coupling methods. Results show that online coupling method is of great value in improving the model simulation. Furthermore,the data exchange frequency of online coupling has some effect on simulation result.展开更多
To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this pap...To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this paper,a coupled BEM(Boundary Element Method)/RANS(Renolds-Averaged Navier−Stokes)solver is used to simulate propeller behind the hull in the self-propulsion test.The motivation of this work is to develop a practical tool to design marine propulsion system without suffering long computational time.An unsteady boundary element method which is also known as panel method is chosen to estimate the propeller forces.Propeller wakes are treated using a time marching wake alignment method.Also,a RANS code coupled with VoF equation is developed to consider the ship motions and wake field effects in the problem.A coupling algorithm is developed to interchange ship wake field to the potential flow solver and propeller thrust to the RANS code.Based on the difference between hull resistance and the propeller thrust,a PI controller is developed to compute the propeller RPM in every time step.Verification of the solver is carried out using the towing tank test report of a 50 m oceanography research vessel.Wake factor and trust deduction coefficient are estimated numerically.Also,the wake rollup pattern of the propeller in open water is compared with the propeller in real wake field.展开更多
We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved e...We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved exactly,yielding a set of Bessel vortices.Additionally,based on linear solution and using variational approximation,the solutions for the full nonlinear system and their ground state phase diagrams are derived,including the vortex states with quantum numbers m=0,1,as well as mixed states.In this work,mixed states in spin-1 spin-orbit coupling(SOC)BEC are interpreted for the first time as weighted superpositions of three vortex states.Furthermore,the results also indicate that under strong Zeeman splitting,the system cannot form localized states.The variational solutions align well with numerical simulations,showing stable evolution and meeting the criteria for long-term observation in experiments.展开更多
A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radi...A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radial and axial magnetic pressure acting on the tube, the magnetic pressure is then used as boundary conditions to model the high velocity deformation of tube with DYNAFORM, The radial magnetic pressure on the tube decreases from the center to the tube end, axial magnetic pressure is greater near the location equal to the coil height and slight in the other region. The radial displacement of deformed workpicces is distributed uniformly near the tube center and decreases from the center to the end; Deformation from the location equal to coil height to the tube end is little. This distribution is consistent with the distribution of radial pressure; Effect of the axial magnetic pressure on deformation can be ignored, The calculated results show well agreements with the experimental results.展开更多
Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the be...Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the best catalytic performance among all samples. In addition, the effects of different addition sequences of Na, W, Mn, Zr, S and P on the catalytic performance were studied. The absence of Na before the addition of Mn and Zr in the catalysts preparation depressed the formation of the active phases of Mn2O3 and ZrO2 and decreased the activities of the catalysts significantly.展开更多
基金financially supported by the National Key Research and Development Plan(2023YFC2811001)the National Natural Science Foundation of China(42206233)the Taishan Scholars Program(tsqn202312280,tsqn202306297)。
文摘Hydraulic fracturing technology has played an important role in the exploitation of unconventional oil and gas resources,however,its application to gas hydrate reservoirs has been rarely studied.Currently,there is still limited understanding of the propagation and extension of fractures around the wellbore during the fracturing process of horizontal wells in hydrate reservoirs,as well as the stress interference patterns between fractures.This study simulates hydraulic fracturing processes in hydrate reservoirs using a fluidsolid coupling discrete element method(DEM),and analyzes the impacts of hydrate saturation and geological and engineering factors on fracture extension and stress disturbance.The results show that hydraulic fracturing is more effective when hydrate saturation exceeds 30%and that fracture pressure increases with saturation.The increase in horizontal stress differential enhances the directionality of fracture propagation and reduces stress disturbance.The distribution uniformity index(DUI)reveals that injection pressure is directly proportional to the number of main fractures and inversely proportional to fracturing time,with fracturing efficiency depending on the spacing between injection points and the distance between wells.This work may provide reference for the commercial exploitation of natural gas hydrates.
基金the National Natural Science Foundation of China(No.U20A20328).
文摘In this paper,the multi-body coupled dynamic characteristics of a semisubmersible platform and an HYSY 229 barge were investigated.First,coupled hydrodynamic analysis of the HYSY 229 barge and the semisubmersible platform was performed.Relevant hydrodynamic parameters were obtained using the retardation function method of three-dimensional frequency-domain potential flow theory.The results of the hydrodynamic analysis were highly consistent with the test findings,verifying the accuracy of the multifloating hydrodynamic coupling analysis,and key hydrodynamic parameters were solved for different water depths and the coupling effect.According to the obtained results,the hydrodynamic influence was the largest in shallow waters when the coupling effect was considered.Furthermore,the coupled motion equation combined with viscous damping,fender system,and mooring system was established,and the hydrodynamics,floating body motion,and dynamic response of the fender system were analyzed.Motion analysis revealed good agreement among the surge,sway,and yaw motions of the two floating bodies.However,when the wave period reached 10 s,the motion of the two floating bodies showed severe shock,and a relative motion was also observed.Therefore,excessive constraints should be added between the two floating bodies during construction to ensure construction safety.The numerical analysis and model test results of the semisubmersible platform and HYSY 229 barge at a water depth of 42 m and sea conditions of 0°,45°,and 90° were in good agreement,and the error was less than 5%.The maximum movement of the HYSY 229 barge reached 2.61 m in the sway direction,whereas that of the semisubmersible platform was 2.11 m.During construction,excessive constraints should be added between the two floating bodies to limit their relative movement and ensure construction safety.
文摘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.
基金a U.S. Federal Railroad Administration (FRA)BAA project,titled “Mitigation of Differential Movement at Railway Transitions for High-Speed Passenger Rail and Joint Passenger/Freight Corridors”the financial support provided by the China Scholarship Council (CSC),which funded Zhongyi Liu’s and Wenjing Li’s time and research efforts for this study
文摘Nonuniform track support and differential settlements are commonly observed in bridge approaches where the ballast layer can develop gaps at crosstie-ballast interfaces often referred to as a hanging crosstie condition.Hanging crossties usually yield unfavorable dynamic effects such as higher wheel loads,which negatively impact the serviceability and safety of railway operations.Hence,a better understanding of the mechanisms that cause hanging crossties and their effects on the ballast layer load-deformation characteristics is necessary.Since the ballast layer is a particulate medium,the discrete element method(DEM),which simulates ballast particle interactions individually,is ideal to explore the interparticle contact forces and ballast movements under dynamic wheel loading.Accurate representations of the dynamic loads from the train and track superstructure are needed for high-fidelity DEM modeling.This paper introduces an integrated modeling approach,which couples a single-crosstie DEM ballast model with a train–track–bridge(TTB)model using a proportional–integral–derivative control loop.The TTB–DEM model was validated with field measurements,and the coupled model calculates similar crosstie displacements as the TTB model.The TTB–DEM provided new insights into the ballast particle-scale behavior,which the TTB model alone cannot explore.The TTB–DEM coupling approach identified detrimental effects of hanging crossties on adjacent crossties,which were found to experience drastic vibrations and large ballast contact force concentrations.
基金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.
基金co-supported by the Equipment Advance Research Project of China(No.50911040803)the National Defense Pre-research Foundation of China(No.2021-JCJQJJ-0805)the Aeronautical Science Foundation of China(No.2024Z006053001)。
文摘The Distributed Propulsion Wing(DPW)presents prominent advantages in terms of energy conservation during flight,but the intense integration of propulsive internal flow with aerodynamic external flow brings significant design challenges.To tackle this issue,this paper undertakes a comprehensive investigation of the aero-propulsive coupling performance of the DPW under both hovering and cruising conditions,and subsequently proposes a multi-level collaboration optimization design method based on the decomposition principle.Specifically,the complex 3D surfaces of DPW are systematically dissociated into simple 2D curves with inherent relationships for design.The decomposition is achieved based on the analysis results of the aero-propulsive coupling characteristics.And a DPW design case is conducted and subsequently analyzed in order to further validate the effectiveness and feasibility of the proposed design method.It is shown that a 115.75%drag reduction of DPW can be achieved at cruise under a specified thrust level.Furthermore,the DPW exhibits inherent characteristics of consistent lift-to-drag ratio with the thrust-drag balance constraint,regardless of variations in incoming flow velocity or total thrust.
基金supported by Shanghai 2021“Science and Technology Innovation Action Plan”:Social Development Science and Technology Research Project(Grant No.21DZ1202703).
文摘The numerical simulation of the fluid flow and the flexible rod(s)interaction is more complicated and has lower efficiency due to the high computational cost.In this paper,a semi-resolved model coupling the computational fluid dynamics and the flexible rod dynamics is proposed using a two-way domain expansion method.The gov-erning equations of the flexible rod dynamics are discretized and solved by the finite element method,and the fluid flow is simulated by the finite volume method.The interaction between fluids and solid rods is modeled by introducing body force terms into the momentum equations.Referred to the traditional semi-resolved numerical model,an anisotropic Gaussian kernel function method is proposed to specify the interactive forces between flu-ids and solid bodies for non-circle rod cross-sections.A benchmark of the flow passing around a single flexible plate with a rectangular cross-section is used to validate the algorithm.Focused on the engineering applications,a test case of a finite patch of cylinders is implemented to validate the accuracy and efficiency of the coupled model.
基金supported by Scientific Research Fund of Institute of Engineering Mechanics,China Earthquake Administration(Grant Nos.2021B06,2021C05)Heilongjiang Natural Science Foundation Joint Guidance Project(Grant No.LH2021E122).
文摘A coupled numerical calculation method combining smooth particle hydrodynamics(SPH)and the finite element method(FEM)was implemented to investigate the seismic response of horizontal storage tanks.Anumericalmodel of a horizontal storage tank featuring a free liquid surface under seismic action was constructed using the SPH–FEM coupling method.The stored liquid was discretized using SPH particles,while the tank and supports were discretized using the FEM.The interaction between the stored liquid and the tank was simulated by using the meshless particle contact method.Then,the numerical simulation results were compared and analyzed against seismic simulation shaking table test data to validate the method.Subsequently,a series of numerical models,considering different liquid storage volumes and seismic effects,were constructed to obtain time history data of base shear and top center displacement,which revealed the seismic performance of horizontal storage tanks.Numerical simulation results and experimental data showed good agreement,with an error rate of less than 18.85%.And this conformity signifies the rationality of the SPH-FEM coupling method.The base shear and top center displacement values obtained by the coupled SPH-FEM method were only 53.3% to 69.1% of those calculated by the equivalent mass method employed in the current code.As the stored liquid volume increased,the seismic response of the horizontal storage tank exhibited a gradual upward trend,with the seismic response increasing from 73% to 388% for every 35% increase in stored liquid volume.The maximum von Mises stress of the tank and the supports remained below the steel yield strength during the earthquake.The coupled SPH-FEM method holds certain advantages in studying the seismic problems of tanks with complex structural forms,particularly due to the representation of the flow field distribution during earthquakes by involving reservoir fluid participation.
文摘Magneto-electro-elastic (MEE) materials, a new type of composite intelligent materials, exhibit excellent multifield coupling effects. Due to the heterogeneity of the materials, it is challenging to use the traditional finite element method (FEM) for mechanical analysis. Additionally, the MEE materials are often in a complex service environment, especially under the influence of the thermal field with thermoelectric and thermomagnetic effects, which affect its mechanical properties. Therefore, this paper proposes the efficient multiscale computational method for the multifield coupling problem of heterogeneous MEE structures under the thermal environment. The method constructs a multi-physics field with numerical base functions (the displacement, electric potential, and magnetic potential multiscale base functions). It equates a single cell of heterogeneous MEE materials to a macroscopic unit and supplements the macroscopic model with a microscopic model. This allows the problem to be solved directly on a macroscopic scale. Finally, the numerical simulation results demonstrate that compared with the traditional FEM, the multiscale finite element method (MsFEM) can achieve the purpose of ensuring accuracy and reducing the degree of freedom, and significantly improving the calculation efficiency.
基金Projects (51175162, 50805045) supported by the National Natural Science Foundation of ChinaProject supported by the Scientific Research Foundation for the Returned Overseas Chinese Scholars,Ministry of Education,China
文摘A mathematical energy coupling model was developed to analyze the light transmission in the keyhole and energy distribution on the keyhole wall.The main characteristics of the model include：1） a prototype of the keyhole and the inverse Bremsstrahlung absorption coefficient in the keyhole plasma are obtained from the experiments;2） instead of using a parallel incident beam,a focused laser beam with real Gaussian intensity distribution is implemented;3） both Fresnel absorption and inverse Bremsstrahlung absorption during multiple reflections are considered.The calculation results show that the distribution of absorbed laser intensity by the keyhole wall is not uniform.The maximum laser energy is absorbed by the bottom of the keyhole,although no rays irradiate directly onto the bottom.According to analysis of beam focusing characteristics,the location of the focal plane plays a more important role in the laser energy absorption by the front wall than by the rear wall.
基金Sponsored by National Science and Technology Support Plan of China (2009AA04Z143)Science and Technology Support Plan of Hebei Province of China (10212101D)Important Natural Science Foundation of Hebei Province of China (E2006001038)
文摘The shape of strip is calculated by iterative method which combines strip plastic deformation model with rolls elastic deformation model through their calculation results, which can be called results coupling method. Be- cause the shape and rolling force distribution are very sensitive to strip thickness transverse distribution% variation, the iterative course is rather unstable and sometimes convergence cannot be achieved. In addition, the calculating speed of results coupling method is low, which restricts its usable range. To solve the problem, a new model cou- pling method is developed, which takes the force distribution between rolls, rolling force distribution and strip＇s exit transverse displacement distribution as basic unknowns, and integrates strip plastic deformation model and rolls elas- tic deformation model as a unified linear equations through their internal relation, so the iterative calculation between the strip plastic deformation model and rolls elastic deformation model can be avoided. To prove the effectiveness of the model coupling method, two examples are calculated by results coupling method and model coupling method re- spectively. The results of front tension stress, back tension stress, strip^s exit gauge, the force between rolls and rolling force distribution calculated by model coupling method coincide very well with results coupling method. How- ever the calculation course of model coupling method is more steady than results coupling method, and its calculating speed is about ten times as much as the maximal speed of results coupling method, which validates its practicability and reliability.
基金Projects(51474251,51874351) supported by the National Natural Science Foundation of China
文摘A new electrical method of conductive carbon-film(with waterproof and anticorrosion ability)was proposed to continuously measure crack propagation rate of brittle rock under THMC coupling condition.A self-designed coupling testing system was used to conduct THMC coupling fracture tests of the pre-cracked red sandstone specimens(where the temperature is only changed)by this new electrical method of conductive carbon-film.Calculation results obtained by the energy method coincide well with the test results.And the higher the temperature is,the earlier the crack is initiated and the larger the crack propagation rate and accelerated velocity are,which can prove the validity of the new electrical method.This new electrical method has advantages of continuously measuring crack propagation rate over the conventional electrical,optical and acoustic methods,and can provide important basis for safety assessment and cracking-arrest design of deep rock mass engineering.
基金Project (No. 50575143) supported by the National Natural Science Foundation of China
文摘An adaptive finite element-element-free Galerkin (FE-EFG) coupling method is proposed and developed for the numerical simulation of bulk metal forming processes. This approach is able to adaptively convert distorted FE elements to EFG domain in analysis. A new scheme to implement adaptive conversion and coupling is presented. The coupling method takes both advantages of finite element method (FEM) and meshless methods. It is capable of handling large deformations with no need of remeshing procedures, while it is computationally more efficient than those full meshless methods. The effectiveness of the proposed method is demonstrated with the numerical simulations of the bulk metal forming processes including forging and extrusion.
基金supported by the National Natural Science Foundation of China (No. 11172192)the College Postgraduate Research and Innovation Project of Jiangsu Province (No. CX10B 029Z)the Nominated Excellent Thesis for PHD Candidates Program of Soochow University (No. 23320957)
文摘One of major difficulties in the implementation of meshfree methods using the mov- ing least square （MLS） approximation, such as element-free Galerkin method （EFG）, is the im- position of essential boundary conditions as the approximations do not pass through the nodal parameter values. Another class of meshfree methods based on the radial basis point interpola- tion can satisfy the essential boundary conditions exactly since its approximation function passes through each node in an influence domain and thus its shape functions possess the properties of delta function. In this paper, a coupled element-free Galerkin（EFG）-radial point interpola- tion method （RPIM） is proposed to enhance their advantages and avoid their disadvantages. Discretized equations of equilibrium are obtained in the RPIM region and the EFG region, respectively. Then a collocation approach is introduced to couple the RPIM and the EFG method. This method satisfies the linear consistency exactly and can maintain the stiffness matrix symmetric. Numerical tests show that this method gives reasonably accurate results consistent with the theory.
基金co-supported by the National Key R&D Program of China(Nos.2018YFF01012401-05)the National Natural Science Foundation of China(No.51975243)+2 种基金Jilin Provincial Department of Education(No.JJKH20180084KJ),Chinathe Fundamental Research Funds for the Central Universities and Jilin Provincial Department of Science&Technology Fund Project,China(Nos.20170101043JC and 20180520072JH)Graduate Innovation Fund of Jilin University,China(No.101832018C184).
文摘A Coupling Magneto-Electro-Elastic(MEE)Node-based Smoothed Radial Point Interpolation Method(CM-NS-RPIM)was proposed to solve the free vibration and transient responses of Functionally Graded Magneto-Electro-Elastic(FGMEE)structures.By introducing the modified Newmark method,the displacement,electrical potential and magnetic potential of the structures under transient mechanical loading were obtained.Based on G space theory and the weakened weak(W2)formulation,the equations of the multi-physics coupling problems were derived.Using triangular background elements,the free vibration and transient responses of three numerical examples were studied.Results proved that CM-NS-RPIM performed better than the standard FEM by reducing the overly-stiff of structures.Moreover,CM-NS-RPIM could reduce the number of nodes while guaranteeing the accuracy.Besides,triangular elements could be generated automatically even for complex geometries.Therefore,the effectiveness and validity of CM-NS-RPIM were demonstrated,which were valuable for the design of intelligence devices,such as energy harvesters and sensors.
基金Supported by the National Natural Science Foundation of China(No.61602477)China Postdoctoral Science Foundation(No.2016M601158)National Key Research and Development Program of China(No.2016YFB0200804)
文摘The future climate dynamical downscaling method is that output of general circulation models( GCMs) is employed to provide initial conditions,lateral boundary conditions,sea surface temperatures,and initial land surface conditions to regional climate models( RCMs). There are two methods of downscaling: offline coupling and online coupling. The two kinds of coupling methods are described in detail by coupling the Weather Research and Forecasting model( WRF) with the Institute of Atmospheric Physics of Chinese Academy of Sciences Atmospheric General Circulation Model Version 4. 0( IAP AGCM4. 0) in the study. And the extreme precipitation event over Beijing on July 212012 is simulated by using the two coupling methods. Results show that online coupling method is of great value in improving the model simulation. Furthermore,the data exchange frequency of online coupling has some effect on simulation result.
文摘To design a propeller for ship power plant,the interaction between ship hull and propeller must be taken into account.The main concern is to apply the wake effect of ship stern on the propeller performance.In this paper,a coupled BEM(Boundary Element Method)/RANS(Renolds-Averaged Navier−Stokes)solver is used to simulate propeller behind the hull in the self-propulsion test.The motivation of this work is to develop a practical tool to design marine propulsion system without suffering long computational time.An unsteady boundary element method which is also known as panel method is chosen to estimate the propeller forces.Propeller wakes are treated using a time marching wake alignment method.Also,a RANS code coupled with VoF equation is developed to consider the ship motions and wake field effects in the problem.A coupling algorithm is developed to interchange ship wake field to the potential flow solver and propeller thrust to the RANS code.Based on the difference between hull resistance and the propeller thrust,a PI controller is developed to compute the propeller RPM in every time step.Verification of the solver is carried out using the towing tank test report of a 50 m oceanography research vessel.Wake factor and trust deduction coefficient are estimated numerically.Also,the wake rollup pattern of the propeller in open water is compared with the propeller in real wake field.
基金supported by the Guangdong Basic and Applied Basic Research Foundation(Grant No.2023A1515110198)the Natural Science Foundation of Guangdong Province,China(Grant Nos.2024A1515030131 and 2021A1515010214)+2 种基金the National Natural Science Foundation of China(Grant Nos.12274077,11905032,and 12475014)the Research Fund of the Guangdong-Hong Kong-Macao Joint Laboratory for Intelligent Micro-Nano Optoelectronic Technology(Grant No.2020B1212030010)the Israel Science Foundation(Grant No.1695/22).
文摘We investigate the ground states of spin-orbit coupled spin-1 Bose-Einstein condensates in the presence of Zeeman splitting.By introducing the generalized momentum operator,the linear version of the system is solved exactly,yielding a set of Bessel vortices.Additionally,based on linear solution and using variational approximation,the solutions for the full nonlinear system and their ground state phase diagrams are derived,including the vortex states with quantum numbers m=0,1,as well as mixed states.In this work,mixed states in spin-1 spin-orbit coupling(SOC)BEC are interpreted for the first time as weighted superpositions of three vortex states.Furthermore,the results also indicate that under strong Zeeman splitting,the system cannot form localized states.The variational solutions align well with numerical simulations,showing stable evolution and meeting the criteria for long-term observation in experiments.
文摘A loose coupling method is used to solve the electromagnetic tube bulging. ANSYS/ EMAG is used to model the time varying electromagnetic field with the discharge current used as excitation, in order to obtain the radial and axial magnetic pressure acting on the tube, the magnetic pressure is then used as boundary conditions to model the high velocity deformation of tube with DYNAFORM, The radial magnetic pressure on the tube decreases from the center to the tube end, axial magnetic pressure is greater near the location equal to the coil height and slight in the other region. The radial displacement of deformed workpicces is distributed uniformly near the tube center and decreases from the center to the end; Deformation from the location equal to coil height to the tube end is little. This distribution is consistent with the distribution of radial pressure; Effect of the axial magnetic pressure on deformation can be ignored, The calculated results show well agreements with the experimental results.
基金supported by the financial support from National Natural Science Foundation of China (20676116)
文摘Na-W-Mn-Zr-S-P/SiO2 catalysts for oxidative coupling of methane （OCM） were prepared by incipient wetness impregnation, sol-gel and mixture slurry methods. The catalyst prepared by mixture slurry method showed the best catalytic performance among all samples. In addition, the effects of different addition sequences of Na, W, Mn, Zr, S and P on the catalytic performance were studied. The absence of Na before the addition of Mn and Zr in the catalysts preparation depressed the formation of the active phases of Mn2O3 and ZrO2 and decreased the activities of the catalysts significantly.
