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.展开更多
In this paper, we obtain gradient estimates for certain nonlinear partial differential equations by coupling methods. First, we derive uniform gradient estimates for certain semi-linear PDEs based on the coupling meth...In this paper, we obtain gradient estimates for certain nonlinear partial differential equations by coupling methods. First, we derive uniform gradient estimates for certain semi-linear PDEs based on the coupling method introduced by Wang in 2011 and the theory of backward SDEs. Then we generalize Wang's coupling to the G-expectation space and obtain gradient estimates for nonlinear diffusion semigroups, which correspond to the solutions of certain fully nonlinear PDEs.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
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.展开更多
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.展开更多
To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Un...To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.展开更多
Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as eleme...Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.展开更多
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.展开更多
In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and t...In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.展开更多
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.展开更多
This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating som...This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.展开更多
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 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.展开更多
Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved ...Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.展开更多
Aiming at obtaining high power density of surface-mounted and interior permanent magnet synchronous motor(SIPMSM),it is important to accurately calculate the temperature field distribution of SIPMSM,and a magnetic-the...Aiming at obtaining high power density of surface-mounted and interior permanent magnet synchronous motor(SIPMSM),it is important to accurately calculate the temperature field distribution of SIPMSM,and a magnetic-thermal coupling method is proposed.The magnetic-thermal coupling mechanism is analyzed.The thermal network model and finite element model are built by this method,respectively.The effects of power frequency on iron losses and temperature fields are analyzed by the magnetic-thermal coupling finite element model under the condition of rated load,and the relationship between the load and temperature field is researched under the condition of the synchronous speed.In addition,the equivalent thermal network model is used to verify the magnetic-thermal coupling method.Then the temperatures of various nodes are obtained.The results show that there are advantages in both computational efficiency and accuracy for the proposed coupling method,which can be applied to other permanent magnet motors with complex structures.展开更多
We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement an...We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is devel- oped using C＋＋ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.展开更多
As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM...As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM） for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with that of both FEM and MPM展开更多
The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load typ...The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.展开更多
An adaptive EFG-FE coupling method is proposed and developed for the numerical simulation of lateral extrusion and forward-backward extrusion. Initially, the simulation has been implemented by using a conventional FE ...An adaptive EFG-FE coupling method is proposed and developed for the numerical simulation of lateral extrusion and forward-backward extrusion. Initially, the simulation has been implemented by using a conventional FE model. During the deforming process, mesh quality is checked at every incremental step. Distorted elements are automatically converted to EFG nodes, whereas, the less distorted elements are reserved. A new algorithm to generate EFG nodes and interface elements is presented. This method is capable of dealing with large deformation and has higher computational efficiency than using an EFG method wholly. Numerical results demonstrate that the adaptive EFG-FE coupling method has reasonable accuracy and is effective for local bulk metal forming such as extrusion processes.展开更多
For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics...For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.展开更多
基金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.
基金supported by NCMIS, National Natural Science Foundation of China (Grant Nos.11871458 and 11688101)Key Research Program of Frontier Sciences, Chinese Academy of Sciences (Grant No.QYZDB-SSW-SYS017)。
文摘In this paper, we obtain gradient estimates for certain nonlinear partial differential equations by coupling methods. First, we derive uniform gradient estimates for certain semi-linear PDEs based on the coupling method introduced by Wang in 2011 and the theory of backward SDEs. Then we generalize Wang's coupling to the G-expectation space and obtain gradient estimates for nonlinear diffusion semigroups, which correspond to the solutions of certain fully nonlinear PDEs.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
文摘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 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 the Fundamental Research Funds for the Central Universities(Grant No.2023YJS053)the National Natural Science Foundation of China(Grant No.52278386).
文摘To fundamentally alleviate the excavation chamber clogging during slurry tunnel boring machine(TBM)advancing in hard rock,large-diameter short screw conveyor was adopted to slurry TBM of Qingdao Jiaozhou Bay Second Undersea Tunnel.To evaluate the discharging performance of short screw conveyor in different cases,the full-scale transient slurry-rock two-phase model for a short screw conveyor actively discharging rocks was established using computational fluid dynamics-discrete element method(CFD-DEM)coupling approach.In the fluid domain of coupling model,the sliding mesh technology was utilized to describe the rotations of the atmospheric composite cutterhead and the short screw conveyor.In the particle domain of coupling model,the dynamic particle factories were established to produce rock particles with the rotation of the cutterhead.And the accuracy and reliability of the CFD-DEM simulation results were validated via the field test and model test.Furthermore,a comprehensive parameter analysis was conducted to examine the effects of TBM operating parameters,the geometric design of screw conveyor and the size of rocks on the discharging performance of short screw conveyor.Accordingly,a reasonable rotational speed of screw conveyor was suggested and applied to Jiaozhou Bay Second Undersea Tunnel project.The findings in this paper could provide valuable references for addressing the excavation chamber clogging during ultra-large-diameter slurry TBM tunneling in hard rock for similar future.
基金supported by the National Nature Science Foundation of China(Grant No.40874055)the Natural Science Foundation of the Hunan Province,China(Grant No.14JJ2012)
文摘Meshfree method offers high accuracy and computational capability and constructs the shape function without relying on predefined elements. We comparatively analyze the global weak form meshfree methods, such as element-free Galerkin method （EFGM）, the point interpolation method （PIM）, and the radial point interpolation method （RPIM）. Taking two dimensional Poisson equation as an example, we discuss the support-domain dimensionless size, the field nodes, and background element settings with respect to their effect on calculation accuracy of the meshfree method. RPIM and EFGM are applied to controlled- source two-dimensional electromagnetic modeling with fixed shape parameters. The accuracy of boundary conditions imposed directly and by a penalty function are discussed in the case of forward modeling of two-dimensional magnetotellurics in a homogeneous medium model. The coupling algorithm of EFG-PIM and EFG-RPIM are generated by integrating the PIM or RPIM and EFGM. The results of the numerical modeling suggest the following. First, the proposed meshfree method and corresponding coupled methods are well-suited for electromagnetic numerical modeling. The accuracy of the algorithm is the highest when the support-domain dimensionless size is 1.0 and the distribution of field nodes is consistent with the nodes of background elements. Second, the accuracy of PIM and RPIM are lower than that of EFGM for the Poisson equation but higher than EFGM for the homogeneous medium MT response. Third, RPIM overcomes the matrix inversion problem of PIM and has a wider selection of support-domain dimensionless sizes as compared to RPIM.
基金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 Science and Technology on Reactor System Design Technology Laboratory,Chengdu,China(LRSDT2020106)
文摘In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.
基金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 National Basic Research Program of China (Grant No 2006CB705500)Chang-Jiang Scholars and Innovative Research Team in University of China (Grant No IRT0605)the National Natural Science Foundation of China (Grant No70631001)
文摘This paper studies the cascading failure on random networks and scale-free networks by introducing the tolerance parameter of edge based on the coupled map lattices methods. The whole work focuses on investigating some indices including the number of failed edges, dynamic edge tolerance capacity and the perturbation of edge. In general, it assumes that the perturbation is attributed to the normal distribution in adopted simulations. By investigating the effectiveness of edge tolerance in scale-free and random networks, it finds that the larger tolerance parameter λ can more efficiently delay the cascading failure process for scale-free networks than random networks. These results indicate that the cascading failure process can be effectively controlled by increasing the tolerance parameter λ. Moreover, the simulations also show that, larger variance of perturbation can easily trigger the cascading failures than the smaller one. This study may be useful for evaluating efficiency of whole traffic systems, and for alleviating cascading failure in such systems.
基金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.
文摘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.
基金support of this work by the National Natural Science Foundation of China(No.51405096)the Fundamental Research Funds for the Central Universities(HEUCF210710).
文摘Based on Kirchhoff plate theory and the Rayleigh-Ritz method,the model for free vibration of rectangular plate with rectangular cutouts under arbitrary elastic boundary conditions is established by using the improved Fourier series in combination with the independent coordinate coupling method(ICCM).The effect of the cutout is taken into account by subtracting the energies of the cutouts from the total energies of the whole plate.The vibration displacement function of the hole domain is based on the coordinate system of the hole domain in this method.From the continuity condition of the vibration displacement function at the cutout,the transition matrix between the two coordinate systems is constructed,and the mass and stiffness matrices are completely obtained.As a result,the calculation is simplified and the computational efficiency of the solution is improved.In this paper,numerical examples and modal experiments are presented to validate the effectiveness of the modeling methods,and parameters related to influencing factors of the rectangular plate are analyzed to study the vibration characteristics.
基金This work was supported by Natural Science Foundation of China(Item number:51777060,U1361109)Natural Science Foundation of Henan province(Item number:162300410117)the he innovative research team plan of Henan Polytechnic University(Item number:T2015-2).
文摘Aiming at obtaining high power density of surface-mounted and interior permanent magnet synchronous motor(SIPMSM),it is important to accurately calculate the temperature field distribution of SIPMSM,and a magnetic-thermal coupling method is proposed.The magnetic-thermal coupling mechanism is analyzed.The thermal network model and finite element model are built by this method,respectively.The effects of power frequency on iron losses and temperature fields are analyzed by the magnetic-thermal coupling finite element model under the condition of rated load,and the relationship between the load and temperature field is researched under the condition of the synchronous speed.In addition,the equivalent thermal network model is used to verify the magnetic-thermal coupling method.Then the temperatures of various nodes are obtained.The results show that there are advantages in both computational efficiency and accuracy for the proposed coupling method,which can be applied to other permanent magnet motors with complex structures.
基金the support under A*STAR SERC grant (132-183-0025)
文摘We propose a multi-field implicit finite element method for analyzing the electromechanical behavior of dielectric elastomers. This method is based on a four-field variational principle, which includes displacement and electric potential for the electromechanical coupling analysis, and additional independent fields to address the incompressible constraint of the hyperelastic material. Linearization of the variational form and finite element discretization are adopted for the numerical implementation. A general FEM program framework is devel- oped using C＋＋ based on the open-source finite element library deal.II to implement this proposed algorithm. Numerical examples demonstrate the accuracy, convergence properties, mesh-independence properties, and scalability of this method. We also use the method for eigenvalue analysis of a dielectric elastomer actuator subject to electromechanical loadings. Our finite element implementation is available as an online supplementary material.
基金supported by the National Basic Research Program of China (2010CB832701)
文摘As a Lagrangian meshless method, the material point method （MPM） is suitable for dynamic problems with extreme deformation, but its efficiency and accuracy are not as good as that of the finite element method （FEM） for small deformation problems. Therefore, an algorithm for the coupling of FEM and MPM is proposed to take advantages of both methods. Furthermore, a conversion scheme of elements to particles is developed. Hence, the material domain is firstly discretized by finite elements, and then the distorted elements are automatically converted into MPM particles to avoid element entanglement. The interaction between finite elements and MPM particles is implemented based on the background grid in MPM framework. Numerical results are in good agreement with that of both FEM and MPM
基金This project is supported by National Natural Science Foundation of China (No.50077019).
文摘The coupling iteration (CI) of the finite element method(FEM) is used to simulate the magnetic and mechanical characteristics for a GMM actuator. The convergent ability under different prestress and different load types is investigated. Then the calculated deformations are compared with the experimental values. The results convince that the CI of FEM is suitable for the simulation of energy coupling and transformation mechanism of the GMM. At last, the output deformation properties are studied under different input currents, showing that there is a good compromise between good linearity and large strain under the prestress 6 MPa.
基金National Natural Science Foundation of China(No.50575143)Specialized Research Fund for the Doctoral Program of HigherEducation (No.20040248005).
文摘An adaptive EFG-FE coupling method is proposed and developed for the numerical simulation of lateral extrusion and forward-backward extrusion. Initially, the simulation has been implemented by using a conventional FE model. During the deforming process, mesh quality is checked at every incremental step. Distorted elements are automatically converted to EFG nodes, whereas, the less distorted elements are reserved. A new algorithm to generate EFG nodes and interface elements is presented. This method is capable of dealing with large deformation and has higher computational efficiency than using an EFG method wholly. Numerical results demonstrate that the adaptive EFG-FE coupling method has reasonable accuracy and is effective for local bulk metal forming such as extrusion processes.
文摘For accurate aeroelastic analysis,the unsteady rotor flowfield is solved by computational fluid dynamics(CFD)module based on RANS/Euler equations and moving-embedded grid system,while computational structural dynamics(CSD)module is introduced to handle blade flexibility.In CFD module,dual time-stepping algorithm is employed in temporal discretization,Jameson two-order central difference(JST)scheme is adopted in spatial discretization and B-L turbulent model is used to illustrate the viscous effect.The CSD module is developed based on Hamilton′s variational principles and moderate deflection beam theory.Grid deformation is implemented using algebraic method through coordinate transformations to achieve deflections with high quality and efficiency.A CFD/CSD loose coupling strategy is developed to transfer information between rotor flowfield and blade structure.The CFD and the CSD modules are verified seperately.Then the CFD/CSD loose coupling is adopted in airloads prediction of UH-60A rotor under high speed forward flight condition.The calculated results agree well with test data.Finally,effects of torsional stiffness properties on airloads of rotors with different tip swept angles(from 10° forward to 30° backward)are investigated.The results are evaluated through pressure distribution and airloads variation,and some meaningful conclusions are drawn the moderated shock wave strength and pressure gradient caused by varied tip swept angle and structural properties.
