Proper matching of forestry machinery is important when raising mechanization levels for forestry production. In the matching process, forestry machinery needs not only expertise, but also improved methods for solving...Proper matching of forestry machinery is important when raising mechanization levels for forestry production. In the matching process, forestry machinery needs not only expertise, but also improved methods for solving problems. I propose combination of case-based reasoning (CBR) and rule-based reasoning (RBR) by calculating the similarity of quantitative parameters of various forestry machines in an analytical and hierarchical process. I calculated the similarity of machin-ery used in forest industries to enable better selection and matching of equipment. I propose a weight-value adjusting method based on sums of squares of deviations in which the individual parameter weights were modified in the process of application. During the process of system design, I put forward a design method knowledge base and generated a dynamic web reasoning framework to integrate the processes of forest industry machinery selection and weight-value adjustment. This enables expansion of the scope of the complete system and enhancement of the reasoning efficiency. I demonstrate the validity and practicability of this method using a practical example.展开更多
This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solve...This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.展开更多
For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geomet...For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.展开更多
This paper presents a new concept called Unified and Integrated Method for a shear deformable beam element.In this method,Timoshenko beam theory is unified and integrated in such a way that takes into account the effe...This paper presents a new concept called Unified and Integrated Method for a shear deformable beam element.In this method,Timoshenko beam theory is unified and integrated in such a way that takes into account the effect of transverse shear and maintains the shear locking free condition at the same time to generate proper behavior in the analysis of thin to thick beams.The unified and integrated method is applied to finite element analysis(FEA)and isogeometric analysis(IGA)on two-node beam element.This method will be used to analyze uniformly loaded beams with various boundary conditions.A shear influence factor of f,which is a function of beam thickness ratio(L/h),is expressed explicitly as control of the transverse shear strain effect.The analysis gives interesting results showing that applying the unified and integrated method in FEA and IGA will yield exact values of DOF’s and displacement function even when using only a single element.Numerical examples demonstrate the validity and efficiency of the unified and integrated methods.展开更多
The morphing technology of hypersonic vehicle improved the flight performance by changing aerodynamic characteristics with shape deformations,but the design of guidance and control system with morphing laws remained t...The morphing technology of hypersonic vehicle improved the flight performance by changing aerodynamic characteristics with shape deformations,but the design of guidance and control system with morphing laws remained to be explored.An Integrated of Guidance,Control and Morphing(IGCM)method for Hypersonic Morphing Vehicle(HMV)was developed in this paper.The IGCM method contributed to an effective solution of morphing characteristic to improve flight performance and reject the disturbance for guidance and control system caused by the morphing system for HMV in gliding phase.The IGCM models were established based on the motion models and aerodynamic models of the variable span vehicle.Then the IGCM method was designed by adaptive block dynamic surface back-stepping method with stability proof.The parallel controlled simulations’results showed the effectiveness in accomplishing the flight mission of IGCM method in glide phase with smaller terminal errors.The velocity loss of HMV was reduced by 32.8%which inferred less flight time and larger terminal flight velocity than invariable span vehicle.Under the condition of large deviations of aerodynamic parameters and atmospheric density,the robustness of IGCM method with variable span was verified.展开更多
In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley ...In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley integrated-empowerment benefit-distribution method.First,through literature survey and expert interview to identify the risk factors at various stages of the project,a dynamic risk-factor indicator system is developed.Second,to obtain a more meaningful risk-calculation result,the subjective and objective weights are combined,the weights of the risk factors at each stage are determined by the expert scoring method and entropy weight method,and the interest distribution model based on multi-dimensional risk factors is established.Finally,an example is used to verify the rationality of the method for the benefit distribution of the charging-pile project.The results of the example indicate that the limitations of the Shapley method can be reasonably avoided,and the applicability of the model for the benefit distribution of the charging-pile project is verified.展开更多
A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains.To realize this method,we design a new wavelet basis function,by which we construct a fifth...A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains.To realize this method,we design a new wavelet basis function,by which we construct a fifth-order numerical scheme for the approximation of multi-dimensional functions and their multiple integrals defined in complex domains.In the solution of differential equations,various derivatives of the unknown function are denoted as new functions.Then,the integral relations between these functions are applied in terms of wavelet approximation of multiple integrals.Therefore,the original equation with derivatives of various orders can be converted to a system of algebraic equations with discrete nodal values of the highest-order derivative.During the application of the proposed method,boundary conditions can be automatically included in the integration operations,and relevant matrices can be assured to exhibit perfect sparse patterns.As examples,we consider several second-order mathematics problems defined on regular and irregular domains and the fourth-order bending problems of plates with various shapes.By comparing the solutions obtained by the proposed method with the exact solutions,the new multiresolution method is found to have a convergence rate of fifth order.The solution accuracy of this method with only a few hundreds of nodes can be much higher than that of the finite element method(FEM)with tens of thousands of elements.In addition,because the accuracy order for direct approximation of a function using the proposed basis function is also fifth order,we may conclude that the accuracy of the proposed method is almost independent of the equation order and domain complexity.展开更多
Topology and performance are the two main topics dealt in the development of robotic mechanisms.However,it is still a challenge to connect them by integrating the modeling and design process of both parts in a unified...Topology and performance are the two main topics dealt in the development of robotic mechanisms.However,it is still a challenge to connect them by integrating the modeling and design process of both parts in a unified frame.As the properties associated with topology and performance,finite motion and instantaneous motion of the robot play key roles in the procedure.On the purpose of providing a fundamental preparation for integrated modeling and design,this paper carries out a review on the existing unified mathematic frameworks for motion description and computation,involving matrix Lie group and Lie algebra,dual quaternion and pure dual quaternion,finite screw and instantaneous screw.Besides the application in robotics,the review of the work from these mathematicians concentrates on the description,composition and intersection operations of the finite and instantaneous motions,especially on the exponential-differential maps which connect the two sides.Furthermore,an in-depth discussion is worked out by investigating the algebraical relationship among these methods and their further progress in integrated robotic development.The presented review offers insightful investigation to the motion description and computation,and therefore would help designers to choose appropriate mathematical tool in the integrated design and modeling and design of mechanisms and robots.展开更多
The objective is to develop a model considering demand dependent on selling price and deterioration occurs after a certain period of time, which follows two-parameter Weibull distribution. Shortages are allowed and fu...The objective is to develop a model considering demand dependent on selling price and deterioration occurs after a certain period of time, which follows two-parameter Weibull distribution. Shortages are allowed and fully backlogged. Fuzzy optimal solution is obtained by considering hexagonal fuzzy numbers and for defuzzification Graded Mean Integration Representation Method. A numerical example is provided for the illustration of crisp and fuzzy, both models. To observe the effect of changes in parameters, sensitivity analysis is carried out.展开更多
Infertility is one of the difficult complicated diseases. Many couples suffer from it. The pathogenesis is very complicated. The imbalance or lesson of any link of the reproductive system can cause infertility. This p...Infertility is one of the difficult complicated diseases. Many couples suffer from it. The pathogenesis is very complicated. The imbalance or lesson of any link of the reproductive system can cause infertility. This paper summarizes the treatment of female infertility by integrated Traditional Chinese Medicine (TCM) and Western Medicine (WM) which can not only improve the ovulation rate and pregnancy rate, but also decrease the complications. The effects are better than that by TCM or WM only. Therefore, the coupling method is worth to be used widely in clinical practice.展开更多
The singularities, oscillatory performances and the contributing factors to the 3-'D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field ...The singularities, oscillatory performances and the contributing factors to the 3-'D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.展开更多
In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then th...In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.展开更多
A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not ...This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration (iterative) and the operator-splitting (non-iterative) is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom (SDOF) cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1:2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1:2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely use...As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.展开更多
The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcom...The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.展开更多
In this paper,the field method for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass.Two examples are given to illustrat...In this paper,the field method for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass.Two examples are given to illustrate its application.展开更多
This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structu...This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.展开更多
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展开更多
基金financially supported by the Fundamental Research Funds for the Central Universities Nos.DL12EB01-03the planning subject of "the Twelfth Five-Year-Plan" in National Science and Technology Nos.2012AA102003-2Heilongjiang Natural Science Fund in China Nos.F201116
文摘Proper matching of forestry machinery is important when raising mechanization levels for forestry production. In the matching process, forestry machinery needs not only expertise, but also improved methods for solving problems. I propose combination of case-based reasoning (CBR) and rule-based reasoning (RBR) by calculating the similarity of quantitative parameters of various forestry machines in an analytical and hierarchical process. I calculated the similarity of machin-ery used in forest industries to enable better selection and matching of equipment. I propose a weight-value adjusting method based on sums of squares of deviations in which the individual parameter weights were modified in the process of application. During the process of system design, I put forward a design method knowledge base and generated a dynamic web reasoning framework to integrate the processes of forest industry machinery selection and weight-value adjustment. This enables expansion of the scope of the complete system and enhancement of the reasoning efficiency. I demonstrate the validity and practicability of this method using a practical example.
基金the National Science and Tech-nology Council,Taiwan for their financial support(Grant Number NSTC 111-2221-E-019-048).
文摘This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.
文摘For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.
基金support from the Ministry of Research and Technology/National Research and Inovation Agency(RISTEK-BRIN),Indonesia,through the PDUPT program(Grant No.NKB-1641/UN2.R3.1/HKP.05.00/2019)is gratefully acknowledged.
文摘This paper presents a new concept called Unified and Integrated Method for a shear deformable beam element.In this method,Timoshenko beam theory is unified and integrated in such a way that takes into account the effect of transverse shear and maintains the shear locking free condition at the same time to generate proper behavior in the analysis of thin to thick beams.The unified and integrated method is applied to finite element analysis(FEA)and isogeometric analysis(IGA)on two-node beam element.This method will be used to analyze uniformly loaded beams with various boundary conditions.A shear influence factor of f,which is a function of beam thickness ratio(L/h),is expressed explicitly as control of the transverse shear strain effect.The analysis gives interesting results showing that applying the unified and integrated method in FEA and IGA will yield exact values of DOF’s and displacement function even when using only a single element.Numerical examples demonstrate the validity and efficiency of the unified and integrated methods.
文摘The morphing technology of hypersonic vehicle improved the flight performance by changing aerodynamic characteristics with shape deformations,but the design of guidance and control system with morphing laws remained to be explored.An Integrated of Guidance,Control and Morphing(IGCM)method for Hypersonic Morphing Vehicle(HMV)was developed in this paper.The IGCM method contributed to an effective solution of morphing characteristic to improve flight performance and reject the disturbance for guidance and control system caused by the morphing system for HMV in gliding phase.The IGCM models were established based on the motion models and aerodynamic models of the variable span vehicle.Then the IGCM method was designed by adaptive block dynamic surface back-stepping method with stability proof.The parallel controlled simulations’results showed the effectiveness in accomplishing the flight mission of IGCM method in glide phase with smaller terminal errors.The velocity loss of HMV was reduced by 32.8%which inferred less flight time and larger terminal flight velocity than invariable span vehicle.Under the condition of large deviations of aerodynamic parameters and atmospheric density,the robustness of IGCM method with variable span was verified.
基金Supported by Science and Technology Foundation of SGCC Research and development of key models for decision support of energy internet companies(NO.SGSDJY00GPJS1900057).
文摘In this study,to develop a benefit-allocation model,in-depth analysis of a distributed photovoltaic-powergeneration carport and energy-storage charging-pile project was performed;the model was developed using Shapley integrated-empowerment benefit-distribution method.First,through literature survey and expert interview to identify the risk factors at various stages of the project,a dynamic risk-factor indicator system is developed.Second,to obtain a more meaningful risk-calculation result,the subjective and objective weights are combined,the weights of the risk factors at each stage are determined by the expert scoring method and entropy weight method,and the interest distribution model based on multi-dimensional risk factors is established.Finally,an example is used to verify the rationality of the method for the benefit distribution of the charging-pile project.The results of the example indicate that the limitations of the Shapley method can be reasonably avoided,and the applicability of the model for the benefit distribution of the charging-pile project is verified.
基金Project supported by the National Natural Science Foundation of China(No.11925204)the 111 Project(No.B14044)。
文摘A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains.To realize this method,we design a new wavelet basis function,by which we construct a fifth-order numerical scheme for the approximation of multi-dimensional functions and their multiple integrals defined in complex domains.In the solution of differential equations,various derivatives of the unknown function are denoted as new functions.Then,the integral relations between these functions are applied in terms of wavelet approximation of multiple integrals.Therefore,the original equation with derivatives of various orders can be converted to a system of algebraic equations with discrete nodal values of the highest-order derivative.During the application of the proposed method,boundary conditions can be automatically included in the integration operations,and relevant matrices can be assured to exhibit perfect sparse patterns.As examples,we consider several second-order mathematics problems defined on regular and irregular domains and the fourth-order bending problems of plates with various shapes.By comparing the solutions obtained by the proposed method with the exact solutions,the new multiresolution method is found to have a convergence rate of fifth order.The solution accuracy of this method with only a few hundreds of nodes can be much higher than that of the finite element method(FEM)with tens of thousands of elements.In addition,because the accuracy order for direct approximation of a function using the proposed basis function is also fifth order,we may conclude that the accuracy of the proposed method is almost independent of the equation order and domain complexity.
基金National Key R&D Program of China(Grant No.2018YFB1307800)National Natural Science Foundation of China(Grant Nos.51875391,51675366)Tianjin Science and Technology Planning Project(Grant Nos.18YFS DZC00010,18YFZCSF00590).
文摘Topology and performance are the two main topics dealt in the development of robotic mechanisms.However,it is still a challenge to connect them by integrating the modeling and design process of both parts in a unified frame.As the properties associated with topology and performance,finite motion and instantaneous motion of the robot play key roles in the procedure.On the purpose of providing a fundamental preparation for integrated modeling and design,this paper carries out a review on the existing unified mathematic frameworks for motion description and computation,involving matrix Lie group and Lie algebra,dual quaternion and pure dual quaternion,finite screw and instantaneous screw.Besides the application in robotics,the review of the work from these mathematicians concentrates on the description,composition and intersection operations of the finite and instantaneous motions,especially on the exponential-differential maps which connect the two sides.Furthermore,an in-depth discussion is worked out by investigating the algebraical relationship among these methods and their further progress in integrated robotic development.The presented review offers insightful investigation to the motion description and computation,and therefore would help designers to choose appropriate mathematical tool in the integrated design and modeling and design of mechanisms and robots.
文摘The objective is to develop a model considering demand dependent on selling price and deterioration occurs after a certain period of time, which follows two-parameter Weibull distribution. Shortages are allowed and fully backlogged. Fuzzy optimal solution is obtained by considering hexagonal fuzzy numbers and for defuzzification Graded Mean Integration Representation Method. A numerical example is provided for the illustration of crisp and fuzzy, both models. To observe the effect of changes in parameters, sensitivity analysis is carried out.
文摘Infertility is one of the difficult complicated diseases. Many couples suffer from it. The pathogenesis is very complicated. The imbalance or lesson of any link of the reproductive system can cause infertility. This paper summarizes the treatment of female infertility by integrated Traditional Chinese Medicine (TCM) and Western Medicine (WM) which can not only improve the ovulation rate and pregnancy rate, but also decrease the complications. The effects are better than that by TCM or WM only. Therefore, the coupling method is worth to be used widely in clinical practice.
基金supported by the National Natural Science Foundation of China (Grant No. 50879090)
文摘The singularities, oscillatory performances and the contributing factors to the 3-'D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.
基金The Project is supported by the National Natural Science Foundation of China
文摘In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金National Science Foundation(NSF)under grant No.CMMI-0748111
文摘This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration (iterative) and the operator-splitting (non-iterative) is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom (SDOF) cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1:2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1:2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
基金supported by National Key Basic Research Program (973 Program, Grant No. 2011CB706804)National Natural Science Foundation of China (Grant No. 50805093)Science & Technology Commission of Shanghai Municipality, China (Grant No. 09QH1401500)
文摘As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.
文摘The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems.
基金The project supported by the National Natural Science Foundation of China
文摘In this paper,the field method for solving the equations of motion of holonomic nonconservative systems is extended to nonholonomic systems with constant mass and with variable mass.Two examples are given to illustrate its application.
基金funded by National Natural Science Foundation of China(NSFC)under Grant Nos.11702238,51904202,and 11902212Nanhu Scholars Program for Young Scholars of XYNU.
文摘This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.
基金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