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展开更多
In this paper,the path integral solutions for a general n-dimensional stochastic differential equa-tions(SDEs)withα-stable Lévy noise are derived and verified.Firstly,the governing equations for the solutions of...In this paper,the path integral solutions for a general n-dimensional stochastic differential equa-tions(SDEs)withα-stable Lévy noise are derived and verified.Firstly,the governing equations for the solutions of n-dimensional SDEs under the excitation ofα-stable Lévy noise are obtained through the characteristic function of stochastic processes.Then,the short-time transition probability density func-tion of the path integral solution is derived based on the Chapman-Kolmogorov-Smoluchowski(CKS)equation and the characteristic function,and its correctness is demonstrated by proving that it satis-fies the governing equation of the solution of the SDE,which is also called the Fokker-Planck-Kolmogorov equation.Besides,illustrative examples are numerically considered for highlighting the feasibility of the proposed path integral method,and the pertinent Monte Carlo solution is also calculated to show its correctness and effectiveness.展开更多
This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integ...This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.展开更多
By elastic-plastic finite deformation finite element analysis (FEA), the conservation of J-integral is investigated in detail for a welded joint with an overmatchingweld in plane stress case. It is indicated that J-in...By elastic-plastic finite deformation finite element analysis (FEA), the conservation of J-integral is investigated in detail for a welded joint with an overmatchingweld in plane stress case. It is indicated that J-integral is path dependent under various conditions at least in the cases studied in this paper. Meanwhile, the above conclusions are verified by the hybrid method results in which combined Moire interferometry with FEA.展开更多
The path integral Monte Carlo(PIMC) method is employed to study the thermal properties of C70 with one, two,and three H2 molecules confined in the cage, respectively. The interaction energies and vibrationally average...The path integral Monte Carlo(PIMC) method is employed to study the thermal properties of C70 with one, two,and three H2 molecules confined in the cage, respectively. The interaction energies and vibrationally averaged spatial distributions under different temperatures are calculated to evaluate the stabilities of(H2)n@C70(n = 1, 2, 3). The results show that(H2)2@C70is more stable than H2@C70. The interaction energy slowly changes in a large temperature range,so temperature has little effect on the stability of the system. For H2@C70and(H2)2@C70, the interaction energies keep negative; however, when three H2 molecules are in the cage, the interaction energy rapidly increases to a positive value.This implies that at most two H2 molecules can be trapped by C70. With an increase of temperature, the peak of the spatial distribution gradually shifts away from the center of the cage, but the maximum distance from the center of H2 molecule to the cage center is much smaller than the average radius of C70.展开更多
A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in t...A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.展开更多
In this paper, by means of the maximum circle tensile stress on curve of constant ω and stress intensity factors by a path independent contour integral method, trajectories of maxed mode crack propagation are simulat...In this paper, by means of the maximum circle tensile stress on curve of constant ω and stress intensity factors by a path independent contour integral method, trajectories of maxed mode crack propagation are simulated through numerical manifold method. The crack propagation is traced dynamically by modifying the neighboring connection between the crack-top and nodes within elements in the calculating process. This method has the advantages such as less modified area, easiness of programming, high realizability and so on. Then a single sharp nicked specimen is used to verified the numerical result. It is shown that the provided method is reasonable and effective.展开更多
基金Project supported by the National Natural Science Foundation of China(Nos.10902020 and 10721062)
文摘An improved precise integration method （IPIM） for solving the differential Riccati equation （DRE） is presented. The solution to the DRE is connected with the exponential of a Hamiltonian matrix, and the precise integration method （PIM） for solving the DRE is connected with the scaling and squaring method for computing the exponential of a matrix. The error analysis of the scaling and squaring method for the exponential of a matrix is applied to the PIM of the DRE. Based ,on the error analysis, the criterion for choosing two parameters of the PIM is given. Three kinds of IPIMs for solving the DRE are proposed. The numerical examples machine accuracy solutions. show that the IPIM is stable and gives the
基金This work was supported by the Key International(Regional)Joint Research Program of the National Natural Science Foundation of China(No.12120101002).
文摘In this paper,the path integral solutions for a general n-dimensional stochastic differential equa-tions(SDEs)withα-stable Lévy noise are derived and verified.Firstly,the governing equations for the solutions of n-dimensional SDEs under the excitation ofα-stable Lévy noise are obtained through the characteristic function of stochastic processes.Then,the short-time transition probability density func-tion of the path integral solution is derived based on the Chapman-Kolmogorov-Smoluchowski(CKS)equation and the characteristic function,and its correctness is demonstrated by proving that it satis-fies the governing equation of the solution of the SDE,which is also called the Fokker-Planck-Kolmogorov equation.Besides,illustrative examples are numerically considered for highlighting the feasibility of the proposed path integral method,and the pertinent Monte Carlo solution is also calculated to show its correctness and effectiveness.
基金The project supported by National Natural Science Foundation of China(9713008)Zhejiang Natural Science Foundation Special Funds No. RC.9601
文摘This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.
文摘By elastic-plastic finite deformation finite element analysis (FEA), the conservation of J-integral is investigated in detail for a welded joint with an overmatchingweld in plane stress case. It is indicated that J-integral is path dependent under various conditions at least in the cases studied in this paper. Meanwhile, the above conclusions are verified by the hybrid method results in which combined Moire interferometry with FEA.
基金supported by the National Natural Science Foundation of China(Grant Nos.11474207 and 11374217)
文摘The path integral Monte Carlo(PIMC) method is employed to study the thermal properties of C70 with one, two,and three H2 molecules confined in the cage, respectively. The interaction energies and vibrationally averaged spatial distributions under different temperatures are calculated to evaluate the stabilities of(H2)n@C70(n = 1, 2, 3). The results show that(H2)2@C70is more stable than H2@C70. The interaction energy slowly changes in a large temperature range,so temperature has little effect on the stability of the system. For H2@C70and(H2)2@C70, the interaction energies keep negative; however, when three H2 molecules are in the cage, the interaction energy rapidly increases to a positive value.This implies that at most two H2 molecules can be trapped by C70. With an increase of temperature, the peak of the spatial distribution gradually shifts away from the center of the cage, but the maximum distance from the center of H2 molecule to the cage center is much smaller than the average radius of C70.
基金the National Natural Science Foundation of China(No.11672064)。
文摘A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.
基金Funded by the National Natural Science Foundation of China (No. 10272033) and Guangdong Provincial Natural Science Foundation(Nos.04105386,5300090 and 05001844).
文摘In this paper, by means of the maximum circle tensile stress on curve of constant ω and stress intensity factors by a path independent contour integral method, trajectories of maxed mode crack propagation are simulated through numerical manifold method. The crack propagation is traced dynamically by modifying the neighboring connection between the crack-top and nodes within elements in the calculating process. This method has the advantages such as less modified area, easiness of programming, high realizability and so on. Then a single sharp nicked specimen is used to verified the numerical result. It is shown that the provided method is reasonable and effective.
