Browse ORBi

- What it is and what it isn't
- Green Road / Gold Road?
- Ready to Publish. Now What?
- How can I support the OA movement?
- Where can I learn more?

ORBi

High-performance modeling of concrete ageing Habera, Michal ; Zilian, Andreas in Proceedings in Applied Mathematics and Mechanics (in press) Long-term behaviour of concrete structural elements is very important for evaluation of its health and serviceability range. The phenomena that must be considered are complex and lead to coupled ... [more ▼] Long-term behaviour of concrete structural elements is very important for evaluation of its health and serviceability range. The phenomena that must be considered are complex and lead to coupled multiphysics formulations. Such formulations are difficult not only from physical perspective, but also from computational perspective. In this contribution attention to computational efficiency and effective implementation is payed. Presented model for concrete ageing is based on microprestress-solidification (MPS) theory of Bazant [1], Kunzel’s model for heat and moisture transport [2] and Mazars model for damage [3]. Ageing linear viscoelastic response, which is immanent to MPS theory and concrete creep, leads to ordinary differetial equation for internal variables solved for every quadrature/nodal point. Numerical structure of the finite element discretisation is examined. Few simplifications on physical model lead to a very efficient linear algebra problem for which standard preconditioned Krylov solvers are reviewed. In parallel, weak and strong scaling tests are performed. All results are produced within open-source finite element framework FEniCS [4]. These models are usually a basis for more involved thermo-hygro-chemo-mechanical (THCM) models with migrating chemical species. It is anticipated, that presented results will help practitioners or other structural engineerers with the choice of suitable and efficient methods for long-term concrete modeling. [less ▲] Detailed reference viewed: 88 (9 UL)Non-localised contact between beams with circular and elliptical cross-sections Magliulo, Marco ; Lengiewicz, Jakub ; Zilian, Andreas et al in Computational Mechanics (in press) The key novelty of this contribution is a dedicated technique to e fficiently determine the distance (gap) function between parallel or almost parallel beams with circular and elliptical cross-sections ... [more ▼] The key novelty of this contribution is a dedicated technique to e fficiently determine the distance (gap) function between parallel or almost parallel beams with circular and elliptical cross-sections. The technique consists of parametrizing the surfaces of the two beams in contact, fixing a point on the centroid line of one of the beams and searching for a constrained minimum distance between the surfaces (two variants are investigated). The resulting unilateral (frictionless) contact condition is then enforced with the Penalty method, which introduces compliance to the, otherwise rigid, beams' cross-sections. Two contact integration schemes are considered: the conventional slave-master approach (which is biased as the contact virtual work is only integrated over the slave surface) and the so-called two-half-pass approach (which is unbiased as the contact virtual work is integrated over the two contacting surfaces). Details of the finite element formulation which is suitably implemented using Automatic Di fferentiation techniques are presented. A set of numerical experiments shows the overall performance of the framework and allows a quantitative comparison of the investigated variants. [less ▲] Detailed reference viewed: 7 (1 UL)Contact between shear-deformable beams with elliptical cross-sections Magliulo, Marco ; Zilian, Andreas ; Beex, Lars in Acta Mechanica (in press) Slender constituents are present in many structures and materials. In associated mechanical models, each slender constituent is often described with a beam. Contact between beams is essential to ... [more ▼] Slender constituents are present in many structures and materials. In associated mechanical models, each slender constituent is often described with a beam. Contact between beams is essential to incorporate in mechanical models, but associated contact frameworks are only demonstrated to work for beams with circular cross-sections. Only two studies have shown the ability to treat contact between beams with elliptical cross-sections, but those frameworks are limited to point-wise contact, which narrows their applicability. This contribution presents initial results of a framework for shear-deformable beams with elliptical cross-sections if contact occurs along a line or at an area (instead of at a point). This is achieved by integrating a penalty potential over one of the beams’ surfaces. Simo-Reissner Geometrically Exact Beam (GEB) elements are employed to discretise each beam. As the surface of an assembly of such beam elements is discontinuous, a smoothed surface is introduced to formulate the contact kinematics. This enables the treatment of contact for large sliding displacements and substantial deformations. [less ▲] Detailed reference viewed: 66 (6 UL)On predictive modelling of yield stress increase in fresh cement paste Huang, Haiqin ; Zilian, Andreas in Proceedings in Applied Mathematics and Mechanics (2019), 19(1), Detailed reference viewed: 82 (2 UL)Simultaneous finite element analysis of circuit-integrated piezoelectric energy harvesting from fluid-structure interaction ; Zilian, Andreas in Mechanical Systems & Signal Processing (2019), 114 Flow-driven piezoelectric energy harvesting is a strongly coupled multiphysics phenomenon that involves complex three-way interaction between the fluid flow, the electromechanical effect of the ... [more ▼] Flow-driven piezoelectric energy harvesting is a strongly coupled multiphysics phenomenon that involves complex three-way interaction between the fluid flow, the electromechanical effect of the piezoelectric material mounted on a deformable substrate structure and the controlling electrical circuit. High fidelity computational solution approaches are essential for the analysis of flow-driven energy harvesters in order to capture the main physical aspects of the coupled problem and to accurately predict the power output of a harvester. While there are some phenomenological and numerical models for flow-driven harvesters reported in the literature, a fully three-dimensional strongly coupled model has not yet been developed, especially in the context of flow-driven energy harvesting. The weighted residuals method is applied to establish a mixed integral equation describing the incompressible Newtonian flow, elastic substrate structure, piezoelectric patch, equipotential electrode and attached electric circuit that form the multiphysics fluid-structure interaction problem. A monolithic numerical solution method is derived that provides consistent and simultaneous solution to all physical fields as well as to fluid mesh deformation. The approximate solution is based on a mixed space-time finite element discretization with static condensation of the auxiliary fields. The discontinuous Galerkin method is utilized for integrating the monolithic model in time. The proposed solution scheme is illustrated in the example of a lid driven cavity with a flexible piezoelectric bottom wall, demonstrating quantification of the amount of electrical energy extractable from fluid flow by means of a piezoelectric harvester device. The results indicate that in order to make reliable predictions on the power output under varying operational states, the realization of strong multiphysics coupling is required for the mathematical model as well as the numerical solution scheme to capture the characteristics of flow-driven energy harvesters. [less ▲] Detailed reference viewed: 168 (15 UL)Non-localized Contact Between Beams with Non-Circular Cross Sections Magliulo, Marco ; Zilian, Andreas ; Beex, Lars in Proceedings in Applied Mathematics and Mechanics (2019) In this contribution, we introduce a contact formulation between beams finite elements with (hyper)elliptical cross sections. The contact scheme allows to model scenarios in which the contact area is ... [more ▼] In this contribution, we introduce a contact formulation between beams finite elements with (hyper)elliptical cross sections. The contact scheme allows to model scenarios in which the contact area is finite or the contact area occurs along a line. Although some contact schemes are yet able to do this, they require one of the beams to have a circular cross section. Here however, we focus on non-circular cross-sections. Consequently, new projections are required, in which the beam surfaces are used explicitly to formulate contact kinematics. [less ▲] Detailed reference viewed: 80 (12 UL)Model order reduction applied to ALE‐fluid dynamics Baroli, Davide ; Zilian, Andreas in Proceedings in Applied Mathematics and Mechanics (2019) Detailed reference viewed: 32 (1 UL)Phase-field predictive model for setting of fresh self-compacting concrete Huang, Haiqin ; Zilian, Andreas in Proceedings in Applied Mathematics and Mechanics (2018), 18(1), The initial setting of fresh concrete is mainly caused by the dissolution of cement grains and the precipitation of calcium-silicate-hydrates during cement hydration. Progressing hydration drives the ... [more ▼] The initial setting of fresh concrete is mainly caused by the dissolution of cement grains and the precipitation of calcium-silicate-hydrates during cement hydration. Progressing hydration drives the transition from a dense suspension to a porous solid phase. Fresh mixture of self-compacting concrete (SCC) can be considered as a phase-changing multi-component material and can be described as a continuum at the macro scale, interacting with a set of transport-reaction-diffusion processes which in turn are driven by phenomena at the level of the microstructure. This contribution focuses on a predictive model for the setting of fresh SCC where the liquid-solid phase transition is captured by a phase-field variable using the Ginzburg-Landau type free energy function. Hydration-related chemical reactions together with heat and mass transfer are volume coupled with the mechanical behaviour and determined by the environmental conditions. The weak form of the predictive model is discretised using the finite element method and implemented with the FEniCS computational framework. [less ▲] Detailed reference viewed: 83 (8 UL)Determination of process-dependent rheological properties of fresh concrete during setting using a suspension/phase-field model Huang, Haiqin ; Zilian, Andreas Scientific Conference (2018, June) Detailed reference viewed: 75 (19 UL)Analysis of FSI-driven energy harvesting devices using parameterised reduced-order models Zilian, Andreas ; Baroli, Davide Scientific Conference (2018, June) Detailed reference viewed: 71 (4 UL)Modelling and analysis of flow-driven energy harvesting devices and associated reduced order models Zilian, Andreas ; Baroli, Davide Scientific Conference (2018, June) A specific class of energy harvester devices for renewable energy resources allows conversion of ambient fluid flow energy to electrical energy via flow-induced vibrations of a piezo-ceramic composite ... [more ▼] A specific class of energy harvester devices for renewable energy resources allows conversion of ambient fluid flow energy to electrical energy via flow-induced vibrations of a piezo-ceramic composite structure positioned in the flow field. This energy converter technology simultaneously involves the interaction of a composite structure and a surrounding fluid, the electric charge accumulated in the piezo-ceramic material and a controlling electrical circuit. In order to predict the efficiency and operational properties of such future devices and to increase their robustness and performance, a mathematical and numerical model of the complex physical system is required to allow systematic computational investigation of the involved phenomena and coupling characteristics. The presentation will discuss a monolithic modelling approach that allows simultaneous analysis of the harvester, which involves surface-coupled fluid-structure interaction, volume-coupled electro-mechanics and a controlling energy harvesting circuit. Based on a finite element discretisation of the weighted residual form of the governing equations, time- and frequency-domain analysis enables investigation of different types of structures (plate, shells) subject to exterior/interior flow with varying parameters, and attached electrical circuits with respect to the electrical power output generated. Consequently, options for parametric reduced-order modelling of flow-driven energy harvesters will be discussed. [less ▲] Detailed reference viewed: 115 (7 UL)XDMF and ParaView: checkpointing format Habera, Michal ; Zilian, Andreas ; Hale, Jack et al Scientific Conference (2018, March 21) Checkpointing, i.e. saving and reading results of finite element computation is crucial, especially for long-time running simulations where execution is interrupted and user would like to restart the ... [more ▼] Checkpointing, i.e. saving and reading results of finite element computation is crucial, especially for long-time running simulations where execution is interrupted and user would like to restart the process from last saved time step. On the other hand, visualization of results in thid-party software such as ParaView is inevitable. In the previous DOLFIN versions (2017.1.0 and older) these two functionalities were strictly separated. Results could have been saved via HDF5File interface for later computations and/or stored in a format understood by ParaView - VTK’s .pvd (File interface) or XDMF (XDMFFile interface). This led to data redundancy and error-prone workflow. The problem essentially originated from incompatibilities between both libraries, DOLFIN and ParaView (VTK). DOLFIN’s internal representation of finite element function is based on vector of values of degrees of freedom (dofs) and their ordering within cells (dofmap). VTK’s representation of a function is given by it’s values at some points in cell, while ordering and geometric position of these points is fixed and standardised within VTK specification. For nodal (iso- and super-parametric) Lagrange finite elements (Pk , dPk ) both representations coincide up to an ordering. This allows to extend XDMF specification and introduce intermediate way of storing finite element function - intrinsic to both, ParaView and DOLFIN. The necessary work was done as a part of Google Summer of Code 2017 project Develop XDMF for- mat for visualisation and checkpointing, see https://github.com/michalhabera/gsoc-summary. New checkpointing functionality is exposed via write checkpoint() and read checkpoint() methods. [less ▲] Detailed reference viewed: 213 (29 UL)Multiscale Modeling of Discrete Mesomodels for Dry-Woven Fabrics Magliulo, Marco ; Beex, Lars ; Zilian, Andreas Scientific Conference (2018, March) Detailed reference viewed: 36 (1 UL)Phase-field predictive model for setting of fresh self-compacting concrete Huang, Haiqin ; Zilian, Andreas Scientific Conference (2018, March) Detailed reference viewed: 102 (11 UL)2018 Lab Report Computational Mechanics Legato-team and Z-team Bordas, Stéphane ; Zilian, Andreas Report (2018) Detailed reference viewed: 1292 (149 UL)ECCOMAS Newsletter - Computational and Data Sciences in Luxembourg Beex, Lars ; Bordas, Stéphane ; Hale, Jack et al Report (2018) Detailed reference viewed: 866 (111 UL)Isogeometric analysis of thin Reissner-Mindlin plates and shells: Locking phenomena and generalized local B-bar method Hu, Qingyuan ; ; 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: 211 (15 UL)Numerical prediction of the rheological properties of fresh self-compacting concrete Huang, Haiqin ; Zilian, Andreas Scientific Conference (2017, July 14) Self-Compacting Concrete (SCC) is a high-performance construction material that can simplify classical handling on concrete construction by avoiding the need for additional vibrational compaction ... [more ▼] Self-Compacting Concrete (SCC) is a high-performance construction material that can simplify classical handling on concrete construction by avoiding the need for additional vibrational compaction. Challenges in the use of SCC lie in ensuring optimal operation of the material in terms of properly filled castings in presence of complex reinforcement arrangements, reduction of entrained gas bubbles and limitation of aggregate separation. A major factor influencing the aforementioned aspects is the rheological properties of SCC mixtures under varying conditions (e.g. content composition, mechanical impact, temperature, moisture). This contribution aims at unified constitutive modelling of SCC in the setting stage. Concrete setting describes the transition from fluid-like fresh concrete, which -in presence of time- dependent transport-reaction processes- develops a porous cementitious structure, to hardened concrete showing solid-like behaviour. The constitutive model is implemented using the open-source finite element framework FENICS and applied to a number of benchmark problems. [less ▲] Detailed reference viewed: 131 (24 UL)Time and frequency domain analysis of piezoelectric energy harvesters by monolithic finite element modeling Ravi, Srivathsan ; Zilian, Andreas in International Journal for Numerical Methods in Engineering (2017) The successful design of piezoelectric energy harvesting devices relies upon the identification of optimal geometrical and material configurations to maximize the power output for a specific band of ... [more ▼] The successful design of piezoelectric energy harvesting devices relies upon the identification of optimal geometrical and material configurations to maximize the power output for a specific band of excitation frequencies. Extendable predictive models and associated approximate solution methods are essential for analysis of a wide variety of future advanced energy harvesting devices involving more complex geometries and material distributions. Based on a holistic continuum mechanics modeling approach to the multi-physics energy harvesting problem, this article proposes a monolithic numerical solution scheme using a mixed-hybrid 3-dimensional finite element formulation of the coupled governing equations for analysis in time and frequency domain. The weak form of the electromechanical/circuit system uses velocities and potential rate within the piezoelectric structure, free boundary charge on the electrodes, and potential at the level of the generic electric circuit as global degrees of freedom. The approximation of stress and dielectric displacement follows the work by Pian, Sze, and Pan. Results obtained with the proposed model are compared with analytical results for the reduced-order model of a cantilevered bimorph harvester with tip mass reported in the literature. The flexibility of the method is demonstrated by studying the influence of partial electrode coverage on the generated power output. [less ▲] Detailed reference viewed: 139 (17 UL)Monolithic modeling and finite element analysis of piezoelectric energy harvesters Ravi, Srivathsan ; Zilian, Andreas in Acta Mechanica (2017), 228(6), 2251-2267 This paper is devoted to monolithic modeling of piezoelectric energy harvesting devices. From a modeling perspective, piezoelectric energy harvesting is a strongly coupled phenomenon with two-way coupling ... [more ▼] This paper is devoted to monolithic modeling of piezoelectric energy harvesting devices. From a modeling perspective, piezoelectric energy harvesting is a strongly coupled phenomenon with two-way coupling between the electromechanical effect of the piezoelectric material and the harvesting circuit. Even in applications related to shunt damping, where the attached electrical circuit is passive, accurate modeling of the strong coupling is crucial for proper evaluation of the relevant parameters. The article proposes a monolithic mixed-hybrid finite element formulation for the predictive modeling and simulation of piezoelectric energy harvesting devices. The governing equations of the coupled electromechanical problem are converted into a single integral form with six independent unknown fields. Such a holistic approach provides consistent solution to the coupled field equations which involve structural dynamics, electromechanical effect of the piezoelectric patches and the dynamics of the attached harvesting circuit. This allows accurate computation of the eigenvalues and corresponding mode shapes of a harvester for any finite resistive load coupled to the harvester. The fully three-dimensional mixed-hybrid formulation is capable of analyzing structures with non-uniform geometry and varying material properties. The results of the finite element model are verified against the analytical results of a bimorph harvester with tip mass reported in the literature. [less ▲] Detailed reference viewed: 126 (22 UL) |
||