New Geometrical Modelling for 2D Fabric and 2.5D Interlock Composites

: A new geometrical modeling tool has been developed to predict the elastic stiffness properties of 2D orthogonal and 2.5D woven interlock composites. The model estimates the change in performance due to changes in the ordering weaving parameters of the 2.5D weave architecture. Analysis results were validated compared to other models developed in published articles and the literature. Numerical analysis was performed to evaluate the accuracy of the results from the proposed models. These results demonstrate the effectiveness of the models presented by comparisons with experimental results, showing that the model could replicate the mechanical behaviors of 2D fabric and 2.5D interlock composite laminates for predicting 2D textile structures and 2.5D interlock composites with different types, shapes, and conditions. The model presented in this paper is able to replicate the behavior of woven composites of ﬁber reinforced with various types.


Introduction
Woven compounds are being increasingly considered for lots of applications because they provide ease in making complex geometries, but the mechanical properties of the different weave material supports are less visible than non-woven (angle-ply) laminates [1]. Note that there is no mathematical modelling of structural fiber-reinforced composite structures, though one measurement was obtained here using estimates in which woven fibers were formed from two orthotropic unidirectional fiber results in a curved twist using a weave. Recently, there has been in increased focus on exploring integrated and complete mechanical properties and methods of textile composites tested for uniaxial or biaxial tension, pressure, flexibility, and short-beam cutting [2][3][4][5]. Various analytical techniques have been developed to predict the elastic properties of representative volume elements (RVEs) of textile composites that include 2D and 3D weave composite. Woven composites can provide a potential solution to the basic limitations of traditional laminated composites: delamination and production that requires more workers. The addition of binder yarns provides through-thickness reinforcement, leading to highly advanced interlaminar materials and fabric binding to allow near-net-shape preforms to be woven and handled [6].
However, apart from these benefits, the use of 2.5D woven compounds is very limited in niche applications. One of the main reasons for this is the lack of predictive numerical tools, which limit their use in the early stages of design. Another method of production is to cut a simple three-dimensional weave, consisting of two two-dimensional (2D) fabrics connected by twisted loops of yarn, to form a 'hairy' fabric. These 2.5D fabrics are then coated with epoxy resin in a standard way, laminated, and treated with autoclave [7]. Woven Composites (2.5D-WC) not only have a higher delamination resistance compared to 2D laminated composites, but also have a simpler structure than 3D textile composites.

State of the Art
Various analytical techniques have been developed to predict the expansive properties of volume representation (RVEs) of textile composites that include 2D and 2.5D weave composites. Bogdanovich and Pastore [9], Whitcomb et al. [10], and Potluri and Thammandra [11] have conducted a complete review of the topic. Based on stiffness averaging and rigidity ratio, Ishikawa and Chou [12][13][14] managed a two-orthogonal set of fabric compositions such as laminate with two laminae placed at 0/90 position. They proposed different models for predictable expansion, such as (1) a mosaic model (where fiber continuity and retraction was left), (2) a crimp model (where fiber continuity and undulation geometry were defined using the functions of -trigonometric), and (3) bridge coupling model (rigid or compliant dimensions in two orthogonal directions and suitable for satin weave combinations). Naik et al. [15,16] analyzed fiber retreat to the two orthogonal directions and used detailed trigonometric functions to extend the crimson models of Ishikawa and Chou (i.e., parallel two-dimensional series models and series-parallel models) to test expandable structures different combinations of PWF. Saka and Harding [17] and Carey et al. [18] developed rectilinear and sinusoidal crimp models to predict expansion structures using the classic laminate plate theory applied in the context of an expandable base and using the same power system. Tong et al. [19] take fiber as a curved rod based on stretch elastic or shear base in a woven composite structure and present a curved bark model to investigate the impact of fiber flexibility on the strength of woven composite material based on the concept of classical laminate plate theory. Kollegal and Sridharan [20] employed hierarchic works to determine the original curved profile of the strings and plate displays and introduced the micro-plate model for PWF analysis.
A significant amount of research was performed using the principles of differential stress [11,21] or FE (FEA) analysis [22][23][24] to analyze the expandable tracts of woven compounds and a detailed stress field across the RVE from the limited equipment. Regarding basic properties, in order to produce FE models, precise vertical segments or continuous mathematical functions are used to demonstrate in detail the ideal PWF geometry and to represent the thread path and cross section shape.
Mesoscopic simulations of transverse compaction have attracted recent interest. Geometric modeling methods have been developed using fabric matching software such as TexGen [25,26] or WiseTex [27,28], which have been used in various articles [29][30][31][32][33]. Multi-chain digital element-based methods have also been used [34][35][36][37][38]. Green et al. [39,40] proposed a numerical procedure for creating mesostructural geometric models in simulated fabrics by combining appropriate TexGen models with multi-chain digital materials. Nguyen et al. [41] followed the method developed by Hivet and Boisse [42] to produce mesostructural geometric models of 2D woven fabrics to mimic the opposite combination. Numerical simulations of hypo elastic material based on fiber rotation were well matched to experiments. Alternatively, based on the discrete homogenization method, Goda et al. [43] and Rahali et al. [44] can reproduce the functional properties of 2.5D machines and 3D interlock textiles. Zhang et al. [45] developed a model based on the optical aspect ratio to predict mechanical responses and failure location under uniaxial and biaxial loading.
However, these methods are based on appropriate geometrical models, therefore cannot explain mesostructural changes in fiber structure, i.e., a detailed study of how fibers degenerate under bonding. Consideration and simplicity are often needed in numerical matching, which reduces calculation costs, facilitates the modeling process, and provides computer estimations. However, this simplification may prevent the investigation of other events occurring in real systems, especially in the case of mesoscopic and microscopic disorders.
Specific studies of mesoscopic transformation of fibrous reinforcement can also be performed with minimally invasive imaging such as X-ray micro-computed tomography (Micro-CT). With this process, it is possible to measure, model, and analyze the geometric structure of fiber cables. Desplentere et al. [46] and Schell et al. [47,48] were among the first to demonstrate the power of Micro-CT to show the geometric features of 3D fabrics on a mesoscopic scale. In Pazmino et al.'s study [49], the geometric features are directly measured by the fiber pulse from small tomographic images to read a single layer of 3D non-crimp woven fabric. Naouar et al.'s [50,51] mesostructural geometry models successfully reconstructed 2D woven fabric and 3D orthogonal fabric from small tomographic images. Wang et al. [52] studied longitudinal compression and the Poisson number of fiber cables followed a similar pattern. Badel et al. [53] performed both numerical simulations and tomographic analyzes of 2D woven fabrics under biaxial tension and inplane shear deformation.
Combined with digital volume integration, Mendoza et al. [54] developed a mathematical model for the distortion of woven compounds, which was recently used to measure the flexibility caused by the production process [55]. However, no research studies to date have reported on the analysis of variables based on meso-structural models of fiberglass fabrics separated by these simplified methods at different levels of assembly. Methodology:

2D Fabric
The user-whether a student or a fellow researcher-can themselves design a fabric composite. To construct the whole structure of composite as he or she likes from the beginning and to make the whole process "user friendly", the process took the shape of building a puzzle by a simple "click and drop" game as shown in the figure below ( Figure 1). The user-whether a student or a fellow researcher-can themselves design a fabric composite. To construct the whole structure of composite as he or she likes from the beginning and to make the whole process "user friendly", the process took the shape of building a puzzle by a simple "click and drop" game as shown in the figure below ( Figure 1).

2D Fabric
The user-whether a student or a fellow researcher-can themselves desig composite. To construct the whole structure of composite as he or she likes from ning and to make the whole process "user friendly", the process took the shape of puzzle by a simple "click and drop" game as shown in the figure below ( Figure 1).  As 2D woven fabrics are made by interlacing yarns in a weaving loom, yarns are divided into two components: one called the warp, running along the length of the loom, and the other is the weft, running in the cross direction. The woven structure is characterized by orthogonal linking of two sets of threads, called warp and weft yarn. The warp threads are aligned with the direction of the fabric leaving the weaving equipment, which is also called the warp direction. As shown in Figure 2, 2D woven fabric has two yarn sets as warp (0 • ) and filling (90 • ) and interlaced to each other to form the surface.
Textiles 2022, 2, FOR PEER REVIEW 5 threads are aligned with the direction of the fabric leaving the weaving equipment, which is also called the warp direction. As shown in Figure 2, 2D woven fabric has two yarn sets as warp (0°) and filling (90°) and interlaced to each other to form the surface. To start from the beginning, the user is free to choose the settings of the puzzle in each dimension, i.e., number of columns and rows, and the number of fibers, as shown in Figure 3. Plain weave is a very common and very strong basic weave where each warp fiber passes alternately below and above each weft. The fabric is symmetrical and has good stability and logical porosity. However, they are the most difficult to weave, and the high quality of fiber crimp offers lower mechanical features compared to other weaving styles.  To start from the beginning, the user is free to choose the settings of the puzzle in each dimension, i.e., number of columns and rows, and the number of fibers, as shown in Figure 3.
Textiles 2022, 2, FOR PEER REVIEW 5 threads are aligned with the direction of the fabric leaving the weaving equipment, which is also called the warp direction. As shown in Figure 2, 2D woven fabric has two yarn sets as warp (0°) and filling (90°) and interlaced to each other to form the surface. To start from the beginning, the user is free to choose the settings of the puzzle in each dimension, i.e., number of columns and rows, and the number of fibers, as shown in Figure 3. Plain weave is a very common and very strong basic weave where each warp fiber passes alternately below and above each weft. The fabric is symmetrical and has good stability and logical porosity. However, they are the most difficult to weave, and the high quality of fiber crimp offers lower mechanical features compared to other weaving styles.  Textiles 2022, 2, FOR PEER REVIEW 5 threads are aligned with the direction of the fabric leaving the weaving equipment, which is also called the warp direction. As shown in Figure 2, 2D woven fabric has two yarn sets as warp (0°) and filling (90°) and interlaced to each other to form the surface. To start from the beginning, the user is free to choose the settings of the puzzle in each dimension, i.e., number of columns and rows, and the number of fibers, as shown in Figure 3. Plain weave is a very common and very strong basic weave where each warp fiber passes alternately below and above each weft. The fabric is symmetrical and has good stability and logical porosity. However, they are the most difficult to weave, and the high quality of fiber crimp offers lower mechanical features compared to other weaving styles.  Twill is a type of textile weave with a pattern of corresponding diagonal ribs. It can be seen by observing the presence of diagonal lines that are pronounced along the width of the fabric. It has a higher resistance to cracking than plain weave because it has fewer wires that connect to each area, hence the greater the internal flow rate. In addition, two strings will carry the load when the fabric is torn.

Sateen/Satin [5-end, 8-end]
Sateen is a fabric that usually has a shiny surface and a hidden back, one of the three basic types of fabric weaving, seamless weaving, and twill. Four or more full-length or single-stranded strands floating over a straight rope, and four straight strands floating on a single weft thread indicate satin weaving. Floating strings are missed interfaces, where the warp thread lies on the weft in a satin with a straight face and where the weft thread lies on straight strands on satin with a weft face.

Basket Weave [2/2, 4/4]
Basket weave is known also as Hopsack and Matt Weave. Hopsack weave, a variant of plain weave, uses two or more warp and/or two or more weft fibers joined as a single thread. This weave is obtained by doubling or repeating the combined points of a plain weave in both the direction of the wrap and weft. These fabrics are made of two or more strings placed in the same shed. The stitching pattern is similar to a plain weave, but two or more strands follow the same parallel pattern. Fabrics designed by Matt are flexible and cannot wrinkle as there are a few cutters that are a square inch. Fabrics look flatter than conventional weaving fabrics. However, longer floats are more flexible. Matt fabric has great anti-tear properties. Matt's design tends to offer more smooth fabrics. For a repetitive Matt weave size, the warp numbers and weft threads are equal.

2D Fabric Hybrid Composites
Hybrid yarns are where two or more fibers are brought together to combine the performance and aesthetic of both, as shown in Figure 8. Although environmental concerns are encouraging many manufacturers to focus on the use of a single fiber, there are still and will remain reasons for bringing fibers with different qualities together, particularly for performance and health and safety applications. Plain weave is a very common and very strong basic weave where each warp fiber passes alternately below and above each weft. The fabric is symmetrical and has good stability and logical porosity. However, they are the most difficult to weave, and the high quality of fiber crimp offers lower mechanical features compared to other weaving styles.

Twill [2/1, 3/1]
Twill is a type of textile weave with a pattern of corresponding diagonal ribs. It can be seen by observing the presence of diagonal lines that are pronounced along the width of the fabric. It has a higher resistance to cracking than plain weave because it has fewer wires that connect to each area, hence the greater the internal flow rate. In addition, two strings will carry the load when the fabric is torn.

Sateen/Satin [5-end, 8-end]
Sateen is a fabric that usually has a shiny surface and a hidden back, one of the three basic types of fabric weaving, seamless weaving, and twill. Four or more full-length or single-stranded strands floating over a straight rope, and four straight strands floating on a single weft thread indicate satin weaving. Floating strings are missed interfaces, where the warp thread lies on the weft in a satin with a straight face and where the weft thread lies on straight strands on satin with a weft face.

Basket Weave [2/2, 4/4]
Basket weave is known also as Hopsack and Matt Weave. Hopsack weave, a variant of plain weave, uses two or more warp and/or two or more weft fibers joined as a single thread. This weave is obtained by doubling or repeating the combined points of a plain weave in both the direction of the wrap and weft. These fabrics are made of two or more strings placed in the same shed. The stitching pattern is similar to a plain weave, but two or more strands follow the same parallel pattern. Fabrics designed by Matt are flexible and cannot wrinkle as there are a few cutters that are a square inch. Fabrics look flatter than conventional weaving fabrics. However, longer floats are more flexible. Matt fabric has great anti-tear properties. Matt's design tends to offer more smooth fabrics. For a repetitive Matt weave size, the warp numbers and weft threads are equal.

2D Fabric Hybrid Composites
Hybrid yarns are where two or more fibers are brought together to combine the performance and aesthetic of both, as shown in Figure 8. Although environmental concerns are encouraging many manufacturers to focus on the use of a single fiber, there are still and will remain reasons for bringing fibers with different qualities together, particularly for performance and health and safety applications.

2D Fabric Hybrid Composites
Hybrid yarns are where two or more fibers are brought together to combine th performance and aesthetic of both, as shown in Figure 8. Although environmental con cerns are encouraging many manufacturers to focus on the use of a single fiber, there ar still and will remain reasons for bringing fibers with different qualities together, partic ularly for performance and health and safety applications.

Interlock Composite 2.5D
The 2.5D, or so-called 3D, interlock composite is a special fabric combination in which mechanical properties are closely related to its structure. The geometry of the interlock is complex and the number of possible structures is infinite. Fabric architecture depends on the extraction, crimp, and size of the fiber fibers. It is made up of a system of two woven and interconnected yarns, straight yarns, and weft (or filling) yarns. Interlacing caused by bending is called "tow crimp".
The shapes in the "puzzle" represent the filler element and the bounder (interlock fiber) of the structure, as shown in the figure below (Figure 9), in 2D view.
terlock is complex and the number of possible structures is infinite. Fabric architecture depends on the extraction, crimp, and size of the fiber fibers. It is made up of a system of two woven and interconnected yarns, straight yarns, and weft (or filling) yarns. Interlacing caused by bending is called "tow crimp".
The shapes in the "puzzle" represent the filler element and the bounder (interlock fiber) of the structure, as shown in the figure below (Figure 9), in 2D view. 2.5D angle closing fabrics can be divided into two types: thick fabrics and fabrics from layer to layer. This type of combination is characterized by its intricate structure. The cell unit consists of a warp weaver and weft fibers attached to 90° plane (xy) ( Figures  10 and 11). A 2.5D angle woven joint includes warp that binds to warp straight cords by attaching warp threads. Warp wires can be tied to different depths where different layout arrangements could be used to produce a wide range of these types of compounds. Tight fabric is a multi-layered fabric where warp weavers move from one piece of fabric to another, holding all the layers together. The layer-to-layer fabric as shown in Figure 10 is a multi-layer fabric that the warp weavers move from one layer to the nearest layer and back. A collection of warp woven together holds all the layers of fabric. In addition, complex geometry as shown in Figure 11, fractional volume fraction, cable volume, and inclination angle of warp strands are not able to allow the structural properties of specific applications. In other words, designers can replicate the efficient performance of fabrics for the necessary mechanical properties.

1
Layer to Layer Angle Interlock  2.5D angle closing fabrics can be divided into two types: thick fabrics and fabrics from layer to layer. This type of combination is characterized by its intricate structure. The cell unit consists of a warp weaver and weft fibers attached to 90 • plane (xy) (Figures 10 and 11). A 2.5D angle woven joint includes warp that binds to warp straight cords by attaching warp threads. Warp wires can be tied to different depths where different layout arrangements could be used to produce a wide range of these types of compounds. Tight fabric is a multi-layered fabric where warp weavers move from one piece of fabric to another, holding all the layers together. The layer-to-layer fabric as shown in Figure 10 is a multi-layer fabric that the warp weavers move from one layer to the nearest layer and back. A collection of warp woven together holds all the layers of fabric. In addition, complex geometry as shown in Figure 11, fractional volume fraction, cable volume, and inclination angle of warp strands are not able to allow the structural properties of specific applications. In other words, designers can replicate the efficient performance of fabrics for the necessary mechanical properties.
lacing caused by bending is called "tow crimp".
The shapes in the "puzzle" represent the filler element and the bounder (interlock fiber) of the structure, as shown in the figure below (Figure 9), in 2D view. 2.5D angle closing fabrics can be divided into two types: thick fabrics and fabrics from layer to layer. This type of combination is characterized by its intricate structure. The cell unit consists of a warp weaver and weft fibers attached to 90° plane (xy) ( Figures  10 and 11). A 2.5D angle woven joint includes warp that binds to warp straight cords by attaching warp threads. Warp wires can be tied to different depths where different layout arrangements could be used to produce a wide range of these types of compounds. Tight fabric is a multi-layered fabric where warp weavers move from one piece of fabric to another, holding all the layers together. The layer-to-layer fabric as shown in Figure 10 is a multi-layer fabric that the warp weavers move from one layer to the nearest layer and back. A collection of warp woven together holds all the layers of fabric. In addition, complex geometry as shown in Figure 11, fractional volume fraction, cable volume, and inclination angle of warp strands are not able to allow the structural properties of specific applications. In other words, designers can replicate the efficient performance of fabrics for the necessary mechanical properties. 1 Layer to Layer Angle Interlock

Construction Testing
Tests should be performed as shown in Figures 10 and 11 to validate the sh  These tests obtained by logical geometrical construction of the 2.5D interlock structure are as constrains that define whether the design (structure) from the user is validapplicable-as shown in Figure 10.
As named, the function of the tests are to check if there is a continuity of structure and if there is a possibility to add a shape in the puzzle as shown in Figure 12. For example, Shape 3 cannot be similar to Shape 2 s incident and reflected interlock at the same time. The puzzle should be symmetrical if the user is constructing a symmetrical shape. If not, the test to check symmetry can be deactivated. Shapes 2, 3, and 4 cannot be repeated vertically. Shapes 2 and 4 cannot be repeated horizontally. Shapes 2 and 4 should be separated only by Shape 3. Shape 3 cannot follow Shape 1 or 2. At least one instance of Shape 1 should separate Shapes 2 and 4 horizontally. In addition, there are other structural and geometrical conditions and boundaries related to the shape.

Construction Testing
Tests should be performed as shown in Figures 10 and 11 to v suggested by the user. These tests are: These tests obtained by logical geometrical construction of the 2.5 ture are as constrains that define whether the design (structure) from applicable-as shown in Figure 10.
As named, the function of the tests are to check if there is a conti and if there is a possibility to add a shape in the puzzle as shown in ample, Shape 3 cannot be similar to Shape 2's incident and reflected int time. The puzzle should be symmetrical if the user is constructing a sym not, the test to check symmetry can be deactivated. Shapes 2, 3, and 4 c vertically. Shapes 2 and 4 cannot be repeated horizontally. Shapes 2 separated only by Shape 3. Shape 3 cannot follow Shape 1 or 2. At lea Shape 1 should separate Shapes 2 and 4 horizontally. In addition, the tural and geometrical conditions and boundaries related to the shape. Note that all these conditions/boundaries are tested automatically code created by the author.
Moreover, if there is any structural mistake made by the design will identify the mistake with the relation of the mentioned tests. In a Note that all these conditions/boundaries are tested automatically by the JavaScript code created by the author.
Moreover, if there is any structural mistake made by the designer (user) the code will identify the mistake with the relation of the mentioned tests. In addition, the code will specify exactly the place of the error with the messages: "Pattern problem at column 2-Follow up problem at (2,2)-Invalid structure" for the designer have to reshape its structure to continue and the sentences "No pattern problem-No follow up problem-Valid structure", as shown in Figure 13 will specify exactly the place of the error with the messages: "Pattern problem at column 2-Follow up problem at (2,2)-Invalid structure" for the designer have to reshape its structure to continue and the sentences "No pattern problem-No follow up problem-Valid structure", as shown in Figure 13

Analytical Modelling
First, the user should fill fibers and matrix parameters such as Young's and shear modulus (Eᵢ,Gᵢ) and fiber volume fraction (Vf) as shown in Figure 14 as a user interface to the analytical modelling.

Analytical Modelling
First, the user should fill fibers and matrix parameters such as Young's and shear modulus (E i ,G i ) and fiber volume fraction (Vf) as shown in Figure 14 as a user interface to the analytical modelling. Then, to calculate the characteristics of the key elements and the contribution rate, they make a macroscopic layer and subsequently the whole unit cell. The current approach creates expressions at the micro level with the aim of calculating more representative volume fractions of a group of elements to the layer to improve the elastic Then, to calculate the characteristics of the key elements and the contribution rate, they make a macroscopic layer and subsequently the whole unit cell. The current approach creates expressions at the micro level with the aim of calculating more representative volume fractions of a group of elements to the layer to improve the elastic stiffness speculation compared to existing analytical modeling methods that use the 2.5D woven composite as a composite component containing layers of unidirectional elements (which are fibrous tows encased in resin). The new modelling approach creates expressions that discretize the unit cell into elements To calculate the macroscopic stiffness of the whole unit cell, it must be broken down or discretized to the microscale. The microscale looks at the individual elements that make a cell and then the layers that make up the macroscopic unit cell. Therefore, the unit cell undergoes the first level discretization into layers, then the second level discretization into the individual elements that make up the layer. An element can be an individual, filler, binder, or matrix region within a layer. Having determined the stiffness of each element, the stiffness of whole row can be found. Once the stiffness of all rows is known, the model formulates and calculates the stiffness of the whole unit cell and finally outputs the elastic constants. The prediction of the unit cell or macroscopic stiffness begins with the calculation at the micro scale (the constituent elements within a cell).
The algorithm applied is as follows: Textiles 2022, 2, FOR PEER REVIEW The fiber matrix (S) strength can be easily calculated. In the pocket of a pure matrix, it is generally regarded as isotropic substances. When the elastic modulu Poisson's ratio of resin matrix are given, it is easy to determine its basic relationshi easy to obtain their compliance matrix in the local material coordinate systems (1-2 (x-y-z) by inverting their stiffness matrix (Sij) as shown in Equations (1) and (2). The final step is the summation of the found matrices in column and in rows to find the Overall Stiffness matrix The fiber matrix (S) strength can be easily calculated. In the pocket of a pure resin matrix, it is generally regarded as isotropic substances. When the elastic modulus and Poisson's ratio of resin matrix are given, it is easy to determine its basic relationship. It is easy to obtain their compliance matrix in the local material coordinate systems (1-2-3) to (x-y-z) by inverting their stiffness matrix (Sij) as shown in Equations (1) and (2).
Each type of cell unit consists of two types of elements, namely fiber and a pure resin matrix. First, the lamina compliance matrix was calculated. Like most micro-mechanical models, the fibers and the resin matrix assume transversely isotropic, and both of them assume to be linearly elastic in the model. In order to achieve the elastic properties of the composites, the Chamis-proposed fiber-matrix ROM was selected to compute the engineering elastic constants. The Chamis micro-mechanical is a widely used and reliable model, providing an equation of all five independent structures that stretch like this [56]: ν 23 = V f ν f 23 + V m (2ν m − ν 12 (E 22 /E 11 )) (7) where V f is the fiber volume fraction, E f 11 is the Young's elastic modulus of the fiber in principle axis 1, E f 22 is the Young's elastic modulus of the fiber in principle axis 2, G f 12 is the longitudinal shear modulus of the fiber, G f 23 is the transverse shear modulus of the fiber, ν f 12 is the primary Poisson's ratio of the fiber, and E m , ν m , and G m represent the Young's elastic modulus, Poisson's ratio, and shear modulus of the matrix, respectively.
The compliance matrix of the fiber can be easily calculated. For the pure resin matrix pocket, it is generally regarded as isotropic material. When the elastic modulus and Poisson's ratio of resin matrix are given, it is easy to determine its basic relationship.
After that, we can obtain the stiffness in the global coordinate system by transforming the stresses and strains with the generalized transformation matrix as the next form: where k is the number of unit cells in the Puzzle structure and the angle defined as the cosines of the angle between the axes of the local and global coordinate systems before and after rotation: s = Sin (θ) (12) where θ is the angle of weaving for the yarns. A unidirectional cell falls under the orthotropic material category. As the cell is thin and does not carry any out-of-plane loads, one can assume plane stress conditions for the cell.
Therefore, assuming σ3 = 0, τ23 = 0, and τ31 = 0, calculation of element contribution is required so that the stiffness of the whole cell can be calculated as well as the overall unit cell stiffness. In our case, the contribution is related to the percentage in volume of the bounder (interlock fiber) with reference to the matrix volume (for example 30% bounder, 70% resin). Mij = (p) C bK ij + (1 − p) C K mij (13) where [C b ] and [Cm] are the bounder and resin stiffness matrix in the global coordinate system, respectively, p is the percentage in volume of the bounder and resin in the composite structure, and k is the number of element in the whole structure. Calculating the total matrix/row by summation of matrices (series summation) as shown in Equation (14): where n is the number of rows in the puzzle and [M i ] are the inverse of matrices found in Equation (14). The final step is the summation of the found matrices in column and in rows as shown in Equation (15): where m is the number of columns in the puzzle and [Mt-1 i ] are the matrices found in Equation (15).

2D Fabric
To ensure the effectiveness of this model based on this work. In comparison between current results modeling, experimental data, and previously developed multi-scale modelling, four different examples are generated. As studies, different mechanical features are considered and results are released. Examples were chosen from the specimens studied by J.J. Xiong and his co-workers [57] to calculate the engineering elastic constants of 2D Fabric Woven composites.
The fabric specifications and mechanical properties of the above four kinds of textile composites are listed in Table 1

2.5D Interlock
The composites are made from carbon fibers T300J and resin matrix RTM6 (Table 2). Two examples were chosen from the specimens studied by Hallal and his co-workers on 3SHM (3 stages homogenization method) [58]. The architecture of these composites involves only warp weaver yarns and weft straight yarns, as shown in Figures 10 and 11. The REV of the composite-H2 is composed of 6 warp yarns and 12 weft yarns, while the REV of the composite 71 is composed of 3 warp yarns and 24 weft yarns. The REV of composite 69 contains 6 warp yarns and 6 weft yarns. Warp yarns have a mean inclination angle equal to 29 • for the linear longitudinal part of weft yarns. The warp yarns have a linear plus undulated longitudinal parts and a flattened elliptical cross-section. They are interlocked with weft yarns in two steps. Undulated parts have a mean inclination angle equal to 24 • for H2 and 12 • for 71. The weft yarns have a linear longitudinal part while its cross-section is a flattened ellipse. In addition, the fiber volume fractions in both warp and weft yarns is taken equal to 0.6. Table 2 shows the mechanical properties of carbon fibers and matrix. Carbon fibers assumed transversely isotropic material, which gives the following assumption: Table 2. Mechanical properties of carbon fibers and matrix [58].

Carbon Fibers
Ef 11

Results and Discussion
With the purpose of validating the developed modeling technique, the results were compared with experimental data published in the open literature. A comparison between the results of the present modeling, experimental data, and previously developed multiscale modeling was conducted. As case studies, different mechanical properties were considered, and the results are shown in Tables 3-5, each of which shows a clear comparison between the experimental results and the presented model.  For example, Table 3 shows the comparison of longitudinal and Young's moduli between four different materials and fibers studied by J.J. Xiong [57]. On the other hand, Tables 4-6 show a clear comparison for the section of 2.5D interlock composite with its different types.

2D Fabric
As shown in Table 3, the comparison was focused on tension moduli from the specimens studied by J.J. Xiong and his co-workers [57]. The method present accurate and precise values with the percentages of error varied from 0.44% for longitudinal Young's modulus Ex, in E-glass/epoxy-II experiment up 1.91% for EW220/5284.

5D Interlock
The following table presents the results found by a previous analytical model from the reference validation mentioned earlier The comparison was made to the engineering properties found using the code. A few notes should be taken into consideration: The model does not have curvature shape between Shapes 2, 3, and 4. The structure is in the coordinate system (x, y) by the reference. The dimensions of fiber are considered to be a tube shaped form and not elliptical. By comparing the result, we see that there is a validation to the code and the new model with a marginal error due to the curvature shape missing and the real/theoretical shape relation between the numerical, analytical, and geometrical model. It was observed from Table 4 that the Young's modulus in direction 2 is the most accurate parameter, up to 1% in comparison 3SHM found in Hallal [58], and that the maximum average error is shown for the Young's modulus in direction 1 with 10% of error. From Table 4, the percentage of error of transverse Young's modulus is much higher than for the longitudinal Young's modulus, with a value of 4.13% as the percentage of error of longitudinal Young's modulus was 0.614%. In addition, the percentages of error for shear modulus in Tables 4 and 5 were 1.12% and 1.21%, respectively. Finally, from Table 6, the percentage of error of transverse Young's modulus was also much higher than for the longitudinal Young's modulus, with a value of 5.07% and the percentage of error of longitudinal Young's modulus was 1%. The percentage of error for shear modulus in Table 6 was 1.31%.
In addition to the high accuracy level of the current modelling technique, the very short runtime required of the modelling renders the developed multi-scale modelling as a cost-effective computational tool for estimating the Young's moduli of laminate composites. The results showed that the present model could be used to effectively evaluate the elastic properties for laminate composites.

Conclusions
A geometrical-modelling tool has been presented to predict the elastic stiffness characteristics of 2D fabric and 2.5D woven interlock composites with the ability to assess change in performance as a consequence of altering weaving parameters. The models presented in this paper are able to reproduce the behavior of woven composites formed by different fiberreinforced types. This approach has been validated against experimental data produced independently of this work for orthogonal interlock weaves and compared to existing modelling approaches. The present model performs better in all predictions compared to the existing modelling efforts. Good agreement was observed among numerical results, with small differences in the longitudinal strains, as shown in the results section that the shear modulus was the most accurate parameter up to 1% in the 2.5D interlock section, while the error in 2D fabric section was lower with values such as 0.44% in tension moduli. As for the longitudinal Young's modulus, the maximum average error is shown with 10% but with an average percentage of error of transversal strains of 4%. This may be due to a lack of precision in the mechanical properties of the components, i.e., the fibers, which are used as input data. The percentage is related to the lack of an accurate and precise database structure of the used material, but it shows a great ability to use the algorithm.
The suggested hypotheses are the simplest that can be applied for determining the engineering data with all required information, and specifications from the user are presented in a "user friendly" form that can make studying and searching a new level of excitement and interactivity. As shown in the methodology section, applying the algorithm is easy and simple. All this research can be accessed easily with all its codes and models, along with its regulations and explanations needed. Working on these codes will eliminate the hard work caused by the computational efforts using software such as Abaqus or Ansys, etc.
At any given time, the user can calculate required forces to have a known deformation or to predict the deformation for a known force; in addition, they would have easy access to all engineering properties such as stiffness or compliance matrices. A future scope is now set in each section to have more research and data analysis by using these codes and models without any super computers and losing time, money, and effort. Any researcher can now find the optimum set of materials for all levels of woven composites materials (2D Plain weave, Twill, Satin, or Basket), along with altering materials for 2.5D interlock composites such as layer to layer or through the thickness angle interlock. The model works best on other types of combinations, and further analysis and evaluation of the models is required. The introduced framework may apply to other fiber conditions as part of a fiber volume, but these conditions will require further development.

Conflicts of Interest:
The authors declare no conflict of interest.