The exact minimax penalty function method is used to solve a noncon- vex differentiable optimization problem with both inequality and equality constraints. The conditions for exactness of the penalization for the exac...The exact minimax penalty function method is used to solve a noncon- vex differentiable optimization problem with both inequality and equality constraints. The conditions for exactness of the penalization for the exact minimax penalty function method are established by assuming that the functions constituting the considered con- strained optimization problem are invex with respect to the same function η （with the exception of those equality constraints for which the associated Lagrange multipliers are negative these functions should be assumed to be incave with respect to η）. Thus, a threshold of the penalty parameter is given such that, for all penalty parameters exceeding this threshold, equivalence holds between the set of optimal solutions in the considered constrained optimization problem and the set of minimizer in its associated penalized problem with an exact minimax penalty function. It is shown that coercivity is not suf- ficient to prove the results.展开更多
Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dth...Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dthermoelastic solution has been reported for a uniformly heated penny-shaped crack in an infinite magnetoelectric space,with impermeable electromagnetic conditions assumed on the crack faces.Exact 3Dsolutions for the penny-shaped crack subjected to uniform or point temperature load are further presented here when the crack faces are electrically and magnetically permeable.The solutions,obtained by the potential theory method,are exact in the sense that all field variables are explicitly derived and expressed in terms of elementary functions.Along with the previously reported solution,the limits or bounds of the stress intensity factor at the crack-tip for a practical crack can be identified.展开更多
In this paper, based on the step reduction method, a new method, the exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational principle, it can be appl...In this paper, based on the step reduction method, a new method, the exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a triangle noncompatible element with 6 degrees of freedom is derived to solve the bending of nonhomogeneous plate. The convergence of displacements and stress resultants which have satisfactory numerical precision is proved. Numerical examples are given at the end of this paper, which indicate satisfactory results of stress resultants and displacements can be obtained by the present method.展开更多
The exact analytic method was given by [1] . It can be used for arbitrary variable coefficient differential equations and the solution obtained can have the second order convergent precision. In this paper, a new high...The exact analytic method was given by [1] . It can be used for arbitrary variable coefficient differential equations and the solution obtained can have the second order convergent precision. In this paper, a new high precision algorithm is given based on [1], through a bending problem of variable cross-section beams. It can have the fourth convergent precision without increasing computation work. The present computation method is not only simple but also fast. The numerical examples are given at the end of this paper which indicate that the high convergent precision can be obtained using only a few elements. The correctness of the theory in this paper is confirmed.展开更多
In[1], the exact analytic method for the solution of differential equation with variable coefficients was suggested and an analytic expression of solution was given by initial parameter algorithm. But to some problems...In[1], the exact analytic method for the solution of differential equation with variable coefficients was suggested and an analytic expression of solution was given by initial parameter algorithm. But to some problems such as the bending, free vibration and buckling of nonhomogeneous long cylinders, it is difficult to obtain their solutions by the initial parameter algorithm on computer. In this paper, the substructure computational algorithm for the exact analytic method is presented through the bending of non-homogeneous long cylindrical shell. This substructure algorithm can he applied to solve the problems which can not he calculated by the initial parameter algorithm on computer. Finally, the problems can he reduced to solving a low order system of algehraic equations like the initial parameter algorithm Numerical examples are given and compared with the initial para-algorithm at the end of the paper, which confirms the correctness of the substructure computational algorithm.展开更多
In this paper, based on the step reduction method and exaet analytic method, a new method, theexacl element method for constructing finite element, is presented. Since the near method doesn't need varialional prin...In this paper, based on the step reduction method and exaet analytic method, a new method, theexacl element method for constructing finite element, is presented. Since the near method doesn't need varialional principle, it can he applied to solve nun-positive and positive definite partial differcntial equations with arbitral varutble coefficients. By this method, a triangle noncompatible element with 15 degrees of freedom is derived to solve the bending of nonhomogenous Reissner's plate. Because the displacement parameters at the nodal point only contain deflection and rotation angle, it is convenient to deal with arbitrary boundary conditions. In this paper, the convergcnceof displacement and stress resultants is proved. The element obtained by the present method can be used for thin and thick plates as well, hour numerical examples are given at the end of this paper, which indicates that we can obtain satisfactory results and have higher numerical precision.展开更多
Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations unde...Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations under arbitrary boundary condition. By this method, the general computation formal is obtained. Its convergence in proved. We can get analytic expressions which converge to exact solution and its higher order derivatives uniformy Four numerical examples are given, which indicate that satisfactory results can he obtanedby this method.展开更多
In this paper, based on the step reduction method[1] and exact analytic method[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational ...In this paper, based on the step reduction method[1] and exact analytic method[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a quadrilateral noncompatible element with 8 degrees of freedom is derived for the solution of plane problem. Since Jacobi's transformation is not applied, the present element may degenerate into a triangle element. It is convenient to use the element in engineering. In this paper, the convergence is proved. Numerical examples are given at the endof this paper, which indicate satisfactory results of stress and displacements can be obtained and have higher numerical precision in nodes.展开更多
In this paper, we use the classical Lie group symmetry method to get the Lie point symmetries of the （2＋1）-dimensional hyperbolic nonlinear Schr6dinger （HNLS） equation and reduce the （2＋1）-dimensional HNLS equ...In this paper, we use the classical Lie group symmetry method to get the Lie point symmetries of the （2＋1）-dimensional hyperbolic nonlinear Schr6dinger （HNLS） equation and reduce the （2＋1）-dimensional HNLS equation to some （1 ＋ 1 ）-dimensional partial differential systems. Finally, many exact travelling solutions of the （2＋1）-dimensional HNLS equation are obtained by the classical Lie symmetry reduced method.展开更多
In this paper, the nonlinear axial symmetric deformation problem of nonhomogeneous ring- and stringer-stiffened shells is first solved by the exact analytic method. An analytic expression of displacements and stress r...In this paper, the nonlinear axial symmetric deformation problem of nonhomogeneous ring- and stringer-stiffened shells is first solved by the exact analytic method. An analytic expression of displacements and stress resultants is obtained and its convergence is proved. Displacements and stress resultants converge to exact solution uniformly. Finally, it is only necessary to solve a system of linear algebraic equations with two unknowns. Four numerical examples are given at the end of the paper which indicate that satisfactory results can be obtained by the exact analytic method.展开更多
In this paper, the boundary control problem of a distributed parameter system described by the Schrodinger equation posed on finite interval α≤x≤β：{ iyt ＋yzz＋｜y｜^2y = 0, y（α, t） = h1 （t）, y（β, t）=h2...In this paper, the boundary control problem of a distributed parameter system described by the Schrodinger equation posed on finite interval α≤x≤β：{ iyt ＋yzz＋｜y｜^2y = 0, y（α, t） = h1 （t）, y（β, t）=h2（t） for t〉0 （S）is considered. It is shown that by choosing appropriate control inputs （hi）, （j = 1, 2） one can always guide the system （S） from a given initial state φ∈H^S（α,β）, （s ∈ R） to a terminal state φ∈ H^s（α,β）, in the time period [0, T]. The exact boundary controllability is obtained by considering a related initial value control problem of SchrSdinger equation posed on the whole line R. The discovered smoothing properties of Schrodinger equation have played important roles in our approach; this may be the first step to prove the results on boundary controllability of （semi-linear） nonlinear Schrodinger equation.展开更多
This paper presents the development and implementation of an innovative mixed integer programming based mathematical model for an open pit mining operation with Grade Engineering framework.Grade Engineering comprises ...This paper presents the development and implementation of an innovative mixed integer programming based mathematical model for an open pit mining operation with Grade Engineering framework.Grade Engineering comprises a range of coarse-separation based pre-processing techniques that separate the desirable(i.e.high-grade)and undesirable(i.e.low-grade or uneconomic)materials and ensure the delivery of only selected quantity of high quality(or high-grade)material to energy,water,and cost-intensive processing plant.The model maximizes the net present value under a range of operational and processing constraints.Given that the proposed model is computationally complex,the authors employ a data preprocessing procedure and then evaluate the performance of the model at several practical instances using computation time,optimality gap,and the net present value as valid measures.In addition,a comparison of the proposed and traditional(without Grade Engineering)models reflects that the proposed model outperforms the traditional formulation.展开更多
A new nonlinear partial differential equation （PDE） in 2＋1 dimensions is obtained from the mKP equation by means of an asymptotically exact reduction method based on Fourier expansion and spatio-temporal resealing....A new nonlinear partial differential equation （PDE） in 2＋1 dimensions is obtained from the mKP equation by means of an asymptotically exact reduction method based on Fourier expansion and spatio-temporal resealing. In order to demonstrate integrability property of the new equation, the corresponding Lax pair is obtained by applying the reduction technique to the Lax pair of the mKP equation.展开更多
In this paper, we made a new breakthrough, which proposes a new recursion–transform（RT） method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found ...In this paper, we made a new breakthrough, which proposes a new recursion–transform（RT） method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found the exact potential formulae of arbitrary m × n cobweb and fan networks by the RT method, and the potential formulae of infinite and semi-infinite networks are derived. As applications, a series of interesting corollaries of potential formulae are given by using the general formula, the equivalent resistance formula is deduced by using the potential formula, and we find a new trigonometric identity by comparing two equivalence results with different forms.展开更多
A class of lump solutions of(2+1)-dimensional Boussinesq equation are obtained with the help of Maple by using Hirota bilinear method.Some contour plots with different determinant values are sequentially made to show ...A class of lump solutions of(2+1)-dimensional Boussinesq equation are obtained with the help of Maple by using Hirota bilinear method.Some contour plots with different determinant values are sequentially made to show that the corresponding lump solution tends to zero when the determinant approaches zero.The particular lump solutions with specific values of the involved parameters are plotted,as illustrative examples.展开更多
The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite el...The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite element model for two-dimensional dynamic analysis of armored cable is presented. This model accounts for the geometric nonlinearities of large displacement of the armored cable, and effects of axial load and bending stiffness. The governing equations are derived by consistent linearization and finite element discretization of the total weak form of the armored cable system, and solved by the Newmark time integration method. To make the solution procedure avoid falling into the local extreme points, a simple adaptive stepping strategy is proposed. The presented model is validated via actual measured data. Results for dynamic configurations, motion and tension of both ends of the armored cable, and resonance-zone are presented for two numerical cases, including the dynamic analysis under the case of only ship heave motion and the case of joint action of ship heave motion and ocean current. The dynamics analysis can provide important reference for the design or product selection of the armored cable in a deep-sea ROV system so as to improve the safety of its marine operation under the sea state of 4 or above.展开更多
A new generalized transformation method is differential equation. As an application of the method, we presented to find more exact solutions of nonlinear partial choose the （3＋1）-dimensional breaking soliton equati...A new generalized transformation method is differential equation. As an application of the method, we presented to find more exact solutions of nonlinear partial choose the （3＋1）-dimensional breaking soliton equation to illustrate the method. As a result many types of explicit and exact traveling wave solutions, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic function solutions, and rational solutions, are obtained. The new method can be extended to other nonlinear partial differential equations in mathematical physics.展开更多
Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,litt...Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,little attention has been paid to the central states,which are exponentially close to each other in terms of system size.We propose a delta-Davidson(DELDAV)method to efficiently find such interior(including the central)states in many-spin systems.The DELDAV method utilizes a delta filter in Chebyshev polynomial expansion combined with subspace diagonalization to overcome the nearly degenerate problem.Numerical experiments on Ising spin chain and spin glass shards show the correctness,efficiency,and robustness of the proposed method in finding the interior states as well as the ground states.The sought interior states may be employed to identify many-body localization phase,quantum chaos,and extremely long-time dynamical structure.展开更多
For most of the conventional crystals with low-index surfaces, the hopping between the nearest neighbor （1NN） crystal planes （CPs） is dominant and the ones from the nNN （2 〈 n 〈 ∞） CPs are relatively weak, co...For most of the conventional crystals with low-index surfaces, the hopping between the nearest neighbor （1NN） crystal planes （CPs） is dominant and the ones from the nNN （2 〈 n 〈 ∞） CPs are relatively weak, considered as small perturbations. The recent theoretical analysisIll has demonstrated the absence of surface states at the level of the hopping approximation between the INN CPs when the original infinite crystal has the geometric reflection symmetry （GRS） for each CP. Meanwhile, based on the perturbation theory, it has also been shown that small perturbations from the hopping between the nNN （2 〈 n 〈 ∞） CPs and surface relaxation have no impact on the above conclusion. However, for the crystals with strong intrinsic spin-orbit coupling （SOC）, the dominant terms of intrinsic SOC associate with two INN bond hoppings. Thus SOC will significantly contribute the hoppings from the INN and/or 2NN CPs except the ones within each CP. Here, we will study the effect of the hopping between the 2NN CPs on the surface states in model crystals with three different type structures （Type I： “……P-P-P-P……”, Type II： “……-P-Q-P-Q……” and Type III：“……P=Q-P=Q……” where P and Q indicate CPs and the signs “-” and “=” mark the distance between the INN CPs）. In terms of analytical and numerical calculations, we study the behavior of surface states in three types after the symmetric/asymmetric hopping from the 2NN CPs is added. We analytically prove that the symmetric hopping from the 2NN CPs cannot induce surface states in Type I when each CP has only one electron mode. The numerical calculations also provide strong support for the conclusion, even up to 5NN. However, in general, the coupling from the 2NN CPs （symmetric and asymmetric） is favorable to generate surface states except Type I with single electron mode only.展开更多
This paper presents a new curved quadrilateral plate element with 12-degree freedom by the exact element method[1]. The method can be used to arbitrary non-positive and positive definite partial differential equations...This paper presents a new curved quadrilateral plate element with 12-degree freedom by the exact element method[1]. The method can be used to arbitrary non-positive and positive definite partial differential equations without variation principle. Using this method, the compatibility conditions between element can be treated very easily, if displacements and stress resultants are continuous at nodes between elements. The displacements and stress resultants obtained by the present method can converge to exact solution and have the second order convergence speed. Numerical examples are given at the end of this paper, which show the excellent precision and efficiency of the new element.展开更多
文摘The exact minimax penalty function method is used to solve a noncon- vex differentiable optimization problem with both inequality and equality constraints. The conditions for exactness of the penalization for the exact minimax penalty function method are established by assuming that the functions constituting the considered con- strained optimization problem are invex with respect to the same function η （with the exception of those equality constraints for which the associated Lagrange multipliers are negative these functions should be assumed to be incave with respect to η）. Thus, a threshold of the penalty parameter is given such that, for all penalty parameters exceeding this threshold, equivalence holds between the set of optimal solutions in the considered constrained optimization problem and the set of minimizer in its associated penalized problem with an exact minimax penalty function. It is shown that coercivity is not suf- ficient to prove the results.
基金This work was supported by the National Natural Sci- ence Foundation of China （11321202） and the Specialized Research Fund for the Doctoral Program of Higher Educa- tion （2013010 1110120）.
文摘Exact solutions of three-dimensional(3D)crack problems are much less in number than those of two-dimensional ones,especially for multi-field coupling media exhibiting a certain kind of material anisotropy.An exact3Dthermoelastic solution has been reported for a uniformly heated penny-shaped crack in an infinite magnetoelectric space,with impermeable electromagnetic conditions assumed on the crack faces.Exact 3Dsolutions for the penny-shaped crack subjected to uniform or point temperature load are further presented here when the crack faces are electrically and magnetically permeable.The solutions,obtained by the potential theory method,are exact in the sense that all field variables are explicitly derived and expressed in terms of elementary functions.Along with the previously reported solution,the limits or bounds of the stress intensity factor at the crack-tip for a practical crack can be identified.
文摘In this paper, based on the step reduction method, a new method, the exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a triangle noncompatible element with 6 degrees of freedom is derived to solve the bending of nonhomogeneous plate. The convergence of displacements and stress resultants which have satisfactory numerical precision is proved. Numerical examples are given at the end of this paper, which indicate satisfactory results of stress resultants and displacements can be obtained by the present method.
文摘The exact analytic method was given by [1] . It can be used for arbitrary variable coefficient differential equations and the solution obtained can have the second order convergent precision. In this paper, a new high precision algorithm is given based on [1], through a bending problem of variable cross-section beams. It can have the fourth convergent precision without increasing computation work. The present computation method is not only simple but also fast. The numerical examples are given at the end of this paper which indicate that the high convergent precision can be obtained using only a few elements. The correctness of the theory in this paper is confirmed.
文摘In[1], the exact analytic method for the solution of differential equation with variable coefficients was suggested and an analytic expression of solution was given by initial parameter algorithm. But to some problems such as the bending, free vibration and buckling of nonhomogeneous long cylinders, it is difficult to obtain their solutions by the initial parameter algorithm on computer. In this paper, the substructure computational algorithm for the exact analytic method is presented through the bending of non-homogeneous long cylindrical shell. This substructure algorithm can he applied to solve the problems which can not he calculated by the initial parameter algorithm on computer. Finally, the problems can he reduced to solving a low order system of algehraic equations like the initial parameter algorithm Numerical examples are given and compared with the initial para-algorithm at the end of the paper, which confirms the correctness of the substructure computational algorithm.
文摘In this paper, based on the step reduction method and exaet analytic method, a new method, theexacl element method for constructing finite element, is presented. Since the near method doesn't need varialional principle, it can he applied to solve nun-positive and positive definite partial differcntial equations with arbitral varutble coefficients. By this method, a triangle noncompatible element with 15 degrees of freedom is derived to solve the bending of nonhomogenous Reissner's plate. Because the displacement parameters at the nodal point only contain deflection and rotation angle, it is convenient to deal with arbitrary boundary conditions. In this paper, the convergcnceof displacement and stress resultants is proved. The element obtained by the present method can be used for thin and thick plates as well, hour numerical examples are given at the end of this paper, which indicates that we can obtain satisfactory results and have higher numerical precision.
文摘Many engineering problems can be reduced to the solution of a variable coefficient differential equation. In this paper, the exact analytic method is suggested to solve variable coefficient differential equations under arbitrary boundary condition. By this method, the general computation formal is obtained. Its convergence in proved. We can get analytic expressions which converge to exact solution and its higher order derivatives uniformy Four numerical examples are given, which indicate that satisfactory results can he obtanedby this method.
文摘In this paper, based on the step reduction method[1] and exact analytic method[2] anew method-exact element method for constructing finite element, is presented. Since the new method doesn 't need the variational principle, it can be applied to solve non-positive and positive definite partial differential equations with arbitrary variable coefficient. By this method, a quadrilateral noncompatible element with 8 degrees of freedom is derived for the solution of plane problem. Since Jacobi's transformation is not applied, the present element may degenerate into a triangle element. It is convenient to use the element in engineering. In this paper, the convergence is proved. Numerical examples are given at the endof this paper, which indicate satisfactory results of stress and displacements can be obtained and have higher numerical precision in nodes.
基金Supported by the National Natural Science Foundation of China under Grant No.10875106
文摘In this paper, we use the classical Lie group symmetry method to get the Lie point symmetries of the （2＋1）-dimensional hyperbolic nonlinear Schr6dinger （HNLS） equation and reduce the （2＋1）-dimensional HNLS equation to some （1 ＋ 1 ）-dimensional partial differential systems. Finally, many exact travelling solutions of the （2＋1）-dimensional HNLS equation are obtained by the classical Lie symmetry reduced method.
文摘In this paper, the nonlinear axial symmetric deformation problem of nonhomogeneous ring- and stringer-stiffened shells is first solved by the exact analytic method. An analytic expression of displacements and stress resultants is obtained and its convergence is proved. Displacements and stress resultants converge to exact solution uniformly. Finally, it is only necessary to solve a system of linear algebraic equations with two unknowns. Four numerical examples are given at the end of the paper which indicate that satisfactory results can be obtained by the exact analytic method.
文摘In this paper, the boundary control problem of a distributed parameter system described by the Schrodinger equation posed on finite interval α≤x≤β：{ iyt ＋yzz＋｜y｜^2y = 0, y（α, t） = h1 （t）, y（β, t）=h2（t） for t〉0 （S）is considered. It is shown that by choosing appropriate control inputs （hi）, （j = 1, 2） one can always guide the system （S） from a given initial state φ∈H^S（α,β）, （s ∈ R） to a terminal state φ∈ H^s（α,β）, in the time period [0, T]. The exact boundary controllability is obtained by considering a related initial value control problem of SchrSdinger equation posed on the whole line R. The discovered smoothing properties of Schrodinger equation have played important roles in our approach; this may be the first step to prove the results on boundary controllability of （semi-linear） nonlinear Schrodinger equation.
文摘This paper presents the development and implementation of an innovative mixed integer programming based mathematical model for an open pit mining operation with Grade Engineering framework.Grade Engineering comprises a range of coarse-separation based pre-processing techniques that separate the desirable(i.e.high-grade)and undesirable(i.e.low-grade or uneconomic)materials and ensure the delivery of only selected quantity of high quality(or high-grade)material to energy,water,and cost-intensive processing plant.The model maximizes the net present value under a range of operational and processing constraints.Given that the proposed model is computationally complex,the authors employ a data preprocessing procedure and then evaluate the performance of the model at several practical instances using computation time,optimality gap,and the net present value as valid measures.In addition,a comparison of the proposed and traditional(without Grade Engineering)models reflects that the proposed model outperforms the traditional formulation.
基金supported by National Natural Science Foundation of China under Grant No. 10575087the Natural Science Foundation of Zhejiang Province under Grant No. 102053
文摘A new nonlinear partial differential equation （PDE） in 2＋1 dimensions is obtained from the mKP equation by means of an asymptotically exact reduction method based on Fourier expansion and spatio-temporal resealing. In order to demonstrate integrability property of the new equation, the corresponding Lax pair is obtained by applying the reduction technique to the Lax pair of the mKP equation.
基金Project supported by the Natural Science Foundation of Jiangsu Province,China(Grant No.BK20161278)
文摘In this paper, we made a new breakthrough, which proposes a new recursion–transform（RT） method with potential parameters to evaluate the nodal potential in arbitrary resistor networks. For the first time, we found the exact potential formulae of arbitrary m × n cobweb and fan networks by the RT method, and the potential formulae of infinite and semi-infinite networks are derived. As applications, a series of interesting corollaries of potential formulae are given by using the general formula, the equivalent resistance formula is deduced by using the potential formula, and we find a new trigonometric identity by comparing two equivalence results with different forms.
基金Supported by the National Natural Science Foundation of China under Grant No.10647112the Fund of Science and Technology Commission of Shanghai Municipality under Grant No.ZX201307000014
文摘A class of lump solutions of(2+1)-dimensional Boussinesq equation are obtained with the help of Maple by using Hirota bilinear method.Some contour plots with different determinant values are sequentially made to show that the corresponding lump solution tends to zero when the determinant approaches zero.The particular lump solutions with specific values of the involved parameters are plotted,as illustrative examples.
基金Project(2008AA09Z201)supported by the National High Technology Research and Development Program of China
文摘The armored cable used in a deep-sea remotely operated vehicle(ROV) may undergo large displacement motion when subjected to dynamic actions of ship heave motion and ocean current. A novel geometrically exact finite element model for two-dimensional dynamic analysis of armored cable is presented. This model accounts for the geometric nonlinearities of large displacement of the armored cable, and effects of axial load and bending stiffness. The governing equations are derived by consistent linearization and finite element discretization of the total weak form of the armored cable system, and solved by the Newmark time integration method. To make the solution procedure avoid falling into the local extreme points, a simple adaptive stepping strategy is proposed. The presented model is validated via actual measured data. Results for dynamic configurations, motion and tension of both ends of the armored cable, and resonance-zone are presented for two numerical cases, including the dynamic analysis under the case of only ship heave motion and the case of joint action of ship heave motion and ocean current. The dynamics analysis can provide important reference for the design or product selection of the armored cable in a deep-sea ROV system so as to improve the safety of its marine operation under the sea state of 4 or above.
基金The project supported by National Natural Science Foundation of China and the Natural Science Foundation of Shandong Province of China
文摘A new generalized transformation method is differential equation. As an application of the method, we presented to find more exact solutions of nonlinear partial choose the （3＋1）-dimensional breaking soliton equation to illustrate the method. As a result many types of explicit and exact traveling wave solutions, which contain solitary wave solutions, trigonometric function solutions, Jacobian elliptic function solutions, and rational solutions, are obtained. The new method can be extended to other nonlinear partial differential equations in mathematical physics.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.91836101,U1930201,and 11574239).
文摘Many numerical methods,such as tensor network approaches including density matrix renormalization group calculations,have been developed to calculate the extreme/ground states of quantum many-body systems.However,little attention has been paid to the central states,which are exponentially close to each other in terms of system size.We propose a delta-Davidson(DELDAV)method to efficiently find such interior(including the central)states in many-spin systems.The DELDAV method utilizes a delta filter in Chebyshev polynomial expansion combined with subspace diagonalization to overcome the nearly degenerate problem.Numerical experiments on Ising spin chain and spin glass shards show the correctness,efficiency,and robustness of the proposed method in finding the interior states as well as the ground states.The sought interior states may be employed to identify many-body localization phase,quantum chaos,and extremely long-time dynamical structure.
基金supported by the National Natural Science Foundation of China(Grant No.11447601)the National Basic Research Program of China(Grant No.2011CB921803)
文摘For most of the conventional crystals with low-index surfaces, the hopping between the nearest neighbor （1NN） crystal planes （CPs） is dominant and the ones from the nNN （2 〈 n 〈 ∞） CPs are relatively weak, considered as small perturbations. The recent theoretical analysisIll has demonstrated the absence of surface states at the level of the hopping approximation between the INN CPs when the original infinite crystal has the geometric reflection symmetry （GRS） for each CP. Meanwhile, based on the perturbation theory, it has also been shown that small perturbations from the hopping between the nNN （2 〈 n 〈 ∞） CPs and surface relaxation have no impact on the above conclusion. However, for the crystals with strong intrinsic spin-orbit coupling （SOC）, the dominant terms of intrinsic SOC associate with two INN bond hoppings. Thus SOC will significantly contribute the hoppings from the INN and/or 2NN CPs except the ones within each CP. Here, we will study the effect of the hopping between the 2NN CPs on the surface states in model crystals with three different type structures （Type I： “……P-P-P-P……”, Type II： “……-P-Q-P-Q……” and Type III：“……P=Q-P=Q……” where P and Q indicate CPs and the signs “-” and “=” mark the distance between the INN CPs）. In terms of analytical and numerical calculations, we study the behavior of surface states in three types after the symmetric/asymmetric hopping from the 2NN CPs is added. We analytically prove that the symmetric hopping from the 2NN CPs cannot induce surface states in Type I when each CP has only one electron mode. The numerical calculations also provide strong support for the conclusion, even up to 5NN. However, in general, the coupling from the 2NN CPs （symmetric and asymmetric） is favorable to generate surface states except Type I with single electron mode only.
基金Outstanding Education Fund and Doctor Point Fund of National Education Committee and the National Science Foundation of China
文摘This paper presents a new curved quadrilateral plate element with 12-degree freedom by the exact element method[1]. The method can be used to arbitrary non-positive and positive definite partial differential equations without variation principle. Using this method, the compatibility conditions between element can be treated very easily, if displacements and stress resultants are continuous at nodes between elements. The displacements and stress resultants obtained by the present method can converge to exact solution and have the second order convergence speed. Numerical examples are given at the end of this paper, which show the excellent precision and efficiency of the new element.
