C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolati...C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.展开更多
A hybrid natural element method（HNEM） for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation sy...A hybrid natural element method（HNEM） for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation system of the HNEM for viscoelasticity problems is obtained using the Hellinger–Reissner variational principle. In contrast to the natural element method（NEM）, the HNEM can directly obtain the nodal stresses, which have higher precisions than those obtained using the moving least-square（MLS） approximation. Some numerical examples are given to demonstrate the validity and superiority of this HNEM.展开更多
We present the hybrid natural element method（HNEM） for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements an...We present the hybrid natural element method（HNEM） for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements and incremental stresses. The incremental form of Hellinger–Reissner variational principle for elastoplastic large deformation problems is deduced to obtain the equation system. The total Lagrangian formulation is used to describe the discrete equation system.Compared with the natural element method（NEM）, the HNEM has higher computational precision and efficiency in solving elastoplastic large deformation problems. Some numerical examples are selected to demonstrate the advantage of the HNEM for large deformation elastoplasticity problems.展开更多
This paper proposes a novel numerical solution approach for the kinematic shakedown analysis of strain-hardening thin plates using the C^(1)nodal natural element method(C^(1)nodal NEM).Based on Koiter’s theorem and t...This paper proposes a novel numerical solution approach for the kinematic shakedown analysis of strain-hardening thin plates using the C^(1)nodal natural element method(C^(1)nodal NEM).Based on Koiter’s theorem and the von Mises and two-surface yield criteria,a nonlinear mathematical programming formulation is constructed for the kinematic shakedown analysis of strain-hardening thin plates,and the C^(1)nodal NEM is adopted for discretization.Additionally,König’s theory is used to deal with time integration by treating the generalized plastic strain increment at each load vertex.A direct iterative method is developed to linearize and solve this formulation by modifying the relevant objective function and equality constraints at each iteration.Kinematic shakedown load factors are directly calculated in a monotonically converging manner.Numerical examples validate the accuracy and convergence of the developed method and illustrate the influences of limited and unlimited strain-hardening models on the kinematic shakedown load factors of thin square and circular plates.展开更多
This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yie...This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yield criterion,the nonlinear mathematical programming formulation for upper bound shakedown analysis of thin plates is established.In this formulation,the trail function of residual displacement increment is approximated by using the C^(1) shape functions,the plastic incompressibility condition is satisfied by introducing a constant matrix in the objective function,and the time integration is resolved by using the Konig’s technique.Meanwhile,the objective function is linearized by distinguishing the non-plastic integral points from the plastic integral points and revising the objective function and associated equality constraints at each iteration.Finally,the upper bound shakedown load multipliers of thin plates are obtained by direct iterative and monotone convergence processes.Several benchmark examples verify the good precision and fast convergence of this proposed method.展开更多
By coupling natural boundary element method (NBEM) with FEM based on domain decomposition, the torsion problem of the square cross-sections bar with cracks have been studied, the stresses of the nodes of the cross-sec...By coupling natural boundary element method (NBEM) with FEM based on domain decomposition, the torsion problem of the square cross-sections bar with cracks have been studied, the stresses of the nodes of the cross-sections and the stress intensity factors have been calculated, and some distribution pictures of the stresses have been drawn. During computing, the effect of the relaxed factors to the convergence speed of the iterative method has been discussed. The results of the computation have confirmed the advantages of the NBEM and its coupling with the FEM.展开更多
The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived ...The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.展开更多
基金supported by the SDUST Spring Bud (2009AZZ021)Taian Science and Technology Development (20112001)
文摘C^1 natural element method （C^1 NEM） is applied to strain gradient linear elasticity, and size effects on mi crostructures are analyzed. The shape functions in C^1 NEM are built upon the natural neighbor interpolation （NNI）, with interpolation realized to nodal function and nodal gradient values, so that the essential boundary conditions （EBCs） can be imposed directly in a Galerkin scheme for partial differential equations （PDEs）. In the present paper, C^1 NEM for strain gradient linear elasticity is constructed, and sev- eral typical examples which have analytical solutions are presented to illustrate the effectiveness of the constructed method. In its application to microstructures, the size effects of bending stiffness and stress concentration factor （SCF） are studied for microspeciem and microgripper, respectively. It is observed that the size effects become rather strong when the width of spring for microgripper, the radius of circular perforation and the long axis of elliptical perforation for microspeciem come close to the material characteristic length scales. For the U-shaped notch, the size effects decline obviously with increasing notch radius, and decline mildly with increasing length of notch.
基金Project supported by the Natural Science Foundation of Shanghai,China(Grant No.13ZR1415900)
文摘A hybrid natural element method（HNEM） for two-dimensional viscoelasticity problems under the creep condition is proposed. The natural neighbor interpolation is used as the test function, and the discrete equation system of the HNEM for viscoelasticity problems is obtained using the Hellinger–Reissner variational principle. In contrast to the natural element method（NEM）, the HNEM can directly obtain the nodal stresses, which have higher precisions than those obtained using the moving least-square（MLS） approximation. Some numerical examples are given to demonstrate the validity and superiority of this HNEM.
基金supported by the Natural Science Foundation of Shanghai,China(Grant No.13ZR1415900)
文摘We present the hybrid natural element method（HNEM） for two-dimensional elastoplastic large deformation problems. Sibson interpolation is adopted to construct the shape functions of nodal incremental displacements and incremental stresses. The incremental form of Hellinger–Reissner variational principle for elastoplastic large deformation problems is deduced to obtain the equation system. The total Lagrangian formulation is used to describe the discrete equation system.Compared with the natural element method（NEM）, the HNEM has higher computational precision and efficiency in solving elastoplastic large deformation problems. Some numerical examples are selected to demonstrate the advantage of the HNEM for large deformation elastoplasticity problems.
基金supported by the Chinese Postdoctoral Science Foundation(2013M540934).
文摘This paper proposes a novel numerical solution approach for the kinematic shakedown analysis of strain-hardening thin plates using the C^(1)nodal natural element method(C^(1)nodal NEM).Based on Koiter’s theorem and the von Mises and two-surface yield criteria,a nonlinear mathematical programming formulation is constructed for the kinematic shakedown analysis of strain-hardening thin plates,and the C^(1)nodal NEM is adopted for discretization.Additionally,König’s theory is used to deal with time integration by treating the generalized plastic strain increment at each load vertex.A direct iterative method is developed to linearize and solve this formulation by modifying the relevant objective function and equality constraints at each iteration.Kinematic shakedown load factors are directly calculated in a monotonically converging manner.Numerical examples validate the accuracy and convergence of the developed method and illustrate the influences of limited and unlimited strain-hardening models on the kinematic shakedown load factors of thin square and circular plates.
基金supported by the Chinese Postdoctoral Science Foundation(2013M540934)supported by the National Key Research and Development Program of China(2016YFC0801905,2017YFF0210704).
文摘This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yield criterion,the nonlinear mathematical programming formulation for upper bound shakedown analysis of thin plates is established.In this formulation,the trail function of residual displacement increment is approximated by using the C^(1) shape functions,the plastic incompressibility condition is satisfied by introducing a constant matrix in the objective function,and the time integration is resolved by using the Konig’s technique.Meanwhile,the objective function is linearized by distinguishing the non-plastic integral points from the plastic integral points and revising the objective function and associated equality constraints at each iteration.Finally,the upper bound shakedown load multipliers of thin plates are obtained by direct iterative and monotone convergence processes.Several benchmark examples verify the good precision and fast convergence of this proposed method.
基金the State Key Laboratory of Science and Engineering Computation
文摘By coupling natural boundary element method (NBEM) with FEM based on domain decomposition, the torsion problem of the square cross-sections bar with cracks have been studied, the stresses of the nodes of the cross-sections and the stress intensity factors have been calculated, and some distribution pictures of the stresses have been drawn. During computing, the effect of the relaxed factors to the convergence speed of the iterative method has been discussed. The results of the computation have confirmed the advantages of the NBEM and its coupling with the FEM.
基金subsidized by the National Basic Research Program of China under the grant 2005CB321701the National Natural Science Foundation of China under the grant 10531080the Beijing Natural Science Foundation under the grant 1072009 and the Research Project of Zhejiang Ocean University (X08M013,X08Z04)
文摘The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.
