The harmonic vibration characteristics of a deeply buried spherical methane tank in viscoelastic soil subjected to cyclic loading in the frequency domain are investigated. The dynamic behavior of the soil is described based on the theory of fractional derivatives. By introducing potential functions, the closed-form expressions for the displacement and the stress of the viscoelastic soil surrounding the deeply buried spherical methane tank are obtained. Two die structures are considered: a homogeneous elastic medium and a shell structure. Based on the theory of elastic motion and the Flügge theory, analytic solutions for the dynamic responses of the spherical methane tank in a fractional-derivative viscoelastic soil are derived explicitly. Analytic solution expressions of the undetermined coefficients are determined by using the continuum boundary conditions. The system dynamic responses to the homogeneous elastic medium and the shell structure and the influences of the parameters of the fractional derivative, soil, and die on the dynamic characteristic of the system are compared and analyzed. The results indicate a significant difference between the dynamic responses of the die structures for the two models.
Department of Construction of Zhejiang Province Science and Technology Project2015K31National Natural Science Foundation of China51708175Business Funds of Central Colleges and Universities2018B008141. Introduction
Not only has an underground, deeply buried spherical digester the advantages of a small spherical mouth, good airtightness, and high gas production but also the fermentation of biogas slurry does not produce a crust and the structure is not affected by the external ambient temperature [1]. Therefore, spherical underground digesters are widely used and are promoted in China’s rural areas.
The vibration response of biogas digesters has wide applications in the fields of seismic resistance, explosion-resistant design, and structural stability of the die. Forrestal and Sagartz studied the dynamic responses of thin-walled spherical shells and cylindrical shells in homogeneous elastic soil under the action of heavy loads [2]. Duffey obtained the transient dynamic response of viscoelastic and viscous-plastic spherical die structures under impulsive loading by using a Laplace transform technique [3]. Glenn regarded the soil as a linear elastic medium and deduced the vibration response of a deeply buried spherical digester under heavy loads and pulse loads [4]. The transient response of a spherical cavity has been studied under explosive loading by using displacement as an unknown quantity [5].
Many experimental studies have been carried out to investigate the clay behavior under dynamic loading; Lin Guo studied the deformation behavior of saturated soft soil under long-term cyclic loading and established two equations for the prediction of long-term resilient modulus and permanent strain [6, 7]. Strain and pore pressure development on soft marine soil have been studied under dynamic loadings and an empirical formula is established to predict the long-term permanent strain [8]. Chuan Gu studied the influence of cyclic loading history on small strain shear modulus of saturated soils by triaxial test; test results show that both cyclic strain history and cyclic stress history will induce reduction of Gmax compared to the corresponding Gmax values with noncyclic loading effects at the same effective stress [9]. Dynamic modulus characteristics of saturated soils under variable confining pressure have been studied by cyclic triaxial tests; tests results shows that both the phase differences and ratios between the cyclic deviatoric stress and cyclic confining pressure exerted a constant and pronounced influence on dynamic modulus [10–12].
Studies of the viscoelastic properties of soils have shown that a fractional-derivative constitutive model is well suited for the characterization of the soil’s creep properties [13]. However, based on the fractional-derivative model, the dynamic equation is a singular integral-partial differential equation and its theoretical analysis and numerical calculation present many difficulties [14] and there have been few reports on the use of this approach for the study of soil dynamic characteristics [15]. What is more, the above studies ignore the thickness of the mold structure. In this paper, the soil element is regarded as a viscoelastic element with fractional-derivative properties and the exact expression of the displacement and stress of the infinite fractional-derivative viscoelastic soil surrounding the biogas digester is deduced based on the viscoelastic theory. In our research, the thickness of the mold structure also can be considered.
The gas fermentation and user consumption cause repeated changes in pressure in the pool and this pressure is uniformly applied to the inner side of the pool wall and is equivalent to the periodic load [16]. In this study, the soil body is regarded as a viscoelastic body with fractional-derivative properties and the exact expression of the displacement and stress of the infinite fractional-derivative viscoelastic soil surrounding the biogas digester is deduced based on the viscoelastic theory. Secondly, the die structure is considered a homogeneous elastic medium and thin-walled shell structure and the analytics solution for the vibration response of the die structure is obtained based on the theory of elastic motion [17] and the Flügge theory [18]. Finally, the concrete expression of the relevant undetermined coefficients is obtained according to the boundary conditions of displacement and stress in the mold structure and the contact surface of the soil body. A comparison of the results of this study and existing data demonstrates the correctness of the assumptions of this study. Finally, the influences of the parameters of the fractional derivative on the vibration response of the system are investigated, providing a theoretical basis for the design and construction of underground spherical digesters.
2. Mathematical Model
Because the spherical digester is deeply embedded in the soil, it is assumed that there is a spherical digester with an inner radius R1 and an outer radius R2, as shown in Figure 1; the distance from the center to the middle surface is a. The spherical digester has a circular frequency ѡ of a uniform radial harmony excitation load (q0eiωt), where i is the imaginary unit. The shear modulus parameters of the soil and die structure are GS and GL, respectively, the material density parameters are ρS and ρL, respectively, and the Poisson’s ratios of the soil and the model are νS and νL, respectively. The mold and the soil interface are in complete and close contact to meet the stress and displacement continuity.
A spherical biogas digester in viscoelastic soil.
3. Establish and Solution of the Soil Control Equation
In the case of spherical symmetry, the stress-displacement constitutive relationship of the constitutive model with the fractional derivative is [11](1)1+τεαDασrS=1+τσαDαλS∂urS∂r+2urSr+2μS∂urS∂r1+τεαDασθS=1+τσαDαλS∂urS∂r+2urSr+2μSurSr1+τεαDασφS=1+τσαDαλS∂urS∂r+2urSr+2μSurSrWhere σr is radial stress of soil layer, σθ is circumferential stress of soil layer, σφ is shear stress of soil layer, σ is stress tensor of soil, ε is strain tensor of soil, τε and τσ are constitutive model parameters, λs and μs are Lame constant of soil skeleton and λS=2νS/(1-2νS), and ν is Poisson’s ratio. Dα=dα/dtα is α-step Riemann-Liouville fractional derivative, which is defined as(2)Dαxt=1Γ1-αddt∫0txτt-ταdτin which Γ(u)=∫0∞tu-1e-udt is Gamma Function.
Under the action of dynamic load, according to the theory of elasticity, the dynamic equation [16] of soil is as follows.(3)∂σrS∂r+2σrS-σθS-σφSr=ρS∂2urS∂t2
Substituting (1) into (3), we obtain(4)λS+2μS1+τσαDα1+τεαDα∂∂r∂urS∂r+2urSr=ρS∂2urS∂t2
Introduce potential function(5)urS=∂∂rφSr,tr
For steady-state vibration, let(6)urS=R2UηSeiωt,φS=R23φ-Seiωt
Introduce dimensionless quantity and constant as follows:(7)η=rR2,λ=ωR2VS,Tσ=τσVSR2,Tε=τεVSR2,η0=1-δ,δ=dR2.
Substituting (5), (6), and (7) into (4), we obtain (8) as follows:(8)d2φ-Sdη2-β2φ-S=0
where(9)β2=-1+Tεαiλα1+Tσαiλα1-2vSλ221-vS
According to (8), we obtain (10)φ-S=Ae-βη+Beβη
where A and B are undetermined coefficient; they are determined by boundary conditions. According to boundary conditions r→∞,urS=0, we obtain B=0,(11)φ-S=Ae-βη
Substituting (5) into (11), we obtain (12) as follows:(12)UηS=-1η2+βηe-βηA
Substituting (12) into (1), we obtain (13) as follows:(13)σrS=μS1+Tσαiλ1+Tεαiλeiωt2vS1-2vSβ2η+22η3+2βη2+β2ηe-βηAσθS=μS1+Tσαiλ1+Tεαiλeiωt2vS1-2vSβ2η-21η3+βη2e-βηA
4. Establishment and Solution of Die Structure Motion Equation4.1. Die Structure Equation of Type 1
We considered the mold structure of the spherical digester as a homogeneous elastic medium and assumed that the mold structure and soil interface are in complete and close contact with no relative slip to satisfy the stress and displacement continuity. Under the action of a dynamic load, the motion equation of the die junction is as follows [16].(14)∂σrL∂r+2σrL-σθL-σφLr=ρL∂2urL∂t2
Constitutive model is(15)σrL=λL∂urL∂r+2urLr+2μL∂urL∂rσθL=λL∂urL∂r+2urLr+2μLurLrσφL=λL∂urL∂r+2urLr+2μLurLr
Substituting (15) into (14), we obtain (16) as follows:(16)λL+2μL∂∂r∂urL∂r+2urLr=ρL∂2urL∂t2
For steady-state vibration, let urL=R2UηLeiωt, φL=R23φ-Leiωt. Introduce displacement potential function(17)urL=∂∂rφLr,tr
Substituting (17) into (16), we obtain(18)d2φ-Ldη2-q2φ-L=0
where(19)GSL=GSGL,ρLS=ρLρS,q2=-1-2vLλ2ρLSGSL21-vL
According to (18), we obtain(20)φ-L=C5e-qη+C6eqη
where C5 and C6 are undetermined coefficient; they are determined by boundary conditions. Substituting (20) into (17), we obtain(21)UηL=-1η2+qηe-qηC5-1η2-qηeqηC6
According to (15), we obtain radial stress amplitude(22)σrL=2vLGL1-2vLq2η+2GLϑ1e-qηeiωtC5+2vLGL1-2vLq2η+2GLϑ2eqηeiωtC6
where(23)ϑ1=2η3+2qη2+q2η,ϑ2=2η3-2qη2+q2η
Satisfy the boundary conditions as follows:
(1) Displacement continuous in the mold structure and soil interface (r=R2, η=1); that is,(24)urS=urLr=R2
(2) Stress continuity in the mold structure and soil interface (r=R2, η=1); that is,(25)σrS=σrLr=R2
(3) Internal hemorrhage loading equal to the radial stress in the mold structure and soil interface (r=R1, η=1) of the die structure; that is,(26)σrL=q0eiωtr=R1
According to (24)~(26), we obtain the actual formulas of the undetermined coefficients A,C5,C6,(27)C6=a2b1-a1b2e1a1b3-a3b1+e2a2b1-a1b2QC5=Qe1-e2e1C6A=-a2a1C5-a3a1C6,
Since the thickness of the die structure is much smaller than the radius of the spherical digester, it can be assumed that the mold structure has little effect on the interface r=R2 of the soil and the middle surface radius [6]r=a. According to the Hamilton principle, the motion equation of the die structure [7] is as follows:(29)urL2GLd1-vLR22+ρLd∂2urL∂t2=Qa
where Qa is net pressure.
For steady-state vibration, we assumed(30)urL=R2UηLeiωt,Qa=GSQ-aeiωtGSL=GSGL,ρLS=ρLρS
Meanwhile, satisfy the boundary conditions as follows:
(1) Stress compatibility in the mold structure and soil interface (r=R2, η=1); that is,(31)GSLQ-aeiωt=Qeiωt-σrSGL
(2) Displacement continuous in the mold structure and soil interface (r=R2, η=1); that is,(32)urS=urLr=R2
Substituting (29) and (13) into (31) and (32), we obtain(33)A=QGSL1+Tσαiλα/1+Tεαiλαϑ3-ϑ41+βe-β
In order to verify the correctness of the calculation, we compared the results with existing results, which include the ideal elastic model of [16], the classical elastic constitutive model [17], and classical viscoelastic constitutive model [18].
The ideal elastic model of [16] ignores the influence of the die structure firstly. Besides, it establishes the dynamic characteristics of the spherical digester in a general and ideal elastic soil. Substituting the equation for the classical elastic constitutive method into (3) and assuming that the material parameters of the soil are the same as in the die structure, we determined that the results of this study were in agreement with other published solutions [16]. The parameter values were as follows: η=1, vS=0.35, ρLS=1.5, GSL=0.05, vL=0.25, δ=0.005, Tσ/TƐ=3, TƐ=10, and α=0.5. Figure 2 shows a curve for the radial displacement and hoop stress amplitude with a dimensionless frequency for the two models of a uniform elastic medium. It can be seen that, in [16], the resonance effect is significantly greater for the model than for the homogeneous elastic dielectric model and that the base frequency is also higher than the base frequency of the resonance effect of the system. There are two main reasons for this: firstly, the author uses the fractional-derivative constitutive model to describe the soil viscosity and [16] uses a general and ideal elastic soil body; secondly, the author considers the influence of the die structure which has a greater stiffness than the soil body.
Comparison between present results and [16] results.
Radial displacement amplitude
Hoop stress amplitude
Thirdly, the fractional-derivative constitutive model used in this study is compared with the classical elastic model [17] and the classical viscoelastic constitutive model [18]. Substituting the classical elastic constitutive model and classical viscoelastic constitutive model used in the literature [19] into the soil equation of motion (3), respectively, changes the results of [19, 20]. Figure 3 shows the changes in the curves of the radial displacement amplitude and circumference displacement amplitude as a result of the changing frequency of the three models. It can be seen that, for the elastic model in [20], the peak value of the displacement amplitude of the soil is maximum; the resonance effect is clear and the resonance phenomenon decreases with increasing. The peak value of the classical viscoelastic model used in [19] is the minimum, while, for the fractional-derivative viscoelastic model, it is in the middle.
Comparison analysis of three models.
Radial displacement amplitude
Hoop stress amplitude
6. Analysis of Example and Graphic Analysis
The effects of the shear modulus ratio GSL of the soil and die, the Poisson ratio νs of the soil material, the fractional-derivative order α, the material parameter ratio Tσ/Tε on the radial displacement, and the hoop stress amplitude were investigated.
Figure 4 shows that the values of the radial displacement and circumferential stress amplitude change with the frequency f as the shear modulus ratio GSL of the soil body and die structure variation change. It is evident that when the modulus ratio GSL is comparatively small, the displacement amplitudes are very different for the model and the second type of mold model and the resonance effect produced by the latter system is significant. With the increase in the modulus ratio λ, the difference in the displacement amplitude for the two models is small. The stress amplitude is clearly greater for the spherical shell die structure than for the uniform elastic medium. As the frequency increases, the difference between the two models decreases and the displacement and the stress amplitude decrease with an increase in the modulus ratio.
Influence of modulus ratio on radial displacement and hoop stress amplitude.
Radial displacement amplitude
Hoop stress amplitude
Figure 5 shows the effect of the changes in Poisson’s ratio νs on the radial displacement and the amplitude of the circumferential stress. As Poisson’s ratio νs increases, the radial displacement and amplitude of the hoop stress increase gradually. The difference in the displacement and stress amplitude between the two types of die structures at high frequency is very small.
Influence of Poisson’s ratio on radial displacement and hoop stress amplitude.
Radial displacement amplitude
Hoop stress amplitude
Figure 6 shows the influence of the fractional-derivative order α on the radial displacement and circumferential stress amplitude. It is evident that the amplitude of radial displacement decreases with an increase in the order number α in genera; with the increase in the frequency λ, the difference in the displacement amplitude for the two models is small. The hoop stress changed with the fractional-derivative order number α being interesting. Its change has a certain relationship with the frequency λ: when the frequency λ is between 0 and 3, the hoop stress amplitude decreases with an increase in the order number α; when the frequency λ is greater than 3, the circumferential stress amplitude increases with an increase in the order number α.
Influence of order on radial displacement and hoop stress amplitude.
Radial displacement amplitude
Hoop stress amplitude
Figure 7 shows the changes in the radial displacement and circumferential stress amplitude as a result of changes in the frequency λ when the material parameter ratio Tσ/TƐ of the constitutive model is changed. The data indicate that the radial displacement and circumferential stress amplitude decrease as the material parameter ratio Tσ/Tε increases; this is caused by the increase in the impedance of the soil body.
Influence of material parameter ratio on radial displacement and hoop stress amplitude.
Radial displacement amplitude
Hoop stress amplitude
7. Conclusions
In this paper, the simplified vibration characteristics of a deeply buried spherical digester in viscoelastic soil under uniform internal pressure in the frequency domain are studied using an analytics method; the effects of relative physical properties and geometrical parameters on the dynamic characteristics of the system are considered. The main conclusions are summarized as follows:
(1) Under the condition of Tσ/TƐ =3, the steady-state response of the viscoelastic soil is greater for the fractional-derivative viscoelastic soil than for the classical viscoelastic soil but less than for the classical elastic soil
(2) The steady-state response of the spherical digester is clearly greater for the spherical shell mold structure model than for the homogeneous elastic medium model
(3) When the modulus is comparatively small, the displacement amplitude is very different for the two types of mold type models and the difference between them decreases with increases in GSL
(4) With the increase in Poisson’s ratio νS of the soil material, the displacement amplitude and the hoop stress amplitude increases gradually
(5) With the increase in the order number α, the radial displacement and the hoop stress amplitude decrease and the system resonance effect weakens
(6) With the increase in the material parameter ratio Tσ/TƐ, the radial displacement and the hoop stress amplitude decrease
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The work presented in this paper was supported by the Department of Construction of Zhejiang Province Science and Technology Project (no. 2015K31), National Natural Science Foundation of China under (Grant no. 51708175), and Business Funds of Central Colleges and Universities (Grant no. 2018B00814).
GongW.DingM.XuQ.JiangX.Earth pressure on pool wall and stress analysis on pool cover of underground biogas digester20112742702762-s2.0-79959745370ForrestalM. J.SagartzM. J.Radiated pressure in an acoustic medium produced by pulsed cylindrical and spherical shells1971384105710602-s2.0-001520624210.1115/1.3408916DuffeyT. A.Transient response of viscous-elastic and viscous-elastic shells submerged in fluid media197698137143GlennA. L.KidderR. E.Blast loading of a spherical container surrounded by an infinite elastic medium200022984992QingyuanZ.RenruiZ.The dynamic response of spherical hollow chamber under explosive loading1994142182185GuC.WangJ.CaiY.SunL.WangP.DongQ.Deformation characteristics of overconsolidated clay sheared under constant and variable confining pressure20165634274392-s2.0-8499231196710.1016/j.sandf.2016.04.009GuoL.WangJ.CaiY.LiuH.GaoY.SunH.Undrained deformation behavior of saturated soft clay under long-term cyclic loading20135028372-s2.0-8487579229410.1016/j.soildyn.2013.01.029WangJ.GuoL.CaiY.XuC.GuC.Strain and pore pressure development on soft marine clay in triaxial tests with a large number of cycles2013741251322-s2.0-8488769747210.1016/j.oceaneng.2013.10.005GuC.WangJ.CaiY.GuoL.Influence of cyclic loading history on small strain shear modulus of saturated clays2014661122-s2.0-8490456503710.1016/j.soildyn.2014.06.027GuC.GuZ.CaiY.WangJ.LingD.Dynamic modulus characteristics of saturated clays under variable confining pressure20175457297352-s2.0-8504285327510.1139/cgj-2016-0441CaiY.WuT.GuoL.WangJ.Stiffness degradation and plastic strain accumulation of clay under cyclic load with principal stress rotation and deviatoric stress variation20181445040180212-s2.0-85041171602WuT.CaiY.GuoL.LingD.WangJ.Influence of shear stress level on cyclic deformation behaviour of intact Wenzhou soft clay under traffic loading201722861702-s2.0-8502653751010.1016/j.enggeo.2017.06.013GaoH. X.WenM. J.Dynamic characteristics of a spherical cavity in nearly saturated viscoelastic soil2012444753761MR3026274LiuL.-C.YangX.Analysis of vertical vibrations of a pile in saturated soil described by fractional derivative model20113225265322-s2.0-79953089979LijunH.LingweiK.WenjunW.A description of creep model for soft soil with fractional derivative201132supplement 2239244WeilandP.Biogas production: current state and perspectives201085484986010.1007/s00253-009-2246-72-s2.0-76649102408ChangjieX.ShimingW.Spherical wave propagation in saturated soils199819324325210.1007/BF024533892-s2.0-0032258370ChangjieX.YuanqiangC.Dynamic response of spherical cavity in viscous-elastic saturated soils20013448892ZakoutU.AkkasN.Transient response of a cylindrical cavity with and without a bonded shell in an infinite elastic medium19973513-12120312202-s2.0-003122350310.1016/S0020-7225(97)00011-6Zbl0902.73011TimoshenkoS. P.GoodierJ. N.1977New York, NY, USAMcGraw-HillMR0045547