Numerical methods for fracture/cutting of heterogeneous materials

Presentation (2016, December)

Shape Optimization Directly from CAD: an Isogeometric Boundary Element Approach Using T-splines

Report (2016)

Bayesian inference for material parameter identification in elastoplasticity

Scientific Conference (2016, September 07)

Multi-scale modelling of fracture

Speeches/Talks (2016)

We present recent models on complexity reduction for computational fracture mechanics

Linear elastic fracture simulation directly from CAD: 2D NURBS-based implementation and role of tip enrichment

in International Journal of Fracture (2016)

A Bayesian approach for parameter identification in elastoplasticity

Scientific Conference (2016, June 09)

Weakening the tight coupling between geometry and simulation in isogeometric analysis

Presentation (2016, June 07)

In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry ... [more ▼]

In the standard paradigm of isogeometric analysis, the geometry and the simulation spaces are tightly integrated, i.e. the same non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is employed for the numerical solution of the problem over the domain. However, there are situations where this tight integration is a bane rather than a boon. Such situations arise where, e.g., (1) the geometry of the domain is simple enough to be represented by low order NURBS, whereas the unknown (exact) solution of the problem is sufficiently regular, and thus, the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry, (2) the constraint of using the same space for the geometry and the numerical solution is particularly undesirable, such as in the shape and topology optimization, and (3) the solution of the problem has low regularity but for the curved boundary of the domain one can employ higher order NURBS. Therefore, we propose to weaken this constraint. An extensive study of patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution will be discussed. It will be shown, with concrete reasoning, that why patch test fails in certain cases, and that those cases should be avoided in practice. Thereafter, selective numerical examples will be presented to address some of the above-mentioned situations, and it will be shown that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces, and thus, improved accuracy of the numerical solution. [less ▲]

Generalizing the isogeometric concept: weakening the tight coupling between geometry and simulation in IGA

Presentation (2016, June 02)

In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the ... [more ▼]

In the standard paradigm of isogeometric analysis [2, 1], the geometry and the simulation spaces are tightly integrated, i.e. the non-uniform rational B-splines (NURBS) space, which is used for the geometry representation of the domain, is also employed for the numerical solution of the problem over the domain. However, in certain situations, such as, when the geometry of the domain can be represented by low order NURBS but the numerical solution can be obtained with improved accuracy by using NURBS of order higher than that required for the geometry; or in the shape and topology optimization where the constraint of using the same space for the geometry and the numerical solution is not favorable, this tight coupling is disadvantageous. Therefore, we study the effect of decoupling the spaces for the geometry representation and the numerical solution, though still using the prevalent functions in CAD/CAGD. To begin with, we perform the patch tests on various combinations of polynomial degree, geometry type, and various cases of varying degrees and control variables between the geometry and the numerical solution. This shows that certain cases, perhaps intuitive, should be avoided in practice because patch test fails. The above-mentioned situations are further explored with some numerical examples, which shows that weakening the tight coupling between geometry and simulation offers more flexibility in choosing the numerical solution spaces. [1] J. Cottrell, T.J.R. Hughes, and Y. Bazilevs. Isogeometric Analysis: Toward Integration of CAD and FEA, volume 80. Wiley, Chichester, 2009. [2] T.J.R. Hughes, J. Cottrell, and Y. Bazilevs. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement. Computer Methods in Applied Mechanics and Engineering, 194:4135–4195, 2005. [less ▲]

Efficient propagation of uncertainty through an inverse non-linear deformation model of soft tissue

Scientific Conference (2016, June)

Linear smoothing over arbitrary polytopes for compressible and nearly incompressible linear elasticity

Scientific Conference (2016, June)

We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. In this approach, a volume-averaged nodal projection operator is ... [more ▼]

We present a displacement based approach over arbitrary polytopes for compressible and nearly incompressible linear elastic solids. In this approach, a volume-averaged nodal projection operator is constructed to project the dilatational strain into an approximation space of equal or lower-order than the approximation space for the displacement field, resulting in a locking-free method. The formulation uses the usual Wachspress interpolants over arbitrary polytopes and the stability of the method is ensured by the addition of bubble like functions. The smoothed strains are evaluated based on the linear smoothing procedure. This further softens the bilinear form allowing the procedure to search for a solution satisfying the divergence- free condition. The divergence-free condition of the proposed approach is verified through systematic numerical study. The formulation delivers optimal convergence rates in the energy and L2-norms. Inf-sup tests are presented to demonstrated the stability of the formulation. [less ▲]