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 arresting process of carrier-based aircraft is widely recognized as a challenging task,characterized by the highest accident rate among all carrier-based aircraft operations.Dynamic simulation plays a crucial role...The arresting process of carrier-based aircraft is widely recognized as a challenging task,characterized by the highest accident rate among all carrier-based aircraft operations.Dynamic simulation plays a crucial role in assessing the intricate responses of the arresting process,favoring the design of carrier-based aircraft.An efficient and accurate rigid-flexible coupling model for analyzing the dynamic response of the arresting process is proposed.By combining the dynamic characteristics of airframe,landing gear,arresting hook and arresting gear system,the rigid-flexible coupling dynamic model is established to reflect the relative motion of the coupling parts and arresting load.The dynamic model is verified through simulations of landing gear landing drops and by comparing the arresting simulation results with corresponding data in the US military standard.Additionally,simulations of the arresting process under different off-center distance and aircraft yaw angle are conducted to obtain the dynamic response of the aircraft during the arresting process.The result indicates that the rigid-flexible coupling dynamic model proposed is effective for analyzing the arresting dynamics response of carrier-based aircraft.The axial force of the arresting cable on both sides of the hook engagement point,pitch and yaw angle of aircraft are inconsistent under yaw and off-center arresting.The analysis method and obtained results provide valuable references for assessing the dynamic responses of carrier-based aircraft during arresting process and offer valuable in-sights in the design of carrier-based aircraft.展开更多
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.展开更多
Considering the dynamic variation of roll gap and the transverse distribution of dynamic rolling force along the work roll width direction, the movement and deformation of rolls system, influenced by the coupling of v...Considering the dynamic variation of roll gap and the transverse distribution of dynamic rolling force along the work roll width direction, the movement and deformation of rolls system, influenced by the coupling of vertical chatter and transverse bending vibration, may cause instability and also bring product defect of thickness difference. Therefore, a rigid-flexible coupling vibration model of the rolls system was presented. The influence of dynamic characteristics on the rolling process stability and strip thickness distribution was investigated. Firstly, assuming the symmetry of upper and lower structures of six-high rolling mill, a transverse bending vibration model of three-beam system under simply supported boundary conditions was established, and a semi-analytical solution method was proposed to deal with this model. Then, considering both variation and change rate of the roll gap, a roll vertical chatter model with structure and process coupled was constructed, and the critical rolling speed for self-excited instability was determined by Routh stability criterion. Furthermore, a rigid-flexible coupling vibration model of the rolls system was built by connecting the vertical chatter model and transverse bending vibration model through the distribution of dynamic rolling force, and the dynamic characteristics of rolls system were analyzed. Finally, the strip exit thickness distributions under the stable and unstable rolling process were compared, and the product shape and thickness distribution characteristics were quantitatively evaluated by the crown and maximum longitudinal thickness difference.展开更多
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.展开更多
Robot-automated spraying is widely used in various fields,such as the automotive,metalworking,furniture,and aero-space industries.Spraying quality is influenced by multiple factors,including robot speed,acceleration,e...Robot-automated spraying is widely used in various fields,such as the automotive,metalworking,furniture,and aero-space industries.Spraying quality is influenced by multiple factors,including robot speed,acceleration,end-effector trajectory,and spraying process constraints.To achieve high-quality spraying under the influence of multiple factors,this study proposes a multi-objective optimization method for the spraying trajectory that integrates spraying process constraints into the optimization process.First,a 7-degree-of-freedom rigid-flexible coupling serial spray painting robot system is introduced,which includes a motion decoupling mechanism and a tension amplification mechanism.Subsequently,a paint deposition model for the spray gun was established,and the influence of process constraints on spraying quality was analyzed.Trajectory planning for the spray painting robot,based on the septic B-spline interpolation method,was then performed.Based on this foundation,objective functions and constraint equations for spraying trajectory optimization were established.A multi-objective trajectory optimization method for spraying by the robot is proposed based on the NSGA-Ⅱ,which integrates the spraying process constraints.Finally,a prototype system of a 7-degree-of-freedom rigid-flexible coupling serial spray painting robot was constructed.Simulations and spraying experiments were conducted to verify the effectiveness of the proposed multi-objective trajectory optimization method.This paper presents a multi-objective optimization method for the spraying trajectory of a robot.In the proposed method,the optimized spraying trajectory is generated with the spraying process as the constraint and time,energy consumption,and impact during the spraying operation of the robot as the optimization objectives.展开更多
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.展开更多
Welded Turnout on Large-span Bridge(WTLB)is a complex multi-layer heterogeneous system and can significantly influence the service performance of High-Speed Railway(HSR).Understanding the coupling dynamic response of ...Welded Turnout on Large-span Bridge(WTLB)is a complex multi-layer heterogeneous system and can significantly influence the service performance of High-Speed Railway(HSR).Understanding the coupling dynamic response of the vehicle and WTLB is essential.Previous research did not consider the dynamic behavior of foundations,leading to an underestimation of the vehicle-turnout-foundation coupling dynamic response,particularly when turnouts were laid on large-span bridges.This study proposes a novel modeling method that includes the foundations,to overcome the previous shortcomings by applying a rigid-flexible coupling system.In this approach,the vehicle was modeled as a rigid body sub-model in a Multi-Body Software(MBS),while WTLB was modeled as a flexible bodies sub-model using Finite Element(FE)software.The modal information from the FE model was imported into the MBS software.The two sub-models were coupled by the wheel-rail contact in the MBS environment and then the Vehicle-turnout-bridge Rigid-flexible Coupling Dynamic(VRCD)calculation model was established and it was discovered that the calculation results showed good agreement with the field test data.Through the VRCD model,the safety of the structure,the stability of the vehicle and the comfort of passengers were investigated,as well as several important infrastructure factors.The results demonstrate that this novel method provides accurate calculations and highlights the complex and significant interactions in the vehicle-turnout-bridge system.展开更多
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.展开更多
The aerodynamic design of a rigid-flexible coupling profile is the decisive factor for the flow-field quality of a supersonic free jet wind tunnel nozzle, and its mechanic dynamic features are the key for engineering ...The aerodynamic design of a rigid-flexible coupling profile is the decisive factor for the flow-field quality of a supersonic free jet wind tunnel nozzle, and its mechanic dynamic features are the key for engineering implementation of continuous Mach number regulations. To fulfill the requirements of a free jet inlet/engine compatibility test within a wide simulation envelop, both uniform flow-fields of continuous acceleration and deceleration are necessary. In this paper, the aerodynamic design methods of an expansion wall and machinery implementation plan for the halfflexible single jack nozzle were researched. The profile control in nozzle flexible plate design was studied with a rigid-flexible coupling method. Design and calculations were performed with the help of numerical simulation. The technique of axial free stretching of the flexible plate was used to improve the matching performance between the designed elasticity profile and the theoretical one, and the rigid-flexible coupling structure was calibrated by wind tunnel tests. Results indicate that the flexible plate aerodynamic design method used here is effective and feasible. Via rigidflexible coupling design, the flexible plate agrees with the rigid body very well, and continuous Mach number changes can be achieved during the tests. The nozzle’s exit flow-field uniformity meets the requirements of China Military Standard(GJB).展开更多
In the present research two different whole vehicle multibody models are established respectively, including rigid and rigid-flexible coupling multibody vehicle models. The former is all composed by rigid bodies while...In the present research two different whole vehicle multibody models are established respectively, including rigid and rigid-flexible coupling multibody vehicle models. The former is all composed by rigid bodies while in the later model, the flexible rear suspension is built based on the finite element method （FEM） and mode superposition method, in which the deformations of the components are considered. The ride simulations with different speeds are carried out on a 3D digitalized road, and the weighted root mean square （RMS） of accelerations on the seat surface,backrest and at the feet are calculated. The comparison between the responses of the rigid and rigid-flexible coupling multibody models shows that the flexibility of the vehicle parts significantly affects the accelerations at each position, and it is necessary to take the flexibility effects into account for the assessment of ride comfort. C 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi： 10.1063/2.1301304]展开更多
A complete geometric nonlinear formulation for rigid-flexible coupling dynamics of a flexible beam undergoing large overall motion was proposed based on virtual work principle, in which all the high-order terms relate...A complete geometric nonlinear formulation for rigid-flexible coupling dynamics of a flexible beam undergoing large overall motion was proposed based on virtual work principle, in which all the high-order terms related to coupling deformation were included in dynamic equations. Simulation examples of the flexible beam with prescribed rotation and free rotation were investigated. Numerical results show that the use of the first-order approximation coupling (FOAC) model may lead to a significant error when the flexible beam experiences large deformation or large deformation velocity. However, the correct solutions can always be obtained by using the present complete model. The difference in essence between this model and the FOAC model is revealed. These coupling high-order terms, which are ignored in FOAC model, have a remarkable effect on the dynamic behavior of the flexible body. Therefore, these terms should be included for the rigid-flexible dynamic modeling and analysis of flexible body undergoing motions with high speed.展开更多
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.展开更多
基金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.
基金This study was co-supported by the National Natural Science Foundation of China(No.T2288101)the National Key Research and Development Project,China(No.2020YFC1512500).
文摘The arresting process of carrier-based aircraft is widely recognized as a challenging task,characterized by the highest accident rate among all carrier-based aircraft operations.Dynamic simulation plays a crucial role in assessing the intricate responses of the arresting process,favoring the design of carrier-based aircraft.An efficient and accurate rigid-flexible coupling model for analyzing the dynamic response of the arresting process is proposed.By combining the dynamic characteristics of airframe,landing gear,arresting hook and arresting gear system,the rigid-flexible coupling dynamic model is established to reflect the relative motion of the coupling parts and arresting load.The dynamic model is verified through simulations of landing gear landing drops and by comparing the arresting simulation results with corresponding data in the US military standard.Additionally,simulations of the arresting process under different off-center distance and aircraft yaw angle are conducted to obtain the dynamic response of the aircraft during the arresting process.The result indicates that the rigid-flexible coupling dynamic model proposed is effective for analyzing the arresting dynamics response of carrier-based aircraft.The axial force of the arresting cable on both sides of the hook engagement point,pitch and yaw angle of aircraft are inconsistent under yaw and off-center arresting.The analysis method and obtained results provide valuable references for assessing the dynamic responses of carrier-based aircraft during arresting process and offer valuable in-sights in the design of carrier-based aircraft.
基金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 the National Natural Science Foundation of China(No.51775038).
文摘Considering the dynamic variation of roll gap and the transverse distribution of dynamic rolling force along the work roll width direction, the movement and deformation of rolls system, influenced by the coupling of vertical chatter and transverse bending vibration, may cause instability and also bring product defect of thickness difference. Therefore, a rigid-flexible coupling vibration model of the rolls system was presented. The influence of dynamic characteristics on the rolling process stability and strip thickness distribution was investigated. Firstly, assuming the symmetry of upper and lower structures of six-high rolling mill, a transverse bending vibration model of three-beam system under simply supported boundary conditions was established, and a semi-analytical solution method was proposed to deal with this model. Then, considering both variation and change rate of the roll gap, a roll vertical chatter model with structure and process coupled was constructed, and the critical rolling speed for self-excited instability was determined by Routh stability criterion. Furthermore, a rigid-flexible coupling vibration model of the rolls system was built by connecting the vertical chatter model and transverse bending vibration model through the distribution of dynamic rolling force, and the dynamic characteristics of rolls system were analyzed. Finally, the strip exit thickness distributions under the stable and unstable rolling process were compared, and the product shape and thickness distribution characteristics were quantitatively evaluated by the crown and maximum longitudinal thickness difference.
基金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.
基金Supported by National Natural Science Foundation of China(Grant Nos.52335002,51925502)。
文摘Robot-automated spraying is widely used in various fields,such as the automotive,metalworking,furniture,and aero-space industries.Spraying quality is influenced by multiple factors,including robot speed,acceleration,end-effector trajectory,and spraying process constraints.To achieve high-quality spraying under the influence of multiple factors,this study proposes a multi-objective optimization method for the spraying trajectory that integrates spraying process constraints into the optimization process.First,a 7-degree-of-freedom rigid-flexible coupling serial spray painting robot system is introduced,which includes a motion decoupling mechanism and a tension amplification mechanism.Subsequently,a paint deposition model for the spray gun was established,and the influence of process constraints on spraying quality was analyzed.Trajectory planning for the spray painting robot,based on the septic B-spline interpolation method,was then performed.Based on this foundation,objective functions and constraint equations for spraying trajectory optimization were established.A multi-objective trajectory optimization method for spraying by the robot is proposed based on the NSGA-Ⅱ,which integrates the spraying process constraints.Finally,a prototype system of a 7-degree-of-freedom rigid-flexible coupling serial spray painting robot was constructed.Simulations and spraying experiments were conducted to verify the effectiveness of the proposed multi-objective trajectory optimization method.This paper presents a multi-objective optimization method for the spraying trajectory of a robot.In the proposed method,the optimized spraying trajectory is generated with the spraying process as the constraint and time,energy consumption,and impact during the spraying operation of the robot as the optimization objectives.
文摘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.
基金supported by the National Natural Science Foundation of China(U23A20666)the China National Railway Group Corporation Science and Technology Research and Development Program(N2023G083).
文摘Welded Turnout on Large-span Bridge(WTLB)is a complex multi-layer heterogeneous system and can significantly influence the service performance of High-Speed Railway(HSR).Understanding the coupling dynamic response of the vehicle and WTLB is essential.Previous research did not consider the dynamic behavior of foundations,leading to an underestimation of the vehicle-turnout-foundation coupling dynamic response,particularly when turnouts were laid on large-span bridges.This study proposes a novel modeling method that includes the foundations,to overcome the previous shortcomings by applying a rigid-flexible coupling system.In this approach,the vehicle was modeled as a rigid body sub-model in a Multi-Body Software(MBS),while WTLB was modeled as a flexible bodies sub-model using Finite Element(FE)software.The modal information from the FE model was imported into the MBS software.The two sub-models were coupled by the wheel-rail contact in the MBS environment and then the Vehicle-turnout-bridge Rigid-flexible Coupling Dynamic(VRCD)calculation model was established and it was discovered that the calculation results showed good agreement with the field test data.Through the VRCD model,the safety of the structure,the stability of the vehicle and the comfort of passengers were investigated,as well as several important infrastructure factors.The results demonstrate that this novel method provides accurate calculations and highlights the complex and significant interactions in the vehicle-turnout-bridge system.
基金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.
基金supported by the National Natural Science Foundation of China (Nos. 90916023 and 51176075)
文摘The aerodynamic design of a rigid-flexible coupling profile is the decisive factor for the flow-field quality of a supersonic free jet wind tunnel nozzle, and its mechanic dynamic features are the key for engineering implementation of continuous Mach number regulations. To fulfill the requirements of a free jet inlet/engine compatibility test within a wide simulation envelop, both uniform flow-fields of continuous acceleration and deceleration are necessary. In this paper, the aerodynamic design methods of an expansion wall and machinery implementation plan for the halfflexible single jack nozzle were researched. The profile control in nozzle flexible plate design was studied with a rigid-flexible coupling method. Design and calculations were performed with the help of numerical simulation. The technique of axial free stretching of the flexible plate was used to improve the matching performance between the designed elasticity profile and the theoretical one, and the rigid-flexible coupling structure was calibrated by wind tunnel tests. Results indicate that the flexible plate aerodynamic design method used here is effective and feasible. Via rigidflexible coupling design, the flexible plate agrees with the rigid body very well, and continuous Mach number changes can be achieved during the tests. The nozzle’s exit flow-field uniformity meets the requirements of China Military Standard(GJB).
基金supported by the National Natural Science Foundation of China(51175379)the National Basic Research Program of China(2011CB711200)
文摘In the present research two different whole vehicle multibody models are established respectively, including rigid and rigid-flexible coupling multibody vehicle models. The former is all composed by rigid bodies while in the later model, the flexible rear suspension is built based on the finite element method （FEM） and mode superposition method, in which the deformations of the components are considered. The ride simulations with different speeds are carried out on a 3D digitalized road, and the weighted root mean square （RMS） of accelerations on the seat surface,backrest and at the feet are calculated. The comparison between the responses of the rigid and rigid-flexible coupling multibody models shows that the flexibility of the vehicle parts significantly affects the accelerations at each position, and it is necessary to take the flexibility effects into account for the assessment of ride comfort. C 2013 The Chinese Society of Theoretical and Applied Mechanics. [doi： 10.1063/2.1301304]
基金Project(10772113) supported by the National Natural Science Foundation of China
文摘A complete geometric nonlinear formulation for rigid-flexible coupling dynamics of a flexible beam undergoing large overall motion was proposed based on virtual work principle, in which all the high-order terms related to coupling deformation were included in dynamic equations. Simulation examples of the flexible beam with prescribed rotation and free rotation were investigated. Numerical results show that the use of the first-order approximation coupling (FOAC) model may lead to a significant error when the flexible beam experiences large deformation or large deformation velocity. However, the correct solutions can always be obtained by using the present complete model. The difference in essence between this model and the FOAC model is revealed. These coupling high-order terms, which are ignored in FOAC model, have a remarkable effect on the dynamic behavior of the flexible body. Therefore, these terms should be included for the rigid-flexible dynamic modeling and analysis of flexible body undergoing motions with high speed.
基金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.
