The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is m...The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.展开更多
Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymme...Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymmetric tilt boundaries are hard to be accounted for based on traditional theoretical models,and the corresponding solute segregation is complex.Herein,atomic structures of a specific asymmetric boundary on{1012}TBs were reveled using aberration-corrected high-angle annular dark-field scanning transmission electron microscopy(HAADF-STEM),molecular dynamics(MD)and density functional theory(DFT)simulations.Reaction between<a60>M dislocations and the{1012}TB can generate a~61°/25°asymmetric tilt boundary.The segregation of Gd and Zn atoms is closely related to the aggregateddislocations and the interfacial interstices of the asymmetric tilt boundary,which is energetically favorable in reducing the total system energy.展开更多
To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and ga...To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.展开更多
This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic...This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.展开更多
The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of...The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.展开更多
In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong co...In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.展开更多
Second period elements(B,C,N,and O)usually appear at the grain boundary(GB)and strongly affect the mechanical performance in austenitic stainless steels.Therefore,it is significant to investigate the effect of solute ...Second period elements(B,C,N,and O)usually appear at the grain boundary(GB)and strongly affect the mechanical performance in austenitic stainless steels.Therefore,it is significant to investigate the effect of solute elements(B,C,N,and O)on the GB.The first-principles calculation based on the density function theory was applied to explore the effect of B,C,N,and O onγ-FeΣ5(210)[001]GB.The GB energy,the segregation energy,the Voronoi volume,and the theoretical tensile test were calculated to investigate the segregation behavior and the strengthening effect.The structural change and electronic evolution were also investigated by bond change,charge density distribution,and density of states.The results show that B is favored to segregate at the capped trigonal prism(CTP)position with a large void and has a strengthening effect on the GB strength,while O and N are preferred to locate at the octahedral(OCT)site and have an embrittling effect on GB cohesion.C can segregate at both the CTP site and the OCT location with little energy difference.As C segregates at the OCT site,it is beneficial for GB strength.However,it is detrimental at the CTP position.It can be seen that the influence of solutes is closely related to the element type and segregated position.展开更多
THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between c...THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between crystallographic orientation,grain boundary(GB)proximity,and pore characteristics(size/location).This study compares single-crystal nickel models along[100],[110],and[111]orientations with equiaxed polycrystalline models containing 0,1,and 2.5 nm pores in surface and subsurface configurations.Our results reveal that crystallographic anisotropy manifests as a 24.4%higher elastic modulus and 22.2%greater hardness in[111]-oriented single crystals compared to[100].Pore-GB synergistic effects are found to dominate the deformation behavior:2.5 nm subsurface pores reduce hardness by 25.2%through stress concentration and dislocation annihilation at GBs,whereas surface pores enable mechanical recovery via accelerated dislocation generation post-collapse.Additionally,size-dependent deformation regimes were identified,with 1 nm pores inducing negligible perturbation due to rapid atomic rearrangement,in contrast with persistent softening in 2.5 nm pores.These findings establish atomic-scale design principles for defect engineering in nickel-based aerospace components,demonstrating how crystallographic orientation,pore configuration,and GB interactions collectively govern nanoindentation behavior.展开更多
This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the pred...This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.展开更多
On the stone-paved lanes of Songyang County that date back to ancient times,morning mist lingered as a faint fragrance of tea wafted from a century-old house.Inside,Yang Junjie,a tea maker born in the 1980s,worked def...On the stone-paved lanes of Songyang County that date back to ancient times,morning mist lingered as a faint fragrance of tea wafted from a century-old house.Inside,Yang Junjie,a tea maker born in the 1980s,worked deftly at the stove,his hands moving swiftly over the scorching iron wok as tender green tea leaves dance between his fingers.展开更多
We report the results of the experiment on synthesizing ^(287,288)Mc isotopes (Z=115) using the fusionevaporation reaction ^(243)Am(^(48)Ca,4n,3n)^(287,288)Mc at the Spectrometer for Heavy Atoms and Nuclear Structure-...We report the results of the experiment on synthesizing ^(287,288)Mc isotopes (Z=115) using the fusionevaporation reaction ^(243)Am(^(48)Ca,4n,3n)^(287,288)Mc at the Spectrometer for Heavy Atoms and Nuclear Structure-2(SHANS2),a gas-filled recoil separator located at the China Accelerator Facility for Superheavy Elements(CAFE2).In total,20 decay chains are attributed to ^(288)Mc and 1 decay chain is assigned to ^(287)Mc.The measured oa-decay properties of ^(287,288)Mc as well as its descendants are consistent with the known data.No additional decay chains originating from the 2n or 5n reaction channels were detected.The excitation function of the ^(243)Am(^(48)Ca,3n)^(288)Mc reaction was measured at the cross-section level of picobarn,which indicates the promising capability for the study of heavy and superheavy nuclei at the facility.展开更多
In conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete bounda...In conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete boundary unit cells(i.e.,boundary defects)even though the bulk polarization is zero,which challenges the conventional understanding of HOTIs.Here,based on a Kekul´e-distorted honeycomb lattice with incomplete unit cells,we reveal that incomplete unit cells exhibit fractional charges through the analysis of Wannier centers by developing a compensation method and creating the concept of Wannier center domain(WCD)which is the smallest region that one Wannier center occupies.This method compensates for the missing parts of these boundary incomplete unit cells with additional WCDs to make them complete.The compensated WCDs automatically carry the corresponding charge,and this charge together with that of the incomplete unit cell constitutes the total charge of the complete unit cell after compensation.We conclude that the emergence of corner states is attributed to the filling anomaly,which is a fundamental mechanism.Our results refresh the understanding of HOTIs,especially those with structural discontinuities,and provide a novel design for topological states which have application value in producing optical functional devices.展开更多
The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analy...The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.展开更多
The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-eleme...The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.展开更多
In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the ...In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.展开更多
The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach...The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables,considering the yield limit,rotation speeds and material density to be the fundamental stochastic variables.Through analyzing a numerical example and a turbo-disk of an aeroengine,the results show that the method developed is successful.展开更多
The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the sin...The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.展开更多
A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate...A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate the convergence to FBMD method by means of introducing the residual modes which take into accaunt the quasi -state contributiort of all neglected modes. As an example, the vibration and sound radiation of a submerged spherical shell excited by axisymmetric force are studied in cases of ka=l,2,3 and 4. From the calculated results we see that the FBMMD method shows a significant improvement to the accuracy of surface sound pressure, normal displacement and directivity patterns of radiating sound, especially to the directivity patterns.展开更多
In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach num...In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.展开更多
基金supported by the Science and Technology Research Project of Henan Province(242102241055)the Industry-University-Research Collaborative Innovation Base on Automobile Lightweight of“Science and Technology Innovation in Central Plains”(2024KCZY315)the Opening Fund of State Key Laboratory of Structural Analysis,Optimization and CAE Software for Industrial Equipment(GZ2024A03-ZZU).
文摘The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.
基金supported by the Scientific and Technological Developing Scheme of Jilin Province under grants no.YDZJ202301ZYTS538the Chinese Academy of Sciences Youth Innovation Promotion Association under grants number 2023234+3 种基金the National Natural Science Foundation of China under grants number U21A20323the Scientific and Technological Developing Scheme of Jilin Province under grants no.SKL202302038the Major Scientific and Technological Projects of Hebei Province under grants No.23291001Zthe Scientific and Technology Project of Hanjiang District.
文摘Asymmetric tilt boundaries on conventional twin boundaries(TBs)are significant for understanding the role of twins on coordinating plastic deformation in many metallic alloys.However,the formation modes of many asymmetric tilt boundaries are hard to be accounted for based on traditional theoretical models,and the corresponding solute segregation is complex.Herein,atomic structures of a specific asymmetric boundary on{1012}TBs were reveled using aberration-corrected high-angle annular dark-field scanning transmission electron microscopy(HAADF-STEM),molecular dynamics(MD)and density functional theory(DFT)simulations.Reaction between<a60>M dislocations and the{1012}TB can generate a~61°/25°asymmetric tilt boundary.The segregation of Gd and Zn atoms is closely related to the aggregateddislocations and the interfacial interstices of the asymmetric tilt boundary,which is energetically favorable in reducing the total system energy.
基金Supported by the China National Natural Science Foundation Project(52274048)Beijing Natural Science Foundation Project(3222037)。
文摘To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.
基金supported by the Shanxi Scholarship Council of China(Grant No.2023-036)the Natural Science Foundation of Shanxi Province(Grant No.202303021222020).
文摘This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.
基金funded by research grant OPUS no.2021/41/B/ST8/02432 entitled Probabilistic entropy in engineering computations sponsored by The National Science Center in Polandthe Institute of Structural Analysis of Poznan University of Technology in the framework of the internal research grant 0411/SBAD/0010.
文摘The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.
基金supported by National Scientific Foundation of China with Grant No. 41774067
文摘In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.
基金supported by National Key R&D Program of China(No.2022YFB3705202)National Natural Science Foundation of China(Nos.51831008,52171049 and 52104330).
文摘Second period elements(B,C,N,and O)usually appear at the grain boundary(GB)and strongly affect the mechanical performance in austenitic stainless steels.Therefore,it is significant to investigate the effect of solute elements(B,C,N,and O)on the GB.The first-principles calculation based on the density function theory was applied to explore the effect of B,C,N,and O onγ-FeΣ5(210)[001]GB.The GB energy,the segregation energy,the Voronoi volume,and the theoretical tensile test were calculated to investigate the segregation behavior and the strengthening effect.The structural change and electronic evolution were also investigated by bond change,charge density distribution,and density of states.The results show that B is favored to segregate at the capped trigonal prism(CTP)position with a large void and has a strengthening effect on the GB strength,while O and N are preferred to locate at the octahedral(OCT)site and have an embrittling effect on GB cohesion.C can segregate at both the CTP site and the OCT location with little energy difference.As C segregates at the OCT site,it is beneficial for GB strength.However,it is detrimental at the CTP position.It can be seen that the influence of solutes is closely related to the element type and segregated position.
基金The National Natural Science Foundation of China(Grant No.12462006)Beijing Institute of Structure and Environment Engineering Joint Innovation Fund(No.BQJJ202414).
文摘THE mechanical response and deformation mechanisms of pure nickel under nanoindentation were systematically investigated using molecular dynamics(MD)simulations,with a particular focus on the novel interplay between crystallographic orientation,grain boundary(GB)proximity,and pore characteristics(size/location).This study compares single-crystal nickel models along[100],[110],and[111]orientations with equiaxed polycrystalline models containing 0,1,and 2.5 nm pores in surface and subsurface configurations.Our results reveal that crystallographic anisotropy manifests as a 24.4%higher elastic modulus and 22.2%greater hardness in[111]-oriented single crystals compared to[100].Pore-GB synergistic effects are found to dominate the deformation behavior:2.5 nm subsurface pores reduce hardness by 25.2%through stress concentration and dislocation annihilation at GBs,whereas surface pores enable mechanical recovery via accelerated dislocation generation post-collapse.Additionally,size-dependent deformation regimes were identified,with 1 nm pores inducing negligible perturbation due to rapid atomic rearrangement,in contrast with persistent softening in 2.5 nm pores.These findings establish atomic-scale design principles for defect engineering in nickel-based aerospace components,demonstrating how crystallographic orientation,pore configuration,and GB interactions collectively govern nanoindentation behavior.
基金supported by the National Key R&D Program of China[grant number 2023YFC3008004]。
文摘This study introduces a new ocean surface friction velocity scheme and a modified Thompson cloud microphysics parameterization scheme into the CMA-TYM model.The impact of these two parameterization schemes on the prediction of the movement track and intensity of Typhoon Kompasu in 2021 is examined.Additionally,the possible reasons for their effects on tropical cyclone(TC)intensity prediction are analyzed.Statistical results show that both parameterization schemes improve the predictions of Typhoon Kompasu’s track and intensity.The influence on track prediction becomes evident after 60 h of model integration,while the significant positive impact on intensity prediction is observed after 66 h.Further analysis reveals that these two schemes affect the timing and magnitude of extreme TC intensity values by influencing the evolution of the TC’s warm-core structure.
文摘On the stone-paved lanes of Songyang County that date back to ancient times,morning mist lingered as a faint fragrance of tea wafted from a century-old house.Inside,Yang Junjie,a tea maker born in the 1980s,worked deftly at the stove,his hands moving swiftly over the scorching iron wok as tender green tea leaves dance between his fingers.
基金supported in part by the National Key R&D Program of China (Contract Nos.2023YFA1606500,2024YFE0109800,and 2024YFE0110400)Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB34010000)+5 种基金the Gansu Key Project of Science and Technology (Grant No.23ZDGA014)the Guangdong Major Project of Basic and Applied Basic Research (Grant No.2021B0301030006)the National Natural Science Foundation of China (Grant Nos.12105328,W2412040,12475126,12422507,12035011,12375118,12435008,and W2412043)the Chinese Academy of Sciences Project for Young Scientists in Basic Research(Grant No.YSBR-002)the Youth Innovation Promotion Association of the Chinese Academy of Sciences (Grant Nos.2020409 and 2023439)the Russian Science Foundation (Grant No.25-42-00003)。
文摘We report the results of the experiment on synthesizing ^(287,288)Mc isotopes (Z=115) using the fusionevaporation reaction ^(243)Am(^(48)Ca,4n,3n)^(287,288)Mc at the Spectrometer for Heavy Atoms and Nuclear Structure-2(SHANS2),a gas-filled recoil separator located at the China Accelerator Facility for Superheavy Elements(CAFE2).In total,20 decay chains are attributed to ^(288)Mc and 1 decay chain is assigned to ^(287)Mc.The measured oa-decay properties of ^(287,288)Mc as well as its descendants are consistent with the known data.No additional decay chains originating from the 2n or 5n reaction channels were detected.The excitation function of the ^(243)Am(^(48)Ca,3n)^(288)Mc reaction was measured at the cross-section level of picobarn,which indicates the promising capability for the study of heavy and superheavy nuclei at the facility.
基金supported by the Natural Science Basic Research Program of Shaanxi Province (Grant Nos.2024JC-JCQN-06 and2025JC-QYCX-006)the National Natural Science Foundation of China (Grant No.12474337)Chinese Academy of Sciences Project (Grant Nos.E4BA270100,E4Z127010F,E4Z6270100,and E53327020D)。
文摘In conventional higher-order topological insulators(HOTIs),the emergence of topological states can be explained by using the nonzero bulk polarization index.However,corner states emerge in HOTIs with incomplete boundary unit cells(i.e.,boundary defects)even though the bulk polarization is zero,which challenges the conventional understanding of HOTIs.Here,based on a Kekul´e-distorted honeycomb lattice with incomplete unit cells,we reveal that incomplete unit cells exhibit fractional charges through the analysis of Wannier centers by developing a compensation method and creating the concept of Wannier center domain(WCD)which is the smallest region that one Wannier center occupies.This method compensates for the missing parts of these boundary incomplete unit cells with additional WCDs to make them complete.The compensated WCDs automatically carry the corresponding charge,and this charge together with that of the incomplete unit cell constitutes the total charge of the complete unit cell after compensation.We conclude that the emergence of corner states is attributed to the filling anomaly,which is a fundamental mechanism.Our results refresh the understanding of HOTIs,especially those with structural discontinuities,and provide a novel design for topological states which have application value in producing optical functional devices.
文摘The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.
基金sponsored by the National Natural Science Foundation of China Research(Grant No.41274138)the Science Foundation of China University of Petroleum(Beijing)(No.KYJJ2012-05-02)
文摘The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.
文摘In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.
文摘The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables,considering the yield limit,rotation speeds and material density to be the fundamental stochastic variables.Through analyzing a numerical example and a turbo-disk of an aeroengine,the results show that the method developed is successful.
文摘The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.
文摘A finite element / boundary element-modified modal decomposition method (FBMMD) is presented for predicting the vibration and sound radiation from submerged shell of revolution. Improvement has been made to accelerate the convergence to FBMD method by means of introducing the residual modes which take into accaunt the quasi -state contributiort of all neglected modes. As an example, the vibration and sound radiation of a submerged spherical shell excited by axisymmetric force are studied in cases of ka=l,2,3 and 4. From the calculated results we see that the FBMMD method shows a significant improvement to the accuracy of surface sound pressure, normal displacement and directivity patterns of radiating sound, especially to the directivity patterns.
基金the China Scholarship Council(CSC)the Aeronautics Science Foundation of China(No.20163252037)+2 种基金the China Postdoctoral Science Foundation(No.2017M610325)the Natural Science Foundation of Jiangsu Province of China(No.BK20170771)the Fundamental Research Funds for the Central Universities of China(No.NP2017202)for their support
文摘In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.
