In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
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.展开更多
For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study prop...For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.展开更多
The safety of risers in hang-off states is a vital challenge in offshore oil and gas engineering.A new hang-off system installed on top of risers is proposed for improving the security of risers.This approach leads to...The safety of risers in hang-off states is a vital challenge in offshore oil and gas engineering.A new hang-off system installed on top of risers is proposed for improving the security of risers.This approach leads to a challenging problem:coupling the dynamics of risers with a new hang-off system combined with multiple structures and complex constraints.To accurately analyze the dynamic responses of the coupled system,a coupled dynamic model is established based on the Euler-Bernoulli beam-column theory and penalty function method.A comprehensive analysis method is proposed for coupled dynamic analysis by combining the finite element method and the Newmarkβmethod.An analysis program is also developed in MATLAB for dynamic simulation.The simulation results show that the dynamic performances of the risers at the top part are significantly improved by the new hang-off system,especially the novel design,which includes the centralizer and articulation joint.The bending moment and lateral deformation of the risers at the top part decrease,while the hang-off joint experiences a great bending moment at the bottom of the lateral restraint area which requires particular attention in design and application.The platform navigation speed range under the safety limits of risers expands with the new hang-off system in use.展开更多
The structural transformation from a liquid into a crystalline solid is an important subject in condensed matter physics and materials science. In the present study, first-principles molecular dynamics calculations ar...The structural transformation from a liquid into a crystalline solid is an important subject in condensed matter physics and materials science. In the present study, first-principles molecular dynamics calculations are performed to investigate the structure and properties of aluminum during the solidification which is induced by cooling and compression. In the cooling process and compression process, it is found that the icosahedral short-range order is initially enhanced and then begin to decay, the face-centered cubic short-range order eventually becomes dominant before it transforms into a crystalline solid.展开更多
High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic ...High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic mechanical properties of shale.Dynamic experiments on shale subjected to true triaxial compression at different strain rates are first conducted in this research.The dynamic stress-strain curves,peak strain,peak stress and failure modes of shale are investigated.The results of the study indicate that the intermediate principal stress and the minor principal stress have the significant influence on the dynamic mechanical behaviors,although this effect decreases as the strain rate increases.The characteristics of compression-shear failure primarily occur in shale subjected to triaxial compression at high strain rates,which distinguishes it from the fragmentation characteristics observed in shale under dynamic uniaxial compression.Additionally,a numerical three-dimensional Split Hopkinson Pressure Bar(3D-SHPB),which is established by coupling PFC3D and FLAC3D methods,is validated to replicate the laboratory characteristics of shale.The dynamic mechanical characteristics of shale subjected to different confining stresses are systematically investigated by the coupling PFC3D and FLAC3D method.The numerical results are in good agreement with the experimental data.展开更多
This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology co...This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.展开更多
Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite ano...Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.展开更多
The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scann...The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.展开更多
The sloshing in a group of rigid cylindrical tanks with baffles and on soil foundation under horizontal excitation is studied analytically.The solutions for the velocity potential are derived out by the liquid subdoma...The sloshing in a group of rigid cylindrical tanks with baffles and on soil foundation under horizontal excitation is studied analytically.The solutions for the velocity potential are derived out by the liquid subdomain method.Equivalent models with mass-spring oscillators are established to replace continuous fluid.Combined with the least square technique,Chebyshev polynomials are employed to fit horizontal,rocking and horizontal-rocking coupling impedances of soil,respectively.A lumped parameter model for impedance is presented to describe the effects of soil on tank structures.A mechanical model for the soil-foundation-tank-liquid-baffle system with small amount of calculation and high accuracy is proposed using the substructure technique.The analytical solutions are in comparison with data from reported literature and numerical codes to validate the effectiveness and correctness of the model.Detailed dynamic properties and seismic responses of the soil-tank system are given for the baffle number,size and location as well as soil parameter.展开更多
A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid ...A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.展开更多
Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the productio...Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.展开更多
The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-...The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-tion(DMD).The pseudo-time marching of every real-time step is approximated as an infinite-dimensional linear dynamical system.Thereafter,DMD is utilized to analyze the adjoint vectors sampled from these pseudo-time marching.First-order zero frequency mode is selected to accelerate the pseudo-time marching of unsteady adjoint equations in every real-time step.Through flow past a stationary circular cylinder and an unsteady aerodynamic shape optimization example,the efficiency of solving unsteady adjoint equations is significantly improved.Re-sults show that one hundred adjoint vectors contains enough information about the pseudo-time dynamics,and the adjoint dominant mode can be precisely predicted only by five snapshots produced from the adjoint vectors,which indicates DMD analysis for pseudo-time marching of unsteady adjoint equations is efficient.展开更多
We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatia...We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.展开更多
In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avo...In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.展开更多
The underground hydropower projects in Southwest China is characterized by large excavation sizes,high geostresses,complicated geological conditions and multiple construction processes.Various disasters such as collap...The underground hydropower projects in Southwest China is characterized by large excavation sizes,high geostresses,complicated geological conditions and multiple construction processes.Various disasters such as collapses,large deformations,rockbursts are frequently encountered,resulting in serious casualties and huge economic losses.This review mainly presents some representative results on microseismic(MS)monitoring and forecasting for disasters in hydropower underground engineering.First,a set of new denoising,spectral analysis,and location methods were developed for better identification and location of MS signals.Then,the tempo-spatial characteristics of MS events were analyzed to understand the relationship between field construction and damages of surrounding rocks.Combined with field construction,geological data,numerical simulation and parametric analysis of MS sources,the focal mechanism of MS events was revealed.A damage constitutive model considering MS fracturing size was put forward and feedback analysis considering the MS damage of underground surrounding rocks was conducted.Next,an MS multi-parameter based risk assessment and early warning method for dynamic disasters were proposed.The technology for control of the damage and deformation of underground surrounding rocks was proposed for underground caverns.Finally,two typical underground powerhouses were selected as case studies.These achievements can provide significant references for prevention and control of dynamic disasters for underground engineering with similar complicated geological conditions.展开更多
Dynamic coupling modeling and analysis of rotating beams based on the nonlinear Green-Lagrangian strain are introduced in this work.With the reservation of the axial nonlinear strain,there are more coupling terms for ...Dynamic coupling modeling and analysis of rotating beams based on the nonlinear Green-Lagrangian strain are introduced in this work.With the reservation of the axial nonlinear strain,there are more coupling terms for axial and transverse deformations.The discretized dynamic governing equations are obtained by using the finite element method and Lagrange’s equations of the second kind.Time responses are conducted to compare the proposed model with other previous models.The stretching deformation due to rotating motion is observed and calculated by special formulations under dynamic equilibrium.The stretching deformation and the change of the associated equilibrium position are taken into account to analyze the free vibration and frequency response of the rotating beams.Analytical and numerical comparisons show that the proposed model can provide reliable results,while the previous models may lead to imprecise results,especially in high-speed conditions.展开更多
The solubility of 2,3,4-trichloro-1,5-dinitrobenzene(TCDNB) was measured by a laser dynamic method over the temperature range from 278.15 K to 323.15 K under 0.1 MPa in fifteen mono-solvents(methanol,ethanol,isopropan...The solubility of 2,3,4-trichloro-1,5-dinitrobenzene(TCDNB) was measured by a laser dynamic method over the temperature range from 278.15 K to 323.15 K under 0.1 MPa in fifteen mono-solvents(methanol,ethanol,isopropanol,n-butanol,toluene,dichloromethane,chloroform,tetrachloromethane,1,2-dichloroethane,acetone,ethyl acetate,acetonitrile,N-methylpyrrolidone(NMP),N,Ndimethylformamide dimethyl sulfoxide(DMF),dimethyl sulfoxide(DMSO).The solubility of TCDNB could be increased with increasing temperature in fifteen mono-solvents.TCDNB solubility is in the following order at 298.15 K:NMP>DMF>DMSO>toluene>acetone>ethyl acetate>dichloromethane>1,2-dich loroethane>chloroform>acetonitrile>tetrachloromethane>methanol>ethanol>n-butanol>isopropanol.The KAT-LSER model was used to investigate the solvent effect,which revealed that the hydrogen bond acidity of solvents has a greater effect on TCDNB solubility.The van't Hoff model,the modified Apelblat model,theλh model,and the non-random two liquid(NRTL)model were used to correlate the solubility of TCDNB.The calculated solubility data agreed well with the experimental data,and the modified Apelblat model fit best.Furthermore,the van't Hoff and Gibbs equations were also used to calculate the dissolution thermodynamic properties of TCDNB in various solvents.TCDNB dissolution could be an enthalpy-driven,non-spontaneous,and endothermic process in fifteen mono-solvents.The determination and fitting solubility of TCDNB,as well as the calculation of its thermodynamic properties,would be critical in the purification and crystallization of its preparation process research.展开更多
We present a study on the dynamic stability of porous functionally graded(PFG)beams under hygro-thermal loading.The variations of the properties of the beams across the beam thicknesses are described by the power-law ...We present a study on the dynamic stability of porous functionally graded(PFG)beams under hygro-thermal loading.The variations of the properties of the beams across the beam thicknesses are described by the power-law model.Unlike most studies on this topic,we consider both the bending deformation of the beams and the hygro-thermal load as size-dependent,simultaneously,by adopting the equivalent differential forms of the well-posed nonlocal strain gradient integral theory(NSGIT)which are strictly equipped with a set of constitutive boundary conditions(CBCs),and through which both the stiffness-hardening and stiffness-softening effects of the structures can be observed with the length-scale parameters changed.All the variables presented in the differential problem formulation are discretized.The numerical solution of the dynamic instability region(DIR)of various bounded beams is then developed via the generalized differential quadrature method(GDQM).After verifying the present formulation and results,we examine the effects of different parameters such as the nonlocal/gradient length-scale parameters,the static force factor,the functionally graded(FG)parameter,and the porosity parameter on the DIR.Furthermore,the influence of considering the size-dependent hygro-thermal load is also presented.展开更多
Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System (IAP-DCP) developed by the Institute of Atmospheric Physics, a correction met...Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System (IAP-DCP) developed by the Institute of Atmospheric Physics, a correction method that can account for the dependence of model's systematic biases on SST anomalies is proposed. It is shown that this correction method can improve the hindcast skill of the IAP-DCP for summer rainfall anomalies over China, especially in western China and southeast China, which may imply its potential application to real-time seasonal prediction.展开更多
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金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.
基金National Natural Science Foundation of China under Grant Nos.51978213 and 51778190the National Key Research and Development Program of China under Grant Nos.2017YFC0703605 and 2016YFC0701106。
文摘For real-time dynamic substructure testing(RTDST),the influence of the inertia force of fluid specimens on the stability and accuracy of the integration algorithms has never been investigated.Therefore,this study proposes to investigate the stability and accuracy of the central difference method(CDM)for RTDST considering the specimen mass participation coefficient.First,the theory of the CDM for RTDST is presented.Next,the stability and accuracy of the CDM for RTDST considering the specimen mass participation coefficient are investigated.Finally,numerical simulations and experimental tests are conducted for verifying the effectiveness of the method.The study indicates that the stability of the algorithm is affected by the mass participation coefficient of the specimen,and the stability limit first increases and then decreases as the mass participation coefficient increases.In most cases,the mass participation coefficient will increase the stability limit of the algorithm,but in specific circumstances,the algorithm may lose its stability.The stability and accuracy of the CDM considering the mass participation coefficient are verified by numerical simulations and experimental tests on a three-story frame structure with a tuned liquid damper.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52271300,52071337,and 51809279)the National Key Research and Development Program of China(Grant No.2022YFC2806501)the High-tech Ship Research Projects Sponsored by MIIT(Grant No.CBG2N21-4-2-5).
文摘The safety of risers in hang-off states is a vital challenge in offshore oil and gas engineering.A new hang-off system installed on top of risers is proposed for improving the security of risers.This approach leads to a challenging problem:coupling the dynamics of risers with a new hang-off system combined with multiple structures and complex constraints.To accurately analyze the dynamic responses of the coupled system,a coupled dynamic model is established based on the Euler-Bernoulli beam-column theory and penalty function method.A comprehensive analysis method is proposed for coupled dynamic analysis by combining the finite element method and the Newmarkβmethod.An analysis program is also developed in MATLAB for dynamic simulation.The simulation results show that the dynamic performances of the risers at the top part are significantly improved by the new hang-off system,especially the novel design,which includes the centralizer and articulation joint.The bending moment and lateral deformation of the risers at the top part decrease,while the hang-off joint experiences a great bending moment at the bottom of the lateral restraint area which requires particular attention in design and application.The platform navigation speed range under the safety limits of risers expands with the new hang-off system in use.
基金Project supported by the National Natural Science Foundation of China(Grant No.51701180)the Foundation of the State Key Laboratory of Coal Conversion,China(Grant No.J22-23-103)。
文摘The structural transformation from a liquid into a crystalline solid is an important subject in condensed matter physics and materials science. In the present study, first-principles molecular dynamics calculations are performed to investigate the structure and properties of aluminum during the solidification which is induced by cooling and compression. In the cooling process and compression process, it is found that the icosahedral short-range order is initially enhanced and then begin to decay, the face-centered cubic short-range order eventually becomes dominant before it transforms into a crystalline solid.
基金supported by the National Natural Science Foundation of China(Nos.51839009 and 52027814)the Natural Science Foundation of Hubei Province(No.2023AFB589).
文摘High-energy gas fracturing of shale is a novel,high efficacy and eco-friendly mining technique,which is a typical dynamic perturbing behavior.To effectively extract shale gas,it is important to understand the dynamic mechanical properties of shale.Dynamic experiments on shale subjected to true triaxial compression at different strain rates are first conducted in this research.The dynamic stress-strain curves,peak strain,peak stress and failure modes of shale are investigated.The results of the study indicate that the intermediate principal stress and the minor principal stress have the significant influence on the dynamic mechanical behaviors,although this effect decreases as the strain rate increases.The characteristics of compression-shear failure primarily occur in shale subjected to triaxial compression at high strain rates,which distinguishes it from the fragmentation characteristics observed in shale under dynamic uniaxial compression.Additionally,a numerical three-dimensional Split Hopkinson Pressure Bar(3D-SHPB),which is established by coupling PFC3D and FLAC3D methods,is validated to replicate the laboratory characteristics of shale.The dynamic mechanical characteristics of shale subjected to different confining stresses are systematically investigated by the coupling PFC3D and FLAC3D method.The numerical results are in good agreement with the experimental data.
基金supported by the National Natural Science Foundation of China (Grant No. 41877524, No. 42172320, No. 41971214)。
文摘This study employs the smoothed particle hydrodynamics–finite element method(SPH–FEM) coupling numerical method to investigate the impact of debris flow on reinforced concrete(RC)-frame buildings. The methodology considers the variables of debris flow depth and velocity and introduces the intensity index IDV(IDV = DV) to evaluate three different levels of debris flow impact intensity. The primary focus of this study is to investigate the dynamic response and failure mechanism of RC-frame buildings under debris flow impact, including structural failure patterns, impact force and column displacement. The results show that under a highintensity impact, a gradual collapse process of the RCframe building can be observed, and the damage mode of the frame column reflects shear failure or plastic hinge failure mechanism. First, the longitudinal infill walls are damaged owing to their low out-of-plane flexural capacity;the critical failure intensity index IDV value is approximately 7.5 m2/s. The structure cannot withstand debris flows with an intensity index IDV greater than 16 m2/s, and it is recommended that the peak impact force should not exceed 2100 k N. The impact damage ability of debris flow on buildings mostly originates from the impact force of the frontal debris flow, with the impact force of the debris flow body being approximately 42% lower than that of the debris flow head. Finally, a five-level classification system for evaluating the damage status of buildings is proposed based on the numerical simulation and investigation results of the disaster site.
基金supported by the National Natural Science Foundation of China(T2322015,22209094,22209093,and 22109086)the National Key Research and Development Program(2021YFB2500300)+2 种基金the Open Research Fund of CNMGE Platform&NSCC-TJOrdos-Tsinghua Innovative&Collaborative Research Program in Carbon Neutralitythe Tsinghua University Initiative Scientific Research Program。
文摘Electric double layer(EDL)is a critical topic in electrochemistry and largely determines the working performance of lithium batteries.However,atomic insights into the EDL structures on heteroatom-modified graphite anodes and EDL evolution with electrode potential are very lacking.Herein,a constant-potential molecular dynamics(CPMD)method is proposed to probe the EDL structure under working conditions,taking N-doped graphite electrodes and carbonate electrolytes as an example.An interface model was developed,incorporating the electrode potential and atom electronegativities.As a result,an insightful atomic scenario for the EDL structure under varied electrode potentials has been established,which unveils the important role of doping sites in regulating both the EDL structures and the following electrochemical reactions at the atomic level.Specifically,the negatively charged N atoms repel the anions and adsorb Li~+at high and low potentials,respectively.Such preferential adsorption suggests that Ndoped graphite can promote Li~+desolvation and regulate the location of Li~+deposition.This CPMD method not only unveils the mysterious function of N-doping from the viewpoint of EDL at the atomic level but also applies to probe the interfacial structure on other complicated electrodes.
基金supported by the National Natural Science Foundation of China under Grants Nos.52165013 and 51565021.
文摘The meso-dynamical behaviour of a high-speed rail ballast bed with under sleeper pads(USPs)was studied.The geometrically irregular refined discrete element model of the ballast particles was constructed using 3D scanning techniques,and the 3D dynamic model of the rail-sleeper-ballast bed was constructed using the coupled discrete element method-multiflexible-body dynamics(DEM-MFBD)approach.We analyse the meso-mechanical dynamics of the ballast bed with USPs under dynamic load on a train and verify the correctness of the model in laboratory tests.It is shown that the deformation of the USPs increases the contact area between the sleeper and the ballast particles,and subsequently the number of contacts between them.As the depth of the granular ballast bed increases,the contact area becomes larger,and the contact force between the ballast particles gradually decreases.Under the action of the elastic USPs,the contact forces between ballast particles are reduced and the overall vibration level of the ballast bed can be reduced.The settlement of the granular ballast bed occurs mainly at the shallow position of the sleeper bottom,and the installation of the elastic USPs can be effective in reducing the stress on the ballast particles and the settlement of the ballast bed.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51978336 and 11702117)the Science and Technology Plan Project of Department of Communications of Zhejiang Province(Grant No.2021051)Nantong City Social Livelihood Science and Technology Project(Grant No.MS22022067).
文摘The sloshing in a group of rigid cylindrical tanks with baffles and on soil foundation under horizontal excitation is studied analytically.The solutions for the velocity potential are derived out by the liquid subdomain method.Equivalent models with mass-spring oscillators are established to replace continuous fluid.Combined with the least square technique,Chebyshev polynomials are employed to fit horizontal,rocking and horizontal-rocking coupling impedances of soil,respectively.A lumped parameter model for impedance is presented to describe the effects of soil on tank structures.A mechanical model for the soil-foundation-tank-liquid-baffle system with small amount of calculation and high accuracy is proposed using the substructure technique.The analytical solutions are in comparison with data from reported literature and numerical codes to validate the effectiveness and correctness of the model.Detailed dynamic properties and seismic responses of the soil-tank system are given for the baffle number,size and location as well as soil parameter.
基金supported by the National Natural Science Foundation of China(Grant No.52078010)Beijing Natural Science Foundation(Grant No.JQ19029).
文摘A 5-MW wind turbine has been modeled and analyzed for fluid-structure interaction and aerodynamic performance.In this study, a full-scale model of a 5-MW wind turbine is first developed based on a computational fluid dynamics(CFD) approach, in which the unsteady, noncompressible Reynolds Averaged Navier-Stokes(RANS) method is used. The main focus of the study is to analyze the tower shadow effect on the aerodynamic performance of the wind turbine under different inlet flow conditions. Subsequently, the finite element model is established by considering fluid/structure interactions to study the structural stress, displacement, strain distributions and flow field information of the structure under the uniform wind speed. Finally, the fluid-structure interaction model is established by considering turbulent wind and the tower shadow effect. The variation rules of the dynamic response of the one-way and two-way fluid-structure interaction(FSI) models under different wind speeds are analyzed, and the numerical calculation results are compared with those of the centralized mass model. The results show that the tower shadow effect and structural deformation are the main factors affecting the aerodynamic load fluctuation of the wind turbine, which in turn affects the aerodynamic performance and structural stability of the blades. The structural dynamic response of the coupled model shows significant similarity, while the structural displacement response of the former exhibits less fluctuation compared with the conventional centralized mass model. The one-way fluid-structure interaction(FSI)model shows a higher frequency of stress-strain and displacement oscillations on the blade compared with the two-way FSI model.
基金National Natural Science Foundation of China(U2004176,22008055)Technology Research Project of Henan Province(232102240034)are gratefully acknowledged.
文摘Liquid phase exfoliation(LPE)process for graphene production is usually carried out in stirred tank reactor and the interactions between the solvent and the graphite particles are important as to improve the production efficiency.In this paper,these interactions were revealed by computational fluid dynamics–discrete element method(CFD-DEM)method.Based on simulation results,both liquid phase flow hydrodynamics and particle motion behavior have been analyzed,which gave the general information of the multiphase flow behavior inside the stirred tank reactor as to graphene production.By calculating the threshold at the beginning of graphite exfoliation process,the shear force from the slip velocity was determined as the active force.These results can support the optimization of the graphene production process.
基金the Natural Science Foundation of Jiangsu Province(Grants No.BK20230202)Basic Science(Natural Science)Re-search Project of Colleges and Universities in Jiangsu Province(Grant No.22KJB130005)+3 种基金Changzhou Science and Technology Project(Grant No.CJ20220242)for financial supportJiaqing Kou would like to thank the support of the Alexander von Humboldt Foundation(Ref 3.5-CHN-1227287-HFST-P)Wenkai Yang would like to thank the support of the National Natural Science Foundation of China(Grant No.52205335)supported by Changzhou Sci&Tech Pro-gram(Grant No.CM20223013).
文摘The computational cost of unsteady adjoint equations remains high in adjoint-based unsteady aerodynamic op-timization.In this letter,the solution of unsteady adjoint equations is accelerated by dynamic mode decomposi-tion(DMD).The pseudo-time marching of every real-time step is approximated as an infinite-dimensional linear dynamical system.Thereafter,DMD is utilized to analyze the adjoint vectors sampled from these pseudo-time marching.First-order zero frequency mode is selected to accelerate the pseudo-time marching of unsteady adjoint equations in every real-time step.Through flow past a stationary circular cylinder and an unsteady aerodynamic shape optimization example,the efficiency of solving unsteady adjoint equations is significantly improved.Re-sults show that one hundred adjoint vectors contains enough information about the pseudo-time dynamics,and the adjoint dominant mode can be precisely predicted only by five snapshots produced from the adjoint vectors,which indicates DMD analysis for pseudo-time marching of unsteady adjoint equations is efficient.
文摘We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.
基金the National Natural Science Foundation of China(Grant No.41807223)the Fundamental Research Funds for the Central Universities(Grant No.B210202096)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDA 23090202).
文摘In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.
基金The authors are grateful for the financial support from the National Natural Science Foundation of China(Grant Nos.42177143,42277461)the Science Foundation for Distinguished Young Scholars of Sichuan Province(Grant No.2020JDJQ0011).Thanks to the Chn Energy Dadu River Hydropower Development Co.,Ltd,China Three Gorges Construction Engineering Corporation,Yalong River Hydropower Development Company,Ltd,Power China Chengdu Engineering Co.,Ltd,Power China Northwest Engineering Co.,Ltd,Power China Sinohydro Bureau 7 Co.,Ltd,China Gezhouba Group No.1 Engineering Co.,Ltd.,and the 5th Engineering Co.,Ltd.of China Railway Construction Bridge Engineering Bureau Group for the support and assistance.
文摘The underground hydropower projects in Southwest China is characterized by large excavation sizes,high geostresses,complicated geological conditions and multiple construction processes.Various disasters such as collapses,large deformations,rockbursts are frequently encountered,resulting in serious casualties and huge economic losses.This review mainly presents some representative results on microseismic(MS)monitoring and forecasting for disasters in hydropower underground engineering.First,a set of new denoising,spectral analysis,and location methods were developed for better identification and location of MS signals.Then,the tempo-spatial characteristics of MS events were analyzed to understand the relationship between field construction and damages of surrounding rocks.Combined with field construction,geological data,numerical simulation and parametric analysis of MS sources,the focal mechanism of MS events was revealed.A damage constitutive model considering MS fracturing size was put forward and feedback analysis considering the MS damage of underground surrounding rocks was conducted.Next,an MS multi-parameter based risk assessment and early warning method for dynamic disasters were proposed.The technology for control of the damage and deformation of underground surrounding rocks was proposed for underground caverns.Finally,two typical underground powerhouses were selected as case studies.These achievements can provide significant references for prevention and control of dynamic disasters for underground engineering with similar complicated geological conditions.
基金the National Natural Science Foundation of China(Nos.12232012,12202110,12102191,and 12072159)the Fundamental Research Funds for the Central Universities of China(No.30922010314)the Natural Science Foundation of Guangxi Province of China(No.2020GXNSFBA297010)。
文摘Dynamic coupling modeling and analysis of rotating beams based on the nonlinear Green-Lagrangian strain are introduced in this work.With the reservation of the axial nonlinear strain,there are more coupling terms for axial and transverse deformations.The discretized dynamic governing equations are obtained by using the finite element method and Lagrange’s equations of the second kind.Time responses are conducted to compare the proposed model with other previous models.The stretching deformation due to rotating motion is observed and calculated by special formulations under dynamic equilibrium.The stretching deformation and the change of the associated equilibrium position are taken into account to analyze the free vibration and frequency response of the rotating beams.Analytical and numerical comparisons show that the proposed model can provide reliable results,while the previous models may lead to imprecise results,especially in high-speed conditions.
文摘The solubility of 2,3,4-trichloro-1,5-dinitrobenzene(TCDNB) was measured by a laser dynamic method over the temperature range from 278.15 K to 323.15 K under 0.1 MPa in fifteen mono-solvents(methanol,ethanol,isopropanol,n-butanol,toluene,dichloromethane,chloroform,tetrachloromethane,1,2-dichloroethane,acetone,ethyl acetate,acetonitrile,N-methylpyrrolidone(NMP),N,Ndimethylformamide dimethyl sulfoxide(DMF),dimethyl sulfoxide(DMSO).The solubility of TCDNB could be increased with increasing temperature in fifteen mono-solvents.TCDNB solubility is in the following order at 298.15 K:NMP>DMF>DMSO>toluene>acetone>ethyl acetate>dichloromethane>1,2-dich loroethane>chloroform>acetonitrile>tetrachloromethane>methanol>ethanol>n-butanol>isopropanol.The KAT-LSER model was used to investigate the solvent effect,which revealed that the hydrogen bond acidity of solvents has a greater effect on TCDNB solubility.The van't Hoff model,the modified Apelblat model,theλh model,and the non-random two liquid(NRTL)model were used to correlate the solubility of TCDNB.The calculated solubility data agreed well with the experimental data,and the modified Apelblat model fit best.Furthermore,the van't Hoff and Gibbs equations were also used to calculate the dissolution thermodynamic properties of TCDNB in various solvents.TCDNB dissolution could be an enthalpy-driven,non-spontaneous,and endothermic process in fifteen mono-solvents.The determination and fitting solubility of TCDNB,as well as the calculation of its thermodynamic properties,would be critical in the purification and crystallization of its preparation process research.
基金Project supported by the National Natural Science Foundation of China(No.12172169)the Natural Sciences and Engineering Research Council of Canada(No.NSERC RGPIN-2023-03227)。
文摘We present a study on the dynamic stability of porous functionally graded(PFG)beams under hygro-thermal loading.The variations of the properties of the beams across the beam thicknesses are described by the power-law model.Unlike most studies on this topic,we consider both the bending deformation of the beams and the hygro-thermal load as size-dependent,simultaneously,by adopting the equivalent differential forms of the well-posed nonlocal strain gradient integral theory(NSGIT)which are strictly equipped with a set of constitutive boundary conditions(CBCs),and through which both the stiffness-hardening and stiffness-softening effects of the structures can be observed with the length-scale parameters changed.All the variables presented in the differential problem formulation are discretized.The numerical solution of the dynamic instability region(DIR)of various bounded beams is then developed via the generalized differential quadrature method(GDQM).After verifying the present formulation and results,we examine the effects of different parameters such as the nonlocal/gradient length-scale parameters,the static force factor,the functionally graded(FG)parameter,and the porosity parameter on the DIR.Furthermore,the influence of considering the size-dependent hygro-thermal load is also presented.
文摘Based on the hindcast results of summer rainfall anomalies over China for the period 1981-2000 by the Dynamical Climate Prediction System (IAP-DCP) developed by the Institute of Atmospheric Physics, a correction method that can account for the dependence of model's systematic biases on SST anomalies is proposed. It is shown that this correction method can improve the hindcast skill of the IAP-DCP for summer rainfall anomalies over China, especially in western China and southeast China, which may imply its potential application to real-time seasonal prediction.