CN107014704B - A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave - Google Patents
A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave Download PDFInfo
- Publication number
- CN107014704B CN107014704B CN201710337797.8A CN201710337797A CN107014704B CN 107014704 B CN107014704 B CN 107014704B CN 201710337797 A CN201710337797 A CN 201710337797A CN 107014704 B CN107014704 B CN 107014704B
- Authority
- CN
- China
- Prior art keywords
- wave
- viscoelasticity
- strain
- stress
- viscosity coefficient
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/30—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight
- G01N3/303—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight generated only by free-falling weight
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0075—Strain-stress relations or elastic constants
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The present invention relates to a kind of short rock bar viscosity coefficient test methods that analysis is propagated based on viscoelasticity wave, comprising the following steps: establishes the governing equation of rock bar viscoelasticity constitutive equation and the propagation of one-dimensional viscoelasticity wave；Governing equation, transflector theory and the viscoelasticity defect theory propagated based on one-dimensional viscoelasticity wave design one-dimensional viscoelasticity stress wave analysis program；Impact test is simulated using one-dimensional viscoelasticity stress wave analysis program, input load wave is simultaneously compared the strain wave that impact test obtains is simulated with verifying wave, keep the strain waveform of simulation consistent with verifying waveform by adjusting parameter, and then obtains the viscosity coefficient of rock bar.The present invention is more reasonable by the viscosity coefficient that the method for the continuous trial and error of simulation test obtains；Only a complete waveform need to be measured in dynamic impact test as Loaded contact analysis, another verifying waveform should not Seeking Truth be no superposition and completely, can shorten the length of rock bar in testing, be conducive to the test of natural rock viscosity coefficient.
Description
Technical field
The present invention relates to a kind of rock viscosity coefficient test method more particularly to it is a kind of based on viscoelasticity wave propagate analysis
Short rock bar viscosity coefficient test method.
Background technique
In the rock mass engineering projects such as mining, tunnel excavation and the construction of hydropower stations, driling explosion generate mechanical oscillation and
The Explosive stress wave generated during explosive charge, can all cause the disturbance of surrounding rocks and be propagated outward in the form of wave.It encloses
There is various structural planes and joint etc. in rock, external disturbance is propagated in jointed rock mass can generate energy attenuation and waveform more
It dissipates, some rocks for containing microdefect and moisture property that not only demonstrates flexibility also shows certain adhesion properties, utilizes
Stress wave propagation in jointed rock mass can be equivalent to stress wave and propagated in viscous-elastic material by the method for effect medium, therefore be carried out
The research of viscoelasticity wave propagation law is the mechanism for inquiring into stress wave and propagating in jointed rock mass and then inducing mine power disaster
One ring of key.
It mainly using established viscoelasticity constitutive relation and is moved currently based on the viscosity coefficient test of EFFECTIVE MEDIUM method
The strain waveform of two measuring points (P1 and P2) on the rock bar measured in state impact test calculates stress wave on rock bar according to formula (1)
Then the attenuation coefficient of propagation is solved according to the relationship of viscosity coefficient in attenuation coefficient and constitutive equation.If viscoelasticity sheet
Structure relationship takes the constitutive equation of Kelvin-Voigt model, shown in the relationship such as formula (2) of viscosity coefficient and attenuation coefficient.
It is corresponding that each harmonic wave can be acquired according to the strain wave progress Fourier transformation measured of two measuring point P1 and P2 on rock bar
Attenuation coefficient, as shown in formula (1).
In formula, α is attenuation coefficient, and ω is the harmonic wave angular frequency after Fourier transformation, ε2And ε1Respectively P2 point and P1 point are surveyed
The strain obtained, F indicate that, to strain wave progress Fourier transformation, l is the distance between P1 point and P2 point, and Re indicates to take plural reality
Portion.
Using Kelvin-Voigt model foundation viscoelasticity constitutive equation, constructed then in conjunction with one-dimensional theory of wave propagation one-dimensional
Viscoelasticity wave transmission controe equation, viscosity coefficient can be acquired by equation (2).
The test method of above-mentioned viscosity coefficient exists following insufficient: firstly, being frequency by the attenuation coefficient that equation (1) acquires
Function, therefore it is not a definite value that the viscosity coefficient being finally calculated, which is also the function of frequency, needs to determine according to dominant frequency
The viscosity coefficient of rock bar；Secondly, test rock bar on two o'clock strain waveform, need to guarantee the length of rock bar can measure two it is complete
Whole Loaded contact analysis, it is higher to the length requirement of rock bar, and the natural rock pole length drilled through is shorter, therefore the rock met the requirements
Bar comparison is difficult.
Summary of the invention
The technical problem to be solved by the embodiment of the invention is that proposing the examination that laboratory test and numerical simulation combine
Proved recipe method constantly adjusts viscosity coefficient, one most suitable viscosity coefficient of optimization by the method for trial and error.
In order to solve the above technical problems, providing a kind of short rock bar viscosity coefficient test for propagating analysis based on viscoelasticity wave
Method, comprising the following steps:
The governing equation that the viscoelasticity constitutive equation and one-dimensional viscoelasticity wave for establishing rock bar are propagated；
It is one-dimensional using transflector theory and the design of viscoelasticity defect theory based on the governing equation that one-dimensional viscoelasticity wave is propagated
Viscoelasticity stress wave analysis program；
Impact test is simulated using one-dimensional viscoelasticity stress wave analysis program, input load wave simultaneously obtains simulation impact test
The strain wave obtained is compared with verifying wave, keeps the strain waveform of simulation consistent with verifying waveform by adjusting the parameter of setting,
And then obtain the viscosity coefficient of rock bar.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, the rock bar of establishing
Viscoelasticity constitutive equation specifically:
Using pendulum load in strain rate Hopkinson pressure bar test system carry out rock bar SHPB test, test its dynamic
Load-deformation curve；
Constitutive model is constructed by elastic element and stickiness element, establishes rock in conjunction with dynamic stress-strain curve of rock bar
The dynamic constitutive equation of bar.
It is described to establish one-dimensional stick in the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave
The governing equation of elastic wave propagation specifically:
The dynamic constitutive equation of the continuity equation, the equation of motion and rock bar propagated based on one-dimensional wave, establishes one-dimensional glutinous bullet
Property wave propagate governing equation.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, one-dimensional viscoelasticity is utilized
Stress wave analysis program simulate impact test process the following steps are included:
Pre-treatment step: establishing the geometrical model of rock bar sample, carries out grid dividing for model partition into multiple units, right
The relevant parameter of each unit carries out assignment, and boundary condition and loading environment is arranged；
Stress analysis step: carry out node Stress calculation simultaneously judging unit whether occur damage and to degree of impairment into
Row processing；
Post-processing step: the stress value, strain value and elastic mould value of key node are exported.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, walked in the preceding processing
Carry out assignment to following parameter in rapid: density, elasticity modulus, cross-sectional area, viscosity coefficient, uniaxial compressive strength and single shaft are anti-
Tensile strength；
Loading environment is set in pre-treatment step specifically: to loading form, the load wave being applied on rock bar sample
Shape and waveform parameter are configured.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, stress analysis step tool
Body the following steps are included:
One new boundary stress is applied to rock bar sample；
According to the stress value and strain value of governing equation and transflector theoretical calculation node；
Unit damage is judged whether there is according to failure criteria, if haveing damage, is calculated and is damaged according to viscoelasticity defect theory
Situation, and modify the stress value and strain value that each node is recalculated after relevant physical parameter；
If not damaging, the stress value and strain value of each node are exported.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, stickiness damage is according to such as
Lower formula is calculated:
η=(1-D) η0；
Wherein,
D is damage factor in formula, between 0 to 1, η and η0It is the viscosity coefficient after damaging and before damage respectively；
ε is the strain of unit, εt0For elastic limit strain, εtuFor the limit tensile strain of unit, εc0For elastic limit pair
The compressive strain answered；
E0Elasticity modulus before damage, ftrFor tensile failure residual strength, fcrFor compression failure residual strength.
In the short rock bar viscosity coefficient test method for propagating analysis the present invention is based on viscoelasticity wave, post-processing step is specific
Are as follows:
Stress value, strain value, the elastic mould value of key node are extracted from the calculated result of stress analysis step, and defeated
Stress diagram, strain figure and springform spirogram out.
In natural short rock bar viscosity coefficient test method of the invention, according to the stress value, strain value, elasticity of key node
Modulus value draws strain waveform figure, and strain wave and verifying wave are compared, and makes strain waveform and verifying wave by adjusting parameter
Shape is consistent, and then obtains the viscosity coefficient of rock bar.
In short rock bar viscosity coefficient test method of the invention natural, using pendulum load in strain rate Hopkinson pressure bar
Pilot system carries out the impact test of natural short rock bar, obtains the waveform of two measuring points respectively as required in simulation impact test
Loaded contact analysis and verifying waveform.
Above-mentioned technical proposal at least has the following beneficial effects: that the invention proposes laboratory tests and numerical simulation to combine
Method, measure the dynamic of rock this structure curve by SHPB test and construct corresponding numerical model, then using dynamically
Impact test measures the Loaded contact analysis needed in numerical simulation and verifying waveform, constantly adjusts stickiness system by the method for trial and error
Number, one most suitable viscosity coefficient of optimization.
For the deficiency of existing test method, test method of the invention has following improvement: firstly, measured by test
Dynamic stress-strain curve building dynamic constitutive is more in line with the dynamic characteristic of rock, measures the viscosity coefficient of rock more
Adduction reason；Secondly, the viscosity coefficient obtained by the method for the continuous trial and error of numerical simulation corresponds to viscosity coefficient more compared to dominant frequency
It is representative；Again, only need to measure a complete waveform in dynamic impact test of the invention as numerical model
Loaded contact analysis, another verifying waveform are not required for whether having superposition and length that is complete, therefore can shortening rock bar in test,
Be conducive to the test of natural rock viscosity coefficient.
Detailed description of the invention
Fig. 1 is Kelvin-Voigt model schematic；
Elastic damage constitutive relation schematic diagram when Fig. 2 microscopic element single shaft tension；
Fig. 3 is one-dimensional viscoelasticity stress wave analysis program flow diagram.
Specific embodiment
It is including following the present invention provides a kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave
Step: the governing equation that the viscoelasticity constitutive equation and one-dimensional viscoelasticity wave for establishing rock bar are propagated；It is passed based on one-dimensional viscoelasticity wave
Governing equation, transflector theory and the viscoelasticity defect theory broadcast design one-dimensional viscoelasticity stress wave analysis program；Utilize one
Tie up viscoelasticity stress wave analysis program simulate impact test, input load wave and by simulate impact test obtain strain wave with test
Card wave compares, and keeps the strain waveform of simulation consistent with verifying waveform by adjusting the parameter of setting, and then obtain rock bar
Viscosity coefficient.
Wherein, the viscoelasticity constitutive equation of rock bar is established specifically: utilize strain rate Hopkinson pressure bar in pendulum load
Pilot system carries out the SHPB test of rock bar, tests dynamic stress-strain curve of rock bar sample.By elastic element and stick
Property element construct constitutive model, elastic element can be used spring, and stickiness element is using glutinous pot.Fig. 1 institute can be used when specific implementation
Constitutive model of the Kevlin body component models shown as viscous-elastic material is answered in conjunction with the dynamic of the rock bar of SHPB test acquisition
Force-strain curve establishes the dynamic constitutive equation of rock bar, as shown in formula (3).
In formula, σ, ε are the stress and strain of Kevlin body respectively, and E and η are the elasticity modulus of spring in Kevlin body respectively
With the viscosity coefficient of newton body.
After establishing the dynamic constitutive equation of rock, continuity equation, the equation of motion and rock based on one-dimensional wave propagation
Dynamic constitutive equation establishes the governing equation that one-dimensional viscoelasticity wave is propagated.
From the continuity equation of one-dimensional longitudinal wave:
The equation of motion of one-dimensional longitudinal wave can be obtained by Newton's Theorem:
In formula, v is particle velocity, and ρ is the density of medium, and x is the coordinate of particle.
It can be in the hope of the governing equation of one-dimensional viscoelasticity wave according to formula (3) (4) (5):
The harmonic wave solution of equation (6) are as follows:
U (x, t)=u0 exp(αx)exp[i(ωt-kx)] (7)
In formula, α and k are respectively decay factor and wave number, and α is negative number representation harmonic wave with propagation distance exponential damping, u and ω
Respectively indicate displacement and angular frequency.
Phase velocity c=ω/k of stress wave propagation, is the function of frequency, can be obtained by equation (8)
After constructing the governing equation that one-dimensional viscoelasticity wave is propagated based on the above process, damaged using transflector theory and viscoelasticity
Hurt Theoretical Design one-dimensional wave and propagates analysis program.The design principle that one-dimensional wave propagates analysis program is specifically described below:
1, stress wave is in viscoelastic rod interface transflector
When one-dimensional harmonic wave is propagated in no interface viscoelastic rod, harmonic amplitude only declines with propagation distance according to exponential function
Subtract, harmonic frequency and phase velocity do not change.By formula (10) it is found that the spread speed of harmonic wave is not only by the material properties shadow of bar
It rings and related with harmonic frequency.When one-dimension stress wane propagates the interface for encountering material properties variation in elastic rod, stress
Wave increases and decreases through the waveforms amplitude at interface according to transflector relationship, but frequency remains unchanged, i.e. stress wave waveform figure is longitudinally stretched
Contracting and it is laterally invariant.When harmonic wave is propagated in viscoelastic rod, the material properties of viscoelastic rod have no effect on the frequency of harmonic wave, therefore
When encountering the interface of material properties variation when harmonic wave is propagated in viscoelastic rod, harmonic wave passes through the waveform frequency after interface
Rate also remains unchanged, i.e. harmonic wave forms figure is also longitudinal extension and laterally invariant.
From the condition of continuity:
FI+FR=FT (11)
VI+VR=VT (12)
It is obtained by the conservation of momentum:
FIDt=ρ1A1c1dtVI FI=σIA1 (13)
FRDt=- ρ1A1c1dtVR FR=σRA1 (14)
FTDt=ρ2A2c2dtVT FT=σTA2 (15)
In formula: F, V and σ are respectively power, particle velocity and stress on interface, and subscript I, R and T respectively indicate incidence
Wave, back wave and transmitted wave, ρ1And ρ2The respectively density on the interface left side and the right, c1 and c2 be respectively the interface left side and
The harmonic velocity on the right, dt are unit time, A1And A2The respectively sectional area on the interface left side and the right.
It is obtained by formula (11)-(15):
By formula (16) and (17) it is found that when harmonic wave is by the interface of material and cross-section variation, the transmission coefficient of harmonic wave
It is related with the density of material on interface both sides, harmonic wave phase velocity and cross-sectional area with reflection coefficient.Since harmonic wave is in viscoelastic rod
The phase velocity of middle propagation is influenced by material properties and harmonic frequency, therefore the transmission coefficient of harmonic wave and reflection coefficient are by harmonic wave
The harmonic wave transflector coefficient of frequency influence, different frequency is different.Material is passed through by the stress wave that the harmonic superposition of different frequency is formed
When the interface of matter and cross-section variation, since the transflector coefficient of different frequency harmonic wave is inconsistent, cause to be formed after superposition
The frequency of stress wave will change, i.e., behind the interface that stress wave is changed by wave impedance, stress wave waveform is not only in longitudinal direction
Flexible and flexible laterally having, the transflector rule that the stress wave of this and arbitrary shape is propagated in one-dimensional elastic rod is not yet
Together.
2, viscoelasticity wave loads BORDER PROCESSING
Since strain and stress equation form is the same in one-dimensional viscoelasticity wave governing equation, stress or strain harmonic wave exist
Propagation law is consistent in viscoelasticity stock.But it due to loading the inconsistent of boundary, needs to perform corresponding processing.
(1) variable boundary is answered
The load boundary of model is strain controlling, needs to solve stress by strain, by this structure side of equation (3) Kevlin body
Journey based on the particle stress that strain wave is propagated it is found that solved as follows:
ε in formulat、σt、EtAnd ηtIndicate that strain, stress, elasticity modulus and the viscosity coefficient of t moment particle, Δ t indicate
Incremental time.The stress at particle moment can be solved by the strain at this moment and last moment.
(2) stress boundary
The load boundary of model is stress loading, and model medium wave is propagated in the form strained, and stress wave is needed to turn
Strain wave is changed into, from equation (20)
The initial time of stress-wave loading, t=0, ε=0.
3, stickiness defect theory
In one-dimensional viscoelasticity stock viscous-elastic material unit institute's tension stress or compression be greater than uniaxial tensile strength or
Material damages when compression strength, and damage includes two parts, and a part is elastic damage, by the Young's modulus for changing material
Express the influence of damage, another part is stickiness damage, express the influence damaged to material by changing viscosity coefficient.
(1) elastic damage principle
According to Elastic Damage Theory, after unit generates damage, material parameter will be weakened, and elasticity modulus is according to formula (22)
It is handled.
E=(1-D) E0 (22)
D is damage factor in formula, between 0 to 1；E and E0It is the elasticity modulus after damaging and before damage respectively.
It is assumed that the elastic damage constitutive relation of microscopic element is as shown in Figure 2.The tensile stress suffered by the unit is greater than material
When tensile strength, there is tensile damage in unit, and the strain of unit at this time is ε, it is known that εt0≤ε<εtu, εt0It is corresponding for elastic limit
Strain, that is, elastic limit strain, εtuFor the limit tensile strain of unit.Have centainly since the unit not yet destroys still completely
Bearing capacity, and this bearing capacity can define a residual strength ftrIt characterizes, the calculating of damage factor such as formula (23)
It is shown.
When the compression suffered by the unit is greater than the compression strength of material, there is crushing wound, use and tensile damage in unit
Similar method handles compressive damage, as shown in formula (24).
In formula, εc0For the corresponding compressive strain of elastic limit, fcrFor compression failure residual strength.
(2) stickiness damage principle
The strain when calculating of the elastic damage factor is elasticity modulus and the damage before utilizing residual strength, material damage.
When material is viscous-elastic material, the stress of unit is made of Elastic Term and stickiness item two parts, is damaged by means of elastic material
The calculation method of the factor calculates the damage factor of viscous-elastic material.Tensile damage and the damage factor of compressive damage difference are as follows
Shown in formula (25) and (26).
The viscosity coefficient and damage of intact rock in stickiness damage are provided referring to the calculation formula of Young's modulus in elastic damage
The equation (27) of rock viscosity coefficient after wound.
η=(1-D) η0 (27)
η and η in formula0It is the viscosity coefficient after damaging and before damage respectively.
According to the above-mentioned communication theory about stress wave in one-dimensional viscoelastic rod, transflector is theoretical and viscoelasticity damage is managed
The one-dimensional viscoelasticity stress wave analysis program of the design of opinion, program circuit are as shown in Figure 3.One-dimensional viscoelasticity stress wave analysis is used again
Program simulates the strain wave of impact test test, and strain waveform and verifying waveform are compared, and is calculated using the method for trial and error
The viscosity coefficient of natural rock.Specific implementation process is as follows:
As shown in figure 3, using one-dimensional wave propagate analysis program simulation impact test process include pre-treatment, stress analysis and
Post-process three parts.Pre-treatment step major function is to establish the geometrical model of rock bar sample, carries out grid dividing and draws model
It is divided into multiple units, assignment is carried out to the relevant parameter of each unit, boundary condition and loading environment are set.Stress analysis step master
Wanting function is to carry out node stress to calculate whether judging unit simultaneously occurs to damage and handle degree of impairment.Post-processing step
Rapid major function is stress value, strain value and the elastic mould value for exporting key node.
In pre-treatment step, carry out assignment to the following parameter of each unit first: density, cross-sectional area, is sticked elasticity modulus
Property coefficient, uniaxial compressive strength and uniaxial tensile strength.Again to be applied on rock bar sample loading form, Loaded contact analysis with
And waveform parameter is configured.When it is implemented, first with pendulum load in strain rate Hopkinson pressure bar test system carry out
The impact test of natural short rock bar, and then obtain the waveform of two measuring points.It is impacted one of them complete waveform as simulation
Loaded contact analysis and another waveform in test is as verifying waveform.
After the completion of pre-treatment step, stress analysis step is carried out.Stress analysis step specifically includes the following steps:
One new boundary stress is applied to rock bar sample；
According to the stress value and strain value of governing equation and each node of transflector theoretical calculation；
Unit damage is judged whether there is according to failure criteria, if haveing damage according to stickiness defect theory and Elastic Damage Theory
Degree of impairment is calculated, and modifies the stress value and strain value for recalculating each node after relevant physical parameter；
If not damaging, the stress value and strain value of each node are exported.
Stress value, strain value, the elastic mould value of key node are extracted from the calculated result of stress analysis step, and defeated
Stress diagram, strain figure and springform spirogram out.
Strain waveform figure is drawn according to the stress value, strain value, elastic mould value of key node, by strain wave and verifying wave
It compares, keeps strain waveform consistent with verifying waveform by adjusting parameter, and then obtain the viscosity coefficient of rock bar.
The test method that the present invention uses laboratory test and numerical simulation to combine, initially sets up this structure of the viscoelasticity of rock bar
The governing equation that equation and one-dimensional viscoelasticity wave are propagated；Based on one-dimensional viscoelasticity wave propagate governing equation, transflector theory with
And viscoelasticity defect theory designs one-dimensional viscoelasticity stress wave analysis program；It is simulated using one-dimensional viscoelasticity stress wave analysis program
Impact test, input load wave is simultaneously compared the strain wave that impact test obtains is simulated with verifying wave, by adjusting setting
Parameter make simulation strain waveform with verifying waveform it is consistent, and then acquisition rock bar viscosity coefficient.
The dynamic stress that the present invention is measured by test-strain curve constructs dynamic constitutive and is more in line with the dynamic of rock
Step response, the viscosity coefficient for measuring rock are more reasonable；The viscosity coefficient phase obtained by the method for the continuous trial and error of numerical simulation
It is more representative that viscosity coefficient is corresponded to than dominant frequency；It only needs to measure a complete wave in dynamic impact test of the invention
Loaded contact analysis of the shape as numerical model, another verifying waveform is not required for whether having superposition and complete, therefore can shorten
The length of rock bar in test, is conducive to the test of natural rock viscosity coefficient.
The foregoing is a specific embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (8)
1. a kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave, which is characterized in that including following step
It is rapid:
Using pendulum load in strain rate Hopkinson pressure bar test system carry out rock bar SHPB test, test its dynamic and answer
Force-strain curve；
It is built by elastic element and stickiness element building Kevlin body component models in conjunction with dynamic stress-strain curve of rock bar
The dynamic constitutive equation of vertical rock bar；
The dynamic constitutive equation of the continuity equation, the equation of motion and rock bar propagated based on one-dimensional wave, establishes one-dimensional viscoelasticity wave
The governing equation of propagation；
Based on the governing equation that one-dimensional viscoelasticity wave is propagated, one-dimensional glutinous bullet is designed using transflector theory and viscoelasticity defect theory
Property stress wave analysis program；
Impact test is simulated using one-dimensional viscoelasticity stress wave analysis program, input load wave simultaneously will simulate what impact test obtained
Strain wave is compared with verifying wave, keeps the strain waveform of simulation consistent with verifying waveform by adjusting the parameter of setting, in turn
Obtain the viscosity coefficient of rock bar.
2. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as described in claim 1, feature exists
In, using one-dimensional viscoelasticity stress wave analysis program simulation impact test process the following steps are included:
Pre-treatment step: establishing the geometrical model of rock bar sample, carries out grid dividing for model partition at multiple units, to each list
The relevant parameter of member carries out assignment, and boundary condition and loading environment is arranged；
Stress analysis step: carry out node Stress calculation simultaneously judging unit whether occur damage and to degree of impairment at
Reason；
Post-processing step: the stress value, strain value and elastic mould value of key node are exported.
3. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as claimed in claim 2, feature exists
In to following parameter progress assignment in the pre-treatment step: density, elasticity modulus, cross-sectional area, viscosity coefficient, single shaft
Compression strength and uniaxial tensile strength；
Loading environment is set in pre-treatment step specifically: to be applied on rock bar sample loading form, Loaded contact analysis with
And waveform parameter is configured.
4. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as claimed in claim 2, feature exists
In, stress analysis step specifically includes the following steps:
One new boundary stress is applied to rock bar sample；
According to the stress value and strain value of governing equation and transflector theoretical calculation node；
Unit damage is judged whether there is according to failure criteria, if haveing damage, according to stickiness defect theory and Elastic Damage Theory
Degree of impairment is calculated, and modifies the stress value and strain value for recalculating each node after relevant physical parameter；
If not damaging, the stress value and strain value of each node are exported.
5. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as claimed in claim 4, feature exists
In stickiness damage is calculated according to the following formula:
η=(1-D) η0；
Wherein,
D is damage factor in formula, between 0 to 1, η and η0It is the viscosity coefficient after damaging and before damage respectively；
ε is the strain of unit, εt0For elastic limit strain, εtuFor the limit tensile strain of unit, εc0It is corresponding for elastic limit
Compressive strain；
E0Elasticity modulus before damage, ftrFor tensile failure residual strength, fcrFor compression failure residual strength.
6. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as claimed in claim 2, feature exists
In post-processing step specifically:
Stress value, strain value, the elastic mould value of key node are extracted from the calculated result of stress analysis step, and are exported and answered
Try hard to, strain figure and springform spirogram.
7. the short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave as described in claim 2 or 6, feature
Be, according to the stress value, strain value, elastic mould value of key node draw strain waveform figure, by strain wave and verifying wave into
Row comparison keeps strain waveform consistent with verifying waveform, and then obtains the viscosity coefficient of rock bar by adjusting parameter.
8. propagating the short rock bar viscosity coefficient test method of analysis based on viscoelasticity wave as described in claim 1, feature exists
In, using pendulum load in strain rate Hopkinson pressure bar test system carry out the impact test of natural short rock bar, obtain two
The waveform of measuring point is respectively as Loaded contact analysis needed for simulation impact test and verifying waveform.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710337797.8A CN107014704B (en) | 2017-05-15 | 2017-05-15 | A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710337797.8A CN107014704B (en) | 2017-05-15 | 2017-05-15 | A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107014704A CN107014704A (en) | 2017-08-04 |
| CN107014704B true CN107014704B (en) | 2019-06-25 |
Family
ID=59450105
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710337797.8A Active CN107014704B (en) | 2017-05-15 | 2017-05-15 | A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107014704B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108195725B (en) * | 2017-12-13 | 2020-06-26 | 中国人民解放军陆军工程大学 | Method for measuring one-dimensional stress wave attenuation effect of particulate matter based on SHPB (split harmonic elimination photon quenching) |
| CN110133104B (en) * | 2019-05-30 | 2021-07-30 | 东北大学 | Method for testing dynamic mechanical characteristics of full-age filling body |
| CN114861395A (en) * | 2022-03-30 | 2022-08-05 | 西北核技术研究所 | Method and model for analyzing parameters and energy of elastic-plastic stress wave in one-dimensional rod |
| CN117871293B (en) * | 2024-03-11 | 2024-05-28 | 四川嘉乐地质勘察有限公司 | Self-correction-based high strain force detection method and detection device |
| CN117875096B (en) * | 2024-03-13 | 2024-05-24 | 西安交通大学 | Rock stress wave load propagation characteristic simulation method and related device |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002060347A2 (en) * | 2001-01-31 | 2002-08-08 | Maria-Grazia Ascenzi | Modeling viscoelastic torsional properties of osteons |
| CN102141545A (en) * | 2010-11-27 | 2011-08-03 | 江西理工大学 | Method for testing rock mass mechanics parameters based on explosion seismic wave space-time attenuation law |
| CN103234844A (en) * | 2013-03-29 | 2013-08-07 | 东北大学 | Pendulum load intermediate strain rate split Hopkinson pressure bar test device and method |
| CN103412327A (en) * | 2013-08-01 | 2013-11-27 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | Viscoelasticity parameter extraction method of fractured reservoir |
| CN103994921A (en) * | 2014-06-12 | 2014-08-20 | 江西理工大学 | Testing method for elasticity modulus of rock mass weak intercalated layer based on wavelet waveform change rule |
| CN105136347A (en) * | 2015-05-19 | 2015-12-09 | 中国石油大学(华东) | Visco-elastic strain field ground stress measurement device and measurement method |
| CN106353797A (en) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | High-precision earthquake forward modeling method |
-
2017
- 2017-05-15 CN CN201710337797.8A patent/CN107014704B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2002060347A2 (en) * | 2001-01-31 | 2002-08-08 | Maria-Grazia Ascenzi | Modeling viscoelastic torsional properties of osteons |
| CN102141545A (en) * | 2010-11-27 | 2011-08-03 | 江西理工大学 | Method for testing rock mass mechanics parameters based on explosion seismic wave space-time attenuation law |
| CN103234844A (en) * | 2013-03-29 | 2013-08-07 | 东北大学 | Pendulum load intermediate strain rate split Hopkinson pressure bar test device and method |
| CN103412327A (en) * | 2013-08-01 | 2013-11-27 | 中国石油化工股份有限公司胜利油田分公司地质科学研究院 | Viscoelasticity parameter extraction method of fractured reservoir |
| CN103994921A (en) * | 2014-06-12 | 2014-08-20 | 江西理工大学 | Testing method for elasticity modulus of rock mass weak intercalated layer based on wavelet waveform change rule |
| CN105136347A (en) * | 2015-05-19 | 2015-12-09 | 中国石油大学(华东) | Visco-elastic strain field ground stress measurement device and measurement method |
| CN106353797A (en) * | 2015-07-17 | 2017-01-25 | 中国石油化工股份有限公司 | High-precision earthquake forward modeling method |
Non-Patent Citations (3)
| Title |
|---|
| "Dynamic Brazilian Test of Rock Under Intermediate Strain Rate:Pendulum Hammer-Driven SHPB Test and Numerical Simulation";W.C.Zhu 等;《Rock Mech Rock Eng》;20150930;第48卷(第5期);全文 |
| "岩体中应力波幅值随时空衰减的关系";王观石 等;《岩土力学》;20101130;第31卷(第11期);全文 |
| "岩石试样SHPB实验过程中应力均匀性的一维波动分析";牛雷雷 等;《东北大学学报（自然科学版）》;20150531;第36卷;全文 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107014704A (en) | 2017-08-04 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107014704B (en) | A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave | |
| Zhu et al. | Peridynamic modelling of blasting induced rock fractures | |
| Abedini et al. | Comparison of ALE, LBE and pressure time history methods to evaluate extreme loading effects in RC column | |
| CN105224742B (en) | A kind of blast action of analyzing is got out of a predicament or an embarrassing situation the method for slope stability | |
| Zhao et al. | UDEC modelling on wave propagation across fractured rock masses | |
| Strayer et al. | Numerical modeling of fold initiation at thrust ramps | |
| Wang et al. | Coupled finite element and discrete element method for underground blast in faulted rock masses | |
| Zhou et al. | Peridynamic micro-elastoplastic constitutive model and its application in the failure analysis of rock masses | |
| Wu et al. | A reformulated rate-dependent visco-elastic model for dynamic deformation and fracture of PMMA with peridynamics | |
| CN105589999B (en) | Determination method applied to underground engineering wall rock grouting program | |
| Aliabadian et al. | Numerical modeling of presplitting controlled method in continuum rock masses | |
| Siamaki et al. | Degradation of a discrete infilled joint shear strength subjected to repeated blast-induced vibrations | |
| Dipasquale et al. | A Stress Tensor-based Failure Criterion for Ordinary State-based Peridynamic Models | |
| CN111062113A (en) | Novel stope mining blasting parameter comprehensive optimization method under complex filling body condition | |
| Li et al. | Crack initiation and propagation simulation for polycrystalline-based brittle rock utilizing three dimensional distinct element method | |
| SÜLÜ | Mechanical behavior of mode-I crack propagation between composite adhesive joints | |
| Lin et al. | A data assimilation method for blast load prediction | |
| Ahmed | Models for analysis of shotcrete on rock exposed to blasting | |
| Yang et al. | Dynamic mechanical response and fracture mechanism of rock-like materials with concealed structural planes of different filling thicknesses | |
| Azevedo et al. | Numerical simulation of soft-body impact on GFRP laminate composites: mixed SPH-FE and pure SPH approaches | |
| CN114218831A (en) | General blasting numerical simulation method | |
| Xinping et al. | The damage fracture analysis of a jointed rock mass and its application in engineering | |
| CN116680832B (en) | Tunnel surrounding rock stability evaluation method considering creep characteristics of structural surface | |
| Shterenlikht et al. | Towards mechanism-based simulation of impact damage using exascale computing | |
| Gong et al. | Modelling rock fracturing by a novel implicit continuous to discontinuous method |
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 |