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.展开更多
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.展开更多
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.
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.展开更多
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.展开更多
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.展开更多
扩展有限元(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.
基金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.
摘要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.
摘要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.
基金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.
基金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.
摘要扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。