In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed ...In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.展开更多
The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels,making it challenging to perform hydraulic fracturing effectively.To solve this issue,it is essential to study inter-action b...The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels,making it challenging to perform hydraulic fracturing effectively.To solve this issue,it is essential to study inter-action behavior between hydraulic fractures(HFs)and gravels.A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method.This paper first investigates the dynamic evolution of HFs in glutenite,then analyzes the influences of various factors such as horizontal stress difference(D_(s)),minimum horizontal stress(s_(h)),gravel content(V_(g)),gravel size(d_(g)),and stiffness ratio of gravel to matrix(R_(s))on HF propagation geometries.Results show that pene-trating the gravel is the primary HF-gravel interaction behavior,which follows sequential and staggered initiation modes.Bypassing the gravel is the secondary behavior,which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(s_(H)).An offset along the gravel boundary is usually formed while penetrating gravels,and the offsets may cause fracture widths to decrease by 37.8%-84.4%.Even if stress dominates the direction of HF propagation,HFs still tend to deflect within gravels.The deviation angle from sH decreases with rising D_(s) and increases with the increase of d_(g) and R_(s).Additionally,intra-gravel shear HFs(IGS-HFs)are prone to be generated in coarse-grained glutenite under high D_(s),while more gravel-bypassing shear HFs(GBS-HFs)tend to be created in argillaceous glutenite with high R_(s) than in sandy glutenite with low R_(s).The findings above prompt the emergence of a novel HF propagation pattern in glutenite,which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.展开更多
A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as th...A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.展开更多
Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The b...Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.展开更多
Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication ...Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.展开更多
In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the s...In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.展开更多
In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lomb...In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lombardi(2012), and Trimarchi(2012). We focus on a 3D test case in which we evaluated the sensitivity of several input parameters on the fluid and structural results. In particular, this analysis provides a starting point from which we can look deeper into specific aspects of these simulations and analyze more realistic cases, e.g., in sails design. In this study, using the commercial software ADINATM, we addressed a well-known unsteadiness problem comprising a square box representing the fluid domain with a flexible bottom modeled with structural shell elements. We compared data from previously published work whose authors used the same numerical approach, i.e., a partitioned approach coupling a finite volume solver(for the fluid domain) and a finite element solver(for the solid domain). Specifically, we established several benchmarks and made comparisons with respect to fluid and solid meshes, structural element types, and structural damping, as well as solution algorithms. Moreover, we compared our method with a monolithic finite element solution method. Our comparisons of new and old results provide an outline of best practices for such simulations.展开更多
A pneumatic launcher is theoretically investigated to study its natural transverse vibration in water. Considering the mass effect of the sealing cover, the launcher is simplified as a uniform cantilever beam with a t...A pneumatic launcher is theoretically investigated to study its natural transverse vibration in water. Considering the mass effect of the sealing cover, the launcher is simplified as a uniform cantilever beam with a top point mass. By introducing the boundary and continuity conditions into the motion equation, the natural frequency equation and the mode shape function are derived. An iterative calculation method for added mass is also presented using the velocity potential function to account for the mass effect of the fluid on the launcher. The first 2 order natural frequencies and mode shapes are discussed in external flow fields and both external and internal flow fields. The results show good agreement with both natural frequencies and mode shapes between the theoretical analysis and the FEM studies. Also, the added mass is found to decrease with the increase of the mode shape orders of the launcher. And because of the larger added mass in both the external and internal flow fields than that in only the external flow field, the corresponding natural frequencies of the former are relatively smaller.展开更多
Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynam...Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.展开更多
In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation...In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.展开更多
In the current information society, the dissemination mechanisms and evolution laws of individual or collective opinions and their behaviors are the research hot topics in the field of opinion dynamics. First, in this...In the current information society, the dissemination mechanisms and evolution laws of individual or collective opinions and their behaviors are the research hot topics in the field of opinion dynamics. First, in this paper, a two-layer network consisting of an individual-opinion layer and a collective-opinion layer is constructed, and a dissemination model of opinions incorporating higher-order interactions(i.e. OIHOI dissemination model) is proposed. Furthermore, the dynamic equations of opinion dissemination for both individuals and groups are presented. Using Lyapunov's first method,two equilibrium points, including the negative consensus point and positive consensus point, and the dynamic equations obtained for opinion dissemination, are analyzed theoretically. In addition, for individual opinions and collective opinions,some conditions for reaching negative consensus and positive consensus as well as the theoretical expression for the dissemination threshold are put forward. Numerical simulations are carried to verify the feasibility and effectiveness of the proposed theoretical results, as well as the influence of the intra-structure, inter-connections, and higher-order interactions on the dissemination and evolution of individual opinions. The main results are as follows.(i) When the intra-structure of the collective-opinion layer meets certain characteristics, then a negative or positive consensus is easier to reach for individuals.(ii) Both negative consensus and positive consensus perform best in mixed type of inter-connections in the two-layer network.(iii) Higher-order interactions can quickly eliminate differences in individual opinions, thereby enabling individuals to reach consensus faster.展开更多
Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we ...Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we present a numerical method of studying nonlinear interactions between a long wave and multiple short harmonic waves in infinitely deep water. Specifically, this method is applied to the calculation of the temporal and spatial evolutions of the surface elevations in which a given long wave interacts with several short harmonic waves. Another important application of our method is to quantitatively analyse the nonlinear interactions between an arbitrary short wave train and another short wave train. From simulation results, we obtain that the mechanism for the nonlinear interactions between one short wave train and another short wave train （expressed as wave train 2） leads to the energy focusing of the other short wave train （expressed as wave train 3）. This mechanism occurs on wave components with a narrow frequency bandwidth, whose frequencies are near that of wave train 3.展开更多
This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance ...This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.展开更多
A method for analysing the vehicle-bridge interaction system with enhanced objectivity is proposed in the paper, which considers the time-variant and random characteristics and allows finding the power spectral densit...A method for analysing the vehicle-bridge interaction system with enhanced objectivity is proposed in the paper, which considers the time-variant and random characteristics and allows finding the power spectral densities(PSDs) of the system responses directly from the PSD of track irregularity. The pseudo-excitation method is adopted in the proposed framework, where the vehicle is modelled as a rigid body and the bridge is modelled using the finite element method. The vertical and lateral wheel-rail pseudo-excitations are established assuming the wheel and rail have the same displacement and using the simplified Kalker creep theory, respectively. The power spectrum function of vehicle and bridge responses is calculated by history integral. Based on the dynamic responses from the deterministic and random analyses of the interaction system, and the probability density functions for three safety factors(derailment coefficient, wheel unloading rate, and lateral wheel axle force) are obtained, and the probabilities of the safety factors exceeding the given limits are calculated. The proposed method is validated by Monte Carlo simulations using a case study of a high-speed train running over a bridge with five simply supported spans and four piers.展开更多
To investigate the effect of soil-pile-structure interaction(SPSI effect)on the dynamic response of a baseisolated structure with buried footings on a pile foundation,certain shake table tests are previously conducted...To investigate the effect of soil-pile-structure interaction(SPSI effect)on the dynamic response of a baseisolated structure with buried footings on a pile foundation,certain shake table tests are previously conducted.Based on the test results and the existing related studies,an efficient simplified model and a corresponding calculation method are verified for estimating the dynamic characteristics of a base-isolated structure with buried footings on a pile foundation with the SSI effect.In this method,the solutions by Veletsos and co-workers for a non-isolated structure with the SSI effect are verified and advanced for a base-isolated structure,and the solutions by Maravas and co-workers for a non-isolated structure on a pile foundation are introduced to consider the effect of the piles.By comparison with the shake table test,this work proves that the simplified method can efficiently estimate the dynamic responses of a base-isolated structure with buried footings on a pile foundation.Using parameter analysis,this work also shows that the dynamic characteristics of a non-isolated structure are quite similar to those of the base-isolated structure when the soil foundation is sufficiently soft,which means that the isolation layer gradually loses its isolation function as the soil foundation softens.展开更多
The Balise Transmission Module(BTM)unit of the on-board train control system is a crucial component.Due to its unique installation position and complex environment,this unit has a higher fault rate within the on-board...The Balise Transmission Module(BTM)unit of the on-board train control system is a crucial component.Due to its unique installation position and complex environment,this unit has a higher fault rate within the on-board train control system.To conduct fault prediction for the BTM unit based on actual fault data,this study proposes a prediction method combining reliability statistics and machine learning,and achieves the fusion of prediction results from different dimensions through multi-method interactive validation.Firstly,a method for predicting equipment fault time targeting batch equipment is introduced.This method utilizes reliability statistics to construct a model of the remaining faultless operating time distribution considering uncertainty,thereby predicting the remaining faultless operating probability of the BTM unit.Secondly,considering the complexity of the BTM unit’s fault mechanism,the small sample size of fault cases,and the potential presence of multiple fault features in fault text records,an individual-oriented fault prediction method based on Bayesian-optimized Gradient Boosting Regression Tree(Bayes-GBRT)is proposed.This method achieves better prediction results compared to linear regression algorithms and random forest regression algorithms,with an average absolute error of only 0.224 years for predicting the fault time of this type of equipment.Finally,a multi-method interactive validation approach is proposed,enabling the fusion and validation of multi-dimensional results.The results indicate that the predicted fault time and the actual fault time conform to a log-normal distribution,and the parameter estimation results are basically consistent,verifying the accuracy and effectiveness of the prediction results.The above research findings can provide technical support for the maintenance and modification of BTM units,effectively reducing maintenance costs and ensuring the safe operation of high-speed railway,thus having practical engineering value for preventive maintenance.展开更多
We investigate the geometries and energies of seven electronic states X-1A1, A1B1, a-3B1, B-1A2,b-3A2, C1B2 and c-3B2 of CF2 carbene using internally contracted multireference configuration interaction methods includi...We investigate the geometries and energies of seven electronic states X-1A1, A1B1, a-3B1, B-1A2,b-3A2, C1B2 and c-3B2 of CF2 carbene using internally contracted multireference configuration interaction methods including Davidson correction （icMRCIq-Q） with different basis sets aug-cc-pVXZ （X=T, Q, 5）. For the first time, the potential energy curves of electronic states of CF2 related icMRCI＋Q/aug-cc-pVTZ level. The ab initio results will and dynamics of electronic states of CF2 radical. to the lowest dissociation limit are calculated at the further increase our understanding of the structures展开更多
We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. ...We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.展开更多
We give a detailed examination of potential energy curves of the singlet and triplet states of CFC1 correlated with the lowest three dissociation limits. The calculations are carried out at the internally contracted m...We give a detailed examination of potential energy curves of the singlet and triplet states of CFC1 correlated with the lowest three dissociation limits. The calculations are carried out at the internally contracted multi- reference configuration interaction/cc-pV（T＋d）Z level with the other two geometric parameters fixed at the state equilibrium conformation. The vertical transition energy, the oscillator strength, the main configuration and the electron transition are also investigated at the same level.展开更多
Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to e...Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.展开更多
文摘In order to simulate and analyze the dynamic characteristics of the parachute from advanced tactical parachute system(ATPS),a nonlinear finite element algorithm and a preconditioning finite volume method are employed and developed to construct three dimensional parachute fluid-structure interaction(FSI)model.Parachute fabric material is represented by membrane-cable elements,and geometrical nonlinear algorithm is employed with wrinkling technique embedded to simulate the large deformations of parachute structure by applying the NewtonRaphson iteration method.On the other hand,the time-dependent flow surrounding parachute canopy is simulated using preconditioned lower-upper symmetric Gauss-Seidel(LU-SGS)method.The pseudo solid dynamic mesh algorithm is employed to update the flow-field mesh based on the complex and arbitrary motion of parachute canopy.Due to the large amount of computation during the FSI simulation,massage passing interface(MPI)parallel computation technique is used for all those three modules to improve the performance of the FSI code.The FSI method is tested to simulate one kind of ATPS parachutes to predict the parachute configuration and anticipate the parachute descent speeds.The comparison of results between the proposed method and those in literatures demonstrates the method to be a useful tool for parachute designers.
基金supported by the National Natural Science Foundation of China(Grant No.52304003)the Natural Science Foundation of Sichuan Province(Grant No.2024NSFSC0961)the Postdoctoral Fellowship Program of CPSF(Grant No.GZB20230090).
文摘The glutenite reservoir is strongly heterogeneous due to the random distribution of gravels,making it challenging to perform hydraulic fracturing effectively.To solve this issue,it is essential to study inter-action behavior between hydraulic fractures(HFs)and gravels.A coupled hydro-mechanical model is proposed for HF propagation in glutenite using a grain-based discrete element method.This paper first investigates the dynamic evolution of HFs in glutenite,then analyzes the influences of various factors such as horizontal stress difference(D_(s)),minimum horizontal stress(s_(h)),gravel content(V_(g)),gravel size(d_(g)),and stiffness ratio of gravel to matrix(R_(s))on HF propagation geometries.Results show that pene-trating the gravel is the primary HF-gravel interaction behavior,which follows sequential and staggered initiation modes.Bypassing the gravel is the secondary behavior,which obeys the sequential initiation mode and occurs when the orientation of the gravel boundary is inclined to the maximum horizontal stress(s_(H)).An offset along the gravel boundary is usually formed while penetrating gravels,and the offsets may cause fracture widths to decrease by 37.8%-84.4%.Even if stress dominates the direction of HF propagation,HFs still tend to deflect within gravels.The deviation angle from sH decreases with rising D_(s) and increases with the increase of d_(g) and R_(s).Additionally,intra-gravel shear HFs(IGS-HFs)are prone to be generated in coarse-grained glutenite under high D_(s),while more gravel-bypassing shear HFs(GBS-HFs)tend to be created in argillaceous glutenite with high R_(s) than in sandy glutenite with low R_(s).The findings above prompt the emergence of a novel HF propagation pattern in glutenite,which helps to understand the real HF geometries and to provide theoretical guidance for treatments in the field.
文摘A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.
文摘Fluid-structure interaction （FSI） problems in microchannels play a prominent role in many engineering applications. The present study is an effort toward the simulation of flow in microchannel considering FSI. The bottom boundary of the microchannel is simulated by size-dependent beam elements for the finite element method （FEM） based on a modified cou- ple stress theory. The lattice Boltzmann method （LBM） using the D2Q13 LB model is coupled to the FEM in order to solve the fluid part of the FSI problem. Because of the fact that the LBM generally needs only nearest neighbor information, the algorithm is an ideal candidate for parallel computing. The simulations are carried out on graphics processing units （GPUs） using computed unified device architecture （CUDA）. In the present study, the governing equations are non-dimensionalized and the set of dimensionless groups is exhibited to show their effects on micro-beam displacement. The numerical results show that the displacements of the micro-beam predicted by the size-dependent beam element are smaller than those by the classical beam element.
基金supported by the National Natural Science Foundations of China(Nos.50821063 and 50823004)973 Program(No.2007CB714701)the Fundamental Research Funds for the Central Universities(No.2010XS34)
文摘Based on the train-track coupling dynamics and high-speed train aerodynamics, this paper deals with an improved algorithm for fluid-structure interaction of high-speed trains. In the algorithm, the data communication between fluid solver and structure solver is avoided by inserting the program of train-track coupling dynamics into fluid dynamics program, and the relaxation factor concerning the load boundary of the fluid-structure interface is introduced to improve the fluctuation and convergence of aerodynamic forces. With this method, the fluid-structure dynamics of a highspeed train are simulated under the condition that the velocity of crosswind is 13.8 m/s and the train speed is 350 km/h. When the relaxation factor equals 0.5, the fluctuation of aerodynamic forces is lower and its convergence is faster than in other cases. The side force and lateral displacement of the head train are compared between off-line simulation and co-simulation. Simulation results show that the fluid-structure interaction has a significant influence on the aerodynam- ics and attitude of the head train under crosswind conditions. In addition, the security indexes of the head train worsen after the fluid-structure interaction calculation. Therefore, the fluid-structure interaction calculation is necessary for high-speed trains.
文摘In this paper, an investigation into the propagation of far field explosion waves in water and their effects on nearby structures are carried out. For the far field structure, the motion of the fluid surrounding the structure may be assumed small, allowing linearization of the governing fluid equations. A complete analysis of the problem must involve simultaneous solution of the dynamic response of the structure and the propagation of explosion wave in the surrounding fluid. In this study, a dynamic adaptive finite element procedure is proposed. Its application to the solution of a 2D fluid-structure interaction is investigated in the time domain. The research includes：a） calculation of the far-field scatter wave due to underwater explosion including solution of the time-depended acoustic wave equation, b） fluid-structure interaction analysis using coupled Euler-Lagrangian approach, and c） adaptive finite element procedures employing error estimates, and re-meshing. The temporal mesh adaptation is achieved by local regeneration of the grid using a time-dependent error indicator based on curvature of pressure function. As a result, the overall response is better predicted by a moving mesh than an equivalent uniform mesh. In addition, the cost of computation for large problems is reduced while the accuracy is improved.
文摘In this paper, we report our study on a numerical fluid-structure interaction problem originally presented by Mok et al.(2001) in two dimensions and later studied in three dimensions by Valdés Vazquez(2007), Lombardi(2012), and Trimarchi(2012). We focus on a 3D test case in which we evaluated the sensitivity of several input parameters on the fluid and structural results. In particular, this analysis provides a starting point from which we can look deeper into specific aspects of these simulations and analyze more realistic cases, e.g., in sails design. In this study, using the commercial software ADINATM, we addressed a well-known unsteadiness problem comprising a square box representing the fluid domain with a flexible bottom modeled with structural shell elements. We compared data from previously published work whose authors used the same numerical approach, i.e., a partitioned approach coupling a finite volume solver(for the fluid domain) and a finite element solver(for the solid domain). Specifically, we established several benchmarks and made comparisons with respect to fluid and solid meshes, structural element types, and structural damping, as well as solution algorithms. Moreover, we compared our method with a monolithic finite element solution method. Our comparisons of new and old results provide an outline of best practices for such simulations.
基金Foundation item： Supported by the National Natural Science Foundation of China （51379083） and the Specialized Research Fund for the Doctoral Program of Hiher Education （20120142110051）.
文摘A pneumatic launcher is theoretically investigated to study its natural transverse vibration in water. Considering the mass effect of the sealing cover, the launcher is simplified as a uniform cantilever beam with a top point mass. By introducing the boundary and continuity conditions into the motion equation, the natural frequency equation and the mode shape function are derived. An iterative calculation method for added mass is also presented using the velocity potential function to account for the mass effect of the fluid on the launcher. The first 2 order natural frequencies and mode shapes are discussed in external flow fields and both external and internal flow fields. The results show good agreement with both natural frequencies and mode shapes between the theoretical analysis and the FEM studies. Also, the added mass is found to decrease with the increase of the mode shape orders of the launcher. And because of the larger added mass in both the external and internal flow fields than that in only the external flow field, the corresponding natural frequencies of the former are relatively smaller.
基金supported by the National Natural Science Foundation of China(51879159,51490675,11432009,and 51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for doctoral student
文摘Nowadays,an increasing number of ships and marine structures are manufactured and inevitably operated in rough sea.As a result,some phenomena related to the violent fluid-elastic structure interactions(e.g.,hydrodynamic slamming on marine vessels,tsunami impact on onshore structures,and sloshing in liquid containers)have aroused huge challenges to ocean engineering fields.In this paper,the moving particle semi-implicit(MPS)method and finite element method(FEM)coupled method is proposed for use in numerical investigations of the interaction between a regular wave and a horizontal suspended structure.The fluid domain calculated by the MPS method is dispersed into fluid particles,and the structure domain solved by the FEM method is dispersed into beam elements.The generation of the 2D regular wave is firstly conducted,and convergence verification is performed to determine appropriate particle spacing for the simulation.Next,the regular wave interacting with a rigid structure is initially performed and verified through the comparison with the laboratory experiments.By verification,the MPS-FEM coupled method can be applied to fluid-structure interaction(FSI)problems with waves.On this basis,taking the flexibility of structure into consideration,the elastic dynamic response of the structure subjected to the wave slamming is investigated,including the evolutions of the free surface,the variation of the wave impact pressures,the velocity distribution,and the structural deformation response.By comparison with the rigid case,the effects of the structural flexibility on wave-elastic structure interaction can be obtained.
文摘In consideration of the problem that the effect of conduit structure on water hammer has been ignored in the classical theory,the Poisson coupling between the fluid and the pipeline was studied and a fourteen-equation mathematical model of fluid-structure interaction(FSI)was developed.Then,the transfer matrix method(TMM)was used to calculate the modal frequency,modal shape and frequency response.The results were compared with that in experiment to verify the correctness of the TMM and the results show that the fluid-structure coupling has a greater impact on the modal frequencies than the modal shape.Finally,the influence on the response spectrum of different damping ratios was studied and the results show that the natural frequency under different damping ratios has changed little but there is a big difference for the pressure spectrum.With the decreasing of damping ratio,the damping of the system on frequency spectrum is more and more significant and the dispersion and dissipation is more and more apparent.Therefore the appropriate damping ratio should be selected to minimize the effects of the vibration of the FSI.The results provide references for the theory research of FSI in the transient process.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.72031009 and 61473338)。
文摘In the current information society, the dissemination mechanisms and evolution laws of individual or collective opinions and their behaviors are the research hot topics in the field of opinion dynamics. First, in this paper, a two-layer network consisting of an individual-opinion layer and a collective-opinion layer is constructed, and a dissemination model of opinions incorporating higher-order interactions(i.e. OIHOI dissemination model) is proposed. Furthermore, the dynamic equations of opinion dissemination for both individuals and groups are presented. Using Lyapunov's first method,two equilibrium points, including the negative consensus point and positive consensus point, and the dynamic equations obtained for opinion dissemination, are analyzed theoretically. In addition, for individual opinions and collective opinions,some conditions for reaching negative consensus and positive consensus as well as the theoretical expression for the dissemination threshold are put forward. Numerical simulations are carried to verify the feasibility and effectiveness of the proposed theoretical results, as well as the influence of the intra-structure, inter-connections, and higher-order interactions on the dissemination and evolution of individual opinions. The main results are as follows.(i) When the intra-structure of the collective-opinion layer meets certain characteristics, then a negative or positive consensus is easier to reach for individuals.(ii) Both negative consensus and positive consensus perform best in mixed type of inter-connections in the two-layer network.(iii) Higher-order interactions can quickly eliminate differences in individual opinions, thereby enabling individuals to reach consensus faster.
基金Project supported by the National High Technology Research and Development Program of China (Grant No 2007AA12Z170)the Major Research Plan of the National Natural Science Foundation of China (Grant No 40706058)+1 种基金the Science-Technology Chenguang foundation for Young Scientist of Wuhan,China (Grant No 200850731388)the Canadian Space Agency Government Related Initiatives Program (GRIP) entitled Building Satellite Data into Fisheries and Oceans Operational Systems
文摘Although the nonlinear interactions between a single short gravity wave and a long wave can be solved analytically, the solution is less tractable in more general cases involving multiple short waves. In this work we present a numerical method of studying nonlinear interactions between a long wave and multiple short harmonic waves in infinitely deep water. Specifically, this method is applied to the calculation of the temporal and spatial evolutions of the surface elevations in which a given long wave interacts with several short harmonic waves. Another important application of our method is to quantitatively analyse the nonlinear interactions between an arbitrary short wave train and another short wave train. From simulation results, we obtain that the mechanism for the nonlinear interactions between one short wave train and another short wave train （expressed as wave train 2） leads to the energy focusing of the other short wave train （expressed as wave train 3）. This mechanism occurs on wave components with a narrow frequency bandwidth, whose frequencies are near that of wave train 3.
文摘This paper presents an improved unstructured grid immersed boundary method.The advantages of both immersed boundary method and body fitted grids which are generated by unstructured grid technology are used to enhance the computation efficiency of fluid structure interaction in complex domain.The Navier-Stokes equation was discretized spacially with collocated finite volume method and Euler implicit method in time domain.The rigid body motion was simulated by immersed boundary method in which the fluid and rigid body interface interaction was dealt with VOS(volume of solid) method.A new VOS calculation method based on graph was presented in which both immersed boundary points and cross points were collected in arbitrary order to form a graph.The method is verified with flow past oscillating cylinder.
文摘A method for analysing the vehicle-bridge interaction system with enhanced objectivity is proposed in the paper, which considers the time-variant and random characteristics and allows finding the power spectral densities(PSDs) of the system responses directly from the PSD of track irregularity. The pseudo-excitation method is adopted in the proposed framework, where the vehicle is modelled as a rigid body and the bridge is modelled using the finite element method. The vertical and lateral wheel-rail pseudo-excitations are established assuming the wheel and rail have the same displacement and using the simplified Kalker creep theory, respectively. The power spectrum function of vehicle and bridge responses is calculated by history integral. Based on the dynamic responses from the deterministic and random analyses of the interaction system, and the probability density functions for three safety factors(derailment coefficient, wheel unloading rate, and lateral wheel axle force) are obtained, and the probabilities of the safety factors exceeding the given limits are calculated. The proposed method is validated by Monte Carlo simulations using a case study of a high-speed train running over a bridge with five simply supported spans and four piers.
基金funded by the National Natural Science Foundation of China(No.51778282)the Science Research Foundations of Nanjing Institute of Technology(CKJA201505,JCYJ201618)
文摘To investigate the effect of soil-pile-structure interaction(SPSI effect)on the dynamic response of a baseisolated structure with buried footings on a pile foundation,certain shake table tests are previously conducted.Based on the test results and the existing related studies,an efficient simplified model and a corresponding calculation method are verified for estimating the dynamic characteristics of a base-isolated structure with buried footings on a pile foundation with the SSI effect.In this method,the solutions by Veletsos and co-workers for a non-isolated structure with the SSI effect are verified and advanced for a base-isolated structure,and the solutions by Maravas and co-workers for a non-isolated structure on a pile foundation are introduced to consider the effect of the piles.By comparison with the shake table test,this work proves that the simplified method can efficiently estimate the dynamic responses of a base-isolated structure with buried footings on a pile foundation.Using parameter analysis,this work also shows that the dynamic characteristics of a non-isolated structure are quite similar to those of the base-isolated structure when the soil foundation is sufficiently soft,which means that the isolation layer gradually loses its isolation function as the soil foundation softens.
基金supported by the Integrated Rail Transit Dispatch Control and Intermodal Transport Service Technology Project(Grant No.2022YFB4300500).
文摘The Balise Transmission Module(BTM)unit of the on-board train control system is a crucial component.Due to its unique installation position and complex environment,this unit has a higher fault rate within the on-board train control system.To conduct fault prediction for the BTM unit based on actual fault data,this study proposes a prediction method combining reliability statistics and machine learning,and achieves the fusion of prediction results from different dimensions through multi-method interactive validation.Firstly,a method for predicting equipment fault time targeting batch equipment is introduced.This method utilizes reliability statistics to construct a model of the remaining faultless operating time distribution considering uncertainty,thereby predicting the remaining faultless operating probability of the BTM unit.Secondly,considering the complexity of the BTM unit’s fault mechanism,the small sample size of fault cases,and the potential presence of multiple fault features in fault text records,an individual-oriented fault prediction method based on Bayesian-optimized Gradient Boosting Regression Tree(Bayes-GBRT)is proposed.This method achieves better prediction results compared to linear regression algorithms and random forest regression algorithms,with an average absolute error of only 0.224 years for predicting the fault time of this type of equipment.Finally,a multi-method interactive validation approach is proposed,enabling the fusion and validation of multi-dimensional results.The results indicate that the predicted fault time and the actual fault time conform to a log-normal distribution,and the parameter estimation results are basically consistent,verifying the accuracy and effectiveness of the prediction results.The above research findings can provide technical support for the maintenance and modification of BTM units,effectively reducing maintenance costs and ensuring the safe operation of high-speed railway,thus having practical engineering value for preventive maintenance.
基金Supported by the 2014 Postdoctoral Sustentation Fund of Qingdao under Grant No 01020120517the Natural Science Foundation of Shandong Province under Grant No ZR2014AP001+1 种基金the National Natural Science Foundation of China under Grant No11447226the Scientific Research Foundation of Shandong University of Science and Technology for Recruited Talents under Grant No 2015RCJJ015
文摘We investigate the geometries and energies of seven electronic states X-1A1, A1B1, a-3B1, B-1A2,b-3A2, C1B2 and c-3B2 of CF2 carbene using internally contracted multireference configuration interaction methods including Davidson correction （icMRCIq-Q） with different basis sets aug-cc-pVXZ （X=T, Q, 5）. For the first time, the potential energy curves of electronic states of CF2 related icMRCI＋Q/aug-cc-pVTZ level. The ab initio results will and dynamics of electronic states of CF2 radical. to the lowest dissociation limit are calculated at the further increase our understanding of the structures
文摘We review recent advances in the finite element method (FEM) simulations of interactions between waves and structures. Our focus is on the potential theory with the fully nonlinear or second-order boundary condition. The present paper has six sections. A review of previous work on interactions between waves and ocean structures is presented in Section one. Section two gives the mathematical formulation. In Section three, the finite element discretization, mesh generation and the finite element linear system solution methods are described. Section four presents numerical methods including time marching schemes, computation of velocity, remeshing and smoothing techniques and numerical radiation conditions. The application of the FEM to the wave-structure interactions are presented in Section five followed by the concluding remarks in Section six.
基金Supported by the National Natural Science Foundation of China under Grant No 11447148
文摘We give a detailed examination of potential energy curves of the singlet and triplet states of CFC1 correlated with the lowest three dissociation limits. The calculations are carried out at the internally contracted multi- reference configuration interaction/cc-pV（T＋d）Z level with the other two geometric parameters fixed at the state equilibrium conformation. The vertical transition energy, the oscillator strength, the main configuration and the electron transition are also investigated at the same level.
基金Supported by the National Natural Science Foundation of China(51379125,51490675,11432009,51579145)Chang Jiang Scholars Program(T2014099)+3 种基金Shanghai Excellent Academic Leaders Program(17XD1402300)Program for Professor of Special Appointment(Eastern Scholar)at Shanghai Institutions of Higher Learning(2013022)Innovative Special Project of Numerical Tank of the Ministry of Industry and Information Technology of China(2016-23/09)Lloyd’s Register Foundation for Doctoral Students
文摘Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.
