The problem of multilevel model development for calculating of an elastic property of polymer composite materials with a complex reinforcement structure at high temperatures is considered. It is assumed that thermal destruction processes take place in the matrix and fibers at high temperatures. In order to take into account the change in the elastic properties of the composite depending on the temperature and heating time, a 3-level structural model of the composite is proposed. At the lower level mono-fibers and a matrix consisting of 4 phases, the ratio between which changes when heated are considered. At this level, the analytical relations proposed earlier in the works of Yu.I. Dimitrienko. At the next level of the model, a unidirectional composite is considered, consisting of bundles of monofilaments and a matrix. To calculate elastic properties at this level, the method of asymptotic averaging is used, and a finite element algorithm for solving local problems of the theory of thermoelasticity arising in this method. At the 3rd structural level of the model, composites with complex reinforcement structures, in particular, fabric composites, are considered. The method of asymptotic averaging is also used to calculate the elastic properties of the composite at this level. For the numerical calculation of the elastic characteristics of polymer composites at high temperatures, specialized software has been developed that operates under the control of the SMCM software package created at the Scientific and Educational Center for Supercomputer Engineering Modeling and Development of Software Systems of the Bauman Moscow State Technical University. The article provides examples of the application of the developed multilevel model and software for textile composites based on an epoxy matrix and glass fibers. The values of all components of the tensor of the elastic moduli of the composite are calculated, which vary depending on the heating program of the composite. The microstress fields in the composite are obtained. A comparison is made of the fields of microstresses and effective elastic constants at normal temperatures, with similar values obtained using the ANSYS software package, which has been modified to enable the calculation of effective elastic constants in accordance with the proposed model. A very good agreement was obtained between the calculation results, both of the effective constants and of the microstresses fields, which allows us to speak of the high accuracy of the developed software.
Димитриенко Ю.И., Юрин Ю.В., Сборщиков С.В., Богданов И.О., Яхновский А.Д., Баймурзин Р.Р. Конечно-элементное моделирование упругих свойств тканевых полимерных композитов при высоких температурах. Математическое моделирование и численные методы. 2020. № 1. с. 3–27
The problem of calculating the integral characteristics of the viscoelasticity of composite materials is considered, based on information on similar characteristics of the composite components and its microstructure. An algorithm is proposed for predicting the effective relaxation and creep kernels of composites with an arbitrary reinforcement microstructure. The algorithm is based on the Fourier transform application and the inverse Fourier transform, as well as the method of asymptotic averaging for composites under steady-state polyharmonic vibrations. The algorithm uses exponential relaxation and creep kernels for the initial components of the composite. The basis of the computational procedure of the proposed algorithm is the finite element solution of local viscoelasticity problems over the composite periodicity cell. The result of the algorithm application is the determination of the exponential relaxation and creep kernels parameters for composite materials, which makes it possible to obtain a problem solution in a completely closed form. As an example, a numerical simulation of the viscoelastic-tic characteristics of unidirectionally reinforced carbon /epoxy composites has been carried out. It is shown that the developed algorithm allows one to obtain effective relaxation and creep kernels of the composite with high accuracy, without oscillations, which, as a rule, ac-company the methods of inverting Fourier transforms.
Димитриенко Ю.И., Юрин Ю.В., Сборщиков С.В., Яхновский А.Д., Баймурзин Р.Р. Моделирование эффективных ядер релаксации и ползучести вязко-упругих композитов методом асимптотического осреднения. Математическое моделирование и численные методы, 2020, № 3, с. 22–46.