On étudie l'écoulement turbulent d'un fluide sur une canopée, que l'on modélise comme un milieu poreux déformable. Ce milieu poreux est en fait composé d'un tapis de fibres susceptibles de se courber sous la charge hydrodynamique du fluide, et ainsi de créer un couplage fluide-structure à l'échelle d'une hauteur de fibre (honami). L'objectif de la thèse est de développer un modèle macroscopique de cette interaction fluide-structure, afin d'en réaliser des simulations numériques. Une approche numérique de simulation aux grandes échelles est donc mise en place pour capturer les grandes structures de l'écoulement et leur couplage avec les déformations du milieu poreux. Pour cela nous dérivons les équations régissant la grande échelle, au point de vue du fluide ainsi que de la phase solide. À cause du caractère non-local de la phase solide, une approche hybride est proposée. La phase fluide est décrite d'un point de vue Eulerien, tandis que la description de la dynamique de la phase solide nécessite une représentation Lagrangienne. L'interface entre le fluide et le milieu poreux est traitée de manière continue. Cette approche de l'interface fluide/poreux est justifiée par un développement théorique sous forme de bilan de masse et de quantité de mouvement à l'interface. Ce modèle hybride est implémenté dans un solveur écrit en C$++$, à partir d'un solveur fluide disponible dans la librairie CFD \openfoam. Un préalable nécessaire à la réalisation d'un tel modèle macroscopique est la connaissance des phénomènes de la petite échelle en vue de les modéliser. Deux axes sont explorés concernant cet aspect. Le premier consiste à étudier les effets de l'inertie sur la perte de charge en milieu poreux. Un paramètre géométrique est proposé pour caractériser la sensibilité d'une microstructure poreuse à l'inertie de l'écoulement du fluide dans ses pores. L'efficacité de ce paramètre géométrique est validée sur une diversité de microstructures et le caractère général du paramètre est démontré. Une loi asymptotique est ensuite proposée pour modéliser les effets de l'inertie sur la perte de charge, et comprendre comment celle-ci évolue en fonction de la nature de la microstructure du milieu poreux. Le deuxième axe d'étude de la petite échelle consiste à étudier l'effet de l’interaction fluide-structure à l'échelle du pore sur la perte de charge au niveau macroscopique. Comme les cas présentent de grands déplacements de la phase solide, une approche par frontières immergées est proposée. Ainsi deux méthodes numériques sont employées pour appliquer la condition de non-glissement à l'interface fluid/solide: l'une par interface diffuse, l'autre par reconstitution de l'interface. Cela permet une validation croisée des résultats et d'atteindre des temps de calcul acceptables tout en maîtrisant la précision des résultats numériques. Cette étude permet de montrer que l'interaction fluide-structure à l'échelle du pore a un effet considérable sur la perte de charge effective au niveau macroscopique. Des questions fondamentales sont ensuite abordées, telles que la taille d'un élément représentatif ou la forme des équations de transport dans un milieu poreux souple. / We study the turbulent flow of a fluid over a canopy, that we model as a deformable porous medium. This porous medium is more precisely a carpet of fibres that bend under the hydrodynamic load, hence initiating a fluid-structure coupling at the scale of a fibre's height (honami). The objective of the thesis is to develop a macroscopic model of this fluid-structure interaction in order to perform numerical simulations of this process. The volume averaging method is implemented to describe the large scales of the flow and their interaction with the deformable porous medium. An hybrid approach is followed due to the non-local nature of the solid phase; While the large scales of the flow are described within an Eulerian frame by applying the method of volume averaging, a Lagrangian approach is proposed to describe the ensemble of fibres. The interface between the free-flow and the porous medium is handle with a One-Domain- Approach, which we justify with the theoretical development of a mass- and momentum- balance at the fluid/porous interface. This hybrid model is then implemented in a parallel code written in C$++$, based on a fluid- solver available from the \openfoam CFD toolbox. Some preliminary results show the ability of this approach to simulate a honami within a reasonable computational cost. Prior to implementing a macroscopic model, insight into the small-scale is required. Two specific aspects of the small-scale are therefore studied in details; The first development deals with the inertial deviation from Darcy's law. A geometrical parameter is proposed to describe the effect of inertia on Darcy's law, depending on the shape of the microstructure of the porous medium. This topological parameter is shown to efficiently characterize inertia effects on a diversity of tested microstructures. An asymptotic filtration law is then derived from the closure problem arising from the volume averaging method, proposing a new framework to understand the relationship between the effect of inertia on the macroscopic fluid-solid force and the topology of the microstructure of the porous medium. A second research axis is then investigated. As we deal with a deformable porous medium, we study the effect of the pore-scale fluid-structure interaction on the filtration law as the flow within the pores is unsteady, inducing time-dependent fluidstresses on the solid- phase. For that purpose, we implement pore-scale numerical simulations of unsteady flows within deformable pores, focusing for this preliminary study on a model porous medium. Owing to the large displacements of the solid phase, an immersed boundary approach is implemented. Two different numerical methods are compared to apply the no-slip condition at the fluid-solid interface: a diffuse interface approach and a sharp interface approach. The objective is to find the proper method to afford acceptable computational time and a good reliability of the results. The comparison allows a cross-validation of the numerical results, as the two methods compare well for our cases. This numerical campaign shows that the pore-scale deformation has a significant impact on the pressure drop at the macroscopic scale. Some fundamental issues are then discussed, such as the size of a representative computational domain or the form of macroscopic equations to describe the momentum transport within a soft deformable porous medium.

Para la presente investigación de tesis desarrollada entre los distritos de Imperial y Pampas, provincia de Tayacaja, en el departamento de Huancavelica, se tiene como unidad de análisis a 60 muros de contención, donde se aplicó criterios muestréales de 23 muros f´c 210 kg/cm2. Como resultados alcanzados se logró elaborar concreto hidráulico de calidad usando agregado fino marginal con el cual se pudo satisfacer los parámetros de resistencia y durabilidad. Para demostrar la durabilidad se realizaron distintos ensayos; los que no cumplieron con las especificaciones de la norma EG-2000 fueron el de Equivalente de arena MTC E 114, Material que pasa el Tamiz N°200 MTC E 202 y Análisis granulométrico por tamizado ASTM C 136, por lo cual se efectuaron análisis más rigurosos como: el ensayo de Valor de azul de metileno y el ensayo de Análisis petrográfico macroscópico y microscópico en agregados para concreto ASTM C 295, donde se obtuvo como resultado que la arena no presenta elementos activos en su composición, lo que quiere decir que ante las condiciones ambientales de la zona estos elementos no se expanden ni se contraen y adicionalmente como medida de control se diseñó la mezcla de concreto hidráulico con aditivo incorporador de aire AirMix-200. Del mismo informe petrográfico se concluyó que la composición principal del agregado es la caliza, la cual es materia prima para la elaboración del cemento, que de alguna manera adicional ayuda a cumplir satisfactoriamente el parámetro de resistencia. Así mismo quedó demostrado que a través de ensayos a compresión no confinada de testigos de concreto, se alcanzó resistencias superiores a lo requerido en la norma EG-2000. Research for this thesis developed between the districts of Imperial and Pampas, Tayacaja province, in the department of Huancavelica, it has the analysis unit 60 retaining walls, where 23 sample criteria applied walls f´c 210 kg/cm2. Results obtained are as elaborate hydraulic concrete achievement of marginal quality using fine aggregate with which it was able to satisfy the parameters of strength and durability. To demonstrate the durability various tests were performed; those who do not comply with the specifications of the EG-2000 were rule the sand equivalent MTC E 114, passes what material the sieve No. 200 MTC E 202 and sieve analysis ASTM C 136, here by which they were made Analysis More rigorous like: The test value of methylene blue and testing of macroscopic and microscopic petrographic analysis in paragraph aggregate concrete ASTM C 295, where it resulted the sand without presents elements active in their composition, which means they say to the area Environmental Conditions These elements do not expand or contract and additionally as a control measure of self-design hydraulic concrete mix with entraining admixture AirMix-200 air. Same petrographic report concluded that the composition of the director of aggregate is limestone, which is a raw material for the production of cement, which in some additionally helps one perform satisfactorily resistance parameter. Also it demonstrated that one through essays of unconfined compressive concrete witness, higher strengths than required was reached in the EG-2000 standard.

La simulation des procédés de mise en forme des composites à renforts tissés de type RTM est un enjeu majeur pour les industries de pointe mettant en œuvre ce type de matériaux. Au cours de ces procédés, la préforme tissée est souvent soumise à des déformations importantes. La connaissance et la simulation du comportement mécanique de la préforme à l’échelle macroscopique et à l’échelle mésoscopique s’avère souvent nécessaire pour optimiser la phase de conception de pièces composites formées par de tels procédés. Une analyse du comportement mésoscopique des préformes tissées de composites est d’abord proposée. Une loi de comportement hyperélastique isotrope transverse est développée, permettant de décrire le comportement mécanique de chacun des modes de déformation de la mèche : élongation dans la direction des fibres, compaction et distorsion dans le plan d’isotropie de la mèche, cisaillement le long des fibres. Une méthodologie est proposée pour identifier les paramètres de cette loi de comportement à l’aide d’essais sur la mèche et sur le tissu, et une validation par comparaison avec des essais expérimentaux est présentée. Une analyse du comportement macroscopique des renforts interlocks est ensuite proposée : une loi de comportement hyperélastique orthotrope est développée et implémentée. Cette loi, extension de la loi de comportement pour la mèche, est également basée sur une description phénoménologique des modes de déformation de la préforme. Une méthode d’identification des paramètres de cette loi de comportement est mise en œuvre, utilisant des essais expérimentaux classiques dans le contexte des renforts tissés (tension uniaxiale, compression, bias extension test, flexion). Cette seconde loi de comportement est validée par comparaison avec des essais de flexion et d’emboutissage hémisphérique. / Simulating the preforming stage of RTM-like fabric-reinforced composites manufactoring processes is a major stake for industries which use such materials. During such processes, the woven preform often undergoes finite deformations. Simulation methods are then required to optimize the conception of composites parts formed by RTM. An analysis of the mechanical behaviour of woven preforms at mesoscale is first presented. A transversely isotropic hyperelastic behaviour law is developped in order to describe the mechanical behaviour of each deformation mode of the yarn : elongation in the direction of fibres, compaction and distorsion in the transverse plane and along-fibres shear. An identification method is set up for this behaviour law which allows to compute its parameters by use of simple experimental tests on the yarn and on the fabric. The behaviour law is then validated by comparizon between simulations ans experimental tests. An analysis of the mechanical behaviour of interlock woven preforms at macroscale is the presented. An orthotropic hyperelastic behaviour law is developped and implemented as an extension of the behaviour law for the yarn. A phenomenological approach is also used to describe the mechanical behaviour of each deformation mode of the preform. An identification method is set up and put into place, based on tests well known in the field of fabric reinforcements : tensile test, crushing test, bias extension test, flexure test. A hemispherical stamping simulation is set up and compared to experiment for validation purpose.

Afin de pouvoir prédire le comportement des renforts de composites 3D interlock au cours d'un procédé de mise en forme, il est nécessaire de connaitre la position des mèches dans le renfort durant la phase de préformage du procédé. Les travaux présentés ici traitent de la simulation du préformage de renforts 3D épais à l'aide d'un élément fini hexaédrique semi-discret spécifique. En utilisant le principe des travaux virtuels, on distingue le travail interne virtuel dû à la tension des mèches des autres travaux virtuels. La raideur due aux tensions de mèches, qui constitue la contribution principale de la rigidité du matériau, est prise en compte à l'aide de barres incluses dans les éléments. Les rigidités dues aux autres sollicitations, comme la compression transverse, les cisaillements ou les frottements inter-mèches, sont décrites par un matériau continu additionnel. La combinaison de ce modèle discret du premier ordre et d'un matériau continu hyperélastique anisotrope dit du second ordre, pour formuler le comportement du matériau va permettre la simulation du préformage des renforts tissés épais. Conjointement aux travaux sur la simulation, des travaux expérimentaux pour l'identification des paramètres matériau de la loi de comportement ont été définis et réalisés. Ces paramètres concernent les deux parties de la formulation du comportement. / In order to simulate 3D interlock composite reinforcement behavior during forming process, it is necessary to predict yarns positions in the fabric during the preforming stage of the process. The present work deals with thick 3D interlock fabric forming simulation using specific hexahedral semi-discrete finite elements. Using the virtual work principle, we distinguish the virtual internal work due to tensions in yarns from other internal virtual works. The stiffness relative to yarns tension which is the main part of the rigidity is described by bars within the elements. The other rigidities - like transverse compression, shears or friction between yarns - are depicted by a continuous additional material. A combination of this "first order" discrete model and a continuous orthotropic hyperelastic "second order" material formulation will enable us to simulate the interlock preforming process. Jointly to the simulation work, we also had to specify and perform experimental testing identification of material parameters. These parameters concern both parts of the model.

L'industrie aéronautique doit faire face aux nouvelles exigences environnementales, tout particulièrement concernant la réduction de la consommation des énergies fossiles. L'utilisation de matériaux composites plus léger permet de répondre en partie à cette attente. Pour limiter les coûts lors de la fabrication et du développement des composites à renforts tissés 3D, il est nécessaire d'utiliser des outils de simulation performants. En particulier, les outils existants, qui discrétisent à une échelle mésoscopique l'architecture des tissus 3D, ne tiennent pas compte de l'influence du procédé de fabrication sur la constitution de la structure textile. Si des outils numériques dédiés à la modélisation du procédé de tressage et de tricotage sont disponibles, il n'en est rien concernant le tissage. Cette étude avait donc pour but de s'intéresser plus particulièrement à la simulation du prodécé de tissage pour pouvoir obtenir une structure de tissu sèche déformée numériquement. La production de différentes architectures de tissu en verre E dans notre laboratoire nous a permis d'observer les différents éléments en contact avec le fil ou le tissu sur la machine à tisser, par le biais de l'utilisation d'une caméra rapide par exemple. Le développement d'un modèle numérique par éléments finis reproduisant le procédé de tissage a été réalisé. Une loi de comportement isotrope transverse fut utilisée pour modéliser les fils de verre. Des premières simulations numériques encourageantes pour la fabrication d'un tissu d'armure toile et d'un tissu d'armure croisé 2-2 sont présentées et comparées avec les tissus réels produits correspondants. / The aeronautical industry faces new challenges regarding the reduction of fossil fuel consumption. One way to address this issue is to use lighter composite materials. The ability to predict the geometry and the mechanical properties of the unit cell is necessary in order to develop 3D reinforcements in composite materials for these aeronautical applications. There is a difficulty to get realistic geometries for these unit cells due to the complexity of their architecture. Currently, existing tools which model 3D fabrics at a meso scale don't take into account manufacturing process influence on the shape modification of the textile structure. There is already some numerical tools that can model the braiding or knitting process, but none have been developed for weaving so far. Consequently, this study deals with the numerical simulation of the weaving process to obtain a deformed dry fabric structure. During the weaving process of E-glass fabrics, achieved in our laboratory, it has been observed that large deformations led to the modification of transverse section of meshes, or local density changes, that can modify the fabrics mechanical resistance. For this reason, a numerical tool of the weaving process, based on finite element modelling, has been developped to predict these major deformations and their influences on the final textile structure. The correlation between numerical results and fabrics produced with glass fibres has been achieved for plain weave and 2-2 twill.

Due to failure of the continuum hypothesis for higher Knudsen numbers, rarefied gases and microflows of gases are particularly difficult to model. Macroscopic transport equations compete with particle methods, such as the direct simulation Monte Carlo method (DSMC) to find accurate solutions in the rarefied gas regime. Due to growing interest in micro flow applications, such as micro fuel cells, it is important to model and understand evaporation in this flow regime. To gain a better understanding of evaporation physics, a non-steady simulation for slow evaporation in a microscopic system, based on the Navier-Stokes-Fourier equations, is conducted. The one-dimensional problem consists of a liquid and vapor layer (both pure water) with respective heights of 0.1mm and a corresponding Knudsen number of Kn=0.01, where vapor is pumped out. The simulation allows for calculation of the evaporation rate within both the transient process and in steady state. The main contribution of this work is the derivation of new evaporation boundary conditions for the R13 equations, which are macroscopic transport equations with proven applicability in the transition regime. The approach for deriving the boundary conditions is based on an entropy balance, which is integrated around the liquid-vapor interface. The new equations utilize Onsager relations, linear relations between thermodynamic fluxes and forces, with constant coefficients that need to be determined. For this, the boundary conditions are fitted to DSMC data and compared to other R13 boundary conditions from kinetic theory and Navier-Stokes-Fourier solutions for two steady-state, one-dimensional problems. Overall, the suggested fittings of the new phenomenological boundary conditions show better agreement to DSMC than the alternative kinetic theory evaporation boundary conditions for R13. Furthermore, the new evaporation boundary conditions for R13 are implemented in a code for the numerical solution of complex, two-dimensional geometries and compared to Navier-Stokes-Fourier (NSF) solutions. Different flow patterns between R13 and NSF for higher Knudsen numbers are observed which suggest continuation of this work. / Graduate

Dans cette thèse, nous présentons des résultats d’homogénéisation périodique d’un problème spectral et de l’équation d’ondes de Bloch. Il permet de modéliser les ondes à basse et haute fréquences. La partie modèle à basse fréquence est bien connu et n’est pas donc abordée dans ce travail. A contrario ; la partie à haute fréquence du modèle, qui représente des oscillations aux échelles microscopiques et macroscopiques, est un problème laissé ouvert. En particulier, les conditions aux limites de l’équation macroscopique à hautes fréquences établies dans [36] n’étaient pas connues avant le début de la thèse. Ce dernier travail apporte trois contributions principales. Les deux premières contributions, portent sur le comportement asymptotique du problème d’homogénéisation périodique du problème spectral et de l’équation des ondes en une dimension. La troisième contribution consiste en une extension du modèle du problème spectral posé dans une bande bi dimensionnelle et bornée. Le résultat d’homogénéisation comprend des effets de couche limite qui se produisent dans les conditions aux limites de l’équation macroscopique à haute fréquence. / In this dissertation, we present the periodic homogenization of a spectral problem and the waveequation with periodic rapidly varying coefficients in a bounded domain. The asymptotic behavioris addressed based on a method of Bloch wave homogenization. It allows modeling both the lowand high frequency waves. The low frequency part is well-known and it is not a new point here.In the opposite, the high frequency part of the model, which represents oscillations occurringat the microscopic and macroscopic scales, was not well understood. Especially, the boundaryconditions of the high-frequency macroscopic equation established in [36] were not known prior to thecommencement of thesis. The latter brings three main contributions. The first two contributions, areabout the asymptotic behavior of the periodic homogenization of the spectral problem and waveequation in one-dimension. The third contribution consists in an extension of the model for thespectral problem to a thin two-dimensional bounded strip Ω = (0; _) _ (0; ") _ R2. The homogenizationresult includes boundary layer effects occurring in the boundary conditions of the high-frequencymacroscopic equation.

Dans cette thèse, on considère plusieurs problèmes issus de la modélisation macroscopique des mouvements de foule. Le premier modèle consiste en une loi de conservation avec un flux discontinu, le second est un système mixte hyperbolique-elliptique et le dernier est une équation non-locale. D'abord, on utilise le modèle de Hughes une dimension pour décrire l'évacuation d'un couloir avec deux sorties. Ce modèle couple une loi de conservation avec un flux discontinu à une équation eikonale. On implémente la méthode de suivi de fronts, qui traite explicitement le comportement de la solution non-classique au point de rebroussement, afin d'obtenir des solutions de référence. Elles serviront à tester numériquement la convergence de schémas aux volumes finis classiques. Ensuite, on modélise le croisement de deux groupes marchant dans des directions opposées avec un système de lois de conservation mixte hyperbolique-elliptique dont le flux dépend des deux densités. Le système perd son hyperbolicité pour certainement valeurs de densité. On assiste à l'apparition d'oscillations persistantes mais bornées, ce qui conduit à la reformulation du problème associé dans le cadre des mesures de probabilités. Finalement, on étudie un modèle non-local de trafic piétonnier en deux dimensions. Le modèle consiste en une loi de conservation dont le flux dépend d'une convolution de la densité. Avec ce modèle, on résout un problème d'optimisation pour une évacuation d'une salle avec une méthode de descente, évaluant l'impact du calcul explicite du gradient de la fonction coût avec la méthode de l'état adjoint plutôt que son approximation par différences finies. / In this thesis, we consider nonclassical problems brought out by the macroscopic modeling of pedestrian flow. The first model consists of a conservation law with a discontinuous flux, the second is a mixed hyperbolic-elliptic system of conservation laws and the last one is a nonlocal equation. In the first chapter, we use the Hughes model in one space-dimension to represent the evacuation of a corridor with two exits. The model couples a conservation law with discontinuous flux to an eikonal equation. We implement the wave front tracking scheme, treating explicitly the solution nonclassical behavior at the turning point, to provide a reference solution, which is used to numerically test the convergence of classical finite volume schemes. In the second chapter, we model the crossing of two groups of pedestrians walking in opposite directions with a system of conservation laws whose flux depends on the two densities. This system loses its hyperbolicity for certain density values. We assist to the rising of persistent but bounded oscillations, that lead us to the recast of the problem in the framework of measure-valued solutions. Finally we study a nonlocal model of pedestrian flow in two space-dimensions. The model consists of a conservation law whose flux depends on a convolution of the density. With this model, we solve an optimization problem for a room evacuation with a descent method, evaluating the impact of the explicit computation of the cost function gradient with the adjoint state method rather than approximating it with finite differences.

Dans le cadre du transport de marchandises par voie ferroviaire, un certain nombre de verrous limitent les efforts consentis pour un fret plus "propre". Des améliorations doivent être faites afin de limiter le nombre de locotracteurs assurant les charges/décharges de marchandises en bout de ligne. Dans cette perspective, FEMTO-ST et Alstom Transport ont pour objectif de concevoir et développer un système de gestion d'énergie pour locomotive électrique hybride. Composée d'un groupe électrogène couplé à des accumulateurs électrochimiques et des super-condensateurs, la locomotive électrique hybride permet d'accroître la souplesse et l'efficacité du transport ferroviaire électrique, tout en réduisant encore son impact environnemental. Cette étude a consisté dans un premier temps à développer une modélisation macroscopique de la chaîne de traction électrique hybride et à proposer une structuration de la commande avec l'identification des capteurs matériels et logiciels nécessaires au contrôle optimal de la chaîne de traction. La caractérisation expérimentale des moyens de stockage a permis l'amélioration comportementale et dynamique des modèles correspondant. Dans un second temps, le contrôle de la tension de sortie d'un hacheur dévolteur a permis d'étudier l'application de la logique floue de type-2 (intervalle et générale) dans le cas d'applications industrielles (relativement simples). Enfin, une gestion innovante des flux d'énergie au sein d'un système plus complexe: locomotive électrique hybride, a été développée. Mettant en œuvre la logique floue de type-2 et des algorithmes s'inspirant de la théorie de l'évolution, cette gestion d'énergie optimale opère en temps réel et sans connaissance à priori du cycle de conduite. La gestion fréquentielle ainsi que le contrôle des états de charge des sources secondaires du véhicule contribuent également à leur bon fonctionnement. / To achieve a "greener" freight transport, efforts are still needed to overcome some technological barriers. New improvements must be carried at to limit the shunting locomotives' use intended for the goods' load/unload. Considering this aim, FEMTO-ST and Alstom Transport decided to conceive and develop an energy management strategy system for hybrid electric locomotives. Such locomotives include a diesel driven generator set which is coupled with batteries and ultra-capacitors. The architecture aims at improving the flexibility, the effectiveness of the electrical railway transport and at reducing the environmental impact of these activities again. Firstly, the study consists in implementing a macroscopic modelling of the hybrid electric powertrain. Then, the control is optimally designed by identifying the hardware and software sensors of the powertrain. The dynamics and the behavior of the secondary sources' models are improved thanks to their experimental characterizations. Secondly, the use of Type-2 Fuzzy Logic (interval and general) controllers permits to study their efficiency in the control of a very simple industrial system: the output voltage of a buck converter. Lastly, based on the obtained results, an innovative management of the system's energy flows is developed in the case of the hybrid electrical locomotive. The use of the type-2 fuzzy logic and evolutionary algorithms permit to optimally perform a real time energy management strategy without a priori knowledge of the duty cycle. The frequency approach and the secondary sources' state of charge control contribute to their efficient use.

The diploma thesis deals with the observation and diagnostics of transformations that take place in stone slabs of different mineralogical composition during heating to high temperatures. The theoretical part of the thesis deals with a summary of available information about selected types of rocks. Furthermore, it presents findings from the experimental thermal loading of a granite sample. The principles and instruments of laboratory experimental methods that are used to judge the resistance of rocks are described in detail. In the experimental part, the testing methodology is proposed, the procedure of production of rock test specimens is described and the way in which the thermal loading of stone slabs took place is explained. The results of the experiments are shown in the tables and graphs. A comparison of individual tests and properties of selected samples and evaluation of the obtained results is performed.