# Neutrosophic Functions of the Joint Roughness Coefficient and the Shear Strength: A Case Study from the Pyroclastic Rock Mass in Shaoxing City, China

Authors: Jun Ye; Rui Yong; Qi-Feng Liang; Man Huang; Shi-Gui Du
Journal: Mathematical Problems in Engineering
 (2016)

---

## Abstract

Many studies have been carried out to investigate the scale effect on the shear behavior of rock joints. However, existing methods are difficult to determinate the joint roughness coefficient (JRC) and the shear strength of rock joints with incomplete and indeterminate information; the nature of scale dependency of rock joints is still unknown and remains an ongoing debate. Thus, this paper establishes two neutrosophic functions of the JRC values and the shear strength based on neutrosophic theory to express and handle the incomplete and indeterminate problems in the analyses of the JRC and the shear strength. An example, including four rock joint samples derived from the pyroclastic rock mass in Shaoxing city, China, is provided to show the effectiveness and rationality of the developed method. The experimental results demonstrate that the proposed neutrosophic functions can express and deal with the incomplete and indeterminate problems of the test data caused by geometry complexity of the rock joint surface and sampling bias. They provide a new approach for estimating the JRC values of the different-sized test profiles and the peak shear strength of rock joints.


---

## Body

## 1. Introduction

In mining, civil, hydraulic, and petroleum engineering, engineers often face problems associated with jointed rock mass [1]. One of the most challenging tasks in engineering rock mechanics is comprehensively understanding the characteristics of rock joints, including orientation, extent, planarity, roughness, and strength of wall rock asperities. Roughness refers to the local departures from planarity at both small and large scale [2] and has a direct impact on the shear strength of rock joints. It has been researched in the past four decades since Barton introduced the joint roughness coefficient (JRC) [3].A variety of researches have been carried out to quantitatively describe the surface roughness of fractured surfaces in an attempt to relate the conventional geometric parameters or fractal dimensions of surface roughness to the mechanical behavior of rock joints [4]. The most widely used formula for estimating the shear strength was proposed by Barton [5] as follows:(1)
τ
=
σ
n
tan
⁡
J
R
C
log
10
⁡
J
C
S
σ
n
+
ϕ
b
,where τ is the peak shear strength of rock joints, σ
n is the effective normal stress, ϕ
b is the basic friction angle of the discontinuity surface, and JCS is the joint wall compressive strength.The variation in the mechanical behavior of joints with increasing scale is a well-known phenomenon which makes the mechanical parameters worked out in the laboratory unsuitable for the natural discontinuities located in real condition. To better understand the deformational behavior of rock systems, a large number of studies have been carried out over the last four decades for exploring the scale effect on the shear strength of rock joints [5–9]. However, many of these studies have produced conflicting results and the rock joints of the same lithology may exhibit an extremely large range in mechanical properties [10, 11]. The uncertainties in test results are mainly caused by the complex properties of surface roughness of rock joints. Du [12] found that the rock joint roughness in nature generally has the appearance-like nonuniformity, anisotropy, inhomogeneity, and scale effect. To describe the complex structure of joint surface, Chen et al. [13] proposed a geological statistics method to analyze the anisotropy and the size effects of structural surface of rocks. Meanwhile, the sampling bias and sampling disturbance effects may be responsible for incorrect conclusions concerning some of the apparent scale effects [14]. Time and budget constraints make it unfeasible to obtain sufficient data by lots of field experiments. Due to the special limitation in measurement accuracy and the variability of roughness, it is difficult to find a certain equation for describing the scale effect of different-sized rock joint samples.Neutrosophy means the study of ideas and notions which are not true, nor false, but its study is between true and false, that is, neutral, indeterminate, unclear, vague, ambiguous, incomplete, contradictory, and so forth. Many types of indeterminacies are contained in many real situations. Therefore, the neutrosophic logic, neutrosophic set, neutrosophic probability, neutrosophic statistics, neutrosophic measure, neutrosophic precalculus, neutrosophic calculus, and so forth were born in neutrosophy [15]. Neutrosophy theory is very suitable for describing the incomplete and indeterminate information. Since there is the incomplete and indeterminate information in the JRC values of different-sized samples and the shear strength of rock joints indicated in previous researches, neutrosophic function is helpful to better express these incomplete and indeterminate phenomena and can express the relationship of the scale effect of JRC and the shear strength of rock joints. To do so, we use neutrosophic functions to describe the scale effect of JRC and shear strength and verify the effectiveness and rationality of the neutrosophic functions of the JRC values and the shear strength by real testing data for four large-scale rock joint samples in the field.The rest of the paper is organized as follows. In Section2, we introduce the some concepts of neutrosophic functions. In Section 3, the neutrosophic functions of JRC and the shear strength are established in different sampling lengths. In Section 4, a practical example based on real testing data for four test samples is given to illustrate the application and effectiveness of the developed approach. Finally, some final remarks and further work are given in Section 5.


## 2. Some Concepts of Neutrosophic Functions

Smarandache [15] firstly proposed the neutrosophic functions and neutrosophic calculus to handle various problems of indeterminacy in real world.Let us consider a neutrosophic function (thick function)f
:
R
→
F
(
R
), and then(2)
f
x
=
3
x
,
3
x
+
3
f
o
r
a
n
y
x
∈
R
,where R is all real numbers and F
(
R
) is the set of all subsets of R.Thus, its graph is shown in Figure1.Figure 1
The graph of the neutrosophic function (thick function)f
(
x
).For example, ifx
=
3, f
(
3
)
=
[
3
×
3,3
×
3
+
3
]
=
[
9,12
], which is shown in Figure 1. The grey area indicates the indeterminate range between f
1
(
x
)
=
3
x and f
2
(
x
)
=
3
x
+
3.Let us consider another neutrosophic function (thick function)g
:
R
→
F
(
R
); then(3)
g
x
=
2
e
x
,
2
x
=
2
e
x
,
2
e
2
x
for  any
x
∈
R
.Ifx
=
2, then g
(
2
)
=
[
2
e
2
,
2
e
4
].In general, we may define a neutrosophic function (thick function or interval function)f
:
R
→
F
(
R
) as follows [15]:(4)
f
x
=
f
1
x
,
f
2
x
.For example, ifx
=
1, then f
(
1
)
=
[
f
1
(
1
)
,
f
2
(
1
)
]; the value of the neutrosophic function f
(
x
) is a vertical segment of line at x
=
1, which is the interval value [
f
1
(
1
)
,
f
2
(
1
)
].However, we can also define an open interval functionf
(
x
)
=
(
f
1
(
x
)
,
f
2
(
x
)
) or semiopen/semiclosed interval functions f
(
x
)
=
[
f
1
(
x
)
,
f
2
(
x
)
)
/
f
(
x
)
=
(
f
1
(
x
)
,
f
2
(
x
)
].


## 3. Neutrosophic Functions of JRC and the Shear Strength

As we know in rock mechanics, the joint roughness coefficient (JRC) of rock joints contains the indeterminacy in different sampling lengths and directions of rock joints [16] and then results in the indeterminacy of the shear strength of joints due to the corresponding relationship between JRC and the shear strength.To study surface properties of the rock joints in Fusheng town, Shaoxing city, China, we measured 946 roughness profiles from the large, exposed pyroclastic rock mass by modified mechanical hand profilograph. The length of these samples ranged from 10 cm to 170 cm. The images of these joint roughness profiles were obtained by large format scanner. The roughness profiles were digitized with the interval of 0.5 mm based on the grayscale image processing technique. The JRC values were estimated based on the straight edge method [17]. The maximum and minimum JRC values of different-sized samples were obtained, respectively. Taking the maximum values as an example, Figure 2 shows the distribution of the maximum JRC values with the sample length. It denotes that the maximum values of JRC decrease gradually as the sampling length increases. An exponential function was proposed for describing their correlations as follows:(5)
J
R
C
+
L
=
16.35
+
60.18
×
0.95
6
L
f
o
r
L
∈
10,170
.Figure 2
The correlation between the sampling length and maximum JRC values.Similarly, the correlation between the sampling length and minimum JRC values was obtained as(6)
J
R
C
-
L
=
3.65
+
2.3
×
0.95
8
L
f
o
r
L
∈
10,170
.Then the neutrosophic function for calculating the JRC values was proposed as the following formula:(7)
J
R
C
L
=
J
R
C
-
L
,
J
R
C
+
L
=
3.65
+
2.3
×
0.95
8
L
,
16.35
+
60.18
×
0.95
6
L
f
o
r
L
∈
10,170
,where L is the sampling length and J
R
C
-
(
L
) and J
R
C
+
(
L
) are the lower and upper limits of the neutrosophic function J
R
C
(
L
). Hence, the graph of J
R
C
(
L
) is shown in Figure 2.For example, ifL
=
80 cm, then J
R
C
(
80
)
=
[
3.7243,17.9947
], which is shown in Figure 3.Figure 3
The graph of the neutrosophic function (thick function)J
R
C
(
L
) for L
∈
[
10,170
].Based on the Barton shear strength [5], we propose the following neutrosophic function of shear strength:(8)
τ
J
R
C
=
σ
n
-
,
σ
n
+
tan
⁡
J
R
C
-
L
,
J
R
C
-
L
log
10
⁡
J
C
S
-
,
J
C
S
+
σ
n
-
,
σ
n
+
+
ϕ
b
=
σ
n
-
tan
⁡
J
R
C
-
L
log
10
⁡
J
C
S
-
σ
n
+
+
ϕ
b
,
σ
n
+
tan
⁡
J
R
C
+
L
log
10
⁡
J
C
S
+
σ
n
-
+
ϕ
b
,where τ
(
J
R
C
) is the neutrosophic function of the peak shear strength of rock joints; σ
n
- and σ
n
+ are the lower and upper limits of the effective normal stress; J
R
C
-
(
L
) and J
R
C
+
(
L
) are the lower and upper limits of the neutrosophic function J
R
C
(
L
); ϕ
b is the basic friction angle of the discontinuity surface; and JCS− and JCS+ are the lower and upper limits of the joint wall compressive strength.Obviously, (7) and (8) are two interval functions (indeterminate functions) and their indeterminate values lie within the intervals [J
R
C
-
(
L
), J
R
C
+
(
L
)] and [τ
-
(
J
R
C
), τ
+
(
J
R
C
)
], respectively. Hence, they can effectively describe the indeterminate phenomena of rock joints.


## 4. Experimental Results and Discussion

To demonstrate the effectiveness of the proposed neutrosophic functions, the samples of natural rock joints were selected to evaluate the scale effect of surface roughness in this study. We measured these testing samples from the pyroclastic rock mass that crop out at the east of Shaoxing City in Zhejiang Province, China (Figure4). The pyroclastic rocks are volcaniclastic rocks formed by accumulation of pyroclasts (fragments) during Mesozoic explosive eruption by Yanshan movement. According to the size and abundance of their pyroclasts, these pyroclastic rocks are classified to be the lapillistones (pyroclasts 2–64 mm in size).Figure 4
Location of the test site in simplified geological map (GPS: N 29°57′, E 120°43′).Generally, the rock joints play an important role in mechanical properties and deformation behavior of rock masses, the shear resistance hinges on the frictional resistance, and geometric irregularities along rock joint surfaces. In this study, the rock masses are highly fractured due to the mechanical behavior of these soft rocks. Three dominant joint sets are widely developed in this study area. The safety netting system (SNS) (diamond-shaped steel wire mesh) has been applied for protecting the potential rock slope or rock falls (Figure5). Fifteen profiles of rock joint set (dip 15°, dip direction 95°) were measured in situ using modified mechanical hand profilograph [18, 19]. The length of measured samples ranges from 172.69 cm to 223.94 cm.Figure 5
Photograph showing the fractured rock mass and the measurement area of rock joint set.In the field, the length of the rock joint surface of interest is generally much longer than small-scale laboratory specimen, so JRC value must be estimated for full-scale surface. Barton and Bandis [20] proposed an alternative and simple method for estimating field-scale JRC, as shown in Figure 6. The approximate value of field-scale JRC is approximately determined from the measured values of amplitudes and lengths of joint profiles. In this study, these roughness profiles were digitized in the sampling interval of 0.5 mm, which we then used to determine JRC values. Firstly, each roughness profile was divided into small sized samples, in the length of 10 cm, 20 cm, 30 cm, 40 cm, 50 cm, 60 cm, 70 cm, 80 cm, 90 cm, 100 cm, 110 cm, 120 cm, 130 cm, 140 cm, 150 cm, 160 cm, and 170 cm. Note that the left side of the original rock joint profile along the potential sliding direction was chosen as the testing sample if the sample size is longer than half-length of the original profile. When we take the 176.91 cm original roughness profile as an example, there are 17 samples of 10 cm long, 8 samples of 20 cm long, 5 samples of 30 cm long, 4 samples of 40 cm long, 3 samples of 50 cm long, and two for each 60 cm, 70 cm, and 80 cm sized samples; there is just one sample when the profile length is over 90 cm.Figure 6
Estimating JRC from measurements of surface roughness amplitude from a straight edge [17].Here, original roughness profiles (J1, J2, J3, and J4) were provided as an illustrative example and applied for analyzing the relationship between the joint roughness coefficient and their sample sizes, as shown in Figure7. The test data and mean values of JRC for the various sizes of the rock joint samples J1, J2, J3, and J4 are illustrated in Tables 1–4, and then the graphs between the JRC values and the various sizes of the rock joint samples are shown in Figure 8. Clearly, the mean values of JRC of different-sized samples are decreased with the sample length increasing.Table 1
The test data and mean value of JRC of J1.
Length (cm)
Test data of JRC
Mean value
10
13.608, 8.196, 14.842, 15.747, 15.478, 10.121, 24.452, 20.134, 11.879, 12.280, 15.029, 8.156, 8.374, 19.908, 19.160, 34.872, 15.034
15.722
20
8.884, 15.083, 8.791, 22.140, 10.292, 9.787, 12.900, 17.436
13.164
30
8.364, 7.924, 19.591, 8.385, 12.445
11.342
40
8.172, 13.304, 9.072, 11.256
10.451
50
12.987, 12.770, 4.456
10.602
60
8.581, 10.430
9.506
70
9.473, 10.129
9.801
80
8.486, 6.988
7.737
90
9.169
9.169
100
9.263
9.263
110
9.261
9.261
120
8.490
8.490
130
7.837
7.837
140
7.277
7.277
150
6.791
6.791
160
6.746
6.746
170
7.704
7.704Table 2
The test data and mean value of JRC of J2.
Length (cm)
Test data of JRC
Mean value
10
19.388, 30.850, 9.538, 10.764, 24.210, 8.971, 22.201, 17.600, 7.785, 18.691, 9.024, 10.227, 9.925, 11.114, 19.033, 17.454, 12.663
15.261
20
25.400, 11.854, 14.536, 14.643, 12.103, 5.621, 9.347, 13.679
13.398
30
18.169, 10.002, 14.905, 7.942, 7.4270
11.689
40
13.627, 10.858, 10.810, 11.202
11.624
50
12.987, 12.770, 4.456
10.071
60
11.016, 9.105
10.061
70
10.157, 8.394
9.276
80
8.887, 6.830
7.858
90
8.130
8.130
100
9.392
9.392
110
8.700
8.700
120
8.290
8.290
130
7.746
7.746
140
7.342
7.342
150
6.878
6.878
160
7.351
7.351
170
8.007
8.007Table 3
The test data and mean value of JRC of J3.
Length (cm)
Test data of JRC
Mean value
10
57.464, 13.48, 18.888, 17.990, 17.595, 37.676, 28.627, 12.708, 8.585, 13.930, 10.742, 8.063, 18.024, 18.868, 7.812, 9.189, 19.002
18.744
20
31.402, 11.898, 18.838, 15.541, 11.540, 6.391, 12.870, 6.283
14.345
30
20.935, 13.115, 11.630, 4.936, 9.732
12.070
40
15.701, 13.389, 9.451, 8.704
11.811
50
12.561, 11.487, 7.572
10.540
60
10.744, 7.821
9.282
70
9.209, 7.298
8.254
80
8.0578, 6.182
7.120
90
9.158
9.158
100
9.135
9.135
110
8.304
8.304
120
7.612
7.612
130
7.043
7.043
140
7.009
7.009
150
7.211
7.211
160
7.265
7.265
170
7.955
7.955Table 4
The test data and mean value of JRC of J4.
Length (cm)
Test data of JRC
Mean value
10
33.265, 13.451, 16.894, 7.898, 12.925, 9.519, 17.780, 23.429, 26.439, 13.870, 24.149, 21.935, 13.749, 12.593, 16.588, 15.549, 22.397, 25.772, 14.201, 12.490, 24.323
18.058
20
20.779, 9.821, 9.593, 12.065, 19.604, 15.418, 14.186, 19.551, 17.584, 10.350
14.895
30
15.449, 9.080, 13.102, 13.148, 13.803, 12.591, 8.746
12.274
40
11.587, 11.292, 11.218, 10.786, 13.055
11.587
50
9.847, 13.752, 12.030, 10.533
11.540
60
8.243, 15.787, 11.266
11.765
70
7.144, 11.009, 12.321
10.158
80
6.452, 10.706
8.579
90
8.070, 8.326
8.198
100
7.263, 9.751
8.507
110
9.032
9.032
120
8.688
8.688
130
8.181
8.181
140
7.688
7.688
150
7.175
7.175
160
6.727
6.727
170
6.331
6.331Figure 7
Digitization of original roughness profiles.Figure 8
The relationship between the JRC values and the sizes of the rock joint samples. (a) J1; (b) J2; (c) J3; and (d) J4.
(a)
(b)
(c)
(d)Although the difference between maximum JRC value and the minimum JRC value for 10 cm long samples (small samples) is the biggest of all original roughness profiles (J1, J2, J3, and J4), we can see from Figure8 that the JRC values of the test data fall in the interval value of the neutrosophic function (thick function) J
R
C
(
L
) for L. Clearly, the neutrosophic function J
R
C
(
L
) can indicate the range of the JRC values and show its effectiveness in the JRC analysis between the joint roughness coefficient and their sample sizes. Thus, it will provide an analysis foundation in order to predict the shear strength of rock joints.To analyze the shear strength of rock joints, let us consider that the basic friction angleϕ
b is 21°, σ
n
- is 1 MPa, and σ
n
+ is 2 MPa. The joints are completely unweathered; then JCS will be equal to the unconfined compression strength of the unweathered rock σ
c [5]. The upper limit of the joint wall compressive strength JCS+ is equal to σ
c
=
10 MPa and the lower limit JCS− is equal to 0.5
σ
c
=
5 MPa.According to the neutrosophic function of shear strength in formula (6), the peak shear strength of rock joints can be predicted by the following formula:(9)
τ
L
=
tan
⁡
3.65
+
2.3
×
0.95
8
L
log
10
⁡
5
2
+
21
,
2
tan
⁡
16.35
+
60.18
×
0.95
6
L
log
10
⁡
10
1
+
21
.Then, based on (9), we can give the graph of the neutrosophic function τ
(
L
) for L
∈
[
10, 170
], as shown in Figure 9. The grey area in Figure 9 is its possible range (indeterminate area) of the shear strength of rock joints.Figure 9
The graph of the neutrosophic functionτ
(
L
) for L
∈
[
10
,
17
0
].To predict the possible range of the shear strength of rock joints, we consider the following example.Example 1.
IfL
=
80 cm, then τ
(
80
)
=
[
0.414,1.619
] (MPa), which is shown in Figure 9. Obviously, the possible value of the shear strength of rock joints for L
=
80 cm may be within the interval [
0.414, 1.619
] (MPa). This result may provide a theoretical analysis method for predicting the shear strength of rock joints in the specific district in Shaoxing city.
From the above analyses, we can see that the proposed neutrosophic functions can estimate the JRC values and the shear strength of rock joints and demostrate their effectiveness in the specific district. Therefore, the advantage is that neutrosophic functions can effectively express and analyze these inconplete and indeterminate problems because existing methods cannot express and handle these incomplete and indeterminate problems.


## 5. Conclusion

JRC is one of the most important parameters utilized in calculating shear strength of rough joints in rock masses. However, the incomplete and indeterminate problems generally exist in the scale effect test of rock joints, due to the geometry complexity of the rock joint surface and sampling bias. A quantitative description of roughness, suitable for incorporation into theoretically based models of joint behavior, should be proposed to express and deal with the incomplete and indeterminate problems of the JRC and shear strength of rough joints. This paper developed the two neutrosophic functions for describing the scale effect of JRC values and the shear strength of rock joints. They may provide a theoretical analysis method for predicting the JRC values and the shear strength of rock joints in the specific district in Shaoxing city.By analyzing real test data for the four samples, the JRC values of the different-sized test profiles have greater discreteness in small samples than those in large samples. Then, the mean values of JRC are declined with the sample length increasing. The results of the example showed that all the JRC values fall in the interval values of the neutrosophic function of JRC. Hence, the proposed neutrosophic function of JRC demonstrated its effectiveness in estimating the JRC values. On the basis of the proposed neutrosophic function of the peak shear strength of rock joints, an example was provided to show the effectiveness of a theoretical analysis method for the shear strength of rock joints in the specific district in Shaoxing city.The proposed neutrosophic functions can express and deal with the incomplete and indeterminate problems of the test data, which provide a new way for the analyses of the JRC values of the different-sized test profiles and the peak shear strength of rock joints. In the future, we will further develop and investigate the neutrosophic functions of the JRC and shear strength of rock joints in different districts in Shaoxing city.



---

*Source: 4825709-2016-05-18.xml*