Microsoft Word - numero_52_art_15_2711 A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 181 Finite element modeling of the behavior of polymethyl-methacrylate (PMMA) during high pressure torsion process Ahmed Drai Mechanical Engineering Department, University of Mascara, 29000, LABAB Laboratory of ENPO, Algeria. draiahmed14@yahoo.fr, http://orcid.org/ 0000-0001-6091-9519 Benaoumeur Aour Laboratory of Applied Biomechanics And Biomaterials (LABAB), ENP Oran, BP1523 El Mnaour, Oran, 31000, Algeria. benaoumeur.aour@enp-oran.dz, http://orcid.org/0000-0002-2345-6790 Naima Belayachi Université d’Orléans, Laboratoire de Mécanique Gabriel Lamé Polytech Orléans, 45072 Orléans, France. naima.belayachi@univ-orleans.fr, http://orcid.org/0000-0002-2345-6791 Abderrahim Talha Ecole des Hautes Etudes d’Ingénieur (HEI), 13 rue de Toul, 59046 Lille, Lille Mechanics Unit, University of Lille, France abderrahim.talha@hei.fr, http://orcid.org/0000-0003-2345-6792 Noureddine Benseddiq Lille Mechanics Unit, University of Lille, France noureddine.benseddiq@univ-lille1.fr, http://orcid.org/0000-0003-2345-6793 ABSTRACT. High Pressure Torsion (HPT) is a highly effective super-plastic deformation process for obtaining nano-materials with high performance mechanical properties. In view of its optimization, it is of paramount importance to evaluate the relations between the behavior of the material under the effects of different processing parameters. In this context, this work aims to highlight the plastic strain distribution in the deformed material as a function of the hydrostatic pressure, the torsion angle and the temperature of the material applied during the process. A typical amorphous polymer (Polymethyl-Methacrylate: PMMA) has been tested. Firstly, in order to identify the material parameters of a phenomenological elasto-viscoplastic model compression tests at different temperatures and strain rates have been carried out. Then, the distributions of the effective plastic strain, the equivalent plastic strain rate and the hydrostatic stress were analyzed using finite elements method. Recommendations on process conditions were Citation: Drai, A., Aour, B., Belayachi, N., Talha, A., Benseddiq, N., Numerical modeling of polymethyl-methacrylate (PMMA)behavior during high pressure torsion process, Frattura ed Integrità Strutturale, 52 (2020) 181-196. Received: 20.12.2019 Accepted: 10.02.2020 Published: 01.04.2020 Copyright: © 2020 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. https://youtu.be/30UA-qBjqwY A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 182 proclaimed at the end of this work according to the obtained numerical results. KEYWORDS. HPT; Finite Elements; PMMA; Behavior; Plastic strain. INTRODUCTION here have been extensive investigations over the last two decades into the production of bulk ultrafine-grained (UFG) materials through the application of severe plastic deformation (SPD) [1, 2]. Recently, SPD techniques have been used in the modification of microstructure. They made possible to produce nano-crystalline (NC) microstructures out of metallic materials [3, 4]. Such SPD methods include High Pressure Torsion (HPT), Equal Channel Angular Pressing (ECAP), Accumulative Roll Bonding (ARB), Multiple Forging, Twist Extrusion (TE) and some others. Among the various SPD techniques, Equal-Channel Angular Pressing (ECAP) and High Pressure Torsion (HPT) are the most common [5]. To date, processing by HPT has proven to be the most effective of all the SPD methods in producing bulk nanostructured materials [6]. Its principle consists in putting a disk between two massive anvils under a simultaneously or successively action of high compression and torsion as shown in Fig. 1. The occurrence of superplasticity was reported not only after HPT processing of thin disc specimens [7,8], but also using the ring samples [9,10], bulk samples [11, 12] as well as after continuous high-pressure torsion extrusion [13]. Different materials have been tested experimentally [14] and theoretically [15] through the application of HPT process such as aluminum alloys [8, 16], nickel [17], copper [15], titanium alloys [3, 18], polymers [19], zirconium [20, 21], Magnesium alloys [22], aluminum-zinc alloys [23], intermetallics [24, 25] and others (see [26, 27]). Figure 1: Schematic illustration of HPT process. In the past, FEM simulations of the SPD processes were performed with the purpose of analysis and optimization of the SPD processing variables and prediction of the mechanical and microstructural characteristics of the SPD processed samples [15]. For example, for HPT, the effects of processing parameters, such as geometries of the disc [28], applied pressure [29], temperature distribution [30, 31], rotation rate [31], friction conditions [32], material porosity and properties [32, 33] and slope of the inclined anvil's surface [34] on the plastic strain behavior were reported in the literature extensively. This work is devoted to study the behavior of polymethyl-methacrylate (PMMA) during HPT process. It is important to note that the behavior of this amorphos polymer is sensitive to the strain rate, temperature, hydrostatic pressure and torsion angle [35]. A three-dimensional finite element analysis has been conducted using hexahedral element with 8-nodes. Our particular attention has been focused on the distribution of the plastic strain in the deformed material as a function of the different parameters of the process, namely, the sequence of the different phases, the hydrostatic pressure, the torsion angle and the temperature of the material during the process. T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 183 CONSTITUTIVE MODELING o model the PMMA behavior during HPT process, the elasto-viscoplastic model of Perzyna has been used [36]. The material parameters of this model have been identified using compressive tests on PMMA cylindrical specimens at different temperatures and strain rates. The experimental protocol used for this solicitation mode was also presented in details in the next section. The Perzyna model considers the hypo-elastic relation, and the additive decomposition of the tensor of the strain rates to write [36]: . 4σ H : (D D )vp  (1) with H is the Hooke tensor (tensor of elasticity) given by: 2 H ( ) 2(1 ) 1 2 ik jl il jk ij kl E                  (2) where E and represent respectively the Young's modulus and the Poisson's ratio and  is the Kronecker symbol. The viscoplastic deformation rate is given by the rule of normality. Its direction is normal to the surface flow in the stress space f    , f being the surface flow. By introducing the hypothesis that the function f is independent of the deformations, and supposing the Von Mises criterion [37], the viscoplastic strain is formulated as follows: D φ( )vp f f     (3) with f    = 1 in the case of a perfectly viscoplastic material. where φ(f) is a function generally chosen as a power function: 1 0 ( ) 1 m f           (4) Then, we obtain: 1 0 D 1 m vp          (5) where m and  are respectively the hardening parameter and the viscosity (strain rate sensitivity parameters),  is the flow stress and 0 is the static elastic limit. The Perzyna model describes only the sensitivity of the yield stress to the strain rate, then it is necessary to introduce the post yield softening and hardening. The finite element program used allows us to take into account this post-flow behavior. A multilinear isotropic hardening combined with the Von Mises criterion has been adopted. It should be noted that in this study, we assume that the viscoelastic part of the response does not significantly affect the global results and the pre-yield response is assumed to be linear and elastic. The behavior in large deformations of viscoplastic solids is described by the constitutive equations, for which analytical solutions are difficult to obtain, even for simplified problems. This model has attracted a lot of attention in the literature because it was validated experimentally in the work of Perzyna [38]. In addition, its simplicity of writing made it easier to T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 184 implement in simulation codes. Peric et al. [39] have used this model in their numerical studies to adapt to large deformations the behavior of materials depending on time and the history of deformation. The work of Van Der Sluis et al. [40, 41] on polycarbonate loaded with elastomer particles has brought an interesting approach based on the Perzyna model and a nonlinear hardening which has shown the efficiency of the use of this models type for polymers. EXPERIMENTAL PROCEDURE FOR MATERIAL PARAMETERS IDENTIFICATION o identify the parameters of the elasto-viscoplastic constitutive law presented in the previous paragraph, we used compression tests at different strain rates and different temperatures conducted on polymethyl-methacrylate (PMMA) samples. Noting that compressive loading, unlike traction, provides (or delays as much as possible) the damage mechanisms. Material properties The material studied (PMMA) was supplied by the Goodfellow© company in the form of a cylindrical bar of 8 mm in diameter and 1000 mm in long, with a molar mass of the order of 65 kg.mol-1. The glass transition temperature is about 120°C. The mechanical properties of this material are shown in Tab. 1. Mechanical properties of PMMA Density (g cm-3) 1.19 Coefficient of friction 0.25 – 0.4 Hardness - Rockwell 92-100 Impact resistance - Izod (J m-1) 16-32 Poisson coefficient 0.35 – 0.4 Elongation at fracture (%) 2.5-4 Traction module (GPa) 2.4-3.3 Tensile strength (MPa) 80 Table 1: Mechanical properties of the studied polymer (PMMA). Description of the mechanical test In order to determine the three parameters ( 0 , m,  ) of the elasto-viscoplastic model, the PMMA sample was tested at strain rates ranging from 10-5 s-1 to 10-2 s-1 at a room temperature (25°C) and high temperatures (40°C, 60°C and 80°C). The uniaxial compression tests are performed on an Instron 5867 tensile/compression machine using a 10 KN cell to measure the applied force. The acquisition of experimental data in stresses and strains at each moment during the test was recorded using the Bluehill software. The displacement of the upper plate is controlled with a translational speed of the crossbar calculated from the speed of deformation. All samples were compressed up to 40% of deformation using two parallel compression trays. To study the dependence of the polymer behavior at a high temperature, we used a thermal enclosure allowing tests with temperatures up to 250°C. To ensure the homogeneity of the temperature in the specimen, it was maintained for about 15 minutes at the desired temperature before starting the test. The upper end of the sample is left free during the heating phase to allow thermal expansion to take place freely. To investigate the influence of the temperature increase on the stress-strain curve, the PMMA sample was tested at temperatures ranging from 25°C to 100°C at a constant strain rate. From the registration of the cross-bar displacement and the measurement of the force F during the test, the nominal strain- stress curves are deduced using the Bluehill software integrated in the electromechanical testing machine. To make the necessary transformations in true quantities (true stress, true strain), it is assumed that the material is incompressible and isotropic, and the deformation is supposed to be homogeneous. The extensometer allows to measure the displacement ΔL and by knowing the initial height of the sample L0, the true axial strain is obtained from the following expression: 0 1zz L Ln L         (6) The true stress becomes: T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 185 0 . 1zz N L L          (7) where 0N F S  , is the nominal stress and S0 is the initial section. Fig. 2 shows the axial stress-strain curves resulting from the tests at different temperatures and a speed of 10-3 s-1. From this graph, we can notice that the stress level decreases when the temperature increases. In addition, the curves confirm the very marked nonlinearity of the PMMA behavior, and a weak softening followed by a noticeably observed hardening. Figure 2: True axial stress-strain curves at different temperatures and a strain rate of 10-3 s-1 on a PMMA samples. By plotting on a graph the evolution of the stress as a function of the deformation (Fig. 3), for different strain rates, and for the same temperature, we have been able to highlight the dependence of PMMA behavior on strain rate. These curves show that the elastic limit, the Young's modulus, and the hardening decreases with the decrease of the strain rate. The curves show a peak which is slightly pronounced at room temperature and becomes more significant with the increase of the temperature. PARAMETERS IDENTIFICATION AND VALIDATION OF THE MODEL he identification of viscoplastic parameters for PMMA goes through the classical process using linear least squares regression from PMMA strain-strain curves at different strain rates. We draw the following function: 1 D 1 0 mvp Ln f                     (8) For the three strain rate values 10-4, 10-3, and 10-2 s-1, this evolution is a linear function y=a.x+b which allows us to determine the coefficients a and b. Finally the parameters are deduced as follows: γ = Exp(b) , m = 1/a (9) The values of the parameters obtained for each temperature are given in Tab. 2. T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 186 Parameter T=25°C T=60°C T=80°C E (MPa) 2565 1084 833  0.32 0.33 0.33 m 0.346 0.294 0.295  0.027 0.0137 0.0298 0 (MPa) 89 40 29 Table 2: Material parameters of elastic-viscoplastic model for PMMA at different temperatures. To validate the viscoplastic model, we simulated compression tests performed on PMMA sample at different strain rates and at study temperatures {T = 25, 60 and 80°C}. We have examined the macroscopic behavior of PMMA with the set of parameters reported in Tab. 2. As can be seen in Fig. 3, there is an acceptable agreement between experimental data and simulation results. For high strain rates and at 25°C and 80°C, the predicted flow yield is near to the experimental measurements, but the hardening deviates from it. This difference can be explained by the experimental calculation of the stresses from a rectangular surface, whereas the sample has a barrel shape with a section which changes and becomes more important in the middle. (a) T=25°C (b) T=60°C (c) 80°C Figure 3: Stress-strain curves of the PMMA for different strain rates and at temperatures of: (a) 25°C, (b) 60°C and (c) 80°C. A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 187 SIMULATION OF THE HPT PROCESS he simulations were performed by the finite element code MSC.Marc using eight (8) nodes hexahedral elements. The initial dimensions of the PMMA samples used are 20 mm in diameter and 10 mm in thickness. It should be noted that the mesh size selected is largely sufficient to show accurately the distribution of localized plastic strain within the samples. The upper and lower anvils were considered as rigid bodies in the finite element simulations. The upper anvil is subjected to an imposed displacement of compression, while the lower anvil is subjected to rotation. The effects of the hydrostatic pressure applied and the torsion angle on the evolution of the equivalent plastic strain were highlighted in this numerical investigation of PMMA behavior during HPT process. Effect of the hydrostatic pressure To highlight the effect of hydrostatic compression, four different vertical displacements were imposed by the upper anvil. A constant rotation speed of 0.2618 rad/sec was applied during the torsion. The results obtained for the evolution of the equivalent plastic strain along the radial distance from the center for various imposed compression displacements combined with a 30° of torsion angle are illustrated in Fig. 4. It may be noted that the equivalent plastic strain increases with the increase of the imposed compressive displacement. In addition, it can be seen that the influence of this parameter becomes significant by moving away from the center of the specimen and likewise for the heterogeneity of the plastic strain distribution. Indeed, for an imposed displacement of 1 mm, an equivalent plastic strain (εp) of 1.49 was obtained, whereas, for the one of 2.5 mm, εp=1.91. Figure 4: Distribution of the equivalent plastic strain along the radial distance of the disk for different imposed vertical displacements. The distributions of the equivalent plastic strain in the deformed geometries of the sample at the end of the process for imposed displacements of 1mm, 1.5mm, 2mm and 2.5mm are illustrated in Fig. 5. It can be noted that the highest deformation values are located in the lower and upper regions of the sample. Therefore, in order to improve the level of the plastic strain distribution, it is necessary to increase the torsion angle. Furthermore, it should be interesting to remember that the decrease of the sample thickness contributes also to the increase of the plastic strain as it has been highlighted by Drai and Aour [42]. Effect of torsion angle In order to illustrate the effect of torsion angle on the evolution of equivalent plastic strain a vertical displacement of 1 mm was imposed by the upper anvil on the sample and was maintained during the torsion with different angles of the lower anvil (15°, 30°, 45° and 60°). The simulations were performed with low angular velocity of 0.26 rad/sec in order to conduct the process in isothermal conditions. T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 188 (a) (b) (c) (d) Figure 5: Distribution of the equivalent plastic strain as a function of the imposed displacements on the deformed PMMA sample at the end of HPT process in the case of: (a) 1 mm, (b) 1.5 mm, (c) 2 mm and (d) 2.5 mm. Evolution of equivalent plastic strain Fig. 6 shows the distribution of the equivalent plastic strain along the sample radius for an imposed displacement of 1 mm and different torsion angles. It can be observed that the accumulated plastic deformation increases with the increase of the torsion angle and becomes very important at the edge of the sample. Fig. 7 illustrates the isovalue contours of the equivalent plastic strain distribution in the PMMA samples deformed during the HPT process, in the case of 1mm imposed displacement and a torsion angle of: (a) 15°, (b) 30°, (c) 45° and (d) 60°. From Fig. 8, it can be observed that after the vertical compressive displacement of 1 mm the plastic strain is very low and its maximum value does not exceed 0.3. In addition, at this stage, the plastic strain induced by a simple compression remains very low especially in the middle of the sample. However, after the application of torsion under high pressure, the plastic strain increases with the increase of the torsion angle  and becomes very high in the upper and lower areas, where it reaches a maximum value of 4.72 when  = 60°. On the other hand, in the center it is always weak. Therefore, it is advised to use small A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 189 thicknesses in the case of PMMA samples, to obtain a good homogeneity and a large plastic strain in the sample processed by HPT. Figure 6: Evolution of the equivalent plastic strain as a function of the radial distance for different torsion angles at the end of HPT process. Evolution of the mean normal stress In order to study the state of the stresses in the PMMA sample, the distributions of mean normal stresses along the diameter of the upper face at the end of HPT process using an imposed compressive displacement of 1mm of compression with different torsion angles are presented in Fig. 9. It can be seen that as the torsion angle increases, the compression stress increases at the central portion of the upper surface of the sample. (a) (b) A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 190 (c) (d) Figure 7: Distribution of the equivalent plastic strain in the deformed sample during HPT process in the case of an imposed compressive displacement of 1mm with a torsion angle of: (a) 15°, (b) 30°, (c) 45° and (d) 60°. Figure 8: Evolution of the equivalent plastic strain at the beginning and the end of torsion during HPT process at the upper area of the sample. Figure 9: Evolution of the average normal stress along the diameter of the upper face at the end of HPT process using a compressive displacement of 1mm with different torsion angles. Fig. 10 illustrates the contours plots of the mean normal stress distribution in the PMMA sample deformed during HPT process. It can be observed that the central part of the sample is stressed in compression, however the peripheral region (in yellow color) is under tensile stress and this zone increases with the increase of the torsion angle. This explains the heterogeneity of the equivalent plastic strain distribution. In addition, the maximum value is always located in the central part of the sample-anvil interfaces for the four torsion angles. Effect of temperature From the macroscopic point of view, the behavior of PMMA is highly dependent on temperature. It is the site of behavioural transitions that can be associated with different molecular relaxations. It is the site of behavioral transitions that can be associated with different molecular relaxations, i.e. activation of changes in local conformations, resulting in the most brutal behavior changes in certain temperature ranges. To study the dependence of PMMA behavior during HPT process under A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 191 the effect of temperature, a series of numerical calculations was carried out for three different temperatures T = {25, 60 and 80°C}. Noting that the same constitutive law described in the previous paragraph has been identified at different temperatures under isothermal conditions and the material parameters obtained are presented in Tab. 2. The same operating conditions previously used (compression of 1 mm with a torsion angle of 15°) were applied for the simulation of the HPT process under different temperatures. (a) (b) (c) (d) Figure 10: Distribution of the average normal stress in the deformed sample during HPT process in the case of a 1mm of compression with a torsion angle of: (a) 15°, (b) 30°, (c) 45° and (d) 60°. Distribution of the equivalent plastic strain In order to highlight the distribution of the plastic deformation in the PMMA sample, the iso-value contours of this latter at the end of the process for three different temperatures are presented in Fig. 11. It can be noted that the temperature has a slight influence on the magnitude of the plastic deformation during HPT process. A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 192 (a) (b) (c) Figure 11: Distribution of the equivalent plastic strain in the deformed sample with a compression of 1mm and a torsion angle of 15 ° for different temperatures: (a) 25°C, (b) 60°C and (c) 80°C. This slight influence of temperature can be confirmed by plotting the distribution of the equivalent plastic strain along the radial distance from the upper surface of the sample as shown in Fig. 12. A slight difference was found. In addition, the maximum value of equivalent plastic strain is approximately 2.36. Figure 12: Distribution of the equivalent plastic strain as a function of the radial distance in the upper surface of the sample during the HPT process for different temperatures: 25°C, 60°C and 80°C. Evolution of applied load and torque Fig. 13 shows the simulation results for the evolution of the compression load applied by the upper anvil as a function of the time during HPT process at various temperatures. The load required for each temperature in the ascending order (25°C, 60°C and 80°C) is respectively 11800, 16200 and 36000 N. We can see that when the temperature increases, the load required for HPT process decreases. For example, for a temperature of 25°C, the required load is more than three times of that required at room temperature. Moreover, the load versus time evolution during the HPT process can be subdivided into four main stages. The first stage starts from the origin up to the end of compression phase. This stage corresponds to the elastic deformation of compression. The second stage corresponds to the increase of the load from the elastic limit of A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 193 compression to the maximum value of the load. This stage corresponds to initiation and propagation of the plastic zone between the two upper and lower surfaces of the sample. The maximum force corresponds to the formation of the shear band along the cross-section of the sample, which precedes softening. The third stage corresponds to the decrease of the load from the peak to a minimum and is characterized by the initiation of the plastic flow. Beyond (stage 4), the load remains almost constant. This stage corresponds to the steady-state of the plastic flow when the plastic strain reaches its saturation value. Figure 13: Evolution of the compression load applied by the upper anvil as a function of the time of HPT process for different temperatures. Fig. 14 shows the actual values obtained by numerical simulation of the torque required for torsion phase of HPT process when it carried out at different temperatures (25°C, 60°C and 80°C). It can be noticed that there is an abrupt increase from the torsion phase followed by a plateau (a horizontal line) after about 1.8 sec. This aspect can be attributed to the plastic flow in the direction of rotation (when the elastic limit due to the torsion is exceeded). It should be noted that the conduct of the process with a temperature higher than room temperature and lower than the glass transition temperature of the PMMA allows to reduce the pressure and the torque required for the deformation of the sample during HPT process. Figure 14: Evolution of the torque as a function of time process for different temperatures. A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 194 Three particular properties can be highlighted when the sample is deformed at high temperatures: (a) reduced plastic strain distribution in the sample, (b) decrease of the required load and (c) decrease of the required torque. Furthermore, strain recovery during the cooling stage may affect the process. This point will be the subject of further mechanical and microstructural investigations. CONCLUSION he aim of this work is devoted to study the applicability of HPT process for a typical amorphous polymer which is polymethyl-methacrylate (PMMA). For this purpose, an identification of the parameters of a Perzyna elasto- viscoplastic phenomenological model was performed using compression tests at different temperatures and strain rates. After the identification of the parameters, a numerical investigation was carried out to highlight the effects of the main processing parameters in the case of an imposed compressive displacement and different torsion angles. From the results obtained, we can draw the following conclusions:  The equivalent plastic strain increases with the increase of the imposed compression displacement.  Accumulated plastic deformation increases with the increase of torsion angle and becomes very important at the edge of the sample.  It is advised to use small thicknesses to obtain a good homogeneity of the plastic deformation.  The temperature has a slight effect on the evolution of the plastic deformation. However, it can be used to reduce the force and torque required for HPT process if necessary. ACKNOWLEDGMENTS his research was supported by the General Directorate of Scientific Research and Technological Development (DGRSDT : Direction Générale de la Recherche Scientifique et du Développement Technologique) of Algeria. The authors gratefully acknowledge the scientific support of the two research teams from the Lille Mechanics Unit (France) and the Laboratory of Applied Biomechanics And Biomaterials (LABAB) of Oran (Algeria). REFERENCES [1] Alhajeri, S.N., Khaled, J., Al-Fadhalah, A., Almazrouee, I. and Langdon, T.G. (2016). Microstructure and microhardness of an Al-6061metal matrix composite processed by high-pressure torsion, Mater. Charact. , 118, pp. 270–278. DOI: 10.1016/j.matchar.2016.06.003. [2] Beloshenk, V., Vozniak, Iu., Beygelzimer, Y., Estrin, Y. and Kulagin, R. (2019). Severe Plastic Deformation of Polymers, Mater. Trans., 60(7), pp. 1192-1202. DOI: 10.2320/matertrans.MF201912 [3] Um, H.Y., Park, B.H., Ahn, D.H., AbdElAal, M.I., Park, J., Kim, H.S. (2017). Mechanical and biological behavior of ultrafine-grained Ti alloy aneurysm clip processed using high-pressure torsion, J. Mech. Behav. Biomed. Mater. , 68, pp. 203–209. DOI: 10.1016/j.jmbbm.2017.02.002. [4] Peitang, W., Cheng, L., Kiet T., Lihong, S., Guanyu, D., and Wenbin, H. (2017). A study on the texture evolution mechanism of nickel single crystal deformed by high pressure torsion, Mater. Sci. Eng. A, 684, pp. 239–248. DOI: 10.1016/j.msea.2016.11.098. [5] Yu, I., Kulagin, R., Fedorov, V., Mazilkin, A., Scherer, T., Baretzky, B. and Hahn, H. (2016). High Pressure Torsion Extrusion as a new severe plastic deformation process, Mater. Sci. Eng. A, 664, pp. 247–256. DOI: 10.1016/j.msea.2016.04.008. [6] Figueiredo, R.B., Pereira, P.H.R., Aguilar, M.T.P., Cetlin, P.R. and Langdon, T.G. (2012). Using finite element modeling to examine the temperature distribution in quasi-constrained high-pressure torsion, Acta Mater., 60, pp. 3190–3198. DOI: 10.1016/j.actamat.2012.02.027. [7] Valiev, R.Z., Song, C., McFadden, S.X., Mukherjee, A.K. and Mishra, R.S. (2001). TEM/HREM observations of nanostructured superplastic Ni3Al, Philos. Mag. A, 81, pp. 25–36. DOI: 10.1007/s10853-014-8204-5. T T A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 195 [8] Xu, C., Dobatkin, S.V., Horita, Z. and Langdon, T.G. (2009). Superplastic flow in a nanostructured aluminum alloy produced using high-pressure torsion, Mater. Sci. Eng. A, 500, pp. 170–175. DOI: 10.1016/j.msea.2008.09.049. [9] Harai, Y., Edalati, K., Horita, Z., and Langdon, T.G. (2009). Using ring samples to evaluate the processing characteristics in high-pressure torsion, Acta Mater., 57, pp. 1147–1153. DOI: 10.1016/j.actamat.2008.10.046. [10] Edalati, K., Horita, Z., Furuta, T. and Kuramoto, S. (2013). Dynamic recrystallization and recovery during high-pressure torsion: Experimental evidence by torque measurement using ring specimens, Mater. Sci. Eng. A, 559, pp. 506–509. DOI: 10.1016/j.msea.2012.08.132. [11] Horita, Z. and Langdon, T.G. (2008). Achieving exceptional superplasticity in a bulk aluminum alloy processed by high- pressure torsion, Scr.Mater., 58, pp. 1029–1032. DOI: 10.1016/j.scriptamat.2008.01.043. [12] Faraji, G., Kim, H.S., and Kashi, H.T. (2018). Severe Plastic Deformation Methods for Bulk Samples, Book chapter: Severe Plastic Deformation, pp. 37-112. DOI: 10.1016/B978-0-12-813518-1.00002-3. [13] Ivanisenko, Y., Kulagin, R., Fedorov, V., Mazilkin, A., Scherer, T., Baretzky, B. and Hahn, H. (2016). High Pressure Torsion Extrusion as a new severe plastic deformation process, Mater. Sci. Eng. A , 664, pp. 247-256. DOI: 10.1016/j.msea.2016.04.008. [14] Zhilyaev, A.P., and Langdon, T.G. (2008). Using high-pressure torsion for metal processing: Fundamentals and applications, Prog. Mater. Sci., 53, pp. 893–979. DOI: 10.1016/j.pmatsci.2008.03.002. [15] Lee, D.J., Yoon, E.Y., Ahn, D.H., Park, B.H., Park, H.W., Park, L.J., Estrin, Y. and Kim, H.S. (2014). Dislocation density-based finite element analysis of large strain deformation behavior of copper under high-pressure torsion, Acta Mater., 76, pp. 281-293. DOI: 10.1016/j.actamat.2014.05.027. [16] Vafai, R., Toroghinejad, M.R. and Pippan, R. (2012). Evaluation of mechanical behavior of nano-grained 2024 Al alloy during high pressure torsion (HPT) process at various temperatures, Mater. Sci. Eng. A, 536, pp. 73-81. DOI: 10.1016/j.msea.2011.12.064. [17] Peitang, W., Cheng, L., Kiet, T., Lihong, S.,Guanyu, D. and Wenbin, H. (2017). A study on the texture evolution mechanism of nickel single crystal deformed by high pressure torsion, Mater. Sci. Eng. A, 684, pp. 239-248. DOI: 10.1016/j.msea.2016.11.098. [18] Dimic, I., Cvijovic-Alagic, I., Volker B., Hohenwarter, A., Pippan, R., Veljovic, D., Rakin, M. and Bugarski, B. (2016). Microstructure and metallic ion release of pure titanium and Ti–13Nb–13Zr alloy processed by high pressure torsion, Mater Design., 91, pp. 340-347. DOI: 10.1016/j.matdes.2015.11.088. [19] Draï, A. and Aour, B. (2013). Analysis of plastic deformation behavior of HDPE during high pressure torsion process, Eng. Struct., 46, pp. 87-93. DOI: 10.1016/j.engstruct.2012.06.033. [20] Zhao, Y. and Zhang, J. (2007). Enhancement of yield strength in zirconium metal through high-pressure induced structural phase transition, Appl. Phys. Lett., 91, pp. 201907. DOI: 10.1063/1.2802726. [21] Feng, B., Levitas, V.I. and Kamrani, M. (2018). Coupled strain induced alpha to omega phase transformation and plastic flow flow in zirconium under high pressure torsion in a rotational diamond anvil cell, Mater. Sci. Eng. A, 731, pp. 623- 633. DOI: 10.1016/j.msea.2018.06.061. [22] Seenuvasaperumal, P., Doi, K., Basha, D.A., Singh, A., Elayaperumal, A. and Tsuchiya, K. (2018). Wear behavior of HPT processed UFG AZ31B magnesium alloy, Mater Lett., 227, pp. 194-198. DOI: 10.1016/j.matlet.2018.05.076. [23] Alhamidi, A., Edalati, K., Horita, Z., Hirosawa, S., Matsuda, K. and Terada, D. (2014). Softening by severe plastic deformation and hardening by annealing of aluminum–zinc alloy: Significance of elemental and spinodal decompositions, Mater. Sci. Eng. A, 610, pp. 17-27. DOI: 10.1016/j.msea.2014.05.026. [24] Edalati, K., Matsuo, M., Emami, H., Itano, S., Alhamidi, A., Staykov, A., Smith, D.J., Orimo, S., Akiba, E. and Horita, Z. (2016). Impact of severe plastic deformation on microstructure and hydrogen storage of titanium-iron- manganese intermetallics, Script.Mater., 124, pp. 108-111. DOI: 10.1016/j.scriptamat.2016.07.007. [25] Nikitina, M.A., R. Islamgaliev, K., and Sitdikov, V.D. (2016). Thermal stability of TiAl-based intermetallic alloys subjected to high pressure torsion, Mater. Sci. Eng. A, 651, pp. 306-310. DOI: 10.1016/j.msea.2015.10.121. [26] Matsunoshita, H., Edalati, K., Furui, M., and Horita, Z. (2015). Ultrafine-grained magnesium-lithium alloy processed by high-pressure torsion: Low-temperature superplasticity and potential for hydroforming, Mater. Sci. Eng. A, 640, pp. 443-448. DOI: 10.1016/j.msea.2015.05.103. [27] Edalati, K. and Horita, Z. (2016). A review on high-pressure torsion (HPT) from 1935 to 1988, Mater. Sci. Eng. A, 652, pp. 325-352. DOI: 10.1016/j.msea.2015.11.074. [28] Figueiredo, R.B., Aguilar, M.T.P., Cetlin, P.R. and Langdon, T.G. (2012). Analysis of plastic flow during high-pressure torsion, J. Mater. Sci., 47, pp. 7807–7814. DOI: 10.1007/s10853-012-6506-z. A. Drai et alii, Frattura ed Integrità Strutturale, 52 (2020) 181-196; DOI: 10.3221/IGF-ESIS.52.15 196 [29] Wang, W., Song, Y., Gao, D., Yoon, E.Y., Lee, D.J., Lee, C. and Kim, H.S. (2013). Analysis of stress states in compression stage of high pressure torsion using slab analysis method and finite element method, Met. Mater. Int., 19, pp. 1021-1027. DOI: 10.1007/s12540-013-5014-2. [30] Edalati, K., Miresmaeili, R., Horita, Z., Kanayam, H., and Pippan, R. (2011). Significance of temperature increase in processing by high-pressure torsion, Mater. Sci. Eng. A, 528, pp. 7301–7305. DOI: 10.1016/j.msea.2011.06.031. [31] Pereira, P.H.R., Figueiredo, R.B., Huang, Y., Cetlin, P.R., and Langdon, T.G. (2014). Modeling the temperature rise in high-pressure torsion, Mater. Sci. Eng. A, 593, pp. 185-188. DOI: 10.1016/j.msea.2013.11.015. [32] Kamrani, M., Levitas, V.I., and Feng, B. (2017). FEM simulation of large deformation of copper in the quasi-constrain high pressure-torsion setup, Mater. Sci. Eng. A, 705, pp. 219–230. DOI: 10.1016/j.msea.2017.08.078. [33] Kulagin, R., Zhao, Y., Beygelzimer, Y., Toth, L.S. and Shtern, M. (2017). Modeling strain and density distributions during high-pressure torsion of pre-compacted powder materials, Mater. Res. Lett. 5, pp. 179–186. DOI: 10.1080/21663831.2016.1241318. [34] Lee, D.J. and Kim, H.S. (2014). Finite element analysis for the geometry effect on strain inhomogeneity during high pressure torsion, J. Mater. Sci., 49, pp. 6620-6628. DOI: 10.1007/s10853-014-8283-3. [35] Belayachi, N., Benseddiq, N. and Nait-Abdelaziz, M. (2008). Behavior of the heterogeneous glassy polymers: Computational modeling and experimental approach, Compos. Sci. Technol. 68, pp. 367–375. DOI: 10.1016/j.compscitech.2007.07.002 [36] Perzyna, P. (1985). On Constitutive Modelling of Dissipative Solids for Plastic Flow, Instability and Fracture. In: A. Sawczuk (Eds), (1985), Plasticity Today-Modelling, Methods and Applications, pp. 657-679. [37] Perzyna, P. (1966). Fundamental Problems in Viscoplasticity, In: Advances in Applied Mechanics, Academic press in New York., 9, pp. 243-377. DOI: 10.1016/S0065-2156(08)70009-7. [38] Perzyna, P., and Pecherski, R. (1983). Modified theory of viscoplasticity. Physical foundations and identification of material functions for advanced strains, Arch. Mech., 35, pp. 423-436. [39] Peric, D., Owen, D.R.J. and Honnor, M.E. (1992). A model for finite strain elasto-plasticity based on logarithmic strains: Computational Issues. Computer methods in applied mechanics and engineering, 94, pp. 35-61. DOI: 10.1016/0045-7825(92)90156-E. [40] Sluis, O.V.D., Schreurs, P.J.G, Brekelmans, W.A.M. and Meijer, H.E.H. (2000). Overall behavior of heterogeneous elastoviscoplastic materials: effect of microstructural modeling. Mech. Mat., 32, pp. 449-462. DOI: 10.1016/S0167-6636(00)00019-3. [41] Sluis, O.V.D., Schreurs, P.J., and Meijer, H.E.H. (2001). Homogenisation of structured elastoviscoplastic solids at finite strains, Mech. Mat., 33, pp. 499-522. DOI: 10.1016/S0167-6636(01)00066-7. [42] Draï, A. and Aour, B. (2017). Analysis of the temperature effect on the plastic strain of polymers during high pressure torsion (HPT) process, ICONTES2017: International Conference on Technology, Engineering and Science, October 26 - 29, 2017 Antalya/Turkey, The Eurasia Proceedings of Science, Technology, Engineering & Mathematics (EPSTEM), 1, pp. 40-51. << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice