A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale met...A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.展开更多
Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant i...Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.展开更多
This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material ...This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material for the specimens consists of a single hole model(25% void/cell,16% void/cell and 10% void/cell)and a four-hole model(25%void/cell).Using a representative volume element(RVE),we try to produce the equivalent homogenized properties and work on a homogeneous specimen for the study of fretting fatigue.Next,the fretting fatigue contact problem is performed for 3 new cases of models that consist of a homogeneous and a heterogeneous part(single hole cell)in the contact area.The aim is to analyze the normal and shear stresses of these models and compare them with the results of the corresponding heterogeneous models based on the Direct Numerical Simulation(DNS)method.Finally,by comparing the computational time and%deviations,we draw conclusions about the reliability and effectiveness of the proposed method.展开更多
In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-（ 1/4/λ）,meso-（ 1/100λ and 1/4λ） and micro-（ 1/100λ） scales.Based on the multi-scales rock physics ...In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-（ 1/4/λ）,meso-（ 1/100λ and 1/4λ） and micro-（ 1/100λ） scales.Based on the multi-scales rock physics modelling technique,the seismic azimuthal anisotropy characteristic is analyzed for distinguishing the fractures of meso-scale.Furthermore,by integrating geological core fracture description,image well-logging fracture interpretation,seismic attributes macro-scale fracture prediction and core slice micro-scale fracture characterization,an comprehensive multi-scale fracture prediction methodology and technique workflow are proposed by using geology,well-logging and seismic multi-attributes.Firstly,utilizing the geology core slice observation（Fractures description） and image well-logging data interpretation results,the main governing factors of fracture development are obtained,and then the control factors of the development of regional macro-scale fractures are carried out via modelling of the tectonic stress field.For the meso-scale fracture description,the poststack geometric attributes are used to describe the macro-scale fracture as well,the prestack attenuation seismic attribute is used to predict the meso-scale fracture.Finally,by combining lithological statistic inversion with superposed results of faults,the relationship of the meso-scale fractures,lithology and faults can be reasonably interpreted and the cause of meso-scale fractures can be verified.The micro-scale fracture description is mainly implemented by using the electron microscope scanning of cores.Therefore,the development of fractures in reservoirs is assessed by valuating three classes of fracture prediction results.An integrated fracture prediction application to a real field in Sichuan basin,where limestone reservoir fractures developed,is implemented.The application results in the study area indicates that the proposed multi-scales integrated fracture prediction method and the technique procedureare able to deal with the strong heterogeneity and multi-scales problems in fracture prediction.Moreover,the multi-scale fracture prediction technique integrated with geology,well-logging and seismic multi-information can help improve the reservoir characterization and sweet-spots prediction for the fractured hydrocarbon reservoirs.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multi...An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
The support of catalyst for the direct synthesis of diphenyl carbonate （DPC） by heterogeneous catalytic reaction was prepared by the sol-gel method. Compared with activated charcoal, molecular sieve, porous ceramics...The support of catalyst for the direct synthesis of diphenyl carbonate （DPC） by heterogeneous catalytic reaction was prepared by the sol-gel method. Compared with activated charcoal, molecular sieve, porous ceramics, hopcalite, the support prepared by the sol-gel method has higher activity. The characterization of the support by X-ray diffraction （XRD） and transmission electron microscope （TEM） show that the mare crystal phase is Co2MnO4 and the average particle diameter is about 40 nm. The optimum conditions for synthesis of the support were determined by orthogonal experiments, which indicate that the proportion of Cu, Mn, and Co is the first important factor influencing the yield and selectivity of DPC. Temperature of calcination is the second one. The optimum conditions are： molar proportion of Cu, Mn, and Co being 1 ： 1 ： 1, temperature of calcination 700℃, drying at 100~C, temperature of water bath 85~C. The yield and selectivity of DPC in the process can reach 38% and 99% in the batch operation, respectively. The copper cobalt manganese mixed oxides chosen as the support contribute more to the high catalytic activity than the sol-gel method.展开更多
A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the clas...A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the classical multi-scale method(MSM)and the method of variation of parameters.Assuming that the forced excitation is a constant,one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM.Then,this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation.To certify the correctness and precision of the proposed analytical method,the van der Pol system with two scales subject to slowly periodic excitation is investigated;this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS),and quiescence(QS)responses.The approximate analytical expressions of the three types of responses are given by the MMSM,and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM).The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively,particularly when the excitation frequency is far less than the natural frequency of the system.展开更多
A new method, the dynamic reduction method(DRM) combined with the strain-softening method, was applied to evaluate the possible slip surface of a highly heterogeneous rock slope of the Dagangshan hydropower station in...A new method, the dynamic reduction method(DRM) combined with the strain-softening method, was applied to evaluate the possible slip surface of a highly heterogeneous rock slope of the Dagangshan hydropower station in Southwest China.In DRM, only the strength of the failure elements is reduced and the softening reduction factor K is adopted to calculate the strength parameters. The simulation results calculated by DRM show that the further slip surface on the right slope of the Dagangshan hydropower station is limited in the middle part of the slope, while both SRM(strength reduction method) and LEM(limit equilibrium method) predict a failure surface which extends upper and longer. The observations and analysis from the three recorded sliding events indicate that the failure mode predicted by DRM is more likely the scenario.The results in this study illustrate that for highly heterogeneous slopes with geological discontinuities in different length scales, the proposed DRM can provide a reliable prediction of the location of the slip surface.展开更多
According to the requirement of heterogeneous object modeling in additive manufacturing(AM),the Non-Uniform Rational B-Spline(NURBS)method has been applied to the digital representation of heterogeneous object in this...According to the requirement of heterogeneous object modeling in additive manufacturing(AM),the Non-Uniform Rational B-Spline(NURBS)method has been applied to the digital representation of heterogeneous object in this paper.By putting forward the NURBS material data structure and establishing heterogeneous NURBS object model,the accurate mathematical unified representation of analytical and free heterogeneous objects have been realized.With the inverse modeling of heterogeneous NURBS objects,the geometry and material distribution can be better designed to meet the actual needs.Radical Basis Function(RBF)method based on global surface reconstruction and the tensor product surface interpolation method are combined to RBF-NURBS inverse construction method.The geometric and/or material information of regular mesh points is obtained by RBF interpolation of scattered data,and the heterogeneous NURBS surface or object model is obtained by tensor product interpolation.The examples have shown that the heterogeneous objects fitting to scattered data points can be generated effectively by the inverse construction methods in this paper and 3D CAD models for additive manufacturing can be provided.展开更多
Tight glutenite reservoirs are known for strong heterogeneity,complex wettability,and challenging development.Gas-Assisted Gravity Drainage(GAGD)technology has the potential to significantly improve recovery efficienc...Tight glutenite reservoirs are known for strong heterogeneity,complex wettability,and challenging development.Gas-Assisted Gravity Drainage(GAGD)technology has the potential to significantly improve recovery efficiency in glutenite reservoir.However,there is currently limited research on GAGD processes specifically designed for glutenite reservoirs,and there is a lack of relevant dimensionless numbers for predicting recovery efficiency.In this study,we developed a theoretical model based on the characteristics of glutenite reservoirs and used phase-field method to track the oil-gas interface for numerical simulations of dynamic GAGD processes.To explore the factors influencing gas-driven recovery,we simulated the effects of strong heterogeneity and dynamic wettability on the construction process under gravity assistance.Additionally,we introduced multiple dimensionless numbers(including capillary number,viscosity ratio,and Bond number)and conducted a series of numerical simulations.The results demonstrate that gravity enhances the stability of the oil-gas interface but causes unstable pressure fluctuations when passing through different-sized throat regions,particularly leading to front advancement in smaller throats.Although strong heterogeneity has negative impacts on GAGD,they can be mitigated by reducing injection velocity.Increasing oil-wettability promotes oil displacement by overcoming capillary forces,particularly in narrower pores,allowing residual oils to be expelled.Among the dimensionless numbers,the recovery efficiency is directly proportional to the Bond number and inversely proportional to the capillary number and viscosity ratio.Through sensitivity analysis of the dimensionless numbers’impact on the recovery efficiency,a new dimensionless N_(Glu) considering heterogeneity is proposed to accurately predict GAGD recovery of tight glutenite reservoirs.展开更多
Self-vibrating systems comprised of active materials have great potential for application in the fields of energy harvesting,actuation,bionic instrumentation,and autonomous robotics.However,it is challenging to obtain...Self-vibrating systems comprised of active materials have great potential for application in the fields of energy harvesting,actuation,bionic instrumentation,and autonomous robotics.However,it is challenging to obtain analytical solutions describing these systems,which hinders analysis and design.In this work,we propose a self-vibrating liquid crystal elastomer(LCE)fiber-spring system exposed to spatially-constant gradient light,and determine analytical solutions for its amplitude and period.First,using a dynamic model of LCE,we obtain the equations governing the self-vibration.Then,we analyze two different motion states and elucidate the mechanism of self-vibration.Subsequently,we derive analytical solutions for the amplitude and frequency using the multi-scale method,and compare the solutions with numerical results.The analytical outcomes are shown to be consistent with the numerical calculations,while taking far less computational time.Our findings reveal the utility of the multi-scale method in describing self-vibration,which may contribute to more efficient and accurate analyses of self-vibrating systems.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
In this paper, we present a new algorithm to solve a kind of nonlinear time space-fractional partial differential equations on a finite domain. The method is based on B-spline wavelets approximations, some of these fu...In this paper, we present a new algorithm to solve a kind of nonlinear time space-fractional partial differential equations on a finite domain. The method is based on B-spline wavelets approximations, some of these functions are reshaped to satisfy on boundary conditions exactly. The Adams fractional method is used to reduce the problem to a system of equations. By multiscale method this system is divided into some smaller systems which have less computations. We get an approximated solution which is more accurate on some subdomains by combining the solutions of these systems. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique, also the stability of the method is discussed.展开更多
A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D F...A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.展开更多
This work presents a novel approach to achieve nonlinear vibration response based on the Hamilton principle.We chose the 5-MW reference wind turbine which was established by the National Renewable Energy Laboratory(NR...This work presents a novel approach to achieve nonlinear vibration response based on the Hamilton principle.We chose the 5-MW reference wind turbine which was established by the National Renewable Energy Laboratory(NREL),to research the effects of the nonlinear flap-wise vibration characteristics.The turbine wheel is simplified by treating the blade of a wind turbine as an Euler-Bernoulli beam,and the nonlinear flap-wise vibration characteristics of the wind turbine blades are discussed based on the simplification first.Then,the blade’s large-deflection flap-wise vibration governing equation is established by considering the nonlinear term involving the centrifugal force.Lastly,it is truncated by the Galerkin method and analyzed semi-analytically using the multi-scale analysis method,and numerical simulations are carried out to compare the simulation results of finite elements with the numerical simulation results using Campbell diagram analysis of blade vibration.The results indicated that the rotational speed of the impeller has a significant impact on blade vibration.When the wheel speed of 12.1 rpm and excitation amplitude of 1.23 the maximum displacement amplitude of the blade has increased from 0.72 to 3.16.From the amplitude-frequency curve,it can be seen that the multi-peak characteristic of blade amplitude frequency is under centrifugal nonlinearity.Closed phase trajectories in blade nonlinear vibration,exhibiting periodic motion characteristics,are found through phase diagrams and Poincare section diagrams.展开更多
In order to study the mechanical properties of the heterogeneous core plate of the wind turbine blade,a modeling method of the core plate based on displacement field variables is proposed.Firstly,the wind turbine blad...In order to study the mechanical properties of the heterogeneous core plate of the wind turbine blade,a modeling method of the core plate based on displacement field variables is proposed.Firstly,the wind turbine blade core plate was modeled according to the theory of modeling heterogeneous material characteristics.Secondly,the three-point bending finite element model of the wind turbine blade core plate was solved by the display dynamic equation to obtain the deformation pattern and force-deformation relationship of the core plate.Finally,the three-point bending static test was conducted to compare with the finite element analysis.The test results show that:the damage form of the wind turbine blade core plate includes elasticity,yield,and failure stages.The main failure modes are plastic deformation,core material collapse,and panel-core delamination.The failure load measured by the test is 1.59 kN,which is basically consistent with the load-displacement result obtained by the simulation,with a difference of only 1.9%,which verifies the validity and reliability of the model.It provides data references for wind turbine blade structure design.展开更多
The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and faul...The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and fault are developed,so its development cannot be conducted just based on static and dynamic reservoir description.In order to solve this problem,this paper establishes a numerical well test model of vertical wells based on matrix,fractures and faults(large fractures and small faults)by combining the random generation of natural fracture networks with the unstructured discrete fracture modeling method to break through the traditional continuous medium well test model.In addition,the model is solved by using the finite element method with mixed element,and the typical well test type curves under different random fracture networks are obtained.And the following research results are obtained.First,based on the observed data,the fracture network distribution modes of fractured tight sandstone gas reservoirs are classified into three categories.The influence of random generation of fracture networks on typical well test type curves is discussed.The results of discrete fracture well test model are compared with those of the traditional continuous medium well test model,and the applicable conditions of the traditional continuous medium well test model is determined.Second,there are great differences between the results of discrete fracture model and those of dual porosity medium model.1 The dual porosity medium model is a special case of the discrete fracture model,in which the fractures are evenly distributed within infinitely small spacing.Third,the characteristics of well test type curves under three fracture network distribution modes are discussed.The well test type curves that cannot be interpreted by the conventional dual/triple porosity continuous medium model are successfully interpreted by using the established well test interpretation model of random discrete fracture.The curve matching effect is ideal and the interpreted parameters are reasonable.In conclusion,the new model and the new method reveal the development mechanism of step-by-step production and coordinated gas supply between media of different scales,explain the development characteristics of large inter-well productivity difference and abnormal rapid inter-well pressure response,and provide a reference for the development of similar gas reservoirs.展开更多
基金Supported by the National Natural Science Foundation of China(51105195,51075204)the Aeronautical Science Foundation of China(2011ZB52024)
文摘A straightforward multi-scale boundary element method is proposed for global and local mechanical analysis of heterogeneous material.The method is more accurate and convenient than finite element based multi-scale method.The formulations of this method are derived by combining the homogenization approach and the fundamental equations of boundary element method.The solution gives the convenient formulations to compute global elastic constants and the local stress field.Finally,two numerical examples of porous material are presented to prove the accuracy and the efficiency of the proposed method.The results show that the method does not require the iteration to obtain the solution of the displacement in micro level.
基金supported by the National Natural Science Foundation of China(Grant Nos.42102346,42172301).
文摘Magneto-electro-elastic(MEE)materials are widely utilized across various fields due to their multi-field coupling effects.Consequently,investigating the coupling behavior of MEE composite materials is of significant importance.The traditional finite element method(FEM)remains one of the primary approaches for addressing such issues.However,the application of FEM typically necessitates the use of a fine finite element mesh to accurately capture the heterogeneous properties of the materials and meet the required computational precision,which inevitably leads to a reduction in computational efficiency.To enhance the computational accuracy and efficiency of the FEM for heterogeneous multi-field coupling problems,this study presents the coupling magneto-electro-elastic multiscale finite element method(CM-MsFEM)for heterogeneous MEE structures.Unlike the conventional multiscale FEM(MsFEM),the proposed algorithm simultaneously constructs displacement,electric,and magnetic potential multiscale basis functions to address the heterogeneity of the corresponding parameters.The macroscale formulation of CM-MsFEM was derived,and the macroscale/microscale responses of the problems were obtained through up/downscaling calculations.Evaluation using numerical examples analyzing the transient behavior of heterogeneous MEE structures demonstrated that the proposed method outperforms traditional FEM in terms of both accuracy and computational efficiency,making it an appropriate choice for numerically modeling the dynamics of heterogeneous MEE structures.
文摘This paper deals with modeling of the phenomenon of fretting fatigue in heterogeneous materials using the multi-scale computational homogenization technique and finite element analysis(FEA).The heterogeneous material for the specimens consists of a single hole model(25% void/cell,16% void/cell and 10% void/cell)and a four-hole model(25%void/cell).Using a representative volume element(RVE),we try to produce the equivalent homogenized properties and work on a homogeneous specimen for the study of fretting fatigue.Next,the fretting fatigue contact problem is performed for 3 new cases of models that consist of a homogeneous and a heterogeneous part(single hole cell)in the contact area.The aim is to analyze the normal and shear stresses of these models and compare them with the results of the corresponding heterogeneous models based on the Direct Numerical Simulation(DNS)method.Finally,by comparing the computational time and%deviations,we draw conclusions about the reliability and effectiveness of the proposed method.
基金supported by the national oil and gas major project(No.2011ZX05019-008)National Natural Science Foundation of China(No.41574108 and U1262208)presented at the Exploration Geophysics Symposium 2015 of the EAGE Local Chapter China
文摘In this paper,we implement three scales of fracture integrated prediction study by classifying it to macro-（ 1/4/λ）,meso-（ 1/100λ and 1/4λ） and micro-（ 1/100λ） scales.Based on the multi-scales rock physics modelling technique,the seismic azimuthal anisotropy characteristic is analyzed for distinguishing the fractures of meso-scale.Furthermore,by integrating geological core fracture description,image well-logging fracture interpretation,seismic attributes macro-scale fracture prediction and core slice micro-scale fracture characterization,an comprehensive multi-scale fracture prediction methodology and technique workflow are proposed by using geology,well-logging and seismic multi-attributes.Firstly,utilizing the geology core slice observation（Fractures description） and image well-logging data interpretation results,the main governing factors of fracture development are obtained,and then the control factors of the development of regional macro-scale fractures are carried out via modelling of the tectonic stress field.For the meso-scale fracture description,the poststack geometric attributes are used to describe the macro-scale fracture as well,the prestack attenuation seismic attribute is used to predict the meso-scale fracture.Finally,by combining lithological statistic inversion with superposed results of faults,the relationship of the meso-scale fractures,lithology and faults can be reasonably interpreted and the cause of meso-scale fractures can be verified.The micro-scale fracture description is mainly implemented by using the electron microscope scanning of cores.Therefore,the development of fractures in reservoirs is assessed by valuating three classes of fracture prediction results.An integrated fracture prediction application to a real field in Sichuan basin,where limestone reservoir fractures developed,is implemented.The application results in the study area indicates that the proposed multi-scales integrated fracture prediction method and the technique procedureare able to deal with the strong heterogeneity and multi-scales problems in fracture prediction.Moreover,the multi-scale fracture prediction technique integrated with geology,well-logging and seismic multi-information can help improve the reservoir characterization and sweet-spots prediction for the fractured hydrocarbon reservoirs.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
基金supported by the National Natural Science Foundation(10721062,11072051,90715037,10728205,91015003, 51021140004)the Program of Introducing Talents of Discipline to Universities(B08014)the National Key Basic Research Special Foundation of China(2010CB832704).
文摘An extended multiscale finite element method （EMsFEM） is developed for solving the mechanical problems of heterogeneous materials in elasticity.The underlying idea of the method is to construct numerically the multiscale base functions to capture the small-scale features of the coarse elements in the multiscale finite element analysis.On the basis of our existing work for periodic truss materials, the construction methods of the base functions for continuum heterogeneous materials are systematically introduced. Numerical experiments show that the choice of boundary conditions for the construction of the base functions has a big influence on the accuracy of the multiscale solutions, thus,different kinds of boundary conditions are proposed. The efficiency and accuracy of the developed method are validated and the results with different boundary conditions are verified through extensive numerical examples with both periodic and random heterogeneous micro-structures.Also, a consistency test of the method is performed numerically. The results show that the EMsFEM can effectively obtain the macro response of the heterogeneous structures as well as the response in micro-scale,especially under the periodic boundary conditions.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
文摘The support of catalyst for the direct synthesis of diphenyl carbonate （DPC） by heterogeneous catalytic reaction was prepared by the sol-gel method. Compared with activated charcoal, molecular sieve, porous ceramics, hopcalite, the support prepared by the sol-gel method has higher activity. The characterization of the support by X-ray diffraction （XRD） and transmission electron microscope （TEM） show that the mare crystal phase is Co2MnO4 and the average particle diameter is about 40 nm. The optimum conditions for synthesis of the support were determined by orthogonal experiments, which indicate that the proportion of Cu, Mn, and Co is the first important factor influencing the yield and selectivity of DPC. Temperature of calcination is the second one. The optimum conditions are： molar proportion of Cu, Mn, and Co being 1 ： 1 ： 1, temperature of calcination 700℃, drying at 100~C, temperature of water bath 85~C. The yield and selectivity of DPC in the process can reach 38% and 99% in the batch operation, respectively. The copper cobalt manganese mixed oxides chosen as the support contribute more to the high catalytic activity than the sol-gel method.
基金the National Natural Science Foundation of China(Nos.11672191,11772206,and U1934201)the Hundred Excellent Innovative Talents Support Program in Hebei University(No.SLRC2017053)。
文摘A simple,yet accurate modi?ed multi-scale method(MMSM)for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed.This method depends on the classical multi-scale method(MSM)and the method of variation of parameters.Assuming that the forced excitation is a constant,one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM.Then,this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation.To certify the correctness and precision of the proposed analytical method,the van der Pol system with two scales subject to slowly periodic excitation is investigated;this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS),and quiescence(QS)responses.The approximate analytical expressions of the three types of responses are given by the MMSM,and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM).The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively,particularly when the excitation frequency is far less than the natural frequency of the system.
基金supported by the National Key R&D Program of China (2017YFC1501301)the National Natural Science Foundation of China (Grant Nos. 41521002, 41572283 and 41130745)+2 种基金the Funding of Science and Technology Office of Sichuan Province (Grant Nos. 2015JQ0020 and 2017TD0018)the 1000 Young Talent Program of Chinathe research fund of the State Key Laboratory of Geohazard Prevention and Geoenvironment Protection (Grant No. SKLGP2017Z012)
文摘A new method, the dynamic reduction method(DRM) combined with the strain-softening method, was applied to evaluate the possible slip surface of a highly heterogeneous rock slope of the Dagangshan hydropower station in Southwest China.In DRM, only the strength of the failure elements is reduced and the softening reduction factor K is adopted to calculate the strength parameters. The simulation results calculated by DRM show that the further slip surface on the right slope of the Dagangshan hydropower station is limited in the middle part of the slope, while both SRM(strength reduction method) and LEM(limit equilibrium method) predict a failure surface which extends upper and longer. The observations and analysis from the three recorded sliding events indicate that the failure mode predicted by DRM is more likely the scenario.The results in this study illustrate that for highly heterogeneous slopes with geological discontinuities in different length scales, the proposed DRM can provide a reliable prediction of the location of the slip surface.
文摘According to the requirement of heterogeneous object modeling in additive manufacturing(AM),the Non-Uniform Rational B-Spline(NURBS)method has been applied to the digital representation of heterogeneous object in this paper.By putting forward the NURBS material data structure and establishing heterogeneous NURBS object model,the accurate mathematical unified representation of analytical and free heterogeneous objects have been realized.With the inverse modeling of heterogeneous NURBS objects,the geometry and material distribution can be better designed to meet the actual needs.Radical Basis Function(RBF)method based on global surface reconstruction and the tensor product surface interpolation method are combined to RBF-NURBS inverse construction method.The geometric and/or material information of regular mesh points is obtained by RBF interpolation of scattered data,and the heterogeneous NURBS surface or object model is obtained by tensor product interpolation.The examples have shown that the heterogeneous objects fitting to scattered data points can be generated effectively by the inverse construction methods in this paper and 3D CAD models for additive manufacturing can be provided.
基金supported by the National Natural Science Foundation of China(U22B2075)the Fundamental Research Funds for the Central Universities(2024ZKPYSB03)support from Beijing University of Science and Technology.
文摘Tight glutenite reservoirs are known for strong heterogeneity,complex wettability,and challenging development.Gas-Assisted Gravity Drainage(GAGD)technology has the potential to significantly improve recovery efficiency in glutenite reservoir.However,there is currently limited research on GAGD processes specifically designed for glutenite reservoirs,and there is a lack of relevant dimensionless numbers for predicting recovery efficiency.In this study,we developed a theoretical model based on the characteristics of glutenite reservoirs and used phase-field method to track the oil-gas interface for numerical simulations of dynamic GAGD processes.To explore the factors influencing gas-driven recovery,we simulated the effects of strong heterogeneity and dynamic wettability on the construction process under gravity assistance.Additionally,we introduced multiple dimensionless numbers(including capillary number,viscosity ratio,and Bond number)and conducted a series of numerical simulations.The results demonstrate that gravity enhances the stability of the oil-gas interface but causes unstable pressure fluctuations when passing through different-sized throat regions,particularly leading to front advancement in smaller throats.Although strong heterogeneity has negative impacts on GAGD,they can be mitigated by reducing injection velocity.Increasing oil-wettability promotes oil displacement by overcoming capillary forces,particularly in narrower pores,allowing residual oils to be expelled.Among the dimensionless numbers,the recovery efficiency is directly proportional to the Bond number and inversely proportional to the capillary number and viscosity ratio.Through sensitivity analysis of the dimensionless numbers’impact on the recovery efficiency,a new dimensionless N_(Glu) considering heterogeneity is proposed to accurately predict GAGD recovery of tight glutenite reservoirs.
基金supported by the National Natural Science Foundation of China(No.12172001)the University Natural Science Research Project of Anhui Province(No.2022AH020029)+1 种基金the Anhui Provincial Natural Science Foundation(Nos.2208085Y01 and 2008085QA23)the Housing and Urban-Rural Development Science and Technology Project of Anhui Province(No.2023-YF129),China.
文摘Self-vibrating systems comprised of active materials have great potential for application in the fields of energy harvesting,actuation,bionic instrumentation,and autonomous robotics.However,it is challenging to obtain analytical solutions describing these systems,which hinders analysis and design.In this work,we propose a self-vibrating liquid crystal elastomer(LCE)fiber-spring system exposed to spatially-constant gradient light,and determine analytical solutions for its amplitude and period.First,using a dynamic model of LCE,we obtain the equations governing the self-vibration.Then,we analyze two different motion states and elucidate the mechanism of self-vibration.Subsequently,we derive analytical solutions for the amplitude and frequency using the multi-scale method,and compare the solutions with numerical results.The analytical outcomes are shown to be consistent with the numerical calculations,while taking far less computational time.Our findings reveal the utility of the multi-scale method in describing self-vibration,which may contribute to more efficient and accurate analyses of self-vibrating systems.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
文摘In this paper, we present a new algorithm to solve a kind of nonlinear time space-fractional partial differential equations on a finite domain. The method is based on B-spline wavelets approximations, some of these functions are reshaped to satisfy on boundary conditions exactly. The Adams fractional method is used to reduce the problem to a system of equations. By multiscale method this system is divided into some smaller systems which have less computations. We get an approximated solution which is more accurate on some subdomains by combining the solutions of these systems. Illustrative examples are included to demonstrate the validity and applicability of our proposed technique, also the stability of the method is discussed.
基金supported by the National Natural Science Foundation of China (51109029,51178081,51138001,and 51009020)the State Key Development Program for Basic Research of China (2013CB035905)
文摘A new finite element method （FEM） of B-spline wavelet on the interval （BSWI） is proposed. Through analyzing the scaling functions of BSWI in one dimension, the basic formula for 2D FEM of BSWI is deduced. The 2D FEM of 7 nodes and 10 nodes are constructed based on the basic formula. Using these proposed elements, the multiscale numerical model for foundation subjected to harmonic periodic load, the foundation model excited by external and internal dynamic load are studied. The results show the pro- posed finite elements have higher precision than the tradi- tional elements with 4 nodes. The proposed finite elements can describe the propagation of stress waves well whenever the foundation model excited by extemal or intemal dynamic load. The proposed finite elements can be also used to con- nect the multi-scale elements. And the proposed finite elements also have high precision to make multi-scale analysis for structure.
基金supported by the National Natural Science Foundation of China(No.51965034).
文摘This work presents a novel approach to achieve nonlinear vibration response based on the Hamilton principle.We chose the 5-MW reference wind turbine which was established by the National Renewable Energy Laboratory(NREL),to research the effects of the nonlinear flap-wise vibration characteristics.The turbine wheel is simplified by treating the blade of a wind turbine as an Euler-Bernoulli beam,and the nonlinear flap-wise vibration characteristics of the wind turbine blades are discussed based on the simplification first.Then,the blade’s large-deflection flap-wise vibration governing equation is established by considering the nonlinear term involving the centrifugal force.Lastly,it is truncated by the Galerkin method and analyzed semi-analytically using the multi-scale analysis method,and numerical simulations are carried out to compare the simulation results of finite elements with the numerical simulation results using Campbell diagram analysis of blade vibration.The results indicated that the rotational speed of the impeller has a significant impact on blade vibration.When the wheel speed of 12.1 rpm and excitation amplitude of 1.23 the maximum displacement amplitude of the blade has increased from 0.72 to 3.16.From the amplitude-frequency curve,it can be seen that the multi-peak characteristic of blade amplitude frequency is under centrifugal nonlinearity.Closed phase trajectories in blade nonlinear vibration,exhibiting periodic motion characteristics,are found through phase diagrams and Poincare section diagrams.
基金funded by National Natural Science Foundation of China(Grant No.52075305)Natural Science Foundation of Shandong Province(Grant No.ZR2019-MEE076)Zhoucun District School City Integration Development Project(Grant No.2020ZCXCZH01).
文摘In order to study the mechanical properties of the heterogeneous core plate of the wind turbine blade,a modeling method of the core plate based on displacement field variables is proposed.Firstly,the wind turbine blade core plate was modeled according to the theory of modeling heterogeneous material characteristics.Secondly,the three-point bending finite element model of the wind turbine blade core plate was solved by the display dynamic equation to obtain the deformation pattern and force-deformation relationship of the core plate.Finally,the three-point bending static test was conducted to compare with the finite element analysis.The test results show that:the damage form of the wind turbine blade core plate includes elasticity,yield,and failure stages.The main failure modes are plastic deformation,core material collapse,and panel-core delamination.The failure load measured by the test is 1.59 kN,which is basically consistent with the load-displacement result obtained by the simulation,with a difference of only 1.9%,which verifies the validity and reliability of the model.It provides data references for wind turbine blade structure design.
基金supported by the Major Science and Technology Project of PetroChina Company Limited“Research and application of key technologies for development of deep and ultra-deep gas reservoirs in Kuqa Depression(No.2018E-1803).
文摘The cretaceous gas reservoir in Kelasu Gas Field of the Tarim Basin is a rare ultra-deep and ultra-high pressure fractured tight sandstone gas reservoir where multi-scale discrete fractures of matrix,fracture and fault are developed,so its development cannot be conducted just based on static and dynamic reservoir description.In order to solve this problem,this paper establishes a numerical well test model of vertical wells based on matrix,fractures and faults(large fractures and small faults)by combining the random generation of natural fracture networks with the unstructured discrete fracture modeling method to break through the traditional continuous medium well test model.In addition,the model is solved by using the finite element method with mixed element,and the typical well test type curves under different random fracture networks are obtained.And the following research results are obtained.First,based on the observed data,the fracture network distribution modes of fractured tight sandstone gas reservoirs are classified into three categories.The influence of random generation of fracture networks on typical well test type curves is discussed.The results of discrete fracture well test model are compared with those of the traditional continuous medium well test model,and the applicable conditions of the traditional continuous medium well test model is determined.Second,there are great differences between the results of discrete fracture model and those of dual porosity medium model.1 The dual porosity medium model is a special case of the discrete fracture model,in which the fractures are evenly distributed within infinitely small spacing.Third,the characteristics of well test type curves under three fracture network distribution modes are discussed.The well test type curves that cannot be interpreted by the conventional dual/triple porosity continuous medium model are successfully interpreted by using the established well test interpretation model of random discrete fracture.The curve matching effect is ideal and the interpreted parameters are reasonable.In conclusion,the new model and the new method reveal the development mechanism of step-by-step production and coordinated gas supply between media of different scales,explain the development characteristics of large inter-well productivity difference and abnormal rapid inter-well pressure response,and provide a reference for the development of similar gas reservoirs.
