As a connection between the process and the circuit design, the device model is greatly desired for emerging devices, such as the double-gate MOSFET. Time efficiency is one of the most important requirements for devic...As a connection between the process and the circuit design, the device model is greatly desired for emerging devices, such as the double-gate MOSFET. Time efficiency is one of the most important requirements for device modeling. In this paper, an improvement to the computational efficiency of the drain current model for double-gate MOSFETs is extended, and different calculation methods are compared and discussed. The results show that the calculation speed of the improved model is substantially enhanced. A two-dimensional device simulation is performed to verify the improved model. Furthermore, the model is implemented into the HSPICE circuit simulator in Verilog-A for practical application.展开更多
Hyperparameter tuning is a key step in developing high-performing machine learning models, but searching large hyperparameter spaces requires extensive computation using standard sequential methods. This work analyzes...Hyperparameter tuning is a key step in developing high-performing machine learning models, but searching large hyperparameter spaces requires extensive computation using standard sequential methods. This work analyzes the performance gains from parallel versus sequential hyperparameter optimization. Using scikit-learn’s Randomized SearchCV, this project tuned a Random Forest classifier for fake news detection via randomized grid search. Setting n_jobs to -1 enabled full parallelization across CPU cores. Results show the parallel implementation achieved over 5× faster CPU times and 3× faster total run times compared to sequential tuning. However, test accuracy slightly dropped from 99.26% sequentially to 99.15% with parallelism, indicating a trade-off between evaluation efficiency and model performance. Still, the significant computational gains allow more extensive hyperparameter exploration within reasonable timeframes, outweighing the small accuracy decrease. Further analysis could better quantify this trade-off across different models, tuning techniques, tasks, and hardware.展开更多
Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy...Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time(TET) decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method(CDM), the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ(λ = 4) method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.展开更多
It is important to calculate the reachable domain(RD)of the manned lunar mission to evaluate whether a lunar landing site could be reached by the spacecraft. In this paper, the RD of free return orbits is quickly eval...It is important to calculate the reachable domain(RD)of the manned lunar mission to evaluate whether a lunar landing site could be reached by the spacecraft. In this paper, the RD of free return orbits is quickly evaluated and calculated via the classification and regression neural networks. An efficient databasegeneration method is developed for obtaining eight types of free return orbits and then the RD is defined by the orbit’s inclination and right ascension of ascending node(RAAN) at the perilune. A classify neural network and a regression network are trained respectively. The former is built for classifying the type of the RD, and the latter is built for calculating the inclination and RAAN of the RD. The simulation results show that two neural networks are well trained. The classification model has an accuracy of more than 99% and the mean square error of the regression model is less than 0.01°on the test set. Moreover, a serial strategy is proposed to combine the two surrogate models and a recognition tool is built to evaluate whether a lunar site could be reached. The proposed deep learning method shows the superiority in computation efficiency compared with the traditional double two-body model.展开更多
The constant bubble size modeling approach(CBSM)and variable bubble size modeling approach(VBSM)are frequently employed in Eulerian–Eulerian simulation of bubble columns.However,the accuracy of CBSM is limited while ...The constant bubble size modeling approach(CBSM)and variable bubble size modeling approach(VBSM)are frequently employed in Eulerian–Eulerian simulation of bubble columns.However,the accuracy of CBSM is limited while the computational efficiency of VBSM needs to be improved.This work aims to develop method for bubble size modeling which has high computational efficiency and accuracy in the simulation of bubble columns.The distribution of bubble sizes is represented by a series of discrete points,and the percentage of bubbles with various sizes at gas inlet is determined by the results of computational fluid dynamics(CFD)–population balance model(PBM)simulations,whereas the influence of bubble coalescence and breakup is neglected.The simulated results of a 0.15 m diameter bubble column suggest that the developed method has high computational speed and can achieve similar accuracy as CFD–PBM modeling.Furthermore,the convergence issues caused by solving population balance equations are addressed.展开更多
Piles are long, slender structural elements used to transfer the loads from the superstructure through weak strata onto stiffer soils or rocks. For driven piles, the impact of the piling hammer induces compression and...Piles are long, slender structural elements used to transfer the loads from the superstructure through weak strata onto stiffer soils or rocks. For driven piles, the impact of the piling hammer induces compression and tension stresses in the piles. Hence, an important design consideration is to check that the strength of the pile is sufficient to resist the stresses caused by the impact of the pile hammer. Due to its complexity, pile drivability lacks a precise analytical solution with regard to the phenomena involved.In situations where measured data or numerical hypothetical results are available, neural networks stand out in mapping the nonlinear interactions and relationships between the system’s predictors and dependent responses. In addition, unlike most computational tools, no mathematical relationship assumption between the dependent and independent variables has to be made. Nevertheless, neural networks have been criticized for their long trial-and-error training process since the optimal configuration is not known a priori. This paper investigates the use of a fairly simple nonparametric regression algorithm known as multivariate adaptive regression splines(MARS), as an alternative to neural networks, to approximate the relationship between the inputs and dependent response, and to mathematically interpret the relationship between the various parameters. In this paper, the Back propagation neural network(BPNN) and MARS models are developed for assessing pile drivability in relation to the prediction of the Maximum compressive stresses(MCS), Maximum tensile stresses(MTS), and Blow per foot(BPF). A database of more than four thousand piles is utilized for model development and comparative performance between BPNN and MARS predictions.展开更多
The large and complex structures are divided into hundreds of thousands or millions degrees of freedom(DOF) when they are calculated which will spend a lot of time and the efficiency will be extremely low. The class...The large and complex structures are divided into hundreds of thousands or millions degrees of freedom(DOF) when they are calculated which will spend a lot of time and the efficiency will be extremely low. The classical component modal synthesis method (CMSM) are used extensively, but for many structures in the engineering of high-rise buildings, aerospace systemic engineerings, marine oil platforms etc, a large amount of calculation is still needed. An improved hybrid interface substructural component modal synthesis method(HISCMSM) is proposed. The parametric model of the mistuned blisk is built by the improved HISCMSM. The double coordinating conditions of the displacement and the force are introduced to ensure the computational accuracy. Compared with the overall structure finite element model method(FEMM), the computational time is shortened by23.86%–31.56%and the modal deviation is 0.002%–0.157% which meets the requirement of the computational accuracy. It is faster 4.46%–10.57% than the classical HISCMSM. So the improved HISCMSM is better than the classical HISCMSM and the overall structure FEMM. Meanwhile, the frequency and the modal shape are researched, considering the factors including rotational speed, gas temperature and geometry size. The strong localization phenomenon of the modal shape’s the maximum displacement and the maximum stress is observed in the second frequency band and it is the most sensitive in the frequency veering. But the localization phenomenon is relatively weak in 1st and the 3d frequency band. The localization of the modal shape is more serious under the condition of the geometric dimensioning mistuned. An improved HISCMSM is proposed, the computational efficiency of the mistuned blisk can be increased observably by this method.展开更多
The wave-based method (WBM) has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model...The wave-based method (WBM) has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model to converge. Numerical results show that the prediction tech- nique based on WBM is with higher accuracy and smaller computational effort than the one on FEM, which implies that this new technique on WBM can be applied to higher-frequency range.展开更多
In this paper, a new spatial coherence model of seismic ground motions is proposed by a fitting procedure. The analytical expressions of modal combination (correlation) coefficients of structural response are develo...In this paper, a new spatial coherence model of seismic ground motions is proposed by a fitting procedure. The analytical expressions of modal combination (correlation) coefficients of structural response are developed for multi-support seismic excitations. The coefficients from both the numerical integration and analytical solutions are compared to verify the accuracy of the solutions. It is shown that the analytical expressions of numerical modal combination coefficients are of high accuracy. The results of random responses of an example bridge show that the analytical modal combination coefficients developed in this paper are accurate enough to meet the requirements needed in practice. In addition, the computational efficiency of the analytical solutions of the modal combination coefficients is demonstrated by the response computation of the example bridge. It is found that the time required for the structural response analysis by using the analytical modal combination coefficients is less than 1/20 of that using numerical integral methods.展开更多
This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,mult...This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,multi-patch structure is not easy to achieve high continuity between neighboring patches which will reduce the advantage of isogeometric analysis in a sense.The proposed method can achieve high continuity at surface stitching region with low geometric error,and this technique exploits constructing the approximate surface with several control points are from original surfaces,which guarantees the local feature of the surface can be well-preserved with high precision.With the proposed approximating method,isogeometric analysis results using the new single-patch can be obtained efficiently compared with the original multi-patch structure.Several examples are presented to illustrate the effectiveness,accuracy and efficiency of the proposed method.展开更多
Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes w...Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark's time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.展开更多
In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Conver...In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.展开更多
Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, th...Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.展开更多
Relative dispersion ratio(RDR)can be used to quantify the deviation behavior of a water parcel’s trajectory caused by a disturbance in a hydrodynamic system.It can be calculated by using a standard method for determi...Relative dispersion ratio(RDR)can be used to quantify the deviation behavior of a water parcel’s trajectory caused by a disturbance in a hydrodynamic system.It can be calculated by using a standard method for determining relative dispersion(RD),which accounts for the growth of the deviation of a cluster of particles from a specific initial time.However,the standard method for computing RD is time consuming.It involves numerous computations on tracing many water parcels.In this study,a new method based on the adjoint method is proposed to acquire a series of RDR fields in one round of tracing.Through this method,the continuous variation in the RDR corresponding to a time series of the disturbance time t can be obtained.The consistency and efficiency of the new method are compared with those of the standard method by applying it to a double-gyre flow and an unsteady Arnold-Beltrami-Childress flow field.Results show that the two methods have good consistency in a finite time span.The new method has a notable speedup for evaluating the RDR at multiple t.展开更多
This work investigates the machining temperatures of ultra-fine-grained titanium(UFG Ti),prepared by equal channel angular extrusion,through analytical modeling.UFG Ti has great usefulness in biomedical applications b...This work investigates the machining temperatures of ultra-fine-grained titanium(UFG Ti),prepared by equal channel angular extrusion,through analytical modeling.UFG Ti has great usefulness in biomedical applications because of its high mechanical strength,sufficient manufacturability,and high biocompatibility.The temperatures were predicted using a physics-based predictive model based on material constitutive relation and mechanics of the orthogonal cutting process.The minimization between the stress calculated using Johnson–Cook constitutive model and the same stress calculated using mechanics model yields the estimation of machining temperatures at two deformation zones.Good agreements are observed upon validation to the values reported in the literature.The machinability of UFG Ti is investigated by comparing its machining temperature to that of Ti–6Al–4V alloy under the same cutting conditions.Significantly lower temperatures are observed in machining UFG Ti.The computational efficiency of the presented model is investigated by comparing its average computational time(~0.5 s)to that of a widely used modified chip formation model(8900 s)with comparable prediction accuracy.This work extends the applicability of the presented temperature model to a broader class of materials,specifically ultra-fine-grained metals.The high computational efficiency allows the in situ temperature prediction and optimization of temperature condition with process parameters planning.展开更多
The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a ...The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a need for high resolution. To resolve these difficulties, a multiple nested lattice Boltzmann method(MNLBM) was developed, which contains several overlapped layers with different resolutions in the computational domain. The data transference of flow field on two layers is accomplished by a Filippova procedure which is proved to satisfy the continuity of mass, momentum, and stresses across the interface. The proposed method is based on the standard lattice Boltzmann method, so it is easily performed.By numerical investigation, the result of present method has been agreed with that of literature, but the computation efficiency is higher than the standard lattice Boltzmann method.展开更多
In order to investigate the applicability and performance of the Lattice Boltzmann Method( LBM) and the Monte Carlo Method( MCM) to simulate phonon heat transfer,a problem of phonon heat transfer in square geometry st...In order to investigate the applicability and performance of the Lattice Boltzmann Method( LBM) and the Monte Carlo Method( MCM) to simulate phonon heat transfer,a problem of phonon heat transfer in square geometry structures of silicon dioxide was taken as an example to compare the calculational results and analyze computational efficiency of the two methods. Moreover,this work analyzed the numerical stability for two methods. The results indicate that the MCM takes much more computation time than the LBM in the same condition. In addition,the results of the two methods have a good agreement in diffusive and diffusive-ballistic domain for investigating the phonon heat transfer. So they can be used to verify each other when the experiments of energy transport in these domains faces difficulty. In ballistic domain,duo to the random error,the temperature distribution curve from MCM is fluctuant.展开更多
The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In part...The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical computations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth subsonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the prediction of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higherorder harmonic contents to the local pressure fluctuations,a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method.The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.展开更多
Based on Neumman series and epsilon-algorithm, an efficient computation for dynamic responses of systems with arbitrary time-varying characteristics is investigated. Avoiding the calculation for the inverses of the eq...Based on Neumman series and epsilon-algorithm, an efficient computation for dynamic responses of systems with arbitrary time-varying characteristics is investigated. Avoiding the calculation for the inverses of the equivalent stiffness matrices in each time step, the computation effort of the proposed method is reduced compared with the full analysis of Newmark method. The validity and applications of the proposed method are illustrated by a 4-DOF spring-mass system with periodical time-varying stiffness properties and a truss structure with arbitrary time-varying lumped mass. It shows that good approximate results can be obtained by the proposed method compared with the responses obtained by the full analysis of Newmark method.展开更多
基金Project supported by the National Natural Science Foundation of China (Grant No.60876027)the National Science Foundation for Distinguished Young Scholars of China (Grant No.60925015)+1 种基金the National Basic Research Program of China (Grant No.2011CBA00600)the Fundamental Research Project of Shenzhen Science & Technology Foundation,China (Grant No.JC200903160353A)
文摘As a connection between the process and the circuit design, the device model is greatly desired for emerging devices, such as the double-gate MOSFET. Time efficiency is one of the most important requirements for device modeling. In this paper, an improvement to the computational efficiency of the drain current model for double-gate MOSFETs is extended, and different calculation methods are compared and discussed. The results show that the calculation speed of the improved model is substantially enhanced. A two-dimensional device simulation is performed to verify the improved model. Furthermore, the model is implemented into the HSPICE circuit simulator in Verilog-A for practical application.
文摘Hyperparameter tuning is a key step in developing high-performing machine learning models, but searching large hyperparameter spaces requires extensive computation using standard sequential methods. This work analyzes the performance gains from parallel versus sequential hyperparameter optimization. Using scikit-learn’s Randomized SearchCV, this project tuned a Random Forest classifier for fake news detection via randomized grid search. Setting n_jobs to -1 enabled full parallelization across CPU cores. Results show the parallel implementation achieved over 5× faster CPU times and 3× faster total run times compared to sequential tuning. However, test accuracy slightly dropped from 99.26% sequentially to 99.15% with parallelism, indicating a trade-off between evaluation efficiency and model performance. Still, the significant computational gains allow more extensive hyperparameter exploration within reasonable timeframes, outweighing the small accuracy decrease. Further analysis could better quantify this trade-off across different models, tuning techniques, tasks, and hardware.
基金National Natural Science Foundation of China under Grant Nos.51639006 and 51725901
文摘Finite element(FE) is a powerful tool and has been applied by investigators to real-time hybrid simulations(RTHSs). This study focuses on the computational efficiency, including the computational time and accuracy, of numerical integrations in solving FE numerical substructure in RTHSs. First, sparse matrix storage schemes are adopted to decrease the computational time of FE numerical substructure. In this way, the task execution time(TET) decreases such that the scale of the numerical substructure model increases. Subsequently, several commonly used explicit numerical integration algorithms, including the central difference method(CDM), the Newmark explicit method, the Chang method and the Gui-λ method, are comprehensively compared to evaluate their computational time in solving FE numerical substructure. CDM is better than the other explicit integration algorithms when the damping matrix is diagonal, while the Gui-λ(λ = 4) method is advantageous when the damping matrix is non-diagonal. Finally, the effect of time delay on the computational accuracy of RTHSs is investigated by simulating structure-foundation systems. Simulation results show that the influences of time delay on the displacement response become obvious with the mass ratio increasing, and delay compensation methods may reduce the relative error of the displacement peak value to less than 5% even under the large time-step and large time delay.
基金supported by the National Natural Science Foundation of China (12072365)the Natural Science Foundation of Hunan Province of China (2020JJ4657)。
文摘It is important to calculate the reachable domain(RD)of the manned lunar mission to evaluate whether a lunar landing site could be reached by the spacecraft. In this paper, the RD of free return orbits is quickly evaluated and calculated via the classification and regression neural networks. An efficient databasegeneration method is developed for obtaining eight types of free return orbits and then the RD is defined by the orbit’s inclination and right ascension of ascending node(RAAN) at the perilune. A classify neural network and a regression network are trained respectively. The former is built for classifying the type of the RD, and the latter is built for calculating the inclination and RAAN of the RD. The simulation results show that two neural networks are well trained. The classification model has an accuracy of more than 99% and the mean square error of the regression model is less than 0.01°on the test set. Moreover, a serial strategy is proposed to combine the two surrogate models and a recognition tool is built to evaluate whether a lunar site could be reached. The proposed deep learning method shows the superiority in computation efficiency compared with the traditional double two-body model.
基金the National Natural Science Foundation of China(21625603)for supporting this work。
文摘The constant bubble size modeling approach(CBSM)and variable bubble size modeling approach(VBSM)are frequently employed in Eulerian–Eulerian simulation of bubble columns.However,the accuracy of CBSM is limited while the computational efficiency of VBSM needs to be improved.This work aims to develop method for bubble size modeling which has high computational efficiency and accuracy in the simulation of bubble columns.The distribution of bubble sizes is represented by a series of discrete points,and the percentage of bubbles with various sizes at gas inlet is determined by the results of computational fluid dynamics(CFD)–population balance model(PBM)simulations,whereas the influence of bubble coalescence and breakup is neglected.The simulated results of a 0.15 m diameter bubble column suggest that the developed method has high computational speed and can achieve similar accuracy as CFD–PBM modeling.Furthermore,the convergence issues caused by solving population balance equations are addressed.
文摘Piles are long, slender structural elements used to transfer the loads from the superstructure through weak strata onto stiffer soils or rocks. For driven piles, the impact of the piling hammer induces compression and tension stresses in the piles. Hence, an important design consideration is to check that the strength of the pile is sufficient to resist the stresses caused by the impact of the pile hammer. Due to its complexity, pile drivability lacks a precise analytical solution with regard to the phenomena involved.In situations where measured data or numerical hypothetical results are available, neural networks stand out in mapping the nonlinear interactions and relationships between the system’s predictors and dependent responses. In addition, unlike most computational tools, no mathematical relationship assumption between the dependent and independent variables has to be made. Nevertheless, neural networks have been criticized for their long trial-and-error training process since the optimal configuration is not known a priori. This paper investigates the use of a fairly simple nonparametric regression algorithm known as multivariate adaptive regression splines(MARS), as an alternative to neural networks, to approximate the relationship between the inputs and dependent response, and to mathematically interpret the relationship between the various parameters. In this paper, the Back propagation neural network(BPNN) and MARS models are developed for assessing pile drivability in relation to the prediction of the Maximum compressive stresses(MCS), Maximum tensile stresses(MTS), and Blow per foot(BPF). A database of more than four thousand piles is utilized for model development and comparative performance between BPNN and MARS predictions.
基金Supported by National Natural Science Foundation of China (Grant Nos.51375032,51335003)
文摘The large and complex structures are divided into hundreds of thousands or millions degrees of freedom(DOF) when they are calculated which will spend a lot of time and the efficiency will be extremely low. The classical component modal synthesis method (CMSM) are used extensively, but for many structures in the engineering of high-rise buildings, aerospace systemic engineerings, marine oil platforms etc, a large amount of calculation is still needed. An improved hybrid interface substructural component modal synthesis method(HISCMSM) is proposed. The parametric model of the mistuned blisk is built by the improved HISCMSM. The double coordinating conditions of the displacement and the force are introduced to ensure the computational accuracy. Compared with the overall structure finite element model method(FEMM), the computational time is shortened by23.86%–31.56%and the modal deviation is 0.002%–0.157% which meets the requirement of the computational accuracy. It is faster 4.46%–10.57% than the classical HISCMSM. So the improved HISCMSM is better than the classical HISCMSM and the overall structure FEMM. Meanwhile, the frequency and the modal shape are researched, considering the factors including rotational speed, gas temperature and geometry size. The strong localization phenomenon of the modal shape’s the maximum displacement and the maximum stress is observed in the second frequency band and it is the most sensitive in the frequency veering. But the localization phenomenon is relatively weak in 1st and the 3d frequency band. The localization of the modal shape is more serious under the condition of the geometric dimensioning mistuned. An improved HISCMSM is proposed, the computational efficiency of the mistuned blisk can be increased observably by this method.
基金Project supported by the National Natural Science Foundation of China (No.10472035).
文摘The wave-based method (WBM) has been applied for the prediction of mid-frequency vibrations of fiat plates. The scaling factors, Gauss point selection rule and truncation rule are introduced to insure the wave model to converge. Numerical results show that the prediction tech- nique based on WBM is with higher accuracy and smaller computational effort than the one on FEM, which implies that this new technique on WBM can be applied to higher-frequency range.
基金National Natural Science Foundation of China Under Grant No. 50478112
文摘In this paper, a new spatial coherence model of seismic ground motions is proposed by a fitting procedure. The analytical expressions of modal combination (correlation) coefficients of structural response are developed for multi-support seismic excitations. The coefficients from both the numerical integration and analytical solutions are compared to verify the accuracy of the solutions. It is shown that the analytical expressions of numerical modal combination coefficients are of high accuracy. The results of random responses of an example bridge show that the analytical modal combination coefficients developed in this paper are accurate enough to meet the requirements needed in practice. In addition, the computational efficiency of the analytical solutions of the modal combination coefficients is demonstrated by the response computation of the example bridge. It is found that the time required for the structural response analysis by using the analytical modal combination coefficients is less than 1/20 of that using numerical integral methods.
基金This research was supported by the National Nature Science Foundation of China under Grant Nos.61602138,61772163 and 61761136010the NSFC-Zhejiang Joint Fund for the Integration of Industrialization and Informatization(Grant No.U1909210)Zhejiang Provincial Science and Technology Program in China(2018C01030).
文摘This paper presents a novel approximating method to construct highprecision single-patch representation of B-spline surface from a multi-patch representation for isogeometric applications.In isogeometric analysis,multi-patch structure is not easy to achieve high continuity between neighboring patches which will reduce the advantage of isogeometric analysis in a sense.The proposed method can achieve high continuity at surface stitching region with low geometric error,and this technique exploits constructing the approximate surface with several control points are from original surfaces,which guarantees the local feature of the surface can be well-preserved with high precision.With the proposed approximating method,isogeometric analysis results using the new single-patch can be obtained efficiently compared with the original multi-patch structure.Several examples are presented to illustrate the effectiveness,accuracy and efficiency of the proposed method.
基金supported by the National Natural Science Founda-tion of China(11272118)Open Found of State Key Laboratory of Explosion Science and Technology(KFJJ12-5M)
文摘Condensation technique of degree of freedom is first proposed to improve the computational efficiency of meshfree method with Galerkin weak form for elastic dy- namic analysis. In the present method, scattered nodes with- out connectivity are divided into several subsets by cells with arbitrary shape. Local discrete equation is established over each cell by using moving Kriging interpolation, in which the nodes that located in the cell are used for approxima- tion. Then local discrete equations can be simplified by con- densation of degree of freedom, which transfers equations of inner nodes to equations of boundary nodes based on cells. The global dynamic system equations are obtained by as- sembling all local discrete equations and are solved by using the standard implicit Newmark's time integration scheme. In the scheme of present method, the calculation of each cell is carried out by meshfree method, and local search is imple- mented in interpolation. Numerical examples show that the present method has high computational efficiency and good accuracy in solving elastic dynamic problems.
基金the Natural Science Foundation of China(Grant Nos.61673169,11301127,11701176,11626101,and 11601485)The Natural Science Foundation of Huzhou City(Grant No.2018YZ07).
文摘In this article,we construct the most powerful family of simultaneous iterative method with global convergence behavior among all the existing methods in literature for finding all roots of non-linear equations.Convergence analysis proved that the order of convergence of the family of derivative free simultaneous iterative method is nine.Our main aim is to check out the most regularly used simultaneous iterative methods for finding all roots of non-linear equations by studying their dynamical planes,numerical experiments and CPU time-methodology.Dynamical planes of iterative methods are drawn by using MATLAB for the comparison of global convergence properties of simultaneous iterative methods.Convergence behavior of the higher order simultaneous iterative methods are also illustrated by residual graph obtained from some numerical test examples.Numerical test examples,dynamical behavior and computational efficiency are provided to present the performance and dominant efficiency of the newly constructed derivative free family of simultaneous iterative method over existing higher order simultaneous methods in literature.
基金supported by the National Natural Science Foundation of China (Grant No. U1811464)。
文摘Construction of high-order difference schemes based on Taylor series expansion has long been a hot topic in computational mathematics, while its application in comprehensive weather models is still very rare. Here, the properties of high-order finite difference schemes are studied based on idealized numerical testing, for the purpose of their application in the Global/Regional Assimilation and Prediction System(GRAPES) model. It is found that the pros and cons due to grid staggering choices diminish with higher-order schemes based on linearized analysis of the one-dimensional gravity wave equation. The improvement of higher-order difference schemes is still obvious for the mesh with smooth varied grid distance. The results of discontinuous square wave testing also exhibits the superiority of high-order schemes. For a model grid with severe non-uniformity and non-orthogonality, the advantage of high-order difference schemes is inapparent, as shown by the results of two-dimensional idealized advection tests under a terrain-following coordinate. In addition, the increase in computational expense caused by high-order schemes can be avoided by the precondition technique used in the GRAPES model. In general, a high-order finite difference scheme is a preferable choice for the tropical regional GRAPES model with a quasi-uniform and quasi-orthogonal grid mesh.
文摘Relative dispersion ratio(RDR)can be used to quantify the deviation behavior of a water parcel’s trajectory caused by a disturbance in a hydrodynamic system.It can be calculated by using a standard method for determining relative dispersion(RD),which accounts for the growth of the deviation of a cluster of particles from a specific initial time.However,the standard method for computing RD is time consuming.It involves numerous computations on tracing many water parcels.In this study,a new method based on the adjoint method is proposed to acquire a series of RDR fields in one round of tracing.Through this method,the continuous variation in the RDR corresponding to a time series of the disturbance time t can be obtained.The consistency and efficiency of the new method are compared with those of the standard method by applying it to a double-gyre flow and an unsteady Arnold-Beltrami-Childress flow field.Results show that the two methods have good consistency in a finite time span.The new method has a notable speedup for evaluating the RDR at multiple t.
文摘This work investigates the machining temperatures of ultra-fine-grained titanium(UFG Ti),prepared by equal channel angular extrusion,through analytical modeling.UFG Ti has great usefulness in biomedical applications because of its high mechanical strength,sufficient manufacturability,and high biocompatibility.The temperatures were predicted using a physics-based predictive model based on material constitutive relation and mechanics of the orthogonal cutting process.The minimization between the stress calculated using Johnson–Cook constitutive model and the same stress calculated using mechanics model yields the estimation of machining temperatures at two deformation zones.Good agreements are observed upon validation to the values reported in the literature.The machinability of UFG Ti is investigated by comparing its machining temperature to that of Ti–6Al–4V alloy under the same cutting conditions.Significantly lower temperatures are observed in machining UFG Ti.The computational efficiency of the presented model is investigated by comparing its average computational time(~0.5 s)to that of a widely used modified chip formation model(8900 s)with comparable prediction accuracy.This work extends the applicability of the presented temperature model to a broader class of materials,specifically ultra-fine-grained metals.The high computational efficiency allows the in situ temperature prediction and optimization of temperature condition with process parameters planning.
文摘The standard lattice Boltzmann method utilizes uniform grids to maintain a compact computational procedure. However, it is often less efficient to perform hydrodynamic and aerodynamic flow simulations when there is a need for high resolution. To resolve these difficulties, a multiple nested lattice Boltzmann method(MNLBM) was developed, which contains several overlapped layers with different resolutions in the computational domain. The data transference of flow field on two layers is accomplished by a Filippova procedure which is proved to satisfy the continuity of mass, momentum, and stresses across the interface. The proposed method is based on the standard lattice Boltzmann method, so it is easily performed.By numerical investigation, the result of present method has been agreed with that of literature, but the computation efficiency is higher than the standard lattice Boltzmann method.
基金Sponsored by the National Natural Science Foundation of China(Grant No.51176038)the Foundation for Innovative Research Groups of the National Natural Science Foundation of China(Grant No.51121004)
文摘In order to investigate the applicability and performance of the Lattice Boltzmann Method( LBM) and the Monte Carlo Method( MCM) to simulate phonon heat transfer,a problem of phonon heat transfer in square geometry structures of silicon dioxide was taken as an example to compare the calculational results and analyze computational efficiency of the two methods. Moreover,this work analyzed the numerical stability for two methods. The results indicate that the MCM takes much more computation time than the LBM in the same condition. In addition,the results of the two methods have a good agreement in diffusive and diffusive-ballistic domain for investigating the phonon heat transfer. So they can be used to verify each other when the experiments of energy transport in these domains faces difficulty. In ballistic domain,duo to the random error,the temperature distribution curve from MCM is fluctuant.
基金supported by the State Scholarship Fund of the China Scholarship Council (Grant 2009629129)
文摘The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical computations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth subsonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the prediction of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higherorder harmonic contents to the local pressure fluctuations,a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method.The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.
基金supported by the Foundation of the Science and Technology of Jilin Province (20070541)985-Automotive Engineering of Jilin University and Innovation Fund for 985 Engineering of Jilin University (20080104).
文摘Based on Neumman series and epsilon-algorithm, an efficient computation for dynamic responses of systems with arbitrary time-varying characteristics is investigated. Avoiding the calculation for the inverses of the equivalent stiffness matrices in each time step, the computation effort of the proposed method is reduced compared with the full analysis of Newmark method. The validity and applications of the proposed method are illustrated by a 4-DOF spring-mass system with periodical time-varying stiffness properties and a truss structure with arbitrary time-varying lumped mass. It shows that good approximate results can be obtained by the proposed method compared with the responses obtained by the full analysis of Newmark method.