The Usage of Genetic Methods for Prediction of Fabric Porosity

Advanced fabric production demands developing strategies with regard to new fabric constructions in which sample-production is reduced to a minimum. It is clear that a new fabric construction should have the desired end-usage properties pre-specified as project demands. Achieving such a demand is a complex task based on our knowledge of the relations between the fabric constructional parameters and the predetermined fabric endusage properties that fit the desired quality. Individual fabric properties are difficult to predict when confronting the various construction parameters, which can be separated into the following categories: raw materials, fabric structure, design, and manufacturing parameters.


Introduction
Advanced fabric production demands developing strategies with regard to new fabric constructions in which sample-production is reduced to a minimum. It is clear that a new fabric construction should have the desired end-usage properties pre-specified as project demands. Achieving such a demand is a complex task based on our knowledge of the relations between the fabric constructional parameters and the predetermined fabric endusage properties that fit the desired quality. Individual fabric properties are difficult to predict when confronting the various construction parameters, which can be separated into the following categories: raw materials, fabric structure, design, and manufacturing parameters.
Many attempts have been made to develop predictive models for fabric properties with different modelling tools. There are essentially two types of modelling tools: deterministic (mathematical models, empirical models, computer simulation models) and nondeterministic (models based on genetic methods, neural network models, models based on chaos theory and theory of soft logic), and each of them has its advantages and disadvantages [1].
Deterministic modelling tools present the heart of conventional science and have their basis in first principles, statistical techniques or computer simulations. Mathematical models offer a deep understanding of relations between constructional parameters and predetermined fabric property, but due some simplifying assumptions large prediction errors occur. Empirical models based on statistical techniques show a much better agreement with the real values but the problems with samples preparing, process repeatability, measurements errors and extrapolation occur. They usually refer to the one type of testing method of particular fabric property. The advantage of computer simulation models is their ability to capture the randomness inherent in fabric structure so the predicted values are very near the v V V ε = (1) where, ε is the porosity expressed as coefficient, Vv is the volume of the total void space in cm 3 , and V is the total or bulk body volume in cm 3 . The total volume of the body consists of the volumes of the solid and void components as follows: where, V is the total volume of the body in cm 3 , Vv is the volume of void component in cm 3 , and VS is the volume of solid component in cm 3 . If the volume of void component is exposed from the Equation 2, the Equation 1 can be further written as follows: where, ε is the porosity expressed as coefficient, Vs is the volume of solid component in cm 3 , V is the volume of the body (or bulk volume) in cm 3 , ms is the mass of solid component in g, mb is the mass of the body (or bulk mass) in g, ρb is the bulk density in g/cm 3 , and ρs is the density of solid component in g/cm 3 .
In this way exactly defined porosity of the material is useful parameter, only, when materials with the same porous structure are compared, and gives an indication which material possesses more void space in the bulk volume. It does not give any information about the porous structure of the material, so it is an insufficient parameter for describing fibre assembly characteristics [16]. Namely, the materials with the same porosity could have very different porous structure and consequently, in the case of fabrics, different protection, filtration, sound absorption, etc., properties; so the need to define porous structure and some other porosity parameters is essential. From the theoretical point of view, the porosity parameters could be easily determined on the basis of an ideal geometrical model of the material porous structure. The simpler models consider that all pores, whatever their shape, are the same and regularly arranged in a fibre assembly [16,17]. Ideal models are based also on some other simplifying assumptions depending on the fibre assembly type. Porosity parameters calculated on the basis of ideal models of porous structures are usually not in a good correlation with the real porosity parameters. Real porous media generally have rather complex structures that are relatively difficult to define. But the advantage of ideal geometric models of porous structures is the possibility to understand the influence of porous structure on some end-usage properties of the material, which is crucial by a new product planning.
The fundamental building elements of the material porous structure are pores (also capillaries, channels, holes, free volume) [15,18]. Pores are void spaces within the material which are separated between each other, and are classified [19,20]: 1. according to the position in the material into: a. inter-pores, e.g. pores which lie between the structural elements of the material, b. intra-pores, e.g. pores which lie within the structural element of the material; 2. according to the pore width (the shortest pore diameter) into: a. macropores whose pore-width is greater than 50 nm, b. mesopores whose pore-width lies in the range between 2 and 50 nm, and c. micropores with the pore-width lower than 2 nm; 3. according to the fluid accessibility into ( Figure 1 4. according to the pore shape into ( Figure 2): a. cylindrical pores, b. slit-shape pores, c. cone-shape pores, and d. ink bottle pores; Figure 2. Types of pores according to the pore shape [19] 5. according to the geometry of pore-cross section into ( Figure 3): a. pores with geometrically regular cross-sectional shape and b. pores with geometrically irregular cross-sectional shape, Figure 3. Different shapes of pore cross-sections [20] 6. according to the uniformity of pore cross-section over the pore length into ( Figure 4): a. pores with a permanent cross-section, b. pores with a different cross-sections and for which different diameters are defined (the most constricted, the largest, the mean pore diameters). Four groups of pore descriptors, e.g. size, shape, orientation, and placement, are defined as important parameters [21]. Pores can be mathematically assessed on the basis of known model of pores geometry and constructional parameters of the material with the following parameters: the number of pores, pore size, pore volume, pore surface area, pore length, etc.
On the basis of an ideal geometrical model of porous structure, the pore size distribution which is also an important parameter of material porous structure can not assessed while the pores in geometrical model are usually assumed to be the same sizes. Such situation rarely occurs in the real fabrics. The further considerations of ideal geometrical models of material porous structures and porosity parameters will be focused on different types of fabrics.
Fabrics are flat textile materials which are produced by different manufacturing techniques using different fibrous forms of input material (or structural element), and consequently having different porous structures. Following basic types of fabrics are known ( Figure 5): • woven fabrics which are made by interlacing vertical warp and horizontal weft yarns at right angles to each other, • knitted fabrics which are made by forming the yarn into loops and their interlacing in vertical (warp-knitted fabrics) or horizontal (weft-knitted fabrics) direction, • nonwoven fabrics which are produced from the staple fibres, filaments or yarns by different web-forming, bonding and finishing techniques. While this chapter is focused on the genetic methods in order to predict porosity of woven and nonwoven fabrics, only those types of fabrics and their ideal geometric models of porous structure will be presented.

Woven fabric's ideal geometric model of porous structure
When a woven fabric is treated as a three dimensional formation, different types of pores are detected [22,23,24]: 1. inter-pores, e.g. the pores which are situated between warp and weft yarns (macropores, interyarn pores) and pores which are situated between fibres in the yarns (mesopores, interfiber/intrayarn pores), 2. intra-pores, e.g. the pores which are situated in the fibres (micropores, intrafiber pores). The structure and dimensions of the inter-or intrayarn pores are strongly affected by the yarn structure and the density of yarns in the woven structure [22]. As fibrous materials, woven fabrics have, with regard to knitted fabrics or nonwovens, the most exactly determined an ideal geometrical model of a macroporous structure in the form of a tube-like system, where each macropore has a cylindrical shape with a permanent cross-section over all its length ( Figure 6) [25]. Because the warp density is usually greater than the weft density, the elliptical shape of the pore cross-section is used to represent the situation in Figure 6. Macropores are opened to the external surface and have the same cross-section area. They are separated by warp or weft yarns, and are uniformly distributed over the woven fabric area.
The primary constructional parameters of woven fabrics which alter the porous structure are: • yarn fineness, e.g. the mass of 1000 meter of yarn from which the yarn diameter can be calculated, • type of weave, e.g. the manner how the yarns are interlaced. It has an effect on the pore size as well as on the shape of pore cross-section [26], • the number of yarns in length unit (warp and weft densities), which directly alters the pore size.
When fibre properties (fibre density, dimension, and shape) are different, two woven fabrics with similar woven structures and geometrical configurations can have distinctly different porosity [22]. To compare woven fabrics with porosity, the following porosity parameters can be calculated on the basis of the woven fabric primary constructional parameters and the ideal model of porous structure in the form of a tube-like system: • (total) porosity by using Equation 6 where the bulk density of the material is actually the woven fabric density and the density of solid component is the yarn density. If the fibre volume fraction (yarn packing factor) is exposed from the Equation 7 which represents the yarn diameter calculation, and then inserted in Equation 8 by assuming Equation 9 for woven fabric density at the same time, the porosity of woven fabrics can be then written in the form of where, d is the yarn diameter in cm, T is the yarn fineness in tex, ρyarn is the yarn bulk density in g/cm 3 , ρfib is the fibre density in g/cm 3 , βfib is the fibre volume fraction (or yarn packing factor), ρfab is the woven fabric bulk density in g/cm 3 , m is the woven fabric mass per unit area in g/m 2 , D is the woven fabric thickness in mm, ρb is the body bulk density in g/cm 3 , and ρs is the density of solid component in g/cm 3 . It is worth to mention that in this way defined porosity refers to all types of pores regarding their position in the woven fabric, e.g. inter-and intra-pores; • area of pore cross-section which refers only on macropores in a woven fabric. The ideal model of woven fabric porous structure is based on the assumption that macropores have cylindrical shape with circular cross-section. In real woven fabrics, the macropore cross-section shape is more likely to be irregular rather regular (Figure 7) [26]. The shape of pore cross-section and consequently the area of pore cross-section depend on the type of yarns used. Woven fabrics made from filament yarns have pure macropores with rectangular cross-sections, whilst woven fabrics made from staple yarns have a small percentage of pure macropores, some of partly latticed macropores as well as fully latticed macropores (as the consequence of the phenomenon of latticed pores) with mostly irregular cross-sections. The area of pore cross-section also depends on the phenomenon of changing the position of warp threads according to the longitudinal fabric axis and the phenomenon of thread spacing irregularity [28]. For the theoretical calculations of the macropore cross-section area three types of regular pore crosssection shapes are taken into account, e.g. circular (Equation 11), rectangular (Equation 12) and elliptical (Equation 13) as follows: where, Ap is the area of macropore cross-section in mm 2 , p is the yarn spacing in mm, d is the yarn diameter in mm, g is the number of yarns per unit length in threads/cm, and subscripts 1 and 2 indicate warp and weft yarns, respectively;  • number of macropores in the area unit (pore density). It can be seen from Figure 7, that one macropore belongs to one warp yarn and one weft yarn, so the number of macropores can be calculated on the basis of warp and weft densities using Equation 14: where, Np is the pore density in pores/cm 2 , g1 is the warp density in threads/cm, and g2 is the weft density in threads/cm; • open porosity (open area) which describes the fraction of macropore cross-section area in the area unit of woven fabric. If we assume elliptical macropore cross-section area (Figure 7), the open porosity is calculated as follows: where, εopen is the open porosity, Ap is the macropore cross-section area in mm 2 , Ay is the projection area of warp and weft yarns, which refers to one macropore in mm 2 , p is the yarn spacing in mm, d is the yarn diameter in mm, and subscripts 1 and 2 indicate warp and weft yarns, respectively. Open porosity can be calculated also on the basis of cover factor (Equation 16) where, εopen is the open porosity, K is the woven fabric cover factor, d is the yarn diameter in mm, g is the warp/weft density in threads/cm, Np is the pore density in pores/cm 2 , Ap is the area of macropore cross-section in cm 2 , and subscripts 1 and 2 indicate warp and weft yarns, respectively; • equivalent macropore-diameter. If we assume that macropore has cylindrical shape, then the area of macropore cross-section is equal to the area of circle with radius r (Equation 18). Equivalent macropore diameter is the diameter of macropore with circular cross-section whose area is the same as the area of the macropore with irregular cross-section shape (Equation 19) [30].
where, Acircle is the circular cross-section macropore area in mm 2 , r is the macropore radius in mm, d is the macropore diameter in mm, de is the equivalent macropore diameter in mm, and Ap is the macropore cross-section area of macropore with irregular shape in mm 2 ; • maximal an minimal macropore diameters which refer to the elliptical shape of macropore cross-section. In the case where warp density is greater than weft density the maximal diameter is equal to p2-d2, while minimal diameter is equal to p1-d1 ( Figure 7); • macroporosity which describes the portion of macropore volume in volume unit of woven fabric. In general, it is defined using Equation 20. In the case of the elliptical macropore cross-section shape, the macroporosity, defined with Equation 21, is the same as open porosity: 12 12 where, εmacro is the macroporosity, Vp is the macropore volume in cm 3 , Vy is the volume of warp and weft yarns which refers to one macropore in cm 3 , p is the yarn spacing in mm, d is the yarn diameter in mm, D is the woven fabric thickness in mm, Ap is the macropore area in mm 2 , εopen is the open porosity, and subscripts 1 and 2 indicate warp and weft yarns, respectively.

Nonwoven fabric's ideal geometric model of porous structure
The porous structure of nonwoven fabric is a result of nonwoven construction (the type and properties of fibres or yarns as input materials, fabric mass, fabric thickness, etc.) as well as technological phases, e.g. the type of web production, bonding methods and finishing treatments. According to several different methods to produce non-woven fabrics having consequently very different porous structure, the ideal geometric model of porous structure in the form of tube-like system is partially acceptable only by those nonwovens which are thin and translucence, e.g. light polymer-laid nonwovens and some thin spun-laced or heatbonded nonwovens (Figure 8). Such model is based on the assumptions that fibres having the same diameter are distributed only in the direction of fabric plane and the distance between fibres and the length of individual fibres is much greater than the fibre diameter. Xu [21] found out that in most nonwoven fabrics, pore shape is approximately polygonal and that pores appear more circular when the fabric density increases. Pore orientation to some extent relates to fibre orientation. If pores are elongated and predominantly oriented in one direction, fibres are likely to be oriented in that direction. The variation in pore size is inherently high. Some regions may contain more pores than others or may have larger pores than those in other regions. The primary constructional parameters of nonwoven fabrics which alter the porous structure are: • fibre fineness, e.g. the mass of 1000 meter of fibre, from which the fibre diameter can be calculated, • web mass per unit area and • web thicknesses.
To compare nonwoven fabrics with porosity, the following porosity parameters can be calculated on the basis of the nonwoven fabric primary constructional parameters and the ideal model of porous structure in the form of a tube-like system: • (total) porosity by using Equation 6 where the bulk density of the material is actually the nonwoven fabric density and the density of solid component is the fibre density. The nonwoven fabric density is calculated on the basis of primary nonwoven constructional parameters, e.g. fabric mass and thickness using Equation 9 where index fab in this case refers to the nonwoven fabric. Substituting Equation 9 into Equation 6, final Equation 22 of nonwoven porosity which refers to inter-(pores between fibres in nonwovens) and intra-pores (pores inside the fibres) is obtained: where, ε is the nonwoven fabric porosity, ρb is the body bulk density in g/cm 3 , ρs is the density of solid component in g/cm 3 , ρfab is the nonwoven fabric density in g/cm 3 , ρfib is the fibre density in g/cm 3 , mfab is the nonwoven fabric mass per unit area in g/m 2 , and Dfab is the nonwoven fabric thickness in mm; • opening diameter which is the diameter of the maximum circle that can fit in a pore ( Figure 8). It is predicted on the basis of nonwoven fabric constructional parameters and refers to the heat-bonded nonwoven fabrics, as follows [17,21]: where, d0 is the opening diameter in µm, C is the thickness factor, L is the specific total length of fibres per nonwoven unit area in mm -1 , dfib is the fibre diameter in µm, Dfab is the nonwoven thickness in mm, mfab is the nonwoven fabric mass per unit area in g/m 2 , and ρfib is the fibre density in g/cm 3 ; • average area of pore cross-section which is for un-needled fabrics (e.g. fabrics made of layers of randomly distributed fibres) predicted using Equation 26 [17]: where, Ap is the average area of pore cross-section in mm 2 , ε is the porosity, and dfib is the fibre diameter in µm. On the basis of calculated average area of pore-cross-section, the equivalent pore diameter is then calculated using Equation 19. Needle-punched nonwoven fabric is a sheet of fibres made by mechanical entanglement, penetrating barbed needles into a fibrous mat [31]. Needle-punched nonwovens represent the largest segment of filtration materials used as dust filters [32]. The geometrical model of three-dimensional needle-punched nonwoven fabric proposed by Mao & Rusell [33], is also known from the literature, and it is constructed on a two-dimensional fibre orientation within the fabric plane, with interconnecting fibres oriented in the z-direction ( Figure 10). Such model relies on the following basic assumptions: 1. the fibres in the fabric have the same diameter, and a fraction of the fibres is distributed horizontally in the two-dimensional plane, the rest are aligned in the direction of the fabric thickness, 2. fibre distribution in both the fabric plane and the z-direction is homogeneous and uniform, 3. in each twodimensional plane, the number of fibres oriented in each direction is not the same, but obeys the function of the fibre orientation distribution Ω(α), where α is the fibre orientation angle, 4. the distance between fibres and the length of individual fibres is much greater than the fibre diameter. The basic porosity parameters which are based on the mentioned geometrical model of needle-punched nonwoven fabric are still difficult to define due to the fact that in each fabric planes fibres lie in different direction and in this way produce pores with different orientations, diameters, connectivity and accessibility to fluid flow ( Figure 9). The only porosity parameters that are calculated from such model are: Figure 9. Geometrical models of needle-punched nonwoven fabric and porous structure [14,34] • total porosity (Equation 22) and • mean pore diameter which is deduced from the fibre radius and porosity according to the following relation proposed by White [34]: where, dp is the mean pore diameter in µm, ε is the nonwoven fabric porosity, dfib is the fibre diameter in µm, T is the fibre linear density in tex, and ρfib is the fibre density in g/cm 3 .
Three kinds of pores may be present in needle-punched nonwoven fabrics, namely, closed pores, open pores, and blind pores. The important pore structure characteristics of needlepunched nonwoven fabrics as filter media are the most constricted open pore diameter (smallest detected pore diameter), the largest pore diameter (bubble point pore diameter), and mean pore diameter (mean flow pore diameter) [35].

The usage of genetic programming to predict woven fabric porosity parameters
Porosity parameters based on an ideal geometrical model of porous structure give woven fabric constructor some useful information about porosity by developing a new product, but they are not in a good agreement with the experimental values. In order to balance the difference between the theoretical and experimental values of porosity parameters, genetic programming was used to develop models for predicting the following macro-porosity parameters of woven fabric: the area of macro-pore cross-section, macro-pore density, open porosity, and equivalent macro-pore diameter. The genetic programming is a variant of evolutionary algorithm methods described in many sources (e.g., [2,3,4]). The basic information on the evolutionary algorithms is given at the beginning of the section 4. We implemented Koza's variant of genetic programming [2]. In our research, the independent input variables (the set of terminals) were: yarn fineness T (tex), weave value V, fabric tightness t (%) and denting D (ends/dent in the reed). The set of terminals also included random floating-point numbers between -10 and +10. Variegated reed denting was treated as an average value of treads, dented in the individual reed dent. The dependent output variables were: area of macro-pore cross-section Ap (10 -3 mm 2 ), pore density Np (pores/cm 2 ), and equivalent macro-pore diameter (µm). For all modelling, the initially set of functions included the basic mathematical operations of addition, subtraction, multiplication, and division. In the case of modelling the area of macro-pore cross-section and pore density the initially set of functions also included a power function, whereas the set of functions for modelling of equivalent macro-pore diameter included an exponential function. We then used the genetic programming system to evolve appropriate models consist of abovementioned sets of terminals and functions. Open porosity was calculated on the basis of predicted values of the area of macro-pore cross-section and macro-pore density and Equation 17. The equivalent macro-pore diameter was calculated on the basis of predicted values of the area of macro-pore cross-section using Equation 19. The fitness measure for modelling by genetic programming was exactly the same as defined by Equation 33 in section 4. The goal of the modelling was to find such a predictive model in a symbolic form, that Equation 33 would give as low an absolute deviation as possible.
The evolutionary parameters for modelling by genetic programming were: population size 2000, maximum number of generations to be run 400, probability of reproduction 0.1, probability of crossover 0.8, probability of mutation 0.1, minimum depth for initial random organisms 2, maximum depth for initial random organisms 6, maximum depth of mutation fragment 6, and maximum permissible depth of organisms after crossover 17. The generative method for the initial random population was ramped half-and-half. The method of selection was tournament selection with a group size of 7. For the purpose of this research 100 independent genetic programming runs were executed. Only the results of the best runs (i.e., the models with the smallest error between the measurements and predictions) are presented in the paper.

Materials and porosity measurements
Our experiments involved woven fabrics made from staple yarns with two restrictions: first, only fabrics made from 100% cotton yarns (made by a combing and carding procedure on a ring spinning machine) were used in this research; second, fabrics were measured in the grey state to eliminate the influence of finishing processes. We believe that it is very hard, perhaps even impossible, to include all woven fabrics types to predict individual macroporosity parameters precisely enough, and so we focused our research on unfinished staple yarn cotton fabrics. We would like to show that genetic programming can be used to establish the many relations between woven fabric constructional parameters and particular fabric properties, and that the results are more useful for fabric engineering than ideal theoretical models. The cotton fabrics varied according to yarn fineness (14 tex, 25 tex, and 36 tex), weave type (weave value), fabric tightness (55% -65%, 65% -75%, 75% -85%), and denting. The constructional parameters of woven fabric samples are collected in Table 1. They were woven on a Picanol weaving machine under the same technological conditions. The weave values of plain (0.904), twill (1.188), and satin (1.379) fabrics, as well as fabric tightness, were determined according to Kienbaum's setting theory [36].
We used an optical method to measure porosity parameters of woven fabrics, since it is the most accurate technique for macro-pores with diameters of more than 10 µm. For each fabric specimen, we observed between 50 and 100 macro-pores using a Nikon SMZ-2T computeraided stereomicroscope with special software. We measured the following macro-porosity parameters: area of macro-pore cross-section, pore density, and equivalent macro-pore diameters.

Predictive models of woven fabric porosity parameters
Equations 29 and 30 present predictive models of the area of macro-pore cross-section Ap and macro-pore density Np, respectively [37]. Here V is the weave factor, T is the yarn linear density in tex, t is the fabric tightness in %, and D is the denting in ends per reed dent. The open porosity and equivalent diameter are calculated using Equations 17 and 19, respectively, where for Ap and Np the predicted values are taken into account. Because the model of the area of macro-pore cross-section is more complex, the functions f1, f2,…f10 are not presented here but are written in the appendix. When calculating the values of models, the following rules have to be taken into account: the protected division function returns to 1 if denominator is 0; otherwise, it returns to the normal quotient. The protected power function raises the absolute value of the first argument to the power specified by its second argument.
By a comparison of both GA models (Equations 29 and 30) with the theoretical ones (Equations 11-13 and 14), the complexity of GA models is obvious and derives from the factors involved in the models. Namely, factors involved in GA models don't ignore the irregularity of macro-pores cross-section area as well as the number of pores, due to the phenomenon of latticed pores in the case of staple yarns (which depends on the type of weave -factor V and fabric tightness -factor t) and the phenomenon of thread spacing irregularity (factor D), as theoretical models do. Theoretical model for the macro-pore crosssection area assumes that all macro-pores in woven structure have the same cross-section area regardless the type of used yarns, type of weave, fabric tightness and denting, whilst the theoretical model for the pore density assumes no reduction of the number of pores.  Theoretical values of woven fabric porosity parameters deviate from experimental ones on average by 118.3% (min 8.8%, max 452.9%) for the area of the macro-pore with rectangular cross-section, 111.5% (min 14.5%, max 370.6%) for the area of the macro-pore with circular cross-section, 72.8% (min 0.2%, max 335.3%) for the area of the macro-pore with elliptical cross-section, 37.3% (min 0.0%, max 395.0%) for the macro-pore density, 232.6% (min 19.9%, max 1900.1%) for the open porosity of fabrics with rectangular pore cross-section, 221.0% (min 14.3%, max 1558.0%) for the open porosity of fabrics with circular pore cross-section, 166.3% (min 5.9%, max 1479.0%) for the open porosity of fabrics with elliptical pore cross-section, 43.7% (min 4.3%, max 135.1%) for the equivalent pore diameter where rectangular pore crosssection is taken into account, 43.7% (min 7.0%, max 116.9%) for the equivalent pore diameter where circular cross-section is taken into account, and 28.0% (min 0.1%, max 108.6%) for the equivalent pore diameter where elliptical pore cross-section is taken into account.
The results of woven fabric porosity parameters determined with models based on genetic programming show very good agreement with experimental values (Figure 11) and justify the complexity of GA models. The predicted values deviate from experimental ones on average by 1.5% (min 0.0%, max 10.2%) for the area of the macro-pore cross-section, 2.0% (min 0.0%, max 8.0%) for the macro-pore density, 3.2% (min 0.0%, max 10

The usage of genetic algorithm to predict nonwoven fabric porosity parameters
In this research, the genetic algorithm was used for definition of predictive models of nonwoven fabric porosity parameters. Since needle-punched nonwoven fabrics have completely different porous structure when compared to woven fabrics, it is inappropriate to focus on open porosity through the prediction of the area of macro-pore cross-section and macro-pore density. The most valuable porosity parameters for needle-punched nonwoven porous structure characterisations are total porosity and mean pore diameter, and those parameters were the subjects of our research. Since the basic steps in evolutionary computation are well-known, only a brief description follows. Firstly, the initial population P(t) of the random organisms (solutions) is generated. The variable t represents the generation time. The next step is the evaluation of population P(t) according to the fitness measure. Altering the population P(t) by genetic operations follows. The genetic operations alter one or more parental organism(s); thus, creating their offspring. The evaluation and alteration of population takes place until the termination criterion has been fulfilled. This can be the specified maximum number of generations or a sufficient quality of solutions [38]. More comprehensive information on evolutionary computation can be found in [39].
The independent input variables were fibre fineness -T (dtex), nonwoven fabric area mass -m (g/m 2 ), and nonwoven fabric thickness -D (mm). The dependent output variables were mean pore diameter dp (µm) and total porosity ε (%). Since the GA approach is unsuitable for the evolution of prediction models (organisms) in their symbolic forms, it is necessary to define them in advance [38]. In this study, a quadratic polynominal equation with three variables was used as a prespecified model for the prediction of porosity parameters as follows:

Y c cm cD cT cm cD cT cm D cm T c D T c m D T
where, Y is the dependent output variable, m is the nonwoven fabric mass per unit area in g/m 2 , D is the nonwoven fabric thickness in mm, T is the fibre fineness in dtex, and c1…11 are constants. The main reasons for this selection were as follows: 1. a polynominal model is relatively simple, 2. for the problem studied we did not expect harmonic dependence of the output variables, 3. some preliminary modelling-runs with different types of prespecified models showed that the quadratic polynominal model provides very good selection in terms of prediction quality. In our research, the initial random population P(t) consisted of N prespecified models ( The absolute deviation D(i,t) of individual model i (organism) in generation time t was introduced as a fitness measure. It was defined as: where, E(j) is the experimental value for measurement j, P(i, j) is the predicted value returned by the individual model i for measurement j, and n is the maximum number of measurements. The goal of the optimisation task was to find such a predictive model (defined by Equation 31), that Equation 33 would give as low an absolute deviation as possible. Therefore, the aim was to find out appropriate real-valued constants in Equation 32. However, since it was unnecessary that the smallest values of the above equation also meant the smallest percentage deviation of this model, the average absolute percentage deviation of all measurements for individual model i was defined as: The Equation 33 was not used as a fitness measure for evaluating population, but only for finding the best organism within the population, after completing the run.
The altering of population P(t) was effected by reproduction, crossover, and mutation. For the crossover operation, two parental vectors, e.g., c1 and c2 were randomly selected. Then the crossover took place between two randomly-selected parental genes having the same index. Two offspring genes were created according to the extended intermediate crossover, as considered by Mühlenbeim and Schlierkamp-Voosen [40]. During the mutation operation, one parental vector c was randomly selected. Then, the mutation took place in one randomly selected parental gene. During both the crossover and mutation processes, the numbers of crossover and mutational operations performed on parental vector(s), were randomly selected. The evolutionary parameters for modelling by genetic algorithms were: population size 300, maximum number of generations to be run 5000, probability of reproduction 0.1, probability of crossover 0.8 and probability of mutation 0.1. Tournament selection with a group size of 5 was used. For the purpose of the research 200 independent genetic algorithms runs were carried out. Only the best models are presented in the paper.

Materials and porosity measurements
Bearing in mind the fact that nonwovens have very different structures and, thus, also porosity parameters due to their sequences when web-forming, bonding, as well as finishing methods, the nonwoven fabric samples were limited to one type of nonwoven fabricsthose needle-punched nonwoven fabrics made from a mixture of polyester and viscose staple fibres. Nonwoven multi-layered webs were first obtained from the same manufacturing process by subjecting the fibre mixtures to carding and then orienting the carded webs in a cross-direction by using a cross lapper to achieve web surface mass ranges of 100-150, 150-200, 250-300, and 300-350 g/m 2 , and a web volume mass range of 0.019-0.035 g/cm 3  The porosity parameters of the nonwoven fabric samples were measured using the Pascal 140 computer aided mercury intrusion porosimeter, which measures pores' diameters between 3.8 -120 µm, and operates under low pressure. The mercury intrusion technique is based on the principle that non-wetting liquid (mercury) coming in contact with a solid porous material can not be spontaneously absorbed by the pores of the solid itself because of the surface tension, but can be forced by applying external pressure. The required pressure depends on the pore-size and this relationship is commonly known as the Washburn equation [9]: where, P is the applied pressure, ϒ is the surface tension of mercury, θ is the contact-angle and r is the capillary radius. The distribution of pore size, as well as the total porosity and the specific pore volume can be obtained from the relationship between the pressure necessary for penetration (the pore dimension) and the volume of the penetrated mercury (pore volume). There are certain main assumptions necessary when applying the Washburn equation: the pores are assumed to be of cylindrical shape and the sample is pressure stable.
Each nonwoven sample of known weight was placed in the dilatometer, then the air around the sample was evacuated and finally the dilatometer was filled with mercury by increasing the pressure up to the reference level. The volume and pressure measurements' data were transferred into the computer programme and the following data were detectable or calculated: the specific pore volume (mm 3 /g), the average pore diameter (µm) and the total porosity (%). The volume of penetrated mercury is directly the measure of the sample's pore volume expressed as a specific pore volume in mm 3 /g, and is obtained by means of a capacitive reading system. The average pore diameter is evaluated at 50% of the cumulative volume of mercury.

Predictive models of nonwoven fabric porosity parameters
Equations 36 and 37 present predictive models of the total porosity ε and mean pore diameter dp, respectively. Here T is the fibre fineness in dtex, m is the nonwoven fabric mass per unit area in g/m 2 , and D is the nonwoven fabric thickness in mm. 22 Figure 12 presents a comparison of the experimental, predicted and theoretical values of porosity parameters, e.g. total porosity and mean pore diameter. The theoretical values of total porosity and mean pore diameters were calculated using Equation 22 and 27-28, respectively. In Figure 12, the theoretical values of total porosity and mean pore diameter as well as predicted values of pore diameter are linked with lines while samples (1-3, 4-7, 8-11, and 12-15) are arranged regarding their decreased porosity. The results show that nonwovens with similar porous structure and lower porosity also have lower pore diameter. The experimental values of total porosity are for some samples not in a good agreement with theoretical ones, while samples which should have the highest porosity actually have the lowest (samples No. 1, 8, and 12). The reason may lie in fact, that these samples contain more closed pores which are not detectable with mercury porosimetry.
The results show that the theoretical values of porosity parameters deviate from experimental ones on average by 8.0% (min 0.0%, max 15.4%) for total porosity and by 19.7% (min 2.9, max 57.3%) for pore diameter, whilst the predicted values, calculated using Equations 36-37, are in better agreement with the experimental ones. The mean predicted error is: 1.1% (from 0.0% to 4.4%) for the total porosity and 1.9% (from 0.0% to 12.4%) for the average pore diameter. The correlation coefficients between the predicted and experimental values are 0.9024 and 0.8492 for the total porosity and the average pore diameter, respectively. Scatter plots of the experimental and predicted values for porosity parameters, are depicted in Figure 13.

Conclusion
By a new fabric developing, there is a need to know some relationships between the constructional parameters of fabrics and their predetermined end-usage properties in order to produce fabrics with desired quality. Fabric constructors develop a new fabric construction on the basis of their experiences or predictive models using different modelling tools of which deterministic and nondeterministic are distinguished. In general, the models obtained by deterministic modelling tools are the results of strict mathematical rules while in the case of models obtained by nondeterministic modelling tools, there are no precise, strict mathematical rules. Our study focused on the development of predictive models based on the genetic methods, e.g. genetic programming and genetic algorithms, in order to predict some porosity parameters of woven and nonwoven fabrics. Predictive models of the: 1. area of macro-pore cross-section and macro-pore density of woven fabrics based on the constructional parameters of woven fabrics (yarn linear density, weave factor, fabric tightness, denting), image analysis as testing method of porosity measurements, and genetic programming, and 2. total porosity and mean pore diameter of nonwoven fabrics based on the constructional parameters of nonwoven fabrics (fibre linear density, fabric mass per unit area, fabric thickness), mercury intrusion porosimetry as testing method of porosity measurements, and genetic algorithm, were developed. Open porosity and equivalent pore diameter of woven fabric were also predicted using values calculated on the basis of predictive models of the area of macro-pore cross-section and pore density, and known mathematical relationships. All proposed predictive models were created very precisely and could serve as guidelines for woven/nonwoven engineering in order to develop fabrics with the desired porosity parameters.
In general, for prediction of porosity parameters of woven or nonwoven samples both modelling tools can be used, e.g. GA and GP. Usually, GP method is used for more difficult problems. Our purpose was to show usability and effectiveness of both methods. By woven fabric modelling, the range of porosity parameters' measurements was substantial larger with more input variables when compared to the nonwoven fabrics (and this means more difficult problem), so the GP was used as modelling tool. By GP modelling, the models are developed in their symbolic forms, thus more precise models are developed in regard to the GA modelling, where only coefficients of prespecified models are defined. At the same time, for GP modelling more measurements data are desired for better model accuracy, while by GA modelling good results are achieved by lower number of measurements (in our case 27 measurements were available for woven fabrics and only 15 for nonwoven fabrics). The advantage of GP modelling is its excellent prediction accuracy, while its disadvantage is the complexity of the developed models. In general, by GA modelling, the developed models are simple but less accurate.