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.展开更多
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.展开更多
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.展开更多
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.展开更多
This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with...This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.展开更多
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.展开更多
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.展开更多
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.展开更多
Efficient and accurate strength analysis of bolted connections is essential in analyzing the integral thermal protection system（ITPS） of hypersonic vehicles, since the system bears severe loads and structural failur...Efficient and accurate strength analysis of bolted connections is essential in analyzing the integral thermal protection system（ITPS） of hypersonic vehicles, since the system bears severe loads and structural failures usually occur at the connections. Investigations of composite mechanical properties used in ITPS are still in progress as the architecture of the composites is complex. A new method is proposed in this paper for strength analysis of bolted connections by investigating the elastic behavior and failure strength of three-dimensional C/C orthogonal composites used in ITPS. In this method a multi-scale finite element method incorporating the global–local method is established to ensure high efficiency in macro-scale and precision in meso-scale in analysis.Simulation results reveal that predictions of material properties show reasonable accuracy compared with test results. And the multi-scale method can analyze the strength of connections efficiently and accurately.展开更多
In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and t...In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.展开更多
The degree of spatial similarity plays an important role in map generalization, yet there has been no quantitative research into it. To fill this gap, this study first defines map scale change and spatial similarity d...The degree of spatial similarity plays an important role in map generalization, yet there has been no quantitative research into it. To fill this gap, this study first defines map scale change and spatial similarity degree/relation in multi-scale map spaces and then proposes a model for calculating the degree of spatial similarity between a point cloud at one scale and its gener- alized counterpart at another scale. After validation, the new model features 16 points with map scale change as the x coordinate and the degree of spatial similarity as the y coordinate. Finally, using an application for curve fitting, the model achieves an empirical formula that can calculate the degree of spatial similarity using map scale change as the sole independent variable, and vice versa. This formula can be used to automate algorithms for point feature generalization and to determine when to terminate them during the generalization.展开更多
A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by com...A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by combining the proposed multi-scale crystal plasticity finite element method with experimental measurements.The complex local micro region load was progressively extracted from the simulation results of a macro model and applied to an established three-dimensional multi-grain microstructure model.Subsequently,the evolution histories of the grain shape,size,and orientation inside the adiabatic shear band were quantitatively simulated.The results corresponded closely to the experimental results obtained via transmission electron microscopy and precession electron diffraction.Furthermore,by calculating the grain rotation and temperature rise inside the adiabatic shear band,the microstructural softening and thermal softening effects of typical heavily-deformed α grains were successfully decoupled.The results revealed that the microstructural softening stress was triggered and then stabilized(in general)at a relatively high value.This indicated that the mechanical strength was lowered mainly by the grain orientation evolution or dynamic recrystallization occurring during early plastic deformation.Subsequently,thermal softening increased linearly and became the main softening mechanism.Noticeably,in the final stage,the thermal softening stress accounted for 78.4% of the total softening stress due to the sharp temperature increase,which inevitably leads to the stress collapse and potential failure of the alloy.展开更多
In a multiple voltage source converter(VSC)system,the nonlinear characteristics of phase-locked loops(PLLs)and their interactions have a significant influence on the synchronization stability of converters.In this pap...In a multiple voltage source converter(VSC)system,the nonlinear characteristics of phase-locked loops(PLLs)and their interactions have a significant influence on the synchronization stability of converters.In this paper,these influences are investigated from the perspective of the time domain.First,a novel time-domain model of the multi-VSC system is obtained by using a multi-scale method.On this basis,a stability criterion is proposed to assess the synchronization stability of the system.Then,the accuracy of the time-domain model and its stability criterion in various conditions are discussed.Moreover,the negative impact of the interaction on the system is quantified.Finally,the above theoretical analysis is also verified in the controller hardware-in-the-loop(CHIL)experiments.展开更多
Second-generation high-temperature superconducting(HTS)conductors,specifically rare earth-barium-copper-oxide(REBCO)coated conductor(CC)tapes,are promising candidates for high-energy and high-field superconducting app...Second-generation high-temperature superconducting(HTS)conductors,specifically rare earth-barium-copper-oxide(REBCO)coated conductor(CC)tapes,are promising candidates for high-energy and high-field superconducting applications.With respect to epoxy-impregnated REBCO composite magnets that comprise multilayer components,the thermomechanical characteristics of each component differ considerably under extremely low temperatures and strong electromagnetic fields.Traditional numerical models include homogenized orthotropic models,which simplify overall field calculation but miss detailed multi-physics aspects,and full refinement(FR)ones that are thorough but computationally demanding.Herein,we propose an extended multi-scale approach for analyzing the multi-field characteristics of an epoxy-impregnated composite magnet assembled by HTS pancake coils.This approach combines a global homogenization(GH)scheme based on the homogenized electromagnetic T-A model,a method for solving Maxwell's equations for superconducting materials based on the current vector potential T and the magnetic field vector potential A,and a homogenized orthotropic thermoelastic model to assess the electromagnetic and thermoelastic properties at the macroscopic scale.We then identify“dangerous regions”at the macroscopic scale and obtain finer details using a local refinement(LR)scheme to capture the responses of each component material in the HTS composite tapes at the mesoscopic scale.The results of the present GH-LR multi-scale approach agree well with those of the FR scheme and the experimental data in the literature,indicating that the present approach is accurate and efficient.The proposed GH-LR multi-scale approach can serve as a valuable tool for evaluating the risk of failure in large-scale HTS composite magnets.展开更多
Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for mult...Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.展开更多
We formulate a macroscopic particle modeling analysis of metallic materials (aluminum and copper, etc.) based on theoretical energy and atomic geome<span>tries derivable from their interatomic potential. In fact...We formulate a macroscopic particle modeling analysis of metallic materials (aluminum and copper, etc.) based on theoretical energy and atomic geome<span>tries derivable from their interatomic potential. In fact, particles in thi</span>s framework are presenting a large mass composed of huge collection of atoms and are interacting with each other. We can start from cohesive energy of metallic atoms and basic crystalline unit (e.g. face-centered cubic). Then, we can reach to interparticle (macroscopic) potential function which is presented by the analytical equation with terms of exponent of inter-particle distance, like a Lennard-Jones potential usually used in molecular dynamics simulation. Equation of motion for these macroscopic particles has dissipative term and fluctuation term, as well as the conservative term above, in order to express finite temperature condition. First, we determine the parameters needed in macroscopic potential function and check the reproduction of mechanical behavior in elastic regime. By using the present framework, we are able to carry out uniaxial loading simulation of aluminum rod. The method can also reproduce Young’s modulus and Poisson’s ratio as elastic behavior, though the result shows the dependency on division number of particles. Then, we proceed to try to include plasticity in this multi-scale framework. As a result, a realistic curve of stress-strain relation can be obtained for tensile and compressive loading and this new and simple framework of materials modeling has been confirmed to have certain effectiveness to be used in materials simulations. We also assess the effect of the order of loadings in opposite directions including yield and plastic states and find that an irreversible behavior depends on different response of the particle system between tensile and compressive loadings.展开更多
This paper examined how microstructure influences the homogenized thermal conductivity of cellular structures and revealed a surface-induced size-dependent effect.This effect is linked to the porous microstructural fe...This paper examined how microstructure influences the homogenized thermal conductivity of cellular structures and revealed a surface-induced size-dependent effect.This effect is linked to the porous microstructural features of cellular structures,which stems from the degree of porosity and the distri-bution of the pores.Unlike the phonon-driven surface effect at the nanoscale,the macro-scale surface mechanism in thermal cellular structures is found to be the microstructure-induced changes in the heat conduction path based on fully resolved 3D numerical simulations.The surface region is determined by the microstructure,characterized by the intrinsic length.With the coupling between extrinsic and intrinsic length scales under the surface mechanism,a surface-enriched multiscale method was devel-oped to accurately capture the complex size-dependent thermal conductivity.The principle of scale separation required by classical multiscale methods is not necessary to be satisfied by the proposed multiscale method.The significant potential of the surface-enriched multiscale method was demon-strated through simulations of the effective thermal conductivity of a thin-walled metamaterial struc-ture.The surface-enriched multiscale method offers higher accuracy compared with the classical multiscale method and superior efficiency over high-fidelity finite element methods.展开更多
This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method （IFM） and the homoge...This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method （IFM） and the homogenized initial system method （HISM）. In both methods, nonlinear variable loadings within time intervals are simulated using Chebyshev polynomials of the first kind before a direct integration is performed. Developed on the basis of the integral formula, the recurrence relationship of the integral computation suggested in this paper is combined with the Crout decomposed method to solve linear algebraic equations. In this way, the IFM based on Chebyshev polynomial of the first kind is constructed. Transforming the non-homogenous initial system to the homogeneous dynamic system, and developing a special scheme without dimensional expansion, the HISM based on Chebyshev polynomial of the first kind is able to avoid the matrix inversion operation. The accuracy of the time integration schemes is examined and compared with other commonly used schemes, and it is shown that a greater accuracy as well as less time consuming can be achieved. Two numerical examples are presented to demonstrate the applicability of these new methods.展开更多
The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise w...The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.展开更多
基金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.
基金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 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.
基金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.
基金Project supported by the National Natural Science Foundation of China (No. 10871130)the Ph. D.Programs Foundation of Ministry of Education of China (No. 20093127110005)the Shanghai Leading Academic Discipline Project (No. T0401)
文摘This paper proposes an inexact Newton method via the Lanczos decomposed technique for solving the box-constrained nonlinear systems. An iterative direction is obtained by solving an affine scaling quadratic model with the Lanczos decomposed technique. By using the interior backtracking line search technique, an acceptable trial step length is found along this direction. The global convergence and the fast local convergence rate of the proposed algorithm are established under some reasonable conditions. Furthermore, the results of the numerical experiments show the effectiveness of the pro- posed algorithm.
基金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 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.
文摘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.
基金co-supported by the Priority Academic Program Development of Jiangsu Higher Education Institutionsthe National Natural Science Foundation of China (No. 11302105)
文摘Efficient and accurate strength analysis of bolted connections is essential in analyzing the integral thermal protection system（ITPS） of hypersonic vehicles, since the system bears severe loads and structural failures usually occur at the connections. Investigations of composite mechanical properties used in ITPS are still in progress as the architecture of the composites is complex. A new method is proposed in this paper for strength analysis of bolted connections by investigating the elastic behavior and failure strength of three-dimensional C/C orthogonal composites used in ITPS. In this method a multi-scale finite element method incorporating the global–local method is established to ensure high efficiency in macro-scale and precision in meso-scale in analysis.Simulation results reveal that predictions of material properties show reasonable accuracy compared with test results. And the multi-scale method can analyze the strength of connections efficiently and accurately.
基金supported by Science and Technology on Reactor System Design Technology Laboratory,Chengdu,China(LRSDT2020106)
文摘In this study,a multi-physics and multi-scale coupling program,Fluent/KMC-sub/NDK,was developed based on the user-defined functions(UDF)of Fluent,in which the KMC-sub-code is a sub-channel thermal-hydraulic code and the NDK code is a neutron diffusion code.The coupling program framework adopts the"master-slave"mode,in which Fluent is the master program while NDK and KMC-sub are coupled internally and compiled into the dynamic link library(DLL)as slave codes.The domain decomposition method was adopted,in which the reactor core was simulated by NDK and KMC-sub,while the rest of the primary loop was simulated using Fluent.A simulation of the reactor shutdown process of M2LFR-1000 was carried out using the coupling program,and the code-to-code verification was performed with ATHLET,demonstrating a good agreement,with absolute deviation was smaller than 0.2%.The results show an obvious thermal stratification phenomenon during the shutdown process,which occurs 10 s after shutdown,and the change in thermal stratification phenomena is also captured by the coupling program.At the same time,the change in the neutron flux density distribution of the reactor was also obtained.
基金funded by the Natural Science Foundation Committee,China(41364001,41371435)
文摘The degree of spatial similarity plays an important role in map generalization, yet there has been no quantitative research into it. To fill this gap, this study first defines map scale change and spatial similarity degree/relation in multi-scale map spaces and then proposes a model for calculating the degree of spatial similarity between a point cloud at one scale and its gener- alized counterpart at another scale. After validation, the new model features 16 points with map scale change as the x coordinate and the degree of spatial similarity as the y coordinate. Finally, using an application for curve fitting, the model achieves an empirical formula that can calculate the degree of spatial similarity using map scale change as the sole independent variable, and vice versa. This formula can be used to automate algorithms for point feature generalization and to determine when to terminate them during the generalization.
基金financially supported by the National Natural Science Foundation of China(No.51571031)。
文摘A dynamic compression test was performed on α+β dual-phase titanium alloy Ti20C using a split Hopkinson pressure bar.The formation of adiabatic shear bands generated during the compression process was studied by combining the proposed multi-scale crystal plasticity finite element method with experimental measurements.The complex local micro region load was progressively extracted from the simulation results of a macro model and applied to an established three-dimensional multi-grain microstructure model.Subsequently,the evolution histories of the grain shape,size,and orientation inside the adiabatic shear band were quantitatively simulated.The results corresponded closely to the experimental results obtained via transmission electron microscopy and precession electron diffraction.Furthermore,by calculating the grain rotation and temperature rise inside the adiabatic shear band,the microstructural softening and thermal softening effects of typical heavily-deformed α grains were successfully decoupled.The results revealed that the microstructural softening stress was triggered and then stabilized(in general)at a relatively high value.This indicated that the mechanical strength was lowered mainly by the grain orientation evolution or dynamic recrystallization occurring during early plastic deformation.Subsequently,thermal softening increased linearly and became the main softening mechanism.Noticeably,in the final stage,the thermal softening stress accounted for 78.4% of the total softening stress due to the sharp temperature increase,which inevitably leads to the stress collapse and potential failure of the alloy.
基金supported by the Science and Technology Project of State Grid Corporation of China(5400-202199281A-0-0-00).
文摘In a multiple voltage source converter(VSC)system,the nonlinear characteristics of phase-locked loops(PLLs)and their interactions have a significant influence on the synchronization stability of converters.In this paper,these influences are investigated from the perspective of the time domain.First,a novel time-domain model of the multi-VSC system is obtained by using a multi-scale method.On this basis,a stability criterion is proposed to assess the synchronization stability of the system.Then,the accuracy of the time-domain model and its stability criterion in various conditions are discussed.Moreover,the negative impact of the interaction on the system is quantified.Finally,the above theoretical analysis is also verified in the controller hardware-in-the-loop(CHIL)experiments.
基金Project supported by the National Natural Science Foundation of China(Nos.11932008 and 12272156)the Fundamental Research Funds for the Central Universities(No.lzujbky-2022-kb06)+1 种基金the Gansu Science and Technology ProgramLanzhou City’s Scientific Research Funding Subsidy to Lanzhou University of China。
文摘Second-generation high-temperature superconducting(HTS)conductors,specifically rare earth-barium-copper-oxide(REBCO)coated conductor(CC)tapes,are promising candidates for high-energy and high-field superconducting applications.With respect to epoxy-impregnated REBCO composite magnets that comprise multilayer components,the thermomechanical characteristics of each component differ considerably under extremely low temperatures and strong electromagnetic fields.Traditional numerical models include homogenized orthotropic models,which simplify overall field calculation but miss detailed multi-physics aspects,and full refinement(FR)ones that are thorough but computationally demanding.Herein,we propose an extended multi-scale approach for analyzing the multi-field characteristics of an epoxy-impregnated composite magnet assembled by HTS pancake coils.This approach combines a global homogenization(GH)scheme based on the homogenized electromagnetic T-A model,a method for solving Maxwell's equations for superconducting materials based on the current vector potential T and the magnetic field vector potential A,and a homogenized orthotropic thermoelastic model to assess the electromagnetic and thermoelastic properties at the macroscopic scale.We then identify“dangerous regions”at the macroscopic scale and obtain finer details using a local refinement(LR)scheme to capture the responses of each component material in the HTS composite tapes at the mesoscopic scale.The results of the present GH-LR multi-scale approach agree well with those of the FR scheme and the experimental data in the literature,indicating that the present approach is accurate and efficient.The proposed GH-LR multi-scale approach can serve as a valuable tool for evaluating the risk of failure in large-scale HTS composite magnets.
文摘Seismic energy decays while propagating subsurface, which may reduce the resolution of seismic data. This paper studies the method of seismic energy dispersion compensation which provides the basic principles for multi-scale morphology and the spectrum simulation method. These methods are applied in seismic energy compensation. First of all, the seismic data is decomposed into multiple scales and the effective frequency bandwidth is selectively broadened for some scales by using a spectrum simulation method. In this process, according to the amplitude spectrum of each scale, the best simulation range is selected to simulate the middle and low frequency components to ensure the authenticity of the simulation curve which is calculated by the median method, and the high frequency component is broadened. Finally, these scales are reconstructed with reasonable coefficients, and the compensated seismic data can be obtained. Examples are shown to illustrate the feasibility of the energy compensation method.
文摘We formulate a macroscopic particle modeling analysis of metallic materials (aluminum and copper, etc.) based on theoretical energy and atomic geome<span>tries derivable from their interatomic potential. In fact, particles in thi</span>s framework are presenting a large mass composed of huge collection of atoms and are interacting with each other. We can start from cohesive energy of metallic atoms and basic crystalline unit (e.g. face-centered cubic). Then, we can reach to interparticle (macroscopic) potential function which is presented by the analytical equation with terms of exponent of inter-particle distance, like a Lennard-Jones potential usually used in molecular dynamics simulation. Equation of motion for these macroscopic particles has dissipative term and fluctuation term, as well as the conservative term above, in order to express finite temperature condition. First, we determine the parameters needed in macroscopic potential function and check the reproduction of mechanical behavior in elastic regime. By using the present framework, we are able to carry out uniaxial loading simulation of aluminum rod. The method can also reproduce Young’s modulus and Poisson’s ratio as elastic behavior, though the result shows the dependency on division number of particles. Then, we proceed to try to include plasticity in this multi-scale framework. As a result, a realistic curve of stress-strain relation can be obtained for tensile and compressive loading and this new and simple framework of materials modeling has been confirmed to have certain effectiveness to be used in materials simulations. We also assess the effect of the order of loadings in opposite directions including yield and plastic states and find that an irreversible behavior depends on different response of the particle system between tensile and compressive loadings.
基金supported by the National Key Research and Development Program of China(Grant No.2021YFB1714600)the National Natural Science Foundation of China(Grant No.52175095)the Young Top-Notch Talent Cultivation Program of Hubei Province of China.
文摘This paper examined how microstructure influences the homogenized thermal conductivity of cellular structures and revealed a surface-induced size-dependent effect.This effect is linked to the porous microstructural features of cellular structures,which stems from the degree of porosity and the distri-bution of the pores.Unlike the phonon-driven surface effect at the nanoscale,the macro-scale surface mechanism in thermal cellular structures is found to be the microstructure-induced changes in the heat conduction path based on fully resolved 3D numerical simulations.The surface region is determined by the microstructure,characterized by the intrinsic length.With the coupling between extrinsic and intrinsic length scales under the surface mechanism,a surface-enriched multiscale method was devel-oped to accurately capture the complex size-dependent thermal conductivity.The principle of scale separation required by classical multiscale methods is not necessary to be satisfied by the proposed multiscale method.The significant potential of the surface-enriched multiscale method was demon-strated through simulations of the effective thermal conductivity of a thin-walled metamaterial struc-ture.The surface-enriched multiscale method offers higher accuracy compared with the classical multiscale method and superior efficiency over high-fidelity finite element methods.
基金Hunan Provincial Natural Science Foundation Under Grant No.02JJY2085
文摘This paper introduces two new types of precise integration methods based on Chebyshev polynomial of the first kind for dynamic response analysis of structures, namely the integral formula method （IFM） and the homogenized initial system method （HISM）. In both methods, nonlinear variable loadings within time intervals are simulated using Chebyshev polynomials of the first kind before a direct integration is performed. Developed on the basis of the integral formula, the recurrence relationship of the integral computation suggested in this paper is combined with the Crout decomposed method to solve linear algebraic equations. In this way, the IFM based on Chebyshev polynomial of the first kind is constructed. Transforming the non-homogenous initial system to the homogeneous dynamic system, and developing a special scheme without dimensional expansion, the HISM based on Chebyshev polynomial of the first kind is able to avoid the matrix inversion operation. The accuracy of the time integration schemes is examined and compared with other commonly used schemes, and it is shown that a greater accuracy as well as less time consuming can be achieved. Two numerical examples are presented to demonstrate the applicability of these new methods.
文摘The precise integration method proposed for linear time-invariant homogeneous dynamic systems can provide accurate numerical results that approach an exact solution at integration points. However, difficulties arise when the algorithm is used for non-homogeneous dynamic systems due to the inverse matrix calculation required. In this paper, the structural dynamic equalibrium equations are converted into a special form, the inverse matrix calculation is replaced by the Crout decomposition method to solve the dynamic equilibrium equations, and the precise integration method without the inverse matrix calculation is obtained. The new algorithm enhances the present precise integration method by improving both the computational accuracy and efficiency. Two numerical examples are given to demonstrate the validity and efficiency of the proposed algorithm.
