This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulati...Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.展开更多
Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily sim...Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily simulated because conventional PINNs fail in discontinuous heterogeneous porous media or multi-layer strata when labeled data are missing.This work aims to develop a universal network structure to encode the mass continuity equation and Darcy’s law without labeled data.The finite element approximation,which can decompose a complex heterogeneous domain into simpler ones,is adopted to build the differentiable network.Without conventional DNNs,physics-encoded finite element network(PEFEN)can avoid spectral bias and learn high-frequency functions with sharp/steep gradients.PEFEN rigorously encodes Dirichlet and Neumann boundary conditions without training.Benefiting from its discretized formulation,the discontinuous heterogeneous hydraulic conductivity is readily embedded into the network.Three typical cases are reproduced to corroborate PEFEN’s superior performance over conventional PINNs and the PINN with mixed formulation.PEFEN is sparse and demonstrated to be capable of dealing with heterogeneity with much fewer training iterations(less than 1/30)than the improved PINN with mixed formulation.Thus,PEFEN saves energy and contributes to low-carbon AI for science.The last two cases focus on common geotechnical settings of impermeable sheet pile in singlelayer and multi-layer strata.PEFEN solves these cases with high accuracy,circumventing costly labeled data,extra computational burden,and additional treatment.Thus,this study warrants the further development and application of PEFEN as a novel differentiable network in porous flow of practical geotechnical engineering.展开更多
A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinfo...A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinforced composites(FRCs).The fibers are periodically distributed and unidirectionally aligned in a homogeneous matrix.This framework addresses the static linear elastic micropolar problem through partial differential equations,subject to boundary conditions and perfect interface contact conditions.The mathematical formulation of the local problems and the effective coefficients are presented by the AHM.The local problems obtained from the AHM are solved by the FEM,which is denoted as the SAFEM.The numerical results are provided,and the accuracy of the solutions is analyzed,indicating that the formulas and results obtained with the SAFEM may serve as the reference points for validating the outcomes of experimental and numerical computations.展开更多
Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Eul...Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.展开更多
Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computation...Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.展开更多
Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation a...Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method（FEM） model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters＇ impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.展开更多
The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing matu...The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s^(-1),the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.展开更多
Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,compara...Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.展开更多
Finite element method(FEM) was used to simulate the forming process of shotpeening the wing skin panel. Experiment of shotpeeing the wing skin panel was carried out. The results show that equivalent deformation in sho...Finite element method(FEM) was used to simulate the forming process of shotpeening the wing skin panel. Experiment of shotpeeing the wing skin panel was carried out. The results show that equivalent deformation in shotpeening process can be obtained using the elongation and bending result caused by thermal stress that is induced by applying temperature load on the surface of the part. Deformation of the part in the shotpeeing process can be analyzed using this method. The parameters and their relationships are identified.展开更多
Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid ...Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.展开更多
Simulation of the temperature field of copier paper in copier fusing is very important for improving the fusing property of reprography. The temperature field of copier paper varies with a high gradient when the copie...Simulation of the temperature field of copier paper in copier fusing is very important for improving the fusing property of reprography. The temperature field of copier paper varies with a high gradient when the copier paper is moving through the fusing rollers. By means of conventional shaft elements, the high gradient temperature variety causes the oscillation of the numerical solution. Based on the Daubechies scaling functions, a kind of wavelet based element is constructed for the above problem. The temperature field of the copier paper moving through the fusing rollers is simulated using the two methods. Comparison of the results shows the advantages of the wavelet finite element method, which provides a new method for improving the copier properties.展开更多
In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to pred...In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS) model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.展开更多
This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and...This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.展开更多
A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid techn...A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.展开更多
Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if ...Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if the material contains a large number of voids. Eight node isoparametric elastic plastic finite element method was used to simulate the tensile process of sintered powder material. By setting a number of voids in the analyzed metal cuboid, the initial density was taken into consideration. The material properties of the three dimensional solid for the tensile simulation were defined with reference to the known pure iron material parameters. The load displacement curves during elongation were obtained with a universal testing machine, and then the simulated curves were compared with the experimental results. The factors that cause the stress concentration and strength decrease were analyzed according to the simulated equivalent von Mises stress distribution.展开更多
The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear probl...The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.展开更多
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金financially supported by the Riset Pengembangan dan Penerapan(RPP),Diponegoro University 2023 research scheme with contract number 609-18/UN7.D2/PP/VIII/2023.
文摘Shallow water infrastructure needs to support increased activity on the shores of Semarang.This study chooses several pontoons because of their good stability,rolling motion,and more expansive space.A coupled simulation method consisting of hydrodynamic and structural calculations has been used to evaluate a catamaran pontoon’s motion and structural integrity.Four different space sizes are set for the pontoon system:5 m,5.5 m,6 m,and 6.5 m.The frequency domain shows that the pontoon space affects the RAO in wave periods ranging from 3 s to 5 s.At wave periods of 3 s,4 s,and 5 s,the pontoon space significantly affects the maximum motion and chain tension parameter values,which are evaluated via time domain simulation.The critical stress of the pontoon is shown at a wave period of 5 s for 5 m and 5.5 m of pontoon space,which shows that the stress can reach 248 MPa.
基金supported by the National Natural Science Foundation of China(Grant Nos.42272338 and 41827807)Department of Transportation of Zhejiang Province,China(Grant No.202213).
文摘Physics-informed neural networks(PINNs)have prevailed as differentiable simulators to investigate flow in porous media.Despite recent progress PINNs have achieved,practical geotechnical scenarios cannot be readily simulated because conventional PINNs fail in discontinuous heterogeneous porous media or multi-layer strata when labeled data are missing.This work aims to develop a universal network structure to encode the mass continuity equation and Darcy’s law without labeled data.The finite element approximation,which can decompose a complex heterogeneous domain into simpler ones,is adopted to build the differentiable network.Without conventional DNNs,physics-encoded finite element network(PEFEN)can avoid spectral bias and learn high-frequency functions with sharp/steep gradients.PEFEN rigorously encodes Dirichlet and Neumann boundary conditions without training.Benefiting from its discretized formulation,the discontinuous heterogeneous hydraulic conductivity is readily embedded into the network.Three typical cases are reproduced to corroborate PEFEN’s superior performance over conventional PINNs and the PINN with mixed formulation.PEFEN is sparse and demonstrated to be capable of dealing with heterogeneity with much fewer training iterations(less than 1/30)than the improved PINN with mixed formulation.Thus,PEFEN saves energy and contributes to low-carbon AI for science.The last two cases focus on common geotechnical settings of impermeable sheet pile in singlelayer and multi-layer strata.PEFEN solves these cases with high accuracy,circumventing costly labeled data,extra computational burden,and additional treatment.Thus,this study warrants the further development and application of PEFEN as a novel differentiable network in porous flow of practical geotechnical engineering.
基金Project supported by the National Council of Humanities,Sciences,and Technologies of Mexico(Nos.CF-2023-G-792 and CF-2023-G-1458)the National Council for Scientific and Technological Development of Brazil(No.09/2023)the Research on Productivity of Brazil(No.307188/2023-0)。
文摘A semi-analytical finite element method(SAFEM),based on the two-scale asymptotic homogenization method(AHM)and the finite element method(FEM),is implemented to obtain the effective properties of two-phase fiber-reinforced composites(FRCs).The fibers are periodically distributed and unidirectionally aligned in a homogeneous matrix.This framework addresses the static linear elastic micropolar problem through partial differential equations,subject to boundary conditions and perfect interface contact conditions.The mathematical formulation of the local problems and the effective coefficients are presented by the AHM.The local problems obtained from the AHM are solved by the FEM,which is denoted as the SAFEM.The numerical results are provided,and the accuracy of the solutions is analyzed,indicating that the formulas and results obtained with the SAFEM may serve as the reference points for validating the outcomes of experimental and numerical computations.
文摘Viscoelastic foundation plays a very important role in civil engineering. It can effectively disperse the structural load into the foundation soil and avoid the damage caused by the concentrated load. The model of Euler-Bernoulli beam on viscoelastic Pasternak foundation can be used to analyze the deformation and response of buildings under complex geological conditions. In this paper, we use Hermite finite element method to get the numerical approximation scheme for the vibration equation of viscoelastic Pasternak foundation beam. Convergence and error estimation are rigourously established. We prove that the fully discrete scheme has convergence order O(τ2+h4), where τis time step size and his space step size. Finally, we give four numerical examples to verify the validity of theoretical analysis.
文摘Recently,the finite element method(FEM) has been commonly applied in the engineering analysis of rotor dynamics.Gyroscopic moments,rotary inertia,transverse shear deformation and gravity can be included in computational models of rotor-bearing systems.In this paper,a finite element model and its solution method are presented for the calculation of the dynamics of dual rotor systems.A typical structure with two rotor shafts is discussed and the procedure for obtaining the coupling motion equations of the subsystems is illustrated.A computer program is developed to solve critical speeds and to simulate the transient motion.The influence of gyroscopic moments on co-rotation and counter-rotation is analyzed,and the effect of the speed ratio on critical speed is studied.The dynamic characteristics under different conditions of increasing speed during start-up are demonstrated by comparison with transient nodal displacements.The presented model provides a complete foundation for further investigation of the dynamics of dual rotor systems.
基金supported by National Natural Science Foundation of China(No. 50175034).
文摘Single-point incremental forming （SPIF） is an innovational sheet metal forming method without dedicated dies, which belongs to rapid prototyping technology. In generalizing the SPIF of sheet metal, the deformation analysis on forming process becomes an important and useful method for the planning of shell products, the choice of material, the design of the forming process and the planning of the forming tool. Using solid brick elements, the finite element method（FEM） model of truncated pyramid was established. Based on the theory of anisotropy and assumed strain formulation, the SPIF processes with different parameters were simulated. The resulted comparison between the simulations and the experiments shows that the FEM model is feasible and effective. Then, according to the simulated forming process, the deformation pattern of SPIF can be summarized as the combination of plane-stretching deformation and bending deformation. And the study about the process parameters＇ impact on deformation shows that the process parameter of interlayer spacing is a dominant factor on the deformation. Decreasing interlayer spacing, the strain of one step decreases and the formability of blank will be improved. With bigger interlayer spacing, the plastic deformation zone increases and the forming force will be bigger.
基金financially supported by the National Natural Science Foundation of China(No.51571064)the National Basic Research Program of China(No.2014CB643701)。
文摘The numerical simulation model of single roller melt spinning for rapid quenching process of NdFeB alloy was built,and the vacuum chamber,cooling roller and sample were taken into account as a system.The existing mature technology was in order to verify the correctness of simulation.The rapid quenching ribbons with different roll speeds were used as the simulation objects.The results of the numerical simulation and experiments show that the validity of the model has been testified and the reasons of the formation of complete quenching ribbons and by-product have been explained.The experimental thickness of the ribbons is proportional to the theoretical thickness.In the same spray condition,with the roll speed increasing,the thickness decreases linearly.At the speed range of25-30 m·s^(-1),the simulated calculation date is close to the experimental date,which can be considered as an ideal technological parameter.
基金National Natural Science Foundation of China(No.51674280,51774308,51704033,51722406,51950410591)Shandong Provincial Natural Science Foundation(ZR2019JQ21,JQ201808)+3 种基金the Fundamental Research Funds for the Central Universities(No.20CX02113A)National Science and Technology Major Project(2016ZX05014-000407)Program for Changjiang Scholars and Innovative Research Team in University(IRT_16R69)PetroChina Innovation Foundation(No.2018D-5007-0210)。
文摘Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.
文摘Finite element method(FEM) was used to simulate the forming process of shotpeening the wing skin panel. Experiment of shotpeeing the wing skin panel was carried out. The results show that equivalent deformation in shotpeening process can be obtained using the elongation and bending result caused by thermal stress that is induced by applying temperature load on the surface of the part. Deformation of the part in the shotpeeing process can be analyzed using this method. The parameters and their relationships are identified.
基金supported by the Natural Sciences and Engineering Research Council of Canada through Discovery Grant 341275 (G. Grasselli) and Engage EGP 461019-13
文摘Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.
文摘Simulation of the temperature field of copier paper in copier fusing is very important for improving the fusing property of reprography. The temperature field of copier paper varies with a high gradient when the copier paper is moving through the fusing rollers. By means of conventional shaft elements, the high gradient temperature variety causes the oscillation of the numerical solution. Based on the Daubechies scaling functions, a kind of wavelet based element is constructed for the above problem. The temperature field of the copier paper moving through the fusing rollers is simulated using the two methods. Comparison of the results shows the advantages of the wavelet finite element method, which provides a new method for improving the copier properties.
基金Project supported by the National Natural Science Foundation of China(No.51078230)the Research Fund for the Doctoral Program of Higher Education of China(No.200802480056)the Key Project of Fund of Science and Technology Development of Shanghai(No.10JC1407900),China
文摘In this paper,a stabilized finite element technique,actualized by streamline upwind Petrov-Galerkin(SUPG) stabilized method and three-step finite element method(FEM),for large eddy simulation(LES) is developed to predict the wind flow with high Reynolds numbers.Weak form of LES motion equation is combined with the SUPG stabilized term for the spatial finite element discretization.An explicit three-step scheme is implemented for the temporal discretization.For the numerical example of 2D wind flow over a square rib at Re=4.2×105,the Smagorinsky's subgrid-scale(SSGS) model,the DSGS model,and the DSGS model with Cabot near-wall model are applied,and their results are analyzed and compared with experimental results.Furthermore,numerical examples of 3D wind flow around a surface-mounted cube with different Reynolds numbers are performed using DSGS model with Cabot near-wall model based on the present stabilized method to study the wind field and compared with experimental and numerical results.Finally,vortex structures for wind flow around a surface-mounted cube are studied by present numerical method.Stable and satisfactory results are obtained,which are consistent with most of the measurements even under coarse mesh.
基金supported by the National Nature Science Foundation of China (Grant No 90510017)
文摘This paper proposes a new, simple and efficient method for nonlinear simulation of arch dam cracking from the construction period to the operation period, which takes into account the arch dam construction process and temperature loads. In the calculation mesh, the contact surface of pair nodes is located at places on the arch dam where cracking is possible. A new effective iterative method, the mixed finite element method for friction-contact problems, is improved and used for nonlinear simulation of the cracking process. The forces acting on the structure are divided into two parts： external forces and contact forces. The displacement of the structure is chosen as the basic variable and the nodal contact force in the possible contact region of the local coordinate system is chosen as the iterative variable, so that the nonlinear iterative process is only limited within the possible contact surface and is much more economical. This method was used to simulate the cracking process of the Shuanghe Arch Dam in Southwest China. In order to prove the validity and accuracy of this method and to study the effect of thermal stress on arch dam cracking, three schemes were designed for calculation. Numerical results agree with actual measured data, proving that it is feasible to use this method to simulate the entire process of nonlinear arch dam cracking.
基金by Key Aviation Scientific and Technological Laboratory of High-speed Hydrodynamic under grant MJ-2015-F-028.
文摘A vein model was established to simulate the periodic characteristics of blood flow and valve deformation in blood-induced valve cycles.Using an immersed finite element method which was modified by a ghost fluid technique,the interaction between the vein and blood was simulated.With an independent solid solver,the contact force between vein tissues was calculated using an adhesive contact method.A benchmark simulation of the normal valve cycle validated the proposed model for a healthy vein.Both the opening orifice and blood flow rate agreed with those in the physiology.Low blood shear stress and maximum leaflet stress were also seen in the base region of the valve.On the basis of the healthy model,a diseased vein model was subsequently built to explore the sinus lesions,namely,fibrosis and atrophy which are assumed stiffening and softening of the sinus.Our results showed the opening orifice of the diseased vein was inversely proportional to the corresponding modulus of the sinus.A drop in the transvalvular pressure gradient resulted from the sinus lesion.Compared to the fibrosis,the atrophy of the sinus apparently improved the vein deformability but simultaneously accelerated the deterioration of venous disease and increased the risk of potential fracture.These results provide understandings of the normal/abnormal valve cycle in vein,and can be also helpful for the prosthesis design.
文摘Different material properties leads to different metal fracture behaviors. Even if the powder material is composed of plastic metal, the fracture still does not show macroscopic plastic deformation characteristics if the material contains a large number of voids. Eight node isoparametric elastic plastic finite element method was used to simulate the tensile process of sintered powder material. By setting a number of voids in the analyzed metal cuboid, the initial density was taken into consideration. The material properties of the three dimensional solid for the tensile simulation were defined with reference to the known pure iron material parameters. The load displacement curves during elongation were obtained with a universal testing machine, and then the simulated curves were compared with the experimental results. The factors that cause the stress concentration and strength decrease were analyzed according to the simulated equivalent von Mises stress distribution.
基金supported by the National Natural Science Foundation of China(Nos.51378293,51078199,50678093,and 50278046)the Program for Changjiang Scholars and the Innovative Research Team in University of China(No.IRT00736)
文摘The element energy projection （EEP） method for computation of super- convergent resulting in a one-dimensional finite element method （FEM） is successfully used to self-adaptive FEM analysis of various linear problems, based on which this paper presents a substantial extension of the whole set of technology to nonlinear problems. The main idea behind the technology transfer from linear analysis to nonlinear analysis is to use Newton＇s method to linearize nonlinear problems into a series of linear problems so that the EEP formulation and the corresponding adaptive strategy can be directly used without the need for specific super-convergence formulation for nonlinear FEM. As a re- sult, a unified and general self-adaptive algorithm for nonlinear FEM analysis is formed. The proposed algorithm is found to be able to produce satisfactory finite element results with accuracy satisfying the user-preset error tolerances by maximum norm anywhere on the mesh. Taking the nonlinear ordinary differential equation （ODE） of second-order as the model problem, this paper describes the related fundamental idea, the imple- mentation strategy, and the computational algorithm. Representative numerical exam- ples are given to show the efficiency, stability, versatility, and reliability of the proposed approach.
