Research article

Geotechnical characterization of index and deformation properties of Stockholm clays

  • This paper describes the geotechnical characteristics of index and deformation properties of Stockholm clays. These clays exhibit large variations both geographically and with depth. The clays range from highly organic clays with high plasticity to silty clays with low plasticity. First, the geological conditions of the clays are outlined to qualitatively explain the typical soil stratigraphy encountered in the Stockholm area. Second, two large generic databases are presented, containing 3,500 and 1,600 data points, respectively. The data originates from routine testing and constant rate of strain oedometer tests conducted in commercial projects. The data is analyzed and compared with results from high quality block sampling. It is seen that a common feature of the clays are low undrained shear strengths, and consequently low yield stresses and oedometer moduli. Further, several deformation properties such as preconsolidation pressure and several oedometer moduli is shown to depend on the soils natural water content and its plasticity. Differences in sample quality is shown to highly affect some properties, highlighting the importance of quality sampling and handling of samples. Criteria for sample quality for this type of clay is proposed based on the oedometer moduli before and after the preconsolidation pressure. The paper can hopefully work as a useful reference to engineers working on similar soils worldwide.

    Citation: Solve Hov, David Gaharia. Geotechnical characterization of index and deformation properties of Stockholm clays[J]. AIMS Geosciences, 2023, 9(2): 258-284. doi: 10.3934/geosci.2023015

    Related Papers:

    [1] Rahmad Doni Widodo, Rusiyanto, Wahyudi, Melisa Kartika Sari, Deni Fajar Fitriyana, Januar Parlaungan Siregar, Tezara Cionita, Natalino Fonseca Da Silva Guterres, Mateus De Sousa Da Silva, Jamiluddin Jaafar . Effect of compression molding temperature on the characterization of asbestos-free composite friction materials for railway applications. AIMS Materials Science, 2023, 10(6): 1105-1120. doi: 10.3934/matersci.2023059
    [2] Bauyrzhan Rakhadilov, Lyaila Bayatanova, Sherzod Kurbanbekov, Ravil Sulyubayev, Nurdaulet Shektibayev, Nurbol Berdimuratov . Investigation on the effect of technological parameters of electrolyte-plasma cementation method on phase structure and mechanical properties of structural steel 20X. AIMS Materials Science, 2023, 10(5): 934-947. doi: 10.3934/matersci.2023050
    [3] Leydi J. Cardenas F., Josep Ma. Chimenos, Luis C. Moreno A., Elaine C. Paris, Miryam R. Joya . Enhancing Co3O4 nanoparticles: Investigating the impact of nickel doping and high-temperature annealing on NiCo2O4/CoO heterostructures. AIMS Materials Science, 2023, 10(6): 1090-1104. doi: 10.3934/matersci.2023058
    [4] Ayu Azera Izad, Rosimah Nulit, Che Azurahanim Che Abdullah, A'fifah Abdul Razak, Teh Huey Fang, Mohd Hafiz Ibrahim . Growth, leaf gas exchange and biochemical changes of oil palm (Elaeis guineensis Jacq.) seedlings as affected by iron oxide nanoparticles. AIMS Materials Science, 2019, 6(6): 960-984. doi: 10.3934/matersci.2019.6.960
    [5] Larysa Khomenkova, Mykola Baran, Jedrzej Jedrzejewski, Caroline Bonafos, Vincent Paillard, Yevgen Venger, Isaac Balberg, Nadiia Korsunska . Silicon nanocrystals embedded in oxide films grown by magnetron sputtering. AIMS Materials Science, 2016, 3(2): 538-561. doi: 10.3934/matersci.2016.2.538
    [6] Khang Duy Vu Nguyen, Khoa Dang Nguyen Vo . Magnetite nanoparticles-TiO2 nanoparticles-graphene oxide nanocomposite: Synthesis, characterization and photocatalytic degradation for Rhodamine-B dye. AIMS Materials Science, 2020, 7(3): 288-301. doi: 10.3934/matersci.2020.3.288
    [7] Laís Weber Aguiar, Cleiser Thiago Pereira da Silva, Hugo Henrique Carline de Lima, Murilo Pereira Moises, Andrelson Wellington Rinaldi . Evaluation of the synthetic methods for preparing metal organic frameworks with transition metals. AIMS Materials Science, 2018, 5(3): 467-478. doi: 10.3934/matersci.2018.3.467
    [8] Endi Suhendi, Andini Eka Putri, Muhamad Taufik Ulhakim, Andhy Setiawan, Syarif Dani Gustaman . Investigation of ZnO doping on LaFeO3/Fe2O3 prepared from yarosite mineral extraction for ethanol gas sensor applications. AIMS Materials Science, 2022, 9(1): 105-118. doi: 10.3934/matersci.2022007
    [9] Yana Fajar Prakasa, Sumari Sumari, Aman Santoso, Muhammad Roy Asrori, Ririn Cahyanti . The performance of radar absorption of MnxFe3–xO4/rGO nanocomposites prepared from iron sand beach and coconut shell waste. AIMS Materials Science, 2023, 10(2): 227-248. doi: 10.3934/matersci.2023013
    [10] Tetsuji Saito, Daisuke Nishio-Hamane . Magnetic properties of Sm5Fe17 melt-spun ribbons and their borides. AIMS Materials Science, 2015, 2(4): 392-400. doi: 10.3934/matersci.2015.4.392
  • This paper describes the geotechnical characteristics of index and deformation properties of Stockholm clays. These clays exhibit large variations both geographically and with depth. The clays range from highly organic clays with high plasticity to silty clays with low plasticity. First, the geological conditions of the clays are outlined to qualitatively explain the typical soil stratigraphy encountered in the Stockholm area. Second, two large generic databases are presented, containing 3,500 and 1,600 data points, respectively. The data originates from routine testing and constant rate of strain oedometer tests conducted in commercial projects. The data is analyzed and compared with results from high quality block sampling. It is seen that a common feature of the clays are low undrained shear strengths, and consequently low yield stresses and oedometer moduli. Further, several deformation properties such as preconsolidation pressure and several oedometer moduli is shown to depend on the soils natural water content and its plasticity. Differences in sample quality is shown to highly affect some properties, highlighting the importance of quality sampling and handling of samples. Criteria for sample quality for this type of clay is proposed based on the oedometer moduli before and after the preconsolidation pressure. The paper can hopefully work as a useful reference to engineers working on similar soils worldwide.



    Fe2O3 material is the most common oxide material, very environmentally friendly and used in many fields of science and technology: Welding [1], magnetism [2], supercapacitor electrode [3], biomedicine [4], environment [5], ceramic materials [6,7,8,9,10], environmental contaminant removal [11], photocatalysis [12], chemical sensors [13], electrochemistry [14], magnetic resonance imaging, data storage [15,16,17,18,19], energy [20]. Energy with bandgap (Eg) of Fe2O3 materials with absorption range of 1.9-2.2 eV [21] that corresponds to half of the visible light region of wavelength (λ), λ > 600 µm, absorbing 40% of the sun's spectrum and this is a stable material in solutions with concentrations (pH), pH > 3.0 [22]. Fe2O3 exists in 5 allotropic states: α-Fe2O3 (hematite), β-Fe2O3 [23,24], γ-Fe2O3 (maghemite), ε-Fe2O3, [9,10,25] and ζ-Fe2O3 [25]. In literature, the polymorphic existence of amorphous Fe2O3 was investigated at temperature (T), T = 923 K [26,27,28]. Among them, ε-Fe2O3 was first identified in 1998 [29] and is considered a remarkable phase, with room temperature magnetic resonance T = 300 K [7,10,30] and is applied in the technology of making Jian ceramics of ancient China by heating clay containing Fe, T = ∼1573 K [31]. While, at T = 923 K, Fe2O3 can be converted to Fe3O4 with the energy (E) that decreases very quickly [32,33,34]. On the other hand, the polymorphic phase transition from γ-Fe2O3 → ε-Fe2O3 → α-Fe2O3 has been studied in many previous works. The obtained results show that the material has a thermodynamically unstable structure [10,35]. Wheeler et al. [36] have successfully fabricated single-crystal Fe2O3 nanowires with high density, which have reduced the charge loss by electron-hole recombination. Also, it is known that the magnetic Curie temperature (Tc) of the ferrimagnet material is respectively Tc = 500/585 K [37,38,39]. To study and fabricate Fe2O3 materials based on the experimental methods [40,41], there are the following methods: X-ray diffraction (XRD), Extended X-ray Absorption Fine Structure (EXAFS) method [42], transmission electron microscopy, and High-Resolution Transmission Electron Microscopy (HRTEM) [43]. Schaetz et al. [44] have successfully determined the influence of nano-size, shape and surface elevation effect of nanomaterials on the structure of hematite (α-Fe2O3), and maghemite (γ-Fe2O3). In practical applications, Fe2O3 is also used in the separation of water by electrolysis of α-Fe2O3. Hematite electrodes can generate photons by illumination with a wavelength (λ), λ = 350 μm, without the need for a power source [45,46]. Fe2O3 thin films have been used in photovoltaic cells [47] with hematite nanorods (α-Fe2O3) with thicknesses from 5 nm to 7 nm with lengths from 3 nm to 4 nm [48]. In practical applications [49,50], when it is welded Fe2O3-doped Austenitic steel, better results are obtained than previous conventional methods, due to the self-reversing energy of Marangoni flow at the lower melting temperature of nanowire flux. Similarly, the practical applications of supercapacitor electrodes have also been solved by optimizing the structure of Fe2O3 nanotubes based on polyaniline (PANI) combination [3]. Besides, under the influence of Fe2O3, the positions of defects and oxygen-containing functional groups are improved significantly with band gap (Eg), Eg < 2.8 eV. Because of that, composites Fe2O3 are considered as materials that have visible light with photocatalysts and exhibit photocatalytic activity towards the degradation of cationic dye Rhodamine B (RhB) [12]. O. Iglesias and A. Labarta [51] applied the theoretical method (the Monte Carlo (MC) simulation method combined with periodic boundary conditions) to study the influence of the surface layer on the magnetic properties of nanomaterials at the temperature T = 350 K. VV Hoang and Khanh BTHL used the Molecular Dynamics (MD) method for the determination of the static and thermodynamic properties of liquid and amorphous Fe2O3 nanoparticles and concluded that with amorphous Fe2O3 nano, at T = 3500 K, the Fe-Fe, Fe-O, and O-O corresponding to links lengths are rFe-Fe = 3.1 Å, rFe-O = 1.9 Å, rO-O = 2.7 Å [52]. In other studies, researchers found with the experimental method that Fe2O3 has rFe-Fe = 3.35 Å, rFe-O = 1.91/1.95 Å, rO-O = 2.91 Å [33,53] and rFe-O = 2.08 Å [54]. Among these research methods, the MD simulation method is considered to have many advantages. With a simple method, just by the determination of the interaction between atoms by Newton's second law, together with the averaging statistics method, and applying the computational methods, it serves as a complement to conventional experiments. Based on these obtained results it can be concluded that the nanomaterials have a link length (r) always smaller than that of the bulk material. With the nanomaterials, no pick separation at the first pick of the Radial Distribution Function (RDF) can be observed. With amorphous material, the pick separation occurs at the first vertex of RDF and has 02 Picks, the first Pick has a link length of r11 = 3.11 Å and the second Pick has a link length of r12 = 3.40 Å [55] with link length rFe-Fe = 3.11 Å related to octahedral link pair. Furthermore, the positions of Fe atoms are link to their edges and gradually increase to rFe-Fe = 3.40 Å. Besides, the Coordination Number (CN) of nanomaterials increases when is increased the nano-size. The average coordination number (CN), CNFe-O the range of 4.71-5.28, shows that α-Fe2O3, β-Fe2O3 tend to form an octahedral FeO6 structure in the liquid state with T = 3500 K [54,56,57,58]. Similarly, the bond angle Fe-O-Fe of liquid Fe2O3 nanomaterials at T = 3500 K ranges from 83° to 99°. It is known that for the ideal FeO6 configuration, the O-Fe-O angle is 90°, and for the ideal FeO4 it is 109.5°. However, β-Fe2O3 and ε-Fe2O3 polymorphic structures with only superparamagnetic Fe2O3 nanomaterials [59] have octahedral structural units number of FeO6, FeO4 tetrahedra in ε-Fe2O3 [60,61,62]. The glass temperature (Tg) is an important parameter used for studies about amorphous materials. The glass temperature of Fe2O3 nanomaterials usually is calculated through the intersection of high linear temperature and low temperature extrapolated to the total energy. The values Tg = 2182.18 K, 1897.26 K, 1784.19 K, and 1610.33 K of nanomaterials correspond for the nano-size (D) (2, 3, 4, and 5 nm). Besides, there are also results showing that Tg increases with the increases of D [63,64]. For Fe2O3 materials, Tg = 1400 K has a much smaller value than nanomaterials. Based on the existing data it can be concluded that the results obtained are contradictory, and there are not yet results covering all types of materials such as nanomaterials, thin films, and bulk materials [63,64,65,66]. Besides, in some studies is also determined the phase transition temperature (Tm) [67]. It has been demonstrated that the obtained results depend on the specified research method and applied boundary conditions [66]. Well-documented studies highlighted that there are various factors such as the atomic number (N), temperature (T), pressure (P), surface elevation (h) of the earth's, phase transition process, and heterogeneous kinetic, that can influence the oxide materials CaSiO3 [68], MgSiO3 [69,70]. In these reports, the researchers found a shift in the number of SiOx, MgOy, CaOz structural units with x, y, z = 4/12. However, so far, there have been no specific studies about the structural characteristics of Fe2O3 materials at the high temperature (T) the range is 300-7000 K, and high pressure (P) from P = 0 GPa to P = 360 GPa. To elucidate, the phase transition process can be studied based on the change of the number of structural units FeO4, FeO5, FeO6 in the Fe2O3 material. In our study are given the results of molecular dynamics simulations to determine the structure, phase transition process of Fe2O3 material.

    Initially, a Fe2O3 material with 2000 atoms (800 Fe atoms, 1200 O atoms) was randomly seeded into a model, a cube a nano-size (l), determined by the Eq 1:

    l=(Nρ)13=((2mFenFe+3mOnO)ρ)13 (1)

    where: ρ = 5.19 g/cm3, and ρ is the specific gravity of Fe2O3.

    We use this model by the MD simulation method [51,52,71,72,73,74,75,76,77] with the force field of the Born-Mayer (BM) pair interaction potential calculated with the Eq 2 [52,58]:

    Uij(r)=qiqjr+Bijexp(rRij) (2)

    where: qi, qj are the charges of ions Fe, O; r is the distance from the ith atom to the jth atom, and Bij, Rij are the parameters chosen as: BFe-Fe = 0.0 eV, BFe-O = 2453 eV, BO-O = 1500 eV, Rij = 29 pm, qFe = +3e, qO = −2e [52,58].

    To obtain accurate results, we applied periodic boundary conditions to eliminate surface effects.

    With the content of the MD method, it is only necessary to determine the interaction between atoms according to the content of Newton's II law. With the interaction of the force field, the atoms will move with the displacement (r) and velocity (v). Within the scope of the content of the article, the interactive force field includes the interaction force field between atoms that obeys Newton's II law and the electronic interaction force field between atoms. To determine the quantities, we use Verlet's algorithm to integrate the equations of motion and are combined them with the average statistical method to determine the energy (E) of each atom, according to the following Eq 3:

    E=U+K;U=i>jUij(rij),K=Ni=1miv2i2 (3)

    The rate of increase in temperature (T) and pressure (P) can be calculated with the Eqs 4, 5:

    T(t)=13NkBNi=1mi(vi(t))2 (4)
    P(t)=NVkBT(t)+13Vi<jrij(t)Fij(t) (5)

    where: U, K, N, V are the potential energy, kinetic energy, atomic number, material volume, and kB = 1.381 × 10−23 J/K is the Boltzmann constant.

    The heating process of the model is carried out according to the Eq 6 [52,58]:

    T=T0γt (6)

    where: γ is the cooling coefficient, γ = 4.357 × 1013 Ks−1.

    Initially, it was considered the Fe2O3 material model at a temperature (T), T = 7000 K, and it was run 4 × 104 steps of NVT, with NVP recovery statistics (N is the atomic number, V is volume, T is temperature, and pressure P is constant). The simulation was made for each time step of 0.5 fs so that the model discards the initial configuration and stabilizes at the desired temperature (T), and pressure (P). Then, it was decreased T from T = 7000 K to T = 300 K with a heating rate of dT = 5 K/fs to convert the material model to an amorphous state. This process was followed by an increasing of the temperature T from T = 300 K to T = 900, 1300, 1900, 2300, 2900, 3400, 4500, 5700, 7000 K at P = 0 GPa. Similarly, at T = 300, 2300, 4500, 7000 K, it was increased the pressure P from P = 0 GPa to P = 50, 150, 250, 360 GPa with a turbocharge rate dP = 0.001 GPa/fs with 4 × 104 regression steps that fits all models. After it was obtained the equilibrium in these models, was applied for all models simultaneously running with 107 NVE steps (N is the number of atoms, V is the volume, and E is the energy-being constant). The structural results of the models were determined by an average process after 100 simulation steps. During the computation, for all models, the same conditions (N, V, T, P, and E constant) were applied. To study the structure, phase transition process, and heterogeneous kinetics, data analysis, and visualization methods were used. With the obtained results, the structural features were investigated through the radial distribution function (RDF) with the Eq 7 [68,69,70]:

    g(r)=VN2ini(r)4πr2Δr (7)

    where: r, ni(r), V, N, g(r) are the distance between atoms, the density of atoms, volume, number of atoms, and probability of finding an atom between r to r + Δr.

    In addition, there are other quantities such as nano-size (1), angular distribution, and average coordination number (CN) that are calculated using the following Eq 8 [68,69,70]:

    CN=4πρr10g(r)r2dr (8)

    where: r1 is the first peak position of the RDF, and the phase transition temperature (Tm) is determined through the relationship between E and T.

    For the bond angle, the relationship between the Fe-O-Fe and O-Fe-O bond angle is used for the link lengths (applied for the links O-O, Fe-O, Fe-Fe), which are determined by the Eq 9:

    cosα=2r2FeOr2FeFe2r2FeO (9)

    where: α = Fe-O-Fe or O-Fe-O for the model defined at T and P. During heating and pressure, the model annealing is determined by the Nosé-Hoover formula [78,79]. To confirm the accuracy of the results, it was compared our results with those obtained previously under the same T and P conditions. All these results were written in the C programming language and were computed on the computer central system of the Institute of Physics, the Department of Physics and Astronomy of Zielona Gora University, Poland.

    To study the structure, and phase transition process of Fe2O3 materials with atomic number (N), N = 2000 atoms at temperature (T), T = 300, 900, 1300, 1900, 2300, 3400, 4500, 5700, 7000 K at pressure (P), P = 0 GPa has applied the simulation method of molecular dynamics (MD) with pair interaction potential Born-Mayer (BM) and conditions periodic boundary, according to Eq 2. The results of the shape, nano-size (l), energy (E) of the system, and radial distribution function (RDF) of links Fe-Fe, Fe-O, O-O pairs are different according to Eq 7 (Figure 1, and Table 1).

    Figure 1.  The shape (Fe is red, O is blue) (a), the radial distribution function (b) of materials Fe2O3 material at temperature T = 300 K, pressure P = 0 GPa.
    Table 1.  The structural characteristics of materials Fe2O3 material at different temperatures with the pressure P = 0 GPa.
    T(K) rij(Å) gij(r) 1 (nm) E (eV)
    Fe-Fe Fe-O O-O Fe-Fe Fe-O O-O
    300 3.14 2.00 2.66 3.20 8.75 2.49 2.92 −29.56
    900 3.24 1.98 2.72 2.88 7.23 2.28 2.926 −29.40
    1300 3.32 1.98 2.78 2.83 6.53 2.30 2.932 −29.29
    1900 3.32 1.98 2.82 2.73 5.87 2.23 2.946 −29.13
    2300 3.32 1.96 2.84 2.72 5.44 2.20 2.936 −29.01
    3400 3.36 1.96 2.86 2.57 4.73 2.05 2.98 −28.67
    4500 3.40 1.92 2.94 2.40 4.56 1.88 3.05 −28.24
    5700 3.42 1.94 2.92 2.22 4.26 1.79 3.10 −27.86
    7000 3.46 1.88 2.96 2.09 4.18 1.71 3.19 −27.41
    previous results 3.35 [33,34,53]; 3.11 [33,55,58];
    3.1 [52]
    1.91/1.95 [33,34,53];
    2.08 [54];
    1.925/1.945 [33,53]; 1.9 [33,42,52,55,58]
    2.91 [34,33,53];
    2.7 [33,42,52,55,58]

     | Show Table
    DownLoad: CSV

    The results show that the Fe2O3 material at T = 300 K, P = 0 GPa has a cubic shape, nano-size (1), 1 = 2.92 nm, the energy (E), E = −29.56 eV and it is generated by two atoms Fe and O (Fe in red, O in blue) (Figure 1a). In addition, the Fe2O3 material is made by links Fe-Fe, Fe-O, O-O with the corresponding links lengths rFe-Fe = 3.14 Å, rFe-O = 2.00 Å, rO-O = 2.66 Å (Figure 1b). In this case, the obtained results are completely consistent with experimental results with rFe-Fe = 3.35 Å, rFe-O = 1.93 Å, rO-O = 2.91 Å [34] and rFe-Fe = 3.35 Å, rFe-O = 1.91/1.95 Å, rO-O = 2.91 Å [33,53] and rFe-O = 2.08 Å [54]; whereas the simulation results with rFe-Fe = 3.11 Å, rFe-O = 1.90 Å, rO-O = 2.70 Å [33,55,58] and rFe-Fe = 3.1 Å, rFe-O = 1.9 Å, rO-O = 2.7 Å [52]. It can be noted that the link length of bulk materials is always larger than the link length of nanomaterials with nanomaterials is rFe-Fe = 3.10 Å, rFe-O = 1.9 Å, rO-O = 2.7 Å [33,42,52]. The Fe-Fe link length appears a Pick separation at the first vertex of the RDF into two small peaks, r11 = 3.14 Å, r12 = 3.4 Å, which shows that the Fe2O3 material, at T = 300 K, exists in the amorphous state, and confirms this result as consistent with previous simulation results [33,34,55]. Besides, the Fe-Fe link length of r11 = 3.11 Å is related to the number of FeO4 structural units, whereas r12 = 3.40 Å is related to the number of FeO6 structural units and has a common peak [55]. When it is increased the temperature from T = 300 K to T = 900, 1300, 1900, 2300, 3400, 4500, 5700, 7000 K, l increases from 1 = 2.92 nm to 1 = 2.926, 2.932, 2.946, 2.936, 2.98, 3.05, 3.10, 3.19 nm altitude and the energy (E) increases from E = −29.56 eV to E = −29.40, −29.29, −29.13, −29.01, −28.67, −28.24, −27.86, −27.41 eV. It can be seen that the influence of temperature on the link lengths r of Fe-Fe, Fe-O, O-O of Fe2O3 bulk is very large as rFe-Fe and increases from rFe-Fe = 3.14 Å to rFe-Fe = 3.24, 3.32, 3.32, 3.32, 3.36, 3.40, 3.42, 3.46 Å; rFe-O decreases from rFe-O = 2.00 Å to rFe-O = 1.98, 1.98, 1.98, 1.96, 1.96, 1.92, 1.94, 1.88 Å; rO-O increases from rO-O = 2.66 Å to rO-O = 2.72, 2.78, 2.82, 2.84, 2.86, 2.94, 2.92, 2.96 Å corresponding to the height of the first peak. Also, the radial distribution function g(r) decreases as gFe-Fe increases from gFe-Fe = 3.20 to gFe-Fe = 2.88, 2.83, 2.73, 2.72, 2.57, 2.40, 2.22, 2.09; gFe-O decreases from gFe-O = 8.75 to gFe-O = 7.23, 6.53, 5.87, 5.44, 4.73, 4.56, 4.26, 4.18; gO-O increases from gO-O = 2.49 to gO-O = 2.28, 2.30, 2.23, 2.20, 2.09, 1.88, 1.79, 1.71 (Table 1). On the other hand, during the study, it was found that in the temperature range from T = 300 K to T = 900, 1900, 2300 K and T = 2300 K to T = 3400, 4500, 5700, 7000 K, the g(r) of Fe-O decreases to the smallest value, showing that the Fe2O3 material was a phase transition from the amorphous state to the liquid state. The phase transition of Fe2O3 materials is a type 1 phase transition. During the heating process of Fe2O3 materials, Eq 6 is always applied and the accuracy is evaluated according to the Nosé-Hoover formula [78,79]. To confirm the influence of temperature on structural change, and heterogeneous kinetics, it was determined the average coordination number (CN) according to the Eq 8, and the energy (E) at high temperature (T) (Table 2).

    Table 2.  Coordination number (CN) of link Fe-Fe, Fe-O, O-O and the energy (E) with Fe2O3 material at different temperatures and P = 0 GPa.
    T(K) CN
    Fe-Fe Fe-O O-Fe O-O
    300 11.20 5.44 3.63 15.95
    900 11.42 5.35 3.57 15.90
    1300 11.40 5.33 3.55 15.82
    1900 11.32 5.29 3.53 15.64
    2300 11.67 5.60 3.73 15.93
    3400 11.06 5.11 3.41 15.28
    4500 11.35 5.34 3.56 14.37
    5700 9.63 4.64 3.09 13.78
    7000 8.63 4.12 2.75 3.55

     | Show Table
    DownLoad: CSV

    The results show that, as the temperature increases, the average coordination number (CN) of the link Fe-Fe changes from CNFe-Fe = 11.20 to CNFe-Fe = 11.42, 11.40, 11.32, 11.67, 11.06, 11.35, 9.67, 8.63; Fe-O changes from CNFe-O = 5.44 to CNFe-O = 5.35, 5.33, 5.29, 5.60, 5.11, 5.34, 4.64, 4.12; O-Fe changes from CNO-Fe = 3.63 to CNO-Fe = 3.57, 3.55, 3.53, 3.73, 3.41, 3.56, 3.09, 2.75; O-O changes from CNO-O = 15.95 to CNO-O = 15.90, 15.82, 15.64, 15.93, 15.28, 14.37, 13.78, 3.55 (Table 2), and the energy (E) increases (Table 1).

    In this case, the average coordination number of Fe-O decreases from CN = 5.44 to CN = 5.26; this is the FeO5 structure of Fe2O3 material. When it is increased the temperature from T = 300 K to T = 900, 1900, 2300, 3400, 4500, 5700, 7000 K, the average coordination number (CN), and the bonding density of Fe-O, O-Fe tends to decrease, indicating that the oxide material tends to shift to the liquid state. Based on these results, we can realize that the increase of the temperature leads to an increase in the size (1) and E of the material, which proves that there is a major influence of temperature. With extremes, the distribution number of coordination has the value CN > 5.26 with Fe-O and CN > 3.51 with O-Fe, of which the majority is CN = 5.0 with Fe-O, and at CN = 4.0 with O-Fe. On the other hand, at T = 300 K, the amorphous Fe2O3 material has a distorted FeO6 structure with the average coordination number of the Fe-O pair with CNFe-O value ≈ 6 and it was considered that these FeO6 structural units are as structural defects [58]. In addition, to analyze the structure of the Fe2O3 material oxide material, it was applied the visualization method to study the number of structural units of the material (Figure 2).

    Figure 2.  (a) The shape of the number of the structural units Fe2O3 material as FeO4, (b) FeO5, (c) FeO6, (d) structure shape at T = 300 K, P = 0 GPa.

    The Fe2O3 material at T = 300 K, P = 0 GPa exists with 03 structural units: FeO4 corresponds to one Fe atom link to four O atoms (Figure 2a), similar with FeO5 (Figure 2b), FeO6 (Figure 2c), and the structure shape of Fe2O3 material with the number of the structural unit of FeO4 is black, FeO5 is red, FeO6 is blue at T = 300 K, P = 0 GPa (Figure 2d). The number of structural units in the Fe2O3 material at different temperatures (T) is shown in Table 3.

    Table 3.  Number of atoms with structural units FeO4, FeO5, FeO6 of Fe2O3 material at different temperatures, P = 0 GPa.
    Temperature (K) FeO4 FeO5 FeO6
    300 419 1627 1080
    900 453 1647 1018
    1300 540 1621 1018
    1900 676 1640 850
    2300 789 1552 833
    3400 1175 1501 462
    4500 1551 1225 209
    5700 1634 1949 201
    7000 1634 1049 201

     | Show Table
    DownLoad: CSV

    The obtained results show that, when the temperature is increased, the number of structural units FeO4 increases from 419 to 453, 540, 676, 789, 1175, 1551, 1634, 1634; FeO5 changes about from 1627 to 1647, 1621, 1640, 1552, 1501, 1225, 1949, 1049; FeO6 decrease from 1080 to 1018, 1018, 850, 833, 462, 209, 201, 201 (Table 3), showing that the Fe2O3 material tends to shift to the liquid state, which demonstrates the influence of temperature on link lengths and the number of structural units that are very big. In addition, the angular distributions of Fe-O-Fe and O-Fe-O are shown in Table 4.

    Table 4.  The angular distribution of material Fe2O3 material at different temperatures.
    Temperature (K) Angular distribution
    Fe-O-Fe O-Fe-O
    300 105 95
    900 105 95
    1300 100 95
    1900 100 95
    2300 100 95
    3400 100 95
    4500 95 95
    5700 95 95
    7000 95 95

     | Show Table
    DownLoad: CSV

    The obtained results show that, when increasing the temperature, the angle distribution Fe-O-Fe, O-Fe-O between Fe and O atoms has a negligible change around the value of 105° with Fe-O-Fe and with O-Fe-O is 95°. When it is increased the temperature from T = 300 K to T = 900, 1900, 2300, 3400, 4500, 5700, 7000 K, Fe-O-Fe decreases from 105° to 95° and O-Fe-O has a constant value is 95° (Table 4). The obtained results are completely consistent with the simulation results [55,58].

    For O-Fe-O bond angle nanomaterials, O-Fe-O has major peaks between 99° and 83°. While the ideal uniform octahedral model has a value of 90°, this is considered the ideal tetrahedral value. The value of 109.5°, confirms the reason that there is a moving to the liquid state, and there are structural units of FeO5.

    The results show that the bond angle O-Fe-O reduced from 105° to 90° when the material changes to the liquid state. These results raise the question of whether the bond angle depends on the length of the links. Through the computational process, was applied the Eq 9, and it was highlighted when rFe-Fe, rO-O increases, rFe-O decreases. Besides, there is a relationship between the temperature (T) and the energy (E), (Figure 3).

    Figure 3.  Phase transition temperatures of material Fe2O3 material at P = 0 GPa, with different temperatures.

    The results show that, when it is increased the temperature (T), form 300-2300 K to 2300-7000 K, the energy (E) of Fe2O3 material increases from E = −29.56 eV to E = −29.00 eV and from E = −29.00 eV to E = −27.41 eV. When T increases, T = 300-2300 K then E increases slowly, and between T = 2300 K to T = 7000 K this leads to an E rapid increase. These results show that from T = 300 K to T = 2300 K the Fe2O3 material exists in an amorphous state, and E increases slowly. When there is the Fe2O3 material at state liquid, E increases rapidly, due to the intersection between the fast-rising energy region and the slowly increasing energy region of 2300 K and calling this is as a point of the phase transition or phase transition temperature (Tm), Tm = 2300 K, corresponding to E = −29.00 eV (Figure 3). This leads to an abrupt change in the density of atoms. When it is increased T, this leads to an increase in nano-size (1), E increases, and the results are like those in simulations of CaSiO3 [68], MgSiO3 [69,70] materials.

    The phase transition temperature is one of the most important parameters of amorphous materials, and only the glass temperature (Tg) of Fe2O3 nanomaterials was found in the scientific literature [52]. The obtained results show that Tg depends greatly on the nanoparticle nano-size. When it is increased the nano-size (D) of nanomaterials is D = 2-5 nm, the Tg decreases Tg = 2182-1610 K. These results appear to be incomplete in contradiction to previously published bulk materials. In addition, the characteristics of the phase transition from the liquid state to the crystalline state are still unclear with the statements for previously published papers [65,66]. With the phase transition temperature Tm, there is a lack of consensus on studying the nature of the phase transition [63,64,67]. The obtained results show that there is a great influence of temperature on the structure, phase transition process of Fe2O3 material.

    Similarly, it was studied the effect of pressure (P), P = 0, 50, 150, 250, 360 GPa at temperature (T), T = 300, 2300, 4500, 7000 K, to elucidate the structure, phase transition process of Fe2O3 material.

    The Fe2O3 material at T = 300 K, leads to an P increasing from P = 0 GPa to P = 50, 150, 250, 360 GPa with structural characteristic results (Figure 4).

    Figure 4.  (a) Structural features of material Fe2O3 material as structure shape; the links Fe-Fe, Fe-O, O-O of (b) r, (c) g(r); Fe-Fe, Fe-O, O-Fe, O-O at T = 300 K, different P, (d) 1, (e) E, (g) different P.

    The results show for the material Fe2O3 material at T = 300 K, P = 0 GPa, that it has a structure shape (Figure 4a) and the Fe-Fe, Fe-O, O-O bonding length of Fe-Fe is rFe-Fe = 3.14 Å, rFe-O = 2.0 Å, rO-O = 2.66 Å.

    Also, the first peak height of the radial distribution function (RDF) of the links is gFe-Fe = 3.20, gFe-O = 8.75, gO-O = 2.49. When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, then rFe-Fe decreases from rFe-Fe = 3.14 Å to rFe-Fe = 3.10, 2.86, 2.74, 2.66 Å, and gFe-Fe increases, and decreases from gFe-Fe = 3.20 to gFe-Fe = 3.30, 3.20, 2.83, 2.76; rFe-O decreases from rFe-O = 2.0 Å to rFe-O = 1.92, 1.86, 1.84, 1.84 Å, and gFe-O decreases from gFe-O = 8.75 to gFe-O = 7.28, 6.86, 6.49, 5.99; rO-O decreases from rO-O = 2.66 Å to rO-O = 2.66, 2.50, 2.38, 2.24 Å, and gO-O changes from gO-O = 2.49 to gO-O = 2.72, 2.66, 2, 76, 2.72 (Figure 4b, c); size (1), decreases from 1 = 2.92 nm to 1 = 2.72, 2.55, 2.44, 2.36 nm (Figure 4d), the energy (E) increases from E = -29.56 eV to E = -29.22, -28.34, -27.33, -26.27 eV (Figure 4e) corresponding to the mean coordinate number (CN), CNFe-Fe increases from CNFe-Fe = 11.20 to CNFe-Fe = 13.91, 16.42, 18.69, 20.70; CNFe-O increases from CNFe-Fe = 5.44 to CNFe-Fe = 6.28, 7.25, 8.12, 8.73; CNO-Fe increases from CNFe-Fe = 3.63 to CNFe-Fe = 4.19, 4.83, 5.41, 5.82; CNO-O increases from CNFe-Fe = 15.95 to CNFe-Fe = 18.64, 23.02, 26.80, 30.07 (Figure 4g). In addition, the bonding angle of Fe-O-Fe, O-Fe-O and the number of structural units FeO4, FeO5, FeO6 are determined by Figure 5.

    Figure 5.  (a) The angular distribution, the number of atoms with structural units FeO4, FeO5, (b) FeO6, of Fe2O3 material at T = 300 K with different pressure.

    The results show at P = 0 GPa that the angle of Fe-O-Fe is 105°, and of O-Fe-O is 95° corresponding to the number of structural units FeO4, FeO5, FeO6 is 419, 1627, 1080. When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, then Fe-O-Fe decreases from 105° to 0°, O-Fe-O is 95° (Figure 5a) and FeO4, FeO5 decreases from 419 to 48, 15, 0, 0; FeO5 decreases from 1627 to 1217, 336, 63, 12; FeO6 increases from 1080 to 1638, 1936, 1989, 1998 (Figure 5b). The obtained results are completely consistent with the results from the literature [80]. Besides, it also shows that r, g(r) decreases when Fe2O3 material is at high pressure and low temperature.

    Similarly, with Fe2O3 material at T = 2300 K, when P increases from P = 0 GPa to P = 50, 150, 250, 360 GPa, there are structural characteristic results (Figure 6).

    Figure 6.  (a) Structural features of material Fe2O3 material as structure shape; (b) r, (c) g(r) of the links Fe-Fe, Fe-O, O-O; (d) 1, (e) E, (g) CN of Fe-Fe, Fe-O, O-Fe, O-O at T = 2300 K, different P.

    The results show, when Fe2O3 material at T = 2300 K, P = 0 GPa, there is a structure shape (Figure 6a), with the length of Fe-Fe, Fe-O, O-O links is 3.32 Å, rFe-O = 1.96 Å, respectively rO-O = 2.84 Å, the first peak radial distribution function (RDF) of the links is gFe-Fe = 2.72, gFe-O = 5.44, gO-O = 2.22 (Figure 6a). When it is increased P from P = 0.0 GPa to P = 50, 150, 250, 360 GPa, rFe-Fe decreases from rFe-Fe = 3.32 Å to rFe-Fe = 3.14, 2.96, 2.80, 2.77 Å, and gFe-Fe increased, decreases from gFe-Fe = 2.72 to gFe-Fe = 2.94, 2.72, 2.70, 2.51; rFe-O decreases from rFe-O = 1.96 Å to rFe-O = 1.92, 1.86, 1.84, 1.82, and gFe-O increased, decreases from gFe-O = 5.44 to gFe-O = 6.29, 5.52, 5.42, 5.18; rO-O decreases from rO-O = 2.84 Å to rO-O = 2.64, 2.48, 2.32, 2.22 Å, and gO-O increased, decreases from gO-O = 2.20 to gO-O = 2.44, 2.53, 2.53, 2.48 (Figure 6b, c); nano-size (1), decreases from 1 = 2.94 nm to 1 = 2.72, 2.55, 2.44, 2.36 nm (Figure 6d) and energy (E) increases from E = −29.01 eV to E = −28.78, −27.92, −26.97, −26.00 eV (Figure 6e) corresponding to the mean coordinate number (CN), CNFe-Fe increases from CNFe-Fe = 11.67 to CNFe-Fe = 14.07, 16.61, 18.87, 20.60; CNFe-O increases from CNFe-O = 5.60 to CNFe-O = 6.42, 7.44, 8.32, 8.88; CNO-Fe increases from CNO-Fe = 3.73 to CNO-Fe = 4.28, 4.96, 5.55, 5.92; CNO-O increases from CNO-O = 15.93 to CNO-O = 18.99, 23.32, 26.99, 30.42; (Figure 6g). The same results are also obtained, with T = 2300 K (Figure 7).

    Figure 7.  (a) The angular distribution, (b) the number of atoms with structural units FeO4, FeO5, FeO6 of Fe2O3 material at T = 2300 K with different pressure.

    The obtained results show at P = 0 GPa the bonding angle of Fe-O-Fe is 100°, O-Fe-O is 95° corresponding to the number of structural units (FeO4, FeO5, FeO6 is 789, 1552, 833). When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, then Fe-O-Fe decreases from 100° to 0°, O-Fe-O increases from 95° to 110° (Figure 7a), and FeO4 decreases from 789 to 170, 15, 0, 0; FeO5 decreases 1552 to 1153, 353, 74, 0, and FeO6 increases from 833 to 1661, 1933, 1987, 2000 (Figure 7b). In it, all r, g(r), 1 decrease and E, CN increase but a smaller decrease at T = 300 K; it shows that in the phase transition temperature Tm = 2300 K, Fe2O3 material structure changes to the liquid state, so Fe2O3 material structure changes is more stable, when increasing P, it leads to changes in smaller structural features.

    Similarly, when Fe2O3 material at T = 4500 K, P is increased from P = 0 GPa to P = 50, 150, 250, 360 GPa there are the structural characteristic results (Figure 8).

    Figure 8.  (a) Structural features of material Fe2O3 material as structure shape; (b) r, (c) g(r) of the links Fe-Fe, Fe-O, O-O; (d) 1, (e) E, (g) CN of Fe-Fe, Fe-O, O-Fe, O-O at T = 4500 K, different P.

    The results show when Fe2O3 material at T = 4500 K, P = 0 GPa has a structure shape (Figure 8a), and the length of Fe-Fe, Fe-O, O-O links respectively, rFe-Fe is 3.40 Å, rFe-O is 1.92 Å, rO-O is 2.94 Å, the first peak radial distribution function (RDF) of the links is gFe-Fe = 2.40, gFe-O = 4.56, gO-O = 1.88. When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, then rFe-Fe decreases from rFe-Fe = 3.40 Å to rFe-Fe = 3.10, 2.92, 2.78, 2.74 Å, and gFe-Fe increased, decreases from gFe-Fe = 2.40 down to gFe-Fe = 2.57, 2.52, 2.51, 2.41; rFe-O decreases from rFe-O = 1.92 Å to rFe-O = 1.90, 1.88, 1.84, 1.84, and gFe-O decreases, increased from gFe-O = 4.56 to gFe-O = 4.52, 4.41, 4.55, 4.57; rO-O decreases from rO-O = 2.94 Å to rO-O = 2.80, 2.46, 2.34, 2.22 Å, and gO-O increases from gO-O = 1.88 to gO-O = 2.12, 2.18, 2.26, 2.39 (Figure 8b, c); that corresponds to nano-size (1), that decreases from 1 = 3.05 nm to 1 = 2.72, 2.55, 2.44, 2.36 nm (Figure 8d), and energy (E) increases from E = −28.24 eV to E = −28.17, −27.47, −26.49, −25.45 eV (Figure 8e) corresponding to the mean coordinate number (CN); CNFe-Fe increases from CNFe-Fe = 11.35 to CNFe-Fe = 14.40, 16.87, 18.99, 20.70; CNFe-O increases from CNFe-O = 5.34 to CNFe-O = 6.56, 7.69, 8.40, 8.91; CNO-Fe increases from CNO-Fe = 3.56 to CNO-Fe = 4.37, 5.12, 5.60, 5.94; CNO-O increases from CNO-O = 14.37 to CNO-O = 19.21, 23.56, 27.35, 30.51 (Figure 8g). The same results are obtained with T = 4500 K (Figure 9).

    Figure 9.  (a) The angular distribution, (b) the number of atoms with structural units FeO4, FeO5, FeO6, of Fe2O3 material at T = 4500 K with different pressure.

    The obtained results show when P = 0 GPa the bonding angle of Fe-O-Fe is 100°, O-Fe-O is 95° corresponding to the number of structural units FeO4, FeO5, FeO6 is 1551, 1225, 209. When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, then Fe-O-Fe decreases from 100° to 0°, O-Fe-O increases from 95° to 100° (Figure 9a) and FeO4 decreases from 1551 to 222, 5, 0, 0; FeO5 decreases from 1225 to 1271, 285, 46, 6, and FeO6 increases from 209 to 1583, 1948, 1992, 1999 (Figure 9b).

    It can be observed that r, g(r), 1 decrease and E, CN increases, but a smaller decrease at T = 300, 2300 K show that in the high-temperature region, the Fe2O3 structure in the liquid state has a smaller structural change.

    Similarly, when Fe2O3 material at T = 7000 K, P increased from P = 0 GPa to P = 50, 150, 250, 360 GPa, there are structural characteristic results (Figure 10).

    Figure 10.  (a) Structural features of material Fe2O3 material as structure shape; (b) r, (c) g(r) of the links Fe-Fe, Fe-O, O-O; (d) I, (e) E, (g) CN of Fe-Fe, Fe-O, O-Fe, O-O at T = 7000 K, different P.

    The results show when Fe2O3 material at T = 7000 K, P = 0 GPa there is a structure shape (Figure 10a) and the length of Fe-Fe, Fe-O, O-O links respectively, rFe-Fe is 3.46 Å, rFe-O is 1.88 Å, rO-O is 2.96 Å, the first peak radial distribution function (RDF) of the links is gFe-Fe = 2.09, gFe-O = 4.18, gO-O = 1.71. When it is increased P from P = 0 GPa to P = 50, 150, 250, 360 GPa, rFe-Fe decreases from rFe-Fe = 3.46 Å to rFe-Fe = 3.10, 2.94, 2.78, 2.68 Å, and gFe-Fe increases, decreases from gFe-Fe = 2.09 to gFe-Fe = 2.37, 2.37, 2.36, 2.29; rFe-O decreases from rFe-O = 1.88 Å to rFe-O = 1.90, 1.88, 1.84, 1.82 Å, and gFe-O decreases from gFe-O = 4.18 to gFe-O = 3.62, 3.73, 3.93, 3.96; rO-O decreases from rO-O = 2.96 Å to rO-O = 2.62, 2.42, 2.28, 2.24 Å, and gO-O increases from gO-O = 1.71 to gO-O = 1.97, 2.05, 2.04, 2.15 (Figure 10b, c). Also, nano-size (1), decreases from 1 = 3.19 nm to 1 = 2.72, 2.55, 2.44, 2.36 nm (Figure 10d), energy (E) increases from E = −27.41 eV to E = −27.51, −26.81, −25.87, −24.83 eV (Figure 10e), corresponding to the average number of coordinates (CN); CNFe-Fe increases from CNFe-Fe = 8.63 to CNFe-Fe = 14.58, 17.01, 19.11, 20.83; CNFe-O increases from CNFe-O = 3.55 to CNFe-O = 6.69, 7.79, 8.49, 8.96; CNO-Fe increases from CNO-Fe = 2.75 to CNO-Fe = 4.46, 5.19, 5.66, 5.98; CNO-O increases from CNO-O = 4.12 to CNO-O = 19.36, 23.79, 27.79, 30.52 (Figure 10g). The results are similar with T = 7000 K (Figure 11).

    Figure 11.  (a) The angular distribution, the number of atoms with structural units FeO4, (b) FeO5, FeO6 of Fe2O3 material at T = 7000 K with different pressure.

    The obtained results show, at P = 0 GPa the bonding angle of Fe-O-Fe is 95°, O-Fe-O is 95° corresponding to the number of structural units FeO4, FeO5, FeO6 is 1634, 1049, 201. When increasing P from P = 0 GPa to P = 50, 150, 250, 360 GPa then Fe-O-Fe decreases from 95° to 0°, O-Fe-O increases from 95° to 100° (Figure 11a) and FeO4 decreases from 1634 to 368, 20, 0, 0; FeO5 changes from 1049 to 1259, 397, 47, 0 and FeO6 increases from 201 to 1520, 1925, 1992, 2000 (Figure 11b).

    In it, all r, g(r), 1 decreases and E, CN increase but with a smaller decrease at T = 2300, 4500 K. All obtained results, when an increase of T = 300-7000 K at P = 0.0 GPa and an increasing of P from P = 0 GPa to P = 50, 150, 250, 360 GPa at T = (300, 2300, 4500, 7000 K) then r, g(r) and 1, E all decrease with T > 4500 K, the decrease led to less, the smallest at T = 7000 K, this result corresponding to CN also increases fewer. In addition, the disappearance in the number of structural units appears with T = 2300 K; with T = 4500 K, FeO4 disappears at all material with P ≥ 250 GPa and Fe-O-Fe disappears with P ≥ 150 GPa. The obtained results are used as the basis for future experimental results such as determining the link length, and the bonding angle between atoms.

    The molecular dynamics simulation method has successfully studied the influence of high temperature and high pressure on the structure and phase transition of Fe2O3 materials. The results show that with increasing temperature, the nano-size (1), energy increase, length link decreases, the number of structural units FeO4, FeO5 increases, and FeO6 decreases. In contrast, as pressure increases, then l decreases, E increases, FeO4, FeO5 decreases, and FeO6 increases. In addition, when P ≥ 150 GPa, FeO4 disappears, FeO5 disappears at P ≥ 250 GPa), then only FeO6 exists at T = 2300, 7000 K. It shows that when T increases, FeO4 increases, FeO5 increases, FeO6 decrease when increasing P gives the opposite result. The cause of this phase transition (the phase transition between the amorphous state to the liquid state and from the liquid to the amorphous state) is caused by the nano-size effect. The phase transition follows a type 1 phase transition. With the phase transition, the phase transition temperature (Tm) has been successfully determined, Tm = 2300 K, which is the phase transition temperature from the amorphous state to the state. liquid. The obtained results form the basis for future experimental studies.

    The authors declare no conflict of interest.



    [1] Atterberg A (1911) Clays and its water contents and plasticity limits[In Swedish]. Kungliga Lantbruksakademiens Handlingar och Tidskrift 50: 132–158.
    [2] Atterberg A (1912) Consistency of soils[In Swedish]. Kungliga Lantbruksakademiens Handlingar och Tidskrift 51: 93–123.
    [3] Casagrande A (1932) Research on the Atterberg limits of soils. Public Roads 13: 121–136.
    [4] Statens Järnvägar (1922) Swedish State Railways Geotechnical Commission 1914–1922.[In Swedish].
    [5] Hansbo S (1957) A new approach to the determination of the shear strength of clay by the fall-cone test, Swedish Geotechnical Institute, Stockholm.
    [6] Jakobson B (1954) Influence of sampler type and testing method on shear strength of clay samples, Statens Geotekniska Institut, Stockholm.
    [7] Kallstenius T (1958) Mechanical disturbances in clay samples taken with standard piston sampler, Swedish Geotechnical Institute, Stockholm.
    [8] Kallstenius T (1963) Studies on clay samples taken with standard piston sampler, Swedish Geotechnical Institute, Stockholm.
    [9] Cadling L, Odenstad D (1950) The vane borer. An apparatus for determining the shear strength of clay soils directly in the ground, Swedish Geotechnical Institute, Stockholm.
    [10] Bjerrum L, Flodin N (1960) The development of soil mechanics in Sweden, 1900–1925. Géotechnique 10: 1–18. https://doi.org/10.1680/geot.1960.10.1.1 doi: 10.1680/geot.1960.10.1.1
    [11] Lundin S (2000) Geotechnics in Sweden 1920–1945, Swedish Geotechnical Society Report 1.
    [12] Larsson R (1986) Consolidation of soft soils, Swedish Geotechnical Institute.
    [13] Massarsch R, Aronsson S (2014) The Heritage of Swedish Foundation Engineering. Proceedings, DFI/EFFC International Conference on Piling and Deep Foundations, Stockholm, 1031–1054.
    [14] Phoon KK, Ching J, Cao Z (2022) Unpacking data-centric geotechnics. Underground Space 7: 967–989. https://doi.org/10.1016/j.undsp.2022.04.001 doi: 10.1016/j.undsp.2022.04.001
    [15] Phoon KK, Kulhawy FH (1999) Evaluation of geotechnical property variability. Can Geotech J 36: 625–639. https://doi.org/10.1139/t99-039 doi: 10.1139/t99-039
    [16] Phoon KK, Ching J, Shuku T (2022) Challenges in data-driven site characterization. Georisk Assess Manage Risk Eng Syst Geohazards 16: 114–126. https://doi.org/10.1080/17499518.2021.1896005 doi: 10.1080/17499518.2021.1896005
    [17] TC304, TC304 databases. Available from: http://140.112.12.21/issmge/tc304.htm? = 6.
    [18] Swedish Geotechnical Society, Classification and description of soils. Report 1: 2016.
    [19] Swedish Geological Survey, Interactive online maps. 2022. Available from: https://www.sgu.se/om-sgu/nyheter/2022/juni/sgus-strsndforskjutningsmodell-ger-en-bild-av-den-forntida-fordelningen-mellan-hav-och-land/.
    [20] Bjerrum L (1967) Engineering geology of Norwegian normally-consolidated marine clays as related to settlements of buildings. Géotechnique 17: 83–118. https://doi.org/10.1680/geot.1967.17.2.83 doi: 10.1680/geot.1967.17.2.83
    [21] Hov S, Borgström K, Paniagua P (2022) Full-flow CPT tests in a nearshore organic clay. Cone Penetration Test, 452–458.
    [22] Swedish Geological Survey, Description of the geological map Stockholm. 1964.
    [23] Swedish Geotechnical Society, Geotechnical Field Hand Book[In Swedish]. Report 1: 2013.
    [24] Swedish Geotechnical Institute, Swedish Committe on piston samling, Standard piston sampler, Stockholm. 1961.
    [25] Larsson R (1981) Do we obtain undisturbed samples from the standard piston samplers? A comparison with the Laval sampler[In Swedish]. Swedish Geotechnical Insititute.
    [26] Emdal A, Gylland A, Amundsen H, et al. (2016) Mini-block sampler. Can Geotech J 53: 1235–1245. https://doi.org/10.1139/cgj-2015-0628 doi: 10.1139/cgj-2015-0628
    [27] Hov S, Garcia de Herreros C (2020) Comparison of the mini block sampler and the standard piston sampler in a varved East Swedish clay. Report LabMind.
    [28] Karlsson M, Bergström A, Djikstra J (2015) Comparison of the performance of mini-block and piston sampling in high plasticity clays, Technical Report. Chalmers University of Technology, Gothenburg, Sweden.
    [29] Karlsson M, Emdal A, Dijkstra J (2016) Consequences of sample disturbance for predicting long-term settlements in soft clay. Can Geotech J 53. https://doi.org/10.1139/cgj-2016-0129 doi: 10.1139/cgj-2016-0129
    [30] Hov S, Holmén M (2018) The liquid limit[In Swedish]. Swedish Geotechnical Society.
    [31] Hov S, Holmén M (2018) The fall cone test[In Swedish]. Swedish Geotechnical Society.
    [32] Hov S, Prästings A, Persson E, et al. (2021) On Empirical Correlations for Normalised Shear Strengths from Fall Cone and Direct Simple Shear Tests in Soft Swedish Clays. Geotech Geol Eng 39: 4843–4854. https://doi.org/10.1007/s10706-021-01797-w doi: 10.1007/s10706-021-01797-w
    [33] Karlsson R (1981) Consistency Limits in cooperation with the Laboratory Committee of the Swedish Geotechnical Society. Swed Council Build Res, Stockholm, Sweden, Document D6.
    [34] Karlsson R (1961) Suggested improvements in the liquid limit test, with reference to flow properties of remoulded clays. Proc 5th Int Conf Soil Mech Found Eng 1: 171–184.
    [35] Christensen S (2014) Interpretation of field tests, choice of design cuA based on field and laboratory tests[In Norwegian]. Norway's water resources and energy directorate in cooperation with the Norwegian road and railway administrations report no.
    [36] Sällfors G (1975) Preconsolidation pressure of soft, high-plastic clays. Chalmers University of Technology.
    [37] Larsson R, Sällfors G (1986) Automatic Continuous Consolidation Testing in Sweden, Consolidation of Soils: Testing and Evaluation, ASTM STP, 299–328.
    [38] Tidfors M, Sällfors G (1989) Temperature effects on preconsolidation pressure. Geotech Test J 12: 93–97. doi: 10.1520/GTJ10679J
    [39] Boudali M, Leroueil S, Murthy S (1994) Viscous behaviour of natural clays. XⅢ ICSMFE, 411–416.
    [40] Larsson R, Sällfors G (1981) Calculation of settlements in clay[In Swedish]. Väg-och Vattenbyggaren 3: 39–42.
    [41] Paniagua P, L'Heureux JS, Yang SY, et al. (2016) Study on the practices for preconsolidation stress evaluation from oedometer tests. Proc 17th Nord Geotech Meet, 547–555.
    [42] Karlsrud K, Hernandez-Martinez F (2013) Strength and deformation properties of Norwegian clays from laboratory tests on high-quality block samples. Can Geotech J 50: 1273–1293. https://doi.org/10.1139/cgj-2013-0298 doi: 10.1139/cgj-2013-0298
    [43] Janbu N (1963) Soil compressibility as determined by oedometer and triaxial tests. Proc Eur Conf SMFE 1: 19–25.
    [44] Larsson R, Sällfors G, Bengtsson PE, et al. (2007) Shear strength in cohesive soils[In Swedish]. Swedish Geotechnical Institute.
    [45] Berre T, Lunne T, L'Heureux JS (2022) Quantification of sample disturbance for soft, lightly overconsolidated, sensitive clay samples. Can Geotech J 59: 300–303. https://doi.org/10.1139/cgj-2020-0551 doi: 10.1139/cgj-2020-0551
    [46] Ching J, Li DQ, Phoon KK (2016) Statistical characterization of multivariate geotechnical data. Reliability of geotechnical structures in ISO2394, Balkema, CRC Press.
    [47] Bjerrum L (1972) Embankments on soft ground. State of the art report, ASCE Specialty conference on performance of earth and earth-supported structure.
    [48] D'Ignazio M, Phoon KK, Tan SA, et al. (2016) Correlations for undrained shear strength of Finnish soft clays. Can Geotech J 53: 1628–1645. https://doi.org/10.1139/cgj-2016-003 doi: 10.1139/cgj-2016-003
    [49] Ching J, Phoon KK (2012) Modeling parameters of structured clays as a multivariate normal distribution. Can Geotech J 49: 522–545. https://doi.org/10.1139/t2012-015 doi: 10.1139/t2012-015
    [50] Tavenas F, Jean P, Leblond P, et al. (1983) The permeability of natural soils. Part Ⅱ: permeability characteristics. Can Geot J 20: 645–660. https://doi.org/10.1139/t83-073 doi: 10.1139/t83-073
    [51] Larsson R, Bengtsson PE, Eriksson L (1997) Prediction of settlements of embankments on soft, fine-grained soils. Calculation of settlements and their course with time, Swedish Geotechnical Institute.
    [52] Leroueil S, Bouclin G, Tavenas F, et al. (1990) Permeability anisotropy of natural clays as a function of strain. Can Geotech J 27: 568–579. https://doi.org/10.1139/t90-072 doi: 10.1139/t90-072
    [53] Mesri G, Feng T, Ali S, et al. (1994) Permeability characteristics of soft clays. XⅡ ICSMFE, 187–192.
    [54] Larsson R (2008) Properties of soils[In Swedish]. Swedish Geotechnical Institute.
    [55] Larsson R (1977) Basic behaviour of Scandinavian soft clays, Swedish Geotechnical Institute.
    [56] Massarsch R (1979) Lateral earth pressure in normally consolidated clay. Design Parameters in Geotechnical Engineering. Proc 7th Eur Conf Soil Mech Found Eng 2: 245–249.
    [57] L'Heureux JS, Ozkul Z, Lacasse S, et al. (2017) A revised look at the coefficient of earth pressure at rest for Norwegian Clays. NGF Geoteknikkdagen 35.
    [58] DeGroot DJ, Poirier SE, Landon MM (2005) Sample disturbance—soft clays. Studia Geotechnica et Mechanica, Vol. XXVⅡ, 107–120.
    [59] Hight DW, Boese R, Butcher AP, et al. (1992) Disturbance of the Bothkennar clay prior to laboratory testing. Géotechnique 42: 199–217. https://doi.org/10.1680/geot.1992.42.2.199 doi: 10.1680/geot.1992.42.2.199
    [60] Lacasse S, Berre T, Lefevbre G (1985) Block sampling of sensitive clays. Proceedings of the 11th conference on soil mechanics and foundation engineering, San Francisco, CA, 887–892.
    [61] Lunne T, Berre T, Strandvik S (1997) Sample disturbance effects in soft low plastic Norwegian clay. Proceedings of the conference on recent developments in soil and pavement mechanics, Rio de Janeiro, 81–102.
    [62] Clayton CRI, Siddique A (1999) Tube sampling disturbances—fotgotten truths and new perspectives. ICE Geotech Eng 137: 127–135. https://doi.org/10.1680/gt.1999.370302 doi: 10.1680/gt.1999.370302
    [63] Larsson R (2011) SGIs 200 mm diameter block sampler—Undisturbed sampling in fine-grained soils[In Swedish]. Swedish Geotechnical Institute Report 33.
    [64] Löfroth H (2012) Sampling in normal and high sensitive clay—a comparison of results from specimens taken with the SGI large-diameter sampler and the standard piston sampler St Ⅱ. Swedish Geotechnical Institute.
    [65] Andersson M (2022) Personal communication. Swedish Geotechnical Institute.
    [66] Terzaghi K, Peck R, Mesri G (1997) Soil mechanics in engineering practice, Jon Wiley & Sons.
    [67] Donohue S, Long M (2010) Assessment of sample quality in soft clay using shear wave velocity and suction measurements. Géotechnique 60: 883–889. https://doi.org/10.1680/geot.8.T.007.3741 doi: 10.1680/geot.8.T.007.3741
  • This article has been cited by:

    1. Dung Nguyen Trong, Van Cao Long, Ştefan Ţălu, Umut Saraç, Phu Nguyen Dang, Kien Pham Huu, A New Study on the Structure, and Phase Transition Temperature of Bulk Silicate Materials by Simulation Method of Molecular Dynamics, 2022, 6, 2504-477X, 234, 10.3390/jcs6080234
    2. Dung Nguyen Trong, Van Cao Long, Ştefan Ţălu, Effects of Number of Atoms and Doping Concentration on the Structure, Phase Transition, and Crystallization Process of Fe1-x-yNixCoy Alloy: A Molecular Dynamic Study, 2022, 12, 2076-3417, 8473, 10.3390/app12178473
    3. Nguyen Quang Hoc, Dung Nguyen Trong, Nguyen Chinh Cuong, Bui Duc Tinh, Nguyen Duc Hien, Van Cao Long, Umut Saraç, Ştefan Ţălu, Determination of Young Modulus and Stress-Strain Curve for Metal Fe and Interstitial Alloy FeC, 2022, 6, 2504-477X, 250, 10.3390/jcs6090250
    4. Dung Nguyen Trong, Van Cao Long, Ştefan Ţălu, New insights on the factors affecting the heterogeneous kinetics of bulk Fe2O3 on the Earth: A molecular dynamic simulation, 2022, 12, 2158-3226, 065016, 10.1063/5.0094033
  • Reader Comments
  • © 2023 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(2965) PDF downloads(193) Cited by(0)

Figures and Tables

Figures(20)  /  Tables(3)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog