#### Abstract

This paper focuses on the shear and dilatancy behavior of wheat stored in silos with various densities and normal stresses. The goal is to find a quantitative relationship modeling the peak friction angle and maximum dilatancy angle of wheat stored in silos. A total of 48 direct shear tests were carried out to research the evolution of shear and dilatancy of stored wheat in silos. It is revealed that strength of wheat in bulk attributes to the combination of frictional and dilatant during shearing, in particular attributing to its elliptic shape. An increase in relative density enhances the peak friction angle as well as the dilation. The relationships between relative density, peak friction angle, and dilatancy angle were presented based on the tests data and Bolton’s theory. Then an advanced model is developed to evaluate the peak shear behavior of wheat stored in silos considering the dilatancy of the stored wheat. It is a practical method to predict the strength and dilatancy behavior of wheat stored in silos.

#### 1. Introduction

Grain silos are special structures in industry for storing agriculture materials, which are widely used in the whole world [1]. The failure of silos may not only cause the loss of stored material itself, but also sometimes cause an explosion. Wheat is one of the important raw materials, which provides approximately twenty-one percent of the world’s food [2]. Because of large amount of wheat production and large population in China to serve, most of the wheat in China is stored in silos for three to five years. The long-term storage system is very different from those in other countries [3]. Therefore, the storage of wheat in silos is of great concern to the government, engineers, and scientists. The mechanical properties of the stored materials are primarily important to design in the silos [4].

Janssen’ theory, found in 1895, is still the basis of most internal standards to gain the loads applied on silos [5–10]. It is simple to obtain the loads with some common material properties of the internal friction of stored materials, friction coefficient of agriculture grain to wall, etc. These common properties of grain such as frictional angle, bulk density, and friction angle between concrete and steel have been widely investigated by many researchers [11–16]. However, it is found that the “real silo loads” measuring from silos are more complex than those obtained by Janssen’s theory, which still provides new uncertainties in the design [17]. The stored materials are undergoing various forces, such as gravity loading, friction with the silo walls, which change with depth during silo filling, and discharging [18–25]. For example, the dilatancy phenomenon usually occurred during silo discharging. As a result, the numerical methods are usually adopted to accurately simulate silo loads during the storage conditions. It is necessary to investigate the stress-strain model including friction and dilatancy behavior for agricultural materials at different conditions in silos.

Direct shear tests are widely used to investigate the mechanical properties of materials, especially for soils [26–30]. However, for agriculture material the mechanical properties are still poorly understood [5], especially for the large amount and long-term stored wheat in China. In this work, the strength and dilatancy behavior of wheat is investigated using direct shear tests under different normal stresses and relative densities. The stress-dilatancy model was built based on the relationship of friction and dilatancy angle. The effects of relative densities and normal stresses were also examined.

#### 2. Experimental Procedure and Materials

##### 2.1. Materials

The wheat samples used in this project were uniformly graded winter wheat, grown in Henan Province of China. Different parameters were used to describe the wheat, including:

*(**1) Particle Shape and Particle Size Distribution.* The particle shapes were quantified by Zingg index to measure the particle shape. The length of major axis (L), minor axis (B), and its thickness (T) are used to describe the particle shape shown in Figure 2. The average Zingg index is 1.62, which indicates the rod-shaped particle for wheat. The grain size distribution curve was shown in Figure 3, and the coefficient of uniformity for wheat was obtained. It showed the wheat was very uniform and poorly grade. The main properties of wheat used were summarized in Table 1.

*(**2) Void Ratio.* The granular mass is commonly considered to be consisted of a network of solid particles, enclosing voids of varying sizes. The void ratio of samples is described as where is density of wheat particle with the same quantitative value of the specific gravity =1.34, is the density of wheat mass, and is the void ratio. In this paper, the values of and are 0.668 and 0.450 for wheat, respectively, according to SL237-1999 [31]. The detailed properties of void ratio were listed in Table 2. The density of sample can be evaluated by comparing with the value of and . Here represents the state of void ratio at loosest condition, which is achieved by quickly inverting the wheat in a container. And is measured by vibration without causing wheat particle crushing.

Specimens with different void ratio were prepared by pouring a specified mass of dry wheat into the shear box, by deposition with a spoon for loose samples, in five lifts each tamped by using a steel rammer for dense samples without destroying grain particle, until the required void ratio was achieved for each layer.

*(**3) Moisture Content.* The moisture content was obtained from wheat flour with oven drying method. Some representative crushing wheat flour sample of 10 g was placed in a dish, which was dried at 105°C in an oven for at least 3 hours. The moisture content of wheat particle when crushing into flour was 10.6% by dry basis.

##### 2.2. Testing Apparatus and Procedure

The developed shear box with circular cross section was 100 mm in diameter and 100 mm in depth as shown in Figure 1. The shear cell diameter should be at least 20 times the maximum particle size and not less than 40 times the mean particle size, as recommended by Eurocode 1 Part 4. In the study, the values of the mean particle size and the mean particle size are 5 mm and 4.02 mm, according to the particle distribution curve, respectively, which are partly satisfying the size of cell box in Eurocode 1. However, the size of shear box is not very accurate with the requirement when it exceeds 40 times mean particle size. Most of the shear boxes in the literatures are less than 100 mm because of the limitation of the shear device used. In addition, the results in a large box of 143 mm, considering the ratio of 5% between sample size and grain size, are similar to those by the authors [21]. Regardless, the test results provided by these shear devices show considerable parameters for the agricultures materials.

As shown in Figure 1, the upper part of the shear box is restrained, while the lower part is controlled by a motor to apply the horizontal shear load with displacement-controlled mode. The vertical normal load is imposed by the upper platen, which is fixed with the upper part of shear box. The opening size of the two halves of shear box is 2 mm, which is corresponding to the requirement according to Shibuya et al. [32]. During the tests the loads and displacements at axial and horizontal directions were recorded automatically by a computer-controlled data collection system.

Test samples were prepared by dry pouring method with 10 layers to obtain specified densities. For dense sample it may be vibrated or compaction without grain crushing. In this study, the influences of void ratio and normal stress on stress-dilatancy behavior were examined. The values of void ratio ranged from 0.487 to 0.668 with normal stress of 50 kPa, 100 kPa, and 150 kPa, which usually existed for stored wheat in situ silos. The shearing velocities are set as 2.4 mm/min until the horizontal displacement reaches more than 20 mm.

#### 3. Experimental Results

##### 3.1. Stress-Strain Behavior and Volumetric Response

A series of direct shear tests were carried out with a consideration of different levels of void ratios and normal stresses. The initial void ratio of specimen ranged from 0.487 to 0.668 including loose and dense specimens. The typical curves of three void ratios are shown in Figures 4–6 with normal stresses ranged from 50 kPa to 150 kPa. The stress-strain evolution was plotted in terms of ratio between shear strength and normal stress R= and horizontal displacement . The volumetric response was characterized as vertical normal displacement and horizontal shear displacement. For the dense wheat sample with relative density at normal stress of 50 kPa, increases with the increase of horizontal displacement* u*, while it reached the peak value at about* u*=5 mm. After that the stress ratio keeps constant at a relatively small range until it decreases to a stable state at large horizontal displacement, which is called a critical state. During the shearing process, significant dilative deformation occurs; then a small contraction occurs in the initial stage. The maximum amount of dilation reaches the peak strength state. After the plateau stage of curve in Figure 4, the rate of dilation decreases with continuous shearing and then the critical stage is achieved. With the void ratio increasing, both the peak value of and dilative deformation reduce with the generation of horizontal displacement. Furthermore, dilation nearly vanishes for the loose sample with at 150 kPa.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

Comparing the responses of wheat at different initial void ratios, the loose sample requires larger horizontal displacement to obtain the peak horizontal shear stress than that of dense sample at the same normal stress. However, the post-peak strain softening of loose sample is gentler than that of dense sample. As expected, the critical state is confirmed using the results as shown in Figures 4–6 for all the test samples. The effect of normal stress on stress-strain response is also examined. As shown in Figures 4–6, at the same initial void ratio, the final vertical displacement decreases with the increasing normal stress, and this suppresses the dilation. The samples with different void ratio show the same trend.

##### 3.2. Peak Strength Friction Angle

The peak state represents the maximum shear strength of the wheat sample, which was determined when the horizontal shear stress meets its maximum value. The peak direct shear friction angle is defined with where and denote horizontal shear stress and normal stress.

Two methods can be used to obtain the strength friction angle at the peak state, which are interpreted of a single test (first method) or multiple tests (second method). The first method is the direct measurement from shear stress ratio of mobilized shear stress to normal stress at peak state. The results obtained from the first method are shown in Figure 7. It shows that an increase of peak friction angle performs with the increase of relative density. However, the peak friction angle of wheat decreases with the increase of normal stress. Figure 8 shows an example of the second method to obtain the friction angle with multiple tests curves of shear stress versus normal stress. The intercept equals , in which *φ* is also called internal friction angle. It is one of the useful mechanical properties in silo design as recommended by most of silo codes, which is greatly affected by the density of sample. Figure 9 showed the effect of relative density on internal friction angle at peak state strength. A gradual increase occurs with the increase of relative density of samples, which indicates denser sample exhibits higher internal friction angle and thus higher shear strength.

The shear behavior of wheat was similar to the angular sand, which was also influenced by its angularity particle shape [33]. As is well-known to us all, wheat grain is elliptic shape, which may cause large interparticle locking between the particles when they slide or rotate during shearing. High shear stress occurs to overcome the locking between particles, just before the dilation takes place. Thus, it also implies that the shear resistance is the combined action of particle interlock and dilation. In the initial stage of shearing, the interparticle locking dominates while dilation tends to play prevailing effect after a coupled effect of dilation and interlocking. The plateau curve in Figure 6 describes the consequence of couple effect. More specially, the interlocking action decreases with the development of horizontal displacement. It is noted that these observations are very similar to those of sands [33–35].

##### 3.3. Critical State Friction Angle

In this paper the friction angle at critical state is used to represent the minimum shear strength condition during shearing. It is fundamental to describe the shear behavior at constant volume shearing condition with the theory of dilatancy. In this paper, the critical state friction angle is obtained based on the fitted line of against according to Simoni and Houlsby [26]. Here, is the peak friction, while is measured at the point corresponding to peak strength. The intercept of the fitted line is the derived critical friction angle at the point of dilatancy angle as shown in Figure 10. At least two tests are needed at various densities for this method, which shows more reliable than that obtained from single test results [26]. The result is also compared with Bolton’s dilatancy equation [36, 37]. In this study, is derived from a wide range of densities of wheat direct shear tests. The test results of peak friction and dilatancy angle are regressed as shown in Figure 10. A linear fitted line with the direct shear tests of wheat in bulk is observed. The derived value for wheat at different normal stresses shows substantially different values. The values are , , and for normal stresses 50 kPa, 100 kPa, and 150 kPa, respectively. It can be seen that the critical state friction angle decreases with the increase of imposed normal stress. The observation of derived critical state friction related to the normal stress is probably caused by the effect of deformable wheat particle at large normal stress, which is somewhat similar to that of rockfill materials [38].

**(a)**

**(b)**

**(c)**

#### 4. Dilatancy Behavior of Wheat

As showed in Figures 4–6, the vertical displacement and horizontal displacement are also measured as well as the shear strength. In the direct shear stress condition, the rate of dilation of wheat is directly calculated by the ratio of vertical displacement to the horizontal displacement , which is shown in

The differential value of should be smoothing as the measured displacement is obtained at very short interval of time. The successive readings are tested for several time until the appropriate moving average value of three readings is obtained, which shows a limit scatter and sufficient smooth relation between dilation ratio and shear displacement. It is expected that value of is coincident with the value of . It is also used to calculate the maximum dilation angle .

A group of test results of maximum dilatancy angles vary with densities and normal stresses, as shown in Figure 11. It is revealed that dilatancy angle increases with an increase of relative density, while it decreases with increasing normal stress. As an example, for specimen with , the average dilatancy angle of , and for that is under imposed normal stress of 50 kPa. For the specimen with , the dilatancy angles were , , and under normal stress of 50 kPa, 100 kPa, and 150 kPa, respectively. It can be indicated that the dilation may be suppressed with the increase of imposed normal stress. The stress and dilation behavior of wheat in bulk is consistent with those reported results of various sands [26–28].

#### 5. Discussion

Bolton [37, 38] presented a new empirical equation based on Rowe’s theory, which is widely used for various sands shown in

In this paper, the framework of Bolton’s equation is adopted to interpret the experimental result of direct shear tests of wheat in bulk.

Figure 12 presented the relationship of peak friction angle and the maximum dilation angle in direct shear tests. As seen from Figure 12, the test data of wheat in bulk is fitted using (5) for wheat based on Bolton’s equation.

From the fitted curve* b*=0.85 for normal stress of 50 kPa is matched.

A linear relationship between and is observed in Figure 12 and (5) from experiment investigations. The presentation in (5) is deduced on the framework of Bolton [36, 37], which agrees with (4) by Bolton for the plane strain conditions. As above descriptions, the theory for plane strain conditions proposed by Bolton [36, 37] is also suit for the case of direct shear tests, for example, of wheat direct shear tests. Thus, the equation can be deduced in terms of the friction angle obtained from direct shear tests. The series of scattered test data of wheat is fitted with the deduced equation according to Bolton’s theory.

The influence of density, confining pressure on peak friction angle, was also examined by Bolton [36, 37]. Here, the widely accepted definition of compaction is used in (6), which is named as relative density . In this theory, a new index called relative dilatancy index is developed, which in turn combines and normal stress with its shear resistance and corresponding peak friction angle.

In this work, the effects of relative density and normal stresses are also examined for wheat in bulk. The typical test results are plotted in Figure 13. The quantitative correlation was showed in (7) and (8), where the parameter Q value is 7.45 to 7.64 with an average value of 7.54 for wheat.

#### 6. Conclusions

The shear and dilatancy behavior of wheat stored in silos are presented in this paper. The direct shear tests of wheat are carried out at various levels of densities and normal stresses in order to investigate the influence of relative density and normal stress on shear and dilation behavior of wheat. Both the peak friction angle and internal friction angle increase with the increasing relative density of wheat sample. It is observed that denser sample with higher relative density enhances the dilation by higher dilatancy angle, whereas the dilation is suppressed with increasing normal stress. Based on the framework of Bolton’s theory, the relationship between peak friction angle and maximum dilatancy angle is deduced with a consideration of relative density effect.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

This research was substantially supported by Plan for Youth Core Teachers of Henan University of Technology (2015004), Plan for the Key Projects of High Schools of Henan province (19A560009), and Provincial Key Laboratory for Grain and Oil Storage Facility & Safety, Henan University of Technology (2016KF-B03).