Stress Field Evaluation in Orthotropic Microstructured Composites with Holes as Cosserat Continuum

It is known that the presence of microstructures in solids such as joints and interfaces has an essential influence on the studies of the development of advanced materials, rock mechanics, civil engineering, and so on. However, microstructures are often neglected in the classical local (Cauchy) continuum model, resulting in inaccurate descriptions of the behavior of microstructured materials. In this work, in order to show the impact of microstructures, an implicit ‘non-local’ model, i.e., micropolar continuum (Cosserat), is used to numerically investigate the effects of direction and scale of microstructures on the tension problem of a composite plate with a circular hole. The results show that distributions of field variables (such as displacements and stresses) have an obvious directionality with respect to the microstructures’ direction. As the scale of microstructures increases, such a direction effect becomes more evident. Unlike the isotropic material where stress concentration occurs at the vertex of the hole and the stress concentration factor is close to 3, for the microstructured composite, the stress concentration can be observed at any location depending on the microstructures’ directions, and the concentration factor can exceed 3 to a maximum close to 9 as the increasing scale of microstructures. In addition, differences in the mechanical behavior between Cosserat and Cauchy models can be also observed; such differences are more evident for the material showing a pronounced orthotropic nature.


Introduction
Microstructure is one of the most critical factors that involves many kinds of materials such as rock, ceramic, alloy, human cortical bone, etc. [1][2][3][4]. As an internal structure, microstructure can play a crucial role in determining the gross behavior and mechanical response of materials [5]. However, the microstructures in materials distribute randomly with different lengths and directions, which complicates the understanding of the material's response. In general, materials with microstructure have weaker strength than in intact materials [6,7]. Guo et al. [8] experimentally investigated the effect of bedding angle in phyllite under unloading confining pressures, where the rock bedding joints can be regarded as microstructures. They found that the rock shows different strengths as the bedding angle changes. Numerical research on the tunnel surrounding rocks with different inclination angles also shows a directional effect on the distributions of displacements and stress around the tunnel [9]. By reviewing indentation tests at the micron scale, Bauer et al. [10] demonstrate that an obvious length scale effect (i.e., non-locality) can be found when the material's intrinsic length scales are comparable with the dimension of specimens. Therefore, it is of importance to describe the macroscopic response of these materials by considering the influences of the microstructures.
The existence of microstructure results in the heterogeneity characteristic of materials. There are various methods that can be used to model the behavior of microstructured materials. Discrete modeling with interactions of each constituent in materials is a good option because it can produce an accurate result; however, this approach is often computationally cumbersome [11][12][13]. Alternatively, homogenizing the heterogeneous material as an equivalent continuum could be an efficient approach because it is faster and takes less computational cost [14]. Nevertheless, the application of this approach depends on selections of the homogenization method and macroscopic continuum theory that need to reveal the presence of microstructures. As is known, the classical Cauchy continuum may have disadvantages in describing the gross behavior of microstructured materials since it lacks in internal length descriptions [15,16]. This calls for the application of the non-local continuum theory, as this approach can reveal the presence of internal lengths [17]. In the non-local theory, internal lengths can be represented by adding extra degrees of freedom or parameter as internal variables, corresponding to the so-called implicit and explicit non-local descriptions, respectively [18,19].
The Cosserat continuum theory is a widely used implicit non-local model to investigate the microstructured material's behavior. After the completed mathematical foundations of the micropolar continuum was achieved, this theory became very popular (since the 70s). The Cosserat model introduces to each material point an extra degree of freedom, termed microrotation, which is different from the local rigid rotation (i.e., macrorotation). As a result, the stress and strain fields become asymmetric in this model, which is different from the classical Cauchy model with symmetric measurements. Moraes et al. [1] found that the asymmetrical property in the Cosserat model can be helpful to improve the description of the mechanical behavior of the materials such as rocks. The asymmetric strains also correspond to the relative rotation between the microrotation and macrorotation. Pau and Trovalusci [20] found that the relative rotation is significant in anisotropic materials, whereas it can be negligible in orthotetragonal materials where the internal length trends to vanish. With the advantage of keeping the memory of the microstructure, the Cosserat continuum was used to study many kinds of materials such as layered materials [21,22], fiber-reinforced materials [23,24], granular materials [25], and composites [26][27][28]. Using a homogenization process for the Cosserat continuum, Trovalusci and Masiani [29] numerically and experimentally studied the mechanical behavior of an inclined masonry structure in which microstructures (interfaces) show a different direction from the ordinary masonry structure. However, the results of this literature only showed the micropolar effect resulted from the Cosserat continuum but the effect of microstructures direction was not further discussed.
In this work, the effects of the microstructure's direction and length scale in a composite are studied to extend the understanding of microstructured materials, such as advanced materials with various microstructured and layered rocks with inclined angles. The composite considered here is made of rectangular blocks interacting with each other through their elastic interfaces, and it is homogenized as a Cosserat continuum by an energeticequivalence-based homogenization technique [30]. Thus, the characteristic of non-locality is involved in this study. Six directions and four length scales of the microstructure are investigated by the finite element method (FEM) for a tension problem of a composite plate with a circular hole. Therefore, this work also focuses on the stress concentration problem of microstructured materials. This problem has been widely reported in previous studies [31][32][33][34]. Holes in materials can induce stress concentration around it and hence reduce the mechanical properties [35]. However, a number of solutions have been carried out for holes in isotropic plates [36]. With the increasing research interests on materials, especially with microstructures, it is essential to gain a better understanding in modeling the mechanical behavior of these materials. This paper is structured as follows. After the introduction section, Section 2 introduces the Cosserat theory and its FEM implementation. Section 3 presents the model, parameters, methods, etc. used in the simulations of the tension problem for a plate with a circular hole. In Section 4, numerical simulations are conducted and results of displacements, stresses, and relative rotation are shown. The stress distribution around the hole is discussed and the simulation results are analyzed in Section 5. In the end, conclusions and remarks are drawn in Section 6.

Cosserat Continuum and Its FEM Implementation
The Cosserat continuum is considered to be a multi-scale tool [37] that can be used to investigate the mechanical behavior of materials where the microstructures and, in particular, internal lengths, play a crucial role. As an implicit' non-local continuum, it is equipped with additional degrees of freedom revealing the presence of microstructures. That is, for two-dimensional (2D) Cosserat media, each material point has two translation degrees of freedom u 1 , u 2 and an additional microrotation degree of freedom ω. The microrotation ω is an independent degree of freedom and it is different from the macrorotation θ which is defined as the skew-symmetric part of the gradient of displacement. Thus, a peculiar measurement, the relative rotation θ − ω, can be defined in this continuum. A general displacement vector for the Cosserat material point can be expressed as: where u = u 1 u 2 . Due to the introduction of ω, the tangential strains in the Cosserat model are not reciprocal, i.e., ε 12 = ε 21 , and the microcurvature component is introduced as an additional strain measure; therefore, the linear strain-displacement relation can be expressed as: where ε = ε 11 ε 22 ε 12 ε 21 contains the normal and tangential strains and χ = χ 1 χ 2 contains the microcurvatures. ∇ is the gradient operator, and With the strain measures, the stresses of the Cosserat continuum can be obtained by a linear elastic constitutive equation as: where σ = σ 11 σ 22 σ 12 σ 21 contains the normal and tangential stresses and µ = µ 1 µ 2 contains the couple stresses. The tangential stress components are also not reciprocal (σ 12 = σ 21 ). The constitutive sub-matrices A, B, and D collect the constitutive terms A ijhk , B ijh , and D ij , where i, j, h and k = 1, 2.
A detailed three-dimensional finite element formulation for the Cosserat continuum can be found in [38]. For the sake of simplicity, the displacement-based finite element implementation for 2D Cosserat theory is presented here to model the behavior of microstructured material. Firstly, in the finite element procedure, displacement and microrotation fields should be approximated by the nodal values of an element. In this study, for avoiding the element locking problem, a bi-quadratic (N u ) and a bi-linear (N ω ) shape function are, respectively, used for the displacement and microrotation approximation: (5) whereũ andω are nodal displacement and microrotation values. In the present paper, nine-node quadrangular elements are considered for an element. All nine node values are used to approximate the displacements, whereas values at four corner nodes are used for the microrotation. N u and N ω can be expressed as: Substituting Equation (5) into (2), the strain vectors become: where B ε and B χ are the derivatives of the shape functions.d is the unknown nodal values collectingũ andω. Substituting Equation (7) into Equation (4), the constitutive relations become: Now, the stress and couple stress measures can be obtained from the nodal values. Considering a domain A and boundary Γ, the principle of virtual work can be expressed as: where δ is the variational operator, b is the body force vector. t and m are the traction and couple-traction vectors applied on the boundary Γ. The components (t i and m i ) of t and m should satisfy the equilibrium at external boundary as t i = σ ij n j and m i = µ j n j , where n j is the components of the outward unit normal to the boundary. Substituting Equations (5), (7) and (8) into (9) and excluding body forces, we obtain: where K e and F e are the element stiffness matrix and the element nodal force vector. They can be computed numerically by a Gauss-Legendre integration with 3 × 3 grid. If considering arbitrary δd, we can obtain the standard finite element formulation as: At the end, the unknownd can be obtained by solving this equation. With this solution, in the post-processing stresses and strains are firstly computed at Gauss points for each element and then an extrapolation technique is used to get stresses and strains at element nodal points.
The above implementations are achieved by an updated MATLAB code based on codes of a classical 2D Cauchy continuum as presented in [39].

Numerical Simulation
In this section, we intend to numerically investigate the effects of direction and scale of microstructures for a composite material that can be considered as an assembly made of rigid rectangular blocks in contact with elastic interfaces (Figure 1), where each rectangular block has the width of b and height of h. The assembly is arranged as an interlocking structure and the interfaces of blocks form the microstructures of this composite material. A homogenization procedure presented in [30] can be used to describe the assembly as an equivalent Cosserat continuum. In this work, 7-block representative volume element (RVE) that is highlighted with orange color in Figure 1 is used for the homogenization procedure to produce the Cosserat constitutive matrix in Equation (4). In the highlighted RVE, the blocks' centroids are represented by green crosses, and red lines mean the outward unit normal vectors of the central block's interfaces. The direction of microstructures can be changed by transforming the assembly of an angle β from x − y coordinate system to X − Y coordinate system as shown in Figure 1. In this study, we select 6 values of β (0 • , 30 • , 60 • , 90 • , 120 • , 150 • ). Furthermore, to obtain various scales of microstructures, 4 different block sizes are used by fixing the height of block h = 0.1 m and changing the aspect ratio ρ = b/h = 1.5, 3, 7, and 15, where a greater ρ corresponds to a longer rectangular block, as a consequence, showing more orthotropic nature. Figure 1. Schemes of the considered assembly and RVE.
The blocks interact among themselves through elastic common interfaces. The adopted spring stiffness at the interfaces is: where k n and k t are the normal and tangential stiffness per unit length, respectively. Here we have k n = 576.58 MPa/m and k t = 288.29 MPa/m. The rotation stiffness of interface is computed k r = k n (d/2) 2 , where d is the length of interface. Therefore, the Cosserat constitutive matrices of the reference RVE when β = 0 • can be obtained by the homogenization technique that is based on an equivalence energy criterion between the material's discrete system of and the continuum model [30]. The constitutive matrix of transformed assembly can be obtained as follows: where C 0 is the constitutive matrix when β = 0 • , Q is the usual transformation matrix [40].
For comparison, the Cauchy continuum is also considered here to carry out the same simulations as performed by the Cosserat continuum. Because of the lack in microstructures, the constitutive relation of the Cauchy continuum has the form of: σ =Âε, where the matrixÂ is obtained from A and their relationship can be found in the previous literature as [41]: Tables 1-4 list the constitutive components of all configurations, where components keeping zero for all configurations are not reported. It can be seen that there are more zero components when β = 0 • and 90 • . Actually, only diagonal components of the constitutive matrix exist and an orthotropic nature of the material is observed for these two transformation angles. For other angles, all components of matrices A and D appear. Consequently, materials with these RVEs show a centrosymmetric nature [30]. B = 0 for all configurations, meaning there is no coupling between stresses/microcurvatures and microcouples/strains. As the aspect ratio ρ increases, constitutive components change monotonously except for A 2222 , A 1212 of Cosserat continuum andÂ 2222 of Cauchy continuum when β = 0 • . Because of the fixed height of blocks (h), these components stay the same with increasing ρ.  In the following, a classical tension problem of a square plate with a circular hole is studied for all above-mentioned configurations. Figure 2 shows the sketch of the problem and its finite element meshing. The plate has a side length of 10 m and the radius of the hole is 1.25 m. A total of 1440 elements is used for this model. Due to the singular nature resulting from the presence of the hole, stress concentration is more likely observed around the hole under tension force. To make sure the results are reliable, a finer mesh is applied near the hole. Uniform tensile stress σ 0 = 1 MPa is applied on the right side of the plate. The left side of the plate is fixed symmetrically in the x-direction and the bottom left point is additionally fixed in the y-direction. In the present paper, the direction and scale effect of microstructure is investigated by setting various β and ρ. The simulation results are shown below.

Results
The results of displacements and stresses for the Cosserat and Cauchy models and relative rotation for the Cosserat model are presented in this section to show the effects of microstructure's direction and scale on the behavior of microstructured composite materials. Figure 3 depicts the horizontal displacement results u 1 of Cosserat and Cauchy models. It can be seen that change in direction of microstructure has a significant effect on u 1 for both two models. The smallest u 1 can be observed when β = 0 • . As β changes from 0 • to 150 • , u 1 increases to the greatest at β = 90 • and then decreases but this is expected for the shortest block case ρ = 1.5. In the case of ρ = 1.5, the greatest u 1 occurs at β = 60 • and 120 • but that is also close to u 1 at β = 90 • . For all β representing the directions of microstructure, u 1 has a reduction as ρ increases. For the plate with the shortest blocks, displacement localization can be observed at the middle of the right side of the plate. However, as ρ increases, displacement localization reduces and uniform displacement distribution can be observed at the right side of the plate.
It can be seen from Figure 3 that the difference in u 1 between Cosserat and Cauchy models is not obvious for the orthotropic materials (β = 0 • and 90 • ). Here we take the difference in the maximum u 1 between the two models as ∆u 1 . For various ρ, ∆u 1 is 3.02-8.63 mm when β = 0 • and 2.13-5.85 mm when β = 90 • . However, the difference is more evident for the centrosymmetric materials. ∆u 1 increases from 20 mm to 69 mm as ρ increases when β = 30 • , 60 • , 120 • , and 150 • .   It should be noted that there is a big difference in u 2 between Cosserat and Cauchy models. u 2 of these two models has a similar distribution only when ρ = 1.5. As ρ increases, u 2 distributions of two models become different. When β = 30 • and 150 • , the Cosserat continuum shows high-intensity u 2 distribution on the right plate as ρ increases, e.g., the maximum u 2 is up to 117 mm when β = 150 • ; however, u 2 from the Cauchy continuum has no clear high-intensity distribution with ρ and the maximum u 2 is just 77 mm when β = 150 • . An opposite difference between two models can be observed when β = 60 • and 120 • , that is high-intensity u 2 distribution can be observed by the Cauchy but not the Cosserat model.  The horizontal stresses σ 11 of the two models are shown in Figure 5. The results show directionality of distribution of σ 11 with respect to the direction β. For orthotropic materials, the high-stress area is parallel to the x-direction and the peak stress happens at the top and bottom points of the hole edge (i.e., x = 0). However, for centrosymmetric materials, the high-stress area is inclined with the x-direction by an angle that can be related to the direction of microstructure, and the location of peak stress changes. It can be seen the directionality of σ 11 is more evident for greater ρ. As the increase of ρ, the high-stress area becomes wider. The exception happens at β = 90 • which shows the opposite behavior. The existence of the hole as a singularity can result in the concentration of stress around the hole edge. In the following, the stress concentration at the hole boundary will be discussed.
Vertical stress σ 22 (Figure 6) also shows directionality of stress distribution with respect to β. There can be seen a difference in σ 22 between two models for various ρ. For relative short blocks (ρ = 1.5 and 3), two models produce closed behavior of σ 22 . However, for longer blocks (ρ = 7 and 15) especially at β = 60 • , 90 • and 120 • , σ 22 from Cauchy model is significantly greater than that from Cosserat model. The relative rotation, defined as the difference between macrorotation (θ) and microrotation (ω), is a peculiar measurement in the Cosserat model. Figure 7 shows the relative rotation θ − ω for the Cosserat model. The directionality of θ − ω distribution can be also observed and that is more evident as ρ increases. For orthotropic materials, θ − ω shows point symmetric behavior with respect to the hole center where positive and negative θ − ω can be both observed, but the value of θ − ω is close to 0. Thus, the Cosserat model is very close to the Cauchy one. For β = 30 • and 60 • , the plate domain mainly undergoes a positive relative rotation. Oppositely, for β = 120 • and 150 • , this domain mainly undergoes a negative θ − ω. The above-mentioned indicates that relative rotation can be affected by the direction of microstructure. There is less relative rotation when the microstructures are arranged along parallel and perpendicular to the direction of force, whereas the microstructures arranged along other directions would result in an obvious relative rotation acting at a certain orientation.

Discussions
The problem of stress concentration has always been focused on due to the presence of singularity [32,42,43]. Under horizontal tension in this study, it can be seen that the stresses σ 11 and σ 22 in the plate (Figures 5 and 6) is mainly concentrated around the boundary of the hole. To better show the stress distribution, in the following, by transforming the stress state from Cartesian coordinate to polar coordinate, the hoop stress σ h at the hole boundary is depicted for two models in the polar coordinate system as shown in Figure 8. In this way, σ h represents σ 22 when the polar angle equals to 0 or π, whereas the hoop stress denotes σ 11 when the polar angle is π/2 or 3π/2. Therefore, the location and magnitude of σ h can be clearly observed. For both Cosserat and Cauchy models, it can be seen that the distribution of σ h at the hole boundary is point symmetric to the center of hole and also shows directionality that depends on the aspect ration ρ and angle β. For the orthotropic materials, the distribution of σ h is symmetrical along the vertical direction (π/2 − 3π/2) in the polar coordinate. The highest σ h can be observed at the top and bottom points of the hole boundary (i.e., polar angle equals to π/2 and 3π/2) for the Cosserat model and Cauchy model when β = 0 • , indicating that the peak stress results from the horizontal stress σ 11 . As for the Cauchy model when β = 90 • , it is consistent with the above results for small ρ. As ρ increases, the highest σ h trend to be occurred at the right and left points of the hole boundary (i.e., polar angle equals to 0 and π). Since we observed a significant σ 22 concentration in the Cauchy model in Figure 6b, the vertical stress σ 22 is able to result in the peak stress for these cases. For centrosymmetric materials, the distribution of σ h is no longer symmetrical along the vertical direction but deviates from it to more or less of an extent because of the transformation angle of rectangular blocks. Thus, the highest σ h does not occur at these special points, i.e., polar angle equals to 0, π/2, π, or 3π/2. This is consistent with the results by early study [9], which investigates the stress distribution of the layered surrounding rock tunnel by considering different angles of rock joints (microstructures). It is also shown that the stress distribution is symmetrical when angle equals to 0 • and 90 • . When the angle is 45 • the stress presents an asymmetric distribution and the tunnel even comes into being eccentric-pressed.
It should be noted that the effect of β is less for smaller ρ. When ρ = 1.5, the distribution of σ h is close to each other, that is, the highest σ h is located near to the top and bottom point of hole boundary and its value close to 3, whereas the lowest σ h near to right and left points and its value is around 1. Such a result is close to the well-known analytical solution for an infinite isotropic plate with a circular hole [33,44]. However, as ρ increases, the directionality of σ h distribution becomes more evident and the extreme values of σ h also vary. For example, when β = 0 • , as increase in ρ the highest value of σ h can increase to 9. In addition, the difference between the Cosserat and Cauchy models gets bigger with increasing ρ. The smallest values of σ h are not lower than −3 for the Cosserat continuum, whereas lower σ h to −9 can be obtained from the Cauchy continuum, especially when β = 60 • , 90 • , and 120 • .
By investigating measurements of displacements, stress, etc, and the stress distribution of the square plate with circular hole under horizontal tension stress, a direction effect of microstructures can be found for the present problem. With different transformation angle β, the plate can produce mechanical behavior showing obvious directionality. Such a direction effect of microstructures can be also found in previous studies on surrounding rock roadway and tunnels with different dip angles of surroundings [9,45], where the joints of surrounding rocks can represent the microstructure interfaces. In this present paper, similar behaviors are often observed when β = 0 • and 90 • . That is because the microstructures in these two cases are both parallel and perpendicular to the x-direction, showing an orthotropic nature. The difference between these two cases is actually due to the different values of ρ (ρ < 1 when β = 90 • whereas ρ > 1 when β = 0 • ); therefore, measurements of these two cases have similar behavior but different intensities under the horizontal tension stress. As for microstructures not parallel or perpendicular to the x-direction, i.e., β = 30 • , 60 • , 120 • and 150 • , it can be seen from the constitutive matrices (Tables 1-4) that more coupling in the constitutive components such as dilatancy components [46] appear for these cases, showing a centrosymmetric nature. The centrosymmetric material can show different behavior from the orthotropic material by coupling different stresses and strains as well as the couple stresses and curvatures. As the transformed constitutive matrix is related to β (Equation (13)), the behavior of the centrosymmetric material depends on β. Under horizontal tension stress, it is shown that u 1 first increases and then decreases with β. Oppositely, as β increases the maximum σ h first decreases and then increases. Such an effect of β can be compared with the previous study [8] where a similar effect was found, that is, the strength of phyllite decreases first and then increases with the increase of bedding angle.
The direction effect of microstructures is more evident for higher ρ. When ρ is small (e.g., ρ = 1.5), the width of the rectangular block b is close to its height h. Thus, the assembly made of such blocks can show a nearly orthotetragonal behavior (close to isotropic) that is less sensitive to change in the microstructure direction. This could also account for the small differences between Cosserat and Cauchy results when ρ is small since it was known that orthotetragonal materials are very close to Cauchy continua [30,47].
With the increase in ρ, the assembly becomes more orthotropic and the measurements can show more obvious directionality with respect to β. As the length of microstructure (i.e., internal length) is introduced to the Cosserat continuum, asymmetries are generated for the shear stress and shear strain fields. Therefore, in the Cosserat model each stress is coupled with asymmetric shear strains through two constitutive components (A ij12 and A ij21 ), and vice-versa. In the Cauchy model, there is just one component (Â ij12 ). When ρ is small, the difference between A ij12 and A ij21 is not evident. As ρ increases, such difference increases rapidly, showing a higher degree of asymmetry. However,Â ij12 of the Cauchy model is an arithmetic mean by A ij12 and A ij21 (Equation (14)), which cannot show the asymmetric behavior of continua. Thus, the difference between the results from Cosserat and Cauchy models may be induced, especially for higher ρ. Moreover, since the microstructures are considered in the Cosserat continuum, the additional sub-matrix D is involved in the Cosserat constitutive relation rather than the Cauchy one. The components of D are negligible for small ρ but become prominent as ρ increases. It has been found that the stress can be re-distributed within the Cosserat continuum [34]. For the smaller ρ, the Cosserat continuum behaves close to the Cauchy continuum as mentioned above. As a result, the re-distribution of stress can be neglected. However, as the scale of the microstructure increases, such a re-distribution can be more prominent for higher ρ. This can be used to account for the difference in σ 22 between the Cosserat and Cauchy continuum, especially for large ρ.

Conclusions
The present paper investigates the mechanical behavior of the microstructured composite treated as Cosserat continuum by considering various microstructure's directions (β) and scales (ρ). According to the constitutive parameters obtained from the Cosserat homogenization procedure, the composite studied here can be classified as orthotropic and centrosymmetric materials depending on the direction β. The simulations are conducted for a tension problem of a microstructured plate with a circular hole, so this paper also focuses on the stress distribution around the hole. The main conclusions are as follows: (1) The mechanical behavior of microstructured composite changes as the microstructure's directions β, thereby showing a directionality of measurement distribution such as stresses. In general, orthotropic materials show similar behaviors but with different intensities, and the behavior of centrosymmetric is related to various β. (2) The increasing microstructure's scale ρ can results in more evident effect of β and difference between the Cosserat and Cauchy models. Such an effect of ρ is clearer for the centrosymmetric materials than orthotropic materials. (3) The Cosserat continuum is able to better describe the direction effect of microstructures due to the relative rotation that not only shows the directionality of distribution but also varies with the microstructure direction. The Cauchy continuum does not have such advantages because there is no relative rotation and tangential strains are symmetric. (4) The extreme value and its location of the hoop stress σ h around the hole depend on β.
For smaller ρ, the highest and smallest σ h are close to 3 and −1, which is similar with the classical result of the isotropic material. As ρ increases, a highest σ h up to 9 can be observed. (5) Difference in the hoop stress σ h between the Cosserat and Cauchy model is mainly in the smallest σ h , especially for greater ρ when β = 60 • , 90 • , and 120 • . All the smallest σ h of the Cosserat model are greater than -3, whereas the Cauchy model can have a σ h as low as −9.
From this present study, the effect of the microstructure's directions on mechanical behavior of microstructured composite can be found, especially for the composite with large scale of the microstructure. The area applying the development of this research can be for microstructured materials with various dimensions (from micromaterials to macromaterials), where the scale of the microstructure should be comparable to the material's dimension. For example, the layered rock with joints has different dip angle due to the geological formation. The stress distribution can be also used for such composite with singularity, not only the circular shape used here, but alternative shapes. The rectangular microstructures with standard interlocking give a basic and important research aspect of the effect of the microstructure's direction and scale. As more and more advanced materials are developed nowadays, various microstructures formations are interesting to be studied in future research.

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