Loading [MathJax]/jax/element/mml/optable/Latin1Supplement.js
Research article Special Issues

Assessing compressibility characteristics of silty soils from CPTU: lessons learnt from the Treporti Test Site, Venetian Lagoon (Italy)

  • This paper explores the capability of cone penetration testing as an effective tool to estimate the compressibility characteristics of silts, silt mixtures and sand mixtures. The study uses part of the large amount of piezocone (CPTU) data and subsoil deformation measurements assembled over approximately 6 years at the Treporti Test Site (Venice, Italy), within an extensive research project aimed at thoroughly analyzing the stress-strain-time response of the complex assortment of sandy and silty sediments forming the Venetian lagoon subsoil. The combined analysis of piezocone data and subsoil deformations measured beneath a full-scale test embankment primarily showed that the existing interpretation approaches developed for 'standard' sands or clays - and thus based on the idealized assumption of fully drained or fully undrained testing conditions - generally result in invalid estimates of soil compressibility parameters. Therefore, site-specific correlations have been developed in order to correctly predict the compressibility of Venetian sediments, both primary and secondary, and consequently obtain the best fit of the measured settlements. By taking this very well documented case study as a base, the proposed paper aims at discussing some key issues on the use of piezocone tests for the geotechnical characterization of silts and other sedimentary soils with very scattered grain size distributions, thus falling in the so-called intermediate permeability range. The idea is basically to provide a critical appraisal of the CPTU-based approaches currently available for the mechanical characterization of natural soils, and explore potential advances in the interpretation methods, which might take account of possible partial drainage phenomena around the advancing cone during the test.

    Citation: Laura Tonni, Guido Gottard. Assessing compressibility characteristics of silty soils from CPTU: lessons learnt from the Treporti Test Site, Venetian Lagoon (Italy)[J]. AIMS Geosciences, 2019, 5(2): 117-144. doi: 10.3934/geosci.2019.2.117

    Related Papers:

    [1] Yiqin Bao, Qiang Zhao, Jie Sun, Wenbin Xu, Hongbing Lu . An edge cloud and Fibonacci-Diffie-Hellman encryption scheme for secure printer data transmission. Mathematical Biosciences and Engineering, 2024, 21(1): 96-115. doi: 10.3934/mbe.2024005
    [2] Takao Komatsu, Haotian Ying . The p-Frobenius and p-Sylvester numbers for Fibonacci and Lucas triplets. Mathematical Biosciences and Engineering, 2023, 20(2): 3455-3481. doi: 10.3934/mbe.2023162
    [3] Zhongxue Yang, Yiqin Bao, Yuan Liu, Qiang Zhao, Hao Zheng, Wenbin Xu . Lightweight blockchain fuzzy decision scheme through MQTT and Fibonacci for sustainable transport. Mathematical Biosciences and Engineering, 2022, 19(12): 11935-11956. doi: 10.3934/mbe.2022556
    [4] Yan Wang, Guichen Lu, Jiang Du . Calibration and prediction for the inexact SIR model. Mathematical Biosciences and Engineering, 2022, 19(3): 2800-2818. doi: 10.3934/mbe.2022128
    [5] Anh Quang Tran, Tien-Anh Nguyen, Van Tu Duong, Quang-Huy Tran, Duc Nghia Tran, Duc-Tan Tran . MRI Simulation-based evaluation of an efficient under-sampling approach. Mathematical Biosciences and Engineering, 2020, 17(4): 4048-4063. doi: 10.3934/mbe.2020224
    [6] Yong Tao, Fan Ren, Youdong Chen, Tianmiao Wang, Yu Zou, Chaoyong Chen, Shan Jiang . A method for robotic grasping based on improved Gaussian mixture model. Mathematical Biosciences and Engineering, 2020, 17(2): 1495-1510. doi: 10.3934/mbe.2020077
    [7] Zhipeng Zhao, Zhenyu Hu, Zhiyu Zhao, Xiaoyu Du, Tianfei Chen, Lijun Sun . Fault-tolerant Hamiltonian cycle strategy for fast node fault diagnosis based on PMC in data center networks. Mathematical Biosciences and Engineering, 2024, 21(2): 2121-2136. doi: 10.3934/mbe.2024093
    [8] Auwalu Saleh Mubarak, Zubaida Said Ameen, Fadi Al-Turjman . Effect of Gaussian filtered images on Mask RCNN in detection and segmentation of potholes in smart cities. Mathematical Biosciences and Engineering, 2023, 20(1): 283-295. doi: 10.3934/mbe.2023013
    [9] Sergio Bermudo, Robinson A. Higuita, Juan Rada . k-domination and total k-domination numbers in catacondensed hexagonal systems. Mathematical Biosciences and Engineering, 2022, 19(7): 7138-7155. doi: 10.3934/mbe.2022337
    [10] Lal Hussain, Wajid Aziz, Ishtiaq Rasool Khan, Monagi H. Alkinani, Jalal S. Alowibdi . Machine learning based congestive heart failure detection using feature importance ranking of multimodal features. Mathematical Biosciences and Engineering, 2021, 18(1): 69-91. doi: 10.3934/mbe.2021004
  • This paper explores the capability of cone penetration testing as an effective tool to estimate the compressibility characteristics of silts, silt mixtures and sand mixtures. The study uses part of the large amount of piezocone (CPTU) data and subsoil deformation measurements assembled over approximately 6 years at the Treporti Test Site (Venice, Italy), within an extensive research project aimed at thoroughly analyzing the stress-strain-time response of the complex assortment of sandy and silty sediments forming the Venetian lagoon subsoil. The combined analysis of piezocone data and subsoil deformations measured beneath a full-scale test embankment primarily showed that the existing interpretation approaches developed for 'standard' sands or clays - and thus based on the idealized assumption of fully drained or fully undrained testing conditions - generally result in invalid estimates of soil compressibility parameters. Therefore, site-specific correlations have been developed in order to correctly predict the compressibility of Venetian sediments, both primary and secondary, and consequently obtain the best fit of the measured settlements. By taking this very well documented case study as a base, the proposed paper aims at discussing some key issues on the use of piezocone tests for the geotechnical characterization of silts and other sedimentary soils with very scattered grain size distributions, thus falling in the so-called intermediate permeability range. The idea is basically to provide a critical appraisal of the CPTU-based approaches currently available for the mechanical characterization of natural soils, and explore potential advances in the interpretation methods, which might take account of possible partial drainage phenomena around the advancing cone during the test.


    The Fibonacci sequence is defined Fn=Fn1+Fn2 for n2 with the initial conditions F0=0,F1=1. The Fibonacci sequence and the golden ratio are used in many fields such as cryptology, coding theory and quantum physics in [1,2,3,4,5,6,7,8,9,10,11]. Horadam defined the Gaussian Fibonacci numbers in [12,13] and gave some general identities about the Gaussian Fibonacci numbers. Jordan generalized the Gaussian Fibonacci numbers with a similar definition in [14]. The Gaussian Fibonacci sequence is defined as GFn=GFn1+GFn2 with n>1 where GF0=i,GF1=1 in [14]. It can also be easily seen Fn that GFn=Fn+iFn1 where Fn is the nth Fibonacci number.

    Asci and Gurel generalized these studies and defined the k-order Gaussian Fibonacci numbers in [15] by the following recurrence relation

    GF(k)n=kj=1GF(k)nj,forn>0andk2

    with boundary conditions for 1kn0,

    GF(k)n={1i,k=1ni,k=2n0,otherwise.

    It can also be seen that GF(k)n=F(k)n+iF(k)n1 where F(k)n is the nth k-order Fibonacci number.

    In [16], Asci and Aydinyuz defined the k-order Gaussian Fibonacci polynomials and gave some important results. The k-order Gaussian Fibonacci polynomials {GFn(k)}n=0 are

    GF(k)n(x)=kj=1xkjGF(k)nj(x)

    for n>0 and k2 with the initial conditions for 1kn0,

    GF(k)n(x)={1ix,k=1ni,k=2n0,otherwise.

    It can be seen that

    GF(k)n(x)=F(k)n(x)+iF(k)n1(x)

    where F(k)n(x) is the nth k-order Fibonacci polynomial.

    In order to ensure information security in terms of data transfer over communication channels, a lot of work has been done on this subject and continues to be done. Therefore, coding/decoding algorithms play an important role to ensure information security. Especially, Fibonacci coding theory is one of the most preferred in this field. We can see examples of these in many studies. For example, in [17], Stakhov gave a new approach to a coding theory using the generalization of the Cassini formula for Fibonacci p-numbers and Qp matrices in 2006. In 2009, Basu and Prasad in [18] presented the generalized relations among the code elements for the Fibonacci coding theory. Also, Basu and Das introduced a new coding theory for Tribonacci matrices in [19] and in 2014 they defined the coding theory for Fibonacci n-step numbers by generalizing the coding theory for Tribonacci numbers in [20]. Esmaeili in [21] also described a Fibonacci polynomial-based coding method with error detection and correction.

    In this paper, we remember the coding theory for k-order Gaussian Fibonacci polynomials given in [16] by taking x=1 for k-order Gaussian Fibonacci numbers. We describe the coding theory for the k-order Gaussian Fibonacci numbers. We give illustrative examples. The interesting relation between the elements of the code-message matrix is analyzed and derived. The most important of this paper is that we deal with error detection and error correction methods and examples of error probability is given.

    When we take x=1 in the coding theory defined in [16], we obtain the k-order Gaussian Fibonacci coding method. First, we define the Qk, Rk and En(k) matrices, which plays an important role in this coding theory. Qk, Rk and En(k) are defined the k×k matrices as the following:

    Qk=[111111000001000000000010]k×k,Rk=[GF(k)k1GF(k)k2GF(k)k3GF(k)2GF(k)10GF(k)k2GF(k)k3GF(k)k4GF(k)100GF(k)k3GF(k)k4GF(k)k5000GF(k)2GF(k)1000GF(k)10000i0000i1i]k×k

    and

    E(k)n=[GF(k)n+k1GF(k)n+k2GF(k)n+1GF(k)nGF(k)n+k2GF(k)n+k3GF(k)nGF(k)n1GF(k)n+1GF(k)nGF(k)nk+3GF(k)nk+2GF(k)nGF(k)n1GF(k)nk+2GF(k)nk+1]k×k

    where GF(k)n is the nth k-order Gaussian Fibonacci number.

    Theorem 2.1. For n1, we get in [15] as follow:

    QnkRk=E(k)n. (2.1)

    Corollary 2.1. For k=2, we get

    QnR=[1110]n[1ii1i]
    =[GFn+1GFnGFnGFn1]

    where GFn is the nth usual Gaussian Fibonacci number in [15].

    Corollary 2.2. For k=3, we get

    Qn3R3=[111100010]n[1+i1010i0i1i]
    =[GTn+2GTn+1GTnGTn+1GTnGTn1GTnGTn1GTn2]

    where GTn is the nth Gaussian Tribonacci number in [23].

    In this section, we redefine k-order Gaussian Fibonacci coding theory using k-order Gaussian Fibonacci numbers and Qk, Rk and En(k) play very important role in the construction of k-order Gaussian Fibonacci coding theory.

    We now obtain the matrix En(k) using Qk, Rk matrices for the k=2,k=3 values and examine the inverse.

    • For k=2, introducing the square matrix Q2 of order 2 as:

    Q2=[1110]

    and the square matrix R2 of order 2 as:

    R2=[1ii1i]

    for n=1, we can use (2.1)

    E(2)1=Q2.R2
    =[1+i11i]
    =[GF(2)2GF(2)1GF(2)1GF(2)0]

    such that

    detE(2)1=det(Q2.R2)
    =(1).(2i)
    =2+i.

    The inverse of E(2)1 is as:

    (E(2)1)1=[1525i25+15i25+15i1535i]
    =1detE(2)1[GF0GF1GF1GF2]
    =12+i[GF0GF1GF1GF2]

    such that

    det(E(2)1)1=12+i
    =2515i.

    Also, by (2.1) for n=2, we can get E(2)2 as follows:

    E(2)2=Q22.R2
    =[2+i1+i1+i1]
    =[GF(2)3GF(2)2GF(2)2GF(2)1]

    such that

    detE(2)2=det(Q22.R2)
    =(1)2.(2i)
    =2i.

    The inverse of E(2)2 is as:

    (E(2)2)1=[25+15i1535i1535i35+45i]
    =1detE(2)2[GF(2)1GF(2)2GF(2)2GF(2)3]
    =12i[GF(2)1GF(2)2GF(2)2GF(2)3]

    such that

    det(E(2)2)1=12i
    =25+15i.

    Theorem 2.1.1. E(2)n=[GF(2)n+1GF(2)nGF(2)nGF(2)n1] where E(2)1=[1+i11i].

    Theorem 2.1.2. (E(2)n)1=1detE(2)n[GF(2)n1GF(2)nGF(2)nGF(2)n+1]

    =1(detQ2)n.detR2[GF(2)n1GF(2)nGF(2)nGF(2)n+1].

    • For k=3, introducing the square matrix Q3 of order 3 as:

    Q3=[111100010]

    and the square matrix R3 of order 3 as:

    R3=[1+i1010i0i1i]

    for n=1, we can use (2.1)

    E(3)1=Q3.R3
    =[2+i1+i11+i1010i]
    =[GF(3)3GF(3)2GF(3)1GF(3)2GF(3)1GF(3)0GF(3)1GF(3)0GF(3)1]

    such that

    detE(3)1=det(Q3.R3)
    =1.2i
    =2i.

    The inverse of E(3)1 is as:

    (E(3)1)1=[121212i12i1212i1+i1212i12i1212i12i]
    =1detE(3)1[GF(3)1GF1(3)(GF(3)0)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)0GF2(3)(GF(3)1)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)3GF1(3)(GF(3)1)2GF(3)1GF2(3)GF(3)0GF3(3)(GF(3)1)2GF(3)0GF2(3)GF(3)1GF2(3)GF(3)0GF3(3)GF(3)1GF3(3)(GF(3)2)2]
    =1detR3[GF(3)1GF1(3)(GF(3)0)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)0GF2(3)(GF(3)1)2GF(3)0GF1(3)GF(3)1GF2(3)GF(3)3GF1(3)(GF(3)1)2GF(3)1GF2(3)GF(3)0GF3(3)(GF(3)1)2GF(3)0GF2(3)GF(3)1GF2(3)GF(3)0GF3(3)GF(3)1GF3(3)(GF(3)2)2]

    such that

    det(E(3)1)1=12i
    =i2.

    Also, by (2.1) for n=2, we can get E(3)2 as follows:

    E(3)2=Q23.R3
    =[4+2i2+i1+i2+i1+i11+i10]
    =[GF(3)4GF(3)3GF(3)2GF(3)3GF(3)2GF(3)1GF(3)2GF(3)1GF(3)0]

    such that

    detE(3)2=det(Q23.R3)
    =12.2i
    =2i.

    The inverse of E(3)2 is as:

    (E(3)2)1=[12i1212i12i1212i112+32i12i12+32i1+12i]
    =1detE(3)2[GF(3)0GF(3)2(GF(3)1)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)1GF(3)3(GF(3)2)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)0GF(3)4(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)1GF(3)3(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)2GF(3)4(GF(3)3)2]
    =1detR3[GF(3)0GF(3)2(GF(3)1)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)1GF(3)3(GF(3)2)2GF(3)1GF(3)2GF(3)0GF(3)3GF(3)0GF(3)4(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)1GF(3)3(GF(3)2)2GF(3)2GF(3)3GF(3)1GF(3)4GF(3)2GF(3)4(GF(3)3)2].

    Theorem 2.1.3. E(3)n=[GF(3)n+2GF(3)n+1GF(3)nGF(3)n+1GF(3)nGF(3)n1GF(3)nGF(3)n1GF(3)n2] where E(3)1=[2+i1+i11+i1010i].

    Theorem 2.1.4.

    (E(3)n)1=1detR3[GF(3)n2GF(3)n(GF(3)n1)2GF(3)n1GF(3)nGF(3)n2GF(3)n+1GF(3)n1GF(3)n+1(GF(3)n)2GF(3)n1GF(3)nGF(3)n2GF(3)n+1GF(3)n2GF(3)n+2(GF(3)n)2GF(3)nGF(3)n+1GF(3)n1GF(3)n+2GF(3)n1GF(3)n+1(GF(3)n)2GF(3)nGF(3)n+1GF(3)n1GF(3)n+2GF(3)nGF(3)n+2(GF(3)n+1)2].

    For arbitrary k-positive integers, the square matrix E(k)n of order k and inverses can be found similarly.

    In this section, we describe a new k-order Gaussian Fibonacci coding theory. We put our message in a matrix of M and let us represent the initial message in the form of the square matrix M of order k. We take the k-order Gaussian Fibonacci matrix E(k)n as a coding matrix and its inverse matrix (E(k)n)1 as a decoding matrix for an arbitrary positive integer k. The transformation M×E(k)n=C is called k-order Gaussian Fibonacci coding and we name the transformation C×(E(k)n)1=M as k-order Gaussian Fibonacci decoding. We define C as a code matrix.

    The given example is solved using the alphabet table below.

    Using the arbitrary value of s0, we write the following alphabet according to mod30. We can extend the characters in the table according to our wishes. We begin the “s” for the first character in Table 1.

    Table 1.  Alphabet table.
    A B C D E F G H I J
    s s+1 s+2 s+3 s+4 s+5 s+6 s+7 s+8 s+9
    K L M N O P Q R S T
    s+10 s+11 s+12 s+13 s+14 s+15 s+16 s+17 s+18 s+19
    U V W X Y Z 0 ! ? .
    s+20 s+21 s+22 s+23 s+24 s+25 s+26 s+27 s+28 s+29

     | Show Table
    DownLoad: CSV

    Example 2.3.1. Let us consider the message matrix for the following message text:

    {\rm{“CODE”}} (2.2)

    Step 1: Let’s create the message matrix using the message text:

    M = {\left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right]_{2 \times 2}}

    Step 2: Let’s write the message matrix M according to the alphabet table for the arbitrary value “s” we choose. For s = 2 ;

    M = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]

    Step 3: For k = 2 , n = 2 , we use (2.1);

    E_2^{(2)} = Q_2^2.{R_2}
    = \left[ {\begin{array}{*{20}{c}} {2 + i}&{1 + i} \\ {1 + i}&1 \end{array}} \right]

    Step 4: The code message is:

    C = M \times E_2^{(2)}
    = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {2 + i}&{1 + i} \\ {1 + i}&1 \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right]

    Step 5: The decode message is:

    M = C \times {\left( {E_2^{(2)}} \right)^{ - 1}}
    = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right]\left[ {\begin{array}{*{20}{c}} {\frac{2}{5} + \frac{1}{5}i}&{ - \frac{1}{5} - \frac{3}{5}i} \\ { - \frac{1}{5} - \frac{3}{5}i}&{\frac{3}{5} + \frac{4}{5}i} \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right] .

    Example 2.3.2. Let us consider the message matrix for the following message text:

    {\rm{"PUBLIC \;KEY"}}

    Step 1: Let’s create the message matrix using the message text:

    M = {\left[ {\begin{array}{*{20}{c}} P&U&B \\ L&I&C \\ K&E&Y \end{array}} \right]_{3 \times 3}}

    Step 2: Let’s write the message matrix M according to the alphabet table for the arbitrary value “s” we choose. For s = 5 :

    M = \left[ {\begin{array}{*{20}{c}} {20}&{25}&6 \\ {16}&{13}&7 \\ {15}&9&{29} \end{array}} \right]

    Step 3: For k = 3 , n = 6 , we use (2.1):

    E_6^{(3)} = Q_3^6.{R_3}
    = \left[ {\begin{array}{*{20}{c}} {44 + 24i}&{24 + 13i}&{13 + 7i} \\ {24 + 13i}&{13 + 7i}&{7 + 4i} \\ {13 + 7i}&{7 + 4i}&{4 + 2i} \end{array}} \right]

    Step 4: The code message is:

    C = M \times E_6^{(3)}
    = \left[ {\begin{array}{*{20}{c}} {1558 + 847i}&{847 + 459i}&{459 + 252i} \\ {1107 + 602i}&{602 + 327i}&{327 + 178i} \\ {1253 + 680i}&{680 + 374i}&{374 + 199i} \end{array}} \right]

    Step 5: The decode message is:

    M = C \times {\left( {E_3^{(6)}} \right)^{ - 1}}
    = \left[ {\begin{array}{*{20}{c}} P&U&B \\ L&I&C \\ K&E&Y \end{array}} \right]
    = \left[ {\begin{array}{*{20}{c}} {20}&{25}&6 \\ {16}&{13}&7 \\ {15}&9&{29} \end{array}} \right] .

    In this section, we consider the k-order Gaussian Fibonacci coding/decoding method for k = 2 . We have an interesting relation among the elements of a code matrix C that has a crucial role in the error-correction process, outlined as follows:

    C = M \times E_n^{(2)} = \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right] \times \left[ {\begin{array}{*{20}{c}} {GF_{n + 1}^{(2)}}&{GF_n^{(2)}} \\ {GF_n^{(2)}}&{GF_{n - 1}^{(2)}} \end{array}} \right] (3.1)
    = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] (3.2)

    and

    M = C \times {\left( {E_n^{(2)}} \right)^{ - 1}} (3.3)
    = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] \times \left( {\frac{1}{{{{\left( {\det \left( {{Q_2}} \right)} \right)}^n}\det {R_2}}}\left[ {\begin{array}{*{20}{c}} {GF_{n - 1}^{(2)}}&{ - GF_n^{(2)}} \\ { - GF_n^{(2)}}&{GF_{n + 1}^{(2)}} \end{array}} \right]} \right) (3.4)
    = \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right]

    For the case for an even integer n = 2m , we obtain the following equation

    \left[ {\begin{array}{*{20}{c}} {{m_1}}&{{m_2}} \\ {{m_3}}&{{m_4}} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] \times \left( {\frac{1}{{2 - i}}\left[ {\begin{array}{*{20}{c}} {GF_{n - 1}^{(2)}}&{ - GF_n^{(2)}} \\ { - GF_n^{(2)}}&{GF_{n + 1}^{(2)}} \end{array}} \right]} \right) (3.5)

    since the \det \left({{Q_2}} \right) = 1 and \det \left({{R_2}} \right) = 2 - i .

    It follows from (3.5) that the elements of the matrix M can be obtained according to the following formulas:

    {m_1} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_1} - GF_n^{(2)}{c_2}} \right) (3.6)
    {m_2} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_1} + GF_{n + 1}^{(2)}{c_2}} \right) (3.7)
    {m_3} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_3} - GF_n^{(2)}{c_4}} \right) (3.8)
    {m_4} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_3} + GF_{n + 1}^{(2)}{c_4}} \right) (3.9)

    Since the s \geqslant 0 , the elements of the matrix M are

    {m_1} \geqslant 0 , {m_2} \geqslant 0 , {m_3} \geqslant 0 \;{\rm{ and}}\; {m_4} \geqslant 0 (3.10)

    Because of the condition (3.10) we can write the equalities (3.6)–(3.9) as:

    {m_1} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_1} - GF_n^{(2)}{c_2}} \right) \geqslant 0 (3.11)
    {m_2} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_1} + GF_{n + 1}^{(2)}{c_2}} \right) \geqslant 0 (3.12)
    {m_3} = \frac{1}{{2 - i}}\left( {GF_{n - 1}^{(2)}{c_3} - GF_n^{(2)}{c_4}} \right) \geqslant 0 (3.13)
    {m_4} = \frac{1}{{2 - i}}\left( { - GF_n^{(2)}{c_3} + GF_{n + 1}^{(2)}{c_4}} \right) \geqslant 0 (3.14)

    From the Eqs (3.11) and (3.12), we can get

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}}{c_2} \leqslant {c_1} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}}{c_2} (3.15)

    or

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}} \leqslant \frac{{{c_1}}}{{{c_2}}} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}} (3.16)

    Similarly, from (3.13) and (3.14), we can obtain:

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}}{c_4} \leqslant {c_3} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}}{c_4} (3.17)

    or

    \frac{{GF_n^{(2)}}}{{GF_{n - 1}^{(2)}}} \leqslant \frac{{{c_3}}}{{{c_4}}} \leqslant \frac{{GF_{n + 1}^{(2)}}}{{GF_n^{(2)}}} (3.18)

    Since two consecutive Gaussian Fibonacci numbers approach the golden ratio, we also obtain the following equations in (3.16) and (3.18) that connect the elements of the code matrix in (3.1):

    {c_1} \approx \tau {c_2} (3.19)
    {c_3} \approx \tau {c_4} (3.20)

    where \tau = \frac{{1 + \sqrt 5 }}{2} is the “golden ratio”.

    Similarly, for an odd integer n = 2m + 1 , we can write the same approximate equalities (3.19) and (3.20) that connect the elements {c_1}, {c_2}, {c_3} and {c_4} of the code matrix in (3.1).

    Thus, we have generated some important identities that connect the elements of the code matrix in (3.1) for the case k = 2 .

    If we consider the coding/decoding method for k = 3 , we can obtain interesting identities as in k = 2 . Since two consecutive Gaussian Tribonacci numbers approach \Phi where \Phi = 1.8393. However, for general case of k we can find mathematical identities that connect the code matrix elements similar to the (3.19) and (3.20).

    Example 3.1. Suppose k = 2 . Thus, we have \tau = \frac{{1 + \sqrt 5 }}{2} \cong 1.618 . Assume that the following message-matrix is to be transmitted in (2.2):

    For s = 2 , we can get

    M = \left[ {\begin{array}{*{20}{c}} C&O \\ D&E \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} 4&{16} \\ 5&6 \end{array}} \right]

    • If n = 2 then

    C = \left[ {\begin{array}{*{20}{c}} {24 + 20i}&{20 + 4i} \\ {16 + 11i}&{11 + 5i} \end{array}} \right] = \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right]

    In this case, we have the following numbers rounded off to their first five digits:

    \frac{{{c_1}}}{{{c_2}}} = 1.3462 + 0.7307i , \frac{{{c_3}}}{{{c_4}}} = 1.5822 + 0.2808i

    • For n = 5 , we have

    C = \left[ {\begin{array}{*{20}{c}} {112 + 68i}&{68 + 44i} \\ {70 + 43i}&{43 + 27i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.6171 - 0.0463i , \frac{{{c_3}}}{{{c_4}}} = 1.6179 - 0.0159i

    • For n = 10 , we have

    C = \left[ {\begin{array}{*{20}{c}} {1236 + 764i}&{764 + 472i} \\ {775 + 479i}&{479 + 296i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.618 + 0.0003i , \frac{{{c_3}}}{{{c_4}}} = 1.618 + 0.0001i

    • For n = 15 , we have

    C = \left[ {\begin{array}{*{20}{c}} {13708 + 8472i}&{8472 + 5236i} \\ {8595 + 5312i}&{5312 + 3283i} \end{array}} \right]

    In this case

    \frac{{{c_1}}}{{{c_2}}} = 1.618 , \frac{{{c_3}}}{{{c_4}}} = 1.618

    These show that for n = 15 the relation (3.19) and (3.20) holds very well.

    Therefore, we select the value of n large enough, as seen in Examples 2.3.1 we get a good approximation of \tau . If we apply similar operation in Example 2.3.2, when n is chosen large enough we get a good approximation of \Phi where \Phi = 1.8393 and \Phi is the ratio between two Gaussian Tribonacci numbers.

    The k-order Gaussian Fibonacci coding/decoding method considered above provides an interesting possibility to detect and correct “errors” in code message C . Error detection and correction for k = 2 depends on the states of the determinant of the matrix C and the connections with the code matrix elements at (3.19) and (3.20). These mathematical relations given by (3.19), (3.20) and the determinant of the code matrix C play a role in checking relations of the k-order Gaussian Fibonacci coding/decoding method. The error-correction algorithm for Fibonacci coding has been explained in [10] and [13]. This defined algorithm is valid and applicable to the k-order Gaussian Fibonacci numbers we have given in this study.

    First, let’s calculate the determinant of our message matrix M and send it to the communication channel right after the code matrix elements. Here, the determinant of M will be the control element for the C code matrix received from the communication channel. After taking the determinant of the code matrix C and its control element M , we can calculate the determinant of the code matrix C and compare the determinant of M with the determinant of C with respect to the control relation. If the operations we have done after this check are exactly appropriate, we can conclude that the elements of the code matrix C are transmitted on the communication channel without errors. If the opposite happens, we can say that there are errors in the elements of the code matrix C or the determinant of the control element M .

    Suppose the elements of the C code matrix errors. This matrix can have a one-fold error, double-fold error or {k^2} - fold error. To explain how to correct these errors, we consider the 2 \times 2 type matrix obtained for k = 2 . These cases are considered.

    Our first hypothesis is that we have the case of a “single” error in the code matrix C received from the communication channel. It is clearly seen that there are four types of single-fold errors in the C code matrix:

    (a) \left[ {\begin{array}{*{20}{c}} x&{{c_2}} \\ {{c_3}}&{{c_4}} \end{array}} \right] (b) \left[ {\begin{array}{*{20}{c}} {{c_1}}&y \\ {{c_3}}&{{c_4}} \end{array}} \right] (c) \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ z&{{c_4}} \end{array}} \right] (d) \left[ {\begin{array}{*{20}{c}} {{c_1}}&{{c_2}} \\ {{c_3}}&t \end{array}} \right] (4.1)

    where x, y, z, t are possible “destroyed” elements.

    Since it is C = M \times E_n^{(2)} , if we take the determinant of both sides, we get the following equations:

    x{c_4} - {c_2}{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a\; possible \;“single \;error” \;is \;in \;the \;element}}\; {c_1} ) (4.2)
    {c_1}{c_4} - y{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a\; possible \;“single \;error” \;is \;in \;the\; element}}\; {c_2} ) (4.3)
    {c_1}{c_4} - {c_2}z = \det (M){( - 1)^n}(2 - i) ({\rm{a \;possible \;“single \;error”\; is \;in \;the \;element}}\; {c_3} ) (4.4)
    {c_1}t - {c_2}{c_3} = \det (M){( - 1)^n}(2 - i) ({\rm{a \;possible \;“single \;error” \;is \;in \;the \;element}}\; {c_4} ) (4.5)

    or equivalently

    x = \frac{{\det (M){{( - 1)}^n}(2 - i) + {c_2}{c_3}}}{{{c_4}}} (a) (4.6)
    y = \frac{{ - \det (M){{( - 1)}^n}(2 - i) + {c_1}{c_4}}}{{{c_3}}} (b) (4.7)
    z = \frac{{ - \det (M){{( - 1)}^n}(2 - i) + {c_1}{c_4}}}{{{c_2}}} (c) (4.8)
    t = \frac{{\det (M){{( - 1)}^n}(2 - i) + {c_2}{c_3}}}{{{c_1}}} (d) (4.9)

    The formulas (4.6)–(4.9) give four possible states of a “single error”. However, we need to choose the correct case from among x, y, z, t integer solutions. We have to choice such solutions, which satisfies to the additional “checking relations” (3.19) and (3.20). If calculations by formulas (4.6)–(4.9) do not give an integer result, we can conclude that our “single error” hypothesis is incorrect or that the determinant of control element M is “error”. In the second case, we can use approximate Eqs (3.19) and (3.20) to check the correctness of the C code matrix.

    Similarly, we can easily check “double errors” of the C code matrix. Suppose that the code matrix C has elements x and y error as shown below:

    \left[ {\begin{array}{*{20}{c}} x&y \\ {{c_3}}&{{c_4}} \end{array}} \right] (4.10)

    Then,

    x{c_4} - y{c_3} = \det (M){( - 1)^n}(2 - i) (4.11)

    However, according to the (3.19) there is the following relation between x and y

    x \approx \tau y (4.12)

    Again, only integer solutions are acceptable. If no integer solution is obtained, then two-fold errors have not occurred.

    Hence, we can show that with this type of approach, there is the possibility to correct all possible triple-fold errors in the C code matrix.

    \left[ {\begin{array}{*{20}{c}} x&y \\ z&{{c_4}} \end{array}} \right] .

    Therefore, our method of error correction depends on the verification of different hypotheses regarding errors in the code matrix using the determinant of the control relationship C matrix, (3.19) and (3.20), and the elements of the code matrix as integers. If all our solutions do not reach integer solutions, it means that the determinant of the control item M is wrong or there is a four-fold error in the C code matrix, and we should reject it. The code matrix C is defective and not correctable.

    As a result, there are 15 error conditions in the C code matrix. According to the method given in [17], it means that the correctable probability of the method is equal since 14 cases between them can be corrected and

    {S_{cor}} = \frac{{14}}{{15}} = 0.9333 = 93.33\%

    If we generalize this equation as in [20], since only {k^2} - fold errors cannot be corrected with this method, the error correction capacity of the method is

    \frac{{{2^{{k^2}}} - 1}}{{{2^{{k^2}}}}}

    where k is the order of the message-matrix. So far sufficiently large values of k , the probability of decoding error is almost zero.

    Consequently, for a sufficiently large value of k , in the case of n = k , the correct ability of this method is

    \frac{{{2^{{k^2}}} - 2}}{{{2^{{k^2}}} - 1}} .

    Therefore, for large value of k , the correct possibility of the method is

    \frac{{{2^{{k^2}}} - 2}}{{{2^{{k^2}}} - 1}} \approx 1 = 100\% .

    In this paper, we obtained the coding theory for k-order Gaussian Fibonacci polynomials given in [16] by taking x = 1 for k-order Gaussian Fibonacci numbers and gave illustrative examples. This coding method differs from classical algebraic coding methods in this respect. We can show this difference with the following features:

    1)      Since the n, k and s \geqslant 0 values are arbitrarily chosen in the k-order Gaussian Fibonacci coding theory, it is difficult to estimate the information transmission between the two channels by a third channel, which increases the reliability of the information transmission. In addition, this method depends on matrix multiplication and can be performed quickly and easily by today’s computers.

    2)      With k-order Gaussian Fibonacci coding theory, we can encrypt and send messages of the desired length by enlarging the k value sufficiently.

    3)      The main practical feature of this method is that matrix elements are error detection and correction objects.

    4)      The Gaussian Fibonacci coding method, which is the simplest coding method obtained for the k = 2 value, has been redefined by handling the 2 \times 2 type matrices with “single-fold, double-fold and triple-fold errors”.

    5)      In the simplest case, for k = 2 , the correct capability of the method is essentially equal to 93.33%, exceeding all well-known correction codes.

    6)      This error correction method given for k = 2 is generalized and the correct ability of the errors of this method increases as k increases, and the correct ability approaches 100% for a sufficiently large value of k .

    7)      This article is just a brief outline of a new coding theory based on the k-order Gaussian Fibonacci matrices of the articles in [17,18,20,21] and [22].

    The authors declare there is no conflict of interest.



    [1] Rowe PW (1973) Soil mechanic aspect of the cores of the deep borehole VE 1 in Venice. A critical analysis and recommended future investigations. Venice (Italy): National Research Council.
    [2] Ricceri G, Butterfield R (1974) An analysis of compressibility data from a deep borehole in Venice. Géotechnique 24: 175–191. doi: 10.1680/geot.1974.24.2.175
    [3] Ricceri G, Simonini P, Cola S (2002) Applicability of piezocone and dilatometer to characterize the soils of the Venice Lagoon.   Geotech Geol Eng 20: 89–121.
    [4] Jamiolkowski M, Ricceri G, Simonini P (2009) Safeguarding Venice from high tides: site characterization and geotechnical problems. In: Proc. 17th ICSMGE, The Netherlands: IOS Press, 3209–3227.
    [5] Tonni L, Gottardi G (2011) Analysis and interpretation of piezocone data on the silty soils of the Venetian lagoon (Treporti test site). Can Geotech J 48: 616–633. doi: 10.1139/t10-085
    [6] Cola S, Simonini P (2002) Mechanical behavior of silty soils of the Venice lagoon as a function of their grading characteristics. Can Geotech J 39: 879–893. doi: 10.1139/t02-037
    [7] Tonni L, García Martínez MF, Simonini P, et al. (2016) Piezocone-based prediction of secondary compression settlements of coastal defence structures on natural silt mixtures. Ocean Eng 116: 101–116. doi: 10.1016/j.oceaneng.2016.02.015
    [8] Simonini P (2004) Characterization of the Venice lagoon silts from in-situ tests and the performance of a test embankment. In: Proc. 2nd Int. Conf. on Site Characterization, ISC'2, Rotterdam (The Netherlands): Millpress, 187–207.
    [9] Gottardi G, Tonni L (2004) Use of piezocone tests to characterize the silty soils of the Venetian lagoon (Treporti Test Site). In: Proc 2nd Int Conf on Site Characterization, ISC'2, Rotterdam (The Netherlands): Millpress, 1643–1650.
    [10] Marchetti S, Monaco P, Calabrese M, et al. (2004) DMT-predicted vs measured settlements under a full-scale instrumented embankment at Treporti (Venice, Italy). In: Proc 2nd Int Conf on Site Characterization, ISC'2, Rotterdam (The Netherlands): Millpress, 1511–1518.
    [11] Monaco P, Amoroso S, Marchetti S, et al. (2014) Overconsolidation and stiffness of Venice lagoon sands and silts from SDMT and CPTU. J Geotech Geoenviron Eng 140: 215–227. doi: 10.1061/(ASCE)GT.1943-5606.0000965
    [12] McGillivray A, Mayne PW (2004) Seismic piezocone and seismic flat dilatometer tests at Treporti. In: Proc 2nd Int Conf on Site Characterization, ISC'2, Rotterdam (The Netherlands): Millpress, 1695–1700.
    [13] Tonni L, Simonini P (2013) Shear wave velocity as function of cone penetration test measurements in sand and silt mixtures. Eng Geol 163: 55–67. doi: 10.1016/j.enggeo.2013.06.005
    [14] Tonni L, Gottardi G (2009) Partial drainage effects in the interpretation of piezocone tests in Venetian silty soils. In: Proc 17th Int Conf on Soil Mechanics and Geotechnical Eng. (ICSMGE),The Netherlands: IOS Press, 1004–1007.
    [15] Tonni L, Gottardi G (2010) Interpretation of piezocone tests in Venetian silty soils and the issue of partial drainage. In: Deep Foundations and Geotechnical In Situ Testing, ASCE Geotechnical Special Publication (205 GSP), Reston (VA): ASCE, 367–374.
    [16] Schnaid F, Lehane BM, Fahey M (2004) Characterisation of unusual geomaterials. In: Proc 2ndInt Conf on Geotechnical and Geophysical Site Characterization (ISC-2), Rotterdam: Millpress, 49–73
    [17] García Martínez MF, Tonni L, Gottardi G (2016) On the interpretation of piezocone tests in natural silt and sand mixtures. In: Proc 1st IMEKO TC-4 International Workshop on Metrology for Geotechnics, Metrogeotechnics 2016, IMEKO-International Measurement Confederation, 63–68.
    [18] García Martínez MF, Tonni L, Gottardi G, et al. (2016) Influence of penetration rate on CPTU measurements in saturated silty soils. In: Proc 5th Int Conf on Geotechnical and Geophysical Site Characterisation (ISC'5), Sidney: Australian Geomechanics Society, 473–478.
    [19] Schneider JA, Randolph MF, Mayne PW, et al. (2008) Analysis of factors influencing soil classification using normalized piezocone tip resistance and pore pressure parameters. J Geotech Geoenviron Eng 134: 1569–1586. doi: 10.1061/(ASCE)1090-0241(2008)134:11(1569)
    [20] Robertson PK (2009) Interpretation of cone penetration tests-a unified approach. Can Geotech J 46: 1337–1355. doi: 10.1139/T09-065
    [21] Robertson PK (2016) Cone penetration test (CPT)-based soil behaviour type (SBT) classification system-an update. Can Geotech J 53: 1910–1927. doi: 10.1139/cgj-2016-0044
    [22] Lunne T, Robertson PK, Powell JJM (1997) Cone penetration testing in geotechnical practice. London: Blackie Academic & Professional.
    [23] Dejong JT, Randolph MF (2012) Influence of partial consolidation during cone penetration on estimated soil behaviour type and pore pressure dissipation measurements. J Geotech Geoenviron Eng 138: 777–788. doi: 10.1061/(ASCE)GT.1943-5606.0000646
    [24] Teh CI, Houlsby GT (1991) An analytical study of the cone penetration test in clay. Géotechnique 41: 17–34. doi: 10.1680/geot.1991.41.1.17
    [25] Randolph MF (2004) Characterisation of soft sediments for offshore applications. In: Proc 2ndInt Conf on Site Characterization, ISC'2, Rotterdam (The Netherlands): Millpress, 209–232.
    [26] Leroueil S, Kabbaj M, Tavenas F, et al. (1985) Stress-strain-strain rate relation for the compressibility of sensitive natural clays. Géotechnique 35: 159–180. doi: 10.1680/geot.1985.35.2.159
    [27] Senneset K, Janbu N, Svanø G (1982) Strength and deformation parameters from cone penetration tests. In: Proc 2nd European Symp on Penetration Testing (ESOPT-II), Rotterdam: Balkema, 863–870.
    [28] Lunne T, Christophersen HP (1983) Interpretation of cone penetrometer data for offshore sands. In: Proc of the Offshore Technology Conference, 4464.
    [29] Senneset K, Sandven R, Lunne T et al. (1988) Piezocone tests in silty soils. In: Proc 1st Int Symp on Penetration Testing (ISOPT-1), Rotterdam: Balkema, 955–966.
    [30] Kulhawy FH, Mayne PH (1990) Manual on estimating soil properties for foundation design. Electric Power Research Institute (EPRI), Palo Alto, CA.
    [31] Simonini P, Tonni L (2014) Primary and secondary compression characteristics of silts from CPTU: lessons learnt from the Venetian Lagoon sediments. In: Proc 3rd Int Symp on Cone Penetration Testing (CPT'14), 745–752.
    [32] Augustesen A, Liingaard M, Lade PV (2004) Evaluation of time-dependent behaviour of soils. Int J of Geomech 4: 137–156. doi: 10.1061/(ASCE)1532-3641(2004)4:3(137)
    [33] Leroueil S (2006) The isotache approach.Where are we 50 years after its development by Professor Šuklje? In: Proc XIII Danube-European Conference on Geotechnical Engineering,Ljubiana, Slovenia, 1: 55–88
    [34] Mesri G, Godlewski PM (1977) Time and stress compressibility interrelationship. J Geotech Eng Div, ASCE 103 (GT5): 417–430.
    [35] Tonni L, Simonini P (2013) Evaluation of secondary compression of sands and silts from CPTU. Geomech Geoeng 8: 141–154. doi: 10.1080/17486025.2012.726748
  • This article has been cited by:

    1. Öznur ÖZTUNÇ KAYMAK, Coding theory for h(x)-Fibonacci polynomials, 2024, 26, 1301-7985, 226, 10.25092/baunfbed.1347379
    2. Bahar Demirtürk, Combinatorial Analysis of k-Oresme and k-Oresme–Lucas Sequences, 2025, 17, 2073-8994, 697, 10.3390/sym17050697
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6466) PDF downloads(1379) Cited by(9)

Figures and Tables

Figures(20)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog