The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can ov...The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.展开更多
The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conf...The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.展开更多
In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a ...In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a generalization of the formaulation of Liu[4-6] for the thin film hydrodynamic lubrication re- gime is presented. This method is then aplied to a strip rolling simulation,using the Arbitrary La- grangian eulerian (ALE) formalism.展开更多
Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its e...Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its effects. The goal of this paper is to get results of stresses and displacements of a coffee harvester structure by using FEM for static simulation. The main parts of the coffee harvester analyzed were: engine frame, body right and left sides, front and rear end, main beam, coffee reservoir, wheels and fuel tank. Two different design concepts of a coffee harvester machine were analyzed (structure with rear wheels aligned and misaligned) and the results were compared. It was observed that the model with rear wheels misaligned showed maximum displacement lower than the model with rear wheels aligned. Although higher stress was found in the rear wheels misaligned, it was observed that average stresses for the misaligned wheels design were lower in most structural components analyzed. Based on FEM results, the coffee harvester machine with misaligned rear wheels was built and subjected to operational tests without showing any structural failure.展开更多
On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node fi...On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.展开更多
A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements f...A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.展开更多
Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular character...Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.展开更多
This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families l...This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.展开更多
This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left an...This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left and right limits for the discontinuous finite element at nodes) has the optimal-order ultraconvergence 2k + 2. For nanlinear Hamiltonian systems (e.g., SchrSdinger equation and Kepler system) with momentum conservation, the discontinuous finite element methods preserve momentum at nodes. These properties are confirmed by numerical experiments.展开更多
In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discr...In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.展开更多
Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 n...Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.展开更多
The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the tw...The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.展开更多
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.展开更多
In this paper,we introduce a stabilizer free weak Galerkin(SFWG)finite element method for second order elliptic problems on rectangular meshes.With a special weak Gradient space,an order two superconvergence for the S...In this paper,we introduce a stabilizer free weak Galerkin(SFWG)finite element method for second order elliptic problems on rectangular meshes.With a special weak Gradient space,an order two superconvergence for the SFWG finite element solution is obtained,in both L 2 and H1 norms.A local post-process lifts such a Pk weak Galerkin solution to an optimal order Pk+2 solution.The numerical results confirm the theory.展开更多
For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of ...For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.展开更多
This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and ...This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.展开更多
In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hy...In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale pr...Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.展开更多
The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling ...The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.展开更多
文摘The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.
文摘The family of Falk-Neilan P_(k)finite elements,combined with the Argyris P_(k+1)finite elements,solves the Reissner-Mindlin plate equation quasi-optimally and locking-free,on triangular meshes.The method is truly conforming or consistent in the sense that no projection/reduction is introduced.Theoretical proof and numerical confirmation are presented.
文摘In this paper, the lubrication problem in nmmerical sindation of forming processes is considered.After enumerationg the difficulties encountered when trying to solve such a problem with the finite element method, a generalization of the formaulation of Liu[4-6] for the thin film hydrodynamic lubrication re- gime is presented. This method is then aplied to a strip rolling simulation,using the Arbitrary La- grangian eulerian (ALE) formalism.
文摘Stress concentration and large displacements are usual problems in the components of the structure of agricultural machinery such harvesters coffee, and that finite element method (FEM) can be a tool to minimize its effects. The goal of this paper is to get results of stresses and displacements of a coffee harvester structure by using FEM for static simulation. The main parts of the coffee harvester analyzed were: engine frame, body right and left sides, front and rear end, main beam, coffee reservoir, wheels and fuel tank. Two different design concepts of a coffee harvester machine were analyzed (structure with rear wheels aligned and misaligned) and the results were compared. It was observed that the model with rear wheels misaligned showed maximum displacement lower than the model with rear wheels aligned. Although higher stress was found in the rear wheels misaligned, it was observed that average stresses for the misaligned wheels design were lower in most structural components analyzed. Based on FEM results, the coffee harvester machine with misaligned rear wheels was built and subjected to operational tests without showing any structural failure.
文摘On the basis of the one-dimension infinite element theory, the coordinate translation and shape function of 3D point-radiate 8-node and 4-node infinite elements are derived. They are coupled with 20-node and 8-node finite elements to compute the compression distortion of the prestressed anchorage segment. The results indicate that when the prestressed force acts on the anchorage head and segment, the stresses and the displacements in the rock around the anchorage head and segment concentrate on the zone center with the anchor axis, and they decrease with exponential forms. Therefore,the stresses and the displacement spindles are formed. The calculating results of the infinite element are close to the theoretical results. This indicates the method is right. This article introduces a new way to study the mechanism of prestressed anchors. The obtained results have an important role in the research of the anchor mechanism and engineering application.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.
文摘Materials with a negative Poisson's ratio(PR)are called auxetics;they are characterized by expansion/contraction when tensioned/compressed.Given this counterintuitive behavior,they present very particular characteristics and mechanical behavior.Geometrical models have been developed to justify and artificiall reproduce such materials' auxetic behavior.The focus of this study is the exploration of a reentrant model by analyzing the variation in the PR of reentrant structures as a function of geometrical and base material parameters.It is shown that,even in the presence of protruding ribs,there may not be auxetic behavior,and this depends on the geometry of each reentrant structure.Values determined for these parameters can be helpful as approximate reference data in the design and fabrication of auxetic lattices using reentrant geometries.
文摘This paper presents the formulation of finite elements based on Deslauriers-Dubuc interpolating scaling functions, also known as Interpolets, for their use in wave propagation modeling. Unlike other wavelet families like Daubechies, Interpolets possess rational filter coefficients, are smooth, symmetric and therefore more suitable for use in numerical methods. Expressions for stiffness and mass matrices are developed based on connection coefficients, which are inner products of basis functions and their derivatives. An example in 1-D was formulated using Central Difference and Newmark schemes for time differentiation. Encouraging results were obtained even for large time steps. Results obtained in 2-D are compared with the standard Finite Difference Method for validation.
基金Project supported by the National Natural Science Foundation of China (No. 10771063)
文摘This paper discusses the k-degree averaging discontinuous finite element solution for the initial value problem of ordinary differential equations. When k is even, the averaging numerical flux (the average of left and right limits for the discontinuous finite element at nodes) has the optimal-order ultraconvergence 2k + 2. For nanlinear Hamiltonian systems (e.g., SchrSdinger equation and Kepler system) with momentum conservation, the discontinuous finite element methods preserve momentum at nodes. These properties are confirmed by numerical experiments.
文摘In this paper, we study mixed finite elements for parabolic integro-differential equations, and introduce a kind of nonclassical mixed projection, its optimal L-2 and h(-s) estimates are obtained. We define semi-discrete and full-discrete mixed finite elements for the equations, and obtain the optimal L-2 error estimates.
基金Project supported by the National Natural Science Foundation of China (Nos. 10771198 and 11071226)the Foundation of International Science and Technology Cooperation of Henan Province
文摘Two new locking-free nonconforming finite elements for the pure displacement planar elasticity problem are presented. Convergence rates of the elements are uniformly optimal with respect to A. The energy norm and L2 norm errors are proved to be O(h2) and O(h3), respectively. Numerical tests confirm the theoretical analysis.
文摘The paper presents our contribution to the full 3D finite element modelling of a hybrid stepping motor using COMSOL Multiphysics software. This type of four-phase motor has a permanent magnet interposed between the two identical and coaxial half stators. The calculation of the field with or without current in the windings (respectively with or without permanent magnet) is done using a mixed formulation with strong coupling. In addition, the local high saturation of the ferromagnetic material and the radial and axial components of the magnetic flux are taken into account. The results obtained make it possible to clearly observe, as a function of the intensity of the bus current or the remanent induction, the saturation zones, the lines, the orientations and the magnetic flux densities. 3D finite element modelling provide more accurate numerical data on the magnetic field through multiphysics analysis. This analysis considers the actual operating conditions and leads to the design of an optimized machine structure, with or without current in the windings and/or permanent magnet.
基金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.
基金Xiu Ye was supported in part by National Science Foundation Grant DMS-1620016.
文摘In this paper,we introduce a stabilizer free weak Galerkin(SFWG)finite element method for second order elliptic problems on rectangular meshes.With a special weak Gradient space,an order two superconvergence for the SFWG finite element solution is obtained,in both L 2 and H1 norms.A local post-process lifts such a Pk weak Galerkin solution to an optimal order Pk+2 solution.The numerical results confirm the theory.
基金supported by National Natural Science Foundation of China(11771257)the Shandong Provincial Natural Science Foundation of China(ZR2023YQ002,ZR2023MA007,ZR2021MA004)。
文摘For singularly perturbed convection-diffusion problems,supercloseness analysis of the finite element method is still open on Bakhvalov-type meshes,especially in the case of 2D.The difficulties arise from the width of the mesh in the layer adjacent to the transition point,resulting in a suboptimal estimate for convergence.Existing analysis techniques cannot handle these difficulties well.To fill this gap,here a novel interpolation is designed delicately for the smooth part of the solution,bringing about the optimal supercloseness result of almost order 2 under an energy norm for the finite element method.Our theoretical result is uniform in the singular perturbation parameterεand is supported by the numerical experiments.
基金Project supported by the National Natural Science Foundation of China(Nos.12072209,U21A2043012192211)+1 种基金the Natural Science Foundation of Hebei Province of China(No.A2020210009)the S&T Program of Hebei Province of China(No.225676162GH)。
文摘This paper theoretically studies the axisymmetric frictionless indentation of a transversely isotropic piezoelectric semiconductor(PSC)half-space subject to a rigid flatended cylindrical indenter.The contact area and other surface of the PSC half-space are assumed to be electrically insulating.By the Hankel integral transformation,the problem is reduced to the Fredholm integral equation of the second kind.This equation is solved numerically to obtain the indentation behaviors of the PSC half-space,mainly including the indentation force-depth relation and the electric potential-depth relation.The results show that the effect of the semiconductor property on the indentation responses is limited within a certain range of variation of the steady carrier concentration.The dependence of indentation behavior on material properties is also analyzed by two different kinds of PSCs.Finite element simulations are conducted to verify the results calculated by the integral equation technique,and good agreement is demonstrated.
文摘In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金Supported by Science Center for Gas Turbine Project of China (Grant No.P2022-B-IV-014-001)Frontier Leading Technology Basic Research Special Project of Jiangsu Province of China (Grant No.BK20212007)the BIT Research and Innovation Promoting Project of China (Grant No.2022YCXZ019)。
文摘Thermal conductivity is one of the most significant criterion of three-dimensional carbon fiber-reinforced SiC matrix composites(3D C/SiC).Represent volume element(RVE)models of microscale,void/matrix and mesoscale proposed in this work are used to simulate the thermal conductivity behaviors of the 3D C/SiC composites.An entirely new process is introduced to weave the preform with three-dimensional orthogonal architecture.The 3D steady-state analysis step is created for assessing the thermal conductivity behaviors of the composites by applying periodic temperature boundary conditions.Three RVE models of cuboid,hexagonal and fiber random distribution are respectively developed to comparatively study the influence of fiber package pattern on the thermal conductivities at the microscale.Besides,the effect of void morphology on the thermal conductivity of the matrix is analyzed by the void/matrix models.The prediction results at the mesoscale correspond closely to the experimental values.The effect of the porosities and fiber volume fractions on the thermal conductivities is also taken into consideration.The multi-scale models mentioned in this paper can be used to predict the thermal conductivity behaviors of other composites with complex structures.
基金The Construction S&T Project of the Department of Transportation of Sichuan Province(Grant No.2023A02)the National Natural Science Foundation of China(No.52109135).
文摘The surrounding geological conditions and supporting structures of underground engineering are often updated during construction,and these updates require repeated numerical modeling.To improve the numerical modeling efficiency of underground engineering,a modularized and parametric modeling cloud server is developed by using Python codes.The basic framework of the cloud server is as follows:input the modeling parameters into the web platform,implement Rhino software and FLAC3D software to model and run simulations in the cloud server,and return the simulation results to the web platform.The modeling program can automatically generate instructions that can run the modeling process in Rhino based on the input modeling parameters.The main modules of the modeling program include modeling the 3D geological structures,the underground engineering structures,and the supporting structures as well as meshing the geometric models.In particular,various cross-sections of underground caverns are crafted as parametricmodules in themodeling program.Themodularized and parametric modeling program is used for a finite element simulation of the underground powerhouse of the Shuangjiangkou Hydropower Station.This complicatedmodel is rapidly generated for the simulation,and the simulation results are reasonable.Thus,this modularized and parametric modeling program is applicable for three-dimensional finite element simulations and analyses.