In this paper,the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing.On the basis of the data taken from the spiral computerized tomography images of...In this paper,the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing.On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome(OSAHS),three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering.Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis.The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate.The results of the pressure and velocity distributions in the airflow field are quantitatively determined,and the displacement of soft palate is presented.Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing.However,the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway,and causes higher displacement of the soft palate.The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.展开更多
Debris flow can cause serious damage,and it is a part of the study of fluid-structure interaction(FSI).FSI analysis was established on the interaction between unsteady flow and the slit-trap dam's vibration,with a...Debris flow can cause serious damage,and it is a part of the study of fluid-structure interaction(FSI).FSI analysis was established on the interaction between unsteady flow and the slit-trap dam's vibration,with a coupling bench which can transfer fluid pressure and structure displacement.Debris flow can be seen as the Bingham body of incompressible.Based on ANSYS and CFX softwares,unidirectional and bidirectional coupling methods were used to study the transient interaction between debris flow and dam.The comparison between lateral fluid pressure states under different velocities and the equivalent stresses of the dam under different coupling conditions was made.The result shows that fluid-structure coupling becomes stronger with the increase of flow velocity.The maximum equivalent stress appears at the dam foundation,while the minimum equivalent stress appears at the dam abutment.With the increase of height,the fluid pressure decreases.The fluid pressure based on unidirectional FSI analysis is larger than that based on bidirectional FSI analysis and the maximum appears on the joint of the dam foundation and channel.The maximum equivalent stress of the dam based on the former is less than that based on the latter.展开更多
A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevo...A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.展开更多
Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for nume...Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for numerical simulation based on computational fluid dynamics and nonlinear finite element analysis methods. The typical results of the numerical simulation and wind tunnel experiment, including the overall lift and deformation, are in good agreement with each other. The results obtained indicate that the effect of fluid-structure interaction is noticeable and should be considered for non-rigid airships. Flow- induced deformation can further intensify the upward lift force and pitching moment, which can lead to a large deformation. Under a wind speed of 15 m/s, the lift force of the non-rigid model is increased to approximatelv 60% compared with that of the rigid model under a high angle of attack.展开更多
Based on the theories of the gas seepage in coal seams and the deformation of the coal-rock medium,the gas seepage field in coal-rock mass is coupled with the deformation field of the coal-rock mass to establish a flu...Based on the theories of the gas seepage in coal seams and the deformation of the coal-rock medium,the gas seepage field in coal-rock mass is coupled with the deformation field of the coal-rock mass to establish a fluidstructure interaction model for the interaction between coal gas and coal-rock masses.The outburst process in coal-rock masses under the joint action of gas pressure and crustal stress is simulated using the material point method.The simulation results show the changes in gas pressure,velocity distribution,maximum principal stress distribution,and damage distribution during the process of the coal and gas outburst,as well as themovement and accumulation of coal-rock masses after the occurrence of the outburst.It was found that the gas pressure gradient was greatest at theworking face after the occurrence of the outburst,the gas pressures and pressure gradients at each location within the coal seamgradually decreased with time,and the damage distribution was essentially the same as the minimum principal stress distribution.The simulation further revealed that the outburst first occurred in themiddle of the tunnel excavation face and that the speed at which particles of coal mass were ejected was highest at the center and decreased toward the upper and lower sides.The study provides a scientific basis for enhancing our understanding of the mechanism behind coal and gas outbursts,as well as their prevention and control.展开更多
In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aquedu...In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aqueduct deck, filled water, cables and support towers, respectively. A multi-element cable formulation is introduced to simulate the cable vibration. The dry （without water） and wet （with water） modes of the aqueduct bridge are both extracted and investigated in detail. The dry modes of the aqueduct bridge are basically similar to those of highway cable-stayed bridges. A dry mode may correspond to two types of wet modes, which are called the in-phase （with lower frequency） and out-of-phase （with higher frequency） modes. When the water-structure system vibrates in the in-phase/out-of-phase modes, the aqueduct deck moves and water sloshes in the same/opposite phase-angle, and the sloshing water may take different surface-wave modes. The wet modes of the system reflect the properties of interaction among the deck, towers, cables and water. The in-phase wet frequency generally decreases as the water depth increases, and the out-of-phase wet frequency may increase or decrease as the water depth increases.展开更多
Ice resistance prediction is a critical issue in the preliminary design of ships navigating brash ice conditions, which is closely related to the safety of a ship to navigate encounter brash ice, and has significant e...Ice resistance prediction is a critical issue in the preliminary design of ships navigating brash ice conditions, which is closely related to the safety of a ship to navigate encounter brash ice, and has significant effects on the kinds of propellers and motor power needed. In research on this topic, model tests and full-scale tests on ships have thus far been the primary approaches. In recent years, the application of the finite element method(FEM) has also attracted interest. Some researchers have conducted numerical simulations on ship–ice interactions using the fluid–structure interaction(FSI) method. This study used this method to predict and analyze the resistance of an ice-going ship, and compared the results with those of model ship tests conducted in a towing tank with synthetic ice to discuss the feasibility of the FEM. A numerical simulation and experimental methods were used to predict the brash ice resistance of an ice-going container ship model in a condition with three concentrations of brash ice(60%, 80%, and 90%). A comparison of the results yielded satisfactory agreement between the numerical simulation and the experiments in terms of both observed phenomena and resistance values, indicating that the proposed numerical simulation has significant potential for use in related studies in the future.展开更多
By using the shear stress transport （SST） model to predict the effect ot random now motion in a fluid zone, and using the Newmark method to solve the oscillation equations in a solid zone, a coupling model of the .p...By using the shear stress transport （SST） model to predict the effect ot random now motion in a fluid zone, and using the Newmark method to solve the oscillation equations in a solid zone, a coupling model of the .powerhouse and its tube water was developed. The effects of fluid-structure interaction are considered through the kinematic and dynamic conditions applied to the fluid-structure interfaces （FSI）. Numerical simulation of turbulent flow through the whole flow passage of the powerhouse and concrete structure vibration analysis in the time domain were carried out with the model. Considering the effect of coupling the turbulence and the powerhouse structure, the time history response of both turbulent flows through the whole flow passage and powerhouse structure vibration were generated. Concrete structure vibration analysis shows that the displacement, velocity, and acceleration of the dynamo floor respond dramatically to pressure fluctuations in the flow passage. Furthermore, the spectrum analysis suggests that pressure fluctuation originating from the static and dynamic disturbances of hydraulic turbine blades in the flow passage is one of the most important vibration sources.展开更多
The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for w...The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.展开更多
The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This m...The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.展开更多
Computational fluid dynamics(CFD)and the finite element method(FEM)are used to investigate the wind-driven dynamic response of cantilever traffic signal support structures as a whole.By building a finite element model...Computational fluid dynamics(CFD)and the finite element method(FEM)are used to investigate the wind-driven dynamic response of cantilever traffic signal support structures as a whole.By building a finite element model with the same scale as the actual structure and performing modal analysis,a preliminary understanding of the dynamic properties of the structure is obtained.Based on the two-way fluid-structure coupling calculation method,the wind vibration response of the structure under different incoming flow conditions is calculated,and the vibration characteristics of the structure are analyzed through the displacement time course data of the structure in the crosswind direction and along-wind direction.The results show that the maximum response of the structure increases gradually with the increase of wind speed under 90°wind direction angle,showing a vibration dispersion state,and the vibration response characteristics are following the vibration phenomenon of galloping;under 270°wind direction angle,the maximum displacement response of the structure occurs at the lower wind speed of 5 and 6m/s,and the vibration generated by the structure is vortex vibration at this time;the displacement response of the structure in along-wind direction increaseswith the increase of wind speed.The along-wind displacement response of the structure will increase with increasing wind speed,and the effective wind area and shape characteristics of the structurewill also affect the vibration response of the structure.展开更多
A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calcu...A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.展开更多
This paper studies the effect of ice resistance on the icebreaking capacity and speed of an icebreaking vessel.We combine an improved Correct Smoothed Particle Method(CSPM)with a material low-speed collision fracture ...This paper studies the effect of ice resistance on the icebreaking capacity and speed of an icebreaking vessel.We combine an improved Correct Smoothed Particle Method(CSPM)with a material low-speed collision fracture model to numerically simulate the continuous icebreaking and rolling process of crushed.Using this model,we investigate the icebreaking resistance and immersion resistance during the icebreaking process,taking into account the fluid(water)as the elastic boundary support and the fluid-solid coupling interaction.We compare the icebreaking resistance and broken ice fracture shapes obtained by the numerical calculation with the theoretical analytical results,and thus validate the improved CSPM method.Further,we compare the immersion resistance results from our simulation against that from Puntigliano[Puntigliano,Hamburgische Schiffbau-Versuchsanstalt GmbH(1995)],and demonstrate that the proposed method can accurately predict ice resistance.展开更多
A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid techn...A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.展开更多
The wind-induced vibration of the front windshield concerns the traffic safety and the aerodynamic characteristics of cars. In this paper, the numerical simulation and the experiment are combined to study the wind-ind...The wind-induced vibration of the front windshield concerns the traffic safety and the aerodynamic characteristics of cars. In this paper, the numerical simulation and the experiment are combined to study the wind-induced vibrations of the front windshield at different speeds of a van-body model bus. The Fluid-Structure Interaction （FSI） model is used for the finite element analysis of the vibration characteristics of the front windshield glass in the travelling process, and the wind-induced vibration response characteristics of the glass is obtained. A wind-tunnel experiment with an eddy current displacement sensor is carried out to study the deformation of the windshield at different wind speeds, and to verify the numerical simulation results. It is shown that the windshield of the model bus windshield undergoes a noticeable deformation as the speed changes, and from the deformation curve obtained, it is seen that in the accelerating process, the deformation of the glass increases as the speed increases, and with the speed being stablized, it also tends to a certain value. The results of this study can provide a scientific basis for the safety design of the windshield and the body.展开更多
Vortex-induced vibration(VIV)of a stay cablesubjected to a wind profile is numerically simulatedthrough combining computational fluid dynamics(CFD)code CFX 10.0 and computational structural dynamics(CSD)code ANSYS 10....Vortex-induced vibration(VIV)of a stay cablesubjected to a wind profile is numerically simulatedthrough combining computational fluid dynamics(CFD)code CFX 10.0 and computational structural dynamics(CSD)code ANSYS 10.0.A stay cable with the inclinedangle of 30°is used as the numerical model.Under aprofile of mean wind speed,unsteady aerodynamic liftcoefficients of the cable have been analyzed in both timedomain and frequency domain when VIV occurs.Theresults indicate that the lift coefficient wave response of thestay cable under a wind profile is different from that of aninfinitely long cable under a uniform flow in water(i.e.,without consideration of profile)obtained by directnumerical simulation.Cable oscillations can severelyaffect the unsteady aerodynamic frequencies,changeflow field distribution near the cable and affect the vortexshedding in the wake.展开更多
For the large deformation of the flexible body may cause the fluid grid distortion,which will make the numerical calculation tedious,even to end,the numerical simulation of the flexible body coupling with the fluid is...For the large deformation of the flexible body may cause the fluid grid distortion,which will make the numerical calculation tedious,even to end,the numerical simulation of the flexible body coupling with the fluid is always a tough problem.In this paper,the flexible body is under two kinds of constrained conditions and the ratio of length-diameter is 1:30.The Reynolds number of the airflow is 513,belonging to the area of low Reynolds number.The control equations of the coupling of flexible body with airflow are built and the adaptive grid control method is adopted to conduct the three-dimensional numerical simulation of the movement of the flexible body.The numerical results show that it is possible to simulate the characteristics of the flexible body's movement in the low Reynolds number airflow when the appropriate control equations are modeled and suitable equation-solving method is adopted.Unconstrained flexible body would turn over forward along the airflow's diffusion direction,while constrained flexible body in the flow field will make periodic rotation motion along the axis of the flexible body,and the bending deformation is more obvious than that of unconstrained flexible body.The preliminary three-dimensional numerical simulation can provide references for further research on the characteristics of the yam movement in high Reynolds number airflow.展开更多
The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifical...The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.展开更多
基金The project supported by the National Natural Science Foundation of China(10672036,10472025 and 10421002)the Natural Science Foundation of Liaoning Province(20032109)
文摘In this paper,the authors present airflow field characteristics of human upper airway and soft palate movement attitude during breathing.On the basis of the data taken from the spiral computerized tomography images of a healthy person and a patient with Obstructive Sleep Apnea-Hypopnea Syndrome(OSAHS),three-dimensional models of upper airway cavity and soft palate are reconstructed by the method of surface rendering.Numerical simulation is performed for airflow in the upper airway and displacement of soft palate by fluid-structure interaction analysis.The reconstructed threedimensional models precisely preserve the original configuration of upper airways and soft palate.The results of the pressure and velocity distributions in the airflow field are quantitatively determined,and the displacement of soft palate is presented.Pressure gradients of airway are lower for the healthy person and the airflow distribution is quite uniform in the case of free breathing.However,the OSAHS patient remarkably escalates both the pressure and velocity in the upper airway,and causes higher displacement of the soft palate.The present study is useful in revealing pathogenesis and quantitative mutual relationship between configuration and function of the upper airway as well as in diagnosingdiseases related to anatomical structure and function of the upper airway.
基金Science and Technology Support Program,China(No.2014BAL05B01)Project of Institute of Mountain Hazards and Environment of Chinese Academy of Sciences,China(No.KZZD-EW-Q5-Q1)
文摘Debris flow can cause serious damage,and it is a part of the study of fluid-structure interaction(FSI).FSI analysis was established on the interaction between unsteady flow and the slit-trap dam's vibration,with a coupling bench which can transfer fluid pressure and structure displacement.Debris flow can be seen as the Bingham body of incompressible.Based on ANSYS and CFX softwares,unidirectional and bidirectional coupling methods were used to study the transient interaction between debris flow and dam.The comparison between lateral fluid pressure states under different velocities and the equivalent stresses of the dam under different coupling conditions was made.The result shows that fluid-structure coupling becomes stronger with the increase of flow velocity.The maximum equivalent stress appears at the dam foundation,while the minimum equivalent stress appears at the dam abutment.With the increase of height,the fluid pressure decreases.The fluid pressure based on unidirectional FSI analysis is larger than that based on bidirectional FSI analysis and the maximum appears on the joint of the dam foundation and channel.The maximum equivalent stress of the dam based on the former is less than that based on the latter.
基金Supported by Marie Curie International Incoming Fellowship (No. PIIF-GA-2009-253453)
文摘A numerical simulation of the interaction between laminar flow with low Reynolds number and a highly flexible elastic sheet is presented. The mathematical model for the simulation includes a three-dimensional finitevolume based fluid solver for incompressible viscous flow and a combined finite-discrete element method for the three-dimensional deformation of solid. An immersed boundary method is used to couple the simulation of fluid and solid. It is implemented through a set of immersed boundary points scattered on the solid surface. These points provide a deformable solid wall boundary for the fluid by adding body force to Navier-Stokes equations. The force from the fluid is also obtained for each point and then applied on the boundary nodes of the solid. The vortex-induced vibration of the highly flexible elastic sheet is simulated with the established mathematical model. The simulated results for both swing pattern and oscillation frequency of the elastic sheet in low Reynolds number flow agree well with experimental data.
基金the National Natural Science Foundation of China (11202215 and 11332011)the Youth Innovation Promotion Association of CAS (2015015)
文摘Fluid-structure interaction is an important issue for non-rigid airships with inflated envelopes. In this study, a wind tunnel test is conducted, and a loosely coupled procedure is correspondingly established for numerical simulation based on computational fluid dynamics and nonlinear finite element analysis methods. The typical results of the numerical simulation and wind tunnel experiment, including the overall lift and deformation, are in good agreement with each other. The results obtained indicate that the effect of fluid-structure interaction is noticeable and should be considered for non-rigid airships. Flow- induced deformation can further intensify the upward lift force and pitching moment, which can lead to a large deformation. Under a wind speed of 15 m/s, the lift force of the non-rigid model is increased to approximatelv 60% compared with that of the rigid model under a high angle of attack.
基金The article received China National Natural Science Found(41601574).
文摘Based on the theories of the gas seepage in coal seams and the deformation of the coal-rock medium,the gas seepage field in coal-rock mass is coupled with the deformation field of the coal-rock mass to establish a fluidstructure interaction model for the interaction between coal gas and coal-rock masses.The outburst process in coal-rock masses under the joint action of gas pressure and crustal stress is simulated using the material point method.The simulation results show the changes in gas pressure,velocity distribution,maximum principal stress distribution,and damage distribution during the process of the coal and gas outburst,as well as themovement and accumulation of coal-rock masses after the occurrence of the outburst.It was found that the gas pressure gradient was greatest at theworking face after the occurrence of the outburst,the gas pressures and pressure gradients at each location within the coal seamgradually decreased with time,and the damage distribution was essentially the same as the minimum principal stress distribution.The simulation further revealed that the outburst first occurred in themiddle of the tunnel excavation face and that the speed at which particles of coal mass were ejected was highest at the center and decreased toward the upper and lower sides.The study provides a scientific basis for enhancing our understanding of the mechanism behind coal and gas outbursts,as well as their prevention and control.
基金National Natural Science Foundation of China Under Grant No.50678121Open Research Fund Program of State key Laboratory of Hydro-science and Engineering
文摘In this paper, a full-scale 3-D finite element model of the Jundushan cable-stayed aqueduct bridge is established with ANSYS Code. The shell, fluid, tension-only spar and beam elements are used for modeling the aqueduct deck, filled water, cables and support towers, respectively. A multi-element cable formulation is introduced to simulate the cable vibration. The dry （without water） and wet （with water） modes of the aqueduct bridge are both extracted and investigated in detail. The dry modes of the aqueduct bridge are basically similar to those of highway cable-stayed bridges. A dry mode may correspond to two types of wet modes, which are called the in-phase （with lower frequency） and out-of-phase （with higher frequency） modes. When the water-structure system vibrates in the in-phase/out-of-phase modes, the aqueduct deck moves and water sloshes in the same/opposite phase-angle, and the sloshing water may take different surface-wave modes. The wet modes of the system reflect the properties of interaction among the deck, towers, cables and water. The in-phase wet frequency generally decreases as the water depth increases, and the out-of-phase wet frequency may increase or decrease as the water depth increases.
基金financially supported by the National Natural Science Foundation of China(Grant No.51679052)the Natural Science Foundation of Heilongjiang Province of China(Grant No.E2018026)the Defense Industrial Technology Development Program(Grant No.JCKY2016604B001)
文摘Ice resistance prediction is a critical issue in the preliminary design of ships navigating brash ice conditions, which is closely related to the safety of a ship to navigate encounter brash ice, and has significant effects on the kinds of propellers and motor power needed. In research on this topic, model tests and full-scale tests on ships have thus far been the primary approaches. In recent years, the application of the finite element method(FEM) has also attracted interest. Some researchers have conducted numerical simulations on ship–ice interactions using the fluid–structure interaction(FSI) method. This study used this method to predict and analyze the resistance of an ice-going ship, and compared the results with those of model ship tests conducted in a towing tank with synthetic ice to discuss the feasibility of the FEM. A numerical simulation and experimental methods were used to predict the brash ice resistance of an ice-going container ship model in a condition with three concentrations of brash ice(60%, 80%, and 90%). A comparison of the results yielded satisfactory agreement between the numerical simulation and the experiments in terms of both observed phenomena and resistance values, indicating that the proposed numerical simulation has significant potential for use in related studies in the future.
基金supported by the National Natural Science Foundation of China (Grant No.90510017)
文摘By using the shear stress transport （SST） model to predict the effect ot random now motion in a fluid zone, and using the Newmark method to solve the oscillation equations in a solid zone, a coupling model of the .powerhouse and its tube water was developed. The effects of fluid-structure interaction are considered through the kinematic and dynamic conditions applied to the fluid-structure interfaces （FSI）. Numerical simulation of turbulent flow through the whole flow passage of the powerhouse and concrete structure vibration analysis in the time domain were carried out with the model. Considering the effect of coupling the turbulence and the powerhouse structure, the time history response of both turbulent flows through the whole flow passage and powerhouse structure vibration were generated. Concrete structure vibration analysis shows that the displacement, velocity, and acceleration of the dynamo floor respond dramatically to pressure fluctuations in the flow passage. Furthermore, the spectrum analysis suggests that pressure fluctuation originating from the static and dynamic disturbances of hydraulic turbine blades in the flow passage is one of the most important vibration sources.
基金supported by the CHN Energy United Power Technology Co.,Ltd.,China(Contract No.2020-75).
文摘The blades of large-scale wind turbines can obviously deform during operation,and such a deformation can affect the wind turbine’s output power to a certain extent.In order to shed some light on this phenomenon,for which limited information is available in the literature,a bidirectional fluid-structure interaction(FSI)numerical model is employed in this work.In particular,a 5 MW large-scale wind turbine designed by the National Renewable Energy Laboratory(NREL)of the United States is considered as a testbed.The research results show that blades’deformation can increase the wind turbine’s output power by 135 kW at rated working conditions.Compared with the outcomes of the simulations conducted using the model with no blade deformation,the results obtained with the FSI model are closer to the experimental data.It is concluded that the bidirectional FSI model can replicate the working conditions of wind turbines with great fidelity,thereby providing an effective method for wind turbine design and optimization.
文摘The development and rapid usage of numerical codes for fluid-structure interaction(FSI) problems are of great relevance to researchers in many engineering fields such as civil engineering and ocean engineering. This multidisciplinary field known as FSI has been expanded to engineering fields such as offshore structures, tall slender structures and other flexible structures applications. The motivation of this paper is to investigate the numerical model of two-way coupling FSI partitioned flexible plate structure under fluid flow. The adopted partitioned method and approach utilized the advantage of the existing numerical algorithms in solving the two-way coupling fluid and structural interactions. The flexible plate was subjected to a fluid flow which causes large deformation on the fluid domain from the oscillation of the flexible plate. Both fluid and flexible plate are subjected to the interaction of load transfer within two physics by using the strong and weak coupling methods of MFS and Load Transfer Physics Environment, respectively. The oscillation and deformation results have been validated which demonstrate the reliability of both strong and weak method in resolving the two-way coupling problem in contribution of knowledge to the feasibility field study of ocean engineering and civil engineering.
基金funded by the National Natural Science Foundation of China(Grant No.51578512)the Cultivating Fund Project for Young Teachers of Zhengzhou University(Grant No.JC21539028).
文摘Computational fluid dynamics(CFD)and the finite element method(FEM)are used to investigate the wind-driven dynamic response of cantilever traffic signal support structures as a whole.By building a finite element model with the same scale as the actual structure and performing modal analysis,a preliminary understanding of the dynamic properties of the structure is obtained.Based on the two-way fluid-structure coupling calculation method,the wind vibration response of the structure under different incoming flow conditions is calculated,and the vibration characteristics of the structure are analyzed through the displacement time course data of the structure in the crosswind direction and along-wind direction.The results show that the maximum response of the structure increases gradually with the increase of wind speed under 90°wind direction angle,showing a vibration dispersion state,and the vibration response characteristics are following the vibration phenomenon of galloping;under 270°wind direction angle,the maximum displacement response of the structure occurs at the lower wind speed of 5 and 6m/s,and the vibration generated by the structure is vortex vibration at this time;the displacement response of the structure in along-wind direction increaseswith the increase of wind speed.The along-wind displacement response of the structure will increase with increasing wind speed,and the effective wind area and shape characteristics of the structurewill also affect the vibration response of the structure.
基金supported by the National Natural Science Foundation of China (Grant Nos. 11822203 and 11532012)。
文摘A three-dimensional(3D) embedded Eulerian-Lagrangian method is proposed to simulate the 3D fluid-structure interaction(FSI) problems subjected to explosion and impact loading. This method achieves a high-quality calculation of fluid and structure deformation by adding Lagrangian particles to Eulerian grids. The overall computational domain is solved by the Eulerian method, and the Lagrangian particles with specified volume and influence domains are used to track structural deformations. The bidirectional mapping of physical quantities is achieved using the weighted average of the influence domain, which are based on the topological relationship between Eulerian grids and Lagrangian particles. Then, the data dependence solution and parallel algorithm realization are presented for the large-scale numerical calculations of explosion and impact problems. Additionally, the corresponding parallel program is developed based on the message passing interface(MPI) standard, and the parallel efficiency of parallel hydrocode are tested. The numerical results of typical explosion and impact problems are compared with corresponding experimental data to verify the effectiveness of the method. These comparisons show that the embedded EulerianLagrangian method successfully combine the advantages of both the Eulerian and Lagrangian methods to efficiently calculate the processes of large deformation and dynamic damage to the materials. The results presented in this work provide a useful reference point for further research on explosion and impact problems.
文摘This paper studies the effect of ice resistance on the icebreaking capacity and speed of an icebreaking vessel.We combine an improved Correct Smoothed Particle Method(CSPM)with a material low-speed collision fracture model to numerically simulate the continuous icebreaking and rolling process of crushed.Using this model,we investigate the icebreaking resistance and immersion resistance during the icebreaking process,taking into account the fluid(water)as the elastic boundary support and the fluid-solid coupling interaction.We compare the icebreaking resistance and broken ice fracture shapes obtained by the numerical calculation with the theoretical analytical results,and thus validate the improved CSPM method.Further,we compare the immersion resistance results from our simulation against that from Puntigliano[Puntigliano,Hamburgische Schiffbau-Versuchsanstalt GmbH(1995)],and demonstrate that the proposed method can accurately predict ice resistance.
基金by Key Aviation Scientific and Technological Laboratory of High-speed Hydrodynamic under grant MJ-2015-F-028.
文摘A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10972123,10802042)the Natural Science Foundation of Shandong Province(Grant No.Y2007A04)
文摘The wind-induced vibration of the front windshield concerns the traffic safety and the aerodynamic characteristics of cars. In this paper, the numerical simulation and the experiment are combined to study the wind-induced vibrations of the front windshield at different speeds of a van-body model bus. The Fluid-Structure Interaction （FSI） model is used for the finite element analysis of the vibration characteristics of the front windshield glass in the travelling process, and the wind-induced vibration response characteristics of the glass is obtained. A wind-tunnel experiment with an eddy current displacement sensor is carried out to study the deformation of the windshield at different wind speeds, and to verify the numerical simulation results. It is shown that the windshield of the model bus windshield undergoes a noticeable deformation as the speed changes, and from the deformation curve obtained, it is seen that in the accelerating process, the deformation of the glass increases as the speed increases, and with the speed being stablized, it also tends to a certain value. The results of this study can provide a scientific basis for the safety design of the windshield and the body.
基金This research was financially supported by the National Excellent Young Scientists Fund(No.50525823)the National Natural Science Foundation of China(Grant Nos.50738002 and 90715015).
文摘Vortex-induced vibration(VIV)of a stay cablesubjected to a wind profile is numerically simulatedthrough combining computational fluid dynamics(CFD)code CFX 10.0 and computational structural dynamics(CSD)code ANSYS 10.0.A stay cable with the inclinedangle of 30°is used as the numerical model.Under aprofile of mean wind speed,unsteady aerodynamic liftcoefficients of the cable have been analyzed in both timedomain and frequency domain when VIV occurs.Theresults indicate that the lift coefficient wave response of thestay cable under a wind profile is different from that of aninfinitely long cable under a uniform flow in water(i.e.,without consideration of profile)obtained by directnumerical simulation.Cable oscillations can severelyaffect the unsteady aerodynamic frequencies,changeflow field distribution near the cable and affect the vortexshedding in the wake.
基金supported by Zhejiang Provincial Natural Science Foundation under Grant No.LZ14E050004,LQ12A02002 etc
文摘For the large deformation of the flexible body may cause the fluid grid distortion,which will make the numerical calculation tedious,even to end,the numerical simulation of the flexible body coupling with the fluid is always a tough problem.In this paper,the flexible body is under two kinds of constrained conditions and the ratio of length-diameter is 1:30.The Reynolds number of the airflow is 513,belonging to the area of low Reynolds number.The control equations of the coupling of flexible body with airflow are built and the adaptive grid control method is adopted to conduct the three-dimensional numerical simulation of the movement of the flexible body.The numerical results show that it is possible to simulate the characteristics of the flexible body's movement in the low Reynolds number airflow when the appropriate control equations are modeled and suitable equation-solving method is adopted.Unconstrained flexible body would turn over forward along the airflow's diffusion direction,while constrained flexible body in the flow field will make periodic rotation motion along the axis of the flexible body,and the bending deformation is more obvious than that of unconstrained flexible body.The preliminary three-dimensional numerical simulation can provide references for further research on the characteristics of the yam movement in high Reynolds number airflow.
文摘The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.
