The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) metho...The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) method to address these limitations. In the vLEC method, the waverider is directly designed from the preassigned leading edge and the variable shock angle distribution based on the Leading-Edge Cone (LEC) concept. Since the vLEC method is an approximate method, two test waveriders are designed and evaluated using numerical simulations to validate the shock design accuracy and the effectiveness of the vLEC method. The results show that the shocks of the test waveriders coincide well with the preassigned positions. Furthermore, four specifically designed application cases are conducted to analyze the performance benefits of the vLEC waveriders. The results of these cases indicate that, due to their variable shock angle distributions, the vLEC waveriders exhibit higher lift-to-drag ratios and better longitudinal static stability than conventional waveriders. Additionally, the vLEC waveriders demonstrate superior volumetric capacities near the symmetry plane, albeit with a minor decrease in volumetric efficiency.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial ...The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial thin film growth models.Under the stepratio condition 0<τ_(n)/τ_(n-1)<4.864,the modified energy dissipation law is proven at the discrete levels with regardless of time step size.Nu‐merical experiments are presented to demonstrate the accuracy and efficiency of the proposed numerical scheme.展开更多
Mixed-variable problems are inevitable in engineering. However, few researches pay attention to discrete variables. This paper proposed a mixed-variable experimental design method (ODCD): first, the design variables w...Mixed-variable problems are inevitable in engineering. However, few researches pay attention to discrete variables. This paper proposed a mixed-variable experimental design method (ODCD): first, the design variables were divided into discrete variables and continuous variables;then, the DVD method was employed for handling discrete variables, the LHD method was applied for continuous variables, and finally, a Columnwise-Pairwise Algorithm was used for the overall optimization of the design matrix. Experimental results demonstrated that the ODCD method outperforms in terms of the sample space coverage performance.展开更多
Government credibility is an important asset of contemporary national governance, an important criterion for evaluating government legitimacy, and a key factor in measuring the effectiveness of government governance. ...Government credibility is an important asset of contemporary national governance, an important criterion for evaluating government legitimacy, and a key factor in measuring the effectiveness of government governance. In recent years, researchers’ research on government credibility has mostly focused on exploring theories and mechanisms, with little empirical research on this topic. This article intends to apply variable selection models in the field of social statistics to the issue of government credibility, in order to achieve empirical research on government credibility and explore its core influencing factors from a statistical perspective. Specifically, this article intends to use four regression-analysis-based methods and three random-forest-based methods to study the influencing factors of government credibility in various provinces in China, and compare the performance of these seven variable selection methods in different dimensions. The research results show that there are certain differences in simplicity, accuracy, and variable importance ranking among different variable selection methods, which present different importance in the study of government credibility issues. This study provides a methodological reference for variable selection models in the field of social science research, and also offers a multidimensional comparative perspective for analyzing the influencing factors of government credibility.展开更多
In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least s...In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.展开更多
Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patte...Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.展开更多
This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of...This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of structural collapse.The ET method involves dynamic analysis of a structure under an intensifying record over time.While conventional ET methods can determine the median of the structural response,they lack the ability to calculate its dispersion.To address this limitation,the present study utilizes ET analysis and single-degree-of-freedom(SDOF)systems to develop a method that considers the record-to-record variability for calculating the probability distribution of structural response.The accuracy of this method is evaluated by comparing it with the incremental dynamic analysis(IDA)method using special moment frames.The results demonstrate that the proposed method achieves a reasonably accurate estimation of dispersion while significantly reducing the computational burden(by approximately 95%)compared to the IDA method.展开更多
In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order a...In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.展开更多
This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence ...This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.展开更多
In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is pres...In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.展开更多
The eigenvalue problem of an infinite-dimensional Hamiltonian operator appearing in the isotropic plane magnetoelectroelastic solids is studied. First, all the eigenvalues and their eigenfunctions in a rectangular dom...The eigenvalue problem of an infinite-dimensional Hamiltonian operator appearing in the isotropic plane magnetoelectroelastic solids is studied. First, all the eigenvalues and their eigenfunctions in a rectangular domain are solved directly. Then the completeness of the eigenfunction system is proved, which offers a theoretic guarantee of the feasibility of variable separation method based on a Hamiltonian system for isotropic plane magnetoelectroelastic solids. Finally, the general solution for the equation in the rectangular domain is obtained by using the symplectic Fourier expansion method.展开更多
Studies show that different geometries of a Variable Cycle Engine(VCE)can be adjusted during the transient stage of the engine operation to improve the engine performance.However,this improvement increases the complex...Studies show that different geometries of a Variable Cycle Engine(VCE)can be adjusted during the transient stage of the engine operation to improve the engine performance.However,this improvement increases the complexity of the acceleration and deceleration control schedule.In order to resolve this problem,the Transient-state Reverse Method(TRM)is established in the present study based on the Steady-state Reverse Method(SRM)and the Virtual Power Extraction Method(VPEM).The state factors in the component-based engine performance models are replaced by variable geometry parameters to establish the TRM for a double bypass VCE.Obtained results are compared with the conventional component-based model from different aspects,including the accuracy and the convergence rate.The TRM is then employed to optimize the control schedule of a VCE.Obtained results show that the accuracy and the convergence rate of the proposed method are consistent with that of the conventional model.On the other hand,it is found that the new-model-optimized control schedules reduce the acceleration and deceleration time by 45%and 54%,respectively.Meanwhile,the surge margin of compressors,fuel–air ratio and the turbine inlet temperature maintained are within the acceptable criteria.It is concluded that the proposed TRM is a powerful method to design the acceleration and deceleration control schedule of the VCE.展开更多
This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational funct...This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational function solutions, Jacobi elliptic doubly periodic wave solutions, and nonsymmetrical kink solution are obtained. It is shown that the extended auxiliary equation method, with the help of a computer symbolic computation system, is reliable and effective in finding exact solutions of variable coefficient nonlinear evolution equations in mathematical physics.展开更多
In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squar...In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.展开更多
In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-f...In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.展开更多
By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact sol...By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact solutions including some new formal solutions are successfully picked up for the mKdV-sinh-Gordon equation by this approach.展开更多
In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved c...In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.展开更多
On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is present...On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.展开更多
Generalized functional separation of variables to nonlinear evolution equations is studied in terms of the extended group foliation method, which is based on the Lie point symmetry method. The approach is applied to n...Generalized functional separation of variables to nonlinear evolution equations is studied in terms of the extended group foliation method, which is based on the Lie point symmetry method. The approach is applied to nonlinear wave equations with variable speed and external force. A complete classification for the wave equation which admits functional separable solutions is presented. Some known results can be recovered by this approach.展开更多
基金supported by grants from the National Natural Science Foundation of China(No.U20B2006)the Guangdong Basic and Applied Basic Research Foundation(No.2022A1515110145)Young Elite Scientists Sponsorship Program by China Association for Science and Technology(No.2022QNRC001).
文摘The optimization of the waverider is constrained by the reversely designed leading edge and the constant shock angle distribution. This paper proposes a design method called the variable Leading-Edge Cone (vLEC) method to address these limitations. In the vLEC method, the waverider is directly designed from the preassigned leading edge and the variable shock angle distribution based on the Leading-Edge Cone (LEC) concept. Since the vLEC method is an approximate method, two test waveriders are designed and evaluated using numerical simulations to validate the shock design accuracy and the effectiveness of the vLEC method. The results show that the shocks of the test waveriders coincide well with the preassigned positions. Furthermore, four specifically designed application cases are conducted to analyze the performance benefits of the vLEC waveriders. The results of these cases indicate that, due to their variable shock angle distributions, the vLEC waveriders exhibit higher lift-to-drag ratios and better longitudinal static stability than conventional waveriders. Additionally, the vLEC waveriders demonstrate superior volumetric capacities near the symmetry plane, albeit with a minor decrease in volumetric efficiency.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
文摘The second-order backward differential formula(BDF2)and the scalar auxiliary variable(SAV)approach are applied to con‐struct the linearly energy stable numerical scheme with the variable time steps for the epitaxial thin film growth models.Under the stepratio condition 0<τ_(n)/τ_(n-1)<4.864,the modified energy dissipation law is proven at the discrete levels with regardless of time step size.Nu‐merical experiments are presented to demonstrate the accuracy and efficiency of the proposed numerical scheme.
文摘Mixed-variable problems are inevitable in engineering. However, few researches pay attention to discrete variables. This paper proposed a mixed-variable experimental design method (ODCD): first, the design variables were divided into discrete variables and continuous variables;then, the DVD method was employed for handling discrete variables, the LHD method was applied for continuous variables, and finally, a Columnwise-Pairwise Algorithm was used for the overall optimization of the design matrix. Experimental results demonstrated that the ODCD method outperforms in terms of the sample space coverage performance.
文摘Government credibility is an important asset of contemporary national governance, an important criterion for evaluating government legitimacy, and a key factor in measuring the effectiveness of government governance. In recent years, researchers’ research on government credibility has mostly focused on exploring theories and mechanisms, with little empirical research on this topic. This article intends to apply variable selection models in the field of social statistics to the issue of government credibility, in order to achieve empirical research on government credibility and explore its core influencing factors from a statistical perspective. Specifically, this article intends to use four regression-analysis-based methods and three random-forest-based methods to study the influencing factors of government credibility in various provinces in China, and compare the performance of these seven variable selection methods in different dimensions. The research results show that there are certain differences in simplicity, accuracy, and variable importance ranking among different variable selection methods, which present different importance in the study of government credibility issues. This study provides a methodological reference for variable selection models in the field of social science research, and also offers a multidimensional comparative perspective for analyzing the influencing factors of government credibility.
基金supported by the Graduate Student Scientific Research Innovation Project through Research Innovation Fund for Graduate Students in Shanxi Province(Project No.2024KY648).
文摘In this paper,we develop an advanced computational framework for the topology optimization of orthotropic materials using meshless methods.The approximation function is established based on the improved moving least squares(IMLS)method,which enhances the efficiency and stability of the numerical solution.The numerical solution formulas are derived using the improved element-free Galerkin(IEFG)method.We introduce the solid isotropic microstructures with penalization(SIMP)model to formulate a mathematical model for topology opti-mization,which effectively penalizes intermediate densities.The optimization problem is defined with the numerical solution formula and volume fraction as constraints.The objective function,which is the minimum value of flexibility,is optimized iteratively using the optimization criterion method to update the design variables efficiently and converge to an optimal solution.Sensitivity analysis is performed using the adjoint method,which provides accurate and efficient gradient information for the optimization algorithm.We validate the proposed framework through a series of numerical examples,including clamped beam,cantilever beam,and simply supported beam made of orthotropic materials.The convergence of the objective function is demonstrated by increasing the number of iterations.Additionally,the stability of the iterative process is analyzed by examining the fluctuation law of the volume fraction.By adjusting the parameters to an appropriate range,we achieve the final optimization results of the IEFG method without the checkerboard phenomenon.Comparative studies between the Element-Free Galerkin(EFG)and IEFG methods reveal that both methods yield consistent optimization results under identical parameter settings.However,the IEFG method significantly reduces computational time,highlighting its efficiency and suitability for orthotropic materials.
基金the project PID2022-139202OB-I00Neural Networks and Optimization Techniques for the Design and Safe Maintenance of Transportation Infrastructures:Volcanic Rock Geotechnics and Slope Stability(IA-Pyroslope),funded by the Spanish State Research Agency of the Ministry of Science,Innovation and Universities of Spain and the European Regional Development Fund,MCIN/AEI/10.13039/501100011033/FEDER,EU。
文摘Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.
文摘This study proposes a method for calculating the probability distribution of structural responses at different intensities using the endurance time(ET)method.The results can be used to calculate the fragility curve of structural collapse.The ET method involves dynamic analysis of a structure under an intensifying record over time.While conventional ET methods can determine the median of the structural response,they lack the ability to calculate its dispersion.To address this limitation,the present study utilizes ET analysis and single-degree-of-freedom(SDOF)systems to develop a method that considers the record-to-record variability for calculating the probability distribution of structural response.The accuracy of this method is evaluated by comparing it with the incremental dynamic analysis(IDA)method using special moment frames.The results demonstrate that the proposed method achieves a reasonably accurate estimation of dispersion while significantly reducing the computational burden(by approximately 95%)compared to the IDA method.
文摘In this paper, we study the solutions for variable-order time-fractional diffusion equations. A three-point combined compact difference (CCD) method is used to discretize the spatial variables to achieve sixth-order accuracy, while the exponential-sum-approximation (ESA) is used to approximate the variable-order Caputo fractional derivative in the temporal direction, and a novel spatial sixth-order hybrid ESA-CCD method is implemented successfully. Finally, the accuracy of the proposed method is verified by numerical experiments.
文摘This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)the Innovation Fund for Graduate Student of Shanghai University of China (Grant No.SHUCX120125)
文摘In this paper, based on the improved complex variable moving least-square （ICVMLS） approximation, a new complex variable meshless method （CVMM） for two-dimensional （2D） transient heat conduction problems is presented. The variational method is employed to obtain the discrete equations, and the essential boundary conditions are imposed by the penalty method. As the transient heat conduction problems are related to time, the Crank-Nicolson difference scheme for two-point boundary value problems is selected for the time discretization. Then the corresponding formulae of the CVMM for 2D heat conduction problems are obtained. In order to demonstrate the applicability of the proposed method, numerical examples are given to show the high convergence rate, good accuracy, and high efficiency of the CVMM presented in this paper.
基金supported by the National Natural Science Foundation of China (Grant No 10562002)the Natural Science Foundation of Inner Mongolia, China (Grant No 200508010103)+1 种基金the Specialized Research Fund for the Doctoral Program of Higher Education of China (Grant No 20070126002)the Inner Mongolia University Doctoral Scientific Research Starting Foundation
文摘The eigenvalue problem of an infinite-dimensional Hamiltonian operator appearing in the isotropic plane magnetoelectroelastic solids is studied. First, all the eigenvalues and their eigenfunctions in a rectangular domain are solved directly. Then the completeness of the eigenfunction system is proved, which offers a theoretic guarantee of the feasibility of variable separation method based on a Hamiltonian system for isotropic plane magnetoelectroelastic solids. Finally, the general solution for the equation in the rectangular domain is obtained by using the symplectic Fourier expansion method.
基金supported by the Aviation Power Foundation of China(6141B09050382)。
文摘Studies show that different geometries of a Variable Cycle Engine(VCE)can be adjusted during the transient stage of the engine operation to improve the engine performance.However,this improvement increases the complexity of the acceleration and deceleration control schedule.In order to resolve this problem,the Transient-state Reverse Method(TRM)is established in the present study based on the Steady-state Reverse Method(SRM)and the Virtual Power Extraction Method(VPEM).The state factors in the component-based engine performance models are replaced by variable geometry parameters to establish the TRM for a double bypass VCE.Obtained results are compared with the conventional component-based model from different aspects,including the accuracy and the convergence rate.The TRM is then employed to optimize the control schedule of a VCE.Obtained results show that the accuracy and the convergence rate of the proposed method are consistent with that of the conventional model.On the other hand,it is found that the new-model-optimized control schedules reduce the acceleration and deceleration time by 45%and 54%,respectively.Meanwhile,the surge margin of compressors,fuel–air ratio and the turbine inlet temperature maintained are within the acceptable criteria.It is concluded that the proposed TRM is a powerful method to design the acceleration and deceleration control schedule of the VCE.
基金Project supported by the Fundamental Research Funds for the Central Universities （Grant No. 2010B17914） and the National Natural Science Foundation of China （Grant No. 10926162）.
文摘This paper applies an extended auxiliary equation method to obtain exact solutions of the KdV equation with variable coefficients. As a result, solitary wave solutions, trigonometric function solutions, rational function solutions, Jacobi elliptic doubly periodic wave solutions, and nonsymmetrical kink solution are obtained. It is shown that the extended auxiliary equation method, with the help of a computer symbolic computation system, is reliable and effective in finding exact solutions of variable coefficient nonlinear evolution equations in mathematical physics.
基金supported by the National Natural Science Foundation of China (Grant No.11026223)the Shanghai Leading Academic Discipline Project,China (Grant No.S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No.SHUCX112359)
文摘In this paper, the improved complex variable moving least-squares （ICVMLS） approximation is presented. The ICVMLS approximation has an explicit physics meaning. Compared with the complex variable moving least-squares （CVMLS） approximations presented by Cheng and Ren, the ICVMLS approximation has a great computational precision and efficiency. Based on the element-free Galerkin （EFG） method and the ICVMLS approximation, the improved complex variable element-free Galerkin （ICVEFG） method is presented for two-dimensional elasticity problems, and the corresponding formulae are obtained. Compared with the conventional EFC method, the ICVEFG method has a great computational accuracy and efficiency. For the purpose of demonstration, three selected numerical examples are solved using the ICVEFG method.
基金Project supported by the National Natural Science Foundation of China (Grant No. 11171208)the Shanghai Leading Academic Discipline Project, China (Grant No. S30106)the Innovation Fund Project for Graduate Student of Shanghai University,China (Grant No. SHUCX112359)
文摘In this paper, based on the element-free Galerkin （EFG） method and the improved complex variable moving least- square （ICVMLS） approximation, a new meshless method, which is the improved complex variable element-free Galerkin （ICVEFG） method for two-dimensional potential problems, is presented. In the method, the integral weak form of control equations is employed, and the Lagrange multiplier is used to apply the essential boundary conditions. Then the corresponding formulas of the ICVEFG method for two-dimensional potential problems are obtained. Compared with the complex variable moving least-square （CVMLS） approximation proposed by Cheng, the functional in the ICVMLS approximation has an explicit physical meaning. Furthermore, the ICVEFG method has greater computational precision and efficiency. Three numerical examples are given to show the validity of the proposed method.
基金Project supported by the National Natural Science Foundation of China (Grant No 10672053)
文摘By introducing a more general auxiliary ordinary differential equation （ODE）, a modified variable separated ODE method is developed for solving the mKdV-sinh-Gordon equation. As a result, many explicit and exact solutions including some new formal solutions are successfully picked up for the mKdV-sinh-Gordon equation by this approach.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11171208 and U1433104)
文摘In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.
基金supported by the National Natural Science Foundation of China (Grant No.10871124)the Innovation Program of Shanghai Municipal Education Commission,China (Grant No.09ZZ99)
文摘On the basis of the reproducing kernel particle method （RKPM）, a new meshless method, which is called the complex variable reproducing kernel particle method （CVRKPM）, for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM.
文摘Generalized functional separation of variables to nonlinear evolution equations is studied in terms of the extended group foliation method, which is based on the Lie point symmetry method. The approach is applied to nonlinear wave equations with variable speed and external force. A complete classification for the wave equation which admits functional separable solutions is presented. Some known results can be recovered by this approach.
