References of "Bordas, Stéphane 50000969"
     in
Bookmark and Share    
Full Text
Peer Reviewed
See detailCalculating the Malliavin derivative of some stochastic mechanics problems
Hauseux, Paul UL; Hale, Jack UL; Bordas, Stéphane UL

in PLoS ONE (2017), 12(12), 0189994

The Malliavin calculus is an extension of the classical calculus of variations from deterministic functions to stochastic processes. In this paper we aim to show in a practical and didactic way how to ... [more ▼]

The Malliavin calculus is an extension of the classical calculus of variations from deterministic functions to stochastic processes. In this paper we aim to show in a practical and didactic way how to calculate the Malliavin derivative, the derivative of the expectation of a quantity of interest of a model with respect to its underlying stochastic parameters, for four problems found in mechanics. The non-intrusive approach uses the Malliavin Weight Sampling (MWS) method in conjunction with a standard Monte Carlo method. The models are expressed as ODEs or PDEs and discretised using the finite difference or finite element methods. Specifically, we consider stochastic extensions of; a 1D Kelvin-Voigt viscoelastic model discretised with finite differences, a 1D linear elastic bar, a hyperelastic bar undergoing buckling, and incompressible Navier-Stokes flow around a cylinder, all discretised with finite elements. A further contribution of this paper is an extension of the MWS method to the more difficult case of non-Gaussian random variables and the calculation of second-order derivatives. We provide open-source code for the numerical examples in this paper. [less ▲]

Detailed reference viewed: 291 (39 UL)
Full Text
Peer Reviewed
See detailMicro-structured materials: inhomogeneities and imperfect interfaces in plane micropolar elasticity, a boundary element approach
Atroshchenko, Elena; Hale, Jack UL; Videla, Javier A. et al

in Engineering Analysis with Boundary Elements (2017), 83

In this paper we tackle the simulation of microstructured materials modelled as heterogeneous Cosserat media with both perfect and imperfect interfaces. We formulate a boundary value problem for an ... [more ▼]

In this paper we tackle the simulation of microstructured materials modelled as heterogeneous Cosserat media with both perfect and imperfect interfaces. We formulate a boundary value problem for an inclusion of one plane strain micropolar phase into another micropolar phase and reduce the problem to a system of boundary integral equations, which is subsequently solved by the boundary element method. The inclusion interface condition is assumed to be imperfect, which permits jumps in both displacements/microrotations and tractions/couple tractions, as well as a linear dependence of jumps in displacements/microrotations on continuous across the interface tractions/couple traction (model known in elasticity as homogeneously imperfect interface). These features can be directly incorporated into the boundary element formulation. The BEM-results for a circular inclusion in an in finite plate are shown to be in excellent agreement with the analytical solutions. The BEM-results for inclusions in finite plates are compared with the FEM-results obtained with FEniCS. [less ▲]

Detailed reference viewed: 273 (14 UL)
Full Text
Peer Reviewed
See detailIsogeometric analysis of thin Reissner-Mindlin plates and shells: Locking phenomena and generalized local B-bar method
Hu, Qingyuan UL; Xia, Yang; Natarajan, Sundararajan et al

E-print/Working paper (2017)

We propose a generalized local $\bar{B}$ framework, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are ... [more ▼]

We propose a generalized local $\bar{B}$ framework, addressing locking in degenerated Reissner-Mindlin plate and shell formulations in the context of isogeometric analysis. Parasitic strain components are projected onto the physical space locally, i.e. at the element level, using a least-squares approach. The formulation is general and allows the flexible utilization of basis functions of different order as the projection bases. The present formulation is much cheaper computationally than the global $\bar{B}$ method. Through numerical examples, we show the consistency of the scheme, although the method is not Hu-Washizu variationally consistent. The numerical examples show that the proposed formulation alleviates locking and yields good accuracy for various thicknesses, even for slenderness ratios of $1 \times 10^5$, and has the ability to capture deformations of thin shells using relatively coarse meshes. From the detailed numerical study, it can be opined that the proposed method is less sensitive to locking and mesh distortion. [less ▲]

Detailed reference viewed: 248 (15 UL)
Full Text
Peer Reviewed
See detailBayesian inference to identify parameters in viscoelasticity
Rappel, Hussein UL; Beex, Lars UL; Bordas, Stéphane UL

in Mechanics of Time-Dependent Materials (2017)

This contribution discusses Bayesian inference (BI) as an approach to identify parameters in viscoelasticity. The aims are: (i) to show that the prior has a substantial influence for viscoelasticity, (ii ... [more ▼]

This contribution discusses Bayesian inference (BI) as an approach to identify parameters in viscoelasticity. The aims are: (i) to show that the prior has a substantial influence for viscoelasticity, (ii) to show that this influence decreases for an increasing number of measurements and (iii) to show how different types of experiments influence the identified parameters and their uncertainties. The standard linear solid model is the material description of interest and a relaxation test, a constant strain-rate test and a creep test are the tensile experiments focused on. The experimental data are artificially created, allowing us to make a one-to-one comparison between the input parameters and the identified parameter values. Besides dealing with the aforementioned issues, we believe that this contribution forms a comprehensible start for those interested in applying BI in viscoelasticity. [less ▲]

Detailed reference viewed: 535 (177 UL)
Full Text
Peer Reviewed
See detailDiscretisation and Model Selection for Interface Problems in Mechanics
Bordas, Stéphane UL

in International Journal of Computational Methods (2017, August 04)

Detailed reference viewed: 98 (7 UL)
Full Text
See detailExperimental and numerical assessment of the mechanics of keloid-skin composites undergoing large deformations
Sensale, Marco UL; Chambert, Jerome; Chouly, Franz et al

Scientific Conference (2017, August)

Detailed reference viewed: 119 (14 UL)
Full Text
Peer Reviewed
See detailError-controlled computational fracture mechanics
Bordas, Stéphane UL

Scientific Conference (2017, July 12)

Detailed reference viewed: 375 (33 UL)
Full Text
See detailTowards a seamless Integration of CAD and Simulation: CISM Course 2017
Bordas, Stéphane UL

Learning material (2017)

Isogeometric analysis relies on the use of the same basis functions as employed in Computer Aided Design (CAD). This offers the possibility to facilitate design and optimisation. The previous course ... [more ▼]

Isogeometric analysis relies on the use of the same basis functions as employed in Computer Aided Design (CAD). This offers the possibility to facilitate design and optimisation. The previous course “Isogeometric methods for numerical simulation” held in 2013 had the aim to give an introduction to isogeometric analysis, its advantages, drawbacks and to the range of its applications. The aim of the proposed new course will be different. The focus will be more on the connection of simulation to CAD systems and how CAD data can be used directly for simulation, leading to a seamless integration. An overview of recent advances and applications will be also presented. The course will start with an introduction to NURBS and their use in describing geometry and in simulation. This will be followed by lectures from a CAD vendor describing the current state of development. Currently available connections to simulation software will also be discussed. Next the use of NURBS for 3D structural analysis, structural optimisation and damage tolerance assessment will be presented, including such advanced topics as the treatment of discontinuities and real-time solvers. It will also be discussed when it might be advantageous to decouple the boundary discretisation from the field variable discretisation, in particular in shape optimisation. Isogeometric methods for the analysis of beam and shell structures, including shape optimisation and fluid structure interaction, will be presented. Lectures on the mathematical and algorithmic foundations of analysis-suitable geometry will follow. This includes an introduction to T-splines and multilevel spline schemes such as hierarchical B- splines. Common analysis-suitable spline algorithms will be presented in the context of Bézier extraction and projection as well as its application as a foundation for integrated engineering design and analysis. An important aspect of analysis-suitable geometry is the ability to locally adapt the smooth spline basis. Several common refinement algorithms will be reviewed as well as their application in several demanding areas of application. The emerging area of weak geometry will be introduced as well as its application to the rapid construction of complex structural assemblies. With the rapid development of isogeometric analysis in recent years, there is an urgent need for volumetric parameterization such as volumetric T-spline model construction. Several volumetric T- spline modeling techniques, that were developed in recent years will be presented. They include converting any quad/ hex meshes to standard and rational T-splines, polycube-based parametric mapping, feature preservation using eigenfunctions, Boolean operations and skeletons, truncated hierarchical Catmull-Clark subdivision, weighted T-splines, conformal T-spline modeling, as well as incorporating T-splines into commercial CAD and FEA software, will be presented. The target audience will be engineers, interested in simulation, software developers and researchers. [less ▲]

Detailed reference viewed: 852 (21 UL)
Full Text
Peer Reviewed
See detailThree-dimensional remeshed smoothed particle hydrodynamics for the simulation of isotropic turbulence
Obeidat, Anas UL; Bordas, Stéphane UL

in International Journal for Numerical Methods in Fluids (2017)

We present a remeshed particle-mesh method for the simulation of three-dimensional compressible turbulent flow. The method is related to the mesh free smoothed particle hydrodynamic (SPH) method, but the ... [more ▼]

We present a remeshed particle-mesh method for the simulation of three-dimensional compressible turbulent flow. The method is related to the mesh free smoothed particle hydrodynamic (SPH) method, but the present method introduces a mesh for efficient calculation of the pressure gradient, and laminar and turbulent diffusion. In addition, the mesh is used to remesh (reorganise uniformly) the particles to ensure a regular particle distribution and convergence of the method. The accuracy of the presented methodology is tested for a number of benchmark problems involving two- and three-dimensional Taylor-Green flow, thin double shear layer, and three-dimensional isotropic turbulence. Two models were implemented, direct numerical simulations, and Smagorinsky model. Taking advantage of the Lagrangian advection, and the finite difference efficiency, the method is capable of providing quality simulations while maintaining its robustness and versatility [less ▲]

Detailed reference viewed: 177 (10 UL)
Full Text
Peer Reviewed
See detailModelling hydraulic fractures in porous media using flow cohesive interface elements
Nguyen, Vinh Phu; Lian, Haojie; Rabczuk, Timon et al

in Engineering Geology (2017), 225

Detailed reference viewed: 141 (2 UL)
Full Text
See detailWhat makes Data Science different? A discussion involving Statistics2.0 and Computational Sciences
Ley, Christophe; Bordas, Stéphane UL

E-print/Working paper (2017)

Data Science is today one of the main buzzwords be it in business, industrial or academic settings. Machine learning, experimental design, data-driven modelling are all, undoubtedly, rising disciplines if ... [more ▼]

Data Science is today one of the main buzzwords be it in business, industrial or academic settings. Machine learning, experimental design, data-driven modelling are all, undoubtedly, rising disciplines if one goes by the soaring number of research papers and patents appearing each year. The prospect of becoming a ``Data Scientist'' appeals to many. A discussion panel organised as part of the European Data Science Conference (European Association for Data Science (EuADS)) asked the question: ``What makes Data Science different?'' In this paper we give our own, personal and multi-facetted view on this question, from an engineering and a statistics perspective. In particular, we compare Data Science to Statistics and discuss the connection between Data Science and Computational Science. [less ▲]

Detailed reference viewed: 981 (36 UL)
Full Text
Peer Reviewed
See detailReduced basis Nitsche-based domain decomposition: a biomedical application
Baroli, Davide UL; Beex, Lars UL; Hale, Jack UL et al

Scientific Conference (2017, March 10)

Nowadays, the personalized biomedical simulations demand real-time efficient and reliable method to alleviate the computational complexity of high-fidelity simulation. In such applications, the necessity ... [more ▼]

Nowadays, the personalized biomedical simulations demand real-time efficient and reliable method to alleviate the computational complexity of high-fidelity simulation. In such applications, the necessity of solving different substructure, e.g. tissues or organs, with different numbers of the degrees of freedom and of coupling the reduced order spaces for each substructure poses a challenge in the on-fly simulation. In this talk, this challenge is taken into account employing the Nitsche-based domain decomposition technique inside the reduced order model [1]. This technique with respect to other domain decomposition approach allows obtaining a solution with the same accuracy of underlying finite element formulation and to flexibly treat interface with non-matching mesh. The robustness of the coupling is determined by the penalty coefficients that is chosen using ghost penalty technique [2]. Furthermore, to reduce the computational complexity of the on-fly assembling it is employed the empirical interpolation approach proposed in [3]. The numerical tests, performed using FEniCS[4], petsc4py and slepc4py [5], shows the good performance of the method and the reduction of computation cost. [1] Baroli, D., Beex L. and Bordas, S. Reduced basis Nitsche-based domain decomposition. In preparation. [2] Burman, E., Claus, S., Hansbo, P., Larson, M. G., & Massing, A. (2015). CutFEM: Discretizing geometry and partial differential equations. International Journal for Numerical Methods in Engineering, 104(7), 472-501. [3] E. Schenone, E., Beex,L., Hale, J.S., Bordas S. Proper Orthogonal Decomposition with reduced integration method. Application to nonlinear problems. In preparation. [4] A. Logg, K.-A. Mardal, G. N. Wells et al. Automated Solution of Differential Equations by the Finite Element Method, Springer 2012. [5] L. Dalcin, P. Kler, R. Paz, and A. Cosimo, Parallel Distributed Computing using Python, Advances in Water Resources, 34(9):1124-1139, 2011. http://dx.doi.org/10.1016/j.advwatres.2011.04.013 [less ▲]

Detailed reference viewed: 282 (10 UL)
Full Text
Peer Reviewed
See detailStable 3D XFEM/vector-level sets for non-planar 3D crack propagation and comparison of enrichment schemes
Agathos, Konstantinos UL; Ventura, Giulio; Chatzi, Eleni et al

in International Journal for Numerical Methods in Engineering (2017)

We present a three-dimensional (3D) vector level set method coupled to a recently developed stable extended finite element method (XFEM). We further investigate a new enrichment approach for XFEM adopting ... [more ▼]

We present a three-dimensional (3D) vector level set method coupled to a recently developed stable extended finite element method (XFEM). We further investigate a new enrichment approach for XFEM adopting discontinuous linear enrichment functions in place of the asymptotic near-tip functions. Through the vector level set method, level set values for propagating cracks are obtained via simple geometrical operations, eliminating the need for solution of differential evolution equations. The first XFEM variant ensures optimal convergence rates by means of geometrical enrichment, i.e., the use of enriched elements in a fixed volume around the crack front, without giving rise to conditioning problems. The linear enrichment approach significantly simplifies implementation and reduces the computational cost associated with numerical integration. The two dicretization schemes are tested for different benchmark problems, and their combination to the vector level set method is verified for non-planar crack propagation problems. [less ▲]

Detailed reference viewed: 360 (28 UL)
Full Text
Peer Reviewed
See detailTrefftz polygonal finite element for linear elasticity: convergence, accuracy, and properties
Hirshikesh; Natarajan, Sundararajan; Ratna Kumar, A. K. et al

in Asia Pacific Journal on Computational Engineering (2017)

In this paper, the accuracy and the convergence properties of Trefftz finite element method over arbitrary polygons are studied. Within this approach, the unknown displacement field within the polygon is ... [more ▼]

In this paper, the accuracy and the convergence properties of Trefftz finite element method over arbitrary polygons are studied. Within this approach, the unknown displacement field within the polygon is represented by the homogeneous solution to the governing differential equations, also called as the T-complete set. While on the boundary of the polygon, a conforming displacement field is independently defined to enforce the continuity of the field variables across the element boundary. An optimal number of T-complete functions are chosen based on the number of nodes of the polygon and the degrees of freedom per node. The stiffness matrix is computed by the hybrid formulation with auxiliary displacement frame. Results from the numerical studies presented for a few benchmark problems in the context of linear elasticity show that the proposed method yields highly accurate results with optimal convergence rates. [less ▲]

Detailed reference viewed: 102 (1 UL)
Full Text
Peer Reviewed
See detailA linear smoothed higher-order CS-FEM for the analysis of notched laminated composites
Wan, Detao; Hu, Dean; Natarajan, Sundararajan et al

in Engineering Analysis with Boundary Elements (2017), 85

Higher-order elements with highly accurate solutions are attractive for stress analysis and stress concentration problems. However, the distorted eight-node serendipity quadrilateral element is known to ... [more ▼]

Higher-order elements with highly accurate solutions are attractive for stress analysis and stress concentration problems. However, the distorted eight-node serendipity quadrilateral element is known to yield inaccurate re- sults and sub-optimal convergence rate. In this paper, we present a higher order CS-FEM to alleviate the effect of distorted mesh and guarantee the quality of solutions by employing a linear smoothing technique over eight-node quadratic serendipity elements. The modified. strain matrix is computed by the divergence theorem between the nodal shape functions and their derivatives using Taylor’s expansion of the weak form. The proposed method eliminates the need for isoparametric mapping and numerical studies demonstrate that the proposed method is insensitive to mesh distortion. The improved accuracy and superior convergence rates are numerically demon- strated with a few benchmark problems. The analysis of the stress concentration around cutouts also proves that the present method has good performance for the laminated composites. [less ▲]

Detailed reference viewed: 129 (1 UL)
Full Text
Peer Reviewed
See detailA fully smoothed XFEM for analysis of axisymmetric problems with weak discontinuities
Wan, Detao; Hu, Dean; Natarajan, Sundararajan et al

in International Journal for Numerical Methods in Engineering (2017), 110(3), 203-226

In this paper, we propose a fully smoothed extended finite element method (SmXFEM) for axisymmetric problems with weak discontinuities. The salient feature of the proposed approach is that all the terms ... [more ▼]

In this paper, we propose a fully smoothed extended finite element method (SmXFEM) for axisymmetric problems with weak discontinuities. The salient feature of the proposed approach is that all the terms in the stiffness and mass matrixes can be computed by smoothing technique. This is accomplished by combining the Green’s divergence theorem with the evaluation of indefinite integral based on smoothing technique, which is used to transform the domain integral into boundary integral. The proposed technique completely eliminates the need for isoparametric mapping and the computing of Jacobian matrix even for the mass matrix. When employed over the enriched elements, the proposed technique does not require sub-triangulation for the purpose of numerical integration. The accuracy and convergence properties of the proposed technique are demonstrated with a few problems in elastostatics and elastodynamics with weak discontinuities. It can be seen that the proposed technique yields stable and accurate solutions and is less sensitive to mesh distortion. [less ▲]

Detailed reference viewed: 133 (2 UL)
Full Text
Peer Reviewed
See detailNumerical evaluation of buckling behaviour induced by compression on patch-repaired composites
Deng, Jian; Zhou, Guangming; Bordas, Stéphane UL et al

in Composite Structures (2017), 168

A progressive damage model is proposed to predict buckling strengths and failure mechanisms for both symmetric and asymmetric patch repaired carbon-fibre reinforced laminates subjected to compression ... [more ▼]

A progressive damage model is proposed to predict buckling strengths and failure mechanisms for both symmetric and asymmetric patch repaired carbon-fibre reinforced laminates subjected to compression without lateral restrains. Solid and cohesive elements are employed to discretize composite and adhesive layers, respectively. Coupling with three dimensional strain failure criteria, an energy-based crack band model is applied to address the softening behaviour in composites with mesh dependency elimination. Both laminar and laminate scaled failure are addressed. Patch debonding is simulated by the cohesive zone model with a trapezoidal traction–separation law applied for the ductile adhesive. Geometric imperfection is introduced into the nonlinear analysis by the first order linear buckling configuration. Regarding strengths and failure patterns, the simulation demonstrates an accurate and consistent prediction compared with experimental observations. Though shearing is the main contributor to damage initiation in adhesive, stress analysis shows that lateral deformation subsequently reverses the distribution of normal stresses which stimulates patch debonding at one of the repair sides. The influence of patch dimensions on strengths and failure mechanisms can be explained by stress distributions in adhesive and lateral deformation of repairs. Comparison between symmetric and asymmetric regarding strength and failure modes shows that structural asymmetry can intensify lateral flexibility. This resulted in earlier patch debonding and negative effects on strengths. [less ▲]

Detailed reference viewed: 112 (2 UL)
Full Text
Peer Reviewed
See detailAn implicit potential method along with a meshless technique for incompressible fluid flows for regular and irregular geometries in 2D and 3D
Bourantas, Georgios UL; Loukopoulos, V. C.; Chowdhury, H. A. et al

in Engineering Analysis with Boundary Elements (2017), 77

We present the Implicit Potential (IPOT) numerical scheme developed in the framework of meshless point collocation. The proposed scheme is used for the numerical solution of the steady state ... [more ▼]

We present the Implicit Potential (IPOT) numerical scheme developed in the framework of meshless point collocation. The proposed scheme is used for the numerical solution of the steady state, incompressible Navier-Stokes (N-S) equations in their primitive variable (u-v-w-p) formulation. The governing equations are solved in their strong form using either a collocated or a semi-staggered type meshless nodal configuration. The unknown field functions and derivatives are calculated using the Modified Moving Least Squares (MMLS) interpolation method. Both velocity-correction and pressure correction methods applied ensure the incompressibility constraint and mass conservation. The proposed meshless point collocation (MPC) scheme has the following characteristics: (i) it can be applied, in a straightforward manner to: steady, unsteady, internal and external fluid flows in 2D and 3D, (ii) it equally applies to regular an irregular geometries, (iii) a distribution of points is sufficient, no numerical integration in space nor any mesh structure are required, (iv) there is no need for pressure boundary conditions since no pressure constitutive equation is solved, (v) it is quite simple and accurate, (vi) results can be obtained using collocated or semi-staggered nodal distributions, (vii) there is no need to compute the velocity potential nor the unit normal vectors and (viii) there is no need for a curvilinear system of coordinates. Simulations of fluid flow in 2D and 3D for regular and irregular geometries indicate the validity of the proposed methodology. [less ▲]

Detailed reference viewed: 135 (2 UL)
Full Text
Peer Reviewed
See detailReal-time Error Control for Surgical Simulation
Phuoc Bui, Huu; Tomar, Satyendra; Courtecuisse, Hadrien et al

in IEEE Transactions on Biomedical Engineering (2017)

To present the first real-time a posteriori error-driven adaptive finite element approach for realtime simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational ... [more ▼]

To present the first real-time a posteriori error-driven adaptive finite element approach for realtime simulation and to demonstrate the method on a needle insertion problem. Methods: We use corotational elasticity and a frictional needle/tissue interaction model. The problem is solved using finite elements within SOFA. The refinement strategy relies upon a hexahedron-based finite element method, combined with a posteriori error estimation driven local h-refinement, for simulating soft tissue deformation. Results: We control the local and global error level in the mechanical fields (e.g. displacement or stresses) during the simulation. We show the convergence of the algorithm on academic examples, and demonstrate its practical usability on a percutaneous procedure involving needle insertion in a liver. For the latter case, we compare the force displacement curves obtained from the proposed adaptive algorithm with that obtained from a uniform refinement approach. Conclusions: Error control guarantees that a tolerable error level is not exceeded during the simulations. Local mesh refinement accelerates simulations. Significance: Our work provides a first step to discriminate between discretization error and modeling error by providing a robust quantification of discretization error during simulations. [less ▲]

Detailed reference viewed: 205 (6 UL)