Next Article in Journal
Sulfonated Polyether Ketone Membranes Embedded with Nalidixic Acid—An Emerging Controlled Drug Releaser
Next Article in Special Issue
Petal-like Patterning of Polylactide/Poly (Butylene Succinate) Thin Films Induced by Phase Separation
Previous Article in Journal
Antimicrobial Wound Dressings based on Bacterial Cellulose and Independently Loaded with Nutmeg and Fir Needle Essential Oils
Previous Article in Special Issue
Experimental Study of Curing Temperature Effect on Mechanical Performance of Carbon Fiber Composites with Application to Filament Winding Pressure Vessel Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Study Influence Analysis of the Mathematical Model Choice for Describing Polymer Behavior

by
Anna A. Kamenskikh
*,
Yuriy O. Nosov
and
Anastasia P. Bogdanova
Department of Computational Mathematics, Mechanics and Biomechanics, Perm National Research Polytechnic University, 614990 Perm, Russia
*
Author to whom correspondence should be addressed.
Polymers 2023, 15(17), 3630; https://doi.org/10.3390/polym15173630
Submission received: 31 July 2023 / Revised: 25 August 2023 / Accepted: 30 August 2023 / Published: 1 September 2023
(This article belongs to the Special Issue Polymer Composites in Engineering: Multiscale/Multiphysics Analyses)

Abstract

:
The article considered the three types of description of the material behavior model: elastic, elastic–plastic, and viscoelastic. The problem is considered in the framework of deformable solid mechanics. The paper considers the possibility of describing modern polymeric and composite materials used as antifriction sliding layers in the viscoelasticity framework. A numerical procedure for finding the coefficients to describe the viscoelastic material behavior using the Prony model has been implemented. Numerical results and experimental data are compared. The model problem of spherical indenter penetration into polymer half-space is realized. The influence of the system discretization on the numerical solution is analyzed. The influence of the polymer behavior description in static and dynamic problem formulations is analyzed.

Graphical Abstract

1. Introduction

1.1. Research Objectives

The investigate different types of mathematical models of behavior describing of a polymer material is the purpose of the work: elasticity; elastoplasticity; viscoelasticity. As well as to evaluate the possibility of applying the model depending on the type of analysis and load.
Research objectives:
  • The behavior description of polymeric materials in viscoelastic terms using the generalized Maxwell model;
  • Comparison of experimental data and numerical solutions of three different mathematical models of material behavior: elastic, plastic, and viscoelastic;
  • Carrying out verification of the numerical Hertz problem;
  • Constructing numerical models of the specimen behavior (elastic–plastic, viscoelastic);
  • The influence study of mathematical models of polymers’ behavior in a dynamic setting (loading and subsequent exposure for 1 h at room temperature).

1.2. Problem Context

The study of various characteristics of materials, as well as their application possibility in this or that field, is an urgent task among scientists all over the world [1,2,3]. One of the modern popular materials is polymers. Such materials are useful because, with the help of certain combinations of their structural links, it is possible to obtain materials with different characteristics [4]. In all production areas, there is a search for a type of polymeric material that will have all possible positive characteristics but, at the same time, will have low cost and a low environmental footprint, both in manufacturing and disposal [5,6]. Scientists from all over the world create certain combinations of monomers and further investigate their characteristics and properties [7]. One important quality of polymers is the combination of various useful properties, such as impact resistance, electrical conductivity, and insulation [8], plasticity and elasticity, lightness, chemical resistance, impermeability, corrosion resistance, etc.
In particular, one of the actively used polymeric materials in industry is PTFE, along with its various modifications [9,10]. It is a material that has unique properties such as high chemical resistance, low thermal conductivity, abrasion, and corrosion resistance, and high dielectric constant. Both pure PTFE and composites based on it are used in various industries, such as mechanical engineering, the oil and gas industry, electrical engineering, medicine, etc. PTFE and its modifications are widely used in the production of electronic components and equipment such as capacitors, resistors, switches, and other components since they have a high dielectric constant and high temperature resistance [11,12,13].
In the chemical industry, manufacturers and researchers use PTFE as a gasket and rod material [14]. PTFE is a biocompatible material and has no allergic reactions. Therefore, it has wide application in medicine for the production of prosthetic joint implants, prostheses, catheters, and other medical devices [15,16]. In mechanical engineering, PTFE and composites based on it are used in the production of bearings, shafts, seals, and other components that must operate under conditions of high temperature and pressure [17,18,19]. In its pure form, PTFE has a number of characteristics; however, its essential feature is its bonding with reinforcing elements. Fluoroplastic is mainly used as a matrix to create a new material. In particular, its gamma radiation modifications are the matrix for a wide range of composite materials [20]. In numerical simulation of the performance of polymeric materials and composites based on them, an important step is the correct choice of the approach to describe both the material’s properties and its behavior [21]. There are different approaches to building qualitative models of composites. One of them is to describe the properties of the matrix and the inclusions or reinforcing elements separately from one another [22,23].
Studies of the properties of polymers and composites based on them have existed since the appearance of polymers themselves. However, much of this research has focused on the study of these materials as elastic or elastic–plastic bodies. It was found that the data obtained in solving such problems are not sufficient to qualitatively describe the behavior of the material. This is due not only to the lack of description of viscous behavior but also to the presence of a description of plasticity. It is noted that polymers are generally divided into two types: viscoelastic–plastic and viscoelastic. In addition, if, in the first type there, is a trace of plastic deformations, the second type will have insignificant plastic deformations. Studies aimed at analyzing both the viscoelastic–plastic behavior of materials [24,25,26,27] and viscoelastic behavior of materials [28,29,30] are gaining momentum. The paper considers a mathematical model of only viscoelastic behavior for the first stage of the study.

1.3. Problem Description

Polymer materials have gained a strong foothold in the international market and are actively used on a large scale in various industrial fields. In particular, the materials are used as sliding layers in various types of friction units. However, there is a lack of research on the deformation behavior of polymers under various design conditions [31,32,33]. In particular, studies are most often directed towards a single configuration of a working structure. However, in order to qualitatively describe the behavior of the structure, it is necessary to consider its operation in a complex [34]. For this purpose, it is necessary to consider the problem not in a static formulation, but in a dynamic one, taking into account such factors as variable cyclic load, variable temperature, and different variations of geometrical configuration. However, to qualitatively describe and consider all these factors, it is necessary to choose the right description of the material behavior model.
This work aims to investigate the influence of the choice of a polymer behavior model. Gamma-modified PTFE is described in three ways: an elastic body, an elastic–plastic body, and a viscoelastic body. The model problem of introducing a spherical indenter into polymer half-space is realized.

2. Materials and Methods

2.1. Materials

Composite materials are widely spread in different spheres of human activity. However, creating and study the matrix [35,36] and composite fillers [37,38] is an important task. Matrices from different materials allow the structure to work in different conditions: increased and decreased temperatures [39], responsible friction nodes [40], aggressive environments [41], etc.
The article presents a study of one of the common matrices of composite materials in bridge-building activities—gamma-modified PTFE. The material has high strength and antifriction properties [42]. A full-scale experiments series were realized to determine the following: elastic characteristics, stress–strain dependence, friction properties, and dynamic characteristics for this material. Experimental studies were performed by Dr. Adamov A.A. using the equipment of the Ural Branch of the Russian Academy of Sciences. The experimental study was conducted using Zwick Z100SN5A (Zwick Roell AG, Ulm, Germany), which allows mechanical, quasi-static tests for uniaxial tension, compression, creep, etc. Test specimens are made in the form of cylindrical bodies (Figure 1) with characteristic dimensions: length, l = 20 ± 0.1 mm; radius, r = 10 ± 0.1 mm.
The specimen is deformed by Δ l up to 2 mm from the original size. The surface S 1 is rigidly fixed. Displacements are applied at a constant rate V Δ l = 0.006 mm/min on the surface S 2 . The paper has established that the material has no barrel-like deformation during experimental investigations. Consequently, only is realized in the z-axis direction within the experiment.
The elastic compression modulus E = 863.8 MPa was determined within the experiment at uniaxial deformation under constrained compression. The Poisson’s coefficient ν = 0.461 was determined using uniaxial deformation experiment.
The specimen was uniaxial compression stress tested ( ε = 10 %) to obtain the stress–strain relationship at (Figure 2).
The figure shows that the PTFE elastic behavior occurs when the sample is deformed up to 1%, after which the material enters plasticity. It was also observed that the material recovered its original geometric dimensions after some time. Therefore, plastic deformation is negligible in the polymer. From this, it can be concluded that the stress reduction is due to the viscoelastic behavior of the material.
The next step was to determine the dynamic characteristics of the material under uniaxial stress. One loading cycle was carried out:
-
The specimen was tested at a strain rate of 0.006 mm/min up to 10%;
-
Constant strain was held for 15 min;
-
The load was removed from the material until normal stresses of 0.1 MPa was reached at a rate of 0.006 mm/min.
Based on this, the following dependencies were obtained: strains on time (Figure 3a), normal stresses on time (Figure 3b), and stresses on strain (Figure 3c).
The obtained data will be used to describe the material behavior via the deformation theory of elastic–plasticity and the model of viscoelastic behavior based on the Prony series.

2.2. Description of Polymer Behavior

Currently, three types of material description are common: an elastic body, an elastic–plastic body [43], and a viscoelastic body [44]. The elastic–plastic body is described by the deformation theory of elastic–plasticity [45].
Many authors consider polymeric materials, matrices, and composites in the form of a Maxwell body [46,47]. The Prony series is the most common model for describing a Maxwell body [48].
The mathematical formulation of the problem includes the equilibrium equation:
div σ ^ = 0 ,
where σ ^ is the stress tensor.
The problem is considered in the framework of large deformations:
ε ^ = [ u ¯ + ( u ¯ ) T + u ¯ ( u ¯ ) T ] / 2 ,
where u ¯ is the displacement vector; and ε ^ is the strain tensor.
Let us write down the stress–strain relationship for each case when describing the polymer material by different models:
-
Elastic body:
σ ^ = λ I 1 ( ε ^ ) I ^ + 2 μ ε ^ ,
where λ and μ are Lame parameters; I 1 ( ε ^ ) is the first invariant of the stress tensor; and I ^ is a unit vector;
-
Elastoplastic body:
σ ^ = 2 σ I [ ε ^ I 1 ( ε ^ ) I ^ / 3 ] / ( 3 ε I ) + K I 1 ( ε ^ ) I ^ ,
where σ I and ε I are the stress and strain intensity, respectively; and K is the bulk modulus of elasticity;
-
Viscoelastic body:
σ ^ = 0 T [ E + E 0 i = 1 k α i exp ( ( T τ ) / β i ) ] d ε ( τ ) ,
where E 0 and E are Young’s moduli at the initial and final moment of time, respectively; α i is weighting coefficient; and β i is relaxation times.
A numerical procedure is created to describe the viscoelastic behavior of a polymer material (the scheme is presented in Figure 4). Experimental data are input to this procedure. The unknown coefficients are found by solving the minimization problem of function (6). The presented functional has a parabolic form. Therefore, it has one global minimum (one solution). They are necessary for a correct description of material behavior.
F = j = 1 n [ σ j E x p σ j N u m ( x ¯ ) ] 2 ,
where n is the number of experimental points; σ j E x p is the experimental stress value at j point; σ j n u m ( x ¯ ) is the numerical stress value at j point; and x ¯ = ( α i , β i ) is the vector of unknowns, which consists of 2 k parameters.
This numerical procedure is used to find the coefficients necessary to describe the viscoelastic behavior of the material. The solution is realized using the finite element method in ANSYS Mechanical APDL 2021R2 (Livermore, California, USA) software package in synergy with Python. During the operation of the numerical procedure, the coefficients are selected, then the experiment is numerically realized, and the results of numerical and in situ experiments are verified. The numerical procedure runs until the error is less than 5% between the experimental and numerical values.
The numerical procedure consists of 3 stages:
-
First (preliminary) stage: input of experimental data in the form of a text file with data, model selection, and generation of the initial vector of unknowns from Equation (5);
-
The second stage is based on the Nelder–Mead optimization algorithm: creation of a script-file describing the numerical experiment in ANSYS Mechanical APDL, conducting the numerical experiment, obtaining the results file, comparison of numerical and experimental data, generation of the vector of unknowns from Equation (5) at step i, and transition to the next iteration;
-
The third (and final) step is performed when the error between the experimental and numerical data reaches 5%: formation of the final vector of unknowns from Equation (5) and exit from the procedure.

2.3. The Hertz Formulation

The Hertz contact problem is one of the common problems of contact interaction between an indenter and a half-space [49]. Figure 5 shows the contact interaction scheme between a spherical indenter and a half-space.
The spherical indenter of radius R = 0.2 m is penetrated with force F = 1000 N into a half-space with geometrical characteristics: length is l p = 0.31 m; height is h p = 0.05 m.
All possible contact states at the site S K are considered within the problem. The contact boundary conditions are of the following forms:
-
Sliding friction: u n 1 = u n 2 , u τ 1 1 u τ 1 2 , u τ 2 1 u τ 2 2 , σ n 1 = σ n 2 , σ n τ 1 1 = σ n τ 1 2 , σ n τ 2 1 = σ n τ 2 2 , when | σ n τ 1 | = μ ( σ n ) | σ n | ;
-
No contact: | u n 1 u n 2 | 0 , σ n τ 1 = σ n τ 2 = σ n = 0 ;
-
Adhesion: u ¯ 1 = u ¯ 2 , σ n 1 = σ n 2 , σ n τ 1 1 = σ n τ 1 2 , σ n τ 2 1 = σ n τ 2 2 ,
where μ ( σ n ) is friction coefficient; τ 1 and τ 2 are the axes designations that lie in the plane tangent to the contact surface; u n are displacements along the normal to the corresponding contact boundary; u τ 1 and u τ 2 are displacements in the tangent plane; σ n is stress along the normal to the contact boundary; σ n τ 1 and σ n τ 2 are tangential stresses at the contact boundary; σ n τ is the value of the tangential contact stress vector; and 1 and 2 are conditional numbers of the contacting surfaces.
The Hertz solution is obtained for the case of a parabolic pressure profile and has the following form:
p ( r ) = p 0 ( 1 r 2 / a 2 ) 1 / 2 ,
where r is the distance for an arbitrary point on the plane; a is the Hertz contact radius; and p 0 is the maximum contact pressure.
This dependence will be used to analyze the convergence of the problem within the static problem formulation in a subsequent study.

3. Results

3.1. Invastigation of Mathematical Models

The mathematical model choice for describing material behavior is an important part of computer engineering as it affects the accuracy of the results obtained in the study. This article compared experimental data and mathematical models at the first stage of the study.
The numerical procedure (Section 2.2) allows us to find the vector of unknowns for the Prony series with an error of less than 5%. The final vector of unknowns is presented in Figure 6.
It can be noted that the weight coefficients α i have the largest values in the values range of relaxation times β i [ 10 2 ; 10 4 ].
Numerical modeling of a cylindrical specimen uniaxial deformation has been conducted—similarly to the field experiment (paragraph 2.1.). Uniaxial deformation occurs up to a strain value equal to 10%. The problem was solved in a static formulation and let us obtain the stress–strain diagram for all considered variants of body behavior (Figure 7).
Elastic and viscoelastic bodies behave linearly over the entire deformation interval of the specimen when solving the static problem. At the same time, the elastic–plastic body describes the deformation of the specimen up to 10%, with an error of less than 5%. Consequently, the use of the mathematical model of the elastic–plastic body allows us to explore problems in analyzing the strength of a structure within the framework of static calculations.
Moreover, predicting the performance of the structure at all stages of its life cycle is an important feature of numerical modeling. Let us perform numerical modeling of the experiment depending on the time of the load impact on the specimen (Figure 8).
The article notes that when comparing the numerical solution with the experimental one, it does not fully describe the behavior of the material and has a linear character over the entire range of deformation in the framework of the theory of elasticity. In the framework of elastic–plasticity and viscoelasticity theory, the time dependencies of stresses present a more qualitative description of material behavior. The difference from the experimental data in the elastic–plastic model is, at the loading stage, 5%, and at the unloading stage, 20%. This is due to significant plastic deformations, due to which there is an extreme drop in stress. In the case of the viscoelastic body the difference from experimental data is, at loading stage, 1%, and at the unloading stage, 0.53%. Consequently, the description of polymer behavior by the viscoelastic model allows for a better description of its dynamic behavior.
For a qualitative assessment of the material behavior, we show the stress–strain dependence in Figure 9.
Similar to the time dependence of stresses, a significant difference between the elastic body and all others is noted. The elastic body behaves linearly during the loading and unloading phases of the numerical model. It should be noted that for an elastic–plastic body at unloading, the stress values decrease linearly, which is due to the accumulated plastic deformation in the body. There is no zone of stress reduction at load holding. When considering the viscoelastic model, it is worth noting the qualitative and quantitative description of material behavior in the dynamic formulation.
The following conclusions can be drawn from the above: it is necessary to use the elastic–plastic model to analyze the structure strength in the static setting, and it is necessary to use the viscoelastic model to predict the performance of the structure during the life cycle.

3.2. Hertz Contact Calculation Model

The model problem of indenter penetration into a half-space is solved to practice the use of mathematical models in the structure operation.
The first study step is to determine the optimal value of the mesh size within the framework of computer engineering. For this purpose, we conduct a series of numerical experiments to determine the optimal finite element partition size (Figure 10). The size of the finite element near the contact is chosen as the variable parameter h e .
When analyzing the finite element partitioning, it can be seen that the largest error occurs in the leftmost and rightmost contact nodes. However, when the finite element size is reduced, the numerical solution approaches (7). An error of less than 1% is achieved with finite element partitioning h e = 0.03125 mm .
Further, we realize the problem in dynamic formulation for two models: the elastic–plastic body and the viscoelastic body. The paper also considers the value of contact pressure at maximum load and its endurance for 1 h (Figure 11).
An elastic–plastic body does not change with time. In a viscoelastic body, the following are observed: a decrease in contact pressure values, and an increase in the contact area between the spherical indenter and the half-space.
Further, the article considers the distribution of maximum strain values as a function of time (Figure 12).
It should be noted that the initial point of the strain intensity plots coincide for elastic–plastic and viscoelastic bodies. However, with the time passage, the values of strain intensity in the elastic–plastic body are constant. At the same time, the value of the viscoelastic body grows nonlinearly. From the above, we can conclude that the creep of the material occurs with the passage of time.
Let us also consider the maximum values distribution of stress intensity (Figure 13).
Similar to deformations, the stress intensity has a different distribution of values over time depending on the body type. An elastic–plastic body has a constant value over the entire time range. At the same time, for a viscoelastic body, the level of stress intensity decreases with time due to the material creep (Figure 12).

4. Discussion

4.1. Limitation Statement

The paper presents the results of numerical identification and simulation of the behavior of gamma-modified PTFE. The work has a number of limitations that are planned to be eliminated in the future:
  • The material behavior is considered at a constant temperature of 20 °C;
  • The model problem of spherical indenter penetration into a half-space is considered;
  • For each material, it is necessary to carry out a separate description of the mathematical model;
  • Long time ranges are not considered, while the material works for a long time.
Further directions for the development of the work:
  • Investigation of the material on a large range of operating temperatures;
  • Study of the material on the dependence on the load impact rate on the polymer material;
  • Study of temperature characteristics of the material;
  • Realization of the problem on the example of a bridge bearing structure under cyclic loading.
In the future, we plan to proceed to consideration of a real structure within the framework of contact interaction problems. The real structure consists of two steel plates: one with a spherical indenter and one with a spherical notch, and there is a sliding layer between them. The Hertz model was chosen as a simplified model for the initial verification of the obtained results on the viscoelastic behavior of the polymer material.
When considering the model of a real structure, we plan to simulate its operation under cyclic loads, as well as at different temperature parameters. In this case, we plan to obtain a description of deformation behavior as close as possible to the behavior of a real structure in which irreversible deformations occur over time.

4.2. On the Choice of a Mathematical Model

There is a certain variety of descriptions for the viscoelastic–plastic behavior of a material; in particular, there are such models as the cooperative-viscoplasticity theory based on overstress (VBO) model [16], combining a nonlinear viscoelastic model with a viscoplastic model using the von Mises yield criterion [12,24,50,51,52], etc. However, in the present work, gamma-modified PTFE is considered. In a number of experiments, it is noted that there is no need to describe the behavior of the material as a viscoelastic–plastic body; a viscoelastic model is sufficient for a qualitative description of its behavior. Within the framework of the first approximation, the Prony viscoelastic model was chosen.
As an example, consider the work [53], in which an experimental and numerical study of a Prony series as the main relation of the matrix description was conducted. It is shown that when the number of unknown terms increases, the numerical data describe the experimental data with an error of 10%.
This model has a number of advantages: a sufficient description of material behavior, the active application of this model by other researchers [29,54], the relative simplicity of the mathematical description, etc. The numerical algorithm has been tested early for the description of lubricant behavioral models [55].

4.3. Applicability of the Research

Gamma-modified PTFE is widely used in bridge bearings as an antifriction material [40,45,56]. In this design, the material operates under large temperature differences (from −40 °C to +40 °C) and cyclic loads directed in different directions; thus, numerical experiments in dynamic formulation are necessary for qualitative prediction of the load-bearing structure performance. It was obtained that the elastic–plastic description of materials is suitable only for describing the performance of the structure in the static formulation; thus, to describe the material behavior over time, it is necessary to describe it in a viscoelastic formulation. In the future, it is planned to transfer to a model of a real bridge support structure with a different set of temperature and cyclic tests.

5. Conclusions

Description of material behavior is an important part of computational engineering research. Its correct description allows us to qualitatively predict the performance of a structure during its life cycle. Within the framework of this work, gamma-modified PTFE has been described in three ways: an elastic body, an elastic–plastic body, and a viscoelastic body. Within the numerical analysis, the following has been established:
-
The use of an elastic–plastic body to describe the material behavior can be used only in static problems to determine the strength properties of the structure;
-
The use of a viscoelastic body to describe the mathematical model of material behavior allows for evaluation of the performance of a structure at the entire stage of its life cycle.
The model problem of introducing a spherically shaped indenter into polymer half-space is realized.
The research presented in the article will draw attention to the necessity of correct and detailed descriptions of mathematical models of material behavior. This, in turn, will allow for improvements to the quality of solved problems within the computer engineering framework. At the same time, at early stages, it will be possible to track problem areas of the structure operation and perform manipulations for their elimination in a timely manner.
The presented study can be used in compression contact nodes; for realization in tensile problems, it is necessary to undertake additional research.

Author Contributions

Conceptualization, A.A.K., Y.O.N. and A.P.B.; methodology, A.A.K.; software, A.A.K., Y.O.N. and A.P.B.; validation, A.A.K.; writing—original draft preparation, Y.O.N. and A.P.B.; writing—review and editing, Y.O.N. and A.P.B.; visualization, A.A.K., Y.O.N. and A.P.B.; funding acquisition, A.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by a grant from the Russian Science Foundation (project No. 22-29-01313).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bernard, C.A.; Correia, J.P.M.; Ahzi, S.; Bahlouli, N. Numerical implementation of an elastic-viscoplastic constitutive model to simulate the mechanical behaviour of amorphous polymers. Int. J. Mater. Form. 2017, 10, 607–621. [Google Scholar] [CrossRef]
  2. Hartmann, S.; Gilbert, R.R. Identifiability of material parameters in solid mechanics. Arch. Appl. Mech. 2018, 88, 3–26. [Google Scholar] [CrossRef]
  3. Varna, J.; Pupure, L. Effect of Material State and Temperature on Nonlinear Viscoelastic Response: 3D Constitutive Model and Incremental Formulation for Numerical Analysis. Mech. Compos. Mater. 2023, 59, 193–218. [Google Scholar] [CrossRef]
  4. Yin, B.; Kaliske, M. Fracture simulation of viscoelastic polymers by the phase-field method. Comput. Mech. 2020, 65, 293–309. [Google Scholar] [CrossRef]
  5. Forintos, N.; Czigany, T. Multifunctional application of carbon fiber reinforced polymer composites: Electrical properties of the reinforcing carbon fibers—A short review. Compos. Part B Eng. 2018, 162, 331–343. [Google Scholar] [CrossRef]
  6. Umar, K.; Yaqoob, A.A.; Ibrahim, M.N.M.; Parveen, T.; Safian, M.T.U. Chapter Thirteen-Environmental applications of smart polymer composites. Smart Polym. Nanocomposites 2021, 295–312. [Google Scholar] [CrossRef]
  7. Bakbak, O.; Colak, O. Modeling wide range of viscoelastic–viscoplastic behavior of Araldite LY 564 epoxy using cooperative viscoplasticity theory based on overstress model. J. Appl. Polym. Sci. 2023, 140, e53557. [Google Scholar] [CrossRef]
  8. Yao, B.; Hong, W.; Chen, T.; Han, Z.; Xu, X.; Hu, R.; Hao, J.; Li, C.; Li, H.; Perini, S.E.; et al. Highly Stretchable Polymer Composite with Strain-Enhanced Electromagnetic Interference Shielding Effectiveness. Adv. Mater. 2020, 32, e1907499. [Google Scholar] [CrossRef]
  9. Van Meter, K.E.; Junk, C.P.; Campbell, K.L.; Babuska, T.F.; Krick, B.A. Ultralow Wear Self-Mated PTFE Composites. Macromolecules 2022, 55, 3924–3935. [Google Scholar] [CrossRef]
  10. Peng, S.; Zhang, L.; Xie, G.; Guo, Y.; Si, L.; Luo, J. Friction and wear behavior of PTFE coatings modified with poly (methyl methacrylate). Compos. Part B Eng. 2019, 172, 316–322. [Google Scholar] [CrossRef]
  11. Mashtalyar, D.V.; Pleshkova, A.I.; Piatkova, M.A.; Nadaraia, K.V.; Imshinetskiy, I.M.; Belov, E.A.; Suchkov, S.N.; Sinebryukhov, S.L.; Gnedenkov, S.V. PTFE-Containing Coating Obtained on Ti by Spraying and PEO Pretreatment. Coatings 2023, 13, 1249. [Google Scholar] [CrossRef]
  12. Gooneratne, R.; Iroh, J.O. Thermomechanical and Pre-Ignition Properties of Multicomponent Poly(Vnylidene Fluoride)/Aluminum Oxide/Single-Walled Carbon Nanotube Hybrid Nanocomposites. J. Compos. Sci. 2022, 6, 380. [Google Scholar] [CrossRef]
  13. Zhang, G.; Yin, J.; Yan, M.; Zhang, H.; Yin, X. Improved Through-Plane Thermal Conductivity and Mechanical Properties of Low-Dielectric FEP/HBN Composites by Adding PTFE Nanoparticles. ACS Appl. Electron. Mater. 2021, 3, 4568–4578. [Google Scholar] [CrossRef]
  14. Stoporev, A.; Kadyrov, R.; Adamova, T.; Statsenko, E.; Nguyen, T.H.; Yarakhmedov, M.; Semenov, A.; Manakov, A. Three-Dimensional-Printed Polymeric Cores for Methane Hydrate Enhanced Growth. Polymers 2023, 15, 2312. [Google Scholar] [CrossRef] [PubMed]
  15. Rezvova, M.A.; Klyshnikov, K.Y.; Gritskevich, A.A.; Ovcharenko, E.A. Polymeric Heart Valves Will Displace Mechanical and Tissue Heart Valves: A New Era for the Medical Devices. Int. J. Mol. Sci. 2023, 24, 3963. [Google Scholar] [CrossRef] [PubMed]
  16. Bembenek, M.; Popadyuk, O.; Shihab, T.; Ropyak, L.; Uhryński, A.; Vytvytskyi, V.; Bulbuk, O. Optimization of Technological Parameters of the Process of Forming Therapeutic Biopolymer Nanofilled Films. Nanomaterials 2022, 12, 2413. [Google Scholar] [CrossRef] [PubMed]
  17. Tretiak, O.; Kritskiy, D.; Kobzar, I.; Arefieva, M.; Nazarenko, V. The Methods of Three-Dimensional Modeling of the Hydrogenerator Thrust Bearing. Computation 2022, 10, 152. [Google Scholar] [CrossRef]
  18. Bakkara, A.; Sadykov, B.; Zhapekova, A.; Oserov, T.; Batkal, A.; Khairullina, A.; Mofa, N. Efficiency and Prospects of the Use of Mechanochemical Treatment to Obtain Innovative Composite Systems. Chemengineering 2022, 6, 90. [Google Scholar] [CrossRef]
  19. Tretiak, O.; Kritskiy, D.; Kobzar, I.; Arefieva, M.; Selevko, V.; Brega, D.; Maiorova, K.; Tretiak, I. Stress-Strained State of the Thrust Bearing Disc of Hydrogenerator-Motor. Computation 2023, 11, 60. [Google Scholar] [CrossRef]
  20. Mylsamy, G.; Krishnasamy, P. A Review on Electrical Properties of Fiber-Reinforced Polymer Material: Fabrication, Measurement, and Performances. Trans. Indian Inst. Met. 2023, 76, 1691–1708. [Google Scholar] [CrossRef]
  21. Yang, G.; Meng, G.; Gao, H.; Lin, Q. Micromorphology and mechanical properties of UHMWPE/CNF composites under accelerated aging. Polym. Compos. 2022, 43, 3716–3730. [Google Scholar] [CrossRef]
  22. Panaitescu, D.M.; Vuluga, Z.; Frone, A.N.; Gabor, A.R.; Nicolae, C.-A.; Uşurelu, C.-D. Complex Effects of Hemp Fibers and Impact Modifiers in Multiphase Polypropylene Systems. Polymers 2023, 15, 409. [Google Scholar] [CrossRef] [PubMed]
  23. Zhu, R.; Li, X.; Wu, C.; Du, L.; Du, X. The In-Situ Mechanical Properties of Carbon Fiber/Epoxy Composite under the Electric-Current Loading. Polymers 2022, 14, 4452. [Google Scholar] [CrossRef] [PubMed]
  24. Doh, J.; Hur, S.H.; Lee, J. Viscoplastic parameter identification of temperature-dependent mechanical behavior of modified polyphenylene oxide polymers. Polym. Eng. Sci. 2019, 59, E200–E211. [Google Scholar] [CrossRef]
  25. Liu, E.; Wu, J.; Li, H.; Liu, H.; Xiao, G.; Shen, Q.; Kong, L.; Lin, J. Research on viscoelastic behavior of semi-crystalline polymers using instrumented indentation. J. Polym. Sci. 2021, 59, 1795–1806. [Google Scholar] [CrossRef]
  26. Alam, I.; Khan, D.; Mittal, Y.; Kumar, S. Effect of crack tip shape on near-tip deformation and fields in plastically compressible solids. J. Braz. Soc. Mech. Sci. Eng. 2019, 41, 441. [Google Scholar] [CrossRef]
  27. O’connor, J.; dos Santos, B.B.; Borges, L.; da Costa, M.F.; Castello, D.A. Computational modeling of viscoplastic polymeric material response during micro-indentation tests. J. Braz. Soc. Mech. Sci. Eng. 2020, 42, 438. [Google Scholar] [CrossRef]
  28. Nunes, S.G.; Saseendran, S.; Joffe, R.; Amico, S.C.; Fernberg, P.; Varna, J. On Temperature-Related Shift Factors and Master Curves in Viscoelastic Constitutive Models for Thermoset Polymers. Mech. Compos. Mater. 2020, 56, 573–590. [Google Scholar] [CrossRef]
  29. Chae, S.-H.; Zhao, J.-H.; Edwards, D.R.; Ho, P.S. Characterization of the Viscoelasticity of Molding Compounds in the Time Domain. J. Electron. Mater. 2010, 39, 419–425. [Google Scholar] [CrossRef]
  30. Freire, R.T.S.; Nunes, S.G.; Amico, S.C.; Al-Ramahi, N.J.; Joffe, R.; Varna, J. On Determination of the Linear Viscoelastic Compliance and Relaxation Functions for Polymers in One Tensile Test. Mech. Compos. Mater. 2023, 58, 765–786. [Google Scholar] [CrossRef]
  31. Gonabadi, H.; Oila, A.; Yadav, A.; Bull, S. Investigation of the Effects of Environmental Fatigue on the Mechanical Properties of GFRP Composite Constituents Using Nanoindentation. Exp. Mech. 2022, 62, 585–602. [Google Scholar] [CrossRef]
  32. Wan, Y.; Yu, Y.; Cao, L.; Zhang, M.; Gao, J.; Qi, C. Corrosion and tribological performance of PTFE-coated electroless nickel boron coatings. Surf. Coatings Technol. 2016, 307, 316–323. [Google Scholar] [CrossRef]
  33. Sajid, M.; Ilyas, M. PTFE-coated non-stick cookware and toxicity concerns: A perspective. Environ. Sci. Pollut. Res. 2017, 24, 23436–23440. [Google Scholar] [CrossRef] [PubMed]
  34. Vahabi, H.; Laoutid, F.; Mehrpouya, M.; Saeb, M.R.; Dubois, P. Flame retardant polymer materials: An update and the future for 3D printing developments. Mater. Sci. Eng. R Rep. 2021, 144, 100604. [Google Scholar] [CrossRef]
  35. El-Sheikh, M.N.; Metwally, B.S.; Mubarak, M.F.; Ahmed, H.A.; Moghny, T.A.; Zayed, A.M. Fabrication of electrospun polyamide–weathered basalt nano-composite as a non-conventional membrane for basic and acid dye removal. Polym. Bull. 2023, 80, 8511–8533. [Google Scholar] [CrossRef]
  36. Murtada, K.; Nazdrajić, E.; Pawliszyn, J. Polybenzimidazole: A novel, fluorocarbon-free, SPME sorbent binder with good thermal and solvent resistance properties for GC and LC analysis. Microchim. Acta 2023, 190, 323. [Google Scholar] [CrossRef] [PubMed]
  37. Mohammed, M.I.; Yahia, I.S.; El-Sayed, F. Study the microstructure, optical, mechanical, and dielectric characteristics of Cu2+/Zn2+-codoped PVA for CUT-OFF laser filters. Polym. Bull. 2023, 80, 8377–8405. [Google Scholar] [CrossRef]
  38. Chulliyote, R.; Hareendrakrishnakumar, H.; Kannan, S.K.; Joseph, M.G. Biomass-derived inherently doped multifunctional hierarchically porous carbon as an efficient electrode material for high-performance supercapacitors. J. Porous Mater. 2023, 30, 1129–1141. [Google Scholar] [CrossRef]
  39. Pavlov, A.A.; Soboleva, E.S. Anti-Frictional Copper-Fluoroplast Coating for Parts of Power Equipment. Chem. Petrol. Eng. 2021, 57, 56–59. [Google Scholar] [CrossRef]
  40. Adamov, A.A.; Kamenskikh, A.A.; Pankova, A.P.; Strukova, V.I. Comparative Analysis of the Work of Bridge Spherical Bearing at Different Antifriction Layer Locations. Lubricants 2022, 10, 207. [Google Scholar] [CrossRef]
  41. Zhu, G.; Zhang, C.; Li, K.; Xi, X.; Zhang, X.; Lei, H. A multifunctional zeolitic imidazolate framework-8@wood aerogel composite intergrating superior performance of dye adsorption capacity and flame-retardant property. J. Porous Mater. 2023, 30, 1171–1182. [Google Scholar] [CrossRef]
  42. Adamov, A.A.; Kamenskih, A.A.; Pankova, A.P. Numerical analysis of the spherical bearing geometric configuration with antifriction layer made of different materials. PNRPU Mech. Bull. 2020, 4, 15–26. [Google Scholar] [CrossRef]
  43. Chebakov, M.I.; Danil’chenko, S.A. Theoretical-Experimental Approach to Determining the Mechanical Characteristics of Thin Polymer Composite Coatings. Mech. Compos. Mater. 2023, 58, 803–810. [Google Scholar] [CrossRef]
  44. Polychronopoulos, N.; Vlachopoulos, J. Polymer Processing and Rheology. Funct. Polym. 2019, 1–47. [Google Scholar] [CrossRef]
  45. Nosov, Y.O.; Kamenskikh, A.A. Influence Analysis of Lubricant Recesses on the Working Capacity of the Bridge Span Spherical Bearing. Lubricants 2022, 10, 283. [Google Scholar] [CrossRef]
  46. Khokhlov, A.V.; Shaporev, A.V.; Stolyarov, O.N. Loading-Unloading-Recovery Curves for Polyester Yarns and Identification of the Nonlinear Maxwell-Type Viscoelastoplastic Model. Mech. Compos. Mater. 2023, 59, 129–146. [Google Scholar] [CrossRef]
  47. Khokhlov, A.V. Generalization of a Nonlinear Maxwell-Type Viscoelastoplastic Model and Simulation of Creep and Recovery Curves. Mech. Compos. Mater. 2023, 59, 441–454. [Google Scholar] [CrossRef]
  48. Lesnikova, Y.I.; Trufanov, A.N.; Kamenskikh, A.A. Analysis of the Polymer Two-Layer Protective Coating Impact on Panda-Type Optical Fiber under Bending. Polymers 2022, 14, 3840. [Google Scholar] [CrossRef]
  49. Liu, T.; Fan, J.; Peng, Z. Central Load-Bearing Control in the Construction Process of the Concrete Spherical Joint Nandu River Swing Bridge: A Case Study. Buildings 2022, 12, 511. [Google Scholar] [CrossRef]
  50. Chen, Y.; Smith, L.V. A nonlinear viscoelastic–viscoplastic constitutive model for adhesives under creep. Mech. Time-Depend. Mater. 2022, 26, 663–681. [Google Scholar] [CrossRef]
  51. Sofiyev, A. On the Solution of Dynamic Stability Problem of Functionally Graded Viscoelastic Plates with Different Initial Conditions in Viscoelastic Media. Mathematics 2023, 11, 823. [Google Scholar] [CrossRef]
  52. Reza, A.; Shishesaz, M.; Sedighi, H.M. Micromechanical approach to viscoelastic stress analysis of a pin-loaded hole in unidirectional laminated PMC. Polym. Polym. Compos. 2021, 29, S1144–S1157. [Google Scholar] [CrossRef]
  53. Reichl, K.K.; Avci, O.; Inman, D.J. Temperature Dependent Damping in Additively Manufactured Polymer Structures. J. Appl. Comput. Mech. 2021, 7, 993–1008. [Google Scholar] [CrossRef]
  54. Endo, V.T.; de Carvalho Pereira, J.C. Linear orthotropic viscoelasticity model for fiber reinforced thermoplastic material based on Prony series. Mech. Time-Depend. Mater. 2017, 21, 199–221. [Google Scholar] [CrossRef]
  55. Nosov, Y.O.; Kamenskikh, A.A. Experimental Study of the Rheology of Grease by the Example of CIATIM-221 and Identification of Its Behavior Model. Lubricants 2023, 11, 295. [Google Scholar] [CrossRef]
  56. Adamov, A.A.; Kamenskikh, A.A.; Pankova, A.P. Influence Analysis of the Antifriction Layer Materials and Thickness on the Contact Interaction of Spherical Bearings Elements. Lubricants 2022, 10, 30. [Google Scholar] [CrossRef]
Figure 1. Experimental studies: (a) experiment scheme; (b) full-scale sample. Sections 1 and 2 are upper and lower steel sections, respectively; 3—polymer.
Figure 1. Experimental studies: (a) experiment scheme; (b) full-scale sample. Sections 1 and 2 are upper and lower steel sections, respectively; 3—polymer.
Polymers 15 03630 g001
Figure 2. The σ ε diagram for gamma-modified PTFE.
Figure 2. The σ ε diagram for gamma-modified PTFE.
Polymers 15 03630 g002
Figure 3. Dynamic characteristics determination of the material: (a) strain–time dependence, (b) stress–time dependence; (c) stress–strain dependence.
Figure 3. Dynamic characteristics determination of the material: (a) strain–time dependence, (b) stress–time dependence; (c) stress–strain dependence.
Polymers 15 03630 g003
Figure 4. The numerical procedure scheme.
Figure 4. The numerical procedure scheme.
Polymers 15 03630 g004
Figure 5. Scheme of the Hertz contact problem.
Figure 5. Scheme of the Hertz contact problem.
Polymers 15 03630 g005
Figure 6. Dependence of weight coefficients α i on relaxation time β i . Red dots are the value of the weight coefficient α i at a certain relaxation time β i .
Figure 6. Dependence of weight coefficients α i on relaxation time β i . Red dots are the value of the weight coefficient α i at a certain relaxation time β i .
Polymers 15 03630 g006
Figure 7. Stress–strain diagram: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; and the dots are the viscoelastic body.
Figure 7. Stress–strain diagram: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; and the dots are the viscoelastic body.
Polymers 15 03630 g007
Figure 8. Stress dependence on time: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; the dots are the viscoelastic body.
Figure 8. Stress dependence on time: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; the dots are the viscoelastic body.
Polymers 15 03630 g008
Figure 9. Stress dependence on deformation: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; and the dots are the viscoelastic body.
Figure 9. Stress dependence on deformation: the red line is experimental data; the black solid line is the elastic body; the dashed line is the elastic–plastic body; and the dots are the viscoelastic body.
Polymers 15 03630 g009
Figure 10. Analysis of finite element partitioning of the Hertz contact.
Figure 10. Analysis of finite element partitioning of the Hertz contact.
Polymers 15 03630 g010
Figure 11. Dependence of contact pressure on the radius of indenter insertion: the red line is the analytical solution; black lines are maximum load; gray lines are 1 h exposure time; the dashed line is the elastic–plastic body; the dots are the viscoelastic body.
Figure 11. Dependence of contact pressure on the radius of indenter insertion: the red line is the analytical solution; black lines are maximum load; gray lines are 1 h exposure time; the dashed line is the elastic–plastic body; the dots are the viscoelastic body.
Polymers 15 03630 g011
Figure 12. Dependence of maximum values of strain intensity on time: the dotted line is the elastic–plastic body; the dots are the viscoelastic body.
Figure 12. Dependence of maximum values of strain intensity on time: the dotted line is the elastic–plastic body; the dots are the viscoelastic body.
Polymers 15 03630 g012
Figure 13. Dependence of maximum values of stress intensity on time: the dotted line is the elastic–plastic body, the dots are the viscoelastic body.
Figure 13. Dependence of maximum values of stress intensity on time: the dotted line is the elastic–plastic body, the dots are the viscoelastic body.
Polymers 15 03630 g013
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.

Share and Cite

MDPI and ACS Style

Kamenskikh, A.A.; Nosov, Y.O.; Bogdanova, A.P. The Study Influence Analysis of the Mathematical Model Choice for Describing Polymer Behavior. Polymers 2023, 15, 3630. https://doi.org/10.3390/polym15173630

AMA Style

Kamenskikh AA, Nosov YO, Bogdanova AP. The Study Influence Analysis of the Mathematical Model Choice for Describing Polymer Behavior. Polymers. 2023; 15(17):3630. https://doi.org/10.3390/polym15173630

Chicago/Turabian Style

Kamenskikh, Anna A., Yuriy O. Nosov, and Anastasia P. Bogdanova. 2023. "The Study Influence Analysis of the Mathematical Model Choice for Describing Polymer Behavior" Polymers 15, no. 17: 3630. https://doi.org/10.3390/polym15173630

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop