Semi_weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of contin...Semi_weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi_weight functions were obtained as virtual displacement and stress fields with eigenvalue?_lambda. Integral expression of fracture parameters, K Ⅰ and K Ⅱ, were obtained from reciprocal work theorem with semi_weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi_weight function method is a simple, convenient and high precision calculation method.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Since the plasticity of soil and the irregular shape of the excavation,the efficiency and stability of the traditional local radial basis function(RBF)collocation method(LRBFCM)are inadequate for analyzing three-dimen...Since the plasticity of soil and the irregular shape of the excavation,the efficiency and stability of the traditional local radial basis function(RBF)collocation method(LRBFCM)are inadequate for analyzing three-dimensional(3D)deformation of deep excavation.In this work,the technique known as the direct method,where the local influence nodes are collocated on a straight line,is introduced to optimize the LRBFCM.The direct method can improve the accuracy of the partial derivative,reduce the size effect caused by the large length-width ratio,and weaken the influence of the shape parameters on the LRBFCM.The mapping technique is adopted to transform the physical coordinates of a quadratic-type block to normalized coordinates,in which the deformation problem can easily be solved using the direct method.The stability of the LRBFCM is further modified by considering the irregular shape of 3D excavation,which is divided into several quadratic-type blocks.The soil’s plasticity is described by the Drucker-Prager(D-P)model.The improved LRBFCM is integrated with the incremental method to analyze the plasticity.Five different examples,including strip excavations and circular excavations,are presented to validate the proposed approach’s efficiency.展开更多
[Objectives]This study was conducted to investigate the physicochemical and functional properties of soluble dietary fiber(SDF)from apple pomace.[Methods]Soluble dietary fiber(SDF)from apple pomace was extracted by di...[Objectives]This study was conducted to investigate the physicochemical and functional properties of soluble dietary fiber(SDF)from apple pomace.[Methods]Soluble dietary fiber(SDF)from apple pomace was extracted by direct water extraction(W),lactic acid bacteria fermentation(F)and steam explosion(SE)respectively,and the extraction methods and physicochemical and functional characteristics were compared and analyzed.[Results]The solubility,water holding capacities,oil holding capacities and swelling capacities of W-SDF,F-SDF and SE-SDF were(2.13,3.95 and 5.13 g/g),(9.02,13.75 and 15.88 g/g),(2.13,4.08 and 5.11 g/g),and(10.82,14.03 and 15.77 ml/g),respectively.Their emulsifying activity,emulsifying stability and least gelation concentration were(30.28,47.95 and 58.72 ml/100 ml),(37.88,45.25 and 57.13 ml/100 ml),and(12.11,11.25 and 9.87%),respectively.The adsorption capacities of W-SDF,F-SDF and SE-SDF for heavy metals(Pb,As and Cu)in the intestinal environment(pH 7)were(162.7,183.5 and 197.3μmol/g),(132.8,156.7 and 168.9μmol/g),and(57.2,63.5 and 89.2μmol/g)respectively.In the gastric environment(pH 2),they were(72.8,110.5,138.9μmol/g),(82.1,112.5,135.7μmol/g),and(38.9,42.7,55.1μmol/g)respectively.[Conclusions]The study can provide a theoretical basis for functional modification and comprehensive utilization of dietary fiber from apple pomace.展开更多
Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant i...Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.展开更多
Current research on rail vehicle system vibrations primarily relies on numerical methods,with vibration transfer functions commonly derived through data fitting.However,the physical mechanisms underlying these vibrati...Current research on rail vehicle system vibrations primarily relies on numerical methods,with vibration transfer functions commonly derived through data fitting.However,the physical mechanisms underlying these vibrations are not well understood.To clarify the vibration transfer function and its characteristics,four basic input vectors are defined,and an analytical method is proposed.The vibration transfer functions of the vehicle system are solved,and their spatial coherence is analyzed.The results show that there are two spatial scales and four coherent modes in the vehicle system.The track irregularity wavelengths are combined with two spatial scales to alter the proportions of basic input vectors and then show the characteristics of spatial coherence.Four coherent modes are involved in wheel-rail force and primary suspension force;two coherent modes are involved in bogie vertical motion;and their dominant modes vary with the input frequency.On the other hand,the coherent modes involved in the bogie pitching motion and vehicle body motion are single and fixed over the whole range of frequency.This study presents an analytical method for the rapid solution of dynamic responses in vehicle systems and systematically analyzes the coherence behavior of vibration transfer functions with respect to tracking irregularity wavelengths.展开更多
In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving th...In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.展开更多
Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential...Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.展开更多
文摘Semi_weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi_weight functions were obtained as virtual displacement and stress fields with eigenvalue?_lambda. Integral expression of fracture parameters, K Ⅰ and K Ⅱ, were obtained from reciprocal work theorem with semi_weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi_weight function method is a simple, convenient and high precision calculation method.
基金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.
基金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.
基金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.
基金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.
文摘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.
基金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.
基金supported by grants from the National Natural Science Foundation of China(Nos.12172159 and 12362019).
文摘Since the plasticity of soil and the irregular shape of the excavation,the efficiency and stability of the traditional local radial basis function(RBF)collocation method(LRBFCM)are inadequate for analyzing three-dimensional(3D)deformation of deep excavation.In this work,the technique known as the direct method,where the local influence nodes are collocated on a straight line,is introduced to optimize the LRBFCM.The direct method can improve the accuracy of the partial derivative,reduce the size effect caused by the large length-width ratio,and weaken the influence of the shape parameters on the LRBFCM.The mapping technique is adopted to transform the physical coordinates of a quadratic-type block to normalized coordinates,in which the deformation problem can easily be solved using the direct method.The stability of the LRBFCM is further modified by considering the irregular shape of 3D excavation,which is divided into several quadratic-type blocks.The soil’s plasticity is described by the Drucker-Prager(D-P)model.The improved LRBFCM is integrated with the incremental method to analyze the plasticity.Five different examples,including strip excavations and circular excavations,are presented to validate the proposed approach’s efficiency.
基金Supported by High-end Talent Project of Hebei ProvinceTangshan Talent Funding Project(B202302006).
文摘[Objectives]This study was conducted to investigate the physicochemical and functional properties of soluble dietary fiber(SDF)from apple pomace.[Methods]Soluble dietary fiber(SDF)from apple pomace was extracted by direct water extraction(W),lactic acid bacteria fermentation(F)and steam explosion(SE)respectively,and the extraction methods and physicochemical and functional characteristics were compared and analyzed.[Results]The solubility,water holding capacities,oil holding capacities and swelling capacities of W-SDF,F-SDF and SE-SDF were(2.13,3.95 and 5.13 g/g),(9.02,13.75 and 15.88 g/g),(2.13,4.08 and 5.11 g/g),and(10.82,14.03 and 15.77 ml/g),respectively.Their emulsifying activity,emulsifying stability and least gelation concentration were(30.28,47.95 and 58.72 ml/100 ml),(37.88,45.25 and 57.13 ml/100 ml),and(12.11,11.25 and 9.87%),respectively.The adsorption capacities of W-SDF,F-SDF and SE-SDF for heavy metals(Pb,As and Cu)in the intestinal environment(pH 7)were(162.7,183.5 and 197.3μmol/g),(132.8,156.7 and 168.9μmol/g),and(57.2,63.5 and 89.2μmol/g)respectively.In the gastric environment(pH 2),they were(72.8,110.5,138.9μmol/g),(82.1,112.5,135.7μmol/g),and(38.9,42.7,55.1μmol/g)respectively.[Conclusions]The study can provide a theoretical basis for functional modification and comprehensive utilization of dietary fiber from apple pomace.
基金supported by the National Natural Science Foundation of China(Grant Nos.42102346,42172301).
文摘Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.
基金Supported by Fundamental Research Funds for the Central Universities(Grant No.2024QYBS031)Fundamental Research Funds for the Central Universities(Grant No.2022JBQY007)。
文摘Current research on rail vehicle system vibrations primarily relies on numerical methods,with vibration transfer functions commonly derived through data fitting.However,the physical mechanisms underlying these vibrations are not well understood.To clarify the vibration transfer function and its characteristics,four basic input vectors are defined,and an analytical method is proposed.The vibration transfer functions of the vehicle system are solved,and their spatial coherence is analyzed.The results show that there are two spatial scales and four coherent modes in the vehicle system.The track irregularity wavelengths are combined with two spatial scales to alter the proportions of basic input vectors and then show the characteristics of spatial coherence.Four coherent modes are involved in wheel-rail force and primary suspension force;two coherent modes are involved in bogie vertical motion;and their dominant modes vary with the input frequency.On the other hand,the coherent modes involved in the bogie pitching motion and vehicle body motion are single and fixed over the whole range of frequency.This study presents an analytical method for the rapid solution of dynamic responses in vehicle systems and systematically analyzes the coherence behavior of vibration transfer functions with respect to tracking irregularity wavelengths.
基金supported by the National Natural Science Foundation of China(12201228,12171047)the Fundamental Research Funds for the Central Universities(3034011102)supported by National Key R&D Program of China(2020YFA0713701).
文摘In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.
基金supported by the National Natural Science Foundation of China (Grant Nos. 12272269, 11972257,11832014 and 11472193)the Shanghai Pilot Program for Basic Researchthe Shanghai Gaofeng Project for University Academic Program Development。
文摘Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.
