Influenza,an acute respiratory infectious disease caused by the influenza virus,exhibits distinct seasonal patterns in China,with peak activity occurring in winter and spring in northern regions,and in winter and summ...Influenza,an acute respiratory infectious disease caused by the influenza virus,exhibits distinct seasonal patterns in China,with peak activity occurring in winter and spring in northern regions,and in winter and summer in southern areas[1].The World Health Organization(WHO)emphasizes that early warning and epidemic intensity assessments are critical public health strategies for influenza prevention and control.Internet-based flu surveillance,with real-time data and low costs,effectively complements traditional methods.The Baidu Search Index,which reflects flu-related queries,strongly correlates with influenza trends,aiding in regional activity assessment and outbreak tracking[2].展开更多
In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This st...In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.展开更多
In the process of large-scale,grid-connected wind power operations,it is important to establish an accurate probability distribution model for wind farm fluctuations.In this study,a wind power fluctuation modeling met...In the process of large-scale,grid-connected wind power operations,it is important to establish an accurate probability distribution model for wind farm fluctuations.In this study,a wind power fluctuation modeling method is proposed based on the method of moving average and adaptive nonparametric kernel density estimation(NPKDE)method.Firstly,the method of moving average is used to reduce the fluctuation of the sampling wind power component,and the probability characteristics of the modeling are then determined based on the NPKDE.Secondly,the model is improved adaptively,and is then solved by using constraint-order optimization.The simulation results show that this method has a better accuracy and applicability compared with the modeling method based on traditional parameter estimation,and solves the local adaptation problem of traditional NPKDE.展开更多
Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatia...Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatial domain using generalized moving least square method （GMLSM） and then, discrete equations of motion based on Lagrange＇s equation are obtained. Maximum deflection and bending moments are considered as the important design parameters. The design parameter spectra in terms of mass weight and velocity of the moving mass are presented for multispan viscoelastic beams as well as various values of relaxation rate and beam span number. A reasonable good agreement is achieved between the results of the proposed solution and those obtained by other researchers. The results indicate that, although the load inertia effects in beams with higher span number would be intensified for higher levels of moving mass velocity, the maximum values of design parameters would increase either. Moreover, the possibility of mass separation is shown to be more critical as the span number of the beam increases. This fact also violates the linear relation between the mass weight of the moving load and the associated design parameters, especially for high moving mass velocities. However, as the relaxation rate of the beam material increases, the load inertia effects as well as the possibility of moving mass separation reduces.展开更多
A size-dependent computational approach for bending,free vibration and buckling analyses of isotropic and sandwich functionally graded(FG)microplates is in this study presented.We consider both shear deformation and s...A size-dependent computational approach for bending,free vibration and buckling analyses of isotropic and sandwich functionally graded(FG)microplates is in this study presented.We consider both shear deformation and small scale effects through the generalized higher order shear deformation theory and modified couple stress theory(MCST).The present model only retains a single material length scale parameter for capturing properly size effects.A rule of mixture is used to model material properties varying through the thickness of plates.The principle of virtual work is used to derive the discrete system equations which are approximated by moving Kriging interpolation(MKI)meshfree method.Numerical examples consider the inclusions of geometrical parameters,volume fraction,boundary conditions and material length scale parameter.Reliability and effectiveness of the present method are confirmed through numerical results.展开更多
With the east section of the Changji sag Zhunger Basin as a case study, both a principal curvature method and a moving least square method are elaborated. The moving least square method is introduced, for the first ti...With the east section of the Changji sag Zhunger Basin as a case study, both a principal curvature method and a moving least square method are elaborated. The moving least square method is introduced, for the first time, to fit a stratum surface. The results show that, using the same-degree base function, compared with a traditional least square method, the moving least square method can produce lower fitting errors, the fitting surface can describe the morphological characteristics of stratum surfaces more accurately and the principal curvature values vary within a wide range and may be more suitable for the prediction of the distribution of structural fractures. The moving least square method could be useful in curved surface fitting and stratum curvature analysis.展开更多
Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian g...Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.展开更多
Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Latti...Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.展开更多
Current multiscale topology optimization restricts the solution space by enforcing the use of a few repetitive microstructures that are predetermined,and thus lack the ability for structural concerns like buckling str...Current multiscale topology optimization restricts the solution space by enforcing the use of a few repetitive microstructures that are predetermined,and thus lack the ability for structural concerns like buckling strength,robustness,and multi-functionality.Therefore,in this paper,a new multiscale concurrent topology optimization design,referred to as the self-consistent analysis-based moving morphable component(SMMC)method,is proposed.Compared with the conventional moving morphable component method,the proposed method seeks to optimize both material and structure simultaneously by explicitly designing both macrostructure and representative volume element(RVE)-level microstructures.Numerical examples with transducer design requirements are provided to demonstrate the superiority of the SMMC method in comparison to traditional methods.The proposed method has broad impact in areas of integrated industrial manufacturing design:to solve for the optimized macro and microstructures under the objective function and constraints,to calculate the structural response efficiently using a reduced-order model:self-consistent analysis,and to link the SMMC method to manufacturing(industrial manufacturing or additive manufacturing)based on the design requirements and application areas.展开更多
In order to improve depth extraction accuracy, a method using moving array lenslet technique(MALT) in pickup stage is proposed, which can decrease the depth interval caused by pixelation. In this method, the lenslet a...In order to improve depth extraction accuracy, a method using moving array lenslet technique(MALT) in pickup stage is proposed, which can decrease the depth interval caused by pixelation. In this method, the lenslet array is moved along the horizontal and vertical directions simultaneously for N times in a pitch to get N sets of elemental images. Computational integral imaging reconstruction method for MALT is taken to obtain the slice images of the 3 D scene, and the sum modulus(SMD) blur metric is taken on these slice images to achieve the depth information of the 3 D scene. Simulation and optical experiments are carried out to verify the feasibility of this method.展开更多
This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of ...This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange＇s equations via reproducing kernel particle method （RKPM）. For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.展开更多
In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the II...In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.展开更多
The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （B...The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.展开更多
Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving m...Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense.As the state-of-the-art moving mesh method,the variational mesh adaptation approach has been introduced to CFD calculation.However,quickly estimating the flow field on the updated meshes during the iterative algorithm is challenging.A mesh optimization method,which embeds a machine learning regression model into the variational mesh adaptation,is proposed.The regression model captures the mapping between the initial mesh nodes and the flow field,so that the variational method could move mesh nodes iteratively by solving the mesh functional which is built from the estimated flow field on the updated mesh via the regression model.After the optimization,the density of the nodes in the high gradient area increases while the density in the low gradient area decreases.Benchmark examples are first used to verify the feasibility and effectiveness of the proposed method.And then we use the steady subsonic and transonic flows over cylinder and NACA0012 airfoil on unstructured triangular meshes to test our method.Results show that the proposed method significantly improves the accuracy of the local flow features on the adaptive meshes.Our work indicates that the proposed mesh optimization approach is promising for improving the accuracy and efficiency of CFD computation.展开更多
In naval engineering and offshore industry,the fluid-structure interaction(FSI)problem is a very common problem,and water entry is a very representative one.The hydroelasticity effects due to slamming are of great int...In naval engineering and offshore industry,the fluid-structure interaction(FSI)problem is a very common problem,and water entry is a very representative one.The hydroelasticity effects due to slamming are of great interest.In this paper,the water entry problem is simulated by the moving particle semi-implicit&finite element method(MPS-FEM)coupled method.The MPS method is used for the fluid because it is very suitable for the violent free-surface flow.The structure domain is solved by the FEM method because of the maturity in solving structural motion and deformation.The water entry of a rigid cylinder is numerically studied first and the results show good agreements with previous published data.After that,variable analysis is conducted in the water entry simulation of an elastic cylinder,including the structural elasticity and impact velocity.展开更多
Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering application...Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.展开更多
A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large edd...A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.展开更多
Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential proble...Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square（MLS） approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin（MLPG） method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.展开更多
In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral c...In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.展开更多
Special transmission 3D model simulation must be based on surface discretization and reconstruction, but special transmission usually has complicated tooth shape and movement, so present software can＇t provide techni...Special transmission 3D model simulation must be based on surface discretization and reconstruction, but special transmission usually has complicated tooth shape and movement, so present software can＇t provide technical support for special transmission 3D model simulation. Currently, theoretical calculation and experimental method are difficult to exactly solve special transmission contact analysis problem. How to reduce calculation and computer memories consume and meet calculation precision is key to resolve special transmission contact analysis problem. According to 3D model simulation and surface reconstruction of quasi ellipsoid gear is difficulty, this paper employes meshless local Petrov-Galerkin （MLPG） method. In order to reduce calculation and computer memories consume, we disperse tooth mesh into finite points--sparseness points cloud or grid mesh, and then we do interpolation reconstruction in some necessary place of the 3D surface model during analysis. Moving least square method （MLSM） is employed for tooth mesh interpolation reconstruction, there are some advantages to do interpolation by means of MLSM, such as high precision, good flexibility and no require of tooth mesh discretization into units. We input the quasi ellipsoid gear reconstruction model into simulation software, we complete tooth meshing simulation. Simulation transmission ratio during meshing period was obtained, compared with theoretical transmission ratio, the result inosculate preferably. The method using curve reconstruction realizes surface reconstruction, reduce simulation calculation enormously, so special gears simulation can be realized by minitype computer. The method provides a novel solution for special transmission 3D model simulation analysis and contact analysis.展开更多
基金supported by the National Key Research and Development Program of China(Project No.2023YFC2307500).
文摘Influenza,an acute respiratory infectious disease caused by the influenza virus,exhibits distinct seasonal patterns in China,with peak activity occurring in winter and spring in northern regions,and in winter and summer in southern areas[1].The World Health Organization(WHO)emphasizes that early warning and epidemic intensity assessments are critical public health strategies for influenza prevention and control.Internet-based flu surveillance,with real-time data and low costs,effectively complements traditional methods.The Baidu Search Index,which reflects flu-related queries,strongly correlates with influenza trends,aiding in regional activity assessment and outbreak tracking[2].
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China(Grant No.10902076)the Natural Science Foundation of Shanxi Province of China(Grant No.2007011009)+1 种基金the Scientific Research and Development Program of the Shanxi Higher Education Institutions(Grant No.20091131)the Doctoral Startup Foundation of Taiyuan University of Science and Technology(Grant No.200708)
文摘In this paper, a meshfree boundary integral equation （BIE） method, called the moving Kriging interpolation- based boundary node method （MKIBNM）, is developed for solving two-dimensional potential problems. This study combines the DIE method with the moving Kriging interpolation to present a boundary-type meshfree method, and the corresponding formulae of the MKIBNM are derived. In the present method, the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker＇s delta property, then the boundary conditions can be imposed directly and easily. To verify the accuracy and stability of the present formulation, three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.
基金supported by Science and Technology project of the State Grid Corporation of China“Research on Active Development Planning Technology and Comprehensive Benefit Analysis Method for Regional Smart Grid Comprehensive Demonstration Zone”National Natural Science Foundation of China(51607104)
文摘In the process of large-scale,grid-connected wind power operations,it is important to establish an accurate probability distribution model for wind farm fluctuations.In this study,a wind power fluctuation modeling method is proposed based on the method of moving average and adaptive nonparametric kernel density estimation(NPKDE)method.Firstly,the method of moving average is used to reduce the fluctuation of the sampling wind power component,and the probability characteristics of the modeling are then determined based on the NPKDE.Secondly,the model is improved adaptively,and is then solved by using constraint-order optimization.The simulation results show that this method has a better accuracy and applicability compared with the modeling method based on traditional parameter estimation,and solves the local adaptation problem of traditional NPKDE.
文摘Dynamic response of multispan viscoelastic thin beams subjected to a moving mass is studied by an efficient numerical method in some detail. To this end, the unknown parameters of the problem are discretized in spatial domain using generalized moving least square method （GMLSM） and then, discrete equations of motion based on Lagrange＇s equation are obtained. Maximum deflection and bending moments are considered as the important design parameters. The design parameter spectra in terms of mass weight and velocity of the moving mass are presented for multispan viscoelastic beams as well as various values of relaxation rate and beam span number. A reasonable good agreement is achieved between the results of the proposed solution and those obtained by other researchers. The results indicate that, although the load inertia effects in beams with higher span number would be intensified for higher levels of moving mass velocity, the maximum values of design parameters would increase either. Moreover, the possibility of mass separation is shown to be more critical as the span number of the beam increases. This fact also violates the linear relation between the mass weight of the moving load and the associated design parameters, especially for high moving mass velocities. However, as the relaxation rate of the beam material increases, the load inertia effects as well as the possibility of moving mass separation reduces.
文摘A size-dependent computational approach for bending,free vibration and buckling analyses of isotropic and sandwich functionally graded(FG)microplates is in this study presented.We consider both shear deformation and small scale effects through the generalized higher order shear deformation theory and modified couple stress theory(MCST).The present model only retains a single material length scale parameter for capturing properly size effects.A rule of mixture is used to model material properties varying through the thickness of plates.The principle of virtual work is used to derive the discrete system equations which are approximated by moving Kriging interpolation(MKI)meshfree method.Numerical examples consider the inclusions of geometrical parameters,volume fraction,boundary conditions and material length scale parameter.Reliability and effectiveness of the present method are confirmed through numerical results.
基金Projects 2007CB209405 and 2002CB412702 supported by the National Basic Research Program of ChinaKZCX2-YW-113 by the Important Directive Item of the Knowledge Innovation Project of Chinese Academy of Sciences 40772100 by the National Natural Science Foundation of China
文摘With the east section of the Changji sag Zhunger Basin as a case study, both a principal curvature method and a moving least square method are elaborated. The moving least square method is introduced, for the first time, to fit a stratum surface. The results show that, using the same-degree base function, compared with a traditional least square method, the moving least square method can produce lower fitting errors, the fitting surface can describe the morphological characteristics of stratum surfaces more accurately and the principal curvature values vary within a wide range and may be more suitable for the prediction of the distribution of structural fractures. The moving least square method could be useful in curved surface fitting and stratum curvature analysis.
文摘Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.
基金WJD,JYZ,CLC,ZX,and ZGY were supported by the National Natural Science Foundation of China(Grant Number 51705143)the Education Department of Hunan Province(Grant Number 22B0464)the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant Number QL20230249).
文摘Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.
文摘Current multiscale topology optimization restricts the solution space by enforcing the use of a few repetitive microstructures that are predetermined,and thus lack the ability for structural concerns like buckling strength,robustness,and multi-functionality.Therefore,in this paper,a new multiscale concurrent topology optimization design,referred to as the self-consistent analysis-based moving morphable component(SMMC)method,is proposed.Compared with the conventional moving morphable component method,the proposed method seeks to optimize both material and structure simultaneously by explicitly designing both macrostructure and representative volume element(RVE)-level microstructures.Numerical examples with transducer design requirements are provided to demonstrate the superiority of the SMMC method in comparison to traditional methods.The proposed method has broad impact in areas of integrated industrial manufacturing design:to solve for the optimized macro and microstructures under the objective function and constraints,to calculate the structural response efficiently using a reduced-order model:self-consistent analysis,and to link the SMMC method to manufacturing(industrial manufacturing or additive manufacturing)based on the design requirements and application areas.
基金supported by the National Natural Science Foundation of China(Nos.11474169 and 61675100)the Tianjin Natural Science Foundation(No.15JCYBJC16900)
文摘In order to improve depth extraction accuracy, a method using moving array lenslet technique(MALT) in pickup stage is proposed, which can decrease the depth interval caused by pixelation. In this method, the lenslet array is moved along the horizontal and vertical directions simultaneously for N times in a pitch to get N sets of elemental images. Computational integral imaging reconstruction method for MALT is taken to obtain the slice images of the 3 D scene, and the sum modulus(SMD) blur metric is taken on these slice images to achieve the depth information of the 3 D scene. Simulation and optical experiments are carried out to verify the feasibility of this method.
文摘This paper deals with the capabilities of linear and nonlinear beam theories in predicting the dynamic response of an elastically supported thin beam traversed by a moving mass. To this end, the discrete equations of motion are developed based on Lagrange＇s equations via reproducing kernel particle method （RKPM）. For a particular case of a simply supported beam, Galerkin method is also employed to verify the results obtained by RKPM, and a reasonably good agreement is achieved. Variations of the maximum dynamic deflection and bending moment associated with the linear and nonlinear beam theories are investigated in terms of moving mass weight and velocity for various beam boundary conditions. It is demonstrated that for majority of the moving mass velocities, the differences between the results of linear and nonlinear analyses become remarkable as the moving mass weight increases, particularly for high levels of moving mass velocity. Except for the cantilever beam, the nonlinear beam theory predicts higher possibility of moving mass separation from the base beam compared to the linear one. Furthermore, the accuracy levels of the linear beam theory are determined for thin beams under large deflections and small rotations as a function of moving mass weight and velocity in various boundary conditions.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)
文摘In this paper, an improved interpolating moving least-square （IIMLS） method is presented. The shape function of the IIMLS method satisfies the property of the Kronecker 5 function. The weight function used in the IIMLS method is nonsingular. Then the IIMLS method can overcome the difficulties caused by the singularity of the weight function in the IMLS method. The number of unknown coefficients in the trial function of the IIMLS method is less than that of the moving least-square （MLS） approximation. Then by combining the IIMLS method with the Galerkin weak form of the potential problem, the improved interpolating element-free Galerkin （IIEFG） method for two-dimensional potential problems is presented. Compared with the conventional element-free Galerkin （EFG） method, the IIEFG method can directly use the essential boundary conditions. Then the IIEFG method has higher accuracy. For demonstration, three numerical examples are solved using the IIEFG method.
基金Project supported by the National Natural Science Foundation of China (Grant No 10871124)Innovation Program of Shanghai Municipal Education Commission (Grant No 09ZZ99)Shanghai Leading Academic Discipline Project (Grant No J50103)
文摘The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
基金co-supported by the Key Laboratory of Aerodynamic Noise Control,China(No.ANCL20190103)the State Key Laboratory of Aerodynamics,China(No.SKLA20180102)the Aeronautical Science Foundation of China(Nos.2018ZA52002 and 2019ZA052011)。
文摘Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation.In light of the introduced large amount of computational costs for many adaptive mesh methods,moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense.As the state-of-the-art moving mesh method,the variational mesh adaptation approach has been introduced to CFD calculation.However,quickly estimating the flow field on the updated meshes during the iterative algorithm is challenging.A mesh optimization method,which embeds a machine learning regression model into the variational mesh adaptation,is proposed.The regression model captures the mapping between the initial mesh nodes and the flow field,so that the variational method could move mesh nodes iteratively by solving the mesh functional which is built from the estimated flow field on the updated mesh via the regression model.After the optimization,the density of the nodes in the high gradient area increases while the density in the low gradient area decreases.Benchmark examples are first used to verify the feasibility and effectiveness of the proposed method.And then we use the steady subsonic and transonic flows over cylinder and NACA0012 airfoil on unstructured triangular meshes to test our method.Results show that the proposed method significantly improves the accuracy of the local flow features on the adaptive meshes.Our work indicates that the proposed mesh optimization approach is promising for improving the accuracy and efficiency of CFD computation.
基金the National Key Research and Development Program of China(Grant No.2019YFB1704200)the National Natural Science Foundation of China(Grant Nos.51879159 and 52131102).
文摘In naval engineering and offshore industry,the fluid-structure interaction(FSI)problem is a very common problem,and water entry is a very representative one.The hydroelasticity effects due to slamming are of great interest.In this paper,the water entry problem is simulated by the moving particle semi-implicit&finite element method(MPS-FEM)coupled method.The MPS method is used for the fluid because it is very suitable for the violent free-surface flow.The structure domain is solved by the FEM method because of the maturity in solving structural motion and deformation.The water entry of a rigid cylinder is numerically studied first and the results show good agreements with previous published data.After that,variable analysis is conducted in the water entry simulation of an elastic cylinder,including the structural elasticity and impact velocity.
基金supported by the National Natural Science Foundation of China(Grant No.50979059)
文摘Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51479116 and 11272213)
文摘A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 11102125)
文摘Based on the complex variable moving least-square（CVMLS） approximation and a local symmetric weak form,the complex variable meshless local Petrov-Galerkin（CVMLPG） method of solving two-dimensional potential problems is presented in this paper.In the present formulation,the trial function of a two-dimensional problem is formed with a one-dimensional basis function.The number of unknown coefficients in the trial function of the CVMLS approximation is less than that in the trial function of the moving least-square（MLS） approximation.The essential boundary conditions are imposed by the penalty method.The main advantage of this approach over the conventional meshless local Petrov-Galerkin（MLPG） method is its computational efficiency.Several numerical examples are presented to illustrate the implementation and performance of the present CVMLPG method.
基金The research for this paper was supported by(1)the National Natural Science Foundation of China(Grants Nos.51708429,51708428)the Open Projects Foundation(Grant No.2017-04-GF)of State Key Laboratory for Health and Safety of Bridge Structures+1 种基金Wuhan Institute of Technology Science Found(Grant No.K201734)the science and technology projects of Wuhan Urban and Rural Construction Bureau(Grants Nos.201831,201919).
文摘In this paper,general interpolating isogeometric boundary node method(IIBNM)and isogeometric boundary element method(IBEM)based on parameter space are proposed for 2D elasticity problems.In both methods,the integral cells and elements are defined in parameter space,which can reproduce the geometry exactly at all the stages.In IIBNM,the improved interpolating moving leastsquare method(IIMLS)is applied for field approximation and the shape functions have the delta function property.The Lagrangian basis functions are used for field approximation in IBEM.Thus,the boundary conditions can be imposed directly in both methods.The shape functions are defined in 1D parameter space and no curve length needs to be computed.Besides,most methods for the treatment of the singular integrals in the boundary element method can be applied in IIBNM and IBEM directly.Numerical examples have demonstrated the accuracy of the proposed methods.
基金supported by National Natural Science Foundation of China (Grant No. 50905049)Heilongjiang Provincial International Cooperation Project of China (WB06A06)+1 种基金Heilongjiang Provincial Programs for Science and Technology Development of China (GC09A524)Heilongjiang Provincial Postdoctoral Science Foundation of China (LBH-Z09189)
文摘Special transmission 3D model simulation must be based on surface discretization and reconstruction, but special transmission usually has complicated tooth shape and movement, so present software can＇t provide technical support for special transmission 3D model simulation. Currently, theoretical calculation and experimental method are difficult to exactly solve special transmission contact analysis problem. How to reduce calculation and computer memories consume and meet calculation precision is key to resolve special transmission contact analysis problem. According to 3D model simulation and surface reconstruction of quasi ellipsoid gear is difficulty, this paper employes meshless local Petrov-Galerkin （MLPG） method. In order to reduce calculation and computer memories consume, we disperse tooth mesh into finite points--sparseness points cloud or grid mesh, and then we do interpolation reconstruction in some necessary place of the 3D surface model during analysis. Moving least square method （MLSM） is employed for tooth mesh interpolation reconstruction, there are some advantages to do interpolation by means of MLSM, such as high precision, good flexibility and no require of tooth mesh discretization into units. We input the quasi ellipsoid gear reconstruction model into simulation software, we complete tooth meshing simulation. Simulation transmission ratio during meshing period was obtained, compared with theoretical transmission ratio, the result inosculate preferably. The method using curve reconstruction realizes surface reconstruction, reduce simulation calculation enormously, so special gears simulation can be realized by minitype computer. The method provides a novel solution for special transmission 3D model simulation analysis and contact analysis.
