The standard formula for geometric stiffness matrix calculation, which is convenient for most engineering applications, is seen to be unsatisfactory for large strains because of poor accuracy, low convergence rate, an...The standard formula for geometric stiffness matrix calculation, which is convenient for most engineering applications, is seen to be unsatisfactory for large strains because of poor accuracy, low convergence rate, and stability. For very large compressions, the tangent stiffness in the direction of the compression can even become negative, which can be regarded as physical nonsense. So in many cases rubber materials exposed to great compression cannot be analyzed, or the analysis could lead to very poor convergence. Problems with the standard geometric stiffness matrix can even occur with a small strain in the case of plastic yielding, which eventuates even greater practical problems. The authors demonstrate that amore precisional approach would not lead to such strange and theoretically unjustified results. An improved formula that would eliminate the disadvantages mentioned above and leads to higher convergence rate and more robust computations is suggested in this paper. The new formula can be derived from the principle of virtual work using a modified Green-Lagrange strain tensor, or from equilibrium conditions where in the choice of a specific strain measure is not needed for the geometric stiffness derivation (which can also be used for derivation of geometric stiffness of a rigid truss member). The new formula has been verified in practice with many calculations and implemented in the RFEM and SCIA Engineer programs. The advantages of the new formula in comparison with the standard formula are shown using several examples.展开更多
We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the ca...We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the case of geometrically localized stiffness,where the convergence rates of iterative methods are degraded in a localized subregion of the mesh.Slower convergence may be caused by a number of factors,including the mesh size,anisotropy,highly variable coefficients,and more challenging physics.The approach taken in this work is to correct well-known preconditioners such as the block Jacobi and the block incomplete LU(ILU)with an adaptive inner subregion iteration.The goal of these preconditioners is to reduce the number of costly global iterations by accelerating the convergence in the stiff region by iterating on the less expensive reduced problem.The tolerance for the inner iteration is adaptively chosen to minimize subregion-local work while guaranteeing global convergence rates.We present analysis showing that the convergence of these preconditioners,even when combined with an adaptively selected tolerance,is independent of discretization parameters(e.g.,the mesh size and diffusion coefficient)in the subregion.We demonstrate significant performance improvements over black-box preconditioners when applied to several model convection-diffusion problems.Finally,we present performance results of several variations of iterative subregion correction preconditioners applied to the Reynolds number 2.25×10^(6)fluid flow over the NACA 0012 airfoil,as well as massively separated flow at 30°angle of attack.展开更多
文摘The standard formula for geometric stiffness matrix calculation, which is convenient for most engineering applications, is seen to be unsatisfactory for large strains because of poor accuracy, low convergence rate, and stability. For very large compressions, the tangent stiffness in the direction of the compression can even become negative, which can be regarded as physical nonsense. So in many cases rubber materials exposed to great compression cannot be analyzed, or the analysis could lead to very poor convergence. Problems with the standard geometric stiffness matrix can even occur with a small strain in the case of plastic yielding, which eventuates even greater practical problems. The authors demonstrate that amore precisional approach would not lead to such strange and theoretically unjustified results. An improved formula that would eliminate the disadvantages mentioned above and leads to higher convergence rate and more robust computations is suggested in this paper. The new formula can be derived from the principle of virtual work using a modified Green-Lagrange strain tensor, or from equilibrium conditions where in the choice of a specific strain measure is not needed for the geometric stiffness derivation (which can also be used for derivation of geometric stiffness of a rigid truss member). The new formula has been verified in practice with many calculations and implemented in the RFEM and SCIA Engineer programs. The advantages of the new formula in comparison with the standard formula are shown using several examples.
文摘We present a class of preconditioners for the linear systems resulting from a finite element or discontinuous Galerkin discretizations of advection-dominated problems.These preconditioners are designed to treat the case of geometrically localized stiffness,where the convergence rates of iterative methods are degraded in a localized subregion of the mesh.Slower convergence may be caused by a number of factors,including the mesh size,anisotropy,highly variable coefficients,and more challenging physics.The approach taken in this work is to correct well-known preconditioners such as the block Jacobi and the block incomplete LU(ILU)with an adaptive inner subregion iteration.The goal of these preconditioners is to reduce the number of costly global iterations by accelerating the convergence in the stiff region by iterating on the less expensive reduced problem.The tolerance for the inner iteration is adaptively chosen to minimize subregion-local work while guaranteeing global convergence rates.We present analysis showing that the convergence of these preconditioners,even when combined with an adaptively selected tolerance,is independent of discretization parameters(e.g.,the mesh size and diffusion coefficient)in the subregion.We demonstrate significant performance improvements over black-box preconditioners when applied to several model convection-diffusion problems.Finally,we present performance results of several variations of iterative subregion correction preconditioners applied to the Reynolds number 2.25×10^(6)fluid flow over the NACA 0012 airfoil,as well as massively separated flow at 30°angle of attack.
