<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">JMMCE</journal-id><journal-title-group><journal-title>Journal of Minerals and Materials Characterization and Engineering</journal-title></journal-title-group><issn pub-type="epub">2327-4077</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/jmmce.2016.42015</article-id><article-id pub-id-type="publisher-id">JMMCE-65162</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Chemistry&amp;Materials Science</subject><subject> Engineering</subject></subj-group></article-categories><title-group><article-title>
 
 
  Numerical Simulation and Experimental Characterization of Clay Paste under Loads for Energy Saving in Clay Materials Processing
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>rahiman</surname><given-names>Traore</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Kokou</surname><given-names>Esso Atcholi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Ouro-Djobo</surname><given-names>Samah</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Samuel</surname><given-names>Gomes</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tibi</surname><given-names>Beda</given-names></name><xref ref-type="aff" rid="aff3"><sup>3</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Université Bourgogne Franche-Comté, IRTES EA7274, Belfort, France</addr-line></aff><aff id="aff2"><addr-line>Unité de Recherche sur les Agroressources et la Santé Environnementale (URASE), Université de Lomé,
Lomé, Togo</addr-line></aff><aff id="aff3"><addr-line>Département de Physique, Faculté des Sciences, Université de Ngaoundéré, Ngaoundéré, Cameroun</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>brahiman.traore@utbm.fr(RT)</email>;<email>kokou-esso.atcholi@utbm.fr(KEA)</email>;<email>odsamah@yahoo.fr(OS)</email>;<email>samuel.gomes@utbm.fr(SG)</email>;<email>bedti@yahoo.fr(TB)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>11</day><month>03</month><year>2016</year></pub-date><volume>04</volume><issue>02</issue><fpage>154</fpage><lpage>163</lpage><history><date date-type="received"><day>23</day>	<month>February</month>	<year>2016</year></date><date date-type="rev-recd"><day>accepted</day>	<month>27</month>	<year>March</year>	</date><date date-type="accepted"><day>30</day>	<month>March</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Requirements for the respect of the environment encourage to reduce the impact of human activity on the nature. Civil engineering answers these requirements by the development of ecological construction materials. This paper deals with the transformation of clay raw materials which enable the processing of environmentally friendly construction materials: in addition to their biodegradability, the alveolar fired clay materials allow energy saving in home heating thanks to their thermal isolation properties. But their manufacturing is a high energy consumption process, in particular during compaction, drying and firing which contribute to the emission of greenhouse gases. The goal of this paper is to study the rheology of clay pastes in order to develop low energy in manufacturing processes. For this purpose, theoretical and experimental approaches were carried out on six clay varieties. In the theoretical approach, a finite element (FE) simulation model has been developed for pressing a non-rigid material predicting deformations and stresses occurring within the clay structure. Experiments have then been carried out to validate the finite element modelling. In this experimental approach, the clay pastes were transformed with water content respecting the Atterberg limits which determine the plasticity of clays. The samples compaction has been carried out under variable loadings in order to determine the suitable low energy consumption loading.
 
</p></abstract><kwd-group><kwd>Simulation Model</kwd><kwd> Environmentally Friendly</kwd><kwd> Construction Materials</kwd><kwd> Rheological Behaviour</kwd><kwd> Clay Materials</kwd><kwd> Energy Saving</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Due to the recent requirements for the respect of the environment, the development of ecological construction materials has become a major challenge. The answer to this challenge is possible through the eco-design of bio- degradable materials, the development of low energy consumption processes in construction material industries. In the particular case of terra cotta materials such as tiles, bricks, floor tiles, the shaping of materials as well as their drying and firing, consume an important amount of energy. So, in spite of the biodegradability of these materials, their manufacturing doesn’t meet industrial requirements regarding the limitation of CO<sub>2</sub> and NO<sub>x</sub> emissions.</p><p>This paper contributes to understand the rheological behavior of clay raw materials when hydrated in order to optimize clay pastes shaping by lowering the compaction energy demand and by increasing manufacturing productivity. In this study, Finite Elements analyzes will help us to understand if there is a piezo-dependence of the rheological behavior in the volume of the paste and the behavior at the interface between the clay structure and the extrusion die. These numerical analyzes will be coupled with an experimental study in order to obtain the properties of pastes flowing and to study effects of the extrusion pressures through axi-symmetric dies.</p><p>The correlation of results of both approaches will allow defining a low clay material shaping pressure that will contribute to reduce the loading energy consumed in clay materials industries.</p></sec><sec id="s2"><title>2. Materials and Methods</title><sec id="s2_1"><title>2.1. Materials</title><sec id="s2_1_1"><title>2.1.1. Identification of Clay Pastes Viscosity</title><p>Clay minerals consist of two crystalline elements: a tetrahedral element (the silica, SiO<sub>2</sub>) and an octahedral element (aluminium hydroxide, Al(OH)<sub>3</sub>) [<xref ref-type="bibr" rid="scirp.65162-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref2">2</xref>] . These elements are piled up to constitute the mineral layers. Their stacking determines the variety of clays.<sub> </sub></p><p><xref ref-type="fig" rid="fig1">Figure 1</xref>(a) shows the main clay minerals according to elements stacking (kaolinite, illite and smectites).</p><p>The clay paste rheological behavior is also related to the clay elements stacking. Due to the presence of water between mineral layers, these layers are sliding and the clay paste has plastic properties (<xref ref-type="fig" rid="fig1">Figure 1</xref>(b)) [<xref ref-type="bibr" rid="scirp.65162-ref1">1</xref>] - [<xref ref-type="bibr" rid="scirp.65162-ref3">3</xref>] .</p><p>Clay pastes are not purely viscous and have a much more complex rheological behavior which is similar to that of viscoelastic fluids of Bingham with a sliding on walls [<xref ref-type="bibr" rid="scirp.65162-ref4">4</xref>] .</p><p>- Below a certain stress threshold, the viscoelastic paste does not undergo a deformation.</p><p>- Above this stress threshold, it becomes deformed and the value of the stress is related to the paste shearing or flowing velocity. This stress threshold depends on the water content of the clay paste.</p><p>Indeed, due to the foliar structure of clay minerals, the paste flows by the sliding of layers between which there is water (<xref ref-type="fig" rid="fig2">Figure 2</xref>(a)). The stress threshold depends on the thickness of water layers.</p><p><xref ref-type="fig" rid="fig2">Figure 2</xref>(b) shows a clay paste which is deformed between two plates. The down plate is fixed and the upper plate moves at a velocity V by an applied force F. The clay structure undergoes a shearing according to the de-</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Mineral structure of main clay types (a) and clay paste macroscopic structure (b)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x6.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> Deformation of clay paste by piled up layers sliding (a); flow velocity gradient in a clay paste (b)</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x7.png"/></fig><p>formation direction. The shearing stress is expressed by Equation [<xref ref-type="bibr" rid="scirp.65162-ref5">5</xref>] :</p><disp-formula id="scirp.65162-formula111"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x8.png"  xlink:type="simple"/></disp-formula><p>where S is the surface of each plate and F the force applied on the upper plate.</p><p>The gradient of the flowing velocity in the clay paste is expressed as [<xref ref-type="bibr" rid="scirp.65162-ref5">5</xref>] :</p><disp-formula id="scirp.65162-formula112"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x9.png"  xlink:type="simple"/></disp-formula><p>where H is the clay paste thickness and V, the moving velocity of the upper plate.</p><p>In these conditions, the apparent viscosity of the clay paste is given by the expression [<xref ref-type="bibr" rid="scirp.65162-ref5">5</xref>] :</p><disp-formula id="scirp.65162-formula113"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x10.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_1_2"><title>2.1.2. Identification of the Clays Used</title><p>The apparent viscosity of the studied clays has been determined experimentally (<xref ref-type="table" rid="table1">Table 1</xref>) for the need of the finite element simulations. These average values reveal a major difference in the viscosity of studied clays.</p><p>The paper deals with six clay varieties which raw materials are processed and compressed to perform con- struction materials structures density. These raw materials coming from Togo are named according to their respective deposit name and their natural colours:</p><p>- ABB: White clay of Bang&#233;li which is a kaolin;</p><p>- ANT: Black clay of Togbl&#233;kop&#233;;</p><p>- ARA: Red clay of Albi;</p><p>- ARG: Red clay of Gu&#233;rin-kouka;</p><p>- ARK: Red clay of Kouv&#233;;</p><p>- AVK: Green clay from Kouv&#233;.</p><p>Clay raw materials are generally a natural mixture of the following main constituents [<xref ref-type="bibr" rid="scirp.65162-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref6">6</xref>] : kaolin clay (which consists of aluminum oxide Al<sub>2</sub>O<sub>3</sub>, silica 2SiO<sub>2</sub>, water 2H<sub>2</sub>O), feldspar which provides sodium ion (Na<sup>+</sup>) and potassium ion (K<sup>+</sup>) that allow the formation of the vitreous phase when fired, a neutral solid constituent that reduces the shrinkage but does not act on the reactions during burning. The mineral elements contained in these clays are identified by x-rays. <xref ref-type="table" rid="table2">Table 2</xref> gives the average contents values of three of the studied clays.</p></sec></sec><sec id="s2_2"><title>2.2. Methods</title><sec id="s2_2_1"><title>2.2.1. Numerical Identification of Clay Paste Flowing Parameters</title><p>Due to their visco-plastic likely behavior under load, clay pastes can be modeled by the law [<xref ref-type="bibr" rid="scirp.65162-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref8">8</xref>] :</p><disp-formula id="scirp.65162-formula114"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x11.png"  xlink:type="simple"/></disp-formula><p>where σ is the normal stress induced by the loading; k defines the paste consistency; ε is the generalized defor- mation of the paste; <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2710385x12.png" xlink:type="simple"/></inline-formula> is the generalized deformation velocity; m is the sensitivity coefficient to the velocity; n</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Average values of apparent viscosity of studied clays pastes</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Clay variety</th><th align="center" valign="middle" >ABB</th><th align="center" valign="middle" >ANT</th><th align="center" valign="middle" >ARA</th><th align="center" valign="middle" >ARG</th><th align="center" valign="middle" >ARK</th><th align="center" valign="middle" >AVK</th></tr></thead><tr><td align="center" valign="middle" >Apparent viscosity (in poises)</td><td align="center" valign="middle" >24</td><td align="center" valign="middle" >82</td><td align="center" valign="middle" >65</td><td align="center" valign="middle" >65</td><td align="center" valign="middle" >65</td><td align="center" valign="middle" >80</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Mineral elements average contents of studied clays</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Clay variety</th><th align="center" valign="middle" >Silica (SiO<sub>2</sub>)</th><th align="center" valign="middle" >Iron oxide (Fe<sub>2</sub>O<sub>3</sub>)</th><th align="center" valign="middle" >Aluminum oxide (Al<sub>2</sub>O<sub>3</sub>)</th><th align="center" valign="middle" >Calcium oxide (CaO)</th><th align="center" valign="middle" >Potassium oxide (K<sub>2</sub>O)</th></tr></thead><tr><td align="center" valign="middle" >ABB</td><td align="center" valign="middle" >69.84</td><td align="center" valign="middle" >------</td><td align="center" valign="middle" >22.80</td><td align="center" valign="middle" >------</td><td align="center" valign="middle" >7.36</td></tr><tr><td align="center" valign="middle" >ANT</td><td align="center" valign="middle" >66.80</td><td align="center" valign="middle" >10.06</td><td align="center" valign="middle" >20.88</td><td align="center" valign="middle" >2.26</td><td align="center" valign="middle" >-----</td></tr><tr><td align="center" valign="middle" >ARA</td><td align="center" valign="middle" >64.76</td><td align="center" valign="middle" >10.58</td><td align="center" valign="middle" >20.57</td><td align="center" valign="middle" >------</td><td align="center" valign="middle" >-----</td></tr><tr><td align="center" valign="middle" >ARG</td><td align="center" valign="middle" >62.08</td><td align="center" valign="middle" >10.59</td><td align="center" valign="middle" >20.58</td><td align="center" valign="middle" >------</td><td align="center" valign="middle" >6.76</td></tr><tr><td align="center" valign="middle" >ARK</td><td align="center" valign="middle" >65.00</td><td align="center" valign="middle" >10.08</td><td align="center" valign="middle" >20.86</td><td align="center" valign="middle" >2.38</td><td align="center" valign="middle" >6.76</td></tr><tr><td align="center" valign="middle" >AVK</td><td align="center" valign="middle" >64.38</td><td align="center" valign="middle" >10.00</td><td align="center" valign="middle" >20.37</td><td align="center" valign="middle" >-----</td><td align="center" valign="middle" >7.00</td></tr></tbody></table></table-wrap><p>is the coefficient of the hardening.</p><p>By fitting the experimental curves of stress-strain σ(ε) with the law of Equation (4), the obtained parameters of this equation have the magnitude order of those typically given by the literature, namely:</p><disp-formula id="scirp.65162-formula115"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x13.png"  xlink:type="simple"/></disp-formula><p>This law is based on the properties of the class of so called viscometric flows. In this case and for these flows, the stress tensor can be written, considering the axes in cylindrical coordinates, in the following form [<xref ref-type="bibr" rid="scirp.65162-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref10">10</xref>] :</p><disp-formula id="scirp.65162-formula116"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x14.png"  xlink:type="simple"/></disp-formula><p>Then, the measurement of the strain time rate permits to define three functions called viscometric functions and that define the behavior of the clay pastes [<xref ref-type="bibr" rid="scirp.65162-ref9">9</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref10">10</xref>] :</p><disp-formula id="scirp.65162-formula117"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x15.png"  xlink:type="simple"/></disp-formula></sec><sec id="s2_2_2"><title>2.2.2. Mathematical Model of Clay Pastes Flowing</title><p>Due to the cylindrical shape of the molded clay samples, the equations are defined according to the cylindrical coordinates (r, q, z). <xref ref-type="fig" rid="fig3">Figure 3</xref> illustrates these geometrical considerations. Due to the symmetrical and revolution shape, the non null stresses applied to a cylindrical element are s<sub>rr</sub>, s<sub>qq</sub>, s<sub>zz</sub>, t<sub>rz</sub> and t<sub>zr</sub> (<xref ref-type="fig" rid="fig3">Figure 3</xref>(b)).</p><p>In this analysis of clay pastes flowing, we neglect the terms related to gravity, because of the high rigidity of the pastes when completely filled in the shaping cylinder. In the same way we consider only the case of sufficiently slow flows in order to neglect the inertial terms.</p><p>Under these conditions, the balance equations are written as following [<xref ref-type="bibr" rid="scirp.65162-ref6">6</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref11">11</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref12">12</xref>] :</p><disp-formula id="scirp.65162-formula118"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x16.png"  xlink:type="simple"/></disp-formula><p>where u<sub>r</sub>, u<sub>q</sub> and u<sub>z</sub> are the velocity components according to the radial, tangent and axial directions, and p the loading pressure during the paste compaction.</p><p>With the assumption that the clay paste is incompressible, the continuity equation is expressed as follows:</p><disp-formula id="scirp.65162-formula119"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x17.png"  xlink:type="simple"/></disp-formula><fig-group id="fig3"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Compaction of clay paste in two parallel plates (a) and stresses on a cylindrical element (b).</title></caption><fig id ="fig3_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x18.png"/></fig></fig-group><p>It is assumed that the law governing the fluid behavior may be written as:</p><disp-formula id="scirp.65162-formula120"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x19.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/5-2710385x20.png" xlink:type="simple"/></inline-formula> is a function of invariants of the strain time rate tensor that is independent of the deformation history.</p><p>The compression were also carried out on dried samples, in the condition of a low rate of deformation (about 0, 5 mm/min, to avoid an abrupt fracture of the sample) in order to locate and follow the materials cracking.</p></sec><sec id="s2_2_3"><title>2.2.3. Finite Element Modeling of a Clay Paste Compression</title><p>The problem is axi-symmetric due to the samples cylindrical shape. The clay structure is supposed to be iso- tropic when compressed at a low displacement speed. The paste compression is simulated with the code ANSYS. To assess the efficiency of the simulation compared to the experimental approach, a numerical three dimension- al example was studied on a cylindrical shape sample loaded by a compression force F<sub>z</sub> equal to 5 kN and 35 kN or a uniform pressure respectively of 7 MPa and 50 MPa. A 30 mm diameter sample is meshed with a 3D solid triangular element. The bottom of the sample is clamped and its top is considered as a rigid surface [<xref ref-type="bibr" rid="scirp.65162-ref13">13</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref14">14</xref>] . The clay sample is compressed in a hollow solid cylinder which surfaces are clamped. <xref ref-type="fig" rid="fig4">Figure 4</xref>(a) shows the meshed sample under compression loading in the hallow cylinder and <xref ref-type="fig" rid="fig4">Figure 4</xref>(b) shows the corresponding ex- perimental device.</p><p>The calculation code doesn’t allow simulating the compression of viscoelasic pastes. In the present case, the simulation concerns the compression of clay structure that is considered to be an elasto-plastic and isotropic material. The modelling parameters are the Young modulus and the coefficient of Poisson associated with an asymmetric option and an incremental loading.</p></sec><sec id="s2_2_4"><title>2.2.4. Samples Elaboration and Compression Tests Procedure</title><p>The raw material aggregates size plays an important role in the mechanical properties of the obtained paste. A paste with smaller particles sizes has a more high plasticity because the particles react more intensively between each other. Thus, studied clay raw materials were grinded. The resulted powders were kept at 60˚C during 24 hours in an oven until a complete dehydration. For the samples characterization, each clay paste is processed by mixing a mass m<sub>0</sub> of anhydrous clay powder with a mass me of water. The resulting clay paste has an average water content v expressed as:</p><disp-formula id="scirp.65162-formula121"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x21.png"  xlink:type="simple"/></disp-formula><p>where m<sub>ws</sub> is the mass of the wetted sample and m<sub>0s</sub> the mass of the same sample when dried.</p><p>The masses were measured using a METTLE PJ 360 Delta Range balance, with an accuracy of 1/100 g.</p><p>For samples processing, an average water content v of about 18% was chosen in relation with the assessment of pastes fluidity, the consistency of pastes and the external appearance of samples after shaping. The water content</p><fig-group id="fig4"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> A clay paste sample loaded by a uniform pressure in a cylindrical rigid container (a); corresponding device for the experimental shaping of samples by compression (b).</title></caption><fig id ="fig4_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x22.png"/></fig><fig id ="fig4_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x23.png"/></fig></fig-group><p>of the six varieties of studied clays ranges from 15% to 20%. These water content values lie in the limits set by ATTERBERG [<xref ref-type="bibr" rid="scirp.65162-ref12">12</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref15">15</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref16">16</xref>] .</p><p>To get a clay paste having satisfactory plasticity properties, the powder-water mixtures were kneaded for about 2 hours until getting a homogeneous paste. The paste is then stored in a hermetically climatic chamber for more than 24 hours, to insure a uniform moisture content while preventing evaporation and increasing the plas- ticity of the paste under the effect of micro-organisms that play an important role in the process.</p><p>Clay pastes were compacted by compression in a duralumin hallow cylinder shown in <xref ref-type="fig" rid="fig4">Figure 4</xref>(b). This cy- linder has an interior diameter of about f30 mm. The loading is made by a f30 piston. To meet the manufactur- ing requirements in tiles industrial factories, the clay structure must be compressed to present a density of about 2 g/cm<sup>3</sup>. For this purpose, to obtain samples having a height of 30 mm and a diameter of 30 mm, the mass of the paste to be introduced into the mould is calculated so to obtain a structure density of 2 g/cm<sup>3</sup> after compression.</p><p>Compression is performed by imposing the piston motion at a low speed of 3 mm/min to avoid resistance re- lated to the viscosity of the paste. Clay pastes are subject to adhesion to the die walls in industrial processes, and this is also the case for the moulding carried out in this study. In most cases, friction is below the threshold at which shear occurs. Therefore, pastes slip instead of warping. This promotes good extrusion without lubricants [<xref ref-type="bibr" rid="scirp.65162-ref17">17</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref18">18</xref>] . Several samples were processed (about 24 samples). Samples were loaded at 5 kN, 20 kN and 35 kN or respectively at compression stress of about 7 MPa, 28 MPa and 50 MPa.</p></sec></sec></sec><sec id="s3"><title>3. Results and Discussion</title><sec id="s3_1"><title>3.1. Numerical Simulation</title><p><xref ref-type="fig" rid="fig5">Figure 5</xref>(a) and <xref ref-type="fig" rid="fig5">Figure 5</xref>(b) give respectively the displacements field and the radial stresses field in the cylin- drical clay structure loaded at 7 MPa (corresponding to a compression force of 5 kN).</p><p><xref ref-type="fig" rid="fig5">Figure 5</xref>(a) shows that the displacements field varies according to the loading direction. On the other hand, the distribution of the axial stresses (<xref ref-type="fig" rid="fig5">Figure 5</xref>(b)) is not uniform. These results show that the displacement and the radial stress decrease according to the loading direction, i.e., from the top to the bottom of the sample. This gradient of displacements and radial stresses distribution can be partially explained by shear stresses on the in- terface between the clay structure and the solid die cylinder resulting in an important friction on this interface. This friction generates micro-cracks that will widen during drying. The non uniform distribution of the dis- placements and the radial stresses generate also a non uniform drying shrinkage, resulting also in the structure cracking [<xref ref-type="bibr" rid="scirp.65162-ref2">2</xref>] [<xref ref-type="bibr" rid="scirp.65162-ref19">19</xref>] .</p><p>To verify the assumption of the influence of the interfacial friction on the non uniform distribution of the radial stresses and displacements, we simulate the case of a parallelepipedic shape sample. This block has the same loading face surface as the previous cylindrical sample. The loading force is increased up to 35 kN, i.e., a compression pressure of 50 MPa. <xref ref-type="fig" rid="fig6">Figure 6</xref>(a) and <xref ref-type="fig" rid="fig6">Figure 6</xref>(b) give respectively the displacements field and the radial stresses field in the parallelepipedic clay structure loaded at 50 MPa.</p><p>The results of <xref ref-type="fig" rid="fig6">Figure 6</xref>(a) and <xref ref-type="fig" rid="fig6">Figure 6</xref>(b) show that, near the bottom of the sample, the displacements and</p><fig-group id="fig5"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Displacement field (a) and radial stress field (b) in a cylindrical clay structure loaded at 7 MPa (5 kN).</title></caption><fig id ="fig5_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x24.png"/></fig><fig id ="fig5_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x25.png"/></fig></fig-group><fig-group id="fig6"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Displacement field (a) and radial stress field (b) in a clay block structure loaded at 50 MPa (35 kN).</title></caption><fig id ="fig6_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x27.png"/></fig><fig id ="fig6_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x26.png"/></fig></fig-group><p>the radial stresses are more or less important in the corners of the parallelepipedic sample, thus confirming the assumption of a gradient of displacements and stresses due to the friction of the clay structure on the solid die walls. Due to the concentration of displacements and stresses in the corners of the samples during the structure drying, the shrinkage on the corners differs from the shrinkage at the other points of the clay structure. This in- duces the structure cracking during drying, cracks starting by the corners.</p><p>Compared to <xref ref-type="fig" rid="fig5">Figure 5</xref>(b), the result of <xref ref-type="fig" rid="fig6">Figure 6</xref>(b) shows that the increase of compression loading lowers the non uniform distribution of the radial stresses. This reduces the structure cracking by shrinkage during drying.</p><p>Finally, following figures which correspond to the radial stress field of samples loaded respectively at 7 MPa (<xref ref-type="fig" rid="fig7">Figure 7</xref>(a)) and at 50 MPa (<xref ref-type="fig" rid="fig7">Figure 7</xref>(b)) show that the increase of the loading pressure decreases the size of the low stress zones (in blue) in the clay structure volume.</p></sec><sec id="s3_2"><title>3.2. Sample Compression</title><p>The compression is performed with a conventional hydraulic traction-compression machine. The sample defor- mation and the loading force are recorded by computer devices during loading. Thus, the loading force Fz can</p><fig-group id="fig7"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Influence of the die shape and the loading pressure on the internal stress field: cylindrical shape sample loaded at 7 MPa (a) and parallelepipedic shape sample loaded at 50 MPa (b).</title></caption><fig id ="fig7_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x29.png"/></fig><fig id ="fig7_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x28.png"/></fig></fig-group><p>be plotted in function of the loading piston displacement Dz.</p><p><xref ref-type="fig" rid="fig8">Figure 8</xref>(a) shows the curves of the loading force in function of the displacement in the case of samples of AVTK clay loaded at 5 kN, 20 kN and 35 kN. All studied clays have a similar mechanical behavior as curves shown by <xref ref-type="fig" rid="fig8">Figure 8</xref>(a). <xref ref-type="fig" rid="fig8">Figure 8</xref>(b) is the superimposed curves of studied clays loaded at 35 kN.</p><p>Referring to <xref ref-type="fig" rid="fig8">Figure 8</xref>(b), one observes that, if all curves have the same global evolution, they differ clearly:</p><p>- The most flat curve corresponds to the AVK clay which presents the most plasticity aspect;</p><p>- One can introduce a shape parameter to characterize Fz(Dz) curves, for instance the spreading factor &#229; of the curve at a given loading force value. The spreading factors of curves in <xref ref-type="fig" rid="fig8">Figure 8</xref>(b) can be compared as follow;</p><disp-formula id="scirp.65162-formula122"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/5-2710385x30.png"  xlink:type="simple"/></disp-formula><p>- The spreading factor decreases from the most plastic clay variety to less plastic one. The compaction rate of clays let suppose that the more plastic a clay paste is, the more resistant is its structure. This means that the spreading factor decreases from the most resistant clay structure to the less resistant one.</p><p>To verify this relation between a clay plasticity properties and its dried structure resistance, the compression tests were carried out on dried clay samples in order to determine their mechanical resistance. <xref ref-type="fig" rid="fig9">Figure 9</xref>(a) gives the stress-strain σ(ε) curves of AVTK dried clay samples previously shaped at 5 kN, 20 kN and 35 kN and <xref ref-type="fig" rid="fig9">Figure 9</xref>(b) gives superimposed stress-strain curves of dried samples of studied clays shaped at 35 kN.</p><p><xref ref-type="fig" rid="fig9">Figure 9</xref>(a) shows that the mechanical resistance of the dried clay structures increases with the shaping loading force. But samples shaped at 20 kN and 35 kN have the same yield stress. This means that there is a shaping loading force limit over which the loading force does not increase the clay structure dry resistance. By considering the shaping loading force curves of <xref ref-type="fig" rid="fig8">Figure 8</xref>(a), one can estimate this limit around 10 kN or 14 MPa because above this value, the loading force increases linearly in function of the displacement. The deformation of the structure over this limit is no more plastic but it becomes reversible. The compaction of solid particles in the clay structure is completed when loading reaches 14 MPa. So, to shape clay based materials, the shaping loading force should not exceed 14 MPa. This allows saving the energy related to overloading.</p><p>The stress-strain curves of <xref ref-type="fig" rid="fig9">Figure 9</xref>(b) confirm partially the relation between the clay paste plasticity and the mechanical resistance of resulted dried structures. Except of the AVTK clay variety, the studied clays verify this ruler. But one should not forget that these samples must be fired in order to acquire a high mechanical resistance.</p></sec></sec><sec id="s4"><title>4. Conclusions</title><p>From the results obtained, it can be concluded as follows.</p><fig-group id="fig8"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> Compaction curves of AVTK clay paste loaded at 5 kN, 20 kN and 35 kN (a) and super imposition of compaction curves of studied clay pastes loaded at 35 kN (b).</title></caption><fig id ="fig8_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x32.png"/></fig><fig id ="fig8_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x31.png"/></fig></fig-group><fig-group id="fig9"><label><xref ref-type="fig" rid="fig9">Figure 9</xref></label><caption><title> Stress-strain curves of dried clay samples: case of AVTK clay samples shaped at 5 kN, 20 kN and 35 kN (a) and case of studied clays shaped at 35 kN (b).</title></caption><fig id ="fig9_1"><label> (b)</label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x33.png"/></fig><fig id ="fig9_2"><label></label><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/5-2710385x34.png"/></fig></fig-group><p>The FE simulation giving the radial stress field in the clay structure loaded at 7 MPa or 5 kN and at 50 MPa or 35 kN is in accordance with the experimental loading curves. Indeed, when loaded at 7 MPa or 5 kN, the clay structure does not yet reach its completed compaction. But when loaded at 50 MPa or 35 kN, the structure compaction is completed and the radial stresses in the major part of the structure reach the loading pressure.</p><p>The shaping loading force was estimated around 10 kN or 14 MPa. Beyond this value the force of compaction has no effect even if it increases. Knowing this value, we can save the energy related to overloading.</p><p>But at this stage, it is difficult to establish a relation between the stresses field distribution and the clay structure cracking during drying. An improved FE model should allow predicting structure cracks occurring during drying.</p><p>The modeling of the present paper will be completed by simulations predicting thermo-chemical transformations occurring in clay structures during burning and inducing the structure hardening.</p></sec><sec id="s5"><title>Cite this paper</title><p>Brahiman Traore,1 1,Kokou Esso Atcholi,Ouro-Djobo Samah,Samuel Gomes,Tibi Beda, (2016) Numerical Simulation and Experimental Characterization of Clay Paste under Loads for Energy Saving in Clay Materials Processing. Journal of Minerals and Materials Characterization and Engineering,04,154-163. doi: 10.4236/jmmce.2016.42015</p></sec></body><back><ref-list><title>References</title><ref id="scirp.65162-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Caillère, S., Henin, S. and Rautureau, M. (1982) Minéralogie des argiles. INRA, Actualités Scientifique, Masson.</mixed-citation></ref><ref id="scirp.65162-ref2"><label>2</label><mixed-citation publication-type="book" xlink:type="simple">Jouenne, C.A. (1990) Traité de céramiques et de matériaux minéraux. Ed. Septima, Paris.</mixed-citation></ref><ref id="scirp.65162-ref3"><label>3</label><mixed-citation publication-type="book" xlink:type="simple">Barthelemy, B. (1980) Notions pratiques de la mécanique de la rupture. Ed. Eyrolles, p 290.</mixed-citation></ref><ref id="scirp.65162-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Roquet, N. and Saramito, P. (2003) An Adaptive Finite Element Method for Bingham Fluid Flows around a Cylinder. Computer Methods in Applied Mechanics and Engineering, 192, 3317-3341.  
&lt;br /&gt;http://dx.doi.org/10.1016/S0045-7825(03)00262-7</mixed-citation></ref><ref id="scirp.65162-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Ragouilliaux, A. (2007) Etude rhéophysique de systèmes émulsions inverses/argile organophile. Applications aux boues de forage pétrolier, PhD of Université Pierre et Marie Curie, France.</mixed-citation></ref><ref id="scirp.65162-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Costet, J. (1983) Cours pratique de mécanique des sols: Plasticité et calcul des tassements. 3rd Edition, Vol. 1, Dunod, 26-41.</mixed-citation></ref><ref id="scirp.65162-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Ashby, M.F. and Jones, D.R.H. (1991) Microstructure et mise en oeuvre. Dunod, 2, 149-194.</mixed-citation></ref><ref id="scirp.65162-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Chiarelli, A.-S., Shao, J.-F. and Hoteit, N. (2003) Modelling of Elastoplastic Damage Behaviour of a Clay Stone. International Journal of Plasticity, 19, 23-45. &lt;br /&gt;http://dx.doi.org/10.1016/S0749-6419(01)00017-1</mixed-citation></ref><ref id="scirp.65162-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Conil, N., Djeran-Maigre, I., Cabrillac, R. and Su, K. (2004) Poroplastic Damage Model for Claystones. Applied Clay Science, 26, 473-487. http://dx.doi.org/10.1016/j.clay.2003.12.019</mixed-citation></ref><ref id="scirp.65162-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Hirata, S., Yao, S. and Nishida, K. (1990) Multiple Regression Analysis between the Mechanical and Physical Properties of Cohesive Soils. Soils and Foundations, 30, 91-108. &lt;br /&gt;http://dx.doi.org/10.3208/sandf1972.30.3_91</mixed-citation></ref><ref id="scirp.65162-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Carretero, M.I., Dondib, M., Fabbrib, B. and Raimondob, M. (2002) The Influence of Shaping and Firing Technologies on Ceramic Properties of Calcareous and Non-Calcareous Illitic-Chloritic Clays. Applied Clay Science, 20, 301-306.  
http://dx.doi.org/10.1016/S0169-1317(01)00076-X</mixed-citation></ref><ref id="scirp.65162-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Couarraze, G. and Grossiord, J.L. (1991) Initiation à la Rhéologie. Les dernières nouveautés Tec &amp; Doc, 51-91, 5-202.</mixed-citation></ref><ref id="scirp.65162-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Feng, Z.-Q., Magnain, B. and Cros, J.-M. (2006) Solution of Large Deformation Impact Problems with Friction between Blatz-Ko Hyperelastic Bodies. International Journal of Engineering Science, 44, 113-126. 
&lt;br /&gt;http://dx.doi.org/10.1016/j.ijengsci.2005.06.006</mixed-citation></ref><ref id="scirp.65162-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Peyraut, F., Feng, Z.-Q. and Labed, N.(2007) A Material-Independent Algorithm for Preserving of the Orientation of the Spatial Basis Attached to Deforming Medium. Computational Mechanics, 40, 1053-1060. 
http://dx.doi.org/10.1007/s00466-007-0163-0</mixed-citation></ref><ref id="scirp.65162-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Lambe, T.W. and Whitman, R.V. (1979) Soil Mechanics. SI Version, Wiley, New York.</mixed-citation></ref><ref id="scirp.65162-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Yalcin, A. (2007) The Effects of Clay on Landslides: A Case Study. Applied Clay Science, 38, 77-85. 
&lt;br /&gt;http://dx.doi.org/10.1016/j.clay.2007.01.007</mixed-citation></ref><ref id="scirp.65162-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Druyanov, B.A. and Nepershin, R.I. (1994) Problem of Technological Plasticity. Studies in applied Mechanics, 38, 357-389.</mixed-citation></ref><ref id="scirp.65162-ref18"><label>18</label><mixed-citation publication-type="book" xlink:type="simple">Laroze, S. and Barrau, J.-J. (1995) Mécanique des structures solides élastiques, plaques et coques. Ed., Ensae, Col. Sup’Aéro, Vol. 1, bis.</mixed-citation></ref><ref id="scirp.65162-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Jomaa, W. (1991) Séchage des matériaux fortement déformables: Prise en compte de la vitesse de retrait. Thèse No. 91 BOR 10607, Bordeaux.</mixed-citation></ref></ref-list></back></article>