Prerpints.org logo
Preprint
Article

Linear Drucker-Prager Yield Criterion Calibration for Polypropylene under Multiaxial Load

This version is not peer-reviewed.

Submitted:

12 August 2024

Posted:

13 August 2024

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
This study focuses on calibrating the linear Drucker-Prager yield criterion for polypropylene under multiaxial loading conditions, using the Arcana machine to achieve a wide range of triaxiality states. The Drucker-Prager criterion has been adapted for polymers, composites, and metals due to its ability to account for the influence of hydrostatic pressure on yield stress. This study employs the Arcana machine, which allows for testing flat plate samples under various angles, thereby simulating different stress states from pure shear to biaxial tension. Additionally, uniaxial compression tests were conducted to extend the range of triaxiality. The material used in this study is Sabic 83MF10, a polypropylene polymer. Samples were prepared by injection molding and cut to specific dimensions. Tests were performed using a universal testing machine with the Arcana handle, and the results were analyzed to determine the yield strength, pressure, and triaxiality for each sample. The results showed a linear relationship between von Mises stress and hydrostatic pressure, with a friction angle (β) of 20.65° and material cohesion (d) of 27.81. The numerical simulations in Abaqus confirmed the validity of the Drucker-Prager model, accurately reflecting the moment of yielding for the tested samples.
Keywords: 
Subject: 
Engineering  -   Mechanical Engineering

1. Introduction

The introduction of polymeric materials, specifically polypropylene, to the market in 1951 [1], represented a pivotal advancement in the field of materials science and engineering. These polymers have since become integral across a diverse array of industrial applications like automotive, medical, construction, electronics and consumer goods [2,3,4]. The preference for polymeric materials, particularly polypropylene, is based on their outstanding mechanical properties, intrinsic resistance to aging, capability to maintain dimensional stability throughout their lifecycle, superior resistance to crack initiation or propagation, and their suitability for high-volume manufacturing processes such as extrusion and injection molding [5,6,7].
Effective design process of designing and validating component made of polymers in industries requires the usage of simulation methods, including the Finite Element Method [8]. To make this possible, it is necessary to accurately represent the behavior of plastics under the different loading regimes and to relative them in the numerical material model [9,10,11]. In addition to the commonly known nonlinear behavior of thermoplastics, influence of strain rate, temperature, stress state and loading type [12,13,14,15], the critical phenomena is the influence of hydrostatic pressure on the yield stress, plastic flow, the destruction criterion [16,17,18]. Studies conducted in the past have clearly demonstrated the occurrence of different elastoplastic properties for loads of different characteristics resulting from tests ranging from compression, through shear, to tension [19,20,21,22].
One of the widely used constitutive model capable to reproduce the plastic behavior is the linear Drucker-Prager criterion, which makes the transition of the material to the plastic range, i.e. the value of the yield stress, dependent on the value of hydrostatic pressure. Thanks to this, it is possible to distinguish the behavior of polymers under different load conditions, including fact that for most of thermoplastics the yield strength for compression is higher than for tension [23,24,25,26]. The method originally comes from soil research but it has found wide application for plastics, composites and metals [27,28,29,30].
This article focuses on the use of the Arcana machine which is widely used when the only possible form of material samples are flat plates of a defined thickness [31,32,33,34,35,36]. It is necessary to properly cut the plates to obtain a multiaxial state of tension using an appropriate fillet radius at the necking point. An additional verification is a uniaxial compression test with extend criteria to compression state. The obtained material model is finally tested in the Abaqus numerical software to confirm the compliance of the obtained parameters for the Drucker-Prager model and the numerical simulations carried out using it with the actual behavior of polypropylene samples.

2. Material Behavior

Stress tensor ( σ i j ) according to the matrix form from Equation 1 could be explain using hydrostatic (mean) ( σ h ) stress tensor responsible for volumetric change and deviatoric stress tensor ( s i j ) cause body shape deviation.
σ i j = σ x x τ x y   τ x z   τ y x   σ y y τ y z   τ z x   τ z y   σ z z = σ h 0   0   0   σ h 0   0   0   σ h + s x x s x y   s x z   s y x   s y y s y z   s z x   s z y   s z z = σ h δ i j + s i j
where Kronecker Delta function is defined by:
δ i j = 1 ,   i f   i = j 0 ,   i f   i ≠ j
Dilatation state defined by pressure ( p ) could be drawn using stress tensor first invariant ( I 1 ):
p = − σ h = − 1 3 t r a c e σ i j = − σ x x + σ y y + σ z z 3 = − I 1 3
Equivalent stress called von Mises ( q ) is defined by second invariant of deviator state ( J 2 ):
q = σ e q = 3 J 2
J 2 = σ x x − σ y y 2 + σ y y − σ z z 2 + σ z z − σ x x 2 6 + τ x y + τ y z + τ z x
Thirst invariant of deviator state ( J 3 ) needed to explain Drucker-Prager criterion could by defined by:
r = J 3 = d e t s i j
d e t s i j = σ x x σ y y σ z z + τ x y τ y z τ z x + τ x z τ y x τ z y − τ x z σ y y τ z x − τ x y τ y x σ z z − σ x x τ y z τ z y
The linear Drucker–Prager criterion commonly used for polymers can be represented in the meridian plane by the yield surface using deviatoric stress factor ( t ) as a function of hydrostatic pressure ( p ) as in Figure 1, where the parameters defining its location are the friction angle ( β ) and material cohesion ( d ) according to Equation 8:
t = d + t a n ⁡ β
where deviatoric stress factor ( t ) could be defined by Equation 9:
t = 1 2 q 1 + 1 K − 1 − 1 K r q 3
using ratio of yield stress for triaxial tension and compression ( K ), what could be described on deviatoric plane in Figure 2. Based on the research on Polypropylene [22] and due to the limited scope of possible tests on the Arcana machine, a K factor of 1 was assumed. Due to this t = q on Figure 1.
Plastic flow direction ( d ε p l ) for polymer could be defined in Equation 10 by dilatation angle ( ψ ) using Poison ration ( ν ) [37,38].
ψ = t a n − 1 3 1 − 2 ν 2 1 + ν
In Equation 8:
  • tangent of friction angle ( β ) represent the sensitivity for hydrostatic pressure ( p ),
  • material cohesion ( d ) reflect isotropic hardening, which could be defined by Equation 11 using uniaxial tension ( σ t ), compression ( σ c ) or shear ( τ ) yield stress,
    d = 1 K + 1 3 t a n β σ t = 1 − 1 3 t a n β σ c = 3 2 1 + 1 K τ
Drucker-Prager yield surface in the reverse principal axes is a cone around hydrostatic axis with a circular cross section for K factor equal to 1 as on Figure 3.
To clearly determine the loading state of the material, the triaxiality ( η ) parameter (Equation 12) should be used, the basic values of which are given in Table 1.
η = − p q

3. Material, Samples and Experimental Techniques

The following tests were performed for thermoplastic Sabic 83MF10 material, universal polypropylene polymer provided by Saudi Basic Industries Corporation (SABIC) and the basic parameters of which are presented in Table 2. Dilatation angle was calculated based on Equation 2.
Due to the fact that despite the lack of reinforcing inclusions in the material such as fibers, noticeable anisotropy of properties is possible [41,42], it was decided to cut samples from the plate far away from the injection points. This allowed to obtain a quasi-random arrangement of molecules.
Tension and shear samples subjected to analysis on a Zwick type universal testing machine with Arcana handle was obtained in the form of flat plates received in injection molding process, with a thickness of 2,2 m m . The samples were machine cut with the radius of 1.5 m m and dimensions described on Figure 4.
The sample for uniaxial compression tests is a single cylinder with a diameter of 3 m m cut from the previously mentioned plates.
To achieve wide range of triaxiality η and according to it different load type in stated between tension and shear, the Arcana machine was used (Figure 5). This type of holder allows to stretch samples at angles from 0 to 90 degrees. It consists of a pair of holders that secure it to the testing machine, two semicircular discs with holes that allow for mounting at 7.5 degrees and holders that grip the tested sample directly. The entire structure is made of aluminum alloy 7075 (Fortal), ensuring high rigidity of the structure.
All test on Arcana machine were performed with constant head travel speed of 1 m m / s which allows to achieve strain rate ( ε ˙ ) between 0.06 and 0.11 1 / s during the tests. Due to the selected machine speed it was possible to avoid, on the one hand, the influence of material hardening at high strain rates described e.g. by Cowper-Symonds model [39,40] and, on the other hand, the effect of polymer relaxation [43,44].

4. Experimental Data Processing and Results

The linear Drucker-Prager criterion was chosen for further calibration due to the impossibility of obtaining results for triaxial tension on flat samples with undercut. Figure 6 shows the results of tension (and shear) tests on the Arcana machine and the result of the compression test of a cylindrical sample. For each variant, three tests were performed, and two extreme cases (max and min values) in terms of yield strength were rejected from further evaluation.
The following tests were performed to determine the yield strength, which according to ISO 527-1 defines yielding as the first deviation in true stress-strain curve from pure elastic response due to the lack of a distinct yield strength. Force ( F ) and displacement ( d i s p ) were measure using Zwick embedded cross-head sensor. The values of triaxiality ( η ), pressure ( p ) and equivalent von Mises stress ( q ) were determined from FEM simulations representing the behavior of the samples based on the calibration of the material curve with the load case as a average value for whole cross section in the critical most strained surface. All results needed for further evaluation were described below in Table 3.
The larger the angle on the Arcana machine, the lower the yield strength seems to be. This does not reflect the actual parameters, due to the complex stress state, different for each sample. The decreasing apparent stiffness of the sample with increasing angle results from the transition from pure tension to pure shear, where instead of Young's modulus ( E ), the basic parameter defining elasticity is the Kirchhoff’s ( G ) modulus, which calculated according to Equation 13 is 2.8 times lower.
G = E 2 1 − ν
All the result in x relationship were described as a linear function of Drucker-Prager yield criteria in Figure 7. The regression line was obtained using solver which has the objective function with scope to minimize the R 2 error value.
As can be seen in the above graph, the results obtained for the Arcana machine correspond to the values from pure shear ( η = 0 ) to biaxial tension ( η = 2 / 3 ), which was possible thanks to the proposed undercuts on flat samples. The disc compressed on the strength machine obtained a result slightly lower ( η < − 1 / 3 ) than for uniaxial compression due to its relatively small thickness and the friction occurring in the test between the upper and lower external surfaces of the sample and the contact linings of the testing machine. Final values for linear Drucker-Prager criteria were described in Table 4.
All points achieved an average R 2 error value of 0.35 , which suggests a very good fit to the linear relationship. An extension of the proposed linear criterion can be the tension cut-off, which determines the limiting value of the tensile pressure ( p t ), beyond which plastic deformation occurs, regardless of the value of von Mises stress ( q ) [45].
In order to confirm the correct validation of the Drucker-Prager model, simulations were performed on previously tested samples and the calculation results in Abaqus were compared to the values obtained on the strength test bench. Implicit simulations were used using Abaqus Standard solver. The discretized models used in the simulations (shown in Figure 8) for Arcana and compression have respectively 16263 C3D10, 10-node quadratic tetrahedron elements with characteristic size of 0.5 m m in critical space and 2200 C3D20, 20-node quadratic hexahedral elements with 0.2 m m for whole geometry.
The finite element simulation was performed in a way that was designed to reproduce the actual test conditions as precisely as possible. For the samples tested on the Arcana machine, the model was completely fixed on the lower side in the locations indicated as mounting place in Figure 4. On the other side, a displacement was imposed in the vertical direction. The Arcana machine itself was not included in the simulation due to the stiffness of the entire structure made of aluminum that was several orders of magnitude higher. For the disc compressed on both sides, contact surfaces were simulated with a friction coefficient of 0.2 − , of which only one, the upper surface, moves towards the lower fixed plane. In order to follow the simulation progress above the yield point, the yield curve shown in Figure 9 was determined using numerical calibration from sample no 1 (Arcana – 0deg) test curve.
The hardening behavior was determined for range from 0 to 0.4 of plastic strain and implemented in Abaqus numerical software in tabular form with 17 points. Based on the above data, a comparison (on Figure 10) between previous laboratory results and the numerical simulation results was made for all proposed samples.
As shown in the above graph, the simulation result shows the moment of yielding with a satisfactory result. The moment of appearance of the first plasticized finite elements and the transition of the entire critical cross-section into the plastic range are clearly visible. Due to the determination of the plasticity curve used for simulation based on sample no. 1, the simulation course outside the elastic range is best shown. For models with an angle of setting above 30°, the laboratory results are characterized by a less clear moment of yielding compared to the simulation results. This may result from the settings in the Abaqus program, which do not take into account the possibility of small angular movements of the samples at the place of attachment, which can occur during tests on the Arcana machine. For the compressed sample, the moment of yielding is shown satisfactorily in the numerical model, however, a comparison of the course in the plastic range shows clearly that the yield curve differs significantly from that obtained from the tensile sample. The anisotropic nature of the elastic range is noticeable, which is demonstrated by the fact that for Arcana machine settings close to 90° introducing a shear state in the sample, the Young Modulus ( E ) values ​​are nearly 10% lower than those originally used in the simulation.

5. Conclusions and Areas for Further Research

The comprehensive study of Drucker-Prager presented in this article has shown a noticeable and significant linear dependence of the yield strength on the material state and, consequently, an unequivocal effect of pressure on the criterion under study. The proposed approach using a specially manufactured Arcana machine with samples with special undercuts allows for the representation of the state from pure shear to biaxial tension, which is not possible using standard samples. The supplementation and extension of the studies to the scope of uniaxial compression is a compression test performed on a slender disc. The obtained results uniformly achieve a linear character of variation.
All previously laboratory tested samples were numerically simulated, which allowed us to confirm the ability of the Drucker-Prager criterion implemented in the Abaqus program to properly reflect the moment of plasticization of samples. Due to the demonstrated variable course of the plasticization curve for different material loading states (Triaxiality), it seems advisable to take this phenomenon into account in future studies. At the moment, there are no ready-made solutions. Due to the stiffening of the numerical model of the sample mounting points in the holder, a small deviation occurred for the elastic range in samples with a machine inclination angle >15°.
Further research should focus on the longer range of Triaxiality, above 2-axis stretching and determining the so-called tension cut-off, which defines the pressure at which the material immediately plasticizes. A certain simplification of this phenomenon may be an attempt to calibrate the exponent or hyperbolic model, which indirectly takes into account the mentioned phenomenon.

Acknowledgments

The article was created as part of the 6th edition of the program of the Minister of Education and Science entitled “Implementation doctorate” Agreement No. DWD/6/0399/2022.

References

  1. Stinson, Stephen (1987). "Discoverers of Polypropylene Share Prize". Chemical & Engineering News. 65 (10): 30. [CrossRef]
  2. Fried, J. R. (2014). Polymer Science and Technology (Illustrated ed.). Prentice Hall. ISBN 0137039557, 9780137039555.
  3. Young, R. J., & Lovell, P. A. (2011). Introduction to Polymers (3rd ed.). CRC Press. [CrossRef]
  4. Vasile, C., & Pascu, M. (2005). Practical Guide to Polyethylene. Smithers Rapra Technology. ISBN 978-1859574935.
  5. Maddah, H. A. (2016). Polypropylene as a Promising Plastic: A Review. American Journal of Polymer Science, 6(1), 1-11. [CrossRef]
  6. Hossain, M. T., Shahid, M. A., Mahmud, N., Habib, A., Rana, M. M., Khan, S. A., & Hossain, M. D. (2024). Research and application of polypropylene: A review. Nano, 19(2). [CrossRef]
  7. Karger-Kocsis, J., & Bárány, T. (2019). Polypropylene Handbook: Morphology, Blends and Composites. SpringerLink. [CrossRef]
  8. Maier, C., & Calafut, T. (1998). Applications. In C. Maier & T. Calafut (Eds.), Polypropylene (pp. 87-107). William Andrew Publishing. ISBN: 9781884207587.
  9. Johri, N., Agarwal, G., Mishra, R. K., & Thakur, H. C. (2022). FEM analysis of polymeric hybrid composites. Materials Today: Proceedings, 57(Part 2), 383-390. [CrossRef]
  10. Du, S., Hamdi, M., & Sue, H.-J. (2020). Experimental and FEM analysis of mar behavior on amorphous polymers. Wear, 444-445, 203155. [CrossRef]
  11. Gariya, N., Prasad, B., & Kumar, P. (2020). FEM based analysis of soft polymer composites with crack. Materials Today: Proceedings, 28(Part 4), 2426-2430. [CrossRef]
  12. Ayoub, G., Zaïri, F., Naït-Abdelaziz, M., & Gloaguen, J. M. (2010). Modelling large deformation behaviour under loading–unloading of semicrystalline polymers: Application to a high density polyethylene. International Journal of Plasticity, 26(3), 329-347. [CrossRef]
  13. Laiarinandrasana, L., Besson, J., Lafarge, M., & Hochstetter, G. (2009). Temperature dependent mechanical behaviour of PVDF: Experiments and numerical modelling. International Journal of Plasticity, 25(7), 1301-1324. [CrossRef]
  14. Arruda, E. M., Boyce, M. C., & Jayachandran, R. (1995). Effects of strain rate, temperature and thermomechanical coupling on the finite strain deformation of glassy polymers. Mechanics of Materials, 19(2-3), 193-212. [CrossRef]
  15. Chou, S. C., Robertson, K. D., & Rainey, J. H. (1973). The effect of strain rate and heat developed during deformation on the stress-strain curve of plastics - Temperature rise developed during deformation can have significant effects on the stress-strain relationship. Four hard plastics are tested at various strain rates, and temperature changes are measured during deformation of the specimen. Experimental Mechanics, 13(10), 422-432. [CrossRef]
  16. Alves, M., & Jones, N. (1999). Influence of hydrostatic stress on failure of axisymmetric notched specimens. Journal of the Mechanics and Physics of Solids, 47(3), 643-667. [CrossRef]
  17. Abdul-Hameed, H., Messager, T., Ayoub, G., Zaïri, F., Naït-Abdelaziz, M., Qu, Z., & Zaïri, F. (2014). A two-phase hyperelastic-viscoplastic constitutive model for semi-crystalline polymers: Application to polyethylene materials with a variable range of crystal fractions. Journal of the Mechanical Behavior of Biomedical Materials, 37, 323-332. [CrossRef]
  18. Hong, K., & Strobl, G. (2008). Characterizing and modeling the tensile deformation of polyethylene: The temperature and crystallinity dependences. Polymer Science - Series A, 50(5), 483-493. [CrossRef]
  19. Hachour, K., Zaïri, F., Naït-Abdelaziz, M., Gloaguen, J. M., Aberkane, M., & Lefebvre, J. M. (2014). Experiments and modeling of high-crystalline polyethylene yielding under different stress states. International Journal of Plasticity, 54, 1-18. [CrossRef]
  20. Rinaldi, R., Gaertner, R., Brunet, M., Chazeau, L., Vidal-Sallé, E., & Gauthier, C. (2011). Modeling of the mechanical behavior of amorphous glassy polymer based on the quasi-point defect theory—Part II: 3D formulation and finite element modeling of polycarbonate. International Journal of Non-Linear Mechanics, 46(3), 507-518. [CrossRef]
  21. Olsson, R. (2011). A survey of test methods for multiaxial and out-of-plane strength of composite laminates. Composites Science and Technology, 71(6), 773-783. [CrossRef]
  22. Manaia, J. P., Pires, F. A., & de Jesus, A. M. P. (2019). Elastoplastic and fracture behaviour of semi-crystalline polymers under multiaxial stress states. Frattura ed Integrità Strutturale, 13(47), 82-103. [CrossRef]
  23. Pae, K. D., & Bhateja, S. K. (1975). The Effects of Hydrostatic Pressure on the Mechanical Behavior of Polymers. Journal of Macromolecular Science, Part C, 13(1), 1–75. [CrossRef]
  24. Olufsen, S., Clausen, A. H., & Hopperstad, O. S. (2019). Influence of stress triaxiality and strain rate on stress-strain behaviour and dilation of mineral-filled PVC. Polymer Testing, 75, 350-357. [CrossRef]
  25. Pae, K. D., & Bhateja, S. K. (1975). The Effects of Hydrostatic Pressure on the Mechanical Behavior of Polymers. Journal of Macromolecular Science, Part C, 13(1), 1–75. [CrossRef]
  26. Ward, I. M. (1971). Review: The yield behaviour of polymers. Journal of Materials Science, 6(11), 1397-1417. [CrossRef]
  27. Drucker, D. C., & Prager, W. (1952). Soil mechanics and plastic analysis for limit design. Quarterly of Applied Mathematics, 10(2), 157–165. [CrossRef]
  28. Pulungan, D., Lubineau, G., Yudhanto, A., Yaldiz, R., & Schijve, W. (2017). Identifying design parameters controlling damage behaviors of continuous fiber-reinforced thermoplastic composites using micromechanics as a virtual testing tool. International Journal of Solids and Structures, 117, 177-190. [CrossRef]
  29. Naya, F., González, C., Lopes, C. S., Van der Veen, S., & Pons, F. (2017). Computational micromechanics of the transverse and shear behavior of unidirectional fiber reinforced polymers including environmental effects. Composites Part A: Applied Science and Manufacturing, 92, 146-157. [CrossRef]
  30. Chen, J., Wan, L., Ismail, Y., Hou, P., Ye, J., & Yang, D. (2021). Micromechanical analysis of UD CFRP composite lamina under multiaxial loading with different loading paths. Composite Structures, 269, 114024.
  31. Popelar, C., Liechti, K. A Distortion-Modified Free Volume Theory for Nonlinear Viscoelastic Behavior. Mech Time-Depend Mater 7, 89–141 (2003). [CrossRef]
  32. Hao, P., Ud Din, I., & Panier, S. (2019). Development of Modified Arcan Fixture for biaxial loading response of fiber-reinforced composites. Polymer Testing, 80, 106148. [CrossRef]
  33. Cognard, J. Y., Sohier, L., & Davies, P. (2011). A modified Arcan test to analyze the behavior of composites and their assemblies under out-of-plane loadings. Composites Part A: Applied Science and Manufacturing, 42(1), 111-121. [CrossRef]
  34. Hung, S.C., Liechti, K.M. An evaluation of the arcan specimen for determining the shear moduli of fiber-reinforced composites. Experimental Mechanics 37, 460–468 (1997). [CrossRef]
  35. Arcan, M., Hashin, Z., & Voloshin, A. (1978). A method to produce uniform plane-stress states with applications to fiber-reinforced materials - A specially designed specimen yields material properties under pure shear or uniform plane-stress conditions. Experimental Mechanics, 18(4), 141-146. [CrossRef]
  36. Voloshin, A., & Arcan, M. (1980). Pure shear moduli of unidirectional fibre-reinforced materials (FRM). Fibre Science and Technology, 13(2), 125-134. [CrossRef]
  37. Morelle, X.P., Chevalier, J., Bailly, C. et al. Mechanical characterization and modeling of the deformation and failure of the highly crosslinked RTM6 epoxy resin. Mech Time-Depend Mater 21, 419–454 (2017). [CrossRef]
  38. Quino, G., Gargiuli, J., Pimenta, S., Hamerton, I., Robinson, P., & Trask, R. S. (2023). Experimental characterisation of the dilation angle of polymers. Polymer Testing, 125, 108137. [CrossRef]
  39. Xu, M.-m.; Huang, G.-y.; Feng, S.-s.; McShane, G.J.; Stronge, W.J. Static and Dynamic Properties of Semi-Crystalline Polyethylene. Polymers 2016, 8, 77. [CrossRef]
  40. Marangoni, A. L., & Massaroppi, E.. (2017). Cowper-Symonds parameters estimation for ABS material using design of experiments with finite element simulation. Polímeros, 27(3), 220–224. [CrossRef]
  41. Schrauwen, B. A. G., van Breemen, L. C. A., Spoelstra, A. B., Govaert, L. E., Peters, G. W. M., & Meijer, H. E. H. (2004). Structure, deformation, and failure of flow-oriented semicrystalline polymers. Macromolecules, 37(23), 8618-8633. [CrossRef]
  42. Katti, S. S., & Schultz, M. (1982). The microstructure of injection-molded semicrystalline polymers: A review. Polymer Engineering & Science, 22(16), 1001-1017. [CrossRef]
  43. Swallowe, G.M. (1999). Relaxations in Polymers. In: Swallowe, G.M. (eds) Mechanical Properties and Testing of Polymers. Polymer Science and Technology Series, vol 3. Springer, Dordrecht. [CrossRef]
  44. Jourdan, C., Cavaille, J. Y., & Perez, J. (1989). Mechanical relaxations in polypropylene: A new experimental and theoretical approach. Journal of Polymer Science Part B: Polymer Physics, 27(11), 2171-2184. [CrossRef]
  45. Krabbenhoft, K., Karim, M. R., Lyamin, A. V., & Sloan, S. W. (2012). Associated computational plasticity schemes for nonassociated frictional materials. International Journal for Numerical Methods in Engineering, 90(9), 1089–1117. [CrossRef]
Preprints 115028 g001
Figure 1. Yield surface and flow direction for linear Drucker-Prager model.
Figure 1. Yield surface and flow direction for linear Drucker-Prager model.
Preprints 115028 g001
Preprints 115028 g002
Figure 2. Yield surfaces for different ratio of yield stress for triaxial tension and compression ( K ) on deviatoric plane
Figure 2. Yield surfaces for different ratio of yield stress for triaxial tension and compression ( K ) on deviatoric plane
Preprints 115028 g002
Preprints 115028 g003
Figure 3. Yield surfaces for in the reverse principal axes for linear Drucker-Prager criterion.
Figure 3. Yield surfaces for in the reverse principal axes for linear Drucker-Prager criterion.
Preprints 115028 g003
Preprints 115028 g004
Figure 4. Sample scheme for tension and shear tests on arcana machine.
Figure 4. Sample scheme for tension and shear tests on arcana machine.
Preprints 115028 g004
Preprints 115028 g005
Figure 5. Arcana machine for universal testing machine in different angle position.
Figure 5. Arcana machine for universal testing machine in different angle position.
Preprints 115028 g005
Preprints 115028 g006
Figure 6. Force vs Head travel for Zwick from Arcana or compression tests with the yield points.
Figure 6. Force vs Head travel for Zwick from Arcana or compression tests with the yield points.
Preprints 115028 g006
Preprints 115028 g007
Figure 7. Linear Drucker-Prager yield von Mises stress as a function of pressure.
Figure 7. Linear Drucker-Prager yield von Mises stress as a function of pressure.
Preprints 115028 g007
Preprints 115028 g008
Figure 8. Numerical models overview for Arcana tension and compression tests.
Figure 8. Numerical models overview for Arcana tension and compression tests.
Preprints 115028 g008
Preprints 115028 g009
Figure 9. Hardening behavior (plasticity curve) for numerical calculation calibrated from sample no 1 (Arcana – 0deg).
Figure 9. Hardening behavior (plasticity curve) for numerical calculation calibrated from sample no 1 (Arcana – 0deg).
Preprints 115028 g009
Preprints 115028 g010
Figure 10. Force vs Head travel comparison between laboratory Zwick from Arcana or compression tests and numerical simulation results.
Figure 10. Force vs Head travel comparison between laboratory Zwick from Arcana or compression tests and numerical simulation results.
Preprints 115028 g010
Table 1. Example of triaxiality values for different stress states.
Table 1. Example of triaxiality values for different stress states.
Stress state Triaxiality   ( η )
Triaxial tension ∞
Biaxial tension 2 / 3
Uniaxial tension 1 / 3
Shear 0
Uniaxial compression − 1 / 3
Biaxial compression − 2 / 3
Triaxial compression − ∞
Table 2. Material data set for Sabic 83MF10.
Table 2. Material data set for Sabic 83MF10.
Parameter Value
Density 905 k g / m 3
Young modulus 1200 M P a
Poisson ration 0.4 −
Dilatation angle 12.1 °
Table 3. Material data set for Sabic 83MF10.
Table 3. Material data set for Sabic 83MF10.
Sample Yield   strength   ( q )
M P a 		Pressure   ( p )
M P a 		Triaxiality   ( η )
−
No 1 – Arcana – 0 deg 22.24 − 13.57 0.61
No 2 – Arcana – 15 deg 23.12 − 12.03 0.52
No 3 – Arcana – 30 deg 23.29 − 10.24 0.44
No 4 – Arcana – 45 deg 24.64 − 8.01 0.33
No 5 – Arcana – 60 deg 26.72 − 5.72 0.21
No 6 – Arcana – 75 deg 26.96 − 3.54 0.13
No 7 – Arcana – 90 deg 28.21 − 0.39 0.01
No 8 – compression 31.83 12.41 − 0.39
Table 4. Drucker-Prager linear model for Polypropylene Sabic 83MF10.
Table 4. Drucker-Prager linear model for Polypropylene Sabic 83MF10.
Parameter Value
β 20.65 °
d 27.81
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Alerts
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2025 MDPI (Basel, Switzerland) unless otherwise stated