A quadrilateral element with smoothed curvatures for Mindlin-Reissner plates is proposed. The curvature at each point is obtained by a non-local approximation via a smoothing function. The bending stiffness matrix is calculated by a boundary integral along the boundaries of the smoothing elements (smoothing cells). Numerical results show that the proposed element is robust, computational inexpensive and simultaneously very accurate and free of locking, even for very thin plates. The most promising feature of our elements is their insentivity to mesh distortion.

This article describes a bridge between POD-based model order reduction techniques and the classical Newton/Krylov solvers. This bridge is used to derive an efficient algorithm to correct, "on-the-fly", the reduced order modelling of highly nonlinear problems undergoing strong topological changes. Damage initiation problems are addressed and tackle via a corrected hyperreduction method. It is shown that the relevancy of reduced order model can be significantly improved with reasonable additional costs when using this algorithm, even when strong topological changes are involved.

No abstract

SUMMARYThis paper proposes a novel technique to reduce the computational burden associated with the simulation of localized failure. The proposed methodology affords the simulation of damage initiation and propagation while concentrating the computational effort where it is most needed, that is, in the localization zones. To do so, a local/global technique is devised where the global (slave) problem (far from the zones undergoing severe damage and cracking) is solved for in a reduced space computed by the classical proper orthogonal decomposition while the local (master) degrees of freedom (associated with the part of the structure where most of the damage is taking place) are fully resolved. Both domains are coupled through a local/global technique. This method circumvents the difficulties associated with model order reduction for the simulation of highly nonlinear mechanical failure and offers an alternative or complementary approach to the development of multiscale fracture simulators.

The concept of isogeometric analysis, whereby the parametric functions that are used to describe CAD geometry are also used to approximate the unknown fields in a numerical discretisation, has progressed rapidly in recent years. This paper advances the field further by outlining an isogeometric Boundary Element Method (IGABEM) that only requires a representation of the geometry of the domain for analysis, fitting neatly with the boundary representation provided completely by CAD. The method circumvents the requirement to generate a boundary mesh representing a significant step in reducing the gap between engineering design and analysis. The current paper focuses on implementation details of 2D IGABEM for elastostatic analysis with particular attention paid towards the differences over conventional boundary element implementations. Examples of Matlab R code are given whenever possible to aid understanding of the techniques used.