This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the...This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.展开更多
Unconditionally stable higher-order accurate time step integration algorithms based on the differential quadrature method (DQM) for second-order initial value problems were applied and the quadrature rules of DQM, com...Unconditionally stable higher-order accurate time step integration algorithms based on the differential quadrature method (DQM) for second-order initial value problems were applied and the quadrature rules of DQM, computing of the weighting coefficients and choices of sampling grid points were discussed. Some numerical examples dealing with the heat transfer problem, the second-order differential equation of imposed vibration of linear single-degree-of-freedom systems and double-degree-of-freedom systems, the nonlinear move differential equation and a beam forced by a changing load were computed, respectively. The results indicated that the algorithm can produce highly accurate solutions with minimal time consumption, and that the system total energy can remain conservative in the numerical computation.展开更多
The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elas...The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.展开更多
This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredepend...This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.展开更多
The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the...The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.展开更多
Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a t...Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.展开更多
An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise in...An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the展开更多
This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we stu...This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.展开更多
This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe...This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe conveying fluid.Correspondingly,the natural frequency of axially FGM pipes conveying fluid is calculated using the differential quadrature method(DQM).A variable sensitivity analysis(VSA)is introduced to measure the effect of each random variable,and a mode sensitivity analysis(MSA)is introduced to acquire the importance ranking of failure modes.Then,an active learning Kriging(ALK)method is established to calculate the resonance failure probability and sensitivity indices,which greatly improves the application of resonance reliability analysis for pipelines in engineering practice.Based on the resonance reliability analysis method,the effects of fluid velocity,volume fraction and fluid density of axially FGM pipe conveying fluid on resonance reliability are analyzed.The results demonstrate that the proposed method has great performance in the anti-resonance analysis of pipes conveying fluid.展开更多
In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.T...In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.The reduced differential transforms method(RDTM)is one of the interesting methods for finding the approximate solutions for nonlinear problems.We apply the RDTM to discuss the analytic approximate solutions to the SI1I2R model for the spread of virus HCV-subtype and SIR childhood disease model.We discuss the numerical results at some special values of parameters in the approximate solutions.We use the computer software package such as Mathematical to find more iteration when calculating the approximate solutions.Graphical results and discussed quantitatively are presented to illustrate behavior of the obtained approximate solutions.展开更多
Based on the“three box”exergy analysis model,a black box-gray box hierarchical exergy analysis and evaluation method is put forward in this paper,which is applied to evaluate the power generation technology of diffe...Based on the“three box”exergy analysis model,a black box-gray box hierarchical exergy analysis and evaluation method is put forward in this paper,which is applied to evaluate the power generation technology of differential pressure produced by natural gas expansion.By using the exergy analysis theory,the black box-gray box hierarchical exergy analysis models of three differential pressure power generation technologies are established respectively.Firstly,the“black box”analysis models of main energy consuming equipment are established,and then the“gray box”analysis model of the total system is established.Based on the calculation results of exergy analysis indexes,the weak energy consumption equipment in the whole power generation process is accurately located.Taking a gas field in southwest China as an example,the comprehensive energy consumption evaluation of the three power generation technologies is carried out,and the technology with the best energy consumption condition among the three technologies is determined.Finally,the rationalization improvement measures are put forward from improving the air tightness,replacing the deflector and reducing the flow loss.展开更多
We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on t...We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method （HAM）. This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.展开更多
In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical s...In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical solution of the model in fractional order. The obtained results are compaired with those obtained by forth order Runge-Kutta method and nonstandard numerical method in the integer case. Finally, we present some numerical results.展开更多
In this paper, the Homotopy Analysis Method (HAM) has been used to solve an eco-epidemic model equation. The algorithm of approximate analytical solution is obtained. HAM contains the auxiliary parameterhwhich provide...In this paper, the Homotopy Analysis Method (HAM) has been used to solve an eco-epidemic model equation. The algorithm of approximate analytical solution is obtained. HAM contains the auxiliary parameterhwhich provides us with a convenient way to adjust and control convergence region and rate of solution series. The results obtained show that these algorithms are accurate and efficient for the model.展开更多
In this paper, we propose a unified differential operator method to study mechanical vibrations, solving inhomogeneous linear ordinary differential equations with constant coefficients. The main advantage of this new ...In this paper, we propose a unified differential operator method to study mechanical vibrations, solving inhomogeneous linear ordinary differential equations with constant coefficients. The main advantage of this new method is that the differential operator D in the numerator of the fraction has no effect on input functions (i.e., the derivative operation is removed) because we take the fraction as a whole part in the partial fraction expansion. The method in various variants is widely implemented in related fields in mechanics and engineering. We also point out that the same mistakes in the differential operator method are found in the related references [1-4].展开更多
In this study,the impacts of internal heat generation on heat transfer enhancement of porous fin is theoretical investigated using differential transform method.The parametric studies reveal that porosity enhances the...In this study,the impacts of internal heat generation on heat transfer enhancement of porous fin is theoretical investigated using differential transform method.The parametric studies reveal that porosity enhances the fin heat dissipating capacity but the internal heat generation decreases the heat enhancement capacity of extended surface.Also,it is established that when the internal heat parameter increases to some certain values,some negative effects are recorded where the fin stores heat rather than dissipating it.This scenario defeats the prime purpose of the cooling fin.Additionally,it is established in the present study that the limiting value of porosity parameter for thermal stability for the passive device increases as internal heat parameter increases.This shows that although the internal heat parameter can help assist higher range and value of thermal stability of the fin,it produces negative effect which greatly defeats the ultimate purpose of the fin.The results in the work will help in fin design for industrial applications where internal heat generation is involved.展开更多
The main aim of this paper is to present and emphasize the contribution of stochastic numerical methods as must tools for the modern econometric modelisation. Indeed, the stochastic numerical methods play an important...The main aim of this paper is to present and emphasize the contribution of stochastic numerical methods as must tools for the modern econometric modelisation. Indeed, the stochastic numerical methods play an important role in mathematical modelling and the econometric analysis because they model uncertainties that govern the real-world data. However these powerful tools are not well-known and understood by many economists and financial econometricians.展开更多
Nonlinear partial differential equations(PDEs)are crucial for modeling complex fluid dynamics and are foundational to many computational fluid dynamics(CFD)applications.However,solving these nonlinear PDEs is challeng...Nonlinear partial differential equations(PDEs)are crucial for modeling complex fluid dynamics and are foundational to many computational fluid dynamics(CFD)applications.However,solving these nonlinear PDEs is challenging due to the vast computational resources they demand,highlighting the pressing need for more efficient computational methods.Quantum computing offers a promising but technically challenging approach to solving nonlinear PDEs.Recently,Liao[arXiv:2406.15821]proposed a framework that leverages quantum computing to accelerate the solution of nonlinear PDEs based on the homotopy analysis method(HAM),a semi-analytical technique that transforms nonlinear PDEs into a series of linear PDEs.However,the no-cloning theorem in quantum computing poses a major limitation,where directly applying quantum simulation to each HAM step results in exponential complexity growth with the HAM truncation order.This study introduces a“quantum-compatible linearization”approach that maps the whole HAM process into a system of linear PDEs,allowing for a one-time solution using established quantum PDE solvers.Our method preserves the exponential speedup of quantum linear PDE solvers while ensuring that computational complexity increases only polynomially with the HAM truncation order.We demonstrate the efficacy of our approach by applying it to the Burgers'equation and the Korteweg-de Vries(KdV)equation.Our approach provides a novel pathway for transforming nonlinear PDEs into linear PDEs,with potential applications to fluid dynamics.This work thus lays the foundation for developing quantum algorithms capable of solving the Navier-Stokes equations,ultimately offering a promising route to accelerate their solutions using quantum computing.展开更多
The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant i...The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method（DDM） has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction（SFSI） models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method（FFD）.The relative importance（RI） of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM（e.g.,accurate,efficient,insensitive to numerical noise）.Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.展开更多
The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody sys...The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody systems. However, the formulation has so many principles in choosing the generalized coordinates that it hinders the implementation of modeling automation, A first order direct sensitivity analysis approach to multibody systems formulated with novel natural coordinates is presented. Firstly, a new selection method for natural coordinate is developed. The method introduces 12 coordinates to describe the position and orientation of a spatial object. On the basis of the proposed natural coordinates, rigid constraint conditions, the basic constraint elements as well as the initial conditions for the governing equations are derived. Considering the characteristics of the governing equations, the newly proposed generalized-ct integration method is used and the corresponding algorithm flowchart is discussed. The objective function, the detailed analysis process of first order direct sensitivity analysis and related solving strategy are provided based on the previous modeling system Finally, in order to verify the validity and accuracy of the method presented, the sensitivity analysis of a planar spinner-slider mechanism and a spatial crank-slider mechanism are conducted. The test results agree well with that of the finite difference method, and the maximum absolute deviation of the results is less than 3%. The proposed approach is not only convenient for automatic modeling, but also helpful for the reduction of the complexity of sensitivity analysis, which provides a practical and effective way to obtain sensitivity for the optimization problems of multibody systems.展开更多
文摘This paper presents a combined application of the finite element method （FEM） and the differential quadrature method （DQM） to vibration and buckling problems of rectangular plates. The proposed scheme combines the geometry flexibility of the FEM and the high accuracy and efficiency of the DQM. The accuracy of the present method is demonstrated by comparing the obtained results with those available in the literature. It is shown that highly accurate results can be obtained by using a small number of finite elements and DQM sample points. The proposed method is suitable for the problems considered due to its simplicity and potential for further development.
文摘Unconditionally stable higher-order accurate time step integration algorithms based on the differential quadrature method (DQM) for second-order initial value problems were applied and the quadrature rules of DQM, computing of the weighting coefficients and choices of sampling grid points were discussed. Some numerical examples dealing with the heat transfer problem, the second-order differential equation of imposed vibration of linear single-degree-of-freedom systems and double-degree-of-freedom systems, the nonlinear move differential equation and a beam forced by a changing load were computed, respectively. The results indicated that the algorithm can produce highly accurate solutions with minimal time consumption, and that the system total energy can remain conservative in the numerical computation.
基金supported by the Iranian Nanotechnology Development Committee and the University of Kashan(No.363452/10)
文摘The size-dependent effect on the biaxial and shear nonlinear buckling analysis of an isotropic and orthotropic micro-plate based on the surface stress, the modified couple stress theory （MCST）, and the nonlocal elasticity theories using the differential quadrature method （DQM） is presented. Main advantages of the MCST over the classical theory （CT） are the inclusion of the asymmetric couple stress tensor and the consideration of only one material length scale parameter. Based on the nonlinear von Karman assumption, the governing equations of equilibrium for the micro-classical plate consid- ering midplane displacements are derived based on the minimum principle of potential energy. Using the DQM, the biaxial and shear critical buckling loads of the micro-plate for various boundary conditions are obtained. Accuracy of the obtained results is validated by comparing the solutions with those reported in the literature. A parametric study is conducted to show the effects of the aspect ratio, the side-to-thickness ratio, Eringen＇s nonlocal parameter, the material length scale parameter, Young＇s modulus of the surface layer, the surface residual stress, the polymer matrix coefficients, and various boundary conditions on the dimensionless uniaxial, biaxial, and shear critical buckling loads. The results indicate that the critical buckling loads are strongly sensitive to Eringen＇s nonlocal parameter, the material length scale parameter, and the surface residual stress effects, while the effect of Young＇s modulus of the surface layer on the critical buckling load is negligible. Also, considering the size dependent effect causes the increase in the stiffness of the orthotropic micro-plate. The results show that the critical biaxial buckling load increases with an increase in G12/E2 and vice versa for E1/E2. It is shown that the nonlinear biaxial buckling ratio decreases as the aspect ratio increases and vice versa for the buckling amplitude. Because of the most lightweight micro-composite materials with high strength/weight and stiffness/weight ratios, it is anticipated that the results of the present work are useful in experimental characterization of the mechanical properties of micro-composite plates in the aircraft industry and other engineering applications.
基金supported by a grant from the National Science and Technology Council of the Republic of China(Grant Number:MOST 112-2221-E-006-048-MY2).
文摘This work develops a Hermitian C^(2) differential reproducing kernel interpolation meshless(DRKIM)method within the consistent couple stress theory(CCST)framework to study the three-dimensional(3D)microstructuredependent static flexural behavior of a functionally graded(FG)microplate subjected to mechanical loads and placed under full simple supports.In the formulation,we select the transverse stress and displacement components and their first-and second-order derivatives as primary variables.Then,we set up the differential reproducing conditions(DRCs)to obtain the shape functions of the Hermitian C^(2) differential reproducing kernel(DRK)interpolant’s derivatives without using direct differentiation.The interpolant’s shape function is combined with a primitive function that possesses Kronecker delta properties and an enrichment function that constituents DRCs.As a result,the primary variables and their first-and second-order derivatives satisfy the nodal interpolation properties.Subsequently,incorporating ourHermitianC^(2)DRKinterpolant intothe strong formof the3DCCST,we develop a DRKIM method to analyze the FG microplate’s 3D microstructure-dependent static flexural behavior.The Hermitian C^(2) DRKIM method is confirmed to be accurate and fast in its convergence rate by comparing the solutions it produces with the relevant 3D solutions available in the literature.Finally,the impact of essential factors on the transverse stresses,in-plane stresses,displacements,and couple stresses that are induced in the loaded microplate is examined.These factors include the length-to-thickness ratio,the material length-scale parameter,and the inhomogeneity index,which appear to be significant.
文摘The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.
基金co-supported by the National Natural Science Foundation of China(No.12072365)the Technology Innovation Team of Manned Space Engineering,China。
文摘Aimed at the demand of contingency return at any time during the near-moon phase in the manned lunar landing missions,a fast calculation method for three-impulse contingency return trajectories is proposed.Firstly,a three-impulse contingency return trajectory scheme is presented by combining the Lambert transfer and maneuver at the special point.Secondly,a calculation model of three-impulse contingency return trajectories is established.Then,fast calculation methods are proposed by adopting the high-order Taylor expansion of differential algebra in the twobody trajectory dynamics model and perturbed trajectory dynamics model.Finally,the performance of the proposed methods is verified by numerical simulation.The results indicate that the fast calculation method of two-body trajectory has higher calculation efficiency compared to the semi-analytical calculation method under a certain accuracy condition.Due to its high efficiency,the characteristics of the three-impulse contingency return trajectories under different contingency scenarios are further analyzed expeditiously.These findings can be used for the design of contingency return trajectories in future manned lunar landing missions.
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
文摘This study was suggested by previous work on the simulation of evolution equations with scale-dependent processes,e.g.,wave-propagation or heat-transfer,that are modeled by wave equations or heat equations.Here,we study both parabolic and hyperbolic equations.We focus on ADI (alternating direction implicit) methods and LOD (locally one-dimensional) methods,which are standard splitting methods of lower order,e.g.second-order.Our aim is to develop higher-order ADI methods,which are performed by Richardson extrapolation,Crank-Nicolson methods and higher-order LOD methods,based on locally higher-order methods.We discuss the new theoretical results of the stability and consistency of the ADI methods.The main idea is to apply a higher- order time discretization and combine it with the ADI methods.We also discuss the dis- cretization and splitting methods for first-order and second-order evolution equations. The stability analysis is given for the ADI method for first-order time derivatives and for the LOD (locally one-dimensional) methods for second-order time derivatives.The higher-order methods are unconditionally stable.Some numerical experiments verify our results.
基金The funding was provided by Laboratory Fund (Grant No.SYJJ200320).
文摘This paper aims to solve the resonance failure probability and develop an effective method to estimate the effects of variables and failure modes on failure probability of axially functionally graded material(FGM)pipe conveying fluid.Correspondingly,the natural frequency of axially FGM pipes conveying fluid is calculated using the differential quadrature method(DQM).A variable sensitivity analysis(VSA)is introduced to measure the effect of each random variable,and a mode sensitivity analysis(MSA)is introduced to acquire the importance ranking of failure modes.Then,an active learning Kriging(ALK)method is established to calculate the resonance failure probability and sensitivity indices,which greatly improves the application of resonance reliability analysis for pipelines in engineering practice.Based on the resonance reliability analysis method,the effects of fluid velocity,volume fraction and fluid density of axially FGM pipe conveying fluid on resonance reliability are analyzed.The results demonstrate that the proposed method has great performance in the anti-resonance analysis of pipes conveying fluid.
文摘In this paper,we study the approximate solutions for some of nonlinear Biomathematics models via the e-epidemic SI1I2R model characterizing the spread of viruses in a computer network and SIR childhood disease model.The reduced differential transforms method(RDTM)is one of the interesting methods for finding the approximate solutions for nonlinear problems.We apply the RDTM to discuss the analytic approximate solutions to the SI1I2R model for the spread of virus HCV-subtype and SIR childhood disease model.We discuss the numerical results at some special values of parameters in the approximate solutions.We use the computer software package such as Mathematical to find more iteration when calculating the approximate solutions.Graphical results and discussed quantitatively are presented to illustrate behavior of the obtained approximate solutions.
基金financially supported by the National Natural Science Foundation of China(52074089 and 51534004)Natural Science Foundation of Heilongjiang Province of China(LH2019E019)。
文摘Based on the“three box”exergy analysis model,a black box-gray box hierarchical exergy analysis and evaluation method is put forward in this paper,which is applied to evaluate the power generation technology of differential pressure produced by natural gas expansion.By using the exergy analysis theory,the black box-gray box hierarchical exergy analysis models of three differential pressure power generation technologies are established respectively.Firstly,the“black box”analysis models of main energy consuming equipment are established,and then the“gray box”analysis model of the total system is established.Based on the calculation results of exergy analysis indexes,the weak energy consumption equipment in the whole power generation process is accurately located.Taking a gas field in southwest China as an example,the comprehensive energy consumption evaluation of the three power generation technologies is carried out,and the technology with the best energy consumption condition among the three technologies is determined.Finally,the rationalization improvement measures are put forward from improving the air tightness,replacing the deflector and reducing the flow loss.
文摘We present a new reliable analytical study for solving the discontinued problems arising in nanotechnology. Such problems are presented as nonlinear differential-difference equations. The proposed method is based on the Laplace trans- form with the homotopy analysis method （HAM）. This method is a powerful tool for solving a large amount of problems. This technique provides a series of functions which may converge to the exact solution of the problem. A good agreement between the obtained solution and some well-known results is obtained.
文摘In this paper, we consider the giving up smoking model. First, we present the giving up smoking model in fractional order. Then the homotopy analysis method (HAM) is employed to compute an approximate and analytical solution of the model in fractional order. The obtained results are compaired with those obtained by forth order Runge-Kutta method and nonstandard numerical method in the integer case. Finally, we present some numerical results.
文摘In this paper, the Homotopy Analysis Method (HAM) has been used to solve an eco-epidemic model equation. The algorithm of approximate analytical solution is obtained. HAM contains the auxiliary parameterhwhich provides us with a convenient way to adjust and control convergence region and rate of solution series. The results obtained show that these algorithms are accurate and efficient for the model.
文摘In this paper, we propose a unified differential operator method to study mechanical vibrations, solving inhomogeneous linear ordinary differential equations with constant coefficients. The main advantage of this new method is that the differential operator D in the numerator of the fraction has no effect on input functions (i.e., the derivative operation is removed) because we take the fraction as a whole part in the partial fraction expansion. The method in various variants is widely implemented in related fields in mechanics and engineering. We also point out that the same mistakes in the differential operator method are found in the related references [1-4].
文摘In this study,the impacts of internal heat generation on heat transfer enhancement of porous fin is theoretical investigated using differential transform method.The parametric studies reveal that porosity enhances the fin heat dissipating capacity but the internal heat generation decreases the heat enhancement capacity of extended surface.Also,it is established that when the internal heat parameter increases to some certain values,some negative effects are recorded where the fin stores heat rather than dissipating it.This scenario defeats the prime purpose of the cooling fin.Additionally,it is established in the present study that the limiting value of porosity parameter for thermal stability for the passive device increases as internal heat parameter increases.This shows that although the internal heat parameter can help assist higher range and value of thermal stability of the fin,it produces negative effect which greatly defeats the ultimate purpose of the fin.The results in the work will help in fin design for industrial applications where internal heat generation is involved.
文摘The main aim of this paper is to present and emphasize the contribution of stochastic numerical methods as must tools for the modern econometric modelisation. Indeed, the stochastic numerical methods play an important role in mathematical modelling and the econometric analysis because they model uncertainties that govern the real-world data. However these powerful tools are not well-known and understood by many economists and financial econometricians.
基金supported by the National Key Research and Development Program of China(Grant No.2023YFB4502500)the National Natural Science Foundation of China(Grant No.12404564)Anhui Province Science and Technology Innovation(Grant No.202423s06050001)。
文摘Nonlinear partial differential equations(PDEs)are crucial for modeling complex fluid dynamics and are foundational to many computational fluid dynamics(CFD)applications.However,solving these nonlinear PDEs is challenging due to the vast computational resources they demand,highlighting the pressing need for more efficient computational methods.Quantum computing offers a promising but technically challenging approach to solving nonlinear PDEs.Recently,Liao[arXiv:2406.15821]proposed a framework that leverages quantum computing to accelerate the solution of nonlinear PDEs based on the homotopy analysis method(HAM),a semi-analytical technique that transforms nonlinear PDEs into a series of linear PDEs.However,the no-cloning theorem in quantum computing poses a major limitation,where directly applying quantum simulation to each HAM step results in exponential complexity growth with the HAM truncation order.This study introduces a“quantum-compatible linearization”approach that maps the whole HAM process into a system of linear PDEs,allowing for a one-time solution using established quantum PDE solvers.Our method preserves the exponential speedup of quantum linear PDE solvers while ensuring that computational complexity increases only polynomially with the HAM truncation order.We demonstrate the efficacy of our approach by applying it to the Burgers'equation and the Korteweg-de Vries(KdV)equation.Our approach provides a novel pathway for transforming nonlinear PDEs into linear PDEs,with potential applications to fluid dynamics.This work thus lays the foundation for developing quantum algorithms capable of solving the Navier-Stokes equations,ultimately offering a promising route to accelerate their solutions using quantum computing.
基金National Key Research and Development Program of China under Grant No.2016YFC0701106Natural Sciences and Engineering Research Council of Canada via Discovery under Grant No.NSERC RGPIN-2017-05556 Li
文摘The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method（DDM） has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction（SFSI） models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method（FFD）.The relative importance（RI） of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM（e.g.,accurate,efficient,insensitive to numerical noise）.Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.
基金supported by National Defense Pre-research Foundation of China during the 12th Five-Year Plan Period(Grant No.51036050107)
文摘The classical natural coordinate modeling method which removes the Euler angles and Euler parameters from the governing equations is particularly suitable for the sensitivity analysis and optimization of multibody systems. However, the formulation has so many principles in choosing the generalized coordinates that it hinders the implementation of modeling automation, A first order direct sensitivity analysis approach to multibody systems formulated with novel natural coordinates is presented. Firstly, a new selection method for natural coordinate is developed. The method introduces 12 coordinates to describe the position and orientation of a spatial object. On the basis of the proposed natural coordinates, rigid constraint conditions, the basic constraint elements as well as the initial conditions for the governing equations are derived. Considering the characteristics of the governing equations, the newly proposed generalized-ct integration method is used and the corresponding algorithm flowchart is discussed. The objective function, the detailed analysis process of first order direct sensitivity analysis and related solving strategy are provided based on the previous modeling system Finally, in order to verify the validity and accuracy of the method presented, the sensitivity analysis of a planar spinner-slider mechanism and a spatial crank-slider mechanism are conducted. The test results agree well with that of the finite difference method, and the maximum absolute deviation of the results is less than 3%. The proposed approach is not only convenient for automatic modeling, but also helpful for the reduction of the complexity of sensitivity analysis, which provides a practical and effective way to obtain sensitivity for the optimization problems of multibody systems.
