In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,wh...In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,when the performance index is complicated,because one may encounter a two-point boundary value problem of nonlinear differential algebraic equations.To be a numerical method,it is hard to exactly preserve all the specified conditions,which might deteriorate the accuracy of numerical solution.With this in mind,we develop a novel algorithm to find the solution of the optimal control problem of nonlinear Duffing oscillator,which can exactly satisfy all the required conditions for the minimality of the performance index.A new idea of shape functions method(SFM)is introduced,from which we can transform the optimal control problems to the initial value problems for the new variables,whose initial values are given arbitrarily,and meanwhile the terminal values are determined iteratively.Numerical examples confirm the high-performance of the iterative algorithms based on the SFM,which are convergence fast,and also provide very accurate solutions.The new algorithm is robust,even large noise is imposed on the input data.展开更多
The penalty function method of continuum shape optimization and its sensitivity analysis technique are presented. A relatively simple integrated shape optimization system is developed and used to optimize the design o...The penalty function method of continuum shape optimization and its sensitivity analysis technique are presented. A relatively simple integrated shape optimization system is developed and used to optimize the design of the inner frame shape of a three-axis test table. The result shows that the method converges well, and the system is stable and reliable.展开更多
Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with...Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.展开更多
Using the method of undetermined coefficients, we construct a set of shape function spaces of nine-node triangular plate elements converging for any meshes, which generalize Spect's element and Veubeke's element.
A nonlinear finite element method is applied to observe how inclusion shape influence the thermal response of a ceramic-metal functionally graded material (FGM). The elastic and plastic behaviors of the layers which a...A nonlinear finite element method is applied to observe how inclusion shape influence the thermal response of a ceramic-metal functionally graded material (FGM). The elastic and plastic behaviors of the layers which are two-phase isotropic composites consisting of randomly oriented elastic spheroidal Inclusions and a ductile matrix are predicted by cc mean field method. The prediction results show that inclusion shape has remarkable influence on the overall behavior of the composite. The consequences of the thermal response analysis of the FGM are that the response is dependent on inclusion shape and its composition profile cooperatively and that the plastic behavior of each layer should be taken into account in optimum design of a ceramic-metal FGM.展开更多
In this paper, the basic formulae for the semi-analytical graded FEM on FGM members are derived. Since FGM parameters vary along three space coordinates, the parameters can be integrated in mechanical equations. There...In this paper, the basic formulae for the semi-analytical graded FEM on FGM members are derived. Since FGM parameters vary along three space coordinates, the parameters can be integrated in mechanical equations. Therefore with the parameters of a given FGM plate, problems of FGM plate under various conditions can be solved. The approach uses 1D discretization to obtain 3D solutions, which is proven to be an effective numerical method for the mechanical analyses of FGM structures. Examples of FGM plates with complex shapes and various holes are presented.展开更多
The state-space method is employed to evaluate the modal parameters of functionally graded, magneto-electro-elastic, and multilayered plates. Based on the assumption that the properties of the functionally graded mate...The state-space method is employed to evaluate the modal parameters of functionally graded, magneto-electro-elastic, and multilayered plates. Based on the assumption that the properties of the functionally graded material are exponential, the state equation of structural vibration which takes the displacement and stress of the structure as state variables is derived. The natural frequencies and modal shapes are calculated based on the general solutions of the state equation and boundary conditions given in this paper. The influence of the functionally graded exponential factor on the elastic displacement, electric, and magnetic fields of the structure are discussed by assuming a sandwich plate model with different stacking sequences.展开更多
In this work,the electronic transport properties of Z-shaped silicene nanoribbon（ZsSiNR） structure are investigated.The calculations are based on the tight-binding model and Green＇s function method in Landauer-Biit...In this work,the electronic transport properties of Z-shaped silicene nanoribbon（ZsSiNR） structure are investigated.The calculations are based on the tight-binding model and Green＇s function method in Landauer-Biittiker formalism,in which the electronic density of states（DOS）,transmission probability,and current-voltage characteristics of the system are calculated,numerically.It is shown that the geometry of the ZsSiNR structure can play an important role to control the electron transport through the system.It is observed that the intensity of electron localization at the edges of the ZsSiNR decreases with the increase of the spin-orbit interaction（SOI） strength.Also,the semiconductor to metallic transition occurs by increasing the SOI strength.The present theoretical results may be useful to design silicene-based devices in nanoelectronics.展开更多
This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetr...This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.展开更多
The differential transformation method (DTM) is applied to investigate free vibration of functionally graded beams supported by arbitrary boundary conditions, including various types of elastically end constraints. Th...The differential transformation method (DTM) is applied to investigate free vibration of functionally graded beams supported by arbitrary boundary conditions, including various types of elastically end constraints. The material properties of functionally graded beams are assumed to obey the power law distribution. The main advantages of this method are known for its excellence in high accuracy with small computational expensiveness. The DTM also provides all natural frequencies and mode shapes without any frequency missing. Fundamental frequencies as well as their higher frequencies and mode shapes are presented. The significant aspects such as boundary conditions, values of translational and rotational spring constants and the material volume fraction index on the natural frequencies and mode shapes are discussed. For elastically end constraints, some available results of special cases for isotropic beams are used to validate the present results. The new frequency results and mode shapes of functionally graded beams resting on elastically end constraints are presented.展开更多
The exact solutions of the Schr6dinger equation with the double ring-shaped Coulomb potential are presented, including the bound states, continuous states on the ＂k/2π scale＂, and the calculation formula of the pha...The exact solutions of the Schr6dinger equation with the double ring-shaped Coulomb potential are presented, including the bound states, continuous states on the ＂k/2π scale＂, and the calculation formula of the phase shifts. The polar angular wave functions are expressed by constructing the so-called super-universal associated Legendre polynomials. Some special cases are discussed in detail.展开更多
扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数...扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。展开更多
文摘In the optimal control problem of nonlinear dynamical system,the Hamiltonian formulation is useful and powerful to solve an optimal control force.However,the resulting Euler-Lagrange equations are not easy to solve,when the performance index is complicated,because one may encounter a two-point boundary value problem of nonlinear differential algebraic equations.To be a numerical method,it is hard to exactly preserve all the specified conditions,which might deteriorate the accuracy of numerical solution.With this in mind,we develop a novel algorithm to find the solution of the optimal control problem of nonlinear Duffing oscillator,which can exactly satisfy all the required conditions for the minimality of the performance index.A new idea of shape functions method(SFM)is introduced,from which we can transform the optimal control problems to the initial value problems for the new variables,whose initial values are given arbitrarily,and meanwhile the terminal values are determined iteratively.Numerical examples confirm the high-performance of the iterative algorithms based on the SFM,which are convergence fast,and also provide very accurate solutions.The new algorithm is robust,even large noise is imposed on the input data.
文摘The penalty function method of continuum shape optimization and its sensitivity analysis technique are presented. A relatively simple integrated shape optimization system is developed and used to optimize the design of the inner frame shape of a three-axis test table. The result shows that the method converges well, and the system is stable and reliable.
文摘Accurately approximating higher order derivatives is an inherently difficult problem. It is shown that a random variable shape parameter strategy can improve the accuracy of approximating higher order derivatives with Radial Basis Function methods. The method is used to solve fourth order boundary value problems. The use and location of ghost points are examined in order to enforce the extra boundary conditions that are necessary to make a fourth-order problem well posed. The use of ghost points versus solving an overdetermined linear system via least squares is studied. For a general fourth-order boundary value problem, the recommended approach is to either use one of two novel sets of ghost centers introduced here or else to use a least squares approach. When using either ghost centers or least squares, the random variable shape parameter strategy results in significantly better accuracy than when a constant shape parameter is used.
文摘Using the method of undetermined coefficients, we construct a set of shape function spaces of nine-node triangular plate elements converging for any meshes, which generalize Spect's element and Veubeke's element.
基金Funded by National Science Foundation of China(Grant:1987205).
文摘A nonlinear finite element method is applied to observe how inclusion shape influence the thermal response of a ceramic-metal functionally graded material (FGM). The elastic and plastic behaviors of the layers which are two-phase isotropic composites consisting of randomly oriented elastic spheroidal Inclusions and a ductile matrix are predicted by cc mean field method. The prediction results show that inclusion shape has remarkable influence on the overall behavior of the composite. The consequences of the thermal response analysis of the FGM are that the response is dependent on inclusion shape and its composition profile cooperatively and that the plastic behavior of each layer should be taken into account in optimum design of a ceramic-metal FGM.
基金Project supported by the National Natural Science Foundation of China (No. 10432030)
文摘In this paper, the basic formulae for the semi-analytical graded FEM on FGM members are derived. Since FGM parameters vary along three space coordinates, the parameters can be integrated in mechanical equations. Therefore with the parameters of a given FGM plate, problems of FGM plate under various conditions can be solved. The approach uses 1D discretization to obtain 3D solutions, which is proven to be an effective numerical method for the mechanical analyses of FGM structures. Examples of FGM plates with complex shapes and various holes are presented.
基金Project supported by the National Natural Science Foundation of China (No. 50575172).
文摘The state-space method is employed to evaluate the modal parameters of functionally graded, magneto-electro-elastic, and multilayered plates. Based on the assumption that the properties of the functionally graded material are exponential, the state equation of structural vibration which takes the displacement and stress of the structure as state variables is derived. The natural frequencies and modal shapes are calculated based on the general solutions of the state equation and boundary conditions given in this paper. The influence of the functionally graded exponential factor on the elastic displacement, electric, and magnetic fields of the structure are discussed by assuming a sandwich plate model with different stacking sequences.
基金Project supported by the Sari Branch,Islamic Azad University,Iran Grant No.1-24850
文摘In this work,the electronic transport properties of Z-shaped silicene nanoribbon（ZsSiNR） structure are investigated.The calculations are based on the tight-binding model and Green＇s function method in Landauer-Biittiker formalism,in which the electronic density of states（DOS）,transmission probability,and current-voltage characteristics of the system are calculated,numerically.It is shown that the geometry of the ZsSiNR structure can play an important role to control the electron transport through the system.It is observed that the intensity of electron localization at the edges of the ZsSiNR decreases with the increase of the spin-orbit interaction（SOI） strength.Also,the semiconductor to metallic transition occurs by increasing the SOI strength.The present theoretical results may be useful to design silicene-based devices in nanoelectronics.
基金National Natural Science Foundation of China Under Grant No.51278382
文摘This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.
文摘The differential transformation method (DTM) is applied to investigate free vibration of functionally graded beams supported by arbitrary boundary conditions, including various types of elastically end constraints. The material properties of functionally graded beams are assumed to obey the power law distribution. The main advantages of this method are known for its excellence in high accuracy with small computational expensiveness. The DTM also provides all natural frequencies and mode shapes without any frequency missing. Fundamental frequencies as well as their higher frequencies and mode shapes are presented. The significant aspects such as boundary conditions, values of translational and rotational spring constants and the material volume fraction index on the natural frequencies and mode shapes are discussed. For elastically end constraints, some available results of special cases for isotropic beams are used to validate the present results. The new frequency results and mode shapes of functionally graded beams resting on elastically end constraints are presented.
基金Project supported by the National Natural Science Foundation of China(Grant No.11275165)the Natural Science Foundation of Jiangsu Province,China(Grant No.BK2010291)partly by Secretaria de Investigacio'ny Posgrado de Instituto Polite'cnico Nacional,Mexico(Grant No.20131150-SIP-IPN)
文摘The exact solutions of the Schr6dinger equation with the double ring-shaped Coulomb potential are presented, including the bound states, continuous states on the ＂k/2π scale＂, and the calculation formula of the phase shifts. The polar angular wave functions are expressed by constructing the so-called super-universal associated Legendre polynomials. Some special cases are discussed in detail.
文摘扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。
