• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Enhanced prediction of anisotropic deformation behavior using machine learning with data augmentation

    2024-04-18 13:44:24SujeongByunJinyeongYuSehoCheonSeongHoLeeSungHyukPrkTekyungLee
    Journal of Magnesium and Alloys 2024年1期

    Sujeong Byun ,Jinyeong Yu ,Seho Cheon ,Seong Ho Lee ,Sung Hyuk Prk ,Tekyung Lee,?

    a School of Mechanical Engineering, Pusan National University, Busan 46241, South Korea

    b School of Materials Science and Engineering, Kyungpook National University, Daegu 41566, South Korea

    Abstract Mg alloys possess an inherent plastic anisotropy owing to the selective activation of deformation mechanisms depending on the loading condition.This characteristic results in a diverse range of flow curves that vary with a deformation condition.This study proposes a novel approach for accurately predicting an anisotropic deformation behavior of wrought Mg alloys using machine learning (ML) with data augmentation.The developed model combines four key strategies from data science: learning the entire flow curves,generative adversarial networks (GAN),algorithm-driven hyperparameter tuning,and gated recurrent unit (GRU) architecture.The proposed model,namely GANaided GRU,was extensively evaluated for various predictive scenarios,such as interpolation,extrapolation,and a limited dataset size.The model exhibited significant predictability and improved generalizability for estimating the anisotropic compressive behavior of ZK60 Mg alloys under 11 annealing conditions and for three loading directions.The GAN-aided GRU results were superior to those of previous ML models and constitutive equations.The superior performance was attributed to hyperparameter optimization,GAN-based data augmentation,and the inherent predictivity of the GRU for extrapolation.As a first attempt to employ ML techniques other than artificial neural networks,this study proposes a novel perspective on predicting the anisotropic deformation behaviors of wrought Mg alloys.

    Keywords: Plastic anisotropy;Compression;Annealing;Machine learning;Data augmentation.

    1.Introduction

    Mg alloys have drawn substantial attention in the aerospace and automotive industries because of their high strengthto-weight ratios.They also have potential applications in biomedical and electronic devices because of their biocompatibility,biodegradability,and low density [1–4].However,in addition to their low strength and formability,Mg alloys exhibit anisotropic deformation behavior,wherein their mechanical properties significantly vary with the loading direction.This anisotropy is a result of the unique hexagonal closepacked structure and subsequent grain orientation of Mg,and it leads to the selective activation of deformation mechanisms depending on the loading condition[5–9].Accurate prediction and control of anisotropic deformation behavior are crucial because they can significantly impact the overall performance and durability of Mg alloy components.

    Typically,constitutive analysis is used to estimate the anisotropic deformation behavior of Mg alloys [10–15].This approach correlates physical quantities based on a fundamental physical model and subsequently adjusts them using numerical constants for empirical regression.Although constitutive analysis provides valuable information for characterizing the deformation mechanisms,it often relies on several assumptions,numerical fittings,and even manual adjustments of constants for practical implementation,as discussed in detail below.These inherent limitations of constitutive approach highlight the necessity of an alternative method for predicting the anisotropic properties of wrought Mg alloys.

    Currently,researchers are paying considerable attention to machine learning(ML)techniques because of their functionality with respect to powerful data analysis and predictivity for multivariant nonlinear problems.Numerous studies have used this technique to interpret the thermomechanical properties of Mg alloys [16–19].For example,Sani et al.[16]conducted a comparative study between an artificial neural network(ANN)and constitutive analysis in predicting the high-temperature deformation behavior of cast Mg-Al-Ca alloys.Wang et al.[17]utilized an ANN to improve processing maps describing the hot-deformation behavior of the AZ91 Mg alloy.Yu et al.[18]employed ML techniques to precisely estimate the radical temperature change caused by a short-duration electropulse in ZK60 Mg alloy.Moreover,they compared ML algorithms in terms of both accuracy and prediction resources,wherein the extreme gradient boosting algorithm achieved a 1000-fold acceleration in the learning rate compared to that of the ANN.Similarly,Xu et al.[19]compared a support vector machine algorithm with an ANN to estimate the tensile properties of an AZ31 Mg alloy.

    To the best of our knowledge,only two studies have used ML to predict the severe anisotropy of deformation behavior for wrought Mg alloys.Our previous model [20]was the first to use an ANN for this purpose.The predictability of this model was better than those of other regression methods,including multiple linear regression (MLR) and shallow ANN.Zhang et al.[21]recently developed another ANNbased model that correlated the texture of the AZ31 Mg alloy with its tensile properties.Although this model provides new insights into texture-based estimation,it exhibits high prediction errors up to 20 % for YS and 10 % for the elongation to failure (EL),respectively.

    There is scope for the further improvement of the accuracy of ML predictions from the viewpoint of data science.This study proposes a novel ML-based predictive model specifically designed to predict the anisotropic deformation behavior of wrought Mg alloys.This model was constructed based on four key strategies to enhance model accuracy.First,it learns the data of the entire flow curve rather than the summarized mechanical properties (e.g.,YS and EL) to increase the size of the training dataset.Second,for the same purpose,a data augmentation technique using generative adversarial networks(GAN)is utilized.Third,hyperparameter tuning is determined by an algorithm,instead of the trial-and-error method adopted in previous studies.Lastly,the developed model substitutes an ANN with a recurrent neural network (RNN),expecting a better estimation of the time-series data of the flow curves.The developed model is extensively evaluated under various scenarios,including interpolation,extrapolation,and datasets of limited size to ensure its robustness and reliability.

    2.Data acquisition and preprocessing

    This study investigated the ZK60 Mg alloy with a chemical composition of 5.74 Zn,0.51 Zr,0.02 Si,0.02 Al,0.03 Cu,and balanced Mg (numbers in mass percentage).The cylindrical ingots were produced using induction melting with a diameter of 26 mm and a length of 500 mm.They were homogenized at 673 K for 10 h,air-cooled,and then reheated to the same temperature for a multi-pass caliber rolling.The caliber rolling consisted of six deformation passes with the total area reduction of 84 %.The caliber-rolled rods were subsequently annealed in an electric furnace at 473–673 K for 0.5–3 h,followed by air cooling to room temperature.The detailed rolling procedures are described in the authors’previous studies [7,22,23].

    Compressive deformation properties of ZK60 Mg alloy were obtained using a universal mechanical testing machine at a constant cross-head speed of 1 mm min-1at room temperature.Cylindrical samples were machined from caliberrolled rods that were 6 mm in diameter and 9 mm in height based on the ASTM-E9 standard.The heights of the samples were aligned with either of the following directions to evaluate the anisotropic deformation behaviors: rolling direction(RD),transverse direction (TD),or normal direction (ND).The sample was rotated by 90° per caliber-rolling pass,leading to an alternation between the two coordinate axes on the cross-section.Accordingly,the terms “TD” and “ND” in this work indicate the coordinate system for the final pass of caliber rolling.

    Fig.1 depicts the specific experimental conditions used to construct the training and test datasets for ML.Three flow curves were acquired from the compressive tests of the asrolled alloy along RD,TD,and ND.Nine RD flow curves were obtained under 3×3 annealing conditions,encompassing temperatures of 473,573,and 673 K,each for durations of 0.5,1,and 3 h.Similar conditions were applied to TD and ND flow curves,except for the annealing at 673 K for 1 h,leading to the accumulation of 16 flow curves.In addition,the RD flow curve used in a prior investigation [24]was included in the database;the result was obtained after annealing at 523 K for 1 h.Consequently,the current database comprises 29 flow curves with 94,240 instances.This asymmetric database provided valuable insights into this research,as further discussed in Section 5.3.

    Fig.1.Annealing conditions for experimental data used in the current study.The empty dots indicate the annealing conditions,under which the sample was compressed along only RD.The blue circles indicate the test dataset for interpolation and extrapolation.The datum at zero duration indicates the as-rolled (i.e.,unannealed) condition.

    Each instance involved five features: annealing temperature,annealing period,loading direction,engineering strain,and engineering stress.Three loading directions were allocated to integers ranging from 1 to 3.Subsequently,the entire dataset was standardized using an average and standard deviation to minimize unnecessary interference by the range of features.

    This study designed two types of predictive scenarios (i.e.,interpolation and extrapolation) to extensively evaluate the performance of the developed model.Each scenario has its own test datasets for the RD and ND flow curves,because they exhibited the most distinct differences and were thus the most challenging to predict (Fig.1).The interpolated test dataset corresponded to the flow curves annealed at 573 K for 1 h,enabling a direct comparison between the present method and the previous ANN model [20].The extrapolated test dataset was set as the flow curve annealed at 673 K for 3 h,which was expected to be the most difficult to predict owing to the largest variation in flow behaviors across all experiments.Additionally,the developed ML model was trained with an abridged training dataset (i.e.,10 % of the original dataset) to create a harsh environment for prediction.To achieve this,the original training dataset was equally divided into ten categories,nine of which were excluded from the subsequent ML process.Table 1 summarizes the prediction conditions used in this study.

    Table 1 Predicting conditions and number of instances used in this study.

    3.Machine learning

    3.1. Selection of ML models

    The authors explored a wide variety of ML models in pursuit of the research objective,including MLR,shallow ANN,DNN,random forest (RF),extreme gradient boosting (XGB),light gradient-boosting machine (LGBM),gated recurrent unit(GRU),and long short-term memory (LSTM).The first three were investigated in our previous work [20]with the conclusion favoring DNN due to its superior predictive performance.As a result,DNN was selected as one of three ML models for this study.Ensemble ML models using a decision tree (e.g.,RF,XGB,and LGBM) are effective for a certain problem.However,they exhibited comparable or even lower performance when applied to predicting sequential data in literature[18,25–29],leading to their exclusion from this research.

    RNN represents an advanced ANN architecture designed to perceive sequential data using recurrent connections working as memory [30].LSTM and GRU are the most widely employed RNN models,of which the latter was incorporated into this study due to its faster and superior performance in estimating Al6061 alloy behavior [31].This model resolves the vanishing gradient issue of traditional RNN by adopting gating units to modulate information flow [32].The authors developed an ML model,namely ‘GAN-aided GRU’,by training an optimized GRU model with data augmented through GAN.Hence,this study compares the GAN-aided GRU model with two advanced ML models (i.e.,DNN and GRU) and traditional constitutive approaches.

    3.2. ML framework

    GAN is a deep-learning architecture that generates synthetic data resembling genuine data [33].It is composed of two competing neural networks,denoted as a generator and a discriminator.False samples are created by the generator and then detected by the discriminator.The repeated opposition improves the overall performance of the GAN model,causing the generated data to be indistinguishable from the original data.

    In this study,a GAN was adopted for data augmentation,in which the strain and stress were set as the output variables of the generator and the input variables of the discriminator.Both models have a multilayer perceptron structure with 60 units within a single hidden layer.The leaky rectified linear unit [34]was employed as the activation function for the input and hidden layers,and an exponential linear unit [35]was used for the output layer,based on a preliminary comparative study.Stochastic optimization was performed using the adaptive moment estimation [36].The learning rate was set to be 5×10-4.The dropout method [37]was employed to avoid overfitting.The GAN generated 104instances from the training dataset for each flow curve.The generated data with a strain higher than EL per curve were excluded from the subsequent ML process to increase the prediction accuracy.As a result,the training dataset approximately quadrupled through the data augmentation,with deviations depending on the predicted conditions.The detailed sizes of the training datasets for each condition are listed in Table 1.

    The GRU model employed in this study had four input features: the annealing temperature,annealing period,loading direction,and strain.It also exhibits a single output feature of stress.The hyperbolic tangent and adaptive moment estimation were used as the activation function and optimization algorithm,respectively.Meanwhile,the DNN model for comparison underwent optimization in the authors’ previous study[20].It shared the identical input and output features with the GRU model.The rectified linear unit,mean squared error,and adaptive moment estimation were adopted for the activation function,loss function,and optimization,respectively.The dropout routine was omitted in both models because it had marginal impact on enhancing a model performance in initial assessments.The hidden structure and learning rate for each condition were determined through the hyperband optimization [38].

    Model accuracy was verified in terms of coefficient of determination (R2) and root mean square error (RMSE),defined as follows:

    where,yis an experimental value,is an average value,is an estimation,andnis the total number of instances.The higher model accuracy is indicated by either a higherR2close to 1 or a lower RMSE value.

    Finally,this study implemented Shapley additive explanations (SHAP) [39]for detailed analysis.SHAP is an ML method inspired by game theory,which measures the responsibility of a given feature by comparing changes in the model output.The SHAP analysis yields a quantified contribution of individual features to the prediction,termed the SHAP value.The feature information and predicted results were used for SHAP analysis using the deep explainer module.All code was generated using Python ver.3.9 and implemented in Tensor-Flow ver.2.7.0 [40].

    4.Results

    Fig.2 depicts the flow curves of the caliber-rolled Mg alloy under RD,TD,and ND compressions,which significantly varied.The RD and TD flow curves exhibited a sigmoidal morphology (Fig.2a and b).The strain hardening rate rapidly increased in the early stage of deformation,as characterized by the concave-upward curve,and then decreased on further straining.The RD flow curves showed a larger curvature compared to that of the TD.By contrast,the ND flow curves demonstrated the concave-downward morphology across the entire analysis range,indicating a monotonic reduction of strain hardening rate (Fig.2c).These differences were attributed to the selective activation of deformation mechanisms as discussed in Section 5.1.

    Fig.2.Experimental data used for ML: 29 flow curves of caliber-rolled and annealed ZK60 Mg alloys under compression along (a) RD,(b) TD,and (c)ND.The legends present the annealing temperature and time for each curve.

    Fig.3 shows the comparison between the experimental data and the interpolation prediction based on the two methods(i.e.,with and without data augmentation using the GAN).The GRU model without data augmentation presented an accurate prediction of the RD flow curve with a low RMSE(2.97 MPa),whereas the predictivity was degraded for the ND flow curve with a higher RMSE (16.72 MPa).TheR2value for ND (0.91) is significantly lower than that for RD(0.99).These results propose overfitting towards RD compression when the prediction relies on the GRU alone.Notably,integrating a GAN with a GRU resolves the overfitting issue.This GAN-aided GRU model remarkably improved the predictivity for both RD and ND compressions compared to the GRU alone.In particular,RMSE for ND prediction was considerably reduced from 16.72 to 1.60 MPa,whileR2increased from 0.91 to 0.99 (Fig.3b).This improvement is discussed in more detail along with other predictive approaches in Section 5.2.

    Fig.3.Predicted flow curves in the plastic regime compared with the experimental data obtained from (a) RD and (b) ND compressions.The test dataset corresponds to the interpolating condition (i.e.,annealed at 573 K for 1 h).The numbers in brackets indicate RMSE and R2.

    Fig.4 compares the predictivity of the interpolating conditions with the 90 %-reduced training dataset.The predictivity of the RD flow curve was minimally modified by the GAN integration (Fig.4a).Both models provided precise RD predictions,with a significantR2of 0.99.However,they exhibited different predictive abilities for ND compression (Fig.4b).The GRU alone yielded an appreciable error in the ND estimation at plastic strain values higher than 0.12,as indicated by the arrows.RMSE increased from 9.61 to 14.3 MPa when focusing on the region of high plastic strain values.By contrast,the GAN-aided GRU exhibited a consistently precise estimation across the entire strains for ND curve,as characterized by RMSE of 6.65 MPa andR2of 0.99.These results correlate well with the interpolation using the original dataset,as aforementioned (Fig.3b).It is also noted that the 90 %-reduction in the dataset did not cause a significant deterioration in predictivity under the interpolating condition.The estimating errors remained limited to 12.6 and 12.8 MPa in the RD and ND compressions,respectively,based on the GAN-aided GRU model.The results indicate the remarkable performance of ML-based interpolation,even with a smallsized dataset.

    Fig.4.Predicted flow curves in the plastic regime compared with the experimental data obtained from (a) RD and (b) ND compressions.The test dataset corresponds to the interpolating condition (i.e.,annealed at 573 K for 1 h).The training dataset was reduced by 90 % from the original dataset.The numbers in brackets indicate RMSE and R2.The arrows in Fig.4b indicate the significant deviation made by the GRU alone.

    Fig.5 shows the results of extrapolating prediction.In general,extrapolation is considered more challenging than interpolation because it estimates values beyond the experimental range.Moreover,extrapolation is more susceptible to errors than interpolation,particularly with a nonlinear relationship among the variables,or with a variation in the relationship out of the experimental range.However,no radical deterioration in predictivity was observed when the developed models were adopted for extrapolation.This satisfactory result could be attributed to the improved predictivity of the GRU algorithm with respect to extrapolation [41].For RD estimation,both approaches yielded similar results up to a plastic strain of approximately 0.09(Fig.5a).The GAN-aided GRU model exhibits a marginally larger error with increased strain.However,it significantly enhanced the predictability of the ND estimation (Fig.5b).It eliminated the large errors (up to 24.8 MPa)generated by the GRU alone in the 0–0.10 strain range,as indicated by the arrows.

    Fig.5.Predicted flow curves in the plastic regime compared with the experimental data obtained from (a) RD and (b) ND compressions.The test dataset corresponds to the extrapolating condition (i.e.,annealed at 673 K for 3 h).The numbers in brackets indicate RMSE and R2.The arrows in Fig.5b indicate the significant deviation made by the GRU alone.

    Fig.6 shows the extrapolation results for the reduced training dataset.The predictivity was significantly degraded regardless of data augmentation,owing to the most extreme estimating conditions induced by the combination of data reduction and extrapolation.For RD prediction,the GAN-aided GRU model yielded an increased RMSE (14.2 MPa) and decreasedR2(0.97) compared with the GRU alone (8.49 MPa and 0.99,respectively) (Fig.6a).However,it is not conclusive whether data augmentation was ineffective.The GRU model presented unreasonable fluctuations for the RD prediction,as indicated by the arrows in Fig.6a.The GAN-based data augmentation removed these fluctuations and modified the flow curve to be reasonably smooth,with an error limited to 24.9 MPa,despite a slight increase in the RMSE.In this study,the developed model only failed to estimate the ND extrapolation of the 90 %-reduced training dataset (Fig.6b).These results are discussed in Section 5.3 with the SHAP analysis.Consequently,the GAN-aided GRU model developed in this study provided remarkable predictability for all the investigated conditions except for one.

    Fig.6.Predicted flow curves in the plastic regime compared with the experimental data obtained from (a) RD and (b) ND compressions.The test dataset corresponds to the extrapolating condition (i.e.,annealed at 673 K for 3 h).The training dataset was reduced by 90 % from the original dataset.The numbers in brackets indicate RMSE and R2.The arrows in Fig.6a indicate the unreasonable fluctuations made by the GRU alone.

    5.Discussion

    5.1. Anisotropic deformation behaviors induced by the caliber rolling

    The distinct flow behaviors depicted in Fig.2 stem from the selective activation of deformation mechanisms depending on the loading direction.This topic has been comprehensively discussed elsewhere [7,22,23]with a microstructural characterization,given its broader scope beyond the present study.This subsection provides a concise summary of the underlying origins of anisotropic flow behaviors based on the authors’previous studies [7,22,23].

    A caliber-rolled Mg alloy exhibits a unique basal texture in which the basal poles are split by approximately 30° from the ND to the TD on the cross-section [22].RD compression is perpendicular to most basal poles under such conditions,thereby resulting in the extensive activation of {10-12} extension twinning.This resulted in the formation of a concaveupward morphology of the RD flow curves in the early stage of deformation until extension twinning was depleted.Although TD compression induced a similar onset of {10-12}extension twins,the texture intensity was lower than that of RD compression because of the presence of basal poles against twinning (i.e.,poles close to the TD before compression).This resulted in the TD flow curves having lower curvatures than the RD flow curves.ND compression inhibited {10-12} extension twinning,as most basal poles were aligned with the compressive loading direction,whereas plastic deformation was primarily accommodated by basal slip.This was proven by three factors: (i) the highest Schmid factor for basal slip under the state,(ii) insufficient resolved shear stress for prismatic slip,and (iii) in-grain misorientation axis analysis,confirming that the majority of grains were deformed by basal slip [7].The inhibited {10-12}extension twinning results in a monotonic reduction in the strain-hardening rate of the ND flow curves,as shown in Fig.2c.

    5.2. Superior performance of the developed approach

    The GAN-aided GRU model successfully estimated the anisotropic flow curves under all conditions(Figs.3,4,5,and 6a),except for the ND extrapolation with a severely reduced training dataset (Fig.6b).Fig.7 depicts the parity plots of the various predictive approaches under interpolating conditions for further comparison.The GAN-aided GRU model provided the best predictivity and superior generalizability (i.e.,resistance to overfitting)among the analyzed approaches.This section explains the reasons for the superior performance of the GAN-aided GRU by comparing it to previous predictive approaches.

    Fig.7.Predictions for(a)RD and(b)ND compressions using the GAN-aided GRU,GRU alone,previous ANN model[20],and constitutive equations[10].The numbers in brackets indicate RMSE and R2.

    The previous ANN model [20]provided a good ND prediction with a low RMSE (3.28 MPa) and highR2(0.99),although both values were degraded for RD estimation(15.0 MPa and 0.98).This ANN model exploited hyperparameter optimization using the hyperband algorithm [38]and exhibited better predictivity compared to traditional ML methods,such as multiple linear regression and ANN without hyperparameter optimization [20].This finding supports the effectiveness of algorithm-based hyperparameter tuning,which is one of the four strategies involved in the proposed approach.

    Another strategy explored in this study was the implementation of the RNN architecture (i.e.,GRU) known for its aptitude in handling sequential data.A few recent studies have investigated RNN’s applicability in predicting mechanical properties [42,43].For instance,LSTM model reduced RMSE by 68 %–80 % compared to the modified Johnson–Cook constitutive equation for estimating a high-temperature deformation behavior of Ni-base superalloy [42].However,these results do not assert the superiority of RNN over other ML approaches.Indeed,in a recent study [44],RNN performed less effectively than ANN in estimating the mechanical strengths of a metallic clad.The current results also corroborate the marginal improvement in predictive performance by substituting ANN with RNN,which is evident in both RMSE andR2.The GRU model enhanced the RD estimation compared to ANN,reducing the RMSE from 15.0 to 2.97 MPa.However,this improvement was attained at the expense of ND prediction with RMSE rising from 3.28 to 16.7 MPa.The imbalance in predictivity between the RD and ND remained unresolved in either case.This is clearly indicated by the deviation from the experimental results in their parity plots,i.e.,RD prediction by the ANN (Fig.7a) and ND prediction by the GRU (Fig.7b).In other words,both models exhibited insufficient generalizability,despite decent predictivity.

    The GAN-based data augmentation strategy substantially bolstered the performance of the GAN-aided GRU model,enhancing both predictability and generalizability in comparison to the recurrent connections in the GRU.Expanding the size of the training dataset generally improves a model performance and reduces the expected variance,particularly in nonlinear cases[45].Accordingly,data augmentation has been actively explored to resolve the problem of limited datasets derived from a restricted number of experiments [46,47].The novelty of this study involves leveraging GAN for generating the synthetic flow curves.GAN and its variants have conventionally been utilized to generate synthetic images,rather than sequential mechanical data [48].The mechanical-data augmentation has relied on statistical measures instead.For example,Suh et al.[49]expanded the mechanical-property dataset of Mg alloys by incorporating averages and standard deviations.Although this method achieved favorable predictions under specific conditions,it is unsuitable for the present scenario where acquiring such statistical values is unfeasible.The current study successfully addresses this challenge by pioneering the application of GAN-based data augmentation.The proposed integration of GAN with GRU also sets it apart from a hybrid GAN-GRU model [50].The GRU architecture is implemented for the generator and discriminator of GAN in the hybrid model,whereas the GANaided GRU model innovatively leverages both architectures in parallel.

    The constitutive equation is a conventional approach for interpreting the deformation behavior of a material.For the constitutive analysis of the ND estimation,this study adopted the power-law equation,which has been widely used for continuous hardening behavior and is expressed as follows:

    where,σis the stress,εis the strain,nis the strain-hardening exponent,andKis a numerical constant.This model exhibited high accuracy with a low RMSE (4.23 MPa) and highR2(0.99),comparable to those of the previous ANN model (Fig.7b),suggesting that a constitutive approach can be as effective as the ML method in a specific case.However,this cannot be generalized to the RD and TD cases because Eq.(3) cannot describe a concave-upward curve because of its nature.

    To attain the generalizability,Barnett et al.[10]introduced a constitutive framework incorporating multiple microstructural and mechanical parameters as follows:

    where,Xis the volume fraction of grains favorable for{10-12} twinning,βindicates the extent to slip,fis the fraction of twinned area,σsis the slip stress,σTis the twinning stress,σ0tis the stress required to initiate twinning,σ0sis the stress measured when the basal axis was arranged perpendicular to the tensile direction,ξis the ratio of the stresses in the twinned and untwinned areas,Eis the elastic modulus,andεtis the strain caused by twinning.Although this equation provides valuable insights into the physical correlations among the parameters,it has several limitations in practical use.Some parameters (nandσ0s) should be determined from a separate experiment using a single crystal with an orientation unfavorable for {10-12} twinning.The other parameters were determined by assumptions (εt),numerical fitting (ξ),or manual adjustment (β).It should also be noted that this framework requires twinning kinetics (Xandf) with straining,which must be either frequently measured or assumed using another empirical equation.The lack of supplementary experiments resulted in poor predictability (Fig.7a).Other equations [11–15]have similar issues because of the inherent problems in the constitutive approach.The GAN-aided GRU,which combines the ML approach with data augmentation,can partially replace these methods because of its simpler,highly accurate,and generalized estimation with fewer experiments.

    5.3. Extrapolation with the reduced training dataset

    Upon revisiting Fig.6,the GRU model predicted unreasonable fluctuations,which were improved by integrating the GAN.Fig.8 depicts the estimations compared with the entire set of experimental data to clarify the reasons for the fluctuations.Notably,the RD compression predicted by GRU alone fell within the band composed of the experimental data obtained at 673 K,namely,the kernel training dataset (Fig.8a).Note that abrupt changes in the estimation,as indicated by the arrows in Fig.6a,were perfectly limited within this band.In other words,the model selected a relatively arbitrary value from this band instead of providing an effective prediction owing to the extreme estimating conditions (i.e.,extrapolation with the 90 %-reduced training dataset).The data augmentation with the GAN mitigated such severe prediction conditions by increasing the dataset size,thereby enhancing the predictive performance of the RD flow curve.Furthermore,this improved performance supports the conclusion that poor extrapolation can be attributed to the insufficient size of the training datasets.

    Fig.8.Extrapolating estimations using the 90 %-reduced training dataset shown with the entire experimental data for (a) RD and (b) ND compressions.The annealing conditions in the brackets indicate the kernel training datasets that constitute the band marked by green areas.

    The ND flow curve predicted by the GRU alone was also within the band for which the boundaries were created by other flow curves (Fig.8b).In contrast to the RD estimates,neither of the proposed models provided an effective ND prediction.This difference stemmed from the significantly smaller size of the reduced ND training dataset (1982 instances) compared with the RD case (3593 instances).Such an inherent lack of ND dataset was hardly resolved even by integrating the GAN.These results led to two conclusions.First,integrating GAN with GRU can simultaneously improve predictability and generalizability.Second,a dataset threshold size exists for this improvement,below which the GAN becomes ineffective.

    The extrapolated test dataset corresponded to the flow curve at 673 K for 3 h.The band for RD compression was formed by the flow curves obtained at 673 K for 0.5 h and 1 h (Fig.8a),whereas those for ND compression were defined by the curves at 673 K for 0.5 h and 573 K for 3 h,respectively (Fig.8b).This difference arose from the absence of ND data at 673 K for 1 h,as depicted in Fig.1.These results imply that the annealing temperature contributed more to the present estimations than the annealing duration.

    SHAP provides deeper insights into the priority of features in a black-box ML model.Fig.9 presents the distribution of the SHAP values and their absolute averages with respect to the features.The mean absolute SHAP value of a given feature quantifies its contribution to the predicted output (i.e.,stress),whereas the colors represent the correlation direction(i.e.,positive or negative correlations).In the current prediction,the flow stress correlated with the input features in the following order: strain>temperature ≈direction>duration.Strain exhibited the highest contribution to stress,with a positive correlation,which was understood by the monotonic increase observed in most cases (Fig.2).The annealing temperature showed a negative correlation with the flow stress,although several unusual cases were confirmed,as indicated by the blue dots in the purple area.These unusual cases indicate a high stress after annealing at elevated temperatures owing to a beneficial loading direction,short annealing duration,or both.The loading direction and annealing duration presented an undefined order of feature tendencies,suggesting that the other features contributed more to the estimated flow stress.In particular,the SHAP-based contributions of the annealing temperature and duration are consistent with the above deductions.Therefore,the reduced extrapolation for RD compression was primarily governed by annealing temperatures,whereas that for ND compression was affected by both temperature and duration,owing to insufficient experiments along the ND at a temperature equal to the test dataset.

    Fig.9.Contribution of each feature to the estimated flow stress based on the SHAP analysis: (a) mean absolute SHAP value and (b) distribution of SHAP values per feature.

    5.4. GAN-aided GRU for material design

    The developed model precisely predicted the plastic deformation behaviors within specific conditions of subsequent heat treatment while accounting for the significant plastic anisotropy.These factors exert a crucial effect on the mechanical properties of wrought Mg alloys.In specific,the investigated alloy exhibits considerable variations in mechanical properties (i.e.,up to 37 % decrease in tensile strength and 50 % increase in EL) in the processing window of 373–673 K and 10–360 min [23].The trade-off relationship between strength and ductility collapses at elevated temperature exceeding 473 K owing to alterations in microstructural mechanism [23].The plastic anisotropy is also ascribed to varying deformation mechanisms depending on the loading direction,as discussed in Section 5.1.Traditional models struggle to capture these intricate variations as it generally assumes static or linearly changing mechanisms,such as the constant YS-EL balance or homogeneous twin activation.On the other hand,the GAN-aided GRU can provide the enhanced predictability and generalizability due to the inherent strength of ML in handling the nonlinear correlations among various factors.

    Finally,it is worth highlighting that the developed model can be easily implemented into finite-element analysis (FEA)to improve a material design due to the FEA-friendly characteristic of ML models [51].Researchers typically incorporate multiple constitutive equations to enable an FEA model to accommodate plastic anisotropy.For example,Trzepieci′nski and Gelgele [52]utilized three constitutive equations along RD,TD,and 45° to RD,respectively,to simulate the FEA of a rectangular deep drawing.Similarly,Kim et al.[53]employed two constitutive equations to model a three-point bending of Mg sheet.The proposed approach is expected not only to streamline the experimentation process but also to simplify the mechanical modeling in FEA by adopting the unified ML-based equation adaptable to various loading conditions.A notable illustration of this capability is seen in the work of Zhao et al.[54]where FEA was employed to analyze a metal-rubber isolator using mechanical properties predicted by ANN subroutine.This simulation not only estimated an elastic modulus of raw material,but also precisely predicted the shock response of the metal-rubber isolator with a complex design.

    6.Conclusions

    The current study represents the first application of ML techniques,besides the simple ANN,to accurately predict the severe anisotropy of compressive deformation behaviors for the wrought Mg alloys.The unique split texture formed by caliber rolling and subsequent annealing resulted in different compressive deformation behaviors depending on the loading direction.The GAN-aided GRU model demonstrated the remarkable predictivity in most scenarios,irrespective of the estimation type (i.e.,interpolation and extrapolation)and dataset size.It only failed to the ND extrapolation with a 90 %-reduced training dataset owing to acute data limitations.Furthermore,the developed model exhibited the improved generalizability compared to the previous DNN model or GRU alone.The superior predictive performance is ascribed to hyperparameter optimization,GAN-based data augmentation,and the intrinsic extrapolation capability of the GRU architecture.In particular,GAN-based data augmentation plays a profound role in enhancing the model performance.SHAP analysis clarified the higher priority of the annealing temperature than that of the duration,which explained the different kernel datasets between the RD and ND compressions.The developed model is a promising alternative to traditional constitutive approaches with the limitations of numerous assumptions,numerical fittings,and manual adjustments.The implementation of GAN-aided GRU into FEA is expected to enhance the structural design of Mg components.Overall,this study offers a novel perspective for predicting the anisotropic deformation behavior of wrought Mg alloys exploiting ML techniques with data augmentation.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    CRediT authorship contribution statement

    Sujeong Byun:Conceptualization,Software,Validation,Visualization,Writing– original draft,Writing– review &editing.Jinyeong Yu:Investigation,Methodology,Validation,Writing– review &editing.Seho Cheon:Data curation,Methodology,Writing– review &editing.Seong Ho Lee:Data curation,Methodology,Writing– review &editing.Sung Hyuk Park:Data curation,Resources,Writing–review&editing,Conceptualization.Taekyung Lee:Conceptualization,Methodology,Project administration,Resources,Software,Supervision,Validation,Visualization,Writing–original draft,Writing– review &editing.

    Acknowledgment

    This work was supported by Korea Institute of Energy Technology Evaluation and Planning (KETEP) grant funded by the Korea government (Grant No.20214000000140,Graduate School of Convergence for Clean Energy Integrated Power Generation),by Korea Basic Science Institute (National Research Facilities and Equipment Center)grant funded by the Ministry of Education (2021R1A6C101A449),and by the National Research Foundation of Korea grant funded by the Ministry of Science and ICT (2021R1A2C1095139),Republic of Korea.

    精品人妻熟女av久视频| 日韩av在线大香蕉| 2021天堂中文幕一二区在线观| 日日啪夜夜撸| 欧美性猛交╳xxx乱大交人| 中国国产av一级| 亚洲国产精品成人久久小说 | 赤兔流量卡办理| 亚洲欧美精品自产自拍| 国产精品,欧美在线| 国产成人a∨麻豆精品| 如何舔出高潮| 99热全是精品| 看黄色毛片网站| 欧美激情国产日韩精品一区| 亚洲成人中文字幕在线播放| 久久久久久伊人网av| 国产精品一区二区三区四区久久| 国产精品久久久久久精品电影小说 | 亚洲国产色片| 国产av麻豆久久久久久久| 不卡视频在线观看欧美| 亚洲精华国产精华液的使用体验 | 村上凉子中文字幕在线| 男女啪啪激烈高潮av片| av在线观看视频网站免费| 国产成人精品久久久久久| 久久热精品热| 国产一区二区在线观看日韩| 欧美又色又爽又黄视频| 黄片无遮挡物在线观看| 久久这里只有精品中国| 亚洲自偷自拍三级| 久久久精品欧美日韩精品| 亚洲精品影视一区二区三区av| 1024手机看黄色片| 一级av片app| 欧美xxxx黑人xx丫x性爽| 桃色一区二区三区在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 美女大奶头视频| 我要搜黄色片| 国产精品久久电影中文字幕| 婷婷精品国产亚洲av| 少妇丰满av| av在线播放精品| 女人被狂操c到高潮| 国产亚洲精品av在线| 人妻夜夜爽99麻豆av| 两个人的视频大全免费| 又粗又爽又猛毛片免费看| 国产av不卡久久| 在现免费观看毛片| 亚洲人成网站在线观看播放| 国产精品无大码| 一级毛片久久久久久久久女| 中文字幕免费在线视频6| 人妻制服诱惑在线中文字幕| 久久久成人免费电影| 性插视频无遮挡在线免费观看| 成人亚洲精品av一区二区| 久久久久国产网址| 我的女老师完整版在线观看| 亚洲av免费高清在线观看| 国产成人a区在线观看| 精品熟女少妇av免费看| 成年女人永久免费观看视频| 国产美女午夜福利| 日韩国内少妇激情av| 成人三级黄色视频| 18禁在线无遮挡免费观看视频| 亚洲精品乱码久久久v下载方式| 色综合色国产| 亚洲18禁久久av| 1024手机看黄色片| 国产成人一区二区在线| 国产 一区精品| 国产精品女同一区二区软件| 中出人妻视频一区二区| 欧美激情久久久久久爽电影| 日韩欧美一区二区三区在线观看| 老司机影院成人| 青春草视频在线免费观看| 两个人的视频大全免费| 久久久久久久午夜电影| 99热这里只有是精品在线观看| 亚洲精华国产精华液的使用体验 | 婷婷色av中文字幕| 国语自产精品视频在线第100页| 亚洲欧洲国产日韩| 久久国内精品自在自线图片| 国产午夜精品一二区理论片| 国产亚洲5aaaaa淫片| 成人漫画全彩无遮挡| 青春草国产在线视频 | 国国产精品蜜臀av免费| 亚洲欧美清纯卡通| 一级av片app| 亚洲成a人片在线一区二区| 国产又黄又爽又无遮挡在线| 黄色一级大片看看| 欧美精品一区二区大全| 国产一区二区亚洲精品在线观看| 国产69精品久久久久777片| 国产免费男女视频| 国产精品国产三级国产av玫瑰| 国产老妇伦熟女老妇高清| 又黄又爽又刺激的免费视频.| 亚洲性久久影院| h日本视频在线播放| 26uuu在线亚洲综合色| 六月丁香七月| av专区在线播放| 12—13女人毛片做爰片一| a级毛片免费高清观看在线播放| 秋霞在线观看毛片| 国产白丝娇喘喷水9色精品| www.色视频.com| 午夜免费激情av| 久久久a久久爽久久v久久| 伊人久久精品亚洲午夜| av.在线天堂| 国产一区亚洲一区在线观看| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 国产成人精品婷婷| 欧美3d第一页| 午夜精品一区二区三区免费看| 精品少妇黑人巨大在线播放 | 国内少妇人妻偷人精品xxx网站| 国产一区二区亚洲精品在线观看| 精品99又大又爽又粗少妇毛片| 国产精品乱码一区二三区的特点| 欧美区成人在线视频| 国产一区二区三区av在线 | 寂寞人妻少妇视频99o| 黑人高潮一二区| 在线免费十八禁| 色5月婷婷丁香| 欧美bdsm另类| 久久久国产成人免费| 日韩欧美在线乱码| 男女做爰动态图高潮gif福利片| 午夜免费激情av| 成人午夜精彩视频在线观看| 亚洲乱码一区二区免费版| 亚洲av熟女| 国产 一区 欧美 日韩| 国产午夜福利久久久久久| 亚洲自拍偷在线| 26uuu在线亚洲综合色| 我的女老师完整版在线观看| 床上黄色一级片| 男插女下体视频免费在线播放| 亚洲国产精品久久男人天堂| 久久久久久久久久成人| 男人舔奶头视频| 国产精品一区二区三区四区久久| videossex国产| 亚洲人成网站在线播| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 波多野结衣高清作品| 国产成人影院久久av| 高清毛片免费看| 欧美不卡视频在线免费观看| 看非洲黑人一级黄片| 日本免费a在线| 成人鲁丝片一二三区免费| 国产精品国产高清国产av| 免费电影在线观看免费观看| 韩国av在线不卡| 亚洲色图av天堂| 国产精品一区二区三区四区久久| 精品99又大又爽又粗少妇毛片| 亚洲欧美精品自产自拍| 欧美zozozo另类| 精品久久久噜噜| 亚洲乱码一区二区免费版| 日韩欧美精品免费久久| 小蜜桃在线观看免费完整版高清| 日韩欧美精品v在线| 97在线视频观看| 亚洲国产精品成人久久小说 | 99九九线精品视频在线观看视频| 国产高清三级在线| 亚洲av中文av极速乱| 女人十人毛片免费观看3o分钟| 成人美女网站在线观看视频| 精品国产三级普通话版| 99热这里只有是精品在线观看| 国产精品一区二区在线观看99 | 欧美三级亚洲精品| 男女下面进入的视频免费午夜| 久久久成人免费电影| 日韩欧美在线乱码| 成人高潮视频无遮挡免费网站| 又粗又爽又猛毛片免费看| 久久中文看片网| a级毛片a级免费在线| 可以在线观看的亚洲视频| 99热这里只有是精品50| 亚洲欧美精品自产自拍| 五月玫瑰六月丁香| 久久久色成人| 国产色爽女视频免费观看| 亚洲第一区二区三区不卡| 丰满的人妻完整版| av视频在线观看入口| 中文字幕制服av| 大香蕉久久网| 搡女人真爽免费视频火全软件| 毛片一级片免费看久久久久| 中文字幕制服av| 看十八女毛片水多多多| 婷婷色综合大香蕉| 亚洲aⅴ乱码一区二区在线播放| 内地一区二区视频在线| 国产精品一区二区在线观看99 | 男的添女的下面高潮视频| 免费看美女性在线毛片视频| 亚洲丝袜综合中文字幕| 亚洲综合色惰| 亚洲av不卡在线观看| 国内精品美女久久久久久| 欧美日韩乱码在线| 国产成人午夜福利电影在线观看| 亚洲自偷自拍三级| 精品久久久久久久久久久久久| 国产精品久久久久久亚洲av鲁大| 18禁裸乳无遮挡免费网站照片| 国产伦在线观看视频一区| 丰满人妻一区二区三区视频av| 色视频www国产| av.在线天堂| 两个人视频免费观看高清| 日韩视频在线欧美| 成人国产麻豆网| 美女 人体艺术 gogo| 在线a可以看的网站| 国产精品爽爽va在线观看网站| 中文在线观看免费www的网站| 日本与韩国留学比较| 一本一本综合久久| 一级黄色大片毛片| 免费看美女性在线毛片视频| 久久久欧美国产精品| 亚洲av免费在线观看| 欧美极品一区二区三区四区| 国产一区二区三区av在线 | av免费观看日本| 国产精品一区www在线观看| 老熟妇乱子伦视频在线观看| 综合色av麻豆| 国产激情偷乱视频一区二区| 亚洲自偷自拍三级| 国产又黄又爽又无遮挡在线| 九九久久精品国产亚洲av麻豆| 色视频www国产| 三级国产精品欧美在线观看| 国产成人影院久久av| 国产成人freesex在线| 国产精品三级大全| 亚洲在线观看片| 成人国产麻豆网| 99热这里只有是精品在线观看| 色哟哟哟哟哟哟| 性插视频无遮挡在线免费观看| 亚洲无线观看免费| 国产国拍精品亚洲av在线观看| 内地一区二区视频在线| 国产精品无大码| 国产美女午夜福利| 99久久久亚洲精品蜜臀av| 国产高清不卡午夜福利| 一个人免费在线观看电影| 国产激情偷乱视频一区二区| 亚洲内射少妇av| 日韩欧美国产在线观看| 久久久成人免费电影| 久久草成人影院| 国产精品综合久久久久久久免费| 亚洲欧美日韩高清专用| 黄色日韩在线| 天天躁日日操中文字幕| www.av在线官网国产| 国产精品爽爽va在线观看网站| 久久久精品大字幕| 国产成人91sexporn| 男人狂女人下面高潮的视频| 两个人的视频大全免费| 青春草亚洲视频在线观看| 好男人在线观看高清免费视频| 国产成人福利小说| 欧美xxxx黑人xx丫x性爽| 99久久人妻综合| 成人毛片60女人毛片免费| 亚洲精品国产成人久久av| 干丝袜人妻中文字幕| 色综合色国产| 国产精品免费一区二区三区在线| 亚洲欧美中文字幕日韩二区| 欧美高清成人免费视频www| 免费观看a级毛片全部| 久久午夜福利片| 成人一区二区视频在线观看| 最好的美女福利视频网| 国产精品久久视频播放| 国产精品麻豆人妻色哟哟久久 | 中出人妻视频一区二区| 少妇人妻一区二区三区视频| 欧美丝袜亚洲另类| 男人舔奶头视频| 麻豆成人午夜福利视频| 一进一出抽搐gif免费好疼| 午夜福利视频1000在线观看| av.在线天堂| 日韩高清综合在线| 一级黄色大片毛片| 一级毛片我不卡| 大香蕉久久网| 久久韩国三级中文字幕| 成年版毛片免费区| 91久久精品国产一区二区成人| 国产成人一区二区在线| 午夜亚洲福利在线播放| 久久鲁丝午夜福利片| 超碰av人人做人人爽久久| 国产精品人妻久久久久久| 国产精品三级大全| 欧美最黄视频在线播放免费| 亚洲欧洲日产国产| 久久久久久久久中文| 变态另类丝袜制服| 亚洲丝袜综合中文字幕| 国产精品日韩av在线免费观看| 精品久久久噜噜| 免费不卡的大黄色大毛片视频在线观看 | av视频在线观看入口| 少妇的逼水好多| 久久国产乱子免费精品| av视频在线观看入口| 亚洲不卡免费看| 亚洲成人久久性| 国模一区二区三区四区视频| 在线天堂最新版资源| 久久午夜福利片| 人人妻人人澡人人爽人人夜夜 | 精品国产三级普通话版| 久久久久久久久久久丰满| 看非洲黑人一级黄片| 国产一级毛片在线| 午夜a级毛片| 99热这里只有精品一区| 草草在线视频免费看| 人妻久久中文字幕网| 18禁黄网站禁片免费观看直播| 男人狂女人下面高潮的视频| 欧美bdsm另类| 久久精品国产亚洲av天美| 亚洲欧洲国产日韩| 高清在线视频一区二区三区 | 99在线人妻在线中文字幕| 男女那种视频在线观看| 亚洲18禁久久av| 欧美一级a爱片免费观看看| 国产一区亚洲一区在线观看| 成人一区二区视频在线观看| 亚洲无线在线观看| 国产成人一区二区在线| 岛国在线免费视频观看| 欧美最新免费一区二区三区| 美女cb高潮喷水在线观看| 国产熟女欧美一区二区| 久久国内精品自在自线图片| 一区二区三区高清视频在线| 久99久视频精品免费| 亚洲电影在线观看av| 嫩草影院新地址| 免费一级毛片在线播放高清视频| 免费观看在线日韩| 波多野结衣高清无吗| 搡老妇女老女人老熟妇| 日本在线视频免费播放| 亚洲av中文av极速乱| 少妇猛男粗大的猛烈进出视频 | 久久精品久久久久久久性| 男人舔奶头视频| 日韩一区二区三区影片| 亚洲欧美日韩无卡精品| 午夜福利成人在线免费观看| 日韩欧美一区二区三区在线观看| 精品久久久久久久久久免费视频| 有码 亚洲区| 波多野结衣巨乳人妻| a级毛片免费高清观看在线播放| 欧美高清成人免费视频www| 欧美性猛交╳xxx乱大交人| 欧美人与善性xxx| 岛国在线免费视频观看| 国产老妇女一区| 国产亚洲av嫩草精品影院| 天堂中文最新版在线下载 | 国产精品久久久久久久电影| 日韩亚洲欧美综合| 2022亚洲国产成人精品| 久久精品久久久久久久性| 免费观看a级毛片全部| 国产91av在线免费观看| 丰满人妻一区二区三区视频av| 国内精品久久久久精免费| 精品久久国产蜜桃| 国产伦在线观看视频一区| 国产黄片视频在线免费观看| 网址你懂的国产日韩在线| 精品久久久久久久末码| 精品一区二区三区人妻视频| 中文字幕精品亚洲无线码一区| 久久99蜜桃精品久久| 老司机福利观看| 日韩欧美一区二区三区在线观看| 日韩大尺度精品在线看网址| 乱人视频在线观看| 伦理电影大哥的女人| 一边摸一边抽搐一进一小说| 寂寞人妻少妇视频99o| 亚洲人与动物交配视频| 国产一区二区三区av在线 | 精品久久久久久久久av| 国产亚洲精品av在线| 最近中文字幕高清免费大全6| 大香蕉久久网| 毛片女人毛片| 日日啪夜夜撸| 国产91av在线免费观看| 黑人高潮一二区| 成人午夜精彩视频在线观看| 免费观看人在逋| 白带黄色成豆腐渣| 国产精品嫩草影院av在线观看| 天天躁夜夜躁狠狠久久av| 成年免费大片在线观看| 久久婷婷人人爽人人干人人爱| 久久精品国产亚洲av天美| 日韩三级伦理在线观看| av天堂在线播放| 一级黄片播放器| 精品久久久久久久久亚洲| 久久久久久久久久久丰满| 国产高清视频在线观看网站| 97超碰精品成人国产| 最近中文字幕高清免费大全6| 在线观看66精品国产| 亚洲欧美日韩高清专用| 欧美日本亚洲视频在线播放| 免费黄网站久久成人精品| 亚洲人与动物交配视频| 国产精品久久久久久精品电影小说 | 天堂√8在线中文| 一本精品99久久精品77| 国产精品无大码| 少妇人妻精品综合一区二区 | 国产精品一及| 九九爱精品视频在线观看| 亚洲自拍偷在线| 久久精品夜色国产| 日本黄大片高清| 日韩欧美精品免费久久| 伦理电影大哥的女人| 哪里可以看免费的av片| av黄色大香蕉| 看非洲黑人一级黄片| 日本黄色片子视频| 亚洲精品乱码久久久久久按摩| 国产亚洲精品久久久com| 老熟妇乱子伦视频在线观看| 国产黄色小视频在线观看| 日韩欧美在线乱码| 天堂√8在线中文| 国产成人a区在线观看| 丰满人妻一区二区三区视频av| 丝袜喷水一区| 婷婷六月久久综合丁香| 国产欧美日韩精品一区二区| 两个人视频免费观看高清| 夜夜看夜夜爽夜夜摸| 午夜激情福利司机影院| 在线观看一区二区三区| 午夜精品一区二区三区免费看| 青春草亚洲视频在线观看| h日本视频在线播放| 国产淫片久久久久久久久| 国产av在哪里看| 色尼玛亚洲综合影院| 国产一区二区激情短视频| 亚洲欧美日韩东京热| 久久国内精品自在自线图片| 自拍偷自拍亚洲精品老妇| 亚洲第一电影网av| 男女边吃奶边做爰视频| 99久久九九国产精品国产免费| 精品欧美国产一区二区三| 黄片无遮挡物在线观看| 久久久a久久爽久久v久久| 噜噜噜噜噜久久久久久91| 人人妻人人澡人人爽人人夜夜 | 床上黄色一级片| 国产探花极品一区二区| 春色校园在线视频观看| 国产精品福利在线免费观看| 老师上课跳d突然被开到最大视频| 九九热线精品视视频播放| 亚洲国产精品合色在线| 麻豆成人午夜福利视频| 久久精品国产自在天天线| 亚洲av成人精品一区久久| 久久综合国产亚洲精品| 精品熟女少妇av免费看| 久久久国产成人精品二区| 日日撸夜夜添| 中国美白少妇内射xxxbb| 一级毛片我不卡| 2022亚洲国产成人精品| 一本精品99久久精品77| 久久综合国产亚洲精品| 日本-黄色视频高清免费观看| 少妇熟女欧美另类| 91麻豆精品激情在线观看国产| 亚洲成av人片在线播放无| 成人无遮挡网站| 99国产极品粉嫩在线观看| 级片在线观看| 狂野欧美激情性xxxx在线观看| 亚洲经典国产精华液单| 亚洲人与动物交配视频| 日韩欧美一区二区三区在线观看| 亚洲欧美成人精品一区二区| 一边摸一边抽搐一进一小说| 自拍偷自拍亚洲精品老妇| 老女人水多毛片| 青青草视频在线视频观看| 成年女人看的毛片在线观看| 国产一级毛片在线| 99在线人妻在线中文字幕| 少妇高潮的动态图| 最近的中文字幕免费完整| 国产伦精品一区二区三区视频9| av专区在线播放| 国产精品久久久久久av不卡| 国产精品久久久久久亚洲av鲁大| 中文在线观看免费www的网站| 一级二级三级毛片免费看| 日本黄色视频三级网站网址| 国产精品伦人一区二区| 久久人人爽人人爽人人片va| 免费无遮挡裸体视频| 少妇人妻精品综合一区二区 | 成人亚洲精品av一区二区| 亚洲精品粉嫩美女一区| 国产黄片美女视频| 久久精品国产亚洲网站| 日韩欧美三级三区| 亚洲最大成人中文| 国产精品美女特级片免费视频播放器| 成人国产麻豆网| 人人妻人人看人人澡| 在线观看午夜福利视频| 91aial.com中文字幕在线观看| 99热这里只有是精品50| 精品免费久久久久久久清纯| 一个人看的www免费观看视频| 能在线免费看毛片的网站| ponron亚洲| 久久精品国产99精品国产亚洲性色| 国产一区二区在线观看日韩| 又黄又爽又刺激的免费视频.| 国产精品国产三级国产av玫瑰| 国产色婷婷99| 中国国产av一级| 成人高潮视频无遮挡免费网站| 午夜老司机福利剧场| av在线观看视频网站免费| 成人美女网站在线观看视频| 日韩av在线大香蕉| 亚洲精品亚洲一区二区| 十八禁国产超污无遮挡网站| 日本爱情动作片www.在线观看| 男的添女的下面高潮视频| 国产av一区在线观看免费| 我要搜黄色片| 午夜福利高清视频| 老师上课跳d突然被开到最大视频| 国产黄色小视频在线观看| 精品人妻偷拍中文字幕| 中文在线观看免费www的网站| 嫩草影院新地址| 日韩一本色道免费dvd| 亚洲天堂国产精品一区在线| а√天堂www在线а√下载| 成人毛片a级毛片在线播放| 尤物成人国产欧美一区二区三区| 综合色丁香网| ponron亚洲| 国产综合懂色| 亚洲av第一区精品v没综合| 又粗又爽又猛毛片免费看| 真实男女啪啪啪动态图| 日韩av在线大香蕉| 2021天堂中文幕一二区在线观| 99热精品在线国产| 最后的刺客免费高清国语| 国产一区二区三区av在线 | 丝袜美腿在线中文| 国产亚洲5aaaaa淫片| 99热精品在线国产|