In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results i...In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.展开更多
This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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 predic...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.展开更多
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.展开更多
The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at th...The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.展开更多
The basic principle of corrode groove on outside of steel pipe during storage was analyzed in this paper, namely the water film on the contacted surface of steel pipe, which gathered from humidity in the air, rain or ...The basic principle of corrode groove on outside of steel pipe during storage was analyzed in this paper, namely the water film on the contacted surface of steel pipe, which gathered from humidity in the air, rain or gel, and the suspended particles in air, and the corrosive composition, such as SO2, CO2, O2 and NaCI, in addition to the inhomogeneity of the organization and composition, which lead to the corrosion cell reaction, so that cause the corrosion initial from the contact surface of the between steel pipes, so as to form the corrosion groove. At the same time, the corrosion groove with depth of 0.125t （t pipe wall thickness） on the pipe of φ 1016 mm×21 mm ×70 API SPEC 5L was simulated using the FEM （finite element method）, and the stress and strain distribution of the defect area near corrosion groove were solved at the inner pressure of 12 MPa, 10 MPa, 8 MPa, 6 MPa, 4 MPa and 2 MPa, respectively, which showed that no matter the pressure values were, the maximum stress and strain were lied at the bottom of corrosion defects groove and were in good linear relationship with the internal pressure increasing from 2 MPa to 6 MPa. When the internal pres- sures were greater than 6 MPa, they felled into the nonlinear model and to be yielded or even to be destroyed. In addition, the residual strength and the limit operation pressure of the corrode pipe with the defects groove of 0.125t were calculated or simulated according to the theoretical calculation, the finite element method based on the stress, the finite element method based on strain, DNV-RP-F101, ASME B31G and experimental methods respectively. The results showed that the residual strength and the limit operation pressure of the defective parts solved by the finite element method based on stress were 424 MPa, and 15.34 MPa, respectively, which was very close to that of experimental method, the residual strength was 410 MPa and the limit operation pressure 14.78 MPa. Besides, the results also showed that it was feasible and effective to simulate the residual strength of the structure with corrosion defects using the finite element method.展开更多
文摘In modern ZnO varistors,traditional aging mechanisms based on increased power consumption are no longer relevant due to reduced power consumption during DC aging.Prolonged exposure to both AC and DC voltages results in increased leakage current,decreased breakdown voltage,and lower nonlinearity,ultimately compromising their protective performance.To investigate the evolution in electrical properties during DC aging,this work developed a finite element model based on Voronoi networks and conducted accelerated aging tests on commercial varistors.Throughout the aging process,current-voltage characteristics and Schottky barrier parameters were measured and analyzed.The results indicate that when subjected to constant voltage,current flows through regions with larger grain sizes,forming discharge channels.As aging progresses,the current focus increases on these channels,leading to a decline in the varistor’s overall performance.Furthermore,analysis of the Schottky barrier parameters shows that the changes in electrical performance during aging are non-monotonic.These findings offer theoretical support for understanding the aging mechanisms and condition assessment of modern stable ZnO varistors.
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
文摘This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
基金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.
文摘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.
基金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.
基金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.
文摘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.
基金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.
基金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.
基金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.
文摘The finite element method （FEM） and particle image velocimetry （PIV） technique are utilized to get the flow field along the inlet passage, the chamber, the metering port and the outlet passage of spool valve at three different valve openings. For FEM numerical simulation, the stream function ψ-vorticity ω forms of continuity and Navier-Stokes equations are employed and FEM is applied to discrete the equations. Homemade simulation codes are executed to compute the values of stream function and vorticity at each node in the flow domain, then according to the correlation between stream function and velocity components, the velocity vectors of the whole field are calculated. For PIV experiment, pulse Nd： YAG laser is exploited to generate laser beam, cylindrical and spherical lenses are combined each other to produce 1.0 mm thickness laser sheet to illuminate the object plane, Polystyrene spherical particle with diameter of 30-50 μm is seeded in the fluid as a tracing particles, Kodak ES 1.0 CCD camera is employed to capture the images of interested, the images are processed with fast Fourier transform （FFT） cross-correlation algorithm and the processing results is displayed. Both results of numerical simulation and PIV experimental show that there are three main areas in the spool valve where vortex is formed. Numerical results also indicate that the valve opening have some effects on the flow structure of the valve. The investigation is helpful for qualitatively analyzing the energy loss, noise generating, steady state flow forces and even designing the geometry structure and flow passage.
基金supported by the National Natural Science Foundation of China(Nos.51101127 and 51171154)Soar Star of Northwestern Polytechnical University(2011)Fundamental Research Foundation of Northwestern Polytechnical University(No.JC201213)
文摘The basic principle of corrode groove on outside of steel pipe during storage was analyzed in this paper, namely the water film on the contacted surface of steel pipe, which gathered from humidity in the air, rain or gel, and the suspended particles in air, and the corrosive composition, such as SO2, CO2, O2 and NaCI, in addition to the inhomogeneity of the organization and composition, which lead to the corrosion cell reaction, so that cause the corrosion initial from the contact surface of the between steel pipes, so as to form the corrosion groove. At the same time, the corrosion groove with depth of 0.125t （t pipe wall thickness） on the pipe of φ 1016 mm×21 mm ×70 API SPEC 5L was simulated using the FEM （finite element method）, and the stress and strain distribution of the defect area near corrosion groove were solved at the inner pressure of 12 MPa, 10 MPa, 8 MPa, 6 MPa, 4 MPa and 2 MPa, respectively, which showed that no matter the pressure values were, the maximum stress and strain were lied at the bottom of corrosion defects groove and were in good linear relationship with the internal pressure increasing from 2 MPa to 6 MPa. When the internal pres- sures were greater than 6 MPa, they felled into the nonlinear model and to be yielded or even to be destroyed. In addition, the residual strength and the limit operation pressure of the corrode pipe with the defects groove of 0.125t were calculated or simulated according to the theoretical calculation, the finite element method based on the stress, the finite element method based on strain, DNV-RP-F101, ASME B31G and experimental methods respectively. The results showed that the residual strength and the limit operation pressure of the defective parts solved by the finite element method based on stress were 424 MPa, and 15.34 MPa, respectively, which was very close to that of experimental method, the residual strength was 410 MPa and the limit operation pressure 14.78 MPa. Besides, the results also showed that it was feasible and effective to simulate the residual strength of the structure with corrosion defects using the finite element method.
