CC BY
2
EDN: BQKMXB УДК 539.3
Computer-aided Design of Polymer-based Composites Using Multi-scale Models
Sergey V. Panin* Boris A. Lyukshin^
Institute of Strength Physics and Materials Science SB RAS National Research Tomsk Polytechnic University Tomsk, Russian Federation
Nataliya Y. Grishaeva*
Tomsk State University of Control Systems and Radioelectronics
Tomsk, Russian Federation
Svetlana A. Bochkareva^ Iliya L. Panov^
Institute of Strength Physics and Materials Science SB RAS Tomsk State University of Control Systems and Radioelectronics National Research Tomsk State University Tomsk, Russian Federation
Received 10.10.2023, received in revised form 07.11.2023, accepted 24.12.2023
Abstract. The study is devoted to computer-aided design of particulate polytetrafluoroethylene-based composites for increasing their thermal conductivity values (as much as possible) with minimal deterioration of the mechanical properties. To solve the issue, an algorithm has been developed and applied using a multi-scale model. It has enabled to fabricate such composites, characterized by both high thermal conductivity and acceptable mechanical properties, via loading the polymer matrix with copper particles of two size ranges, differing by more than an order of magnitude.
Keywords: multi-scale model, polymer composite, particulate filler, stress-strain state, thermal conductivity, polytetrafluoroethylene (PTFE).
Citation: S.V. Panin, B.A. Lyukshin, N.Y. Grishaeva, S.A. Bochkareva, I.L. Panov, Computer-aided Design of Polymer-based Composites Using Multi-scale Models, J. Sib Fed. Univ. Math. Phys., 2024, 17(1), 65-74. EDN: BQKMXB.
Introduction
Polymer-based composites have become widespread because of the possibility to vary their functional characteristics over a wide range [1,2]. As a rule, loading with fillers, including particulate ones, changes their both structure and properties but does not affect the chemical nature of the polymer matrix [3]. For tribological applications, numerous composites are based on polytetrafluoroethylene (PTFE) [4], which is an advanced high-molecular polymer, suitable for manufacturing of various products [5,6]. In particular, PTFE-based composites are used for the fabrication of seals in flange connections, since they meet all requirements for such components. They include prolonged service life, great resistance to wear, negative environmental impacts, as
* svp@ispms.ru https://orcid.org/0000-0001-7623-7360 t lba2008@ yandex.ru
* anohina@mail2000.ru
§ svetlanab7@yandex.ru ^ panov.iliya@mail.ru © Siberian Federal University. All rights reserved
well as chemical, aggressive and abrasive substances. At the same time, PTFE has a very low friction coefficient [7,8]. Its loading with graphite enhances the seal durability by many times due to the formation of a strengthened transfer film on the surface of tribological counterparts (shafts, rods, gears, fittings, etc.) [9-11]. The functional characteristics of PTFE are determined by its specific chemical and supermolecular structure [12-15], among which the most outstanding are high thermal resistance at temperatures up to 260 °C, as well as flexibility and elasticity in the range from -70 to +270 °C [16]. The PTFE-based composites are also applied for the manufacturing of piston rings for gas compressors, sealing gaskets, rings and cuffs for connecting equipment units, operating both stationary and at relatively low sliding speeds in the chemical, as well as oil and gas industries. The required set of the functional characteristics of the PTFE-based composites for various operating conditions is ensured by variations in their composition and production routes [17]. As most polymers, PTFE possesses low thermal conductivity that does not enable to dissipate heat properly when seals are operated at elevated temperatures. One of the methods to solve this issue is the addition of fillers, primarily metal powders. In this case, rising the dimensions of their particles substantially increase this parameter for such composites [18], but simultaneously deteriorates their mechanical properties, since such inclusions act as stress concentrators under applied loads. Ceteris paribus, stress peaks near inclusions contribute to the initiation of cracks in the composites, leading to their preliminary failure. It should be noted that the stress peaks enhance with increasing the dimensions of particulate filler. Thus, the problem arises of searching for a compromise between increasing thermal conductivity of the PTFE-based composites and reducing their mechanical properties. This research proposes an approach for solving this problem, considered in relation to the design of ones, characterized by both increased thermal conductivity and acceptable mechanical properties.
1. Problem statement and research method
The objects of the study were the PTFE-based composites that included the F-4 polymer matrix (according to the Russian state standard GOST 10007-80) and loaded with two industrially available fine copper powders with particle sizes of 1.5 pm (met the technical specification TU 1793-011-50316079-2004, designated as "filler 1") and 60 pm (met the Russian state standard GOST 4960-2009, designated as "filler 2") [19]. The key physical properties of the composites components in Tab. 1.
Table 1. The key physical properties of the composites components
Filler 1 Filler 2
Parameter Matrix (PTFE) (1.5 pm (60 pm
copper powder) copper powder)
Elastic modulus, MPa 410 129800 110000
Poisson's ratio 0.37 0.32 0.32
Yield point, MPa 11.8 - -
Ultimate tensile strength, MPa 20 - -
Stress at break, MPa 27 314 -
Elongation at break, % 250-500 - -
Density, kg/m3 2200 8960 8930
Specific heat capacity, J/(kg °C) 1040 - 385
Thermal conductivity, W/(m °C) 0.200 - 401
Computer-aided design of the PTFE-based composites with the required properties was carried out as follows. Initially, it was necessary to find a combination of the control parameters that would ensure the maximum thermal conductivity value (> 0.3 W/m °C) and elongation at break
of > 0.02. To solve this problem, an approach proposed the authors of [20-22] was modified in relation to loading the polymer matrix with particles whose sizes differed by more than an order of magnitude. An algorithm of computer-aided design of the composites with the required set of properties is shown in Fig. 1. At the first stage of computer-aided design, both mechanical and thermophysical properties were assessed for the composites. In this way, both stress-strain state (SSS) parameters and thermal conductivity values were calculated for various applied loads and thermal fluxes using the equations of mechanics of a deformable solid and solving thermal conductivity problems. In general, the functional characteristics, the filler sizes and shapes, the filling degree and other parameters determining both composition and properties of the composites could be optimized. In the studied cases, the filling degree was varied for both loaded copper powders. The problems were solved using the Abaqus CAE software package based on the finite element method.
Fig. 1. The algorithm of computer-aided design of the PTFE-based composites with the required set of properties
At the second stage, inverse problems of computer-aided design were solved, generally aimed at manufacturing the PTFE-based composites with the specified properties. A solution was based on constructing response surfaces of the functional characteristics to the control parameters in a space of states. To achieve this goal, a number of reference points, obtained from solving linear equations, were used. The presence of particles with sizes differing by more than an order of magnitude complicated the considered algorithm (Fig. 1) when solving the direct problems, since they involved several structural levels. Their number (two for the studied cases) depended on the size of fillers [22,23]. At the first level, solving the direct problems made it possible to calculate the functional characteristics of the matrix loaded with filler 1 (1.5 pm). Then, such a composite was considered as the matrix loaded with filler 2 (60 pm) at the second level. As a result, the functional characteristics were assessed for the composites as a whole. An implemented multi-scale algorithm is presented in Fig. 2.
Fig. 3 shows computational domains implemented for determining the SSS parameters under uniaxial tension and for solving the thermal conductivity problems in relation to the PTFE-based composites loaded with copper particles of different sizes.
At the first structural level, the rectangular computational domain was characterized by a side of 25 pm and particles with a radius of 1.5 pm (Fig. 3, right). The bottom edge of the domain
The multi-scale determination of the equivalent characteristics of the composites
Solving a number of direct problems of both SSS and thermal conductivity parameters (by varying tlie functional characteristics, the filler sizes and shapes, the volume fraction, etc.) at the first and second structural levels
V v
Calculation of SSS at the first structural level Calculation of thermal conductivity at the first
(by varying the volume fraction for filler 1) structural level (by varying the volume fraction for filler 1)
v v
Calculation of a stress-strain curve, the elastic Calculation of thermal conductivity
modulus and elongation at break at the first structural level
at the first structural level depending on the volume fraction for filler 1
depending on the volume fraction for filler 1
V V
Calculation of SSS at the second structural level Solution of the thermal conductivity problem at the
(by varying the SSS parameters obtained at the first structural level and the volume fraction for filler 2) second structural level (by varying the thermal conductivity values obtained at the first structural level
and the volume fraction for filler 2)
v v
Calculation elastic moduli values Calculation thermal conductivity values
for the composites and elongation at break for the composites
Fig. 2. The multi-scale algorithm for solving the direct problems to determine the physical and mechanical properties of the PTFE-based composites
was rigidly fixed. On the top edge, a tensile load was applied along the normal axis, while the side edges were stress-free. At the "matrix-particle" boundaries, ideal adhesion conditions were assumed. The problems were solved in a flat statement. Upon loading, the non-linear PTFE properties were taken into account. The calculated results, obtained according to the algorithm shown in Fig. 2, are given below.
2. Results and discussion
The calculated stress-strain properties. After solving the problems of determining the SSS parameters via averaging procedures [22,23], the mechanical properties of the composites were calculated at the first structural level. Tab. 2 presents dependences of their values on the filling degree for the composite loaded with filler 1 (1.5 pm).
Table 2. The mechanical properties calculated at the first structural level for the composite loaded with filler 1 (1.5 pm)
Volume fraction, <fm 10 vol.% 20 vol.% 30 vol.%
Elastic modulus, MPa 471 487 495
Poisson's ratio 0.370 0.381 0.376
Elongation at break 0.068 0.064 0.045
Rising the volume fraction expectedly led to an increase in the elastic modulus and a decrease in the average stresses over the computational domain (Fig. 4), corresponding to the crack initiation stage (marked by an elliptical contour).
Fig. 5 shows stress-strain curves obtained by stretching composites with different degrees of filling, characterizing the effective properties of PTFE loaded with filler 1 (1.5 pm).
Further, the mechanical properties obtained at the first structural level were used to determine the SSS parameters for the composite loaded with filler 2 (60 pm). In this way, they were
Fig. 3. The computational domains for determining the physical and mechanical propertiesat the first and second structural levels (dimensions are in microns)
Fig. 4. A tensile stress distribution for the composite loaded with 10 vol.% filler 1 at the crack initiation stage
designated as values for the "modified matrix". Fig. 6 shows a stress distribution in the computational domain during uniaxial tension of the composite loaded with 30 vol.% filler 1 (60 pm) at the crack initiation stage (represented by a horizontal segment inside the elliptical contour).
Tab. 3 shows dependences of elongation at break versus the volume fraction and the elastic modulus of the modified matrix. It should be noted that the data obtained at the second structural level were similar to those for the first one. In particular, rising the volume fraction was accompanied by enhancing the elastic modulus and lowering the strain level, corresponding to the beginning of fracture under the tensile load.
Calculated thermophysical properties of the composites. An approach applied previously for particulate composites [23,24] was to a certain extent similar to the method of determining the mechanical properties in the studied cases. The static thermal conductivity problem was solved in relation to the non-uniform computational domain, taking into account the specific characteristics of both matrix and filler materials. As a result, a temperature distribution over the computational domain was calculated, on the basis of which thermal conductivity was assessed according to the method proposed in [24]. The computational domains were the same as those used to determine the SSS parameters shown in Fig. 3. The lowest preset temperature level was 0 °C, while the highest one was 20 °C. The obtained thermal conductivity values are presented in Tab. 4 for the modified matrix loaded with filler 2 (60 pm) at its different volume fractions. With rising the contents of thermally conductive copper particles, the corresponding parameter enhanced for the composites as a whole. At the second structural level, the thermal conductivity
15
0,00 0,02 0,04 0,06 0,08
Strain
Fig. 5. The stress-strain curves for the composite loaded with filler 1 (1.5 pm) at its different filling degrees
Fig. 6. The stress distribution in the computational domain during uniaxial tension of the composite loaded with 30 vol.% filler 1 (60 pm) at the crack initiation stage
problem was solved in the same way.
As an example, Fig. 7 shows a temperature distribution in the composite loaded with 20 vol.% filler 2 (60pm)
Tab. 5 presents the thermal conductivity values for the composite loaded with filler 2 (60 pm), obtained at the second structural level depending on both its volume fraction and thermal conductivity of the modified matrix. It could be concluded based on these data that an increase in the volume fraction expectedly caused rising the thermal conductivity values.
Considering that the thermal conductivity value of PTFE was taken as 0.20 W/(m °C), the data presented in Tab. 5 indicated the possibility of its rising up to 0.33 W/(m °C), i.e. by more than 60%. In this case, as shown above, the limiting elongation at break level decreased, but not as significantly as after loading with the only filler 2 (60 pm). So, loading with both fillers, corresponding to two structural levels that enabled to improve the functional characteristics of the PTFE-based composites, gave the much greater positive effect. The data presented in Tabs. 3 and 5 are also shown in Figs. 8 and 9, respectively, as surfaces and ranges of the control parameters that satisfied the specified boundary conditions. Since the thermal conductivity and elastic modulus values were obtained at the first structural level for different volume fractions of
Table 3. The dependences of the elongation at break values on the Em elastic modulus of the modified matrix and the <m volume fraction for the filler 2 (60 pm)
Elastic modulus Em, MPa Volume fraction , <fm
10 vol.% 20 vol.% 30 vol.%
471 0.0220 0.0226 0.0203
487 0.0279 0.0260 0.0245
495 0.0191 0.0187 0.0176
Table 4. The thermal conductivity values for the modified matrix obtained at the first structural level
Volume fraction, <m 10 vol.% 20 vol.% 30 vol.%
Thermal conductivity, W/(m°C) 0.240 0.254 0.277
Fig. 7. The temperature distribution in the composite loaded with 20 vol.% filler 2 (60 pm)
Table 5. The dependence of thermal conductivity W/(m °C) of the composites on both fillers volume fraction and thermal conductivity of the modified matrix
Thermal conductivity of the modified matrix, W/(m °C) Volume fraction, 10% 20% 30%
0.240 0.250 0.266 0.290
0.254 0.266 0.280 0.305
0.277 0.290 0.307 0.330
filler 1, the horizontal scale corresponds to this parameter, while the axis values are normalized for drawing the continuous surfaces in Figs. 8 and 9.
Fig. 10 shows the overlapped areas that determined the control parameters and the required functional characteristics, i.e. thermal conductivity of > 0.3 W/(m °C) and elongation at break > 0.02.
Fig. 8. The temperature distribution in the composite loaded with 20 vol.% filler 2 (60 pm)
Fig. 9. The temperature distribution in the composite loaded with 20 vol.% filler 2 (60 pm)
Fig. 10. The temperature distribution in the composite loaded with 20 vol.% filler 2 (60 pm)
Conclusions
In the study, the advanced approach was implemented to increase thermal conductivity of the PTFE-based composites with the minimal decrease in their mechanical properties. Loading with copper powders characterized by different particle sizes made it possible to solve this issue. Due to the large difference in the particle sizes, it was done at several structural levels. In this way, the sophisticated algorithm was developed and applied by implementing the multi-scale model.
The work was performed according to the government research assignment for ISPMS SB RAS, project FWRW-2021-0010.
References
[1] Multifunctionality of Polymer Composites: Challenges and New Solutions, Edited by K.Friedrich and U.Breuer, Elsevier, 2015.
[2] Handbook of multiphase polymer systems, Edited by A.Boudenne, L.Ibos, Y. Candau and S.Thomas, Wiley & Sons, 2011.
[3] G.Wypych, Handbook of Fillers, ChemTec, Publishing, 2016.
[4] S.Ebnesajjad, Expanded PTFE Applications Handbook, Technology, Manufacturing and Applications, 1st Edition, Imprint: W. Andrew, 2016.
[5] S.Kolhe, A.Deshpande, K.Wangikar, Wear Behavior of Polytetrafluoroethylene Composites: A Review, In: M.Kolhe, P. Labhasetwar, H. Suryawanshi (eds), Smart Technologies for Energy, Environment and Sustainable Development, Lecture Notes on Multidisciplinary Industrial Engineering, Springer, Singapore, 2019.
[6] Teflon PTFE Fluoropolymer Resin: Properties Handbook, DuPont Fluoroproducts, 1996.
[7] J.K.Lancaster, Polymer-based bearing materials: the role of fillers and fibre reinforcement, Tribol. Int. 5(1972), 249-255.
[8] B.J.Briscoe, Wear of polymers: an assay on fundamental aspects, Tribol. Int., 14(1981), 231-243.
[9] F.Yan, Q.Xue, IR study of the interaction between graphite and PTFE during mixing process, Chin. Sci. Bull, 42(1997), 1164-1168.
10] T.Onodera, M.Park, K.Souma, N.Ozawa, M.Kubo, Transfer-Film Formation Mechanism of Polytetrafluoroethylene: A Computational Chemistry Approach, J. Phys. Chem. C, 117(2013), 10464-10472. DOI: 10.1021/jp400515j
11] Yu.Wang, F.Yan, Tribological properties of transfer films of PTFE-based composites, Wear, 261(2006), 1359-1366. DOI: 10.1016/J.WEAR.2006.03.050
12] D.Gong, Q.Xue, H.Wang, Physical models of adhesive wear of polytetrafluoroethylene and its composites, Wear, 147(1991), 9-24.
13] Y.Yang, M.Strobel, S.Kirk, M.J.Kushner, Fluorine plasma treatments of poly(propylene) films, 2 - Modeling reaction mechanisms and scaling, Plasma Processes and Polymers, 7(2010), no. 2, 123-150. DOI: 10.1002/ppap.200900114123
14] G.Liu, C.Gao, X.Li, C.Guo, Y.Chen, J.Lv. Preparation and properties of porous polyte-trafluoroethylene hollow fiber membrane through mechanical operations, Journal of Applied Polymer Science, 132(2015), no. 43, 42696-42705. DOI: 10.1002/app.42696
15] S.Feng, Z.Zhong, Y.Wang, W.Xing, E.Drioli, Progress and perspectives in PTFE membrane: preparation, modification, and applications, Journal of Membrane Science, 549(2018), 332-349.
16] S.Ebnesajjad. Fluoroplastics, Non-Melt Processible Fluoropolymers - The Definitive User's Guide and Data Book, Vol. 1, Imprint: W. Andrew, 2014.
17] H.Wang, A.Sun, X.Qi, Y.Dong, B.Fan, Wear properties of textured. lubricant films filled with graphite and polytetrafluoroethylene (PTFE) via laser surface texturing (LST), Tri-bology International, 167(2022), 107414.
[18] M.F.Ashby, Materials Selection in Mechanical Design, 5th ed.; Butterworth-Heinemann: Oxford, UK, 2016.
[19] A.A.Adamov, I.E.Keller, S.G.Ostrer, D.V.Seletkov, Evaluation of the performance of antifriction PTFE composites at a pressure over 60 MPa. I. Comparison of hardness and deformation properties under conditions of free and constrained compression, Mech Compos Mater, 58(2022), no. 5, 965-986 (in Russian).
[20] N.Yu.Anokhina, N.Yu.Matolygina, B.A.Lyukshin, P.A. Lyukshin, Computer design of a filled polymer composition with the required deformation and strength properties. Mechanics of composite materials and structures, 15(2009), no. 4, 600-609 (in Russian).
[21] S.A.Bochkareva, N.Y.Grishaeva, B.A.Lyukshin et al., Design of Composites with a Specified Set of Physicomechanical Properties Using Three Control Parameters, Phys Mesomech, 24(2021), 196-201.
[22] B.A.Lyukshin, S.V.Panin et al., Computer modeling and construction of filled compositions, Editor S. A. Zelepugin, Novosibirsk: Publishing House of SB RAS Nauka, 2015 (in Russian).
[23] S.A.Bochkareva, N.Y.Grishaeva et al., A Unified Approach to Determining the Effective Physicomechanical Characteristics of Filled Polymer Composites Based on Variational Principles, Mech Compos Mater, 54(2019), 775-788.
[24] N.Y.Grishaeva, P.A.Lyukshin et al., Modification of thermophysical characteristics of polymers by the introduction of microfillers / Matolygina, Mechanics of composite materials and structures, 22(2016), no. 3, 342-361 (in Russian).
Компьютерное конструирование полимерных композитов с использованием многоуровневых моделей
Сергей В. Панин Борис А. Люкшин
Институт физики прочности и материаловедения СО РАН Национальный исследовательский Томский политехнический университет
Томск, Российская Федерация
Наталья Ю. Гришаева
Томский государственный университет систем управления и радиоэлектроники
Томск, Российская Федерация
Светлана А. Бочкарева Илья Л. Панов
Институт физики прочности и материаловедения СО РАН Томский государственный университет систем управления и радиоэлектроники Национальный исследовательский Томский государственный университет
Томск, Российская Федерация
Аннотация. В работе решалась задача компьютерного конструирования дисперсно-наполненного полимерного композита на основе политетрафторэтилена (ПТФЭ) с целью повышения его теплопроводности (получения максимально возможного значения его коэффициента) с учетом ограничений на деформационно-прочностные свойства. Для решения задачи разработан и применен алгоритм с использованием многоуровневой модели. Показано, что за счет использования в качестве наполнителей частиц меди двух размерных диапазонов, отличающихся размерами более чем на порядок, можно получить материалы, которые характеризуются высокой теплопроводностью при сохранении приемлемых деформационно-прочностных свойств.
Ключевые слова: многоуровневое моделирование, полимерный композит, дисперсное наполнение, напряженно-деформированное состояние, теплопроводность, ПТФЭ.
