In this study,an improved integrated radial basis function with nonuniform shape parameter is introduced.The proposed shape parameter varies in each support domain and is defined byθ=1/d_(max),where d_(max)is the max...In this study,an improved integrated radial basis function with nonuniform shape parameter is introduced.The proposed shape parameter varies in each support domain and is defined byθ=1/d_(max),where d_(max)is the maximum distance of any pair of nodes in the support domain.The proposed method is verified and shows good performance.The results are stable and accurate with any number of nodes and an arbitrary nodal distribution.Notably,the support domain should be large enough to obtain accurate results.This method is then applied for transient analysis of curved shell structures made from functionally graded materials with complex geometries.Through several numerical examples,the accuracy of the proposed approach is demonstrated and discussed.Additionally,the influence of various factors on the dynamic behavior of the structures,including the power-law index,different materials,loading conditions,and geometrical parameters of the structures,was investigated.展开更多
Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the ...Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the establishment of relationships between the station’s function and the damage state,function loss,and failure probability of components and professional equipment in each layer.Also,the“4 stages-6 sequences”post-earthquake repair method is present,taking into account the functional and structural characteristics of railway stations.Additionally,a novel piecewise function for the post-earthquake functional dynamic recovery of railway stations is developed.A case study is conducted on a typical railway station to demonstrate the analysis procedure.Results indicate that under fortification,rare,and extremely rare earthquake scenarios,the interlayer drift ratio(IDR)of the railway station were 1/276,1/143,and 1/52,respectively,and corresponding peak floor acceleration(PFA)were 6.31 m/s^(2),7.82 m/s^(2),and 8.57 m/s^(2),respectively.The post-earthquake function of the railway station was 93.21%,82.33%,and 64.16%of its initial function.The repair times were 6.66 days,18.65 days,and 37.42 days.The displacement-sensitive,non-structural components were identified as the most vulnerable to damage.And the first repair stage(R_(1))which was mainly used to repair structural components and non-structural transport components,accounted for the highest proportion of total repair time.展开更多
Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of thes...Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.展开更多
This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the acc...This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.展开更多
A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the re...A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.展开更多
We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground...We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.展开更多
Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instr...Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.展开更多
Lightweight porous materials with high load-bearing,damage tolerance and energy absorption(EA)as well as intelligence of shape recovery after material deformation are beneficial and critical for many applications,e.g....Lightweight porous materials with high load-bearing,damage tolerance and energy absorption(EA)as well as intelligence of shape recovery after material deformation are beneficial and critical for many applications,e.g.aerospace,automobiles,electronics,etc.Cuttlebone produced in the cuttlefish has evolved vertical walls with the optimal corrugation gradient,enabling stress homogenization,significant load bearing,and damage tolerance to protect the organism from high external pressures in the deep sea.This work illustrated that the complex hybrid wave shape in cuttlebone walls,becoming more tortuous from bottom to top,creates a lightweight,load-bearing structure with progressive failure.By mimicking the cuttlebone,a novel bionic hybrid structure(BHS)was proposed,and as a comparison,a regular corrugated structure and a straight wall structure were designed.Three types of designed structures have been successfully manufactured by laser powder bed fusion(LPBF)with NiTi powder.The LPBF-processed BHS exhibited a total porosity of 0.042% and a good dimensional accuracy with a peak deviation of 17.4μm.Microstructural analysis indicated that the LPBF-processed BHS had a strong(001)crystallographic orientation and an average size of 9.85μm.Mechanical analysis revealed the LPBF-processed BHS could withstand over 25000 times its weight without significant deformation and had the highest specific EA value(5.32 J·g^(−1))due to the absence of stress concentration and progressive wall failure during compression.Cyclic compression testing showed that LPBF-processed BHS possessed superior viscoelastic and elasticity energy dissipation capacity.Importantly,the uniform reversible phase transition from martensite to austenite in the walls enables the structure to largely recover its pre-deformation shape when heated(over 99% recovery rate).These design strategies can serve as valuable references for the development of intelligent components that possess high mechanical efficiency and shape memory capabilities.展开更多
Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the...Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the joint effect of grain size(GS)and texture on the functional properties of NiTi SMAs and the corresponding microscopic mechanisms have not been clearly understood yet.In this work,based on the phase field method,the effect of texture on the GS-dependent functional properties of NiTi SMAs,including super-elasticity(SE),one-way shape memory effect(OWSME),and stress-assisted two-way shape memory effect(SATWSME),is investigated,and the corresponding microscopic mechanisms are revealed.Moreover,the samples with discrete geometrical gradients and/or texture gradients are designed to achieve graded functional properties.The simulation results indicate that the dependence of functional properties on texture is due to the effect of crystallographic orientation on martensite transformation and reorientation,which can lead to different inelastic strains.In the designed samples with texture gradients,the stress–strain responses of sheets with various textures are different,allowing for the coordination of overall deformation of the sample by combining such sheets,with varying inelastic deformation degrees.Thus,the overall response of the sample differs from that without texture gradient,leading to the achievement of graded functional properties.The simulation results and new findings in this work contribute to a deeper understanding of the effects of texture,GS,and their interaction on the functional properties of SMAs,and provide valuable reference for the design and development of SMA-based devices with desired functional properties.展开更多
In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,wh...In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,when the performance index is complicated,because one may encounter a two-point boundary value problem of nonlinear differential algebraic equations.To be a numerical method,it is hard to exactly preserve all the specified conditions,which might deteriorate the accuracy of numerical solution.With this in mind,we develop a novel algorithm to find the solution of the optimal control problem of nonlinear Duffing oscillator,which can exactly satisfy all the required conditions for the minimality of the performance index.A new idea of shape functions method(SFM)is introduced,from which we can transform the optimal control problems to the initial value problems for the new variables,whose initial values are given arbitrarily,and meanwhile the terminal values are determined iteratively.Numerical examples confirm the high-performance of the iterative algorithms based on the SFM,which are convergence fast,and also provide very accurate solutions.The new algorithm is robust,even large noise is imposed on the input data.展开更多
Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,...Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.展开更多
We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projec...We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node...This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.展开更多
With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying micr...With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.展开更多
Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore th...Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.展开更多
This paper aims to seek expedited fatigue analysis methods using the infrared self-heating technique.The fatigue analysis of NiTi shape memory alloys is obtained through a hybrid approach:fatigue tests to failure yiel...This paper aims to seek expedited fatigue analysis methods using the infrared self-heating technique.The fatigue analysis of NiTi shape memory alloys is obtained through a hybrid approach:fatigue tests to failure yield relatively shorter fatigue lives,while determining the fatigue limit,normally involving extremely high cycles approaching 107 cycles,is directly achieved via self-heating tests.This methodology significantly reduces testing cycles,costing only a fraction of the several-thousand-cycle tests typically required.The validity of this approach is successfully demonstrated through fatigue testing of 18Ni steel:the entire S–N curve is examined using the traditional fatigue test until a life of up to 10^(7) cycles,and the indicated fatigue limit agrees well with the one directly determined through the self-heating method.Subsequently,this developed approach is applied to the fatigue analysis of shape memory alloys under complex loading,enabling the concurrent estimation of the limits of phase transformation-dominated low-cycle fatigue and high-cycle fatigue in the elastic regime on a single specimen.The results obtained align well with other supporting evidence.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on...To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on postearthquake loss and recovery time is improved.A threelevel function tree model is established,which can consider the dynamic changes in weight coefficients of different category of components relative to their functional losses.Bayesian networks are utilized to quantify the impact of weather conditions,construction technology levels,and worker skill levels on component repair time.A method for determining the real-time functional recovery curve of buildings based on the component repair process is proposed.Taking a three-story teaching building as an example,the seismic resilience indices under basic earthquakes and rare earthquakes are calculated.The results show that the seismic resilience grade of the teaching building is comprehensively judged as GradeⅢ,and its resilience grade is more significantly affected by postearthquake loss.The proposed method can be used to predict the seismic resilience of buildings prior to earthquakes,identify weak components within buildings,and provide guidance for taking measures to enhance the seismic resilience of buildings.展开更多
The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wav...The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wave functions.Two plasma mediums,namely,the Debye plasma and quantum plasma environments are considered.The screened Coulomb potential(SCP)obtained from Debye-Hückel model is used to represent Debye plasma environments and the exponential cosine screened Coulomb potential(ECSCP)obtained from a modified Debye-Hückel model is used to represent quantum plasma environments.The resonance parameters(resonance positions and widths)are presented in terms of the screening parameters.展开更多
基金Ho Chi Minh City University of Technology (HCMUT), VNU-HCM for supporting this study
文摘In this study,an improved integrated radial basis function with nonuniform shape parameter is introduced.The proposed shape parameter varies in each support domain and is defined byθ=1/d_(max),where d_(max)is the maximum distance of any pair of nodes in the support domain.The proposed method is verified and shows good performance.The results are stable and accurate with any number of nodes and an arbitrary nodal distribution.Notably,the support domain should be large enough to obtain accurate results.This method is then applied for transient analysis of curved shell structures made from functionally graded materials with complex geometries.Through several numerical examples,the accuracy of the proposed approach is demonstrated and discussed.Additionally,the influence of various factors on the dynamic behavior of the structures,including the power-law index,different materials,loading conditions,and geometrical parameters of the structures,was investigated.
基金National Natural Science Foundation of China under Grant No.52278534the Sichuan Provincial Natural Science Foundation of China under Grant No.2022NSFSC0423。
文摘Quantifying the post-earthquake functional recovery of railway stations presents significant challenges.This paper first establishes a post-earthquake function calculation method for railway stations,encompassing the establishment of relationships between the station’s function and the damage state,function loss,and failure probability of components and professional equipment in each layer.Also,the“4 stages-6 sequences”post-earthquake repair method is present,taking into account the functional and structural characteristics of railway stations.Additionally,a novel piecewise function for the post-earthquake functional dynamic recovery of railway stations is developed.A case study is conducted on a typical railway station to demonstrate the analysis procedure.Results indicate that under fortification,rare,and extremely rare earthquake scenarios,the interlayer drift ratio(IDR)of the railway station were 1/276,1/143,and 1/52,respectively,and corresponding peak floor acceleration(PFA)were 6.31 m/s^(2),7.82 m/s^(2),and 8.57 m/s^(2),respectively.The post-earthquake function of the railway station was 93.21%,82.33%,and 64.16%of its initial function.The repair times were 6.66 days,18.65 days,and 37.42 days.The displacement-sensitive,non-structural components were identified as the most vulnerable to damage.And the first repair stage(R_(1))which was mainly used to repair structural components and non-structural transport components,accounted for the highest proportion of total repair time.
基金supported by the Open Funds for Hubei Key Laboratory of Marine Geological Resources,China University of Geosciences(No.MGR202308)the Natural Science Foundation of Shandong Province(No.ZR2020MD085)+3 种基金the National Natural Science Foundation of China(No.41821004)the Taishan Scholar Program(No.tstp2022114)the Shandong Provincial Natural Science Foundation(No.DKXZZ202206)the National Key Research and Development Program of China(No.2016YFC1402404).
文摘Ocean remote sensing satellites provide observations with high spatiotemporal resolution.However,the influence of clouds,fog,and haze frequently leads to significant data gaps.Accurate and effective estimation of these missing data is highly valuable for engineering and scientific research.In this study,the radial basis function(RBF)method is used to estimate the spatial distribution of total suspended matter(TSM)concentration in Hangzhou Bay using remote sensing data with severe data gaps.The estimation precision is validated by comparing the results with those of other commonly used interpolation methods,such as the Kriging method and the basic spline(B-spline)method.In addition,the applicability of the RBF method is explored.Results show that the estimation of the RBF method is significantly close to the observation in Hangzhou Bay.The average of the mean absolute error,mean relative error,and root mean square error in all the experiments is evidently smaller than those of the Kriging and B-spline interpolations,indicating that the proposed method is more appropriate for estimating the spatial distribution of the TSM in Hangzhou Bay.Finally,the TSM distribution in the blank observational area is predicted.This study can provide some reference values for handling watercolor remote sensing data.
基金supported by the Advance Research Project of Civil Aerospace Technology(Grant No.D020304)National Nat-ural Science Foundation of China(Grant Nos.52205257 and U22B2083).
文摘This paper proposes a new step-by-step Chebyshev space-time spectral method to analyze the force vibration of functionally graded material structures.Although traditional space-time spectral methods can reduce the accuracy mismatch between tem-poral low-order finite difference and spatial high-order discre tization,the ir time collocation points must increase dramatically to solve highly oscillatory solutions of structural vibration,which results in a surge in computing time and a decrease in accuracy.To address this problem,we introduced the step-by-step idea in the space-time spectral method.The Chebyshev polynomials and Lagrange's equation were applied to derive discrete spatial goverming equations,and a matrix projection method was used to map the calculation results of prev ious steps as the initial conditions of the subsequent steps.A series of numerical experiments were carried out.The results of the proposed method were compared with those obtained by traditional space-time spectral methods,which showed that higher accuracy could be achieved in a shorter computation time than the latter in highly oscillatory cases.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.12174048 and 12204128)。
文摘A Wentzel-Kramers-Brillouin(WKB)method is introduced for obtaining a uniform asymptotic solution for underwater sound propagation at very low frequencies in deep ocean.The method utilizes a mode sum and employs the reference functions method to describe the solution to the depth-separated wave equation approximately using parabolic cylinder functions.The conditions for the validity of this approximation are also discussed.Furthermore,a formula that incorporates waveguide effects for the modal group velocity is derived,revealing that boundary effects at very low frequencies can have a significant impact on the propagation characteristics of even low-order normal modes.The present method not only offers improved accuracy compared to the classical WKB approximation and the uniform asymptotic approximation based on Airy functions,but also provides a wider range of depth applicability.Additionally,this method exhibits strong agreement with numerical methods and offers valuable physical insights.Finally,the method is applied to the study of very-low-frequency sound propagation in the South China Sea,leading to sound transmission loss predictions that closely align with experimental observations.
基金supported by the National Key R&D Program of China(No.2023YFA1606701)the National Natural Science Foundation of China(Nos.12175042,11890710,11890714,12047514,12147101,and 12347106)+1 种基金Guangdong Major Project of Basic and Applied Basic Research(No.2020B0301030008)China National Key R&D Program(No.2022YFA1602402).
文摘We employed random distributions and gradient descent methods for the Generator Coordinate Method(GCM)to identify effective basis wave functions,taking halo nuclei ^(6)He and ^(6)Li as examples.By comparing the ground state(0^(+))energy of ^(6)He and the excited state(0^(+))energy of 6 Li calculated with various random distributions and manually selected generation coordinates,we found that the heavy tail characteristic of the logistic distribution better describes the features of the halo nuclei.Subsequently,the Adam algorithm from machine learning was applied to optimize the basis wave functions,indicating that a limited number of basis wave functions can approximate the converged values.These results offer some empirical insights for selecting basis wave functions and contribute to the broader application of machine learning methods in predicting effective basis wave functions.
基金supported by the Center for Advanced Systems Understanding(CASUS),financed by Germany’s Federal Ministry of Education and Research(BMBF)and the Saxon State Government out of the State Budget approved by the Saxon State Parliamentfunding from the European Research Council(ERC)under the European Union’s Horizon 2022 research and innovation programme(Grant Agreement No.101076233,“PREXTREME”)funding from the European Union’s Just Transition Fund(JTF)within the project Röntgenlaser-Optimierung der Laserfusion(ROLF),Contract No.5086999001,co-financed by the Saxon State Government out of the State Budget approved by the Saxon State Parliament.
文摘Ab initio modeling of dynamic structure factors(DSF)and related density response properties in the warm dense matter(WDM)regime is a challenging computational task.The DSF,convolved with a probing X-ray beam and instrument function,is measured in X-ray Thom-son scattering(XRTS)experiments,which allow the study of electronic structure properties at the microscopic level.Among the various ab initio methods,linear-response time-dependent density-functional theory(LR-TDDFT)is a key framework for simulating the DSF.The standard approach in LR-TDDFT for computing the DSF relies on the orbital representation.A significant drawback of this method is the unfavorable scaling of the number of required empty bands as the wavenumber increases,making LR-TDDFT impractical for modeling XRTS measurements over large energy scales,such as in backward scattering geometry.In this work,we consider and test an alternative approach to LR-TDDFT that employs the Liouville–Lanczos(LL)method for simulating the DSF of WDM.This approach does not require empty states and allows the DSF at large momentum transfer values and over a broad frequency range to be accessed.We compare the results obtained from the LL method with those from the solution of Dyson’s equation using the standard LR-TDDFT within the projector augmented-wave formalism for isochorically heated aluminum and warm dense hydrogen.Additionally,we utilize exact path integral Monte Carlo results for the imaginary-time density-density correlation function(ITCF)of warm dense hydrogen to rigorously benchmark the LL approach.We discuss the application of the LL method for calculating DSFs and ITCFs at different wavenumbers,the effects of pseudopotentials,and the role of Lorentzian smearing.The successful validation of the LL method under WDM conditions makes it a valuable addition to the ab initio simulation landscape,supporting experimental efforts and advancing WDM theory.
基金supported by the National Natural Science Foundation of China(Grant No.52225503)National Key Research and Development Program of China(Grant No.2022YFB3805701)+1 种基金Development Program of Jiangsu Province(Grant Nos.BE2022069 and BE2022069-1)Postgraduate Research&Practice Innovation Program of Jiangsu Province(Grant No.KYCX21-0207).
文摘Lightweight porous materials with high load-bearing,damage tolerance and energy absorption(EA)as well as intelligence of shape recovery after material deformation are beneficial and critical for many applications,e.g.aerospace,automobiles,electronics,etc.Cuttlebone produced in the cuttlefish has evolved vertical walls with the optimal corrugation gradient,enabling stress homogenization,significant load bearing,and damage tolerance to protect the organism from high external pressures in the deep sea.This work illustrated that the complex hybrid wave shape in cuttlebone walls,becoming more tortuous from bottom to top,creates a lightweight,load-bearing structure with progressive failure.By mimicking the cuttlebone,a novel bionic hybrid structure(BHS)was proposed,and as a comparison,a regular corrugated structure and a straight wall structure were designed.Three types of designed structures have been successfully manufactured by laser powder bed fusion(LPBF)with NiTi powder.The LPBF-processed BHS exhibited a total porosity of 0.042% and a good dimensional accuracy with a peak deviation of 17.4μm.Microstructural analysis indicated that the LPBF-processed BHS had a strong(001)crystallographic orientation and an average size of 9.85μm.Mechanical analysis revealed the LPBF-processed BHS could withstand over 25000 times its weight without significant deformation and had the highest specific EA value(5.32 J·g^(−1))due to the absence of stress concentration and progressive wall failure during compression.Cyclic compression testing showed that LPBF-processed BHS possessed superior viscoelastic and elasticity energy dissipation capacity.Importantly,the uniform reversible phase transition from martensite to austenite in the walls enables the structure to largely recover its pre-deformation shape when heated(over 99% recovery rate).These design strategies can serve as valuable references for the development of intelligent components that possess high mechanical efficiency and shape memory capabilities.
基金The National Natural Science Foundation of China(12202294 and 12022208)the Project funded by China Postdoctoral Science Foundation(2022M712243)the Fundamental Research Funds for the Central Universities(2023SCU12098)are acknowledged.
文摘Texture is inevitably introduced during the manufacturing of most NiTi shape memory alloys(SMAs),and the textured nanocrystalline NiTi has been extensively employed in engineering.However,the effect of texture,and the joint effect of grain size(GS)and texture on the functional properties of NiTi SMAs and the corresponding microscopic mechanisms have not been clearly understood yet.In this work,based on the phase field method,the effect of texture on the GS-dependent functional properties of NiTi SMAs,including super-elasticity(SE),one-way shape memory effect(OWSME),and stress-assisted two-way shape memory effect(SATWSME),is investigated,and the corresponding microscopic mechanisms are revealed.Moreover,the samples with discrete geometrical gradients and/or texture gradients are designed to achieve graded functional properties.The simulation results indicate that the dependence of functional properties on texture is due to the effect of crystallographic orientation on martensite transformation and reorientation,which can lead to different inelastic strains.In the designed samples with texture gradients,the stress–strain responses of sheets with various textures are different,allowing for the coordination of overall deformation of the sample by combining such sheets,with varying inelastic deformation degrees.Thus,the overall response of the sample differs from that without texture gradient,leading to the achievement of graded functional properties.The simulation results and new findings in this work contribute to a deeper understanding of the effects of texture,GS,and their interaction on the functional properties of SMAs,and provide valuable reference for the design and development of SMA-based devices with desired functional properties.
文摘In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,when the performance index is complicated,because one may encounter a two-point boundary value problem of nonlinear differential algebraic equations.To be a numerical method,it is hard to exactly preserve all the specified conditions,which might deteriorate the accuracy of numerical solution.With this in mind,we develop a novel algorithm to find the solution of the optimal control problem of nonlinear Duffing oscillator,which can exactly satisfy all the required conditions for the minimality of the performance index.A new idea of shape functions method(SFM)is introduced,from which we can transform the optimal control problems to the initial value problems for the new variables,whose initial values are given arbitrarily,and meanwhile the terminal values are determined iteratively.Numerical examples confirm the high-performance of the iterative algorithms based on the SFM,which are convergence fast,and also provide very accurate solutions.The new algorithm is robust,even large noise is imposed on the input data.
基金supported by the Office of Naval Research (Grant No.N000141010778) to the University of Oklahomathe National Natural Sciences Foundation of China (Grant Nos. 40930950,41075043,and 4092116037) to the Institute of Atmospheric Physicsprovided by NOAA/Office of Oceanic and Atmospheric Research under NOAA-University of Oklahoma Cooperative Agreement No. (NA17RJ1227),U.S. Department of Commerce
文摘Built on the integral formulas in Part I,numerical methods are developed for computing velocity potential and streamfunction in a limited domain.When there is no inner boundary（around a data hole） inside the domain,the total solution is the sum of the internally and externally induced parts.For the internally induced part,three numerical schemes（grid-staggering,local-nesting and piecewise continuous integration） are designed to deal with the singularity of the Green＇s function encountered in numerical calculations.For the externally induced part,by setting the velocity potential（or streamfunction） component to zero,the other component of the solution can be computed in two ways：（1） Solve for the density function from its boundary integral equation and then construct the solution from the boundary integral of the density function.（2） Use the Cauchy integral to construct the solution directly.The boundary integral can be discretized on a uniform grid along the boundary.By using local-nesting（or piecewise continuous integration）,the scheme is refined to enhance the discretization accuracy of the boundary integral around each corner point（or along the entire boundary）.When the domain is not free of data holes,the total solution contains a data-hole-induced part,and the Cauchy integral method is extended to construct the externally induced solution with irregular external and internal boundaries.An automated algorithm is designed to facilitate the integrations along the irregular external and internal boundaries.Numerical experiments are performed to evaluate the accuracy and efficiency of each scheme relative to others.
基金financial support from the National Natural Science Foundation of China (Grant No. 12227901)the financial support from the National Natural Science Foundation of China (Grant Nos. 11974263 and 12174291)。
文摘We report a linear-scaling random Green's function(rGF) method for large-scale electronic structure calculation. In this method, the rGF is defined on a set of random states and is efficiently calculated by projecting onto Krylov subspace. With the rGF method, the Fermi–Dirac operator can be obtained directly, avoiding the polynomial expansion to Fermi–Dirac function. To demonstrate the applicability, we implement the rGF method with the density-functional tight-binding method. It is shown that the Krylov subspace can maintain at small size for materials with different gaps at zero temperature, including H_(2)O and Si clusters. We find with a simple deflation technique that the rGF self-consistent calculation of H_(2)O clusters at T = 0 K can reach an error of~ 1 me V per H_(2)O molecule in total energy, compared to deterministic calculations. The rGF method provides an effective stochastic method for large-scale electronic structure simulation.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金Anhui Provincial Natural Science Foundation(2308085QD124)Anhui Province University Natural Science Research Project(GrantNo.2023AH050918)The University Outstanding Youth Talent Support Program of Anhui Province.
文摘This study proposes an effective method to enhance the accuracy of the Differential Quadrature Method(DQM)for calculating the dynamic characteristics of functionally graded beams by improving the form of discrete node distribution.Firstly,based on the first-order shear deformation theory,the governing equation of free vibration of a functionally graded beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam axial displacement,transverse displacement,and cross-sectional rotation angle by considering the effects of shear deformation and rotational inertia of the beam cross-section.Then,ignoring the shear deformation of the beam section and only considering the effect of the rotational inertia of the section,the governing equation of the beam is transformed into the eigenvalue problem of ordinary differential equations with respect to beam transverse displacement.Based on the differential quadrature method theory,the eigenvalue problem of ordinary differential equations is transformed into the eigenvalue problem of standard generalized algebraic equations.Finally,the first several natural frequencies of the beam can be calculated.The feasibility and accuracy of the improved DQM are verified using the finite element method(FEM)and combined with the results of relevant literature.
基金the National Key Research and Development Program of China(Grant Number 2021YFB1714600)the National Natural Science Foundation of China(Grant Number 52075195)the Fundamental Research Funds for the Central Universities,China through Program No.2172019kfyXJJS078.
文摘With the continuous advancement in topology optimization and additive manufacturing(AM)technology,the capability to fabricate functionally graded materials and intricate cellular structures with spatially varying microstructures has grown significantly.However,a critical challenge is encountered in the design of these structures–the absence of robust interface connections between adjacent microstructures,potentially resulting in diminished efficiency or macroscopic failure.A Hybrid Level Set Method(HLSM)is proposed,specifically designed to enhance connectivity among non-uniform microstructures,contributing to the design of functionally graded cellular structures.The HLSM introduces a pioneering algorithm for effectively blending heterogeneous microstructure interfaces.Initially,an interpolation algorithm is presented to construct transition microstructures seamlessly connected on both sides.Subsequently,the algorithm enables the morphing of non-uniform unit cells to seamlessly adapt to interconnected adjacent microstructures.The method,seamlessly integrated into a multi-scale topology optimization framework using the level set method,exhibits its efficacy through numerical examples,showcasing its prowess in optimizing 2D and 3D functionally graded materials(FGM)and multi-scale topology optimization.In essence,the pressing issue of interface connections in complex structure design is not only addressed but also a robust methodology is introduced,substantiated by numerical evidence,advancing optimization capabilities in the realm of functionally graded materials and cellular structures.
基金the National Natural Science Foundation of China(Nos.12302007,12372006,and 12202109)the Specific Research Project of Guangxi for Research Bases and Talents(No.AD23026051)。
文摘Functionally graded materials(FGMs)are a novel class of composite materials that have attracted significant attention in the field of engineering due to their unique mechanical properties.This study aims to explore the dynamic behaviors of an FGM stepped beam with different boundary conditions based on an efficient solving method.Under the assumptions of the Euler-Bernoulli beam theory,the governing differential equations of an individual FGM beam are derived with Hamilton’s principle and decoupled via the separation-of-variable approach.Then,the free and forced vibrations of the FGM stepped beam are solved with the transfer matrix method(TMM).Two models,i.e.,a three-level FGM stepped beam and a five-level FGM stepped beam,are considered,and their natural frequencies and mode shapes are presented.To demonstrate the validity of the method in this paper,the simulation results by ABAQUS are also given.On this basis,the detailed parametric analyses on the frequencies and dynamic responses of the three-level FGM stepped beam are carried out.The results show the accuracy and efficiency of the TMM.
基金National Natural Science Foundation of China,12272305,Yahui Zhang,12372123,Xiaojun GuBasic Research Program,JCKY2022603C016,Xiaojun Gu。
文摘This paper aims to seek expedited fatigue analysis methods using the infrared self-heating technique.The fatigue analysis of NiTi shape memory alloys is obtained through a hybrid approach:fatigue tests to failure yield relatively shorter fatigue lives,while determining the fatigue limit,normally involving extremely high cycles approaching 107 cycles,is directly achieved via self-heating tests.This methodology significantly reduces testing cycles,costing only a fraction of the several-thousand-cycle tests typically required.The validity of this approach is successfully demonstrated through fatigue testing of 18Ni steel:the entire S–N curve is examined using the traditional fatigue test until a life of up to 10^(7) cycles,and the indicated fatigue limit agrees well with the one directly determined through the self-heating method.Subsequently,this developed approach is applied to the fatigue analysis of shape memory alloys under complex loading,enabling the concurrent estimation of the limits of phase transformation-dominated low-cycle fatigue and high-cycle fatigue in the elastic regime on a single specimen.The results obtained align well with other supporting evidence.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
基金The National Key Research and Development Program of China(No.2023YFC3805003)。
文摘To quantify the seismic resilience of buildings,a method for evaluating functional loss from the component level to the overall building is proposed,and the dual-parameter seismic resilience assessment method based on postearthquake loss and recovery time is improved.A threelevel function tree model is established,which can consider the dynamic changes in weight coefficients of different category of components relative to their functional losses.Bayesian networks are utilized to quantify the impact of weather conditions,construction technology levels,and worker skill levels on component repair time.A method for determining the real-time functional recovery curve of buildings based on the component repair process is proposed.Taking a three-story teaching building as an example,the seismic resilience indices under basic earthquakes and rare earthquakes are calculated.The results show that the seismic resilience grade of the teaching building is comprehensively judged as GradeⅢ,and its resilience grade is more significantly affected by postearthquake loss.The proposed method can be used to predict the seismic resilience of buildings prior to earthquakes,identify weak components within buildings,and provide guidance for taking measures to enhance the seismic resilience of buildings.
基金Supported by the Natural Science Foundation of Heilongjiang Province(LH2024A025)。
文摘The effects of plasma screening on the ^(1)P^(o) resonance states of H-and He below the n=3 and n=4 thresholds of the respective subsystemsare investigated using the stabilization method and correlated exponential wave functions.Two plasma mediums,namely,the Debye plasma and quantum plasma environments are considered.The screened Coulomb potential(SCP)obtained from Debye-Hückel model is used to represent Debye plasma environments and the exponential cosine screened Coulomb potential(ECSCP)obtained from a modified Debye-Hückel model is used to represent quantum plasma environments.The resonance parameters(resonance positions and widths)are presented in terms of the screening parameters.
