CN108846172B - Mud-rock-flow-impact load function generation method and device - Google Patents
Mud-rock-flow-impact load function generation method and device Download PDFInfo
- Publication number
- CN108846172B CN108846172B CN201810527748.5A CN201810527748A CN108846172B CN 108846172 B CN108846172 B CN 108846172B CN 201810527748 A CN201810527748 A CN 201810527748A CN 108846172 B CN108846172 B CN 108846172B
- Authority
- CN
- China
- Prior art keywords
- mud
- rock
- flow
- impact load
- amplitude
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 62
- 239000011435 rock Substances 0.000 claims abstract description 72
- 230000001133 acceleration Effects 0.000 claims description 24
- 230000005484 gravity Effects 0.000 claims description 23
- 240000002853 Nelumbo nucifera Species 0.000 claims description 13
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 13
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 13
- 238000004422 calculation algorithm Methods 0.000 claims description 13
- 210000001367 artery Anatomy 0.000 claims description 2
- 230000002045 lasting effect Effects 0.000 claims description 2
- 210000003462 vein Anatomy 0.000 claims description 2
- 241000251468 Actinopterygii Species 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 58
- 230000000875 corresponding effect Effects 0.000 description 21
- 230000008569 process Effects 0.000 description 14
- FQZYTYWMLGAPFJ-OQKDUQJOSA-N tamoxifen citrate Chemical compound [H+].[H+].[H+].[O-]C(=O)CC(O)(CC([O-])=O)C([O-])=O.C=1C=CC=CC=1C(/CC)=C(C=1C=CC(OCCN(C)C)=CC=1)/C1=CC=CC=C1 FQZYTYWMLGAPFJ-OQKDUQJOSA-N 0.000 description 14
- 238000010586 diagram Methods 0.000 description 12
- 230000003068 static effect Effects 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 5
- 239000012530 fluid Substances 0.000 description 5
- 238000013178 mathematical model Methods 0.000 description 5
- 238000004080 punching Methods 0.000 description 5
- 230000004044 response Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000033001 locomotion Effects 0.000 description 4
- 238000005457 optimization Methods 0.000 description 4
- 230000002093 peripheral effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 3
- 239000007788 liquid Substances 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005253 cladding Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 2
- 238000004590 computer program Methods 0.000 description 2
- 230000002596 correlated effect Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012797 qualification Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 239000004575 stone Substances 0.000 description 2
- 238000013316 zoning Methods 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 244000131316 Panax pseudoginseng Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003698 anagen phase Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000009514 concussion Effects 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000012010 growth Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000010348 incorporation Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000013439 planning Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000008439 repair process Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000002889 sympathetic effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The embodiment of the invention provides a kind of mud-rock-flow-impact load function generation method and devices, are related to geotechnical engineering monitoring field.The method includes the actual information based on the target mud-rock flow case got, after the multiple associated eigenvalues for determining mud-rock-flow-impact load, it is based on the multiple associated eigenvalue, determines multiple undetermined coefficients in mud-rock-flow-impact Load Model;The ripple frequency of mud-rock-flow-impact load is obtained again, obtains the duration of mud-rock-flow-impact load and obtains the overburden depth range of mud-rock-flow-impact load;Then the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the overburden depth range are substituted into the mud-rock-flow-impact Load Model, generates mud-rock-flow-impact load function.The mud-rock-flow-impact load function for realizing generation covers the multiple elements of dynamic load, thus closer in objective reality.Effective new way is provided for simulation mud-rock-flow-impact load.
Description
Technical field
The present invention relates to geotechnical engineerings to monitor field, in particular to a kind of mud-rock-flow-impact load function generation side
Method and device.
Background technique
Mud-rock flow is one of the multiple geological disaster in Southwestern China area.The impulsive force of mud flow is mud-rock flow dynamic characteristics
The significant design foundation of the means of transportation such as important component and mud-rock flow branch mark structure, tunnel, bridge, culvert.Due to each
The limitation of kind of condition, the actual measurement impulsive force of mud flow be it is very difficult, this is also that mud-rock flow is prevented and reduced natural disasters one of the difficult point of work,
Especially for still nonevent mud-rock flow in planning prevention and treatment, it is unpractical for obtaining impact force by the means of actual measurement.One
Feasible scheme is to propose to be applicable on the basis of related physical simulation test according to the observational record for having mud-rock flow case
, simplify computation model.The experts and scholars of computation model in relation to the impulsive force of mud flow, countries in the world set out from different perspectives
Establish a variety of impulsive force of mud flow computation models of power.It can be divided into from the angle of technical method: static(al) fluid model, motive fluid mould
Type, momentum-momentum model, energy model etc.;If according to whether the angle for mutually separately considering solid, liquid, and can be divided into one-dimensional
Fluid model and two dimensional fluid model.
As typical two-phase solid-liquid fluid, the impact process of mud-rock flow has significant pulse characteristic, is substantially
One typical dynamics problem.Obvious drawback existing for above-mentioned conventional impulsive force of mud flow computation model is: using quasi-
While static manner carries out highly simplified to the problem of script dynamic, it is also inevitably lost mud-rock flow physical mechanics mistake
The dynamics details of journey.And these details either study itself for abundant dynamics of debris flow, or deeply promote ground
Engineering minute design is all vital.That is, existing solution is simply, rambunctiously with static
Mode avoids dynamic response analysis.To the shortcoming that mud-rock-flow-impact load is held, necessarily cause subsequent geotechnical structure design knot from
The requirement of fining greatly differs from each other, fruit or overly conservative, there is design waste or safety stock is insufficient, there are engineering wind
Danger.
Although these versions are various, solve the problems, such as a common mode is: using static viewpoint come
The problem for the treatment of script dynamic assesses the maximum impact load or average load of mud-rock flow in a manner of Quintic system, and with this
As mud-rock flow supporting file engineering, protection works and other ground, the load design foundation of Structural Engineering.Therefore, this to power
The tupe of problem is coarse.
Summary of the invention
The purpose of the present invention is to provide a kind of mud-rock-flow-impact load function generation method and devices, to improve above-mentioned ask
Topic.To achieve the goals above, the technical solution adopted by the present invention is as follows:
In a first aspect, the embodiment of the invention provides a kind of mud-rock-flow-impact load function generation method, the method packet
The actual information based on the target mud-rock flow case got is included, determines multiple associated eigenvalues of mud-rock-flow-impact load;Base
In the multiple associated eigenvalue, multiple undetermined coefficients in mud-rock-flow-impact Load Model are determined;Obtain mud-rock-flow-impact lotus
The ripple frequency of load;Obtain the duration of mud-rock-flow-impact load;Obtain the overburden depth range of mud-rock-flow-impact load;It will
The multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the overburden depth model
It encloses and substitutes into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load function.
Second aspect, the embodiment of the invention provides a kind of mud-rock-flow-impact load function generating means, described device packets
Include: the first determination unit, the second determination unit, first acquisition unit, second acquisition unit, third acquiring unit and generation are single
Member.First determination unit determines mud-rock-flow-impact load for the actual information based on the target mud-rock flow case got
Multiple associated eigenvalues.Second determination unit determines mud-rock-flow-impact Load Model for being based on the multiple associated eigenvalue
In multiple undetermined coefficients.First acquisition unit, for obtaining the ripple frequency of mud-rock-flow-impact load.Second acquisition unit,
For obtaining the duration of mud-rock-flow-impact load.Third acquiring unit, the covering for obtaining mud-rock-flow-impact load are deep
Spend range.Generation unit, for by the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, described hold
Continuous time and the overburden depth range substitute into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load function.
The embodiment of the invention provides a kind of mud-rock-flow-impact load function generation method and devices, and the method includes bases
After the actual information of the target mud-rock flow case got, the multiple associated eigenvalues for determining mud-rock-flow-impact load, it is based on
The multiple associated eigenvalue determines multiple undetermined coefficients in mud-rock-flow-impact Load Model;Mud-rock-flow-impact lotus is obtained again
The ripple frequency of load, the duration for obtaining mud-rock-flow-impact load and the overburden depth model for obtaining mud-rock-flow-impact load
It encloses;Then by the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, duration and described
Overburden depth range substitutes into the mud-rock-flow-impact Load Model, generates mud-rock-flow-impact load function.Realize generation
Mud-rock-flow-impact load function covers the multiple elements of dynamic load, thus closer in objective reality.To simulate mud-rock-flow-impact
Load provides effective new way.
Other features and advantages of the present invention will be illustrated in subsequent specification, also, partly be become from specification
It is clear that by implementing understanding of the embodiment of the present invention.The objectives and other advantages of the invention can be by written theory
Specifically noted structure is achieved and obtained in bright book, claims and attached drawing.
Detailed description of the invention
In order to illustrate the technical solution of the embodiments of the present invention more clearly, below will be to needed in the embodiment attached
Figure is briefly described, it should be understood that the following drawings illustrates only certain embodiments of the present invention, therefore is not construed as pair
The restriction of range for those of ordinary skill in the art without creative efforts, can also be according to this
A little attached drawings obtain other relevant attached drawings.
Fig. 1 is a kind of structural block diagram that can be applied to electronic equipment provided in an embodiment of the present invention;
Fig. 2 is the flow chart of mud-rock-flow-impact load function generation method provided in an embodiment of the present invention;
Fig. 3 (a) is that lotus is impacted in mud-rock flow vertical section during mud-rock-flow-impact load function provided in an embodiment of the present invention generates
Carry theoretical generalization schematic diagram;
Fig. 3 (b) is that lotus is impacted in mud-rock flow cross section during mud-rock-flow-impact load function provided in an embodiment of the present invention generates
Carry theoretical generalization schematic diagram;
Fig. 4 is that mud-rock-flow-impact load acceleration of gravity additional coefficient k (t) changes over time waveform in the embodiment of the present invention
Example schematic diagram;
Fig. 5 is in the embodiment of the present invention so that in the expressway network project of Guizhou, mud-rock flow is evaluated as case by route
And the mud-rock-flow-impact load acceleration of gravity additional coefficient k (t) generated changes over time waveform;
Fig. 6 is the structural block diagram of mud-rock-flow-impact load function generating means provided in an embodiment of the present invention.
Specific embodiment
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.The present invention being usually described and illustrated herein in the accompanying drawings is implemented
The component of example can be arranged and be designed with a variety of different configurations.Therefore, below to the reality of the invention provided in the accompanying drawings
The detailed description for applying example is not intended to limit the range of claimed invention, but is merely representative of selected implementation of the invention
Example.Based on the embodiments of the present invention, obtained by those of ordinary skill in the art without making creative efforts
Every other embodiment, shall fall within the protection scope of the present invention.
It should also be noted that similar label and letter indicate similar terms in following attached drawing, therefore, once a certain Xiang Yi
It is defined in a attached drawing, does not then need that it is further defined and explained in subsequent attached drawing.Meanwhile of the invention
In description, term " first ", " second " etc. are only used for distinguishing description, are not understood to indicate or imply relative importance.
Fig. 1 shows a kind of structural block diagram of electronic equipment 100 that can be applied in the embodiment of the present invention.As shown in Figure 1,
Electronic equipment 100 may include memory 102, storage control 104, one or more (one is only shown in Fig. 1) processors
106, Peripheral Interface 108, input/output module 110, audio-frequency module 112, display module 114, radio-frequency module 116 and mud-rock flow punching
Hit load function generating means.
Memory 102, storage control 104, processor 106, Peripheral Interface 108, input/output module 110, audio mould
Block 112, display module 114 are directly or indirectly electrically connected between each element of radio-frequency module 116, with realize data transmission or
Interaction.For example, can realize electrical connection by one or more communication bus or signal bus between these elements.Mud-rock flow punching
Storage can be stored in the form of software or firmware (firmware) by hitting load function generation method and respectively including at least one
Software function module in device 102, for example, the mud-rock-flow-impact load function generating means software function module that includes or
Computer program.
Memory 102 can store various software programs and module, such as mud-rock-flow-impact provided by the embodiments of the present application
Load function generation method and the corresponding program instruction/module of device.Processor 106 is stored in the memory 102 by operation
Software program and module, the mud thereby executing various function application and data processing, i.e., in realization the embodiment of the present application
Rock glacier impact load function generation method.
Memory 102 can include but is not limited to random access memory (Random Access Memory, RAM), only
It reads memory (Read Only Memory, ROM), programmable read only memory (Programmable Read-Only
Memory, PROM), erasable read-only memory (Erasable Programmable Read-Only Memory, EPROM),
Electricallyerasable ROM (EEROM) (Electric Erasable Programmable Read-Only Memory, EEPROM) etc..
Processor 106 can be a kind of IC chip, have signal handling capacity.Above-mentioned processor can be general
Processor, including central processing unit (Central Processing Unit, abbreviation CPU), network processing unit (Network
Processor, abbreviation NP) etc.;It can also be digital signal processor (DSP), specific integrated circuit (ASIC), field-programmable
Gate array (FPGA) either other programmable logic device, discrete gate or transistor logic, discrete hardware components.It can
To realize or execute disclosed each method, step and the logic diagram in the embodiment of the present application.General processor can be micro-
Processor or the processor are also possible to any conventional processor etc..
Various input/output devices are couple processor 106 and memory 102 by the Peripheral Interface 108.Some
In embodiment, Peripheral Interface 108, processor 106 and storage control 104 can be realized in one single chip.Other one
In a little examples, they can be realized by independent chip respectively.
Input/output module 110 is used to be supplied to the interaction that user input data realizes user and electronic equipment 100.It is described
Input/output module 110 may be, but not limited to, mouse and keyboard etc..
Audio-frequency module 112 provides a user audio interface, may include one or more microphones, one or more raises
Sound device and voicefrequency circuit.
Display module 114 provides an interactive interface (such as user interface) between electronic equipment 100 and user
Or it is referred to for display image data to user.In the present embodiment, the display module 114 can be liquid crystal display or touching
Control display.It can be the touching of the capacitance type touch control screen or resistance-type of support single-point and multi-point touch operation if touch control display
Control screen etc..Single-point and multi-point touch operation is supported to refer to that touch control display can sense on the touch control display one or more
The touch control operation generated simultaneously at a position, and the touch control operation that this is sensed transfers to processor 106 to be calculated and handled.
Radio-frequency module 116 is used to receive and transmit electromagnetic wave, realizes the mutual conversion of electromagnetic wave and electric signal, thus with
Communication network or other equipment are communicated.
It is appreciated that structure shown in FIG. 1 is only to illustrate, electronic equipment 100 may also include it is more than shown in Fig. 1 or
Less component, or with the configuration different from shown in Fig. 1.Each component shown in Fig. 1 can using hardware, software or its
Combination is realized.
In the embodiment of the present invention, electronic equipment 100 can be used as user terminal, or as server.User terminal
It can be PC (personal computer) computer, tablet computer, mobile phone, laptop, smart television, set-top box, vehicle-mounted
The terminal devices such as terminal.
Referring to Fig. 2, the embodiment of the invention provides a kind of mud-rock-flow-impact load function generation method, the method packet
It includes: step S200, step S210, step S220, step S230, step S240 and step S250.
Step S200: the actual information based on the target mud-rock flow case got determines the more of mud-rock-flow-impact load
A associated eigenvalue.
Theoretical principle is described as follows:
The mathematical model of mud-rock-flow-impact load waveform generates the mathematical formulae system of mud-rock-flow-impact load, be this
The core and theoretical basis of application, are described below:
P (t, h)=k (t) ρ gh (1)
Wherein, P (t, h) indicates the mud-rock-flow-impact pressure function that t and depth location h changes at any time.Such as Fig. 3 (a) and
Fig. 3 (b), h are the depth that mud-rock flow covers under current location, and variation range, i.e. depth coverage area are 0~hmax.ρ is mud
Rock glacier density, g are acceleration of gravity, and k (t) is acceleration of gravity additional coefficient, are the functions an of variation of t at any time, by
(2) formula determines:
The effect of above-mentioned (2) formula is: when being negative by w (t) waveform, zero is taken, to meet mud-rock flow and Rock And Soil side
Boundary's contact surface only applies pressure, does not apply the restrictive condition of pulling force.Wherein, w (t) is the continually changing sinusoidal waveform letter of an amplitude
Number, is determined as follows by (3) formula:
Wherein, α1, α2, β1, β2, γ1, γ2, k, b are waveform control parameter, and R is residual value, and e is the bottom of natural logrithm, and π is
Pi, f are the vibration frequency of mud-rock-flow-impact wave,WithAbove-mentioned point is controlled respectively
The size of section function waveform amplitude, i.e. some moment tiCorresponding amplitude AiIt is the coefficient entry of SIN function in formula (3), i.e.,
Have:
On figure, point (ti,Ai) be function w (t) waveform envelope line on point.
Above-mentioned function has the property that
(4) formula of analysis, to the amplitude of piecewise functionWithFirst derivative is sought respectively,
And enabling it is zero, the maximum value for obtaining amplitude (is denoted as Amax) and maximum value corresponding time, that is, maximum value arrival time (be denoted as
tAmax), since piecewise function is equal at amplitude maximum, there is following formula establishment:
It is zero that the necessary condition of point of inflexion on a curve (the concave-convex separation of curve), which is its second dervative,.Therefore, respectively to point
The amplitude of section functionWithSecond dervative is sought respectively, and enabling it is zero, obtains amplitude
Two inflection points, are located at tAmaxThe left and right sides of (i.e. peak swing reaches the time), is denoted as tLAnd tr。
According to the property of inflection point: (0, tL) in section, amplitude curve fovea superior, amplitude from scratch, is at this time
The rapid growth phase;In (tL,tAmax) in section, amplitude curve is convex, at this time amplitude slowdown in growth, and progressivelyes reach maximum value;?
(tAmax,tr) in section, amplitude curve is convex, and amplitude is quickly fallen after rise by maximum value at this time;In (tr, ∞) and in section, amplitude is fallen after rise
Speed gradually slows down, and is finally intended to residual value R.
In conclusion the impact load mathematical model established by formula (1)~(3).By giving α1, α2, β1, β2, γ1,
γ2, k, b, these undetermined coefficients of R can determine waveform.As an example, table 1 gives one group of parameter, corresponding waveform is such as
Shown in Fig. 4, waveform technology index is shown in Table 2.
1 mud-rock-flow-impact load equation parameter value example of table
The Specifications example of 2 mud-rock-flow-impact load equation of table
Based on step S200, further, under normal circumstances, known to the impact load associated eigenvalue conduct of target case
Amount provides.Multiple associated eigenvalues include: the maximum value A of mud-rock-flow-impact load amplitudemax, maximum value arrival time tAmax, vibration
A little corresponding time of turning left of the time i.e. amplitude of turning left of width, amplitude the time t that turns rightrI.e. amplitude turn right it is a little corresponding
Time, amplitude the time t that turns rightR, be greater than the time t that turns right on impact load residual value R and impact load tail portionR's
One moment corresponding amplitude A1.Detailed description are as follows for associated eigenvalue:
Peak swing AmaxPhysical meaning be: mud-rock flow bring impact force maximum can be equivalent to debris flow body at this
The multiple of depth location static state pressure, that is, the attainable maximum value of k (t) institute, are denoted as k (t)max。Lichtenhalhn
(1973) suggest k (t)maxIt takes between 2.8~4.4, Armanini (1997) then thinks k (t)maxIt can get between 2.8~5;
Scotton and Deganutti (1997) pass through experimental study, it is believed that k (t)maxValue range 2.7~7.5.The application is built
View value, specific value between 2.8~8.0 can be provided according to the specific characteristic of mudstone flow object investigated.
Mud-rock-flow-impact load amplitude maximum value arrival time tAmax, this time has very big difference according to different cases
It is different.Specific to a certain specific mud-rock-flow-impact overall process, the overall process of mud-rock-flow-impact may continue tens of seconds to hundreds of seconds
Deng, and the time for reaching peak value is the several seconds to tens of seconds etc..Specific value needs determine than selecting according to real case.
The time t that turns left of amplitudeL, turn right a time trRespectively in tAmaxThe left and right sides (i.e. tL<tAmax, tr>tAmax)。
Impact load residual value R is the size for controlling impact load tail portion pinching value.For example, as R=1.0, formula
(3) in the tail portion w (t) gradually pinching to 1.0.1.0 are generally taken for typical mud-rock-flow-impact load.
The a certain moment t in impact load tail portion1(it is required that t1>tr) and its corresponding amplitude A1.Introduce t1And A1It is to control
The convergence form of impact load tail portion, for example, converging on impact load residual value R then desirable t as far as possible to be limited in 18 seconds1=
18 seconds, A1The desirable value more slightly larger than 0.
Step S210: it is based on the multiple associated eigenvalue, determines multiple systems undetermined in mud-rock-flow-impact Load Model
Number.
Further, the multiple associated eigenvalue includes the maximum value of mud-rock-flow-impact load amplitude, maximum value arrival
Greater than described on time, the time of turning left of amplitude, the time of turning right of amplitude, impact load residual value and impact load tail portion
It turns right a moment corresponding amplitude of time, step S210 includes:
Maximum value, maximum value arrival time and amplitude based on the mud-rock-flow-impact load amplitude are turned left the time,
Mud-rock-flow-impact Load Model is solved by solution methods of nonlinear equations, is determined in the multiple undetermined coefficient according to solving result
A part of coefficient;
The time of turning right of maximum value, maximum value arrival time, amplitude based on the mud-rock-flow-impact load amplitude, punching
It hits on load residual value and impact load tail portion in a moment corresponding amplitude for being greater than the time of turning right, passes through group's intelligence
Energy algorithm solves mud-rock-flow-impact Load Model, and another part system in the multiple undetermined coefficient is determined according to solving result
Number.
A part of coefficient in the multiple undetermined coefficient is determined according to solving result, comprising:
It is based onDetermine a part of coefficient in the multiple undetermined coefficient, AmaxFor the mud-rock flow
The maximum value of impact load amplitude, tAmaxFor maximum value arrival time, tLFor turning left the time for amplitude, α1、β1、γ1It is described
A part of coefficient.
Another part coefficient in the multiple undetermined coefficient is determined according to solving result, comprising:
It is based onDetermine another part in the multiple undetermined coefficient
Coefficient, AmaxFor the maximum value of the mud-rock-flow-impact load amplitude, tAmaxFor maximum value arrival time, trFor turning right a little for amplitude
Time, R are impact load residual value, t1For on impact load tail portion in the moment for being greater than the time of turning right, A1To impact lotus
It carries on tail portion in the moment t for being greater than the time of turning right1Corresponding amplitude, α2, β2, γ2, k, b are described another part
Coefficient.
Specifically, it is determined that the purpose of undetermined coefficient, seeks to the specific impact load number for establishing formula (1)~(3)
It learns model and has specified correlated characteristic, to be realized with the language of mathematics to specific mud-rock-flow-impact load a certain in reality
The description of case.The method for solving undetermined coefficient has very much, mainly uses two methods herein: the first is using non-linear
Solving equations undetermined coefficient, second is to regard the solution of undetermined coefficient as nonlinear optimal problem, is calculated with swarm intelligence
Method solves.Steps are as follows:
Using solution methods of nonlinear equations, α in (3) formula is determined1, β1, γ1These three unknown numbers.
The maximum value A of known mud-rock-flow-impact load amplitudemaxAnd maximum value arrival time tAmax, amplitude turns left a little
tL, following equation is constructed according to (5)~(7) formula, as shown in (9) formula.
Above-mentioned equation group has α1, β1, γ1Three unknown numbers are made of three equations, and analytic solutions are as follows:
Using Swarm Intelligence Algorithm, α in (3) formula is determined2, β2, γ2, this five unknown numbers of k, b.Known mud-rock-flow-impact lotus
Carry the maximum value A of amplitudemaxAnd maximum value arrival time tAmax, amplitude the point t that turns rightr, impact load residual value R, impact load
The a certain moment t in tail portion1(meet t1>tr) when corresponding amplitude A1.Following equation is constructed according to (5), (6), (8), (4) formula, such as
(11) shown in formula.
Above-mentioned equation group has α2, β2, γ2, five unknown numbers of k, b, but only there are four equations to constitute, therefore there are more solutions.Often
The solution methods of nonlinear equations of rule cannot obtain reasonable numerical solution if not newly-increased qualifications.However, from nonlinear optimization
From the perspective of, (11) formula is actually a constrained optimization problem (constraint condition k > 0).In general, to this
The constrained optimization problem of class, uses the Swarm Intelligence Algorithm using particle swarm algorithm as representative, obtains one group with can be convenient
Or multiple groups numerical solution, without artificially increasing qualifications.Related particle swarm algorithm and other Swarm Intelligence Algorithm solution sides
The correlative detail of journey group is general technology, is not developed in details herein.
In the embodiment of the present application, the Swarm Intelligence Algorithm may include, but be not limited to population, ant colony, Artificial Fish
Group or bee colony.In general, regardless of which method is applied, can be obtained to this kind of constrained optimization problem
Approximate several groups numerical solution, due to theoretical framework difference, numerical value is had a certain difference.Above-mentioned any one swarm intelligence
Algorithm is applicable in the application.
Step S220: the ripple frequency of mud-rock-flow-impact load is obtained.
The ripple frequency f of mud-rock-flow-impact load i.e. the reference frequency of dynamic load.The selection and earthquake of the frequency
Dynamic Site characteristic is related with mud-rock flow inherent characteristic.It is generally believed that the excellent frequency of Rock And Soil is lower, it is distributed in the left side 2-10Hz
It is right.It is inverted to eigenperiod (predominant period) to be see also data such as " Earthquake In China characteristic periods of response spectra zoning maps "
Excellent frequency can be obtained, and using this frequency as the reference frequency of mud-rock flow dynamic load.In view of high-frequency energy decaying faster, and
The harm of low frequency oscillation is bigger, therefore should take low frequency, the mud-rock-flow-impact of the application as far as possible in the case where making least favorable analysis
The ripple frequency of load can be chosen in 2-10Hz or so.
Step S230: the duration of mud-rock-flow-impact load is obtained.
S200- step S220 through the above steps, acceleration of gravity additional coefficient k (t) are it has been determined that the Function feature
It is that over time, gradually pinching converges to definite value R for tail portion.Although can theoretically be selected in the selection of attack time
0~infinity is selected, but in practical kinematic analysis calculating, both can be time-consuming, also unnecessarily.K (t) figure can be made, when function tail
Portion obviously restrains, and the time at this is taken to be truncated.When duration obviously restrains corresponding to k (t) function tail portion for 0
It carves.
Step S240: the overburden depth range of mud-rock-flow-impact load is obtained.
Determine the density p of mud-rock flow, gravity acceleration g, depth h variation range.The bulk density and water content of mud-rock flow has
It closes, can be determined according to specific object.For example, block stone density takes 2650kg/m3, water density takes 1000kg/m3, mudstone current density
Consider according to 70% block of stone and 30% water, then mud-rock flow density p is 2155kg/m3.Gravity acceleration g takes 9.8m/s2.
Step S250: by the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, described lasting
Time and the overburden depth range substitute into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load function.
The mud-rock-flow-impact Load Model is P (t, h)=k (t) ρ gh, and Step S250 includes:
By the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, duration and described
Overburden depth range substitutes into P (t, h)=k (t) ρ gh, and Generate mud-rock-flow-impact load function, P
(t, h) indicates that the mud-rock-flow-impact load pressure functional value that t and depth location h changes at any time, h are mud-rock flow under current location
Overburden depth, ρ is the mudstone current density obtained in advance, and g is acceleration of gravity, and k (t) is acceleration of gravity additional coefficient, w
It (t) is the continually changing sinusoidal waveform function of an amplitude, α1, α2, β1, β2, γ1, γ2, k, b are the multiple undetermined coefficient, and R is
Impact load residual value, e are the bottom of natural logrithm, and π is pi, and f is ripple frequency, and overburden depth range is the maximum value of 0 and h
Between, at the time of the duration obviously restrains corresponding for 0 to k (t) function tail portion.
S200- step S240 through the above steps has obtained all necessary parameters, substitutes them in (1)~(3) formula,
Rock glacier impact load mathematical model can be obtained.The model is the function that t and depth z changes at any time, i.e., different
There is different impact forces at moment at different mud-rock flow cladding thickness.
The mud-rock-flow-impact load waveform of generation has following several features:
First, the concept of acceleration of gravity additional coefficient peak value (maximum value) is remained, it is simultaneous with the analytical model of Quintic system
Hold.
Second, the associated eigenvalue of mud-rock-flow-impact load can be specified according to the actual conditions of target mud-rock flow case,
Then the specific waveforms met the requirements are generated.These features include: the maximum value A of mud-rock-flow-impact load amplitudemaxAnd this is most
Big value arrival time tAmax, time of turning left of amplitude, amplitude the time t that turns rightL, turn right a time tr, impact load residual value
R, a certain moment t in impact load tail portion1(meet t1>tr) when corresponding amplitude A1。
Third, meets mud-rock flow and massif contact surface only applies pressure, does not apply the restrictive condition of pulling force.
4th, power load and evanishment are compatible with static characteristic.I.e. from 0~tAmaxPower amplitude is gradually increased,
In tAmaxReach peak value, then dribbles, the acceleration of gravity of pinching to R times.R value generally takes 1, at this time 1 times of weight
Power acceleration mean mud-rock flow present stationary state, to rock mass apply pressure no longer with time correlation, only covered with mud-rock flow
Buried depth degree (thickness) is related.In the case where considering partially dangerous, R value is desirable is greater than 1, the pressure that mud-rock flow applies rock mass at this time
Greater than according to static liquid pressure consideration value.It is docked in conclusion dynamic process is realized with static(al) process.
5th, the selection of power reference frequency is compatible with Site characteristic parameter required by existing earthquake resistant code.Pass through local area
Ground motion response spectroscopic eigenvalue is can be in the hope of corresponding place excellent frequency, in this, as mud-rock flow dynamic load frequency in accordance with anti-
Shake code requirement.
It is special in the expressway network project of Guizhou in order to further illustrate beneficial effects of the present invention, by route
Mud-rock flow is evaluated as case, and illustrate mud-rock-flow-impact Load Model determines ginseng process.According to " Earthquake In China characteristic periods of response spectra
Zoning map " etc. related datas, it is known that one's respective area earthquake motion peak acceleration be 0.05g, sympathetic earthquakes basic intensity be VI degree, ground
Vibration reaction spectroscopic eigenvalue is 0.35s.The density p of mud-rock flow takes 2155kg/m3.Set such a extreme operating condition: in heavy rain
Triggering under, mud-rock flow racing current and under, climb from the lowest point, and whole process speed is exceedingly fast, and only expends tens of seconds.Mud-rock flow punching
Hit the displaying that the theoretical generalization schematic diagram of load meets Fig. 3 (a), Fig. 3 (b).
Implementation process is as follows:
Specify the correlation of mud-rock-flow-impact load special according to the actual conditions of target mud-rock flow case based on step S200
Value indicative.According to the actual conditions of analysis case, incorporation engineering experience chooses following parameter:
Choose the maximum value A of mud-rock-flow-impact load amplitudemax=7.0.To meet Scotton and Deganutti (1997)
Pass through the recommended value 2.8~8.0 of experimental study.
Choose mud-rock-flow-impact load amplitude maximum value arrival time tAmax=10 seconds.The time is shorter, and impact force is stronger
Strength selects lesser tAmax, it is intended to make the kinematic analysis under unfavorable operating condition.
Choose the time t that turns left of amplitudeL=5.0 seconds, turned right a time tr=12.0 seconds.
Choose impact load residual value R=1.0.
Choose a certain moment t in impact load tail portion1=18.0 seconds corresponding amplitude As1=0.1, it can make in this way at about 20 seconds
Period in, impact load decays to residual value R.
Above-mentioned parameter, which is chosen to conclude, is shown in Table 3.
The associated eigenvalue of mud-rock-flow-impact load in 3 embodiment of table
Mud-rock-flow-impact load wave is solved according to the correlated characteristic of specified mud-rock-flow-impact load based on step S210
Undetermined coefficient in figurate number mathematic(al) function.Steps are as follows for subitem:
Using solution methods of nonlinear equations, α in (3) formula is determined1, β1, γ1These three unknown numbers.Known mud-rock-flow-impact lotus
Carry the maximum value A of amplitudemaxAnd maximum value arrival time tAmax, amplitude the point t that turns leftL, solve as follows according to (10) formula:
4 institute of typing table is as follows.
Mud-rock-flow-impact load equation parameter value (3 parameters of first part) in 4 embodiment of table
Using Swarm Intelligence Algorithm, α in (3) formula is determined2, β2, γ2, this five unknown numbers of k, b.Known mud-rock-flow-impact lotus
Carry the maximum value A of amplitudemaxAnd maximum value arrival time tAmax, amplitude the time t that turns rightr, impact load residual value R, impact
The a certain moment t in load tail portion1(it is required that) t1>tRWhen corresponding amplitude A1.Following equation group is constructed according to (11) formula.
In the case where constraint condition is k > 0, using standard particle group's algorithm, taking population scale is 100, by 100,000 times
The one group of numerical solution obtained after iteration, mean square error are lower than 2.45e-04.It is as shown in table 5 below to retain preceding 7 effective digitals.
Mud-rock-flow-impact load equation parameter value (5 parameters of second part) in 5 embodiment of table
Based on step S220, the ripple frequency of mud-rock-flow-impact load is obtained, due to local area Ground motion response spectroscopic eigenvalue
For 0.35s, inverse is respective frequencies 2.85714Hz.Therefore, power reference frequency selection 2.857Hz meets place spy
Sign, also in accordance with earthquake resistant code requirement.
Based on step S230, the duration for obtaining mud-rock-flow-impact load is 20 seconds.According to existing literature to known mud
The observation of rock glacier case is recorded, and the dynamic impulsion process of typical mud-rock flow generally continues to differ to hundreds of seconds for tens seconds.In conjunction with
Parameter determined by the item of above-mentioned (1)~(4), acceleration of gravity additional coefficient waveform at t=20 seconds, be 1.02, obviously to
1.0 convergences, therefore selection kinematic analysis time, that is, duration is 20 seconds.
Based on step S240, determine the density p of mud-rock flow, gravity acceleration g, depth h variation range.By known item
Part, the density p of mud-rock flow are 2155kg/m3.Acceleration of gravity takes 9.8m/s2.Overburden depth range, that is, mud-rock flow burying depth h
The range of variation is 0-85m.
Based on step S250, mud-rock-flow-impact load function is generated.S200- step S240 through the above steps, has obtained
All necessary parameters substitute them in (1)~(3) formula, rock glacier impact load mathematical model can be obtained.The model is one
At any time that is, at different times the function of t and depth z variation at different mud-rock flow cladding thickness, there is different impacts
Power.
So far, dynamic load waveform is decided completely.The parameter of mud-rock-flow-impact load equation is shown in Table 3, corresponding phase
It closes technical indicator and is shown in Table 6.Impact load acceleration of gravity additional coefficient k (t) changes over time waveform and sees Fig. 5.
Mud-rock-flow-impact load equation parameter in 6 embodiment of table
It should be noted that wave function took in the application formula (3) is SIN function, cosine function is actually taken also
It can.All sin () functions can be replaced cos () function i.e. in the application, without influencing entire steps flow chart and invention effect
Fruit.
The mud-rock-flow-impact Load Model that the application proposes meets the typical module of mud-rock-flow-impact in real world, just like
Lower five features: first, remain the concept of acceleration of gravity additional coefficient peak value (maximum value), the analytical model with Quintic system
It is compatible.Second, the associated eigenvalue of mud-rock-flow-impact load can be specified according to the actual conditions of target mud-rock flow case, after
And generate the specific waveforms met the requirements.Third, meets mud-rock flow and massif contact surface only applies pressure, does not apply the limit of pulling force
Condition processed.4th, power load and evanishment are compatible with static characteristic.The mathematical model is objective to have reproduced typical mud-rock flow
Dynamic impact mode, it may be assumed that " from power concussion load be gradually increased, reach peak value, then dribble, pinching to residual value " this
One physical process of sample.Pinching means mud-rock flow and stationary state is presented, to the pressure of rock mass application to acceleration of gravity residual value
It is only related with mud-rock flow burying depth (thickness) no longer with time correlation.Therefore, under the above conditions, kinematic analysis process with
Static analysis realizes docking.5th, the selection of power reference frequency is simultaneous with Site characteristic parameter required by existing earthquake resistant code
Hold.Effective new way is provided more subtly quantitatively to portray mud-rock flow dynamic impact effect.
The embodiment of the invention provides a kind of mud-rock-flow-impact load function generation methods, and the method includes based on acquisition
The actual information of the target mud-rock flow case arrived, after the multiple associated eigenvalues for determining mud-rock-flow-impact load, based on described more
A associated eigenvalue determines multiple undetermined coefficients in mud-rock-flow-impact Load Model;The arteries and veins of mud-rock-flow-impact load is obtained again
Dynamic frequency, the duration for obtaining mud-rock-flow-impact load and the overburden depth range for obtaining mud-rock-flow-impact load;Then
By the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the overburden depth
Range substitutes into the mud-rock-flow-impact Load Model, generates mud-rock-flow-impact load function.Realize the mud-rock flow punching of generation
It hits load function and covers the multiple elements of dynamic load, thus closer in objective reality.It is provided for simulation mud-rock-flow-impact load
Effective new way.
Referring to Fig. 6, present embodiments providing a kind of mud-rock-flow-impact load function generating means 600, described device 600
It may include the first determination unit 610, the second determination unit 620, first acquisition unit 630, second acquisition unit 640, third
Acquiring unit 650 and generation unit 660.
First determination unit 610 determines that mud-rock flow rushes for the actual information based on the target mud-rock flow case got
Hit multiple associated eigenvalues of load.
Second determination unit 620 determines in mud-rock-flow-impact Load Model for being based on the multiple associated eigenvalue
Multiple undetermined coefficients.
First acquisition unit 630, for obtaining the ripple frequency of mud-rock-flow-impact load.
Second acquisition unit 640, for obtaining the duration of mud-rock-flow-impact load.
Third acquiring unit 650, for obtaining the overburden depth range of mud-rock-flow-impact load.
Generation unit 660 is used for the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, institute
It states the duration and the overburden depth range substitutes into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load letter
Number.
The multiple associated eigenvalue includes the maximum value of mud-rock-flow-impact load amplitude, maximum value arrival time, amplitude
Time of turning left, time of turning right of amplitude, on impact load residual value and impact load tail portion when being greater than described turn right
Between a moment corresponding amplitude, second determination unit 620 is used for: based on the mud-rock-flow-impact load amplitude
Maximum value, maximum value arrival time and amplitude are turned left the time, are passed through solution methods of nonlinear equations and are solved mud-rock-flow-impact lotus
Model is carried, a part of coefficient in the multiple undetermined coefficient is determined according to solving result;Based on the mud-rock-flow-impact load
Big on the maximum value of amplitude, maximum value arrival time, the time of turning right of amplitude, impact load residual value and impact load tail portion
In a moment corresponding amplitude of the time of turning right, mud-rock-flow-impact Load Model is solved by Swarm Intelligence Algorithm,
Another part coefficient in the multiple undetermined coefficient is determined according to solving result.
Second determination unit 620, for being based onDetermine one in the multiple undetermined coefficient
Part coefficient, AmaxFor the maximum value of the mud-rock-flow-impact load amplitude, tAmaxFor maximum value arrival time, tLFor a left side for amplitude
Inflection point time, α1、β1、γ1For a part of coefficient.
Second determination unit 620, for being based onDescribed in determination
Another part coefficient in multiple undetermined coefficients, AmaxFor the maximum value of the mud-rock-flow-impact load amplitude, tAmaxFor maximum value
Arrival time, trFor turning right the time for amplitude, R is impact load residual value, t1To be greater than described turn right on impact load tail portion
The moment of point time, A1For on impact load tail portion in the moment t for being greater than the time of turning right1Corresponding amplitude,
α2, β2, γ2, k, b are another part coefficient.
The mud-rock-flow-impact Load Model is P (t, h)=k (t) ρ gh, and Generation unit 660, being used for will be the multiple
Associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the overburden depth range substitute into P
(t, h)=k (t) ρ gh, and Generate mud-rock-flow-impact load function, P
(t, h) indicates that the mud-rock-flow-impact load pressure functional value that t and depth location h changes at any time, h are mud-rock flow under current location
Overburden depth, ρ is the mudstone current density obtained in advance, and g is acceleration of gravity, and k (t) is acceleration of gravity additional coefficient, w
It (t) is the continually changing sinusoidal waveform function of an amplitude, α1, α2, β1, β2, γ1, γ2, k, b are the multiple undetermined coefficient, and R is
Impact load residual value, e are the bottom of natural logrithm, and π is pi, and f is ripple frequency, and overburden depth range is the maximum value of 0 and h
Between, at the time of the duration obviously restrains corresponding for 0 to k (t) function tail portion.
The Swarm Intelligence Algorithm includes population, ant colony, artificial fish-swarm or bee colony.
The above each unit can be by software code realization, at this point, above-mentioned each unit can be stored in memory 102.
The above each unit can equally be realized by hardware such as IC chip.
The skill of mud-rock-flow-impact load function generating means 400 provided in an embodiment of the present invention, realization principle and generation
Art effect is identical with preceding method embodiment, and to briefly describe, Installation practice part does not refer to place, can refer to preceding method
Corresponding contents in embodiment.
In several embodiments provided herein, it should be understood that disclosed device and method can also pass through
Other modes are realized.The apparatus embodiments described above are merely exemplary, for example, flow chart and block diagram in attached drawing
Show the device of multiple embodiments according to the present invention, the architectural framework in the cards of method and computer program product,
Function and operation.In this regard, each box in flowchart or block diagram can represent the one of a module, section or code
Part, a part of the module, section or code, which includes that one or more is for implementing the specified logical function, to be held
Row instruction.It should also be noted that function marked in the box can also be to be different from some implementations as replacement
The sequence marked in attached drawing occurs.For example, two continuous boxes can actually be basically executed in parallel, they are sometimes
It can execute in the opposite order, this depends on the function involved.It is also noted that every in block diagram and or flow chart
The combination of box in a box and block diagram and or flow chart can use the dedicated base for executing defined function or movement
It realizes, or can realize using a combination of dedicated hardware and computer instructions in the system of hardware.
In addition, each functional module in each embodiment of the present invention can integrate one independent portion of formation together
Point, it is also possible to modules individualism, an independent part can also be integrated to form with two or more modules.
It, can be with if the function is realized and when sold or used as an independent product in the form of software function module
It is stored in a computer readable storage medium.Based on this understanding, technical solution of the present invention is substantially in other words
The part of the part that contributes to existing technology or the technical solution can be embodied in the form of software products, the meter
Calculation machine software product is stored in a storage medium, including some instructions are used so that a computer equipment (can be a
People's computer, server or network equipment etc.) it performs all or part of the steps of the method described in the various embodiments of the present invention.
And storage medium above-mentioned includes: that USB flash disk, mobile hard disk, read-only memory (ROM, Read-Only Memory), arbitrary access are deposited
The various media that can store program code such as reservoir (RAM, Random Access Memory), magnetic or disk.It needs
Illustrate, herein, relational terms such as first and second and the like be used merely to by an entity or operation with
Another entity or operation distinguish, and without necessarily requiring or implying between these entities or operation, there are any this realities
The relationship or sequence on border.Moreover, the terms "include", "comprise" or its any other variant are intended to the packet of nonexcludability
Contain, so that the process, method, article or equipment for including a series of elements not only includes those elements, but also including
Other elements that are not explicitly listed, or further include for elements inherent to such a process, method, article, or device.
In the absence of more restrictions, the element limited by sentence "including a ...", it is not excluded that including the element
Process, method, article or equipment in there is also other identical elements.
The foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, for the skill of this field
For art personnel, the invention may be variously modified and varied.All within the spirits and principles of the present invention, made any to repair
Change, equivalent replacement, improvement etc., should all be included in the protection scope of the present invention.
The above description is merely a specific embodiment, but scope of protection of the present invention is not limited thereto, any
Those familiar with the art in the technical scope disclosed by the present invention, can easily think of the change or the replacement, and should all contain
Lid is within protection scope of the present invention.Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (9)
1. a kind of mud-rock-flow-impact load function generation method, which is characterized in that the described method includes:
Based on the actual information of the target mud-rock flow case got, multiple associated eigenvalues of mud-rock-flow-impact load are determined;
The multiple associated eigenvalue includes that the maximum value of mud-rock-flow-impact load amplitude, maximum value arrival time, amplitude are turned left a little
In a period of time for being greater than the time of turning right on time, the time of turning right of amplitude, impact load residual value and impact load tail portion
Carve corresponding amplitude;
Based on the multiple associated eigenvalue, multiple undetermined coefficients in mud-rock-flow-impact Load Model are determined;
Obtain the ripple frequency of mud-rock-flow-impact load;
Obtain the duration of mud-rock-flow-impact load;
Obtain the overburden depth range of mud-rock-flow-impact load;
By the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the covering
Depth bounds substitute into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load function;
Wherein, it is based on the multiple associated eigenvalue, determines that multiple undetermined coefficients in mud-rock-flow-impact Load Model include:
Turning left the time for maximum value, maximum value arrival time and amplitude based on the mud-rock-flow-impact load amplitude, passes through
Solution methods of nonlinear equations solves the mud-rock-flow-impact Load Model, determines the mud-rock-flow-impact load according to solving result
A part of coefficient in multiple undetermined coefficients of model;
The time of turning right of maximum value, maximum value arrival time, amplitude based on the mud-rock-flow-impact load amplitude, impact lotus
It carries on residual value and impact load tail portion in a moment corresponding amplitude for being greater than the time of turning right, is calculated by swarm intelligence
Method solves the mud-rock-flow-impact Load Model, and another part system in the multiple undetermined coefficient is determined according to solving result
Number.
2. the method according to claim 1, wherein being determined in the multiple undetermined coefficient according to solving result
A part of coefficient, comprising:
It is based onDetermine a part of coefficient in the multiple undetermined coefficient, AmaxFor the mud-rock-flow-impact
The maximum value of load amplitude, tAmaxFor maximum value arrival time, tLFor turning left the time for amplitude, α1、β1、γ1It is described one
Divide coefficient.
3. the method according to claim 1, wherein being determined in the multiple undetermined coefficient according to solving result
Another part coefficient, comprising:
It is based onDetermine another part coefficient in the multiple undetermined coefficient,
AmaxFor the maximum value of the mud-rock-flow-impact load amplitude, tAmaxFor maximum value arrival time, trFor turning right the time for amplitude,
R is impact load residual value, t1For on impact load tail portion in the moment for being greater than the time of turning right, A1For impact load tail
In the moment t for being greater than the time of turning right in portion1Corresponding amplitude, α2, β2, γ2, k, b are another part system
Number.
4. according to the method described in claim 3, it is characterized in that, the mud-rock-flow-impact Load Model is P (t, h)=k (t)
ρ gh, and By the multiple associated eigenvalue, described
Multiple undetermined coefficients, the ripple frequency, the duration and the overburden depth range substitute into the mud-rock-flow-impact lotus
It carries in model, generates mud-rock-flow-impact load function, comprising:
By the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, the duration and the covering
Depth bounds substitute into P (t, h)=k (t) ρ gh, and Generate mud-rock-flow-impact load function, P
(t, h) indicates that the mud-rock-flow-impact load pressure functional value that t and depth location h changes at any time, h are mud-rock flow under current location
Overburden depth, ρ is the mudstone current density obtained in advance, and g is acceleration of gravity, and k (t) is acceleration of gravity additional coefficient, w
It (t) is the continually changing sinusoidal waveform function of an amplitude, α1, α2, β1, β2, γ1, γ2, k, b are the multiple undetermined coefficient, and R is
Impact load residual value, e are the bottom of natural logrithm, and π is pi, and f is ripple frequency, and overburden depth range is the maximum value of 0 and h
Between, at the time of the duration obviously restrains corresponding for 0 to k (t) function tail portion.
5. the method according to claim 1, wherein the Swarm Intelligence Algorithm include population, it is ant colony, artificial
The shoal of fish or bee colony.
6. a kind of mud-rock-flow-impact load function generating means, which is characterized in that described device includes:
First determination unit determines mud-rock-flow-impact load for the actual information based on the target mud-rock flow case got
Multiple associated eigenvalues;The multiple associated eigenvalue includes the maximum value of mud-rock-flow-impact load amplitude, maximum value arrival
Greater than described on time, the time of turning left of amplitude, the time of turning right of amplitude, impact load residual value and impact load tail portion
It turns right a moment corresponding amplitude of time;
Second determination unit, for be based on the multiple associated eigenvalue, determine in mud-rock-flow-impact Load Model it is multiple to
Determine coefficient;
First acquisition unit, for obtaining the ripple frequency of mud-rock-flow-impact load;
Second acquisition unit, for obtaining the duration of mud-rock-flow-impact load;
Third acquiring unit, for obtaining the overburden depth range of mud-rock-flow-impact load;
Generation unit, for by the multiple associated eigenvalue, the multiple undetermined coefficient, the ripple frequency, described lasting
Time and the overburden depth range substitute into the mud-rock-flow-impact Load Model, generate mud-rock-flow-impact load function;
Wherein, it is based on the multiple associated eigenvalue, determines multiple undetermined coefficients in mud-rock-flow-impact Load Model, comprising:
Turning left the time for maximum value, maximum value arrival time and amplitude based on the mud-rock-flow-impact load amplitude, passes through
Solution methods of nonlinear equations solves the mud-rock-flow-impact Load Model, determines the mud-rock-flow-impact load according to solving result
A part of coefficient in multiple undetermined coefficients of model;
The time of turning right of maximum value, maximum value arrival time, amplitude based on the mud-rock-flow-impact load amplitude, impact lotus
It carries on residual value and impact load tail portion in a moment corresponding amplitude for being greater than the time of turning right, is calculated by swarm intelligence
Method solves the mud-rock-flow-impact Load Model, and another part system in the multiple undetermined coefficient is determined according to solving result
Number.
7. device according to claim 6, which is characterized in that second determination unit, for being based onDetermine a part of coefficient in the multiple undetermined coefficient, AmaxFor the mud-rock-flow-impact load amplitude
Maximum value, tAmaxFor maximum value arrival time, tLFor turning left the time for amplitude, α1、β1、γ1For a part of coefficient.
8. device according to claim 6, which is characterized in that second determination unit, for being based onDetermine another part coefficient in the multiple undetermined coefficient, AmaxFor institute
State the maximum value of mud-rock-flow-impact load amplitude, tAmaxFor maximum value arrival time, trFor turning right the time for amplitude, R is impact
Load residual value, t1For on impact load tail portion in the moment for being greater than the time of turning right, A1For on impact load tail portion big
In the moment t of the time of turning right1Corresponding amplitude, α2, β2, γ2, k, b are another part coefficient.
9. device according to claim 8, which is characterized in that the mud-rock-flow-impact Load Model is P (t, h)=k (t)
ρ gh, and By the multiple associated eigenvalue, described
Multiple undetermined coefficients, the ripple frequency, the duration and the overburden depth range substitute into the mud-rock-flow-impact lotus
It carries in model, generates mud-rock-flow-impact load function, P (t, h) indicates mud-rock-flow-impact load pressure t and depth location at any time
The functional value of h variation, h are the overburden depth of mud-rock flow under current location, and ρ is the mudstone current density obtained in advance, and g adds for gravity
Speed, k (t) are acceleration of gravity additional coefficient, and w (t) is the continually changing sinusoidal waveform function of an amplitude, α1, α2, β1, β2,
γ1, γ2, k, b are the multiple undetermined coefficient, and R is impact load residual value, and e is the bottom of natural logrithm, and π is pi, and f is arteries and veins
Dynamic frequency, overburden depth range between 0 and the maximum value of h, obviously restrain for 0 to k (t) function tail portion corresponding by the duration
Moment.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810527748.5A CN108846172B (en) | 2018-05-28 | 2018-05-28 | Mud-rock-flow-impact load function generation method and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810527748.5A CN108846172B (en) | 2018-05-28 | 2018-05-28 | Mud-rock-flow-impact load function generation method and device |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108846172A CN108846172A (en) | 2018-11-20 |
CN108846172B true CN108846172B (en) | 2019-08-30 |
Family
ID=64207980
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810527748.5A Expired - Fee Related CN108846172B (en) | 2018-05-28 | 2018-05-28 | Mud-rock-flow-impact load function generation method and device |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108846172B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111310623B (en) * | 2020-02-03 | 2023-04-07 | 中国地质大学(武汉) | Method for analyzing debris flow sensitivity map based on remote sensing data and machine learning |
CN112818574B (en) * | 2021-01-27 | 2022-10-14 | 江西理工大学 | Numerical method for simulating start-up formation, flow development and re-siltation of debris flow |
CN115575004A (en) * | 2022-08-18 | 2023-01-06 | 西南交通大学 | Impact performance calculation method and device for rock falling impacting buried pipeline and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026085A (en) * | 2010-07-20 | 2012-02-09 | Maruei Concrete Industries Co Ltd | Buffer block |
CN102819650A (en) * | 2012-08-16 | 2012-12-12 | 同济大学 | Computational simulation method of flow slide catastrophe of rock and soil material |
CN106652361A (en) * | 2017-03-16 | 2017-05-10 | 中国科学院、水利部成都山地灾害与环境研究所 | Mud-rock flow early warning method based on rainfall-probability |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000144692A (en) * | 1998-11-17 | 2000-05-26 | Nippon Kokan Light Steel Kk | Steel-made frame dam |
CN107766613A (en) * | 2017-09-09 | 2018-03-06 | 中铁二院工程集团有限责任公司 | The anti-mud-rock flow boulder impact design method of bridge pier |
CN108468305B (en) * | 2018-04-04 | 2019-12-10 | 中国科学院、水利部成都山地灾害与环境研究所 | Method for measuring and calculating impact force based on viscous debris flow acting on dam body |
CN108733942B (en) * | 2018-05-28 | 2020-03-10 | 中国地质大学(武汉) | Tunnel surrounding rock reliability assessment method and device |
-
2018
- 2018-05-28 CN CN201810527748.5A patent/CN108846172B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012026085A (en) * | 2010-07-20 | 2012-02-09 | Maruei Concrete Industries Co Ltd | Buffer block |
CN102819650A (en) * | 2012-08-16 | 2012-12-12 | 同济大学 | Computational simulation method of flow slide catastrophe of rock and soil material |
CN106652361A (en) * | 2017-03-16 | 2017-05-10 | 中国科学院、水利部成都山地灾害与环境研究所 | Mud-rock flow early warning method based on rainfall-probability |
Non-Patent Citations (3)
Title |
---|
水石流冲击信号能量分布试验研究;陈洪凯 等;《振动与冲击》;20121231;第56-59页 |
泥石流冲击荷载的视频分析方法及应用;陈洪凯 等;《防灾减灾工程学报》;20110630;第31卷(第3期);第255-260页 |
雅康路五里沟右线大桥抗泥石流冲击性能研究;江杨 等;《国防交通工程与技术》;20131231;第11-15,27页 |
Also Published As
Publication number | Publication date |
---|---|
CN108846172A (en) | 2018-11-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108846172B (en) | Mud-rock-flow-impact load function generation method and device | |
Qiu et al. | Generation of the North Equatorial Undercurrent jets by triad baroclinic Rossby wave interactions | |
EP2370957B1 (en) | System and method of grid generation for discrete fracture modeling | |
Lamb et al. | Two-dimensional numerical simulations of shoaling internal solitary waves at the ASIAEX site in the South China Sea | |
Adcock et al. | Did the Draupner wave occur in a crossing sea? | |
Grimshaw et al. | Combined effect of rotation and topography on shoaling oceanic internal solitary waves | |
CN108733942A (en) | The appraisal procedure and device of tunnel surrounding reliability | |
Lemieux et al. | Improving the simulation of landfast ice by combining tensile strength and a parameterization for grounded ridges | |
Gao et al. | Influence of offshore topography on the amplification of infragravity oscillations within a harbor | |
CN102914793A (en) | Nonequivalent modeling method and device for heterogeneous reservoir stratum | |
CN109856674A (en) | Engineering dessert evaluating method and device | |
CN110309518A (en) | Fractured-cavernous carbonate reservoir corrosion hole classification modeling method | |
Montoya et al. | Tsunami versus infragravity surge: comparison of the physical character of extreme runup | |
Mörner | Our Oceans-Our Future: New evidence-based sea level records from the Fiji Islands for the last 500 years indicating rotational eustasy and absence of a present rise in sea level | |
CN104391325B (en) | Discontinuous heterogeneous geologic body detection method and device | |
Huang et al. | Initial experiments on improving seismic data inversion with deep learning | |
Wuxi et al. | Effects of tide-surge interaction and wave set-up/set-down on surge: case studies of tropical cyclones landing China's Zhe-Min coast | |
Kulkarni | Numerical Modelling of Coastal Erosion using MIKE21 | |
CN108804785B (en) | A kind of impact load calculation method and device | |
Waseda | Nonlinear processes | |
CN106443773B (en) | A kind of method and device for suppressing coal seam screen effect in seismic profile data | |
Xiao et al. | Modifications to WRF's dynamical core to improve the treatment of moisture for large‐eddy simulations | |
CN106338768A (en) | Processing method, device and system for generating reservoir prediction attribute data | |
Cervantes et al. | Numerical model simulations of continental shelf flows off northern California | |
WO2022153984A1 (en) | Learning data generation method, model generation method, and learning data generation device |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190830 |