With the advantage of simple operation, environmental protection, land saving, less loss and the like, Silo structure has become a universal structure widely used in industry, enterprises and transportation departments [1, 2, 3]. At present, domestic and foreign scholars have done a lot of theoretical and experimental research on static characteristics of silo structure [4, 5, 6, 7]. Pieraccini et al. studied the buckling of a silo with different structure types caused by the impact load during the storage loading and unloading process, and several different models were used to carry out experimental research and numerical simulation analysis. K Knebel et al.  studied the bending buckling problem of full bunker. The results showed that the compressive stress of bulk can greatly increase the buckling critical load of thin shell wall. C.Butenweg et al.  compared two methods of silo force analysis under static and dynamic equivalent load. The interaction between silos and storage materials had been considered in the analysis. Through geometric nonlinear analysis, Ji  researched systematically on the influence of cylindrical support steel silo’s geometric parameters such as diameter and thickness of columns, steel stall, column height, radius and ratio of pillar to width on structural stability and strength. On the basis of the interaction of silos and bulk materials, the dynamic behavior of the silo system under the earthquake action is evaluated by Durmuş, A. et al. , and a simplified method for the approximation of the finite element model is proposed for the analysis. The results show that the method has high accuracy. However, most of these studies only relate to the engineering mechanical properties of a single silo structure[13, 14, 15, 16], lack of the research on silo storage foundation interaction system under dynamic load. and the related mechanism is not clear, but seismic action is one of the main factors that cause instability and failure of silo structure [17, 18, 19]. Therefore, it is of great practical significance to study the dynamic characteristics of silo storage foundation interaction system under earthquake load.
In this paper, the silo-storage-foundation system is taken as the research object. Through the Shaking table simulation test, the dynamic characteristics of the silo-storage-foundation system under different storage conditions are tested, The seismic design of the existing silo-storage-foundation system is carried out, and the effectiveness of the aseismic device is verified by numerical calculation. The research results can provide some theoretical basis and reference value for the design of safe operation of silo structures under special conditions.
2 Basic equations of dynamic model
Silo storage foundation interaction system can be regarded as a composite thin-walled cylindrical shell with one end fixed and another end free. In the role of external periodic load changes, the cylindrical shell diagram and its stratification shown in Figure 1 cylindrical shell in the movement process will have a displacement of each point in the cylindrical coordinate system with u, v, w and the three displacement components are simultaneously functions of the spatial coordinates x, θ, r and time t.
2.1 Geometric equation
The strain component at any point of a thin-walled cylindrical shell has the following relationship with the mid-plane strain, mid-plane bending deflection, and mid-plane distortion:
In the formula: are the film strain components of the middle curve, χx, χθ is the mid-plane bending deflection, χxθ is the mid-plane twist, and z is the distance from any point on the shell to the middle.
According to the Donnell shell theory [20, 21], the first-order derivative nonlinearity of the normal deflection is taken into account in the relationship between the mid-surface strain and the mid-surface displacement.
Of these the underlined items represent the nonlinear terms.
2.2 Physical equations
Taking into account the elastic modulus of the composite material with the vibration frequency changes, the two have the following relationship:
Physical equation of layer K of isotropic laminated shell:
Among them is reduced stiffness matrix, the elements given by:
Among them: Ek (ω) is k-th elastic modulus, μk is k-th layer of Poisson’s ratio.
2.3 Dynamic balance equation
According to the principle of D’Alembert, the dynamic equilibrium equation of a laminated composite cylindrical shell can be established:
3 Model test of silo-storage-foundation dynamic interaction system
3.1 Project examples
Shenyang Jinshan Thermal Power “To the small” heating project is located in Sujiatun District, Shenyang City Hunhe South Bank. Its coal handling system needs 5 silos of 10000 tons, with a diameter of 22m and a height of 39.730 m. According to the specifications of the proposed silo construction project for the warehouse. Coal is lignite, bulk density by 10 kN/m3. Seismic fortification intensity of 7 degrees, the design of the basic seismic acceleration 0.10 g, Seismic structural measures at 7 degrees fortification, the first group of seismic sub-group design, seismic grade two; structural safety rating of two; fire rating of two; silo foundation design grade B, design life of 50 years. The base bearing layer is medium sand, and the silo wall is made of reinforced concrete.
3.2 Test materials
The model test silo material used particulate concrete and galvanized wire mesh instead of the reinforced concrete material of the prototype silo. In order to simulate the interaction between the soil and the structure, a foundation soil model was prepared using spun yarn, calcium carbonate, and gypsum to ensure a reliable connection between the model structure and the vibration table. The configuration of the foundation soil was similar to the elastic modulus and density parameters of the original soil. The original storage was lignite, and similar fine-grained coal particles were used as similar materials for the storage.
3.3 Model test similarity criteria
The similarity ratio of the experimental model is selected as 1:20, the total height of the prototype is 40 m, and the total height of the model is 2.0 m. During the construction of the model, three particle concrete cubic test blocks (70.7 mm×70.7 mm×70.7 mm) and three prism test blocks (100 mm×100 mm×300 mm) were simultaneously produced. After 28 days of curing, material performance tests were conducted to test the actual compressive strength and elastic modulus of the particulate concrete material that determine the model similarity coefficient.
The test measured the elastic modulus of the micro-particle concrete was 9.33 GPa, the bulk density of the micro-particle concrete was 24.9 kN/m3, and the geometric similarity coefficient was 1:20. Based on these three parameters, similarity coefficients of other major parameters can be derived based on the similarity criteria. The following relationship can be derived:
Finally, the model similarity coefficient is determined as shown in the Table 1:
3.4 Test equipment
The vibration table used in this test is a 3 m×3 m electro-hydraulic servo shaker with a holding time of 5 minutes, Equipped with GS2000 servo controller, digitally controlled dynamic data acquisition system. Text model and equipment shown in Figure 2.
The test and measurement instruments used in this test are mainly acceleration sensors. According to the mechanical characteristics of silo-storage-foundation interaction, a total of 14 accelerometers are arranged in the key parts of the structure, as shown in Figure 3.
3.5 Selection and adjustment of seismic waves
The three factors of ground motion are frequency, amplitude and duration. When the seismic waves are selected, the characteristics of the selected seismic waves should be met according to the category of the building site, even if the main period of the design of the earthquake is close to the excellent period of the site. The original site category is class II, and the seismic fortification intensity is 7 degrees.
Considering the comprehensive consideration, 2 natural seismic waves (EL Centro wave and Lanzhou wave) and a artificial waves are selected as experimental input seismic waves. The seismic time history curves and spectral curves are shown in Figure 4.
After selecting the actual seismic wave, the amplitude of the selected seismic wave must also be adjusted. According to the regulations, the maximum value of ground motion acceleration with 7 degrees frequently earthquake is 35 cm/s2, the 7 degree fortification earthquake is 100 cm/s2 and the 7 degree rare earthquake is 220 cm/s2.
3.6 Test program
In this experiment, a multiple loading scheme is used to input two-way horizontal seismic waves at each loading. The acceleration ratio of bidirectional input seismic waves is X : Y = 1 : 0.8. Designed three conditions of empty storage, 1/2 storage and full storage. The peak acceleration of each load was gradually increased, corresponding to 7 degrees frequently earthquake, 7 degree fortification earthquake and 7 degree rare earthquake three seismic levels.
3.7 Model test results
The maximum measured acceleration value of the acceleration sensor in the base position of the silo model is taken as the reference value, and the magnification coefficient of the model under the same working condition can be obtained by comparing the measured maximum value of the acceleration sensor at the different heights of the model.
The change law of structural dynamic amplification factor is relatively complex, which is related to the characteristics of the model structure, the development of inelastic deformation, and the input seismic spectrum characteristics of the vibratory table, which can reflect the changes and destruction of the dynamic response of the model. In order to investigate the changes of the model’s dynamic amplification along its height direction, the change law of the acceleration magnification factor along the height direction of the model is analyzed. The acceleration peaks and the changing trend of dynamic amplification factor under various operating conditions are shown in Figures 5–7.
The peak acceleration measured from the foundation model shows that the input seismic load is reduced in the foundation, and the acceleration peak to the base surface is reduced to about 65%∼ 75% of the input peak. It is necessary to consider the interaction of silo foundation and consider the interaction of silo foundation, which is more suitable for engineering practice.
Under the action of seismic intensity at the design levels, the peak acceleration of the model gradually enlarges from the base to the top, and the maximum acceleration response occurs at the top of the model structure, indicating that the model structure is the first mode of vibration. At the same time, the top of the top has certain whiplash effect, which indicates the existence of higher mode shapes.
Under the condition of the same material, the acceleration magnification coefficient decreases with the increase of seismic intensity, indicating that the structure gradually enters the nonlinear stage, and the stiffness of the structure gradually decreases. Under the action of 7 degree rare earthquake intensity wave, the acceleration magnification coefficient has a larger decrease, indicating that the damage inside the silo model begins to occur.
With the increase of storage materials, the acceleration amplification coefficient of the model clearly decreases. The X acceleration magnification factor of the 7 degree rare warehouse top is 4.60, 4.01 and 2.56 respectively in empty storage, 1/2 storage and full storage. This indicates that the movement of the storage material is not synchronous with the movement of the silo body, and there is phase difference. The interaction between the silo and the storage tank plays a role of damping energy dissipation for the model structure.
This research was substantially funded by the National Natural Science Foundation of China (51504123).
Xuecheng J.X.D., The Physical Model of Pile-Soil and Structure- Pile-Soil Interaction System, J. Vibr. Eng., 1993, 2, 005. Google Scholar
Zongda J.X.Y., Earthquake Response Analysis of Building Foundation Building Interaction System, J. Vibr. Eng., 1998, 1, 23-28. Google Scholar
Ayuga F., Some unresolved problems in the design of steel cylindrical silos, Int. Conf. on structures and granular solids-from scientific principles to engineering applications, The Royal Society of Edinburgh, Scotland, UK, 2008, 6, 123-133. Google Scholar
Toma A.E., Atanasiu G.M., Specific Behavior of Circular Metal Silo, Advanced Engineering Forum, Trans Tech Publications, 2017, 21, 65-72. Google Scholar
Yu X., Raeesi A., Ghaednia H., Heydariha J., Das S., Xie S., Behavior of a Large Steel Field Silo Structure Subject to Grain Loading, J. Perf. Constr. Fac., 2017, 31(5), 401-407. Google Scholar
Pieraccini L., Silvestri S., Trombetti T., Refinements to the Silvestri’s theory for the evaluation of the seismic actions in flat-bottom silos containing grain-like material, Bulletin Earthquake Eng., 2015, 13(11), 3493-3525. Web of ScienceCrossrefGoogle Scholar
Yu J., Study on stability behavior and strength of column supported steel silos with engaged column, 2005, Zhejiang University, (in Chinese) Google Scholar
Durmus A., Livaoglu R., A simplified 3 DOF model of A FEM model for seismic analysis of a silo containing elastic material accounting for soil–structure interaction, Soil Dynam. Earthquake Eng., 2015, 77, 1-14. CrossrefGoogle Scholar
Wang X., Yang Z., Shu X., Feng J., The static contact statuses between granularmaterials and flat-bottomed steel silos, Powder Techn., 2013, 35, 1053-1059. Google Scholar
Kanyilmaz A., Castiglioni C.A., Reducing the seismic vulnerability of existing elevated silos by means of base isolation devices, Eng. Struct., 2017, 143, 477-497. Web of ScienceCrossrefGoogle Scholar
Khellat F., Khormizi M.B., A global solution for a reaction diffusion equation on bounded domains, Appl. Math. Nonlin. Sci., 2018, 3(1), 15-22. Google Scholar
Palermo M., Stefano S., Trombetti T., On the Fundamental Periods of Vibration of Flat-Bottom Ground-Supported Circular Silos containing Gran-like Material, Procedia Eng., 2017, 199, 248-253. CrossrefGoogle Scholar
Simón C.B., López J.C.C., Altuzar L.V., Micó R.J.V., Mean square calculus and random linear fractional differential equations: Theory and applications, Appl. Math. Nonlin. Sci., 2017, 2(2), 317-328. Google Scholar
Dewasurendra M., Vajravelu K., On the Method of Inverse Mapping for Solutions of Coupled Systems of Nonlinear Differential Equations Arising in Nanofluid Flow, Heat and Mass Transfer, Appl. Math. Nonlin. Sci., 2018, 3(1), 1-14. Google Scholar
About the article
Published Online: 2018-08-02
Citation Information: Open Physics, Volume 16, Issue 1, Pages 435–440, ISSN (Online) 2391-5471, DOI: https://doi.org/10.1515/phys-2018-0060.
© 2018 Li Shuwei et al., published by De Gruyter. This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License. BY-NC-ND 4.0