CN107014704A - A kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple - Google Patents
A kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple Download PDFInfo
- Publication number
- CN107014704A CN107014704A CN201710337797.8A CN201710337797A CN107014704A CN 107014704 A CN107014704 A CN 107014704A CN 201710337797 A CN201710337797 A CN 201710337797A CN 107014704 A CN107014704 A CN 107014704A
- Authority
- CN
- China
- Prior art keywords
- viscoelasticity
- ripple
- strain
- rock bar
- 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.)
- Granted
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
Abstract
The present invention relates to a kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple, comprise the following steps:Set up the governing equation that rock bar viscoelasticity constitutive equation and one-dimensional viscoelasticity ripple are propagated;Governing equation, transflector theory and the viscoelasticity defect theory propagated based on one-dimensional viscoelasticity ripple design one-dimensional viscoelasticity stress wave analysis program;Impact test is simulated using one-dimensional viscoelasticity stress wave analysis program, input loading ripple is simultaneously contrasted the strain wave for simulating impact test acquisition with checking ripple, make the strain waveform of simulation consistent with checking waveform by adjusting parameter, and then obtain the viscosity coefficient of rock bar.The viscosity coefficient that the present invention is obtained by the method for the continuous trial and error of simulated test is more reasonable;A complete waveform only need to be measured as Loaded contact analysis in dynamic impact test, and another checking waveform should not Seeking Truth be no has superposition and complete, can shorten the length of rock bar in experiment, be conducive to the test of natural rock viscosity coefficient.
Description
Technical field
The present invention relates to a kind of rock viscosity coefficient method of testing, more particularly to a kind of analysis is propagated based on viscoelasticity ripple
Short rock bar viscosity coefficient method of testing.
Background technology
In the rock mass engineering projects such as mining, tunnel excavation and the construction of hydropower stations, driling explosion produce mechanical oscillation and
The Explosive stress wave produced during explosive charge, can all cause the disturbance of surrounding rocks and outwards be propagated in the form of ripple.Enclose
Various structural planes and joint etc. are there is in rock, external disturbance is propagated in jointed rock mass can produce energy attenuation and waveform more
Dissipate, 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 be propagated in viscous-elastic material by the method for effect medium, therefore be carried out
The research of viscoelasticity ripple propagation law is to inquire into the mechanism that stress wave propagated and then induced mine power disaster in jointed rock mass
A crucial ring.
The viscosity coefficient test for being currently based on EFFECTIVE MEDIUM method mainly uses the viscoelasticity constitutive relation set up and moved
The strain waveform of two measuring points (P1 and P2) on the rock bar measured in state impact test, stress wave is calculated on rock bar according to formula (1)
The attenuation coefficient of propagation, is then solved according to the relation of viscosity coefficient in attenuation coefficient and constitutive equation.If viscoelasticity sheet
Structure relation takes the constitutive equation of Kelvin-Voigt models, shown in the relation such as formula (2) of viscosity coefficient and attenuation coefficient.
It is corresponding that each harmonic wave can be tried to achieve according to the strain wave progress Fourier transformation that two measuring point P1 and P2 on rock bar are measured
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 points and P1 points are surveyed
The strain obtained, F represents that to strain wave progress Fourier transformation l is the distance between P1 points and P2 points, and Re represents to take the reality of plural number
Portion.
Viscoelasticity constitutive equation is set up using Kelvin-Voigt models, is built then in conjunction with one-dimensional theory of wave propagation one-dimensional
Viscoelasticity ripple transmission controe equation, viscosity coefficient can be tried to achieve by equation (2).
The method of testing of above-mentioned viscosity coefficient exists following not enough:First, the attenuation coefficient tried to achieve by equation (1) is frequency
Function, therefore finally calculate obtained viscosity coefficient be also frequency function not a definite value, it is necessary to according to dominant frequency determine
The viscosity coefficient of rock bar;Secondly, test rock bar on 2 points strain waveform, it is necessary to ensure rock bar length 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 meet desired rock
Bar comparison is difficult.
The content of the invention
Technical problem to be solved of the embodiment of the present invention is, it is proposed that the examination that laboratory test and numerical simulation are combined
Proved recipe method, viscosity coefficient, one most suitable viscosity coefficient of optimization are constantly adjusted by the method for trial and error.
In order to solve the above technical problems, being tested there is provided a kind of short rock bar viscosity coefficient for propagating analysis based on viscoelasticity ripple
Method, comprises the following steps:
Set up the governing equation that the viscoelasticity constitutive equation and one-dimensional viscoelasticity ripple of rock bar are propagated;
The governing equation propagated based on one-dimensional viscoelasticity ripple, it is theoretical one-dimensional with the design of viscoelasticity defect theory using transflector
Viscoelasticity stress wave analysis program;
Impact test is simulated using one-dimensional viscoelasticity stress wave analysis program, input loading ripple simultaneously obtains simulation impact test
The strain wave obtained is contrasted with checking ripple, makes the strain waveform of simulation consistent with checking waveform by adjusting the parameter set,
And then obtain the viscosity coefficient of rock bar.
In the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple, the rock bar of setting up
Viscoelasticity constitutive equation is specially:
Strain rate Hopkinson pressure bar test system carries out the SHPB experiments of rock bar in being loaded using pendulum, tests its dynamic
Load-deformation curve;
Constitutive model is built by flexible member and stickiness element, rock is set up with reference to dynamic stress-strain curve of rock bar
The dynamic constitutive equation of bar.
It is described to set up one-dimensional glutinous in the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple
The governing equation of elastic wave propagation is specially:
The dynamic constitutive equation of continuity equation, the equation of motion and the rock bar propagated based on one-dimensional wave, sets up one-dimensional glutinous bullet
Property ripple propagate governing equation.
In the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple, one-dimensional viscoelasticity is utilized
The process of stress wave analysis program simulation impact test comprises the following steps:
Pre-treatment step:The geometrical model of rock bar sample is set up, mesh generation is carried out by model partition into multiple units, it is right
The relevant parameter of each unit carries out assignment, sets boundary condition and loading environment;
Stress analysis step:The Stress calculation of node is carried out while whether judging unit occurs to damage and enter degree of impairment
Row processing;
Post-processing step:Export stress value, strain value and the elastic mould value of key node.
In the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple, in the preceding processing step
Assignment is carried out to following parameter in rapid:Density, modulus of elasticity, cross-sectional area, viscosity coefficient, uniaxial compressive strength and single shaft are anti-
Tensile strength;
Setting loading environment is specially in pre-treatment step:To loading form, the loading ripple being applied on rock bar sample
Shape and waveform parameter are configured.
In the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple, stress analysis step tool
Body comprises the following steps:
One new boundary stress is applied to rock bar sample;
According to governing equation and the stress value and strain value of transflector theoretical calculation node;
Determine whether that unit is damaged according to failure criteria, if there is damage, calculated and damaged according to viscoelasticity defect theory
Situation, and change 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 method of testing that the present invention propagates analysis based on viscoelasticity ripple, 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;
E0Modulus of elasticity before damage, ftrFor tensile failure residual strength, fcrFor compression failure residual strength.
In the short rock bar viscosity coefficient method of testing that the present invention propagates analysis based on viscoelasticity ripple, post-processing step is specific
For:
Stress value, strain value, the elastic mould value of key node are extracted from the result of calculation of stress analysis step, and it is defeated
Go out stress diagram, strain figure and springform spirogram.
In natural short rock bar viscosity coefficient method of testing of the invention, according to the stress value, strain value, elasticity of key node
Modulus value draws strain waveform figure, and strain wave is contrasted with checking ripple, makes strain waveform and checking ripple by adjusting parameter
Shape is consistent, and then obtains the viscosity coefficient of rock bar.
In natural short rock bar viscosity coefficient method of testing of the invention, strain rate Hopkinson pressure bar in being loaded using pendulum
Pilot system carries out the impact test of natural short rock bar, obtains the waveform of two measuring points respectively as needed for simulation impact test
Loaded contact analysis and checking waveform.
Above-mentioned technical proposal at least has the advantages that:The present invention proposes laboratory test and numerical simulation is combined
Method, dynamic this structure curve for measuring rock is tested by SHPB and then corresponding numerical model is built, dynamic is then utilized
Impact test measures the Loaded contact analysis needed in numerical simulation and checking waveform, and stickiness system is constantly adjusted 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:First, measured by experiment
Dynamic stress-strain curve builds the dynamic characteristic that dynamic constitutive more conforms to rock, measures the viscosity coefficient of rock more
Adduction is managed;Secondly, the viscosity coefficient obtained by the method for the continuous trial and error of numerical simulation is compared to dominant frequency correspondence viscosity coefficient more
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 checking waveform is not required for whether having superposition and complete, therefore can shorten the length of rock bar in experiment,
Be conducive to the test of natural rock viscosity coefficient.
Brief description of the drawings
Fig. 1 is Kelvin-Voigt model schematics;
Elastic damage constitutive relation schematic diagram during Fig. 2 microscopic element single shaft tensions;
Fig. 3 is one-dimensional viscoelasticity stress wave analysis program flow diagram.
Embodiment
It is including following the invention provides a kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple
Step:Set up the governing equation that the viscoelasticity constitutive equation and one-dimensional viscoelasticity ripple of rock bar are propagated;Passed based on one-dimensional viscoelasticity ripple
Governing equation, transflector theory and the viscoelasticity defect theory broadcast design one-dimensional viscoelasticity stress wave analysis program;Utilize one
Viscoelasticity stress wave analysis program simulation impact test is tieed up, the strain wave that input loads ripple and obtains simulation impact test is with testing
Card ripple is contrasted, and makes the strain waveform of simulation consistent with checking waveform by adjusting the parameter set, and then obtain rock bar
Viscosity coefficient.
Wherein, the viscoelasticity constitutive equation for setting up rock bar is specially:Strain rate Hopkinson pressure bar in being loaded using pendulum
Pilot system, carries out the SHPB experiments of rock bar, tests dynamic stress-strain curve of rock bar sample.By flexible member and glutinous
Property element build constitutive model, flexible member can use spring, and stickiness element is using glutinous kettle.Fig. 1 institutes can be used during specific implementation
The Kevlin body component models shown, should with reference to the SHPB dynamics for testing the rock bar obtained as the constitutive model of viscous-elastic material
Force-strain curve sets up 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 modulus of elasticity of spring in Kevlin body respectively
With the viscosity coefficient of newton body.
Set up after the dynamic constitutive equation of rock, continuity equation, the equation of motion and the rock propagated based on one-dimensional wave
Dynamic constitutive equation, sets up the governing equation that one-dimensional viscoelasticity ripple is propagated.
From the continuity equation of one-dimensional compressional wave:
The equation of motion of one-dimensional compressional 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.
According to formula (3) (4) (5) can in the hope of one-dimensional viscoelasticity ripple governing equation:
The harmonic wave solution of equation (6) is:
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 ω
Displacement and angular frequency are represented respectively.
Phase velocity c=ω/k of stress wave propagation, is the function of frequency, can be obtained by equation (8)
Built after the governing equation that one-dimensional viscoelasticity ripple is propagated, damaged using transflector is theoretical with viscoelasticity based on said process
Hinder Theoretical Design one-dimensional wave and propagate analysis program.The design principle that one-dimensional wave propagates analysis program is specifically described below:
1st, stress wave is in viscoelastic rod interface transflector
One-dimensional harmonic wave in without interface viscoelastic rod when propagating, and harmonic amplitude only declines with propagation distance according to exponential function
Subtract, harmonic frequency and phase velocity do not change.From formula (10), the spread speed of harmonic wave is not only by the material properties shadow of bar
Ring and relevant with harmonic frequency.When one-dimension stress wane propagates the interface for running into material properties change in elastic rod, stress
Ripple increases and decreases through the waveforms amplitude at interface according to transflector relation, but frequency keeps constant, 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 the interface of material properties change is run into when harmonic wave is propagated in viscoelastic rod, harmonic wave passes through the waveform frequency after interface
Rate also keeps constant, 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)
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 represent incident respectively
Ripple, 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 is unit time, A1And A2The respectively sectional area on the interface left side and the right.
Obtained by formula (11)-(15):
From formula (16) and (17), when harmonic wave is by the interface of material and cross-section variation, the transmission coefficient of harmonic wave
It is relevant with the density of material on interface both sides, harmonic wave phase velocity and cross-sectional area with reflectance factor.Because harmonic wave is in viscoelastic rod
The phase velocity of middle propagation is influenceed by the harmonious wave frequency rate of material properties, thus harmonic wave transmission coefficient and reflectance factor by harmonic wave
Frequency influence, the harmonic wave transflector coefficient of different frequency is different.The stress wave formed by the harmonic superposition of different frequency passes through material
During the interface of matter and cross-section variation, because 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 laterally having flexible, the transflector rule that the stress wave of this and arbitrary shape is propagated in one-dimensional elastic rod is not yet
Together.
2nd, viscoelasticity ripple loading BORDER PROCESSING
Because strain and stress equation form is the same in one-dimensional viscoelasticity ripple governing equation, therefore stress or strain harmonic wave exist
Propagation law is consistent in viscoelasticity stock.But it is due to the inconsistent, it is necessary to be handled accordingly of loading border.
(1) variable boundary is answered
The loading border of model is strain controlling, it is necessary to stress be solved by strain, by this structure side of equation (3) Kevlin body
Cheng Kezhi, is solved as follows based on the particle stress that strain wave is propagated:
ε in formulat、σt、EtAnd ηtStrain, stress, modulus of elasticity and the viscosity coefficient of t particle are represented, Δ t is represented
Incremental time.The stress at particle moment can be solved by the strain at this moment and last moment.
(2) stress boundary
The loading border of model is stress loading, and model medium wave is propagated in the form of straining, it is necessary to which stress wave is turned
Strain wave is changed into, from equation (20)
The initial time of stress-wave loading, t=0, ε=0.
3rd, stickiness defect theory
In one-dimensional viscoelasticity stock viscous-elastic material unit institute's tension stress or compression be more than uniaxial tension test or
Material is damaged during compression strength, and damage includes two parts, and a part is elastic damage, by the Young's modulus for changing material
To express the influence of damage, another part is stickiness damage, and influence of the damage to material is expressed by changing viscosity coefficient.
(1) elastic damage principle
According to Elastic Damage Theory, unit is produced after damage, and material parameter will be weakened, and modulus of elasticity is according to formula (22)
Handled.
E=(1-D) E0 (22)
D is damage factor in formula, between 0 to 1;E and E0It is the modulus of elasticity after damaging and before damage respectively.
It is assumed that the elastic damage constitutive relation of microscopic element is as shown in Figure 2.When the tension suffered by unit is more than material
During tensile strength, there is tensile damage in unit, and the now strain of unit is ε, it is known that εt0≤ε<εtu, εt0For elastic limit correspondence
Strain be elastic limit strain, εtuFor the limit tensile strain of unit.Still have certain because the unit is not yet destroyed completely
Bearing capacity, and this bearing capacity can define a residual strength ftrTo characterize, the calculating such as formula (23) of damage factor
It is shown.
When the compression suffered by unit be more than material compression strength when, there is crushing wound in unit, using and tensile damage
Similar method processing 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 calculating of the elastic damage factor is strain when utilizing modulus of elasticity and damage before residual strength, material damage.
When material is viscous-elastic material, the stress of unit is made up of Elastic Term and stickiness two parts, is damaged by means of elastomeric material
The computational methods of the factor calculate 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 with reference 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 above-mentioned communication theory, transflector theory and viscoelasticity damage reason on stress wave in one-dimensional viscoelastic rod
The one-dimensional viscoelasticity stress wave analysis program of design of opinion, program circuit is as shown in Figure 3.One-dimensional viscoelasticity stress wave analysis is used again
The strain wave of program simulation impact test test, strain waveform is contrasted with checking waveform, 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 set up the geometrical model of rock bar sample, carries out mesh generation 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
It is to carry out node stress to calculate whether judging unit simultaneously occurs to damage and handle degree of impairment to want function.Post processing step
Rapid major function is stress value, strain value and the elastic mould value for exporting key node.
In pre-treatment step, the following parameter first to each unit carries out assignment:It is density, modulus of elasticity, cross-sectional area, glutinous
Property coefficient, uniaxial compressive strength and uniaxial tension test.Again to be applied on rock bar sample loading form, Loaded contact analysis with
And waveform parameter is configured.When it is implemented, strain rate Hopkinson pressure bar test system is carried out in being loaded first with pendulum
The impact test of natural short rock bar, and then obtain the waveform of two measuring points.Impacted one of them complete waveform as simulation
Loaded contact analysis and another waveform in experiment are used as checking waveform.
After the completion of pre-treatment step, stress analysis step is carried out.Stress analysis step specifically includes following steps:
One new boundary stress is applied to rock bar sample;
According to governing equation and the stress value and strain value of each node of transflector theoretical calculation;
Determine whether that unit is damaged according to failure criteria, if having damage according to stickiness defect theory and Elastic Damage Theory
Degree of impairment is calculated, and changes 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.
Stress value, strain value, the elastic mould value of key node are extracted from the result of calculation of stress analysis step, and it is defeated
Go out stress diagram, strain figure and springform spirogram.
Strain waveform figure is drawn according to the stress value, strain value, elastic mould value of key node, by strain wave and checking ripple
Contrasted, make strain waveform consistent with checking waveform by adjusting parameter, and then obtain the viscosity coefficient of rock bar.
The method of testing that the present invention is combined using laboratory test and numerical simulation, initially sets up the viscoelasticity of rock bar this structure
The governing equation that equation and one-dimensional viscoelasticity ripple are propagated;The governing equation propagated based on one-dimensional viscoelasticity ripple, transflector it is theoretical with
And viscoelasticity defect theory designs one-dimensional viscoelasticity stress wave analysis program;Simulated using one-dimensional viscoelasticity stress wave analysis program
Impact test, input loading ripple is simultaneously contrasted the strain wave for simulating impact test acquisition with checking ripple, is set by adjusting
Parameter make the strain waveform of simulation consistent with checking waveform, and then the viscosity coefficient of acquisition rock bar.
The present invention more conforms to the dynamic of rock by testing the dynamic stress measured-strain curve structure dynamic constitutive
Step response, the viscosity coefficient for measuring rock is more reasonable;The viscosity coefficient phase obtained by the method for the continuous trial and error of numerical simulation
It is more representative than dominant frequency correspondence viscosity coefficient;Only need to measure a complete ripple in the dynamic impact test of the present invention
Shape is as the Loaded contact analysis of numerical model, and another checking waveform is not required for whether having superposition and complete, therefore can shorten
The length of rock bar in experiment, is conducive to the test of natural rock viscosity coefficient.
Described above is the embodiment of the present invention, it is noted that for those skilled in the art
For, under the premise without departing from the principles of the invention, some improvements and modifications can also be made, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (10)
1. a kind of short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple, it is characterised in that including following step
Suddenly:
Set up the governing equation that the viscoelasticity constitutive equation and one-dimensional viscoelasticity ripple of rock bar are propagated;
The governing equation propagated based on one-dimensional viscoelasticity ripple, 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 loading ripple simultaneously will simulate what impact test was obtained
Strain wave is contrasted with checking ripple, makes the strain waveform of simulation consistent with checking waveform by adjusting the parameter set, and then
Obtain the viscosity coefficient of rock bar.
2. natural short rock bar viscosity coefficient method of testing as claimed in claim 1, it is characterised in that described to set up the glutinous of rock bar
Elastic constitutive model equation is specially:
Strain rate Hopkinson pressure bar test system carries out the SHPB experiments of rock bar in being loaded using pendulum, tests it and dynamically should
Force-strain curve;
Constitutive model is built by flexible member and stickiness element, rock bar is set up with reference to dynamic stress-strain curve of rock bar
Dynamic constitutive equation.
3. the short rock bar viscosity coefficient method of testing as claimed in claim 1 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In described to set up the governing equation that one-dimensional viscoelasticity ripple propagates and be specially:
The dynamic constitutive equation of continuity equation, the equation of motion and the rock bar propagated based on one-dimensional wave, sets up one-dimensional viscoelasticity ripple
The governing equation of propagation.
4. the short rock bar viscosity coefficient method of testing as claimed in claim 1 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In the process for simulating impact test using one-dimensional viscoelasticity stress wave analysis program comprises the following steps:
Pre-treatment step:The geometrical model of rock bar sample is set up, mesh generation is carried out by model partition into multiple units, to each list
The relevant parameter of member carries out assignment, sets boundary condition and loading environment;
Stress analysis step:Carry out node Stress calculation simultaneously judging unit whether occur damage and to degree of impairment at
Reason;
Post-processing step:Export stress value, strain value and the elastic mould value of key node.
5. the short rock bar viscosity coefficient method of testing as claimed in claim 4 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In to following parameter progress assignment in the pre-treatment step:Density, modulus of elasticity, cross-sectional area, viscosity coefficient, single shaft
Compression strength and uniaxial tension test;
Setting loading environment is specially in pre-treatment step:To be applied on rock bar sample loading form, Loaded contact analysis with
And waveform parameter is configured.
6. the short rock bar viscosity coefficient method of testing as claimed in claim 4 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In stress analysis step specifically includes following steps:
One new boundary stress is applied to rock bar sample;
According to governing equation and the stress value and strain value of transflector theoretical calculation node;
Determine whether that unit is damaged according to failure criteria, if there is damage, according to stickiness defect theory and Elastic Damage Theory
Degree of impairment is calculated, and changes 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.
7. the short rock bar viscosity coefficient method of testing as claimed in claim 6 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In stickiness damage is calculated according to equation below:
η=(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;
E0Modulus of elasticity before damage, ftrFor tensile failure residual strength, fcrFor compression failure residual strength.
8. the short rock bar viscosity coefficient method of testing as claimed in claim 4 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In post-processing step is specially:
Stress value, strain value, the elastic mould value of key node are extracted from the result of calculation of stress analysis step, and output should
Try hard to, strain figure and springform spirogram.
9. the short rock bar viscosity coefficient method of testing that analysis is propagated based on viscoelasticity ripple as described in claim 4 or 8, its feature
It is, strain waveform figure is drawn according to the stress value, strain value, elastic mould value of key node, strain wave is entered with checking ripple
Row contrast, makes strain waveform consistent with checking waveform, and then obtain the viscosity coefficient of rock bar by adjusting parameter.
10. the short rock bar viscosity coefficient method of testing as claimed in claim 1 that analysis is propagated based on viscoelasticity ripple, its feature is existed
In strain rate Hopkinson pressure bar test system carries out the impact test of natural short rock bar in being loaded using pendulum, obtains two
The waveform of measuring point is respectively as the Loaded contact analysis needed for simulation impact test and checking 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 true CN107014704A (en) | 2017-08-04 |
CN107014704B 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) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108195725A (en) * | 2017-12-13 | 2018-06-22 | 中国人民解放军陆军工程大学 | Particulate matter one-dimension stress wane attenuation effect assay method based on SHPB |
CN110133104A (en) * | 2019-05-30 | 2019-08-16 | 东北大学 | A kind of full age dynamic characteristics test method of obturation |
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 |
---|
W.C.ZHU 等: ""Dynamic Brazilian Test of Rock Under Intermediate Strain Rate:Pendulum Hammer-Driven SHPB Test and Numerical Simulation"", 《ROCK MECH ROCK ENG》 * |
牛雷雷 等: ""岩石试样SHPB实验过程中应力均匀性的一维波动分析"", 《东北大学学报(自然科学版)》 * |
王观石 等: ""岩体中应力波幅值随时空衰减的关系"", 《岩土力学》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108195725A (en) * | 2017-12-13 | 2018-06-22 | 中国人民解放军陆军工程大学 | Particulate matter one-dimension stress wane attenuation effect assay method based on SHPB |
CN110133104A (en) * | 2019-05-30 | 2019-08-16 | 东北大学 | A kind of full age dynamic characteristics test method of obturation |
CN110133104B (en) * | 2019-05-30 | 2021-07-30 | 东北大学 | Method for testing dynamic mechanical characteristics of full-age filling body |
Also Published As
Publication number | Publication date |
---|---|
CN107014704B (en) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhu et al. | Peridynamic modelling of blasting induced rock fractures | |
Fan et al. | A Peridynamics-SPH modeling and simulation of blast fragmentation of soil under buried explosive loads | |
Abedini et al. | Comparison of ALE, LBE and pressure time history methods to evaluate extreme loading effects in RC column | |
Hu et al. | Model-based simulation of the synergistic effects of blast and fragmentation on a concrete wall using the MPM | |
Viola et al. | Soft core plane state structures under static loads using GDQFEM and cell method | |
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 | |
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 | |
CN111159949B (en) | Calculation method and system for small clear distance tunnel excavation damaged area | |
CN106055769A (en) | Performance recognition method of bolt interfaces under different tightening torques | |
CN107014704B (en) | A kind of short rock bar viscosity coefficient test method for propagating analysis based on viscoelasticity wave | |
CN106326567B (en) | Steel-steel-polypropylene hybrid fiber concrete Plastic Constitutive parameter obtaining value method in ABAQUS | |
CN110031304A (en) | A kind of rock deformation and failure mode prediction method | |
Zhao et al. | A coupled distinct lattice spring model for rock failure under dynamic loads | |
CN105787182B (en) | A kind of spacecraft firer impacts source function near field experiment correction of peak value method | |
Yu et al. | A dynamic FEM-DEM multiscale modeling approach for concrete structures | |
CN107545078A (en) | The construction method of fluid structure interaction mode and the analysis method of storage fracture hole deformation | |
CN111062113A (en) | Novel stope mining blasting parameter comprehensive optimization method under complex filling body condition | |
Wang et al. | Numerical simulation of fracture failure characteristics of rock-mass with multiple nonparallel fractures under seepage stress coupling | |
Lin et al. | A data assimilation method for blast load prediction | |
Ahmed | Models for analysis of shotcrete on rock exposed to blasting | |
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 | |
Yang et al. | Dynamic mechanical response and fracture mechanism of rock-like materials with concealed structural planes of different filling thicknesses |
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 |