Research article

Scrapping of PEKK-based thermoplastic composites retaining long fibers and their use for compression molded recycled parts

  • In this work, a novel method for a more sustainable recycling and cost-efficient manufacturing technique of polyether ketone ketone (PEKK) based thermoplastic composite materials is proposed to recover and reprocess waste and end-of-life materials in the aerospace industry. For the recycling of carbon fiber reinforced thermoplastics (CFrTP), an innovative scrapping process based on mechanical cutting was developed and the properties of the obtained scrap and the recycled panel were analyzed. Thus, a cutting tool was developed for the delamination of the input material so that long fibers can be retained in the resulting scrap. Different processing approaches of material scrapping were evaluated, aiming to obtain manageable scrap that can be subsequently used for a compression molding process. Additionally, an automatic process was evaluated to manage the scrap and perform the corresponding lay-up to manufacture high-quality thermoplastic composite products with recycled materials.

    Citation: Alejandro Sandá, Rocío Ruiz, Miguel Ángel Mafé, Jon Ander Sarasua, Antonio González-Jiménez. Scrapping of PEKK-based thermoplastic composites retaining long fibers and their use for compression molded recycled parts[J]. AIMS Materials Science, 2023, 10(5): 819-834. doi: 10.3934/matersci.2023044

    Related Papers:

    [1] Lorenza D'Elia . Γ-convergence of quadratic functionals with non uniformly elliptic conductivity matrices. Networks and Heterogeneous Media, 2022, 17(1): 15-45. doi: 10.3934/nhm.2021022
    [2] Andrea Braides, Valeria Chiadò Piat . Non convex homogenization problems for singular structures. Networks and Heterogeneous Media, 2008, 3(3): 489-508. doi: 10.3934/nhm.2008.3.489
    [3] Patrick Henning . Convergence of MsFEM approximations for elliptic, non-periodic homogenization problems. Networks and Heterogeneous Media, 2012, 7(3): 503-524. doi: 10.3934/nhm.2012.7.503
    [4] Leonid Berlyand, Volodymyr Rybalko . Homogenized description of multiple Ginzburg-Landau vortices pinned by small holes. Networks and Heterogeneous Media, 2013, 8(1): 115-130. doi: 10.3934/nhm.2013.8.115
    [5] Manuel Friedrich, Bernd Schmidt . On a discrete-to-continuum convergence result for a two dimensional brittle material in the small displacement regime. Networks and Heterogeneous Media, 2015, 10(2): 321-342. doi: 10.3934/nhm.2015.10.321
    [6] Fabio Camilli, Claudio Marchi . On the convergence rate in multiscale homogenization of fully nonlinear elliptic problems. Networks and Heterogeneous Media, 2011, 6(1): 61-75. doi: 10.3934/nhm.2011.6.61
    [7] Brahim Amaziane, Leonid Pankratov, Andrey Piatnitski . Homogenization of variational functionals with nonstandard growth in perforated domains. Networks and Heterogeneous Media, 2010, 5(2): 189-215. doi: 10.3934/nhm.2010.5.189
    [8] Antonio DeSimone, Martin Kružík . Domain patterns and hysteresis in phase-transforming solids: Analysis and numerical simulations of a sharp interface dissipative model via phase-field approximation. Networks and Heterogeneous Media, 2013, 8(2): 481-499. doi: 10.3934/nhm.2013.8.481
    [9] Liselott Flodén, Jens Persson . Homogenization of nonlinear dissipative hyperbolic problems exhibiting arbitrarily many spatial and temporal scales. Networks and Heterogeneous Media, 2016, 11(4): 627-653. doi: 10.3934/nhm.2016012
    [10] Peter Bella, Arianna Giunti . Green's function for elliptic systems: Moment bounds. Networks and Heterogeneous Media, 2018, 13(1): 155-176. doi: 10.3934/nhm.2018007
  • In this work, a novel method for a more sustainable recycling and cost-efficient manufacturing technique of polyether ketone ketone (PEKK) based thermoplastic composite materials is proposed to recover and reprocess waste and end-of-life materials in the aerospace industry. For the recycling of carbon fiber reinforced thermoplastics (CFrTP), an innovative scrapping process based on mechanical cutting was developed and the properties of the obtained scrap and the recycled panel were analyzed. Thus, a cutting tool was developed for the delamination of the input material so that long fibers can be retained in the resulting scrap. Different processing approaches of material scrapping were evaluated, aiming to obtain manageable scrap that can be subsequently used for a compression molding process. Additionally, an automatic process was evaluated to manage the scrap and perform the corresponding lay-up to manufacture high-quality thermoplastic composite products with recycled materials.



    Natural hazards such as landslides and subsidence have been recognized as important impediments to developing nations' sustainable development [1,2,3,4]. For example, on December 20, 2015, in Guangming, Shenzhen, China, a catastrophic landfill slope landslide claimed the lives of 69 individuals [5]. Natural hazard risk assessment and management will have short-term benefits in terms of severity reduction and long-term benefits in terms of achieving sustainable development goals [1].

    Slope stability evaluation is essential for analyzing and mitigating natural hazards in mountainous environments. Many attempts have been made to evaluate slope stability [6,7,8]. Owing to the inherent complexity and uncertainty, assessing slope stability for circular mode failure, a common problem, remains a challenge for practitioners and researchers [9]. Several methods for evaluating slope stability have been presented, with the limit equilibrium approach and numerical simulation method being the two most commonly used [10]. Limited equilibrium methods, such as the simplified Bishop, Spencer and Morgenstern-Price methods, are frequently implemented in practice. In general, soil material properties (unit weight, cohesion and friction angle) and the pore pressure ratio are required for limited equilibrium methods [11,12]. Numerical methods (e.g., finite element methods) have been extensively used to analyze slope stability. However, the major drawback is that their input parameters need to be back analyzed using in-situ measurements, which is not available in many cases [13]. Both methods have pros and cons. Finding the critical slip surface using the limit equilibrium method is difficult due to the large number of potential slip surfaces [14]. The numerical simulation method's accuracy is greatly influenced by the choice of constitutive models, mechanical parameters and boundary conditions, and it is frequently necessary to have a great deal of engineering expertise and to conduct on-site back analysis in order to make a reasonable choice and obtain reasonable results [15,16]. Consequently, predicting slope stability still presents considerable challenges.

    In recent years, machine learning (ML) models have gained attention for solving very complex, nonlinear and multivariable geotechnical problems [17,18,19]. Assessments of slope stability circular failure using methods based on soft computing are summarized in Table 1. Despite their reliable and precise outputs, however, most algorithms are not readily applicable in practice owing to their complicated training and modeling procedures and "black box" aspects, i.e., these models do not demonstrate a transparent and understandable relationship between inputs and output. Quinlan [20] developed the model tree algorithm to get over these limitations. It integrates principles from decision trees and linear regression. In addition to the widespread application of soft computing techniques, several studies have only been conducted on a limited amount of data, which might restrict the classifier's ability to generalize. In the current study, an updated database of 627 cases comprising unit weight, cohesion, internal friction angle, slope angle and height, pore pressure ratio and stability status of circular mode failure has been compiled. To predict slope stability, a new logistic model tree (LMT) model is developed. This algorithm is an intelligent choice for classification and decision-making since it can solve the classification problem by combining a tree model and a logistic regression (LR) technique. Adding LR to the leaves of the tree allows for a probabilistic interpretation of the model's output, making it easier to explain, as they represent a series of if-then-else rules. The LMT has been employed in predicting pillar stability in geotechnical engineering [21], but it has not yet been used to predict slope stability.

    Table 1.  A summary of the ML-based circular mode failure assessment of slope stability.
    Dataset (Stable/Failed) Input parameters Data Preprocessing ML techniques Reference
    422 (226/196) γ, c, ϕ, β, H, ru Data normalization MDMSE Zhang et al. [22]
    444 (224/220) γ, c, ϕ, β, H, ru Data normalization AdaBoost, GBM Bagging, XRT, RF, HGB Voting Stacked Lin et al. [23]
    19 (13/6) γ, c, ϕ, β, H, ru Data normalization K-means cluster Haghshenas et al. [24]
    153 (83/70) γ, c, ϕ, β, H, ru Data normalization and outlier removing KNN, SVM, SGD, GP, QDA, GNB, DT, ANN, Bagging ensemble, Heterogeneous ensemble Pham et al. [25]
    257 (123/134) γ, c, ϕ, β, H, ru XGB, RF, LR, SVM, BC, LDA, KNN, DT, MLP, GNB, XRT, Stacked ensemble Kardani et al. [26]
    87 (42/45) γ, c, ϕ, β, H, ru J48 Amirkiyaei and Ghasemi [27]
    221 (115/106) γ, c, ϕ, β, H, ru Data normalization ANN, SVM, RF, GBM Zhou et al. [28]
    148 (78/70) γ, c, ϕ, β, H, ru Data normalization LR, DT, RF, GBM, SVM, BP Qi and Tang [29]
    168 (84/84) γ, c, ϕ, β, H, ru Data normalization GP, QDA, SVM, ADB-DT, ANN, KNN Classifier ensemble Qi and Tang [15]
    107 (48/59) γ, c, ϕ, β, H, ru RF, SVM, Bayes, GSA Lin et al. [30]
    82 (49/33) γ, c, ϕ, β, H, ru NB Feng et al. [31]
    168 (84/84) γ, c, ϕ, β, H, ru Data normalization RBF, LSSVM, ELM Hoang and Bui [32]
    168 (84/84) γ, c, ϕ, β, H, ru Data normalization LSSVM Hoang and Pham [33]
    46 (17/29) γ, c, ϕ, β, H, ru Data normalization SVM Xue et al. [34]
    32 (14/18) γ, c, ϕ, β, H, ru ANN Lu and Rosenbaum [35]
    82 (38/44) γ, c, ϕ, β, H, ru BP Feng [36]
    Note: MDMSE-margin distance minimization selective ensemble, AdaBoost-adaptive boosting, GBM-gradient boosting machine, XRT-extremely randomized tree, RF-random forest, HGB-hist gradient boosting classifier, KNN-k-nearest neighbors, SVM-support vector machine, SGD-stochastic gradient descent, GP-Gaussian process, QDA-quadratic discriminant analysis, GNB-Gaussian naive Bayes, DT-decision tree, ANN-artificial neural network, XGB-extreme gradient boosting, BC-bagging classifier, LDA-linear discriminant analysis, MLP-multilayer perceptron, BP-back-propagation, QDA-quadratic discriminant analysis, ADB-DT-adaptive boosted decision tree, GSA-gravitational search algorithm, NB-naive Bayes, RBF-radial basis function, LSSVM-least squares support vector machine, ELM-extreme learning machine, γ-unit weight, c-cohesion, ϕ-angle of internal friction, β-slope angle, H-slope height, ru-pore pressure ratio.

     | Show Table
    DownLoad: CSV

    LR is a straightforward method with features such as stability, low variance and time-efficient training [37], but its prediction outputs are frequently biased. Decision trees are another ML technique for searching a less confined space of candidate models and obtaining nonlinear patterns from a database; they have low bias but high variance and instability, making them susceptible to overfitting. Landwehr et al. [38] presented the LMT methodology as a result. It is based on Quinlan's model tree approach [20] for dealing with regression problems by combining linear regression with decision tree models, and it is extended to classification problems. This section provides a basic introduction to LMT, whereas the seminal work by Landwehr et al. [38] provides a more complete description.

    A LR tree with LR functions built at the leaves constitutes an LMT. It has a set of leaves or terminal nodes T and an inner or non-terminal node set N. Each leaf t of the LMT model has correlated LR functions instead of classification labels or linear regression functions. The output vector is Y, and the input vector X is (X1, X2), for instance. S represents the complete instance space, which can be further partitioned into numerous subspaces St. Figure 1 displays a simple input space that has been partitioned into seven subspaces.

    Figure 1.  The input space is separated into seven subspaces (solid circles and square points represent different classes within the same case, and each subspace has its own function).
    S=tTSt,StSt= (1)

    The model determines LR functions for the seven subspaces represented in Figure 1. Figure 2 depicts the structure of the tree.

    Figure 2.  Example of a simplified LMT structure corresponding to the sample data shown in Figure 1 (adapted from Landwehr et al. [38]).

    In contrast to standard forms of LR, the LogitBoost technique for fitting additive LR models proposed by Friedman et al. [39] is employed for model construction here. The prediction probability is given by Eq (2).

    Pr(G=j|X=x)=eFj(x)Jk=1eFk(x) (2)

    where G denotes the output, J denotes the class labels, X denotes the inputs, and Fj(x) denotes the functions that the LMT will train in the tree's leaves, as follows:

    Fj(x)=Mm=1fmj(x)=αj0+SStαjSS (3)

    where m is the number of iterations, fmj represents the functions of input variables, α represents the intercepts and coefficients of the linear function, and S represents the variables of the subset St at the leaf t.

    A LMT can be established by the following steps: initial tree growth, tree splitting and stopping, and tree pruning. This section presents the basic idea; the reader is referred to Landwehr et al. [38] for additional detail. The M5P technique, which is commonly used for tree growth, can first construct a standard tree, after which a LR model can be established at each node [40,41]. This technique merely trains the model using case histories at each node in isolation, without taking the surrounding tree structure into account; therefore, in order for the LogitBoost algorithm to iteratively change Fj(x) to increase the fit in a natural way, another technique-one that can incrementally refine logistic model fit at high levels-is used [38]. The function fmj is introduced to Fj by changing one of the function's coefficients or introducing another variable (see Eq (3)). As a result, in the initial growing process, a LR tree is formed in the root using proper iteration numbers. The tree then begins to grow by resembling specified subsets (t) from the database (S) to the child nodes, utilizing the C4.5 splitting law [20] to increase the accuracy of the classification variable. The LR functions are generated in the child nodes by running the LogitBoost algorithm with the logistic model, weights and probability estimate from the previous iteration at the parent node. The splitting process is then repeated. When it comes to model fitting, the tree stops splitting when a node has less than 15 cases. After the tree is constructed, the tree pruning approach is used to trade off tree size and model complexity while maintaining predictive accuracy. After experimenting with several pruning initiatives, Landwehr et al. [38] employed the classification and regression trees pruning approach [42] to make pruning decisions while taking training error and model complexity into account. These three processes can be used to create a LMT.

    An updated database with 627 instances was obtained from previous studies [22,23,25,29,33,43] and can be found in Table S1 of the supplementary information file. The database includes the unit weight, cohesion, angle of internal friction, slope angle and height, pore pressure ratio and slope stability status. There are 311 positive (stable) and 316 negative (failed) samples. The statistics of the input features are summarized in Table 2. The data normalization is not carried out owing to the fact that tree-based methods are insensitive to feature scaling; they make decisions based on relative feature values and splits [44]. The database box plot is shown in Figure 3, where solid black spots represent "outliers", and the statistical features are shown in Table 2. The bottom and top quartiles are shown by horizontal lines, while the median values are represented by bold lines inside the boxes. Slopes with "failed" and "stable" instances are also demonstrated separately.

    Table 2.  Statistical values for features in the database.
    Input parameter Unit Min. Max. Mean Std. Dev.
    Unit weight (γ) kN/m3 0.492 33.16 20.185 7.044
    Cohesion (c) kPa 0 300 25.6 31.036
    Angle of internal friction (ϕ) 0 49.5 25.308 12.331
    Slope angle (β) 0.302 65 32.605 13.711
    Slope height (H) m 0.018 565 90.289 120.14
    Pore pressure ratio (ru) 0 1 0.254 0.26

     | Show Table
    DownLoad: CSV
    Figure 3.  Box plots show the variation among data points in the database.

    The classification metrics include accuracy evaluation indices (accuracy (Acc), Matthews correlation coefficient (Mcc), precision (Prec), recall (Rec) and F-score) that are calculated from the confusion matrix (see Figure 4).

    Figure 4.  Confusion matrix (2 × 2) for classification problem.

    Each row of the matrix represents the instances in an actual class, while each column represents the instances in a predicted class [45]. In predictive analytics, a confusion matrix is described as a table with two rows and two columns that provides the numbers of true positives (TPs), true negatives (TNs), false positives (FPs) and false negatives (FNs). The classification evaluation metrics derived from the confusion matrix results were used to compare the prediction performances of the models [46,47,48]:

    Acc=TN+TPTN+FP+FN+TP (4)
    Prec=TNFN+TNorTPTP+FP (5)
    Rec=TNTN+FPorTPTP+FN (6)
    Fscore=21Precision+1Recall (7)
    Mcc=TN×TPFP×FN(FN+TP)(FP+TN)(FP+TP)(FN+TN) (8)

    Acc is the sum of TP and TN and represents the percentage of correctly classified instances in the data as a whole. This metric measures the model's total prediction accuracy. If the data set is unbalanced, that is, the numbers of observations in different classes vary substantially, accuracy will be deceiving [49]. As a result, further assessment metrics such as precision, recall, F-score and Mcc were utilized to examine the model's performance further. Precision is known as positive predicted value, and recall is known as true positive rate (TPR). F-score is a generalized index that evaluates the performance of both recall and precision and ranges from 0 (worst value) to 1 (best value). Mcc denotes the degree of agreement between observed and predicted classes of failed and stable instances [48]. It is a standard metric used by statisticians that accepts values ranging from −1 to 1. Mcc value of −1 indicates complete disagreement (strong negative association), a value of 1 indicates complete agreement (strong positive association), and a value of 0 indicates that the prediction was unrelated to the ground truth (very weak or no correlation between dependent and independent variables). Additionally, another succinct metric, i.e., area under the receiver operating characteristic curve (AUC), was employed for effective evaluation. This curve is expressed in a graphical plot that shows a binary classification system's diagnostic capacity as its discrimination threshold changes. The AUC is plotted by comparing the TPR vs. the false positive rate (FPR) at various threshold levels. An acceptable classification model should have an AUC near to 1. Table 3 shows the main rule for defining discrimination based on AUC value [25].

    Table 3.  Rule for classifying discrimination based on AUC value.
    Discrimination category AUC value
    No discrimination AUC = 0.5
    Acceptable 0.7 ≤ AUC < 0.8
    Excellent 0.8 ≤ AUC < 0.9
    Outstanding 0.9 ≤ AUC

     | Show Table
    DownLoad: CSV

    Waikato Environment for Knowledge Analysis (WEKA) software was used for developing a model based on the acquired data set. WEKA is a collection of ML algorithms that supports data mining tasks by providing a wide range of tools that can be used for data pre-processing, classification, clustering, regression, association and visualization [50]. The database was divided into two parts: training (80%) and testing (20%). The training set contained 249 stable instances and 253 failed cases, while the test set contained 125 instances, 62 of which were stable and 63 of which were failed instances. The stable-to-failed instance ratios in the training and testing sets were close to one, indicating that the distribution of these two instances does not necessitate a cost-sensitive technique to address problem due to imbalance [51,52].

    Every terminal node (or leaf) of the tree model was trained and updated using LR models during training (see Section 2). The minimal number of instances for the LMT model was determined to be 15 based on predictive performance, readily applicable tree structure and total training data. The size of the tree was 45, and it had 23 leaves. LogitBoost had a weight trimming value of 0.2 and a number of iterations of −1. Figure 5 depicts a representation of a tree generated by LMT. The LMT model has 23 logistic functions (LMs), and their detailed expressions are shown in Table 4. It is important to notice that some of the functions in Table 4 do not incorporate all of the parameters that have been chosen. The LM1 function for stable slopes in Table 4, for example, does not account for unit weight (γ). The simple logistic technique is used in the LMT training phase [38]. The goal of the simple logistic method is to control the parameter numbers as simply and straightforwardly as possible. New parameters are gradually introduced during training in order to improve the performance of each function at each node of the tree (see Section 2). This can also help to avoid the issue of model significance in LR, especially when using multiple logistic functions to build a full logistic model with all parameters. However, just a few of the functions had fewer parameters than those selected, showing that most of the factors we chose have an influence on predictive performance.

    Figure 5.  A LMT's structure.
    Table 4.  Logistic functions in the LMT model.
    No. Regression model for the slope stability condition No. Regression model for the slope stability condition
    LM1 FS=0.3+14.53c+4.62ϕ5.83β9.17H0.69ru
    FF=0.314.53c4.62ϕ+5.83β+9.17H+0.69ru
    LM13 FS=1.870.01c0.03ϕ+0.04β+0H+1.3ru
    FF=1.87+0.01c+0.03ϕ0.04β0H1.3ru
    LM2 FS=0.360.04γ+0.02c0.04ϕ+0.13β0.05H
    FF=0.36+0.04γ0.02c+0.04ϕ0.13β+0.05H
    LM14 FS=6.73+0.22γ0.01c0.03ϕ+0.05β0H+1.3ruFF=6.730.22γ+0.01c+0.03ϕ0.05β+0H1.3ru
    LM3 FS=0.730.07γ0.02c0.13ϕ+0.11β+0.08H+3.03ru
    FF=0.73+0.07γ+0.02c+0.13ϕ0.11β0.08H3.03ru
    LM15 FS=10.690.01c0.09ϕ+0.28β+0H+2.97ruFF=10.69+0.01c+0.09ϕ0.28β0H2.97ru
    LM4 FS=5.08+0.02γ0.02c0.1ϕ+0.28β+0H8.66ru
    FF=5.080.02γ+0.02c+0.1ϕ0.28β0H+8.66ru
    LM16 FS=8.310.01c+0.02ϕ+0.02β+0H1.59ruFF=8.31+0.01c0.02ϕ0.02β0H+1.59ru
    LM5 FS=7.18+0.01γ0.01c0.04ϕ+0β
    FF=7.180.01γ+0.01c+0.04ϕ0β
    LM17 FS=2.050.01c+0.02ϕ+0.02β+0H1.67ru
    FF=2.05+0.01c0.02ϕ0.02β0H+1.67ru
    LM6 FS=4.40.51γ+0.03ϕ+0.09β+0.01H
    FF=4.4+0.51γ0.03ϕ0.09β0.01H
    LM18 FS=7.6+0.2γ0.04c0.1ϕ0.01β0.08H1.42ruFF=7.60.2γ+0.04c+0.1ϕ+0.01β+0.08H+1.42ru
    LM7 FS=27.081.49γ0.09c0.17ϕ0.22β0.02H+23.63ru
    FF=27.08+1.49γ+0.09c+0.17ϕ+0.22β+0.02H23.63ru
    LM19 FS=4.4+0.58γ0.04c0.07ϕ0.01β+0.01H1.72ru
    FF=4.40.58γ+0.04c+0.07ϕ+0.01β0.01H+1.72ru
    LM8 FS=9.31+0.08γ0.03c+0.01ϕ+0.03β0H+0.31ru
    FF=9.310.08γ+0.03c0.01ϕ0.03β+0H0.31ru
    LM20 FS=5.13+0.25γ0.11c0.05ϕ+0.08β0.01H+5.74ru
    FF=5.130.25γ+0.11c0.05ϕ0.08β+0.01H5.74ru
    LM9 FS=11.33+0.03γ+0.02c0.04ϕ0.31β+0.02H23.05ru
    FF=11.330.03γ0.02c+0.04ϕ+0.31β0.02H+23.05ru
    LM21 FS=10.75+0.74γ0c0.16ϕ+0.09β0H2.68ru
    FF=10.750.74γ+0c+0.16ϕ0.09β+0H+2.68ru
    LM10 FS=4.750.21γ0.01c0.06ϕ+0.07β+0H2.67ru
    FF=4.75+0.21γ+0.01c+0.06ϕ0.07β0H+2.67ru
    LM22 FS=90.05+1.93γ0.01c+0.41ϕ+0.48β+0H+0.58ru
    FF=90.051.93γ+0.01c0.41ϕ0.48β0H0.58ru
    LM11 FS=0.170.01γ+0c+0.02ϕ+0β1.67ru
    FF=0.17+0.01γ0c0.02ϕ0β+1.67ru
    LM23 FS=98.1+2.52γ+0.86c1.19ϕ0β+0.01H+3.15ru
    FF=98.12.52γ0.86c+1.19ϕ+0β0.01H3.15ru
    LM12 FS=43.040.01γ+0c+0.73ϕ0.2β+43.7ru
    FF=43.04+0.01γ0c0.73ϕ+0.2β43.7ru
    Note: Angle of internal friction (ϕ) and slope angle (β) are in degrees.

     | Show Table
    DownLoad: CSV

    The proposed LMT model was quantified using several performance metrics based on confusion matrices. The confusion matrices of the model in training and test sets were then obtained, as shown in Table 5.

    Table 5.  Confusion matrices and performance metric results for training and testing datasets.
    Dataset Actual Predicted Acc (%) Mcc AUC Prec Rec F-score
    Stable Failed
    Training Stable 237 12 92.23 0.846 0.974 0.898 0.952 0.924
    Failed 27 226 0.950 0.893 0.921
    Testing Stable 51 11 85.60 0.713 0.907 0.879 0.823 0.850
    Failed 7 56 0.836 0.889 0.862

     | Show Table
    DownLoad: CSV

    The proposed LMT model predicts two classes (stable and failed). Table 5 displays the confusion matrices that illustrate training and prediction results. The number of cases that could have been correctly predicted is indicated by the values along the major diagonal. The results of Table 5 show that the LMT correctly classified the majority of the cases.

    The prediction outcomes of classification problems can be assessed using a variety of metrics, such as Acc, Mcc, AUC, recall, precision and F-score. The LMT model performed well in both the training set (Acc = 92.23%, Mcc = 0.846, AUC = 0.974, F-score for stable state = 0.924 and F-score for failed state = 0.921) and testing set (Acc = 85.60%, Mcc = 0.713, AUC = 0.907, F-score for stable state = 0.850 and F-score for failed state = 0.862). The employed metrics show that the prediction's results are accurate and acceptable.

    The LMT model achieved better prediction performance (Acc = 0.856 and AUC = 0.907) in comparison to the SVM model (Acc = 0.812 and AUC = 0.796), SGD model (Acc = 0.640 and AUC = 0688), QDA model (Acc = 0.788 and AUC = 0.817), GNB model (Acc = 0.812 and AUC = 0.775), DT model (Acc = 0.788 and AUC = 0.829) and RT model (Acc = 0.788 and AUC = 0.904), reported by Pham et al. [25], respectively, for the test data. When using ML techniques to predict slope stability, Acc generally ranges from 0.640 to 0.812, according to the results of previous study by Pham et al. [25]. Meanwhile, in the present study, it is 0.856 for the test data set. However, due to the use of different data sets, a comparison between these results is unwarranted. A project that uses different data sets is needed to give a generalized model to predict slope stability. Additionally, the results for relatively small sample sizes (less than 100) are not presented or compared. These comparative results demonstrate conclusively that the LMT model is capable of improved generalization performance compared to other models in the literature. The potential reasons underlying this observation are that the LMT model can capture interactions between features effectively, the branching structure of the decision tree can identify feature interactions, and LR at the leaves can model the relationships between these interactions and the target variable.

    Two case studies are analyzed using our proposed LMT-based slope stability prediction model to determine its efficacy and applicability in engineering practice. The LMT model predicted ten slope stability status events in two different projects. The field data was obtained from the available literature, which includes the Shao Jiazhuang slope failure [53] and the Daguangbao landslide [54].

    The slope at Shao Jiazhuang village in Guizhou province, China (see Figure 6 (a)) [53], is utilized as a real-world example to analyze slope stability with the proposed LMT model and validate its efficacy and practicality in engineering practice application. As depicted in Figure 6(c)(e), during the survey, there were localized shallow surface collapses and new cracks on the slope's left side [53].

    Figure 6.  Slope failure area at Shao Jiazhuang, China: (a) location map, (b) slope collapse, (c), (d) collapse zone, and (e) crack.

    The input data shown in Table 6 were used by the LMT to predict the stability of the case slope based on the material parameters and geometry feature parameters of the case slope. The slope is failed, according to the prediction's outcome. This study demonstrates that the LMT model is reliable for predicting slope stability and that it can be applied to a variety of geotechnical applications.

    Table 6.  Application of the proposed LMT model in real-world slope stability prediction.
    S. No. γ (kN/m3) c (kPa) ϕ (°) β (°) H (m) ru Actual status LMT predicted result
    1 20 22.4 28 39.47 14 0 Failed Failed

     | Show Table
    DownLoad: CSV

    The Daguangbao landslide, is one of the few extremely big landslides known to exist worldwide, with a size of over 1 billion m3 [55]. It is also the most extensive and largest landslide ever recorded in Chinese historical records. Because of its massive volume, unusual genetic mechanism and complex movement process [56,57,58], the Daguangbao landslide has drawn a lot of attention and interest.

    The stability evaluation is carried out on 9 slopes, approved by specialists, from the Daguangbao landslide, of which 5 are stable, and 4 failed [54]. The index values of the samples and predicted results of the LMT model are compared to the fuzzy discriminant method and the unascertained measure method in Table 7.

    Table 7.  Predicted results using LMT, fuzzy discriminant and unascertained measure methods.
    S. No. γ (kN/m3) c (kPa) ϕ (°) β (°) H (m) ru Actual status LMT proposed method Fuzzy discriminant method [54] Unascertained measure method [54]
    1 27 32 33 42.4 289 0.25 Stable Stable Stable Stable
    2 20.41 33 11 16 46 0.2 Failed Failed Stable* Failed
    3 21.43 0 20 20 61 0.5 Failed Failed Failed Failed
    4 19.63 11.97 20 22 21.19 0.405 Failed Failed Failed Failed
    5 18.68 26.34 15 35 8.23 0.25 Failed Failed Failed Stable*
    6 27 50 40 42 407 0.25 Stable Stable Stable Stable
    7 27.3 14 31 50 92 0.25 Stable Stable Stable Failed*
    8 21.4 10 30.34 30 20 0.25 Stable Stable Stable Stable
    9 25 46 35 46 393 0.25 Stable Stable Stable Stable
    Note: (*) wrongly predicted.

     | Show Table
    DownLoad: CSV

    The prediction outcomes in Table 7 indicate that the slope stability was predicted correctly for all cases. The theory of unascertained measure and its mathematical processing were first put forward by Wang [59] in 1990. On this basis, Liu et al. [60,61] established the unascertained mathematical theory and applied it to decision-making problems. Unascertained information is a new kind of uncertainty information which is different from fuzzy information, random information and grey information. Compared with other evaluation methods, the unascertained measure method has the advantages of non-negativity, normalization and additivity, which also ensures the order of evaluation space. A fuzzy discriminant method [62] is a decision classification to construct a numerical tabular knowledge base from historical cases, and it derives inferences from particular case histories using discrimination and connectivity analyses which are based on a theory of fuzzy relations. For further details regarding these methods, readers may refer to [54,59,60,61,62].

    To demonstrate a use of the proposed LMT model for risk analysis, we use the case history data in Table 7. The following lists typical input data for the Shao Jiazhuang, China, slope failure area: γ = 20 kN/m3, c = 22.4 kPa, ϕ = 28°, β = 39.47°, H = 14 m, and ru = 0. Then, using the tree structure from Figure 5, we use the appropriate functions to compute the probability of failure (PoF). For example, to make a prediction for a new input instance, start at the root node of the LMT model (e.g., the value of c in this case is 22.4 kPa) and follow the path through the tree according to the splitting rules. Then, proceed to the end function at a leaf node (for this instance, LM8). The function values (FS and FF) are first calculated using the LM8 equations in Table 4. Then, using Eq (2), the probability of slope stability is given by Eq (11):

    FS=9.31+0.08γ0.03c+0.01ϕ+0.03β0H+0.31ru (9)
    FF=9.310.08γ+0.03c0.01ϕ0.03β+0H0.31ru (10)
    PS=eFSeFS+eFF,PF=eFFeFS+eFF (11)

    Finally, we can get the results: PS = 1%, and PF = 99%, implying that this instance has a PoF equal to 99%. This also accords with the actual facts. Such probability results can then be incorporated into risk analysis with the associated failure cost estimate.

    For the assessment of slope stability, a novel application of LMT is proposed. The tree structure and corresponding functions are used to assess slope stability, given information on several features, such as slope height (H), slope angle (β), cohesion (c), internal friction angle (ϕ), unit weight (γ) and pore pressure ratio (ru). LogitBoost learns the LMT by utilizing a larger database obtained from the literature. The results show that the LMT model can effectively predict slope stability. A testing set was used to validate the trained LMT model. Furthermore, real-world application to new cases—which had not previously been used for training-was used to validate the proposed LMT model. Comparative study and engineering application results show that the LMT model has the best prediction effect and is the best, optimal model. Furthermore, the results indicate that the LMT technique can provide useful information regarding the probability of slope failure, allowing it to be used for risk analysis of slope stability. Finally, the main advantages of LMT over other soft computing models commonly used for slope stability prediction are that it can be trained easily (even with more input parameters) and that its tree structure, with an LR function in each leaf, explicitly demonstrates the relationship between the inputs and predictive output. Furthermore, the LMT has the potential to be used to solve other geotechnical problems in the future due to its intuitive features and ease of implementation. It is recommended for future work that considering tangents (i.e., tan ϕ and tan β) could improve performance.

    The authors declare they have not used artificial intelligence (AI) tools in the creation of this article.

    This work was supported by the National Key Research and Development Plan of China under Grant No. 2021YFB2600703. The authors also wish to thank CEU San Pablo University Foundation for the funds dedicated to the ARIE Research Group, through Project Ref. EC01/0720-MGI23RGL, provided by the CEU San Pablo University.

    The authors declare there is no conflict of interest.



    [1] Fontaine P, Weiss-Hortala E, Botaro V, et al. (2021) Impact of atmosphere on recovered carbon fibers from poly ether ether ketone (PEEK) based composites during thermoconversion. Waste Biomass Valor 12: 6389–6402. https://doi.org/10.1007/s12649-021-01445-7 doi: 10.1007/s12649-021-01445-7
    [2] Borjan D, Knez Ž, Knez M (2021) Recycling of carbon fibre-reinforced composites—difficulties and future perspectives. Materials 12: 4191–4204. https://doi.org/10.3390/ma14154191 doi: 10.3390/ma14154191
    [3] Directive 2008/98/EC of the European Parliament and of the Council of 19 November 2008 on waste and repealing certain Directives. Available from: http://data.europa.eu/eli/dir/2008/98/2018-07-05.
    [4] Chawla KK (2012) Composite Materials Science and Engineering, 3Eds., New York: Springer New York.
    [5] Howarth J, Mareddy SSR, Mativenga PT (2014) Energy intensity and environmental analysis of mechanical recycling of carbon fibre composite. J Clean Prod 81: 46–50. https://doi.org/10.1016/j.jclepro.2014.06.023 doi: 10.1016/j.jclepro.2014.06.023
    [6] Pickering SJ, Yip H, Kennerley JR, et al. (2000) The recycling of carbon fibre composites using a fluidised bed process, In: Gibson AG, FRC 2000–Composites for the Millennium, Cambridge: Woodhead Publishing.
    [7] Lam KL, Oyedun AO, Cheung KY, et al. (2011) Modelling pyrolysis with dynamic heating. Chem Eng Sci 66: 6505–6514. https://doi.org/10.1016/j.ces.2011.09.013 doi: 10.1016/j.ces.2011.09.013
    [8] Oliveux G, Dandy L, Leeke GA (2014) A step change in the recycling of composite materials. 16th European Conference on Composite Materials, ECCM 2014, 2–5.
    [9] Yang Y, Boom R, Irion B, et al. (2012) Recycling of composite materials. Chem Eng Process 51: 53–68. https://doi.org/10.1016/j.cep.2011.09.007 doi: 10.1016/j.cep.2011.09.007
    [10] Li X, Bai R, McKechnie J (2016) Environmental and financial performance of mechanical recycling of carbon fibre reinforced polymers and comparison with conventional disposal routes. J Clean Prod 127: 451–460. https://doi.org/10.1016/j.jclepro.2016.03.139 doi: 10.1016/j.jclepro.2016.03.139
    [11] Vincent GA (2019) Recycling of thermoplastic composite laminates: The role of processing. University of Twente. https://doi.org/10.3990/1.9789036548526 doi: 10.3990/1.9789036548526
    [12] Vincent GA, de Bruijn TA, Wijskamp S, et al. (2019) Shredding and sieving thermoplastic composite scrap: Method development and analyses of the fibre length distributions. Compos Part B-Eng 176: 107197. https://doi.org/10.1016/j.compositesb.2019.107197 doi: 10.1016/j.compositesb.2019.107197
    [13] Li H, Englund K (2017) Recycling of carbon fibre-reinforced thermoplastic composite wastes from the aerospace industry. J Compos Mater 51: 1265–1273. https://doi.org/10.1177/0021998316671796 doi: 10.1177/0021998316671796
    [14] Roux M, Eguémann N, Dransfeld C, et al. (2017) Thermoplastic carbon fibre-reinforced polymer recycling with electrodynamical fragmentation: From cradle to cradle. J Thermoplast Compos 30: 381–403. https://doi.org/10.1177/0892705715599431 doi: 10.1177/0892705715599431
    [15] Li H, Qin X, He G, et al. (2016) Investigation of chip formation and fracture toughness in orthogonal cutting of UD-CFRP. Int J Adv Manuf Technol 82: 1079–1088. https://doi.org/10.1007/s00170-015-7471-x doi: 10.1007/s00170-015-7471-x
    [16] An Q, Cai C, Cai X, et al. (2019) Experimental investigation on the cutting mechanism and surface generation in orthogonal cutting of UD-CFRP laminates. Compos Struct 230: 111441. https://doi.org/10.1016/j.compstruct.2019.111441 doi: 10.1016/j.compstruct.2019.111441
    [17] Herzog D, Schmidt-Lehr M, Oberlander M, et al. (2016) Laser cutting of carbon fibre reinforced plastics of high thickness. Mater Design 92: 742–749. https://doi.org/10.1016/j.matdes.2015.12.056 doi: 10.1016/j.matdes.2015.12.056
    [18] MM IW, Azmi A, Lee C, et al. (2018) Kerf taper and delamination damage minimization of FRP hybrid composites under abrasive water-jet machining. Int J Adv Manuf Technol 94: 1727–1744. https://doi.org/10.1007/s00170-016-9669-y doi: 10.1007/s00170-016-9669-y
    [19] de Bruijn TA, Vincent GA, van Hattum F (2017) Recycling C/PPS laminates into long fibre thermoplastic composites by low shear mixing. 21st International Conference on Composite Materials 2017, 20–25. Available from: http://www.iccm-central.org/Proceedings/ICCM21proceedings/papers/4122.pdf.
    [20] Vincent GA, de Bruijn TA, Wijskamp S, et al. (2021) Characterisation and improvement of the quality of mixing of recycled thermoplastic composites. Compos Part C-Open 4: 100108. https://doi.org/10.1016/j.jcomc.2021.100108 doi: 10.1016/j.jcomc.2021.100108
    [21] Feraboli P, Cleveland T, Ciccu M, et al. (2010) Defect and damage analysis of advanced discontinuous carbon/epoxy composite materials. Compos Part A-Appl S 41: 888–901. https://doi.org/10.1016/j.compositesa.2010.03.002 doi: 10.1016/j.compositesa.2010.03.002
    [22] Landry B, Hubert P (2015) Experimental study of defect formation during processing of randomly-oriented strand carbon/PEEK composites. Compos Part A-Appl S 77: 301–309. https://doi.org/10.1016/j.compositesa.2015.05.020 doi: 10.1016/j.compositesa.2015.05.020
  • matersci-10-05-044-Supplementary.pdf
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(1909) PDF downloads(116) Cited by(1)

Figures and Tables

Figures(5)  /  Tables(5)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog