Next Article in Journal
A Trust Region Reduced Basis Pascoletti-Serafini Algorithm for Multi-Objective PDE-Constrained Parameter Optimization
Next Article in Special Issue
Microwave Characterization and Modelling of PA6/GNPs Composites
Previous Article in Journal
Multivariable Panel Data Cluster Analysis of Meteorological Stations in Thailand for ENSO Phenomenon
Previous Article in Special Issue
Theoretical and Computational Results of a Memory-Type Swelling Porous-Elastic System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interval-Based Computation of the Uncertainty in the Mechanical Properties and the Failure Analysis of Unidirectional Composite Materials

by
Dimitris G. Sotiropoulos
1,2 and
Konstantinos Tserpes
3,*
1
Department of Electrical and Computer Engineering, University of Peloponnese, 26334 Patras, Greece
2
Laboratory of Applied Mathematics, School of Science and Technology, Hellenic Open University, 26335 Patras, Greece
3
Laboratory of Technology & Strength of Materials, Department of Mechanical Engineering & Aeronautics, University of Patras, 26504 Patras, Greece
*
Author to whom correspondence should be addressed.
Math. Comput. Appl. 2022, 27(3), 38; https://doi.org/10.3390/mca27030038
Submission received: 28 March 2022 / Revised: 26 April 2022 / Accepted: 27 April 2022 / Published: 29 April 2022

Abstract

:
An interval-based method is presented to evaluate the uncertainty in the computed mechanical properties and the failure assessment of composite unidirectional (UD) laminates. The method was applied to two composite laminates: a carbon/epoxy and a glass/epoxy. The mechanical properties of the UD lamina were derived using simplified micromechanical equations. An uncertainty level of ±5% was assumed for the input properties of the constituents. The global minimum and maximum values of the properties were computed using an interval branch-and-bound algorithm. Interval arithmetic operations were used to evaluate the uncertainty in the Hashin-type failure criteria in a closed form. Using the closed-form uncertainties of intervals and sets of stresses obtained by finite element analysis, the uncertainty in the failure assessment was quantified for the two composite laminates. For the assumed uncertainty level of ±5%, the computed uncertainty for the mechanical properties ranges from 6.64% to 10.63% for the carbon/epoxy material and from 6.72% to 12.28% for the glass/epoxy material. For evaluating the uncertainty effect on the efficiency of failure criteria, a probability of failure function, which employs interval boundaries, was defined and proved capable of evaluating the whole spectrum of stresses.

1. Introduction

Fiber-reinforced composites, especially carbon-fiber-reinforced plastics (CFRPs), have become the primary structural material of lightweight structures [1]. At the same time, the numerical design of composite structures finds an increasing use in aiming to reduce the time and cost of the development and design phases.
Micromechanics of composites is the branch of mechanics of composites that predicts the elastic properties and strengths of a unidirectional lamina by using the properties of the fibers and the matrix. Due to the anisotropy and inhomogeneity of composites, the complete experimental characterization of their mechanical properties requires many tests that must be conducted by using different lay ups and loading conditions. There have been proposed several sets of micromechanical analytical relations in the literature. The most widely used set is the one proposed by Chamis [2]. For the failure analysis of UD composites, the several sets of failure criteria are mainly due to the diversity of composite materials [3]. Each set is suitable for a specific type of composite materials and specific loading conditions. A widely used set is the Hashin-type failure criteria [4]. They are of a polynomial type, are stress based, and can distinguish between the different failure modes. Moreover, they have high accuracy and are highly compatible with the finite element (FE) method. Due to the above characteristics, they have been extensively used in progressive-damage models [5,6,7,8].
The computational efficiency of micromechanical equations and the failure criteria depends on the experimental error. The scatter of the results is higher for composites than for monolithic materials due to the material’s complexity, the involvement of different types of tests, and the presence of manufacturing defects. It is, therefore, of great importance to quantify the uncertainty of the predictions of micromechanical equations because it reflects the computation of stresses and strains and the predicted damage, thus compromising the performance of composite components, leading to the use of more conservative designs that do not fully exploit the performance and environmental opportunities offered by composites. Mainly, manufacturing uncertainties can be summarized as follows: (a) fibre architecture variations that are usually generated during production, handling, or storage of prepregs, dry textiles, and per-forms, (b) matrix material uncertainties caused by variations in storage conditions or uncertainties in resin composition and formulation, and (c) variations in environmental parameters and process conditions.
The majority of the works reported in the literature for quantifying uncertainty in the mechanical properties of UD composites are based on probabilistic approaches [9]. In addition, Monte-Carlo methods have also been used based on assumed probability distributions [10]. However, Monte-Carlo methods are computationally demanding and unsuitable for complex material systems and parts of a larger scale. On the other hand, interval analysis has been used in multidisciplinary engineering applications such as in the formulation of fuzzy finite elements [11], the time-dependent reliability analysis of fatigue crack growth [12], and structural optimization [13].
In recent works, Alazwari and Rao [14,15] performed uncertainty quantification of the micromechanical properties of composite materials. Moreover, in [16], two different versions of the interval-based uncertainty failure model, based on the maximum stress failure theory and the Tsai–Wu failure theory, using the universal gray system theory and the truncation-based interval analysis, were presented. These are the only published works on the interval approach to quantify uncertainty in the micromechanical properties of UD composites; however, they do not include the quantification of uncertainty in the failure analysis.
Although there are numerous numerical methods for approximating the range of real functions, none of them guarantee the required accuracy of the result, and the majority of them (e.g., Monte Carlo) face computational-time constraints. However, the main advantage of interval arithmetic is inclusion isotonicity, a valuable property in connection with the fundamental theorem [17,18] of interval analysis that guarantees the production of reliable numerical solutions with mathematical correctness. Since the uncertainty problem is equivalent to bounding the range of a multivariate function over a multidimensional rectangle (box), interval mathematics is the most straightforward technique for modeling uncertainty and providing guaranteed enclosures. When the interval parameters of a multivariate function are large, it can be challenging to determine its exact range. This paper proposes an efficient interval algorithm with mathematical correctness to solve this problem.
The article is organized into five sections. In Section 2, we briefly introduce the basics of interval arithmetic, the mean value form, and the corresponding optimal center. In Section 3, we present our computational method based on a rigorous interval branch-and-bound algorithm for computing the range of a multivariate function with large-range uncertain parameters. Section 4 presents the computed uncertainties in the elastic properties and strengths of two UD composite laminae and the sensitivity analysis of the elastic properties. Section 5 describes the methodology and results of the uncertainty in the failure assessment of the composite laminae, performed using the Hashin-type failure criteria. The paper is completed with three Appendices in which we exhibit the Chamis micromechanical equations, the Hashin-type failure criteria, and equations for the computation of uncertainty intervals of the criteria.

2. Interval Analysis

This section presents the interval mathematics required for the rest of the paper. We provide only a summary of basic interval definitions and properties; for a more in-depth discussion, see, for example, [17,18,19]. Interval arithmetic [17] generalizes ordinary arithmetic to closed intervals of the real line. Let I ( ) = [ a , b ] | a b , a , b be the set of real closed intervals and I ( n ) the n -dimensional interval vectors (also called boxes). Given a pair of intervals X = [ x _ , x ¯ ] and Y = [ y _ , y ¯ ] with 0 Y , the four basic interval arithmetic operations are defined by
X + Y = [ x _ + y _ , x ¯ + y ¯ ] , X Y = [ x _ y ¯ , x ¯ y _ ] , X × Y = [ min ( x _ y _ , x _ y ¯ , x ¯ y _ , x ¯ y ¯ ) , max ( x _ y _ , x _ y ¯ , x ¯ y _ , x ¯ y ¯ ) ] , X / Y = X [ 1 / y ¯ , 1 / y _ ] , if 0 Y .
It is worth noting that the above elementary operations are inclusion isotonic, which means: if X Y and X Y then X Y X Y for each binary operation + , , × , / . In practice, rounded interval arithmetic is utilized to bound roundoff error instead of the exact form above. Rounded intervals may be somewhat larger, but they always contain the exact result [20]. The same interval X I ( ) can be defined in three ways: either by its endpoints, or by its midpoint m ( X ) = ( x _ + x ¯ ) / 2 and radius r ( X ) = ( x ¯ x _ ) / 2 as X = m X r X , m X + r X , and, lastly, in terms of its midpoint and its uncertainty level λ as X = m X [ 1 λ , 1 + λ ] , where λ = r ( X ) / | m ( X ) | . Sometimes, for clarity, the lower bound of X is also denoted as inf X , and the upper bound is denoted sup X .
Let f r g ( X ) denote the range of the function f : n over an interval vector X I ( n ) , i.e., f r g ( X ) = min x X f ( x ) , max x X f ( x ) . A function F : I ( n ) I ( ) is called an inclusion function of f if f r g ( X ) F ( X ) for any X I ( n ) . It is well-known that any intersection of inclusion functions is an inclusion function. Natural extension, mean value forms, and Taylor expansion create inclusion functions [17] that are inclusion isotones, i.e., X , Y I ( n ) , X Y F ( X ) F ( Y ) . Each of these forms has slightly different properties and convergence orders. For a thorough discussion on these issues, see [19].

Optimal Mean Value Form

Assuming that f is continuously differentiable, the gradient’s inclusion function f may be computed (for almost any function specified by a finite procedure) employing automatic differentiation [21] and interval arithmetic. The most important centered form is the mean value form (MVF), which is obtained by taking the first order of Taylor’s expression:
F M V F X , c = f c + X c T F X = f ( c ) + i = 1 n G i X X i c i
where X is an n -dimensional vector of intervals X i = [ x _ i , x ¯ i ] , i = 1 , , n c = m ( X ) = c 1 , , c n is the midpoint of the box X , and G i X = [ g _ i , g ¯ i ] is the interval extension of the derivative f / x i . A second-order approximation to the range of f is proved in [22,23] and the isotonicity property when the center c is the midpoint. The mean value form F M V F ensures tighter enclosures than the natural interval extension F when the interval X is narrow. Nevertheless, on the other hand, when the diameter of X is large, F M V F may drastically overestimate the range of f . Consequently, it is common to use the intersection F M V F X , c F ( X ) to obtain a more accurate range estimate.
In [24], Baumann introduced the notion of optimal mean value forms and proved a weaker property called one-sided isotonicity for the produced forms, a valuable property that in the framework of global minimization algorithms suffices. The point c , which yields the greatest lower bound, and the point c + , which yields the lowest upper bound among all centers c X , are called optimal centers. Baumann provided asymmetric formulas for these centers; while, in [25], the authors determined in a straightforward manner equivalent symmetric expressions for the optimal centers, and they proved the conjecture that the width of (1) is minimal for all centers between c and c + . Furthermore, exact formulas for the lower and upper bounds of optimal MVFs are provided, which show that there is no need for the box’s width to tend to zero to have sharp bounds. The lower bound of the MVF attains its maximum at the center c = c 1 , , c n T n with components [25]:
c i = x ¯ i , i f g ¯ i 0 , x _ i , i f g _ i 0 , m ( X i ) r ( X i ) g ¯ i + g _ i g ¯ i g _ i , o t h e r w i s e .
In order to determine the optimal center c + for the best upper bound of MVF, the optimal point c should be mirrored about the midpoint m ( X ) , resulting in c + = m ( X ) + [ m ( X ) c ] , since the optimal centers c and c + are always symmetric concerning the midpoint (see [25]).
For simplicity of notation, the left endpoint of the interval F M V F ( X , c ) will be denoted by inf F M V F ( X , c ) . Let I 0 be the subset of the index set I = { 1 , , n } , where i I 0 if 0 G i ( X ) = [ g _ i , g ¯ i ] , then it can be proved [26] that the lower bound of the MVF (1) with center c is given by:
inf F M V F ( X , c ) = f ( c ) + i I 0 w ( X i ) g ¯ i g _ i g ¯ i g _ i
where w ( X i ) = x ¯ i x _ i denotes the width of the interval X i [25]. It is evident that when the function f is strictly monotone over an interval vector X , then I 0 = and, therefore, inf F M V F ( X , c ) is the exact lower bound; in a different case, inf F M V F ( X , c ) is overestimated, and we must determine a coordinate direction to bisect the box X .
The adopted selection rule for subdividing the box is inspired by Equation (3), and our goal is to maximize the lower bound of the inclusion function after each subdivision. We observe that, for all i I 0 , the summation terms are all negative; therefore, to maximize the lower bound of MVF, the optimal component k I 0 is determined by the largest absolute summation term, which subtracts the most. This selection rule is very similar to Ratz’s rule C [26,27], where the goal was to minimize the width of the inclusion in Equation (2).
Our algorithm, presented in Section 3, uses first-order information through an interval gradient evaluation and takes advantage of the valuable properties of the optimal center (2) of the MVF for both the bounding and the branching process.

3. The Proposed Branch-and-Bound Algorithm

This section introduces an interval branch-and-bound algorithm for computing the range
F = min x X f ( x ) , max x X f ( x )
of a continuously differentiable function f : n over a box X to solve uncertainty problems with large-range interval parameters. As the maximization problem can be transformed into a minimization problem, precisely, max f ( x ) min ( f ( x ) ) , the proposed algorithm aims to compute verified enclosures for the global minimum of functions.
Algorithm 1 formulates the main procedure in pseudo-code. In Step 1, we first find a guaranteed lower bound, F min , on the range of values of f over the search box X 0 , and then apply the same procedure in Step 2, to find a lower bound on ( f ) over X 0 to obtain the upper bound, F max . Finally, in Step 3, we determine the enclosure F of the range from the previous two steps.
Algorithm 1. ComputeRange f , X 0 , ε , F
Input: objective function f , starting box X 0 , and tolerance ε
Output: enclosure F of the range
1: GlobalMinimize f , X 0 , ε , F min ;
2: GlobalMinimize f , X 0 , ε , F max ;
3: F = inf F min , sup F max ;
4: return F ;
Algorithm 2 describes the branch-and-bound interval algorithm we implemented, employing two well-known accelerating devices: the cut-off test and the monotonicity test [21,26]. The algorithm uses first-order information of the objective function utilizing an interval gradient evaluation, firstly, to check monotonicity and then apply the optimal mean value form for bounding the range of the function. For bounding the objective function, we take the intersection of the natural interval extension F ( X ) with optimal mean value form F M V F ( X , c ) to obtain the best lower bound for the range of f . Moreover, as discussed in the previous, Section Optimal Mean Value Form, the calculated derivative bounds are exploited to determine an optimal component k (subdivision direction) to bisect the box during the branching phase of the algorithm.
In general, the search tree is incrementally expanded by iterating the following steps: (i) the initial search box X 0 is subdivided into smaller sub-boxes, (ii) the objective function and its derivatives are bounded over the sub-boxes, and (iii) sub-boxes that cannot contain a global minimizer with certainty are removed. We next summarize the basic interval techniques that accelerate the search process:
  • Function range test: A box Y is discarded from further consideration when the lower bound inf F Y is greater than the current upper bound f max . When the range test fails to remove it, it is stored in the working list with candidate sub-boxes for further investigation.
  • Cut-off test: The function range test is applied for all candidate sub-boxes in the working list when f max is improved. Of course, the greater the improvement in f max is, the more influential the cut-off test is.
  • Monotonicity test: Determines whether the objective function f is strictly monotone in an entire sub-box Y X 0 or at least one coordinate direction, in which case Y cannot contain a global minimizer. Therefore, the whole sub-box is discard or its dimension reduced as much as possible when Y X 0 .
We now give a detailed algorithmic description of our proposed method (Algorithm 2), implemented in C++ using the C-XSC-2.5 library (http://www.math.uni-wuppertal.de/~xsc/ (accessed on 25 March 2022).
Algorithm 2. GlobalMinimize f , X 0 , ε , F
Input: objective function f , starting box X 0 , and tolerance ε
Output: enclosure F for the global minimum value
01: Y = X 0 ; L = { } ;   L = { } ;                /* initialize working and result list */
02: MonotonicityTest X 0 , Y , F Y , d e l ;            /* apply monotonicity test */
03: c = OptimalCenter Y , F ( Y ) ;        /* compute optimal center using Equation (2) */
04: f max = sup F ( c ) ;                /* initialize upper bound for global minimum */
05: F Y = F ( c ) + F ( Y ) ( Y c ) F ( Y ) ;             /* optimal mean value form */
06: k = OptimalComp Y , F ( Y ) ;        /* optimal component for subdividing the box */
07: if  k = 0  then  L = L Y , inf F Y , k ;   /* no subdivision required, append to the result list */
08: else  L = L Y , inf F Y , k ;                    /* append to working list */
09: while  L { }  do
10:   Y , inf F Y , k = PopHead L ;         /* remove a triple from the head of the working list */
11:  if  k 0  then Bisect Y , k , U 1 , U 2 ;             /* subdivide Y into two subboxes */
12:  for  i = 1 to NBoxes do                    /* number of boxes is 1 if k = 0 */
13: if  inf F U i > f max  then next i ;       /* box U i is discarded due to function range test */
14:  MonotonicityTest X 0 , U i , F U i , d e l ;
15: if  d e l = = t r u e  then next i ;          /* box U i is discarded due to monotonicity test */
16:  c = OptimalCenter U i , F U i ;      /* compute optimal center using Equation (2) */
17: if  sup F ( c ) < f max then                /* update upper bound f max if possible */
18:   f max = sup F ( c ) ; CutOffTest L , f max ;            /* and apply cut-off test */
19:  F U = F ( c ) + F U i ( U i c ) F ( U i ) ;      /* optimal mean value form for box U i */
20: if inf F U > f max then next i ;          /* box U i is discarded due to function range test */
21:  k = OptimalComp U i , F U i ;          /* optimal component for subdividing U i */
22: if w r e l [ inf F U , f max ] ε or  k = 0  then  L = L U i , inf F U , k ;
23: else  L = L U i , inf F U , k ;
24: Y , inf F Y , k = Head L ;            /* set Y to the first element of the result list */
25: F = inf F Y , f max ;               /* construct the global minimum enclosure. */
26: return F ;
Initially, the working list L and the result list L are empty. Following the evaluation of F ( Y ) , the monotonicity test is called to possibly reduce the dimension of the original box, and then compute the optimal center c using Equation (2). The upper bound f max is initialized by evaluating F ( c ) using interval arithmetic to bound all rounding errors. In Step 5, we bound the range of f over Y using the intersection F M V F Y , c F ( Y ) . While, in Step 6, the optimal component k is determined by finding the largest absolute summation term of the expression in Equation (4). The value k = 0 indicates that f is strictly monotone over the entire box Y or, in other words, the subset I 0 = . Then inf F M V F ( Y , c ) is the exact lower bound, and, therefore, no subdivision is required and the box is stored in the result list L . In Step 8, the box is stored in the working list L , in a different case. It is essential to notice that, when k = 0 , the total cost of the algorithm is two function evaluations and one gradient evaluation, the working list remains empty, and the algorithm jumps to Step 24.
The search is carried out by prioritizing the sub-boxes Y using the best-first strategy [21,26]. The rationale for this strategy is to concentrate on the most promising box; that is, the one with the lowest lower bound because we are looking for the global minimum. This task can be accomplished quickly by using and maintaining a priority queue or list sorted in a non-decreasing order relative to the lower bounds inf F Y .
Steps 10–23 are applied to each candidate sub-box Y until the working list L is empty. The algorithm removes the first element from L and subdivides the box Y in the coordinate direction k . Special care is taken to handle unnecessary subdivisions if it has previously been proved that f is strictly monotone. For each sub-box (one or two), first, the function range test (Step 13) is applied, and then the monotonicity test algorithm as proposed by Ratz [26]. When the upper bound f max is improved, Steps 17–18, the cut-off test is applied to the working list L to discard sub-boxes that definitely do not contain the global minimum. After bounding the range of f over U i , the range test is performed in Step 20 and U i is either added to the result list L or the working list L , depending on the termination criterion (Steps 22–23). A box is stored in the result list when the relative diameter [21] of F is smaller than a prescribed accuracy ε (defined by the user). When no candidate sub-boxes are contained in the working list, the algorithm terminates by returning an enclosure for the global minimum F .
Several examples have been investigated to test the algorithm’s correctness and efficiency. The inclusion functions have been implemented using the C-XSC routines while the gradients are computed with forward automatic differentiation [22]. All computations, including those involving floating-point numbers, were carried out using rounded interval arithmetic. Although we present the results with fewer digits, the bounds of the range were calculated with 14 significant digits accuracy. Therefore, we have set ε = ( 1 / 2 ) × 10 14 .
To summarize, Algorithm 2 makes three uses of first-order information. Firstly, the monotonicity test can eliminate or reduce the current sub-box. Second, determining optimal MVF centers while simultaneously attempting to improve the global minimum’s current upper bound. Thirdly, using optimal MVF to improve the function range’s enclosure and use a more sophisticated subdivision rule when necessary.

4. Uncertainty Computation in the Elastic Properties and Strengths

The algorithm of Section 3 was applied for the uncertainty computation of the elastic properties and strengths of two different composite plies that find extensive use in the aeronautical and automotive industries, namely, the carbon/epoxy AS4/3501-6 composite ply and the glass/epoxy 21xK43-LY556 composite ply.
The elastic properties of the UD lamina considered are:
  • E11 = longitudinal Young’s modulus (fibers direction);
  • E22 = transverse Young’s modulus (normal to fibers);
  • E33 = normal Young’s modulus (normal to lamina);
  • G12 = in-plane shear modulus;
  • G13 = out-of-plane shear modulus;
  • G23 = out-of-plane shear modulus;
  • ν12 = in-plane Poisson’s ratio;
  • ν13 = out-of-plane Poisson’s ratio;
  • ν23 = out-of-plane Poisson’s ratio.
The strengths of the UD lamina considered are:
  • F1t = longitudinal tensile strength (fibers direction);
  • F1c = longitudinal compressive strength (fibers direction);
  • F2t = transverse tensile strength (normal to fibers);
  • F2c = transverse compressive strength (normal to fibers);
  • F12 = in-plane shear strength;
  • F13 = out-of-plane shear strength;
  • F23 = out-of-plane shear strength.
The mechanical properties of the composite plies have been calculated using the Chamis micromechanical equations [2] reported in Appendix A. The nominal values and the corresponding input uncertainty intervals of the elastic properties and strengths of the fibers and the matrix for the two composite plies are listed in Table 1, Table 2, Table 3 and Table 4. Due to a lack of experimental data, a fundamental assumption for mechanical properties has a ±5% uncertainty interval, as a confidence interval cannot be established without knowing the shape of the uncertainty distribution. For this reason, we assumed that all properties had a rectangular probability distribution for numerical results.
The computed uncertainties in the elastic properties of the composite plies in terms of the interval range, the mid/rad, and the uncertainty level are presented in Table 5 and Table 6 for the two composite materials. The uncertainty level ranges from 6.64% to 10.63% for the AS4-3501/6 material and 6.72% to 12.28% for the 21xK43-LY556 material. Uncertainty levels for each property are influenced by the mathematical relationship between the input and output uncertainty. For a given input uncertainty of ±5% for the properties of the constituent materials, the percentage of the output uncertainties of the elastic properties is significantly higher and, in some cases, reaches 12%.
For the purposes of testing the algorithm, we conducted an artificial experiment with large-range interval parameters (uncertainty level 60%). For each elastic property (objective function), we report the function evaluations (FE), the number of gradient evaluations (GE), the number of bisections (NB), and the maximum list length (LL).
According to Table 7, the most difficult problem seems to be v23, which is modeled as an objective function with seven (7) interval parameters (see Appendix A). While the most accessible elastic property is E11, where the optimal center determines the exact range, and no further refinement is needed.
The level of computed uncertainty for both composite materials is depicted in Figure 1, when all input parameters varied between 1% and 60% of their nominal values. For an input level of 60%, the output level ranges from 71.2% to 89.6% for the AS4-3501/6 material and 71.8% to 94.7% for the 21xK43-LY556 material. We observe from Figure 1b that the output uncertainty of G12 coincides with G23 since they share the same mathematical expression and Gf12 = Gf23.
The computed uncertainties in the strengths of the composite plies in terms of the interval range, the mid/rad, and the uncertainty level are presented in Table 8 and Table 9 for the two composite materials. Note that subscript “t” refers to tensile and “c” to compressive in the symbols. The uncertainty level ranges from 5.00% to 9.98% for both composite plies. It is observed that, for a given input uncertainty of ±5% for the properties of the constituent materials, the percentage of the output uncertainties in strengths is significantly increased. The variation level of the uncertainty of each property is governed by the mathematical relation between the input and output uncertainties.

Sensitivity Analysis of the Elastic Properties

The uncertainties presented in Table 5 and Table 6 have been derived by considering an uncertainty of 5% for all properties of the constituent materials simultaneously. A piece of beneficial information is how the uncertainty of each input property contributes to the overall uncertainty of the elastic properties and strengths. In order to obtain this information, a sensitivity study of the elastic properties was performed by running analyses in which an uncertainty of 5% was used for a specific input property while the other input properties were kept equal to their nominal value. The results from the sensitivity study are shown in the bar charts of Figure 2 and Figure 3 for the two composite materials.
Although, qualitatively, the contribution is the one expected from the mechanic’s viewpoint, it is very interesting to quantify it. For the AS4/3501-6 material, to the uncertainty of E11, Ef11 contributes by 50.2% and Vf by 49.2%, with Em by the minor 0.6%. To the uncertainty of E22, Vf contributes by 38.6%; Ef22 by 30.1%; and Em by 31.3%. To the uncertainty of G12, Vf contributes by 53.1%; Gf12 by 12.4%; and Gm by 34.5%. To the uncertainty of G23, Vf contributes by 43%; Gf23 by 24.8%; and Gm by 32.2%. To the uncertainty of ν12, Vf contributes by 27%, νf12 by 35.1%; and νm by 37.8%. To the uncertainty of νf23, Vf contributes by 36.4%; νm by 36.4%; νf23 by 16.9%; and νf12, Em, Ef11, and Ef22 by 2.6%. Similar values were obtained for the 21xK43-LY556 material (Figure 3). Similarly, for the strengths, to the uncertainty of F1t, Vf contributes by 29.1%; Fft by 50.2%; and Fmt by 0.7%; and to the uncertainty of F1c, Vm contributes by 49.9% and Vm by 50.1.%. Similar values have been obtained for the 21xK43-LY556 material (Figure 3).

5. Uncertainty Computation in the Failure Assessment

Failure analysis of the UD ply was performed using the Hashin-type failure criteria [4]. Seven different failure modes were considered, namely, matrix tensile and compressive failure, fiber tensile and compressive failure, fiber-matrix debonding, and delamination in tension and compression. The mathematical expressions of the failure criteria are listed in Appendix B.

5.1. Methodology

This section exploits the computed bounds of the properties in connection with the Hashin-type failure criteria [4]. We can guarantee that, although there is uncertainty when the upper bound of a criterion is less than 1, there is no failure (100% safe). In different cases, we can quantify the failure probability. The same procedure can be applied to any other failure criteria, such as the Hoffman, Tsai-Wu, and Tsai-Hill criteria [16].
Here we demonstrate our approach using the matrix tensile failure criterion of Hashin
M T F σ 2 2 F 2 t 2 + σ 12 2 + σ 23 2 F 12 2 = 1 , σ 2 > 0
If we already know the ranges of F 2 t and F 12 , we can express the left hand of Equation (5) using interval arithmetic operations, assuming that 0 F 2 t I = [ F 2 t ¯ , F 2 t ¯ ] and 0 F 12 I = [ F 12 ¯ , F 12 ¯ ] , where the attached superscript “I” denotes the corresponding interval quantity. By taking the natural interval extension of M T F and substituting the corresponding interval numbers F 2 t I and F 12 I , we obtain
M T F I = σ 2 2 F 2 t ¯ 2 + σ 12 2 + σ 13 2 F 12 ¯ 2 , σ 2 2 F 2 t ¯ 2 + σ 12 2 + σ 13 2 F 12 ¯ 2
in the endpoint coordinate system. Therefore, the lower bound is a = inf M T F I and the upper bound b = sup M T F I . In addition, we can reformulate the interval M T F I in the midpoint-radius coordinate system as follows
m M T F I = 1 2 σ 2 2 1 F 2 t ¯ 2 + 1 F 2 t ¯ 2 + σ 12 2 + σ 13 2 1 F 12 ¯ 2 + 1 F 12 ¯ 2 r M T F I = 1 2 σ 2 2 1 F 2 t ¯ 2 1 F 2 t ¯ 2 + σ 12 2 + σ 13 2 1 F 12 ¯ 2 1 F 12 ¯ 2
where r M T F I represents the uncertainty of interval M T F I since M T F I = m M T F I ± r M T F I . After some algebraic manipulations, it is straightforward to derive the following formulae for the uncertainty of M T F I .
r ( M T F I ) = 2 σ 2 2 r ( F 2 t I ) m ( F 2 t I ) m ( F 2 t I ) 2 r ( F 2 t I ) 2 2 + 2 σ 12 2 + σ 13 2 r ( F 12 I ) m ( F 12 I ) m ( F 12 I ) 2 m ( F 12 I ) 2 2 ,
while its midpoint can also be calculated from the expression
m ( M T F I ) = σ 2 2 m ( F 2 t I ) 2 + r ( F 2 t I ) 2 m ( F 2 t I ) 2 r ( F 2 t I ) 2 2 + σ 12 2 + σ 13 2 m ( F 12 I ) 2 + r ( F 12 I ) 2 m ( F 12 I ) 2 r ( F 12 I ) 2 2 .
Hence, we derived the combined uncertainty in closed form explicitly. The expressions of the uncertainty of intervals for all Hashin-type failure criteria are reported in Appendix C.
Let [ a , b ] be the range of values of a criterion under uncertainty; then, we define the probability of failure as:
P ( [ a , b ] 1 ) = 1 , a 1 b 1 b a , a < 1 b 0 , b < 1

5.2. Results

The effect of the uncertainty in the strengths of the UD lamina on the efficiency of the failure analysis has been studied based on a UD composite laminated plate subjected to quasi-static tension. The plate dimensions are 250 mm × 25 mm × 2 mm, and the lay-up is [0/45/90-45]2S. The stress analysis of the plate was performed using the ANSYS FE software and the SOLID185 element. In order to model the tensile load, an axial displacement of 2.5 mm (1% strain) was applied at the nodes of one end while the nodes of the other end were kept fully fixed. The computed stresses for each lamina, used as input to the failure criteria, are listed in Table 10. They were taken at the middle point of the plate. It is noted that, since the input stresses were used only for demonstrating the computational method for the uncertainty of the failure analysis, the FE model was neither optimized nor validated.
The uncertainty in the assessment of the first ply failure was computed. Seven failure modes have been considered, namely:
  • Matrix Tensile Failure (MTF);
  • Matrix Compressive Failure (MCF);
  • Fiber Tensile Failure (FTF);
  • Fiber Compressive Failure (FCF);
  • Fiber-Matrix Shear-Out (FMS);
  • Delamination in Tension (DT);
  • Delamination in Compression (DC).
Given that neither material nor geometrical nonlinearity is present, a linear variation of stresses with the applied load is assumed. Therefore, for the developed stresses to cause failure, a stress multiplication factor was applied instead of load increase.
Figure 4 and Figure 5 plot the variation in the lower and upper bounds and the nominal value of the six failure criteria (Equations (A13), (A14), (A15), (A17), (A18), (A19)) as well as the probability of failure (Equation (6)) of each failure mode with regards to the stress multiplication factor for the two UD composite plies. Note that the fiber compressive failure mode has been excluded since it is not a feasible failure mode for a ply under tension; σ1 is always positive. The very large values of the stress multiplication factor appear for delamination because, for the specific loading (uniaxial tension), delamination is very difficult to develop; it is almost not feasible. As long as the upper bound of the failure criteria remains below 1.0, the probability of failure is 0. When the upper bound exceeds 1.0, the probability of failure starts to increase linearly at a high slope.
Table 11 lists the values plotted in the diagram of Figure 4 for the matrix tensile failure for a range of stress multiplication factors between 4.35 and 5.32. It is shown that, for nominal values above 0.9, for which classical analysis means no failure, there is a probability of failure of more than 30%. For nominal values more than 0.95, the probability of failure exceeds 40%.

6. Conclusions

In the present paper, an interval-based method was presented to compute the uncertainty in the computed mechanical properties and the failure analysis of composite unidirectional (UD) laminates. The method was applied to two composite materials: the carbon/epoxy AS4-3501/6 and the glass/epoxy 21xK43-LY556. An uncertainty level of 5% was assumed for the input properties of the constituents (fibers and matrix).
The computed uncertainty level for the mechanical properties ranges from 6.64% to 10.63% for the AS4-3501/6 material and from 6.72% to 12.28% for the 21xK43-LY556 material. For evaluating the uncertainty effect on the efficiency of failure criteria, a probability of failure function has been defined. The probability is zero until the upper bound of the failure criteria becomes 1.0. From that point, it increases linearly at a high slope. Indicative of the uncertainty in the efficiency of failure criteria is that, for nominal values equal to 0.9, the probability of failure is above 30%, and for values close to 0.990 it is above 50%.
The present work demonstrates the potential of interval arithmetic for problems from materials science and mechanics. The authors intend to perform the uncertain assessment using realistic input uncertainty values coming from tests as a future step.

Author Contributions

Conceptualization, D.G.S. and K.T.; methodology, D.G.S. and K.T.; software, D.G.S.; validation, D.G.S. and K.T.; formal analysis, D.G.S. and K.T.; investigation, D.G.S. and K.T.; resources, D.G.S. and K.T.; data curation, D.G.S. and K.T.; writing—original draft preparation, D.G.S. and K.T.; writing—review and editing, D.G.S. and K.T.; visualization, D.G.S. and K.T; supervision, D.G.S. and K.T.; project administration, D.G.S. and K.T. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Micromechanics Relations

The Chamis [2] micromechanical equations used to compute the mechanical properties of the composite plies are depicted below.
E 11 = V f E f 11 + V m E m
E 22 = E 33 = E m 1 V f 1 E m E f 22
G 12 = G 13 = G m 1 V f 1 G m G f 12
G 23 = G m 1 V f 1 G m G f 23
v 12 = v 13 = V f ν f 12 + V m ν m
ν 23 = V f v f 23 + V m 2 ν m ν 12 E 22 E 11
F 1 t = V f F f t + V m F m t
F 2 t = F 3 t = 1 k σ F m t
F 1 c G m 1 V f
F 2 c = F 3 c = F m c k c
F 23 = F m s k τ
F 12 = F 13 = F m s

Appendix B. The Mathematical Expressions of the Hashin-Type Failure Criteria

M T F = σ 2 2 F 2 t 2 + σ 12 2 F 12 2 + σ 23 2 F 23 2 = 1   ( matrix   tensile   failure   for   σ 2 > 0 )
M C F = σ 2 2 F 2 c 2 + σ 12 2 F 12 2 + σ 23 2 F 23 2 = 1 ( matrix   compressive   failure   for   σ 2 < 0 )
F T F = σ 1 2 F 1 t 2 + σ 12 2 F 12 2 + σ 13 2 F 13 2 = 1 ( fiber   tensile   failure   for   σ 1 > 0 )
F C F = σ 1 F 1 c = 1 ( fiber   compressive   failure   for   σ 1 < 0 )
F M S = σ 1 F 1 c + σ 12 2 + σ 13 2 F 12 = 1 ( fiber-matrix   shear   out   for   σ 1 < 0 )
D T = σ 33 2 F 3 t 2 + σ 13 2 F 13 + σ 23 2 F 23 = 1 ( delamination   in   tension   for   σ 3 > 0 )
D C = σ 33 2 F 3 c 2 + σ 13 2 F 13 + σ 23 2 F 23 = 1 ( delamination   in   compression   for   σ 3 < 0 )

Appendix C. Analytic Expressions of the Uncertainty for Hashin-Type Failure Criteria

The superscript “I” denotes the corresponding interval quantity while m ( ) and r ( ) are the midpoint and radius, respectively.
r M T F I = 2 σ 2 2 r ( F 2 t I ) m ( F 2 t I ) m ( F 2 t I ) 2 r ( F 2 t I ) 2 2 + 2 σ 12 2 + σ 13 2 r ( F 12 I ) m ( F 12 I ) m ( F 12 I ) 2 m ( F 12 I ) 2 2
r M C F I = 2 σ 2 2 r ( F 2 t I ) m ( F 2 t I ) m ( F 2 t I ) 2 r ( F 2 t I ) 2 2 + 2 σ 12 2 r ( F 12 I ) m ( F 12 I ) m ( F 12 I ) 2 r ( F 12 I ) 2 2 + + 2 σ 23 2 r ( F 23 I ) m ( F 23 I ) m ( F 23 I ) 2 r ( F 23 I ) 2 2
r F T F I = 2 σ 1 2 r ( F 1 t I ) m ( F 1 t I ) m ( F 1 t I ) 2 r ( F 1 t I ) 2 2 + 2 σ 12 2 + σ 13 2 r ( F 12 I ) m ( F 12 I ) m ( F 12 I ) 2 r ( F 12 I ) 2 2
r F C F I = σ 1 r F 1 c I m ( F 1 c I ) 2 r ( F 1 c I ) 2
r F M S I = σ 1 r F 1 c I m ( F 1 c I ) 2 r ( F 1 c I ) 2 + 2 σ 12 2 + σ 13 2 r ( F 12 I ) m ( F 12 I ) m ( F 12 I ) 2 r ( F 12 I ) 2 2
r D T I = 2 σ 3 2 r ( F 3 t I ) m ( F 3 t I ) m ( F 3 t I ) 2 r ( F 3 t I ) 2 2 + 2 σ 13 2 r ( F 13 I ) m ( F 13 I ) m ( F 13 I ) 2 r ( F 13 I ) 2 2 + + 2 σ 23 2 r ( F 23 I ) m ( F 23 I ) m ( F 23 I ) 2 r ( F 23 I ) 2 2
r D C I = 2 σ 3 2 r ( F 3 c I ) m ( F 3 c I ) m ( F 3 c I ) 2 r ( F 3 c I ) 2 2 + 2 σ 13 2 r ( F 13 I ) m ( F 13 I ) m ( F 13 I ) 2 r ( F 13 I ) 2 2 + + 2 σ 23 2 r ( F 23 I ) m ( F 23 I ) m ( F 23 I ) 2 r ( F 23 I ) 2 2

References

  1. Pantelakis, S.; Tserpes, K. (Eds.) Revolutionizing Aircraft Materials and Processes; Springer International Publishing: Cham, Switzerland, 2020; ISBN 978-3-030-35345-2. [Google Scholar]
  2. Chamis, C. Simplified Composite Micromechanics Equations for Strength, Fracture Toughness, Impact Resistance and Environmental Effects; Lewis Research Center, National Aeronautics and Space Administration: Cleveland, OH, USA, 1984.
  3. Hinton, M.; Soden, P.; Kaddour, A.-S. Failure Criteria in Fibre-Reinforced-Polymer Composites; Elsevier: Amsterdam, The Netherlands, 2004; ISBN 978-0-08-044475-8. [Google Scholar]
  4. Hashin, Z. Failure Criteria for Unidirectional Fiber Composites. J. Appl. Mech. 1980, 47, 329–334. [Google Scholar] [CrossRef]
  5. Tserpes, K.I.; Papanikos, P.; Kermanidis, T. A Three-Dimensional Progressive Damage Model for Bolted Joints in Composite Laminates Subjected to Tensile Loading: Failure of tensile-loaded bolted composite joints. Fatigue Fract. Eng. Mater. Struct. 2001, 24, 663–675. [Google Scholar] [CrossRef]
  6. Tserpes, K.I.; Labeas, G.; Papanikos, P.; Kermanidis, T. Strength Prediction of Bolted Joints in Graphite/Epoxy Composite Laminates. Compos. Part B Eng. 2002, 33, 521–529. [Google Scholar] [CrossRef]
  7. Kolks, G.; Tserpes, K.I. Efficient Progressive Damage Modeling of Hybrid Composite/Titanium Bolted Joints. Compos. Part Appl. Sci. Manuf. 2014, 56, 51–63. [Google Scholar] [CrossRef]
  8. Stamopoulos, A.G.; Tserpes, K.I.; Dentsoras, A.J. Quality Assessment of Porous CFRP Specimens Using X-Ray Computed Tomography Data and Artificial Neural Networks. Compos. Struct. 2018, 192, 327–335. [Google Scholar] [CrossRef]
  9. Murthy, P.L.N.; Mital, S.K.; Shah, A.R. Probabilistic Micromechanics/Macromechanics for Ceramic Matrix Composites. J. Compos. Mater. 1998, 32, 679–699. [Google Scholar] [CrossRef]
  10. Stock, T.A.; Bellini, P.X.; Murthy, P.L.; Chamis, C.C. A Probabilistic Approach to Composite Micromechanics. In Proceedings of the 29th Structures, Structural Dynamics and Materials Conference, Williamsburg, VA, USA, 18–20 April 1988. [Google Scholar]
  11. Chen, N.; Yu, D.; Xia, B.; Liu, J.; Ma, Z. Interval and Subinterval Homogenization-Based Method for Determining the Effective Elastic Properties of Periodic Microstructure with Interval Parameters. Int. J. Solids Struct. 2017, 106–107, 174–182. [Google Scholar] [CrossRef]
  12. Wang, L.; Wang, X.; Su, H.; Lin, G. Reliability Estimation of Fatigue Crack Growth Prediction via Limited Measured Data. Int. J. Mech. Sci. 2017, 121, 44–57. [Google Scholar] [CrossRef]
  13. Wang, L.; Liu, D.; Yang, Y.; Wang, X.; Qiu, Z. A Novel Method of Non-Probabilistic Reliability-Based Topology Optimization Corresponding to Continuum Structures with Unknown but Bounded Uncertainties. Comput. Methods Appl. Mech. Eng. 2017, 326, 573–595. [Google Scholar] [CrossRef]
  14. Alazwari, M.A.; Rao, S.S. Modeling and Analysis of Composite Laminates in the Presence of Uncertainties. Compos. Part B Eng. 2019, 161, 107–120. [Google Scholar] [CrossRef]
  15. Alazwari, M.A.; Rao, S.S. Interval-Based Uncertainty Models for Micromechanical Properties of Composite Materials. J. Reinf. Plast. Compos. 2018, 37, 1142–1162. [Google Scholar] [CrossRef]
  16. Rao, S.S.; Alazwari, M.A. Failure Modeling and Analysis of Composite Laminates: Interval-Based Approaches. J. Reinf. Plast. Compos. 2020, 39, 817–836. [Google Scholar] [CrossRef]
  17. Moore, R.E.; Kearfott, R.B.; Cloud, M.J. Introduction to Interval Analysis; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2009; ISBN 978-0-89871-669-6. [Google Scholar]
  18. Hansen, E.R.; Walster, G.W. Global Optimization Using Interval Analysis. In Monographs and Textbooks in Pure and Applied Mathematics, 2nd ed.; Revised and Expanded; Marcel Dekker: New York, NY, USA, 2004; ISBN 978-0-8247-4059-7. [Google Scholar]
  19. Ratschek, H.; Rokne, J. Computer Methods for the Range of Functions; Ellis Horwood Series in Mathematics and Its, Applications; Horwood, E., Ed.; Halsted Press: Chichester, UK; New York, NY, USA, 1984; ISBN 978-0-85312-703-1. [Google Scholar]
  20. Krämer, W. High Performance Verified Computing Using C-XSC. Comput. Appl. Math. 2013, 32, 385–400. [Google Scholar] [CrossRef]
  21. Kulisch, U.; Hammer, R.; Hocks, M.; Ratz, D. C++ Toolbox for Verified Computing I; Springer: Berlin/Heidelberg, Germany, 1995; ISBN 978-3-642-79653-1. [Google Scholar]
  22. Caprani, O.; Madsen, K. Mean Value Forms in Interval Analysis. Computing 1980, 25, 147–154. [Google Scholar] [CrossRef]
  23. Rall, L.B. Mean Value and Taylor Forms in Interval Analysis. SIAM J. Math. Anal. 1983, 14, 223–238. [Google Scholar] [CrossRef]
  24. Baumann, E. Optimal Centered Forms. BIT 1988, 28, 80–87. [Google Scholar] [CrossRef]
  25. Sotiropoulos, D.G.; Grapsa, T.N. Optimal Centers in Branch-and-Prune Algorithms for Univariate Global Optimization. Appl. Math. Comput. 2005, 169, 247–277. [Google Scholar] [CrossRef]
  26. Ratz, D. Automatische Ergebnisverifikation Bei Globalen Optimierungsproblemen. Ph.D. Dissertation, Universität Karlsruhe (TH), Institute for Applied Mathematics, Karlsruhe, Germany, 1992. [Google Scholar] [CrossRef]
  27. Ratz, D.; Csendes, T. On the Selection of Subdivision Directions in Interval Branch-and-Bound Methods for Global Optimization. J. Glob. Optim. 1995, 7, 183–207. [Google Scholar] [CrossRef]
  28. Soden, P. Lamina Properties, Lay-up Configurations and Loading Conditions for a Range of Fibre-Reinforced Composite Laminates. Compos. Sci. Technol. 1998, 58, 1011–1022. [Google Scholar] [CrossRef]
Figure 1. Output uncertain level versus input level for (a) AS4-3501/6 and (b) 21xK43-LY556 composite ply.
Figure 1. Output uncertain level versus input level for (a) AS4-3501/6 and (b) 21xK43-LY556 composite ply.
Mca 27 00038 g001
Figure 2. The effect of input parameters on the range of (a) the E11, (b) the E22, (c) the G12, (d) the G23, (e) the ν12, and (f) the ν23 for the AS4-3501/6 composite ply. The main parameter varied within ±5% of its nominal value while the other parameters were kept equal to the nominal value.
Figure 2. The effect of input parameters on the range of (a) the E11, (b) the E22, (c) the G12, (d) the G23, (e) the ν12, and (f) the ν23 for the AS4-3501/6 composite ply. The main parameter varied within ±5% of its nominal value while the other parameters were kept equal to the nominal value.
Mca 27 00038 g002aMca 27 00038 g002b
Figure 3. The effect of input parameters on the range of (a) the E11, (b) the E22, (c) the G12, (d) the G23, (e) the ν12, and (f) the ν23 for the 21xK43-LY556 composite ply. The main parameter varied within ±5% of its nominal value while the other parameters were kept equal to the nominal value.
Figure 3. The effect of input parameters on the range of (a) the E11, (b) the E22, (c) the G12, (d) the G23, (e) the ν12, and (f) the ν23 for the 21xK43-LY556 composite ply. The main parameter varied within ±5% of its nominal value while the other parameters were kept equal to the nominal value.
Mca 27 00038 g003aMca 27 00038 g003b
Figure 4. Plots of the upper and lower bounds, the nominal value of failure criteria, and the probability of failure regarding the stress multiplication factor for the AS4-3501/6 UD ply.
Figure 4. Plots of the upper and lower bounds, the nominal value of failure criteria, and the probability of failure regarding the stress multiplication factor for the AS4-3501/6 UD ply.
Mca 27 00038 g004
Figure 5. Plots of the upper and lower bounds, the nominal value of failure criteria, and the probability of failure regarding the stress multiplication factor for the 21xK43-LY556 UD ply.
Figure 5. Plots of the upper and lower bounds, the nominal value of failure criteria, and the probability of failure regarding the stress multiplication factor for the 21xK43-LY556 UD ply.
Mca 27 00038 g005
Table 1. Nominal values and uncertainty intervals of the elastic properties of the fibers and the matrix of the AS4-3501/6 composite ply. The data have been taken from [28].
Table 1. Nominal values and uncertainty intervals of the elastic properties of the fibers and the matrix of the AS4-3501/6 composite ply. The data have been taken from [28].
PropertySymbolNominal ValueUncertainty Interval (Level 5%)
Fiber volume fractionVf0.60[0.569, 0.631]
Longitudinal modulus of the fibers (GPa)Ef11225[213.750, 236.250]
Transverse modulus of the fibers (GPa)Ef2215[14.250, 15.750]
In-plane shear modulus of the fibers (GPa)Gf1215[14.250, 15.750]
Transverse shear modulus of the fibers (GPa)Gf237[6.649, 7.351]
Poisson’s ratio of the fibersνf120.20[0.190, 0.211]
Poisson’s ratio of the fibersνf230.20[0.190, 0.211]
Matrix volume fractionVm0.40[0.380, 0.421]
Young’s modulus of the matrix (GPa)Em4.20[3.989, 4.411]
Shear modulus of the matrix (GPa)Gm1.567[1.488, 1.646]
Poisson’s ratio of the matrixνm0.34[0.323, 0.358]
Table 2. Nominal values and uncertainty intervals of the elastic properties of the fibers and the matrix of the 21xK43-LY556 composite ply. The data have been taken from [28].
Table 2. Nominal values and uncertainty intervals of the elastic properties of the fibers and the matrix of the 21xK43-LY556 composite ply. The data have been taken from [28].
PropertySymbolNominal ValueUncertainty Interval (Level 5%)
Fiber volume fractionVf0.60[0.569, 0.631]
Longitudinal modulus of the fibers (GPa)Ef1180[76.000, 84.000]
Transverse modulus of the fibers (GPa)Ef2280[76.000, 84.000]
In-plane shear modulus of the fibers (GPa)Gf1233.33[31.663, 34.997]
Transverse shear modulus of the fibers (GPa)Gf2333.33[31.663, 34.997]
Poisson’s ratio of the fibersνf120.20[0.190, 0.211]
Poisson’s ratio of the fibersνf230.20[0.190, 0.211]
Matrix volume fractionVm0.40[0.380, 0.421]
Young’s modulus of the matrix (GPa)Em3.25[3.087, 3.413]
Shear modulus of the matrix (GPa)Gm1.24[1.177, 1.303]
Poisson’s ratio of the matrixνm0.35[0.332, 0.368]
Table 3. Nominal values and uncertainty intervals of the strengths of the fibers and the matrix of the AS4/3501-6 composite ply. The data have been taken from [28].
Table 3. Nominal values and uncertainty intervals of the strengths of the fibers and the matrix of the AS4/3501-6 composite ply. The data have been taken from [28].
PropertySymbolNominal
Value
Uncertainty
Interval (Level 5%)
Longitudinal tensile strength (MPa)Fft3350[3182.5, 3517.5]
Longitudinal compressive strength (MPa)Ffc2500[2375, 2625]
Tensile strength (MPa)Fmt69[65.549, 72.451]
Compressive strength (MPa)Fmc250[237.5, 262.5]
Shear strength (MPa)Fms50[47.5, 52.5]
Stress concentration factorkσ1.4[1.329, 1.470]
Stress intensity factorkτ1.0[0.949, 1.051]
Table 4. Nominal values and uncertainty intervals of the strengths of the fibers and the matrix of the 21xK43-LY556 composite ply. The data have been taken from [28].
Table 4. Nominal values and uncertainty intervals of the strengths of the fibers and the matrix of the 21xK43-LY556 composite ply. The data have been taken from [28].
PropertySymbolNominal ValueUncertainty Interval (Level 5%)
Longitudinal tensile strength (MPa)Fft21,500[2042.5, 2257.5]
Longitudinal compressive strength (MPa)Ffc1450[1377.5, 1522.5]
Tensile strength (MPa)Fmt80[76.0, 84.0]
Compressive strength (MPa)Fmc120[114.0, 126.0]
Shear strength (MPa)Fms60[57.0, 63.0]
Stress concentration factorkσ1.82[1.728, 1.912]
Stress intensity factorkτ1.0[0.949, 1.051]
Table 5. Nominal values and computed uncertainty intervals of the AS4/3501/6 composite ply elastic properties.
Table 5. Nominal values and computed uncertainty intervals of the AS4/3501/6 composite ply elastic properties.
Elastic PropertyNominal ValueInterval Rangemid ± radComputed
Uncertainty Level (%)
E11136.680[123.553, 150.470]137.011 ± 13.4589.82
E22 = E339.496[8.742, 10.292]9.517 ± 0.7758.14
G12 = G135.116[4.596, 5.690]5.143 ± 0.54710.63
G233.929[3.595, 4.286]3.940 ± 0.3458.75
v12 = v130.256[0.239, 0.274]0.256 ± 0.0176.64
v230.385[0.351, 0.420]0.386 ± 0.0348.79
Table 6. Nominal values and computed uncertainty intervals of the elastic properties of the 21xK43-LY556 composite ply.
Table 6. Nominal values and computed uncertainty intervals of the elastic properties of the 21xK43-LY556 composite ply.
Elastic PropertyNominal ValueInterval Rangemid ± radComputed Uncertainty Level (%)
E1149.300[44.647, 54.183]49.415 ± 4.7689.65
E22 = E3312.652[11.199, 14.307]12.753 ± 1.55412.18
G12 = G134.878[4.313, 5.522]4.917 ± 0.60412.28
G234.878[4.313, 5.522]4.917 ± 0.60412.28
v12 = v130.260[0.242, 0.278]0.260 ± 0.0186.72
v230.373[0.338, 0.410]0.374 ± 0.0369.55
Table 7. The computational effort with large-range interval parameters (level 60%).
Table 7. The computational effort with large-range interval parameters (level 60%).
Elastic
Property
AS4/3501/621xK43-LY556
FEGENBLLFEGENBLL
E1142004200
E22 = E33422210410621
G12 = G1316104210621
G2330167310621
v12 = v131062110621
v236234164369198987
Table 8. Nominal values and computed uncertainty intervals of the AS4/3501/6 composite ply strengths.
Table 8. Nominal values and computed uncertainty intervals of the AS4/3501/6 composite ply strengths.
StrengthsNominal ValueInterval Rangemid ± radComputed Uncertainty Level (%)
F1t2037.600[1842.211, 2242.832]2042.521 ± 200.3109.81
F2t49.286[44.591, 54.474]49.533 ± 4.9419.98
F1c3.917[3.544, 4.330]3.937 ± 0.3939.98
F2c = F3c178.571[161.564, 197.369]179.467 ± 17.9029.98
F12 = F1350.000[47.500, 52.500]50.000 ± 2.5005.00
F2350.000[45.238, 55.264]50.251 ± 5.0139.98
Table 9. Nominal values and computed uncertainty intervals of strengths of the 21xK43-LY556 composite ply.
Table 9. Nominal values and computed uncertainty intervals of strengths of the 21xK43-LY556 composite ply.
StrengthsNominal ValueInterval Rangemid ± radComputed Uncertainty Level (%)
F1t1322.0[1196.904, 1453.306]1325.105 ± 128.2009.67
F2t43.956[39.769, 48.583]44.176 ± 4.4079.98
F1c3.10[2.804, 3.427]3.116 ± 0.3119.98
F2c = F3c65.934[59.654, 72.875]66.265 ± 6.6109.98
F12 = F1360.00[57.0, 63.0]60.000 ± 3.0005.00
F2360.00[54.285, 66.316]60.301 ± 6.0159.98
Table 10. Input stresses (MPa) are computed by the FE model.
Table 10. Input stresses (MPa) are computed by the FE model.
Layer σ 1 σ 2 σ 3 σ 12 σ 23 σ 13
1138.4600.5100.026−18.250−0.2100.030
45°231.886114.404−0.054128.6840.029−0.030
90°102.376−341.169−0.004−18.2500.029−0.044
−45°231.820114.404−0.014−165.1850.013−0.044
1138.4600.5100.026−18.2500.042−0.044
45°231.886114.404−0.054128.6840.042−0.044
90°102.376−341.169−0.004−18.2500.0420.004
−45°2318.820114.404−0.014−165.1850.0420.013
Table 11. Variation in the interval parameters and the probability of failure regarding the stress multiplication factor for the MTF. The bold number in the 3rd column corresponds to the first value of the failure criterion above 1.0 (failure) and the coloring of the values in the last column correspond to different intensity ranges of the probability of failure; the intensity increases from green to dark red.
Table 11. Variation in the interval parameters and the probability of failure regarding the stress multiplication factor for the MTF. The bold number in the 3rd column corresponds to the first value of the failure criterion above 1.0 (failure) and the coloring of the values in the last column correspond to different intensity ranges of the probability of failure; the intensity increases from green to dark red.
Stress Multiplication FactorLower BoundNominal Value of the Failure CriterionUpper BoundProbability of Failure
4.3515350.668870.8170960.998170.0%
4.4001020.6838840.8354371.0205756.1%
4.5002210.715360.8738881.06754819.2%
4.6003180.7475370.9131961.11556631.4%
4.7002910.780380.9533181.16457942.8%
4.8000350.8138520.9942081.21453153.5%
4.9004270.8482511.0362291.26586563.7%
5.0004170.8832211.0789481.3180573.1%
5.1004070.9188961.1225291.37128982.1%
5.2003160.9552481.1669371.42553890.5%
5.3000610.9922441.2121321.48074898.4%
5.3207711.0000131.2216231.492343100.0%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sotiropoulos, D.G.; Tserpes, K. Interval-Based Computation of the Uncertainty in the Mechanical Properties and the Failure Analysis of Unidirectional Composite Materials. Math. Comput. Appl. 2022, 27, 38. https://doi.org/10.3390/mca27030038

AMA Style

Sotiropoulos DG, Tserpes K. Interval-Based Computation of the Uncertainty in the Mechanical Properties and the Failure Analysis of Unidirectional Composite Materials. Mathematical and Computational Applications. 2022; 27(3):38. https://doi.org/10.3390/mca27030038

Chicago/Turabian Style

Sotiropoulos, Dimitris G., and Konstantinos Tserpes. 2022. "Interval-Based Computation of the Uncertainty in the Mechanical Properties and the Failure Analysis of Unidirectional Composite Materials" Mathematical and Computational Applications 27, no. 3: 38. https://doi.org/10.3390/mca27030038

APA Style

Sotiropoulos, D. G., & Tserpes, K. (2022). Interval-Based Computation of the Uncertainty in the Mechanical Properties and the Failure Analysis of Unidirectional Composite Materials. Mathematical and Computational Applications, 27(3), 38. https://doi.org/10.3390/mca27030038

Article Metrics

Back to TopTop