CN111008418B - Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test - Google Patents
Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test Download PDFInfo
- Publication number
- CN111008418B CN111008418B CN201911124882.1A CN201911124882A CN111008418B CN 111008418 B CN111008418 B CN 111008418B CN 201911124882 A CN201911124882 A CN 201911124882A CN 111008418 B CN111008418 B CN 111008418B
- Authority
- CN
- China
- Prior art keywords
- building
- impact
- model
- equation
- numerical simulation
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test, which is characterized in that weak positions of a building are identified by simulating dynamic characteristics on a surrounding flow field and an incident flow surface of the building under different situations, so that the building design is optimized by adopting targeted measures; by means of the arrangement mode of the building groups, the stress characteristics of the building under the action of flood impact are calculated, so that the design of the building groups is optimized, and the safety performance of more buildings under the action of flood impact is guaranteed.
Description
Technical Field
The invention belongs to the field of building design tests, and relates to a building group energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model tests.
Background
Flood disaster effects are manifested by loss of life and property of personnel caused by flooding or flushing buildings. The impact of floods is a major driving factor in house collapse in floods. During the impact, when the wave crest just contacts the structure flow-facing surface, impact pressure with short duration and extremely high strength exists, and the extremely high impact load can cause instability of the building or cause local damage. Experience shows that building safety is not damaged, 30% of efficiency of post-disaster reconstruction work is improved, and 60% of personnel and property loss is reduced. Therefore, aiming at the characteristics of the flood impact effect, the analysis of the dynamic response and the damage process of the building under the flood impact effect is carried out, the impact damage mechanism and the characteristics of the residential building under the flood impact effect are revealed, the vitality of the residential building after the residential building is subjected to the flood impact effect is evaluated, and the method is very important for improving the flood control safety of the residential building in the flood multiple areas.
Currently, response studies of buildings under the effect of flood impact have become mature. The impact load damages the characteristics of the building, such as forward direction, mainly impact, wave change along the flow direction, shielding effect and the like. The Sho cloud, ge Xueli and the like develop model tests to research the effect of floods of different quantity levels on houses of different opening rates, the impact pressure characteristics of flood impact on building structures are depicted, the relation among impact load, water flow load and hydrostatic pressure is analyzed, the change rules of water flow load, bending moment and resultant force borne by the model by different water heads and different opening rates are discussed, and meanwhile, the structural damage process is reproduced by numerical simulation, so that technical guidance is provided for the flood control and disaster reduction design of village buildings. The grand et al set up the model test and study the damage form of the brick-concrete residential building of small scale under the effect of flood impact, have guiding effect on improving the design of building safety and stability. Hu et al take the instantaneous effect of impact load into consideration, establish a numerical model and perform physical test to study the interaction of dam-break flood and structures in a flow field, which shows that when high-speed water flow violently impacts the surface of a building, the building is easily damaged locally due to strong impact force. However, the study of the curtailment effect of different arrangements of building groups on floods remains blank. When the mountain floods are out, the arrangement modes of the buildings are different, so that the front-row buildings can resist flood impact to consume energy, and a protective barrier is provided for the following buildings, thereby maximizing the safety of the protection building group.
The invention adopts a mode of combining numerical simulation and physical model test to study the response of a single building and a building group under the action of flood impact. Through the impact test of flood on the single building, the impact force and distribution rule of the flood on the building and the water flow pressure on the surface of the house after the impact are researched. Meanwhile, typical observation points are selected, and the change of physical quantities of flow fields such as water depth, flow velocity, pressure intensity and the like in the flood impact process is measured. Impact tests of building groups aim at researching the resistance effect of different arrangement modes of buildings on flood, and maximizing the safety protection of the buildings.
Disclosure of Invention
The invention aims to overcome the defects of the prior art, and provides a building energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test, which can identify weak positions of a building so as to optimize the design of the building by adopting targeted measures; ensuring the safety performance of more buildings under the action of flood impact.
The invention solves the technical problems by the following technical proposal:
the building energy dissipation and impact prevention optimization design method based on the combination of numerical simulation and physical model test is characterized by comprising the following steps of: the design method comprises the following steps:
s1, a physical model experiment operation step;
s2, performing numerical simulation;
the physical model experiment operation steps in the step S1 comprise:
s11, manufacturing a water tank and organic glass house test model: the water tank is 14m multiplied by 1.5m, the gradient is 1%, wherein the water tank head section of 5m multiplied by 1.5m is used as an upstream reservoir, and is separated from the downstream by a gate, the rest of the water tank simulates residential areas, and house models are respectively installed at positions 3m, 5m and 7m away from the gate. According to a model similarity theory, taking a typical village and town building as a house prototype (4.5 m multiplied by 4.2m multiplied by 3 m), and manufacturing an organic glass house model according to a model geometric scale of 1/6;
s12, instrument layout: the impact pressure measurement adopts a DS-30 type data acquisition system, the measuring range of a pressure sensor is 30kPa and 10kPa, the self-oscillation frequency is 500Hz, the sampling interval is 0.023s, 30 pressure sensors are adopted for measurement, wherein 10 pressure sensors and 20 pressure sensors are arranged on the upstream face of a house model to capture the impact pressure time course change; at the beginning of the experiment, the gate was opened by gravity-controlled rapid flip-open in the flow direction, and the impact procedure was photographed from the side using a high-speed digital camera capable of acquiring 100 frames per second at a resolution of 1280 x 720. Measuring flow field information near the building by using a flow velocity meter, and adopting repeated tests for a plurality of times to ensure the measurement accuracy of impact load;
s13, carrying out a model test: the gate is closed, the water tank is stored to the water depths of 0.6m, 0.9m and 1.2m respectively at the first section, after the water level is stable, the gate control switch is immediately started, the bottom of the gate is quickly turned over under the lifting of gravity, water flow quickly impacts the house model, and meanwhile, each measuring device captures and records flow field data. The tail water flows into the reservoir again through the channel;
the numerical simulation step in the step S2 comprises the following steps:
s21, adopting three groups of test schemes for numerical simulation
Scheme one: according to the scheme that the transverse shafts of the buildings are perpendicular to the incoming flow direction, 3 house models are arranged at positions 3m, 5m and 7m away from the gate, impact mechanical characteristics of the water flow directly-rushing building acting on the building and protection effects of the water flow directly-rushing building on subsequent buildings are discussed, 3 groups of comparison tests are respectively arranged, the house models are independently arranged at positions 3m, 5m and 7m away from the gate, and the deformation characteristics of the building under the action of flood impact of the same magnitude are independently born under the protection effect of the front-row-free building are verified;
scheme II: the arrangement mode that the transverse axis of the building is parallel to the incoming flow direction is also provided with 3 groups of comparison tests, and the shielding effect of the front row of buildings on the rear row of buildings under the impact of different water heads is researched;
scheme III: the arrangement mode of the building group has an influence on the mechanical property of water flow impact, and the blocking effect of the arrangement mode of the building on the water flow is researched by adopting two typical arrangement modes from physical and mechanical characteristics, so that references are provided for optimizing the arrangement mode of the downstream building of the dam, improving the response of the downstream building to the flood impact and maximally guaranteeing the safety of the building;
s22, constructing a flood impact numerical model
The Reynolds average Navier-Stokes equation and the RNG k- ε turbulence model are used to solve the over-dam water flow motion and turbulence, and the continuous equation is as follows:
wherein ui represents average speed, xi is dimension, t is time, p is pressure, ρ is fluid density, gi is gravitational acceleration component, v is molecular kinematic viscosity, vt is vortex kinematic viscosity;
the reynolds average Navier-Stokes equation is a control equation for the flow field average variable, assuming that the flow field variable consists of a time average and a pulsation, by introducing the Boussinesq assumption that the turbulent reynolds stress is proportional to strain, the turbulent calculation is attributed to the calculation of the proportionality coefficient between reynolds stress and strain (i.e. turbulent viscosity coefficient). The turbulent vortex viscosity is adopted in the chapter and calculated by turbulent energy k and turbulent dissipation rate epsilon in the following form:
wherein, C is a constant, and the equation can not be closed due to the Reynolds stress introduced after homogenizing in the flow field, namely, the number of unknown numbers is more than the number of equation, so that various turbulence models are generated, RNG k-epsilon and model adopt corrected turbulence consumption and dissipation to represent the turbulence intensity, and the transportation equations of k, epsilon and model are as follows:
wherein R is ε For turbulent shear performance, G is the rate of turbulent kinetic energy, and the coefficients in the model are as follows:
C μ ＝0.085，C ε1 ＝1.42，C ε2 ＝1.68，σ k ＝σ ε = 0.7194, β=0.012, η=4.38, these coefficients being obtained by data fusion with a large number of statistics, widely applicable to turbulent flow characteristics;
VOF (Volume of Fluid) method is used to solve for free surface motion, tracking of gas and liquid free surfaces is solved by continuous equations of the form:
any finite control body being filled with gas or liquid or a mixture of both, depending on the volume fraction alpha of the liquid w The volume fraction within a given control volume is a constant value;
the time step is determined from the Courant-Friedrichs-Lewy (CFL) condition, which in the three-dimensional model has the form:
wherein, C is a dimensionless constant, the value of which depends on a specific equation to be solved, a display time pushing format is adopted, and the value of C is set to be 1;
s23, solving structural deformation
S231, solid phase control equation
Conservation of the solid fraction may be according to newton's second law:
balance equation:
geometric equation:
the physical equation: sigma (sigma) ij ＝λσ ij ε kk +2με il
Wherein ki is dynamic pressure load component, m/s2; epsilon ij Is Cauchy strain constant, s-1; lambda, mu are Lame constants of the carbon steel bent pipes respectively; e is the elastic modulus of the protective film, kg m -1 s - 2, v is Poisson's ratio；
S232 fluid-solid coupling control equation
The fluid-solid coupling interface should satisfy the displacement, stress, temperature, etc. of the fluid and the solid, and also follow the corresponding basic conservation law, and in the fluid-solid coupling interface, the conservation of the variables such as the stress, displacement, temperature, heat flow, etc. of the fluid and the solid should be satisfied, namely the following 4 control equations are satisfied:
τ f ·n f ＝τ s ·n s
d f ＝d s
T f ＝T s
q f ＝q s
wherein the subscript f represents fluid, the subscript s represents solid, τ represents solid stress, d represents displacement, T represents temperature, q represents heat flow
In addition, in the step S11, the organic glass house model is a single-open building, so that the dynamic conditions of the flood impact time course and the impact process can be accurately captured.
And in the step S11, the bottom slope of the water tank is 1% to simulate the impact effect of the dam-breaking flood.
In addition, in the step S11, the opening speed of the gate is controlled by a counterweight, so that the gate can be opened instantaneously, the process of impacting the house by a large amount of flood after dam-break flood is simulated fully, and the consistency of the model test and the actual situation is ensured.
Further, the control test in the step S21 of the step S2 numerical simulation is the same as the physical model test set for verifying the accuracy of the numerical simulation.
The invention has the advantages and beneficial effects that:
1. according to the invention, weak positions of the building are identified by simulating dynamic characteristics on the surrounding flow field and the flow face of the building under different situations, so that the design of the building is optimized by adopting targeted measures; by means of the arrangement mode of the building groups, the stress characteristics of the building under the action of flood impact are calculated, so that the design of the building groups is optimized, and the safety performance of more buildings under the action of flood impact is guaranteed.
Drawings
FIG. 1 is a numerical sink simulation scheme setup diagram;
FIG. 2 is a schematic diagram of a group layout simulation scheme;
FIG. 3 is a diagram of a building structure model;
FIG. 4 is a flow-face pressure sensor layout;
FIG. 5 is a graph of the pressure change over time at different points of the head-on surface of a 3m house from a gate (a: 1.2m head; b:0.9m head; c:0.6m head);
fig. 6 is a structural deformation diagram under the impact of water flow.
Detailed Description
The invention is further illustrated by the following examples, which are intended to be illustrative only and not limiting in any way.
The building energy dissipation and impact prevention optimization design method based on the combination of numerical simulation and physical model test is characterized by comprising the following steps of: the design method comprises the following steps:
s1, a physical model experiment operation step;
the physical model experiment operation steps in the step S1 comprise:
s11, manufacturing a water tank and organic glass house test model: the water tank is 14m multiplied by 1.5m, the gradient is 1%, wherein the water tank head section of 5m multiplied by 1.5m is used as an upstream reservoir, and is separated from the downstream by a gate, the rest of the water tank simulates residential areas, and house models are respectively installed at positions 3m, 5m and 7m away from the gate. According to a model similarity theory, taking a typical village and town building as a house prototype (4.5 m multiplied by 4.2m multiplied by 3 m), and manufacturing an organic glass house model according to a model geometric scale of 1/6; the size of a single-bay house is 0.74m multiplied by 0.69m multiplied by 0.52m, the house model is shown in fig. 3, the physical model test scheme setting is shown in the comparison test of fig. 1, and the invention only aims at developing a model test of a single building in a flow field. The relationship between the physical model and the prototype main scale is shown in table 1.
TABLE 1 model and prototype Primary Scale relationship
S12, instrument layout: the impact pressure measurement adopts a DS-30 type data acquisition system, the measuring range of a pressure sensor is 30kPa and 10kPa, the self-oscillation frequency is 500Hz, the sampling interval is 0.023s, 30 pressure sensors are adopted for measurement, wherein 10 pressure sensors and 20 pressure sensors are arranged on the upstream face of a house model to capture the impact pressure time course change; at the beginning of the experiment, the gate was opened by gravity-controlled rapid flip-open in the flow direction, and the impact procedure was photographed from the side using a high-speed digital camera capable of acquiring 100 frames per second at a resolution of 1280 x 720. Measuring flow field information near the building by using a flow velocity meter, and adopting repeated tests for a plurality of times to ensure the measurement accuracy of impact load;
s13, carrying out a model test: the gate is closed, the water tank is stored to the water depths of 0.6m, 0.9m and 1.2m respectively at the first section, after the water level is stable, the gate control switch is immediately started, the bottom of the gate is quickly turned over under the lifting of gravity, water flow quickly impacts the house model, and meanwhile, each measuring device captures and records flow field data. The tail water flows into the water storage pool again through the channel;
the numerical simulation step in the step S2 comprises the following steps:
s21, adopting three groups of test schemes for numerical simulation
Scheme one: according to the scheme that the transverse shafts of the buildings are perpendicular to the incoming flow direction, 3 house models are arranged at positions 3m, 5m and 7m away from the gate, impact mechanical characteristics of the water flow directly-rushing building acting on the building and protection effects of the water flow directly-rushing building on subsequent buildings are discussed, 3 groups of comparison tests are respectively arranged, the house models are independently arranged at positions 3m, 5m and 7m away from the gate, and the deformation characteristics of the building under the action of flood impact of the same magnitude are independently born under the protection effect of the front-row-free building are verified;
scheme II: the arrangement mode that the transverse axis of the building is parallel to the incoming flow direction is also provided with 3 groups of comparison tests, and the shielding effect of the front row of buildings on the rear row of buildings under the impact of different water heads is researched;
scheme III: the arrangement mode of the building group has an influence on the mechanical property of water flow impact, and the blocking effect of the arrangement mode of the building on the water flow is researched by adopting two typical arrangement modes from physical and mechanical characteristics, so that references are provided for optimizing the arrangement mode of the downstream building of the dam, improving the response of the downstream building to the flood impact and maximally guaranteeing the safety of the building.
In addition, in the step S11, the organic glass house model is a single-open building, so that the dynamic conditions of the flood impact time course and the impact process can be accurately captured.
And in the step S11, the bottom slope of the water tank is 1% to simulate the impact effect of the dam-breaking flood.
In addition, in the step S11, the opening speed of the gate is controlled by a counterweight, so that the gate can be opened instantaneously, the process of impacting the house by a large amount of flood after dam-break flood is simulated fully, and the consistency of the model test and the actual situation is ensured.
Further, the control test in the step S21 of the step S2 numerical simulation is the same as the physical model test set for verifying the accuracy of the numerical simulation.
S22, constructing a flood impact numerical model
The Reynolds average Navier-Stokes equation and the RNG k- ε turbulence model are used to solve the over-dam water flow motion and turbulence, and the continuous equation is as follows:
where ui represents the average velocity, xi is the dimension, t is the time, p is the pressure, ρ is the fluid density, gi is the gravitational acceleration component, v is the molecular kinematic viscosity, and vt is the vortex motion viscosity.
The reynolds average Navier-Stokes equation is a control equation for the flow field average variable, assuming that the flow field variable consists of a time average and a pulsation, by introducing the Boussinesq assumption that the turbulent reynolds stress is proportional to strain, the turbulent calculation is attributed to the calculation of the proportionality coefficient between reynolds stress and strain (i.e. turbulent viscosity coefficient). The turbulent vortex viscosity is adopted in the chapter and calculated by turbulent energy k and turbulent dissipation rate epsilon in the following form:
wherein, C is a constant, and the equation can not be closed due to the Reynolds stress introduced after homogenizing in a flow field, namely, the number of unknown numbers is more than the number of equation, so that various turbulence models are generated, RNG k-epsilon and model adopt corrected turbulence dissipation rates to represent turbulence intensity, and the transportation equations of k, epsilon and model are as follows:
wherein R is ε For turbulent shear performance, G is the rate of turbulent kinetic energy, and the coefficients in the model are as follows:
C u ＝0.085，C ε1 ＝1.42，C ε2 ＝1.68，σ k ＝σ ε = 0.7194, β=0.012, η=4.38. These coefficients are obtained by data fusion using a large amount of statistical data and are widely applicable to turbulent flow characteristics.
VOF (Volume ofFluid) method is used to solve for free surface motion. The tracking of the gas and liquid free surfaces is solved by the following form of continuous equation:
any finite control body being filled with gas or liquid or a mixture of both, depending on the volume fraction alpha of the liquid w The volume fraction within a given control volume is a constant value. Once the liquid volume fraction is determined, the volume fraction of air can be determined, with the free surface at a location where both volume fractions are 0.5.
Based on solving the control equation, dam flood discharge and flood evolution numerical simulation can be smoothly carried out, and the control equation is subjected to numerical dispersion by adopting a finite volume method. In each control volume, calculation is performed according to the volume average of the variables, the pressure, the volume fraction, the density, the viscosity, the turbulence energy and the turbulence dissipation rate are solved one by one, the surface flux, the surface stress and the volume force on the control body are solved by the conservation equation, and the surface flux, the surface stress and the volume force on the control body are also solved by one. The time step is determined from the Courant-Friedrichs-Lewy (CFL) condition, which in the three-dimensional model has the form:
wherein, C is a dimensionless constant, the value of which depends on a specific equation to be solved, a display time pushing format is adopted, and the value of C is set to be 1.
S23, solving structural deformation
S231, solid phase control equation
Conservation of the solid fraction may be according to newton's second law:
balance equation:
geometric equation:
the physical equation: sigma (sigma) ij ＝λσ ij ε kk +2με ij
Wherein ki is dynamic pressure load component, m/s2; epsilon ij Is Cauchy strain constant, s-1; lambda, mu are Lame constants of the carbon steel bent pipes respectively; e is the elastic modulus of the protective film, kg m -1 s - 2, v is poisson's ratio;
s232 fluid-solid coupling control equation
The fluid-solid coupling interface should satisfy the displacement, stress, temperature, etc. of the fluid and the solid, and also follow the corresponding basic conservation law, and in the fluid-solid coupling interface, the conservation of the variables such as the stress, displacement, temperature, heat flow, etc. of the fluid and the solid should be satisfied, namely the following 4 control equations are satisfied:
τ f ·n f ＝τ s ·n s
d f ＝d s
T f ＝T s
q f ＝q s
wherein the subscript f represents fluid, the subscript s represents solid, τ represents solid stress, d represents displacement, T represents temperature, and q represents heat flow.
Fig. 5 shows the head-on pressure profile of a house at a position 3m from the gate under different head impacts. The deformation of the building is determined by integrating the building pressure to determine the force on the building surface.
Figure 6 shows the deformation characteristics of each building group under the action of a 1.2m head impact. In general, the closer to the crumple, the greater the impact force is received by the building, while the highest impact pressure is received by the bottom of the head-on of the building, and the lower the impact pressure is received by the top, but the more easily deformed is unstable. The deformation characteristics of the building group fully prove the different capability of different arrangement modes to cope with flood impact. The front row of buildings can cut down a large amount of energy, thereby protecting the rear row of buildings. In a flood-prone area, the design mode of reducing the flood impact effect to protect the safety of building groups by building guide walls is worth recommending. The diversion wall bears a large amount of energy and changes the direction of flood, so that the influence of the flood on the residential area is reduced.
Although the embodiments of the present invention and the accompanying drawings have been disclosed for illustrative purposes, those skilled in the art will appreciate that: various substitutions, changes and modifications are possible without departing from the spirit and scope of the invention and the appended claims, and therefore the scope of the invention is not limited to the embodiments and the disclosure of the drawings.
Claims (5)
1. The building energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test is characterized by comprising the following steps of: the design method comprises the following steps:
s1, a physical model experiment operation step;
s2, performing numerical simulation;
the physical model experiment operation steps in the step S1 comprise:
s11, manufacturing a water tank and organic glass house test model: the water tank is 14m multiplied by 1.5m, the gradient is 1%, wherein the water tank head section of 5m multiplied by 1.5m is used as an upstream reservoir, the water tank head section and the downstream water tank are separated by a gate, residential areas are simulated by the rest part of the water tank, house models are respectively installed at positions 3m, 5m and 7m away from the gate, a typical village building is taken as a house prototype of 4.5m multiplied by 4.2m multiplied by 3m according to a model similarity theory, and a plexiglass house model is manufactured according to a model geometric scale of 1/6;
s12, instrument layout: the impact pressure measurement adopts a DS-30 type data acquisition system, the measuring range of the pressure sensor is 30kPa and 10kPa, the self-oscillation frequency is 500Hz, the sampling interval is 0.023s, 30 pressure sensors are adopted for measurement, wherein 10 pressure sensors and 20 pressure sensors are arranged on the windward side of the house model to capture the time course change of the impact pressure; when an experiment starts, the gate is controlled by gravity to quickly turn over and open along the flow direction, a high-speed digital camera is used for shooting the impact process from the side, the camera acquires 100 frames of images per second with the resolution of 1280 multiplied by 720, a flow meter is used for measuring flow field information near a building, and repeated experiments are adopted to ensure the measurement accuracy of impact load;
s13, carrying out a model test: the gate is closed, the water tank is stored to the water depths of 0.6m, 0.9m and 1.2m at the first section, after the water level is stable, the gate control switch is immediately started, the bottom of the gate is quickly turned over under the lifting of gravity, water flow quickly impacts a house model, meanwhile, each measuring device captures and records flow field data, and tail water flows into the reservoir again through the channel;
the numerical simulation step in the step S2 comprises the following steps:
s21, adopting three groups of test schemes for numerical simulation;
scheme one: according to the scheme that the transverse shafts of the buildings are arranged vertically to the incoming flow direction, 3 house models are arranged at positions 3m, 5m and 7m away from the gate, impact mechanical characteristics of the water flow direct-impact building acting on the building and protection effects of the water flow direct-impact building on subsequent buildings are discussed, 3 groups of comparison tests are respectively arranged, the house models are independently arranged at positions 3m, 5m and 7m away from the gate, and the deformation characteristics of the building under the effect of flood impact of the same magnitude under the effect of no front-row building protection are verified;
scheme II: arranging the horizontal axis of the building in parallel with the incoming flow direction, and setting 3 groups of comparison tests to study the shielding effect of the front row of buildings on the rear row of buildings under the impact of different water heads;
scheme III: the arrangement mode of the building group has an influence on the mechanical property of water flow impact, and the blocking effect of the arrangement mode of the building on the water flow is researched by adopting two arrangement modes from physical and mechanical characteristics, so that references are provided for optimizing the arrangement mode of the downstream building of the dam, improving the response of the downstream building to the flood impact and maximally guaranteeing the safety of the building;
s22, constructing a flood impact numerical model
The Reynolds average Navier-Stokes equation and the RNG k- ε turbulence model are used to solve the over-dam water flow motion and turbulence, and the continuous equation is as follows:
wherein u is i Represents the average speed x i Is the dimension, t is the time, p is the pressure, ρ is the fluid density, g i Is the gravity acceleration component, v is the molecular kinematic viscosity, v t Is vortex motion viscosity;
the reynolds average Navier-Stokes equation is a control equation for flow field average variables consisting of a time average and a pulsation, by introducing the Boussinesq assumption that the reynolds stress of turbulence is proportional to strain, the turbulence calculation is attributed to the calculation of the proportionality coefficient between reynolds stress and strain, the turbulence viscosity coefficient, by using the vortex motion viscosity calculated by the turbulence energy k and turbulence dissipation rate epsilon of the form:
wherein C is μ As a constant, the equation can not be closed due to the introduction of Reynolds stress after homogenization in a flow field, namely, the number of unknown numbers is more than the number of equation, the RNGk-epsilon turbulence model adopts the corrected turbulence dissipation rate to represent the turbulence intensity, and the transportation equations of k and epsilon are as follows:
wherein R is ε For turbulent shear performance, G is the rate of turbulent kinetic energy, and the coefficients in the model are as follows:
C μ ＝0.085，C ε1 ＝1.42，C ε2 ＝1.68，σ k ＝σ ε = 0.7194, β=0.012, η=4.38, these coefficients are obtained by data fusion using statistical data,
the VOF method is used to solve for free surface motion, and tracking of gas and liquid free surfaces is solved by continuous equations of the form:
any finite control body being filled with gas or liquid or a mixture of both, depending on the volume fraction alpha of the liquid w The volume fraction within a given control volume is a constant value;
the time step is determined from the CFL conditions in the three-dimensional model having the form:
wherein, C is a dimensionless constant, a display time pushing format is adopted, and the value of C is set to be 1;
s23, solving structural deformation;
s231, solid phase control equation
Balance equation:
geometric equation:
the physical equation: sigma (sigma) ij ＝λσ ij ε kk +2με ij
Wherein k is i Is dynamic pressure load component, and has the unit of m/s 2 ；ε ij Is Cauchy strain constant, and has the unit of s -1 The method comprises the steps of carrying out a first treatment on the surface of the Lambda, mu are Lame constants of the carbon steel bent pipes respectively; e is the elastic modulus of the protective film, and the unit is kg.m -1 ·s -2 Omega is poisson's ratio;
s232, a fluid-solid coupling control equation;
the fluid-solid coupling interface satisfies the following 4 control equations:
τ w ·n w ＝τ s ·n s
d w ＝d s
T w ＝T s
q w ＝q s
wherein the subscript w represents fluid, the subscript s represents solid, τ represents solid stress, d represents displacement, T represents temperature, and q represents heat flow.
2. The building group energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test according to claim 1, wherein the method is characterized by comprising the following steps of: and the step S1 of the physical model experiment operation is that the organic glass house experiment model in the step S11 is a single-bay building, so that the flood impact time course and the dynamic condition of the impact process can be accurately captured.
3. The building group energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test according to claim 1, wherein the method is characterized by comprising the following steps of: and (3) in the step S11, the bottom slope of the water tank is 1%, and the impact effect of dam-break flood is simulated.
4. The building group energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test according to claim 1, wherein the method is characterized by comprising the following steps of: and the step S1 of the physical model experiment operation step is that the opening speed of the gate in the step S11 is controlled by a counterweight, so that the gate can be opened instantaneously, the process of impacting a house by a large amount of flood instantaneous gushes after dam-break flood is fully simulated, and the consistency of a model experiment and actual conditions is ensured.
5. The building group energy dissipation and impact prevention optimization design method based on combination of numerical simulation and physical model test according to claim 1, wherein the method is characterized by comprising the following steps of: the control test in the step S21 is the same as the physical model test, and is used for verifying the accuracy of the numerical simulation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911124882.1A CN111008418B (en) | 2019-11-18 | 2019-11-18 | Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911124882.1A CN111008418B (en) | 2019-11-18 | 2019-11-18 | Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111008418A CN111008418A (en) | 2020-04-14 |
| CN111008418B true CN111008418B (en) | 2023-06-06 |
Family
ID=70112854
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911124882.1A Active CN111008418B (en) | 2019-11-18 | 2019-11-18 | Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111008418B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117874903B (en) * | 2024-03-13 | 2024-05-28 | 江苏省建筑工程质量检测中心有限公司 | Engineering quality assessment diagnosis method and system integrating Internet of things and artificial intelligence |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101017516A (en) * | 2007-02-12 | 2007-08-15 | 三峡大学 | Water flow elaborate simulation method of hydraulic of hydraulic engineering |
| CN102663827A (en) * | 2012-03-02 | 2012-09-12 | 天津大学 | Three-dimensional dynamic whole-process simulation method for storm surge and flood routing in complex flooding areas |
| CN104156525A (en) * | 2014-08-05 | 2014-11-19 | 天津大学 | Method for improving storm surge disaster risk predication accuracy |
| CN110083866A (en) * | 2019-03-26 | 2019-08-02 | 中国核电工程有限公司 | A kind of simulation system and method for interior of building water logging |
-
2019
- 2019-11-18 CN CN201911124882.1A patent/CN111008418B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101017516A (en) * | 2007-02-12 | 2007-08-15 | 三峡大学 | Water flow elaborate simulation method of hydraulic of hydraulic engineering |
| CN102663827A (en) * | 2012-03-02 | 2012-09-12 | 天津大学 | Three-dimensional dynamic whole-process simulation method for storm surge and flood routing in complex flooding areas |
| CN104156525A (en) * | 2014-08-05 | 2014-11-19 | 天津大学 | Method for improving storm surge disaster risk predication accuracy |
| CN110083866A (en) * | 2019-03-26 | 2019-08-02 | 中国核电工程有限公司 | A kind of simulation system and method for interior of building water logging |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111008418A (en) | 2020-04-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mason et al. | Pulsed wall jet simulation of a stationary thunderstorm downburst, Part A: Physical structure and flow field characterization | |
| Hsiao et al. | On the evolution and run-up of breaking solitary waves on a mild sloping beach | |
| Lin et al. | A numerical study of breaking waves in the surf zone | |
| Ma et al. | Shock-capturing non-hydrostatic model for fully dispersive surface wave processes | |
| Wood et al. | Physical and numerical modelling of thunderstorm downbursts | |
| Cheng et al. | Numerical simulation of air-water two-phase flow over stepped spillways | |
| Ko et al. | Hydraulic experiments on impact forces from tsunami-driven debris | |
| Tuan et al. | A numerical model of wave overtopping on seadikes | |
| Habibzadeh et al. | Performance of baffle blocks in submerged hydraulic jumps | |
| Qu et al. | Numerical investigation on hydrodynamic load of coastal bridge deck under joint action of solitary wave and wind | |
| Kihara et al. | A physical model of tsunami inundation and wave pressures for an idealized coastal industrial site | |
| CN111008418B (en) | Building group energy dissipation anti-impact optimization design method based on combination of numerical simulation and physical model test | |
| Marino et al. | Influence of wind–waves energy transfer on the impulsive hydrodynamic loads acting on offshore wind turbines | |
| Peng et al. | Numerical investigation of the impact of a dam-break induced flood on a structure | |
| Duarte et al. | An enhanced physically based scour model for considering jet air entrainment | |
| Ma et al. | Effects of unstable thermal stratification on the flow characteristics in an idealized rural-to-urban transition region: A large-eddy simulation study | |
| Robb et al. | Numerical simulation of dam-break flows using depth-averaged hydrodynamic and three-dimensional CFD models | |
| Dutta et al. | Large eddy simulation (LES) of flow and bedload transport at an idealized 90-degree diversion: Insight into Bulle effect | |
| Minor et al. | " Barbs" for river bend bank protection: application of a three-dimensional numerical model | |
| Lin et al. | Large eddy simulations of strong wind mechanisms at pedestrian level around square-section buildings with same aspect ratios and different sizes | |
| Meselhe et al. | Numerical modeling for fish diversion studies | |
| Li et al. | Lagrangian modeling of tidal bores passing through bridge piers | |
| CN110057488B (en) | Method for calculating maximum pressure intensity and vertical distribution of tidal bore impact cylindrical pile | |
| Wisse et al. | Effects of joint opening and block protrusion on the hydraulic parameters affecting rock block erosion in unlined spillways using a reduced-scale model | |
| Wang et al. | Hydraulic jumps and breaking bores: modelling and analysis |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |