Search

While solid mechanics codes are now conventional tools both in industry and research, the increasingly more exigent requirements of both sectors are fuelling the need for more computational power and more advanced algorithms. For obvious reasons, commercial codes are lagging behind academic codes often dedicated either to the implementation of one new technique, or the upscaling of current conventional codes to tackle massively large scale computational problems. Only in a few cases, both approaches have been followed simultaneously. In this article, a solid mechanics simulation strategy for parallel supercomputers based on a hybrid approach is presented. Hybrid parallelization exploits the thread-level parallelism of multicore architectures, combining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multi-physics code. Hybrid parallelization is specially well suited for the current trend of supercomputers, namely large clusters of multicores. The strategy is assessed through transient non-linear solid mechanics problems, both for explicit and implicit schemes, running on thousands of cores. In order to demonstrate the flexibility of the proposed strategy under advance algorithmic evolution of computational mechanics, a non-local parallel overset meshes method (Chimera-like) is implemented and the conservation of the scalability is demonstrated., QC 20151104

Because of their widely generalized use in many industries, composites are the subject of many research campaigns. More particularly, the development of both accurate and flexible numerical models able to capture their intrinsically multiscale modes of failure is still a challenge. The standard finite element method typically requires intensive remeshing to adequately capture the geometry of the cracks and high accuracy is thus often sacrificed in favor of scalability, and vice versa. In an effort to preserve both properties, we present here an extended finite element method (XFEM) for large scale composite fracture simulations. In this formulation, the standard FEM formulation is partially enriched by use of shifted Heaviside functions with special attention paid to the scalability of the scheme. This enrichment technique offers several benefits since the interpolation property of the standard shape function still holds at the nodes. Those benefits include (i) no extra boundary condition for the enrichment degree of freedom, and (ii) no need for transition/blending regions; both of which contribute to maintaining the scalability of the code.Two different cohesive zone models (CZM) are then adopted to capture the physics of the crack propagation mechanisms. At the intralaminar level, an extrinsic CZM embedded in the XFEM formulation is used. At the interlaminar level, an intrinsic CZM is adopted for predicting the failure. The overall framework is implemented in ALYA, a mechanics code specifically developed for large scale, massively parallel simulations of coupled multi-physics problems. The implementation of both intrinsic and extrinsic CZM models within the code is such that it conserves the extremely efficient scalability of ALYA while providing accurate physical simulations of computationally expensive phenomena. The strong scalability provided by the proposed implementation is demonstrated. The model is ultimately validated against a full experimental campaign, QC 20150602

Multiscale modelling and simulation play an important role in sheet metal forming analysis, since the overall material responses at macroscopic engineering scales, e.g. formability and anisotropy, are strongly influenced by microstructural properties, such as grain size and crystal orientations (texture). In the present report, multiscale analysis on deep drawing of dual-phase steels is performed using an efficient grain cluster-based homogenization scheme. The homogenization scheme, called relaxed grain cluster (RGC), is based on a generalization of the grain cluster concept, where a (representative) volume element consists of p x q x r (hexahedral) grains. In this scheme, variation of the strain or deformation of individual grains is taken into account through the, so-called, interface relaxation, which is formulated within an energy minimization framework. An interfacial penalty term is introduced into the energy minimization framework in order to account for the effects of grain boundaries. The grain cluster-based homogenization scheme has been implemented and incorporated into the advanced material simulation platform DAMASK, which purposes to bridge the macroscale boundary value problems associated with deep drawing analysis to the micromechanical constitutive law, e.g. crystal plasticity model. Standard Lankford anisotropy tests are performed to validate the model parameters prior to the deep drawing analysis. Model predictions for the deep drawing simulations are analyzed and compared to the corresponding experimental data. The result shows that the predictions of the model are in a very good agreement with the experimental measurement., QC 20150612

A continuum material model is developed for simulating the mechanical response of high-density cellulose-based materials subjected to stationary and transient loading. The model is formulated in an infinitesimal strain framework, where the total strain is decomposed into elastic and plastic parts. The model adopts a standard linear viscoelastic solid model expressed in terms of Boltzmann hereditary integral form, which is coupled to a rate-dependent viscoplastic formulation to describe the irreversible plastic part of the overall strain. An anisotropic hardening law with a kinematic effect is particularly adopted in order to capture the complex stress-strain hysteresis typically observed in polymeric materials. In addition, the present model accounts for the effects of material densification associated with through-thickness compression, which are captured using an exponential law typically applied in the continuum description of elasticity in porous media. Material parameters used in the present model are calibrated to the experimental data for high-density (press)boards. The experimental characterization procedures as well as the calibration of the parameters are highlighted. The results of the model simulations are systematically analyzed and validated against the corresponding experimental data. The comparisons show that the predictions of the present model are in very good agreement with the experimental observations for both stationary and transient load cases., QC 20151027

The tensile response of dense fiber-based materials like paper or paperboard is mainly dependent of the properties of the fibers, which store most of the elastic energy. In this paper, we investigate the influence of geometrical and material parameters on the mechanical response of the pulp fibers used in paper manufacturing. We developed a three-dimensional finite element model of the fiber, which accounts for microfibril orientation of cellulose fibril, and the presence of lignin in the secondary cell wall. The results showed that the change in the microfibril orientation upon axial straining is mainly a geometrical effect, and is independent of the material properties of the fiber, as long as the deformations are elastic. Plastic strain accelerates the change in microfibril orientation and thus makes it material-dependent. The results also showed that the elastic modulus of the fiber has a non-linear dependency on microfibril angle, with elastic modulus being more sensitive to the change of microfibril angle around small initial values of microfibril angles. Based on numerical results acquired from a 3D fiber model supported by available experimental evidence, we propose an anisotropic-kinematic hardening plasticity model for a fiber within a beam framework. The proposed fiber model is capable of reproducing the main features of the cyclic tensile response of a pulp fiber, such as stiffening due to changing microfibril angle. The constitutive model of the fiber was implemented in a finite-element model of the fiber network. By using the fiber network model, we estimated the level of strain that fiber segments accumulate before the typical failure strain of the entire network is reached., QC 20151113

Because of their widely generalized use in many industries, composites are the subject of many research campaigns. More particularly, the development of both accurate and flexible numerical models able to capture their intrinsically multiscale modes of failure is still a challenge. The standard finite element method typically requires intensive remeshing to adequately capture the geometry of the cracks and high accuracy is thus often sacrificed in favor of scalability, and vice versa. In an effort to preserve both properties, we present here an extended finite element method (XFEM) for large scale composite fracture simulations. In this formulation, the standard FEM formulation is partially enriched by use of shifted Heaviside functions with special attention paid to the scalability of the scheme. This enrichment technique offers several benefits since the interpolation property of the standard shape function still holds at the nodes. Those benefits include (i) no extra boundary condition for the enrichment degree of freedom, and (ii) no need for transition/blending regions; both of which contribute to maintaining the scalability of the code. Two different cohesive zone models (CZM) are then adopted to capture the physics of the crack propagation mechanisms. At the intralaminar level, an extrinsic CZM embedded in the XFEM formulation is used. At the interlaminar level, an intrinsic CZM is adopted for predicting the failure. The overall framework is implemented in ALYA, a mechanics code specifically developed for large scale, massively parallel simulations of coupled multi-physics problems. The implementation of both intrinsic and extrinsic CZM models within the code is such that it conserves the extremely efficient scalability of ALYA while providing accurate physical simulations of computationally expensive phenomena. The strong scalability provided by the proposed implementation is demonstrated. The model is ultimately validated against a full experimental campaig, A.J., A.M., D.T., L.N. and L.W. acknowledge funding through the SIMUCOMP ERA-NET MATERA + project financed by the Fonds National de la Recherche (FNR) of Luxembourg, the Consejería de Educación y Empleo of the Comunidad de Madrid, the Walloon region (agreement no 1017232, CT-EUC 2010–10-12), and by the European Unions Seventh Framework Programme (FP7/2007–2013)., Peer Reviewed, Postprint (author's final draft)

© 2015 Elsevier Ltd. Because of their widely generalized use in many industries, composites are the subject of many research campaigns. More particularly, the development of both accurate and flexible numerical models able to capture their intrinsically multiscale modes of failure is still a challenge. The standard finite element method typically requires intensive remeshing to adequately capture the geometry of the cracks and high accuracy is thus often sacrificed in favor of scalability, and vice versa. In an effort to preserve both properties, we present here an extended finite element method (XFEM) for large scale composite fracture simulations. In this formulation, the standard FEM formulation is partially enriched by use of shifted Heaviside functions with special attention paid to the scalability of the scheme. This enrichment technique offers several benefits since the interpolation property of the standard shape function still holds at the nodes. Those benefits include (i) no extra boundary condition for the enrichment degree of freedom, and (ii) no need for transition/blending regions; both of which contribute to maintaining the scalability of the code.Two different cohesive zone models (CZM) are then adopted to capture the physics of the crack propagation mechanisms. At the intralaminar level, an extrinsic CZM embedded in the XFEM formulation is used. At the interlaminar level, an intrinsic CZM is adopted for predicting the failure. The overall framework is implemented in ALYA, a mechanics code specifically developed for large scale, massively parallel simulations of coupled multi-physics problems. The implementation of both intrinsic and extrinsic CZM models within the code is such that it conserves the extremely efficient scalability of ALYA while providing accurate physical simulations of computationally expensive phenomena. The strong scalability provided by the proposed implementation is demonstrated. The model is ultimately validated against a full

While solid mechanics codes are now conventional tools both in industry and research, the increasingly more exigent requirements of both sectors are fuelling the need for more computational power and more advanced algorithms. For obvious reasons, commercial codes are lagging behind academic codes often dedicated either to the implementation of one new technique, or the upscaling of current conventional codes to tackle massively large scale computational problems. Only in a few cases, both approaches have been followed simultaneously. In this article, a solid mechanics simulation strategy for parallel supercomputers based on a hybrid approach is presented. Hybrid parallelization exploits the thread-level parallelism of multicore architectures, combining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multi-physics code. Hybrid parallelization is specially well suited for the current trend of supercomputers, namely large clusters of multicores. The strategy is assessed through transient non-linear solid mechanics problems, both for explicit and implicit schemes, running on thousands of cores. In order to demonstrate the flexibility of the proposed strategy under advance algorithmic evolution of computational mechanics, a non-local parallel overset meshes method (Chimera-like) is implemented and the conservation of the scalability is demonstrated.

High density (HD) pressboard is an essential element in power transformers combining good electrical insulation properties with effective mechanical characteristics that well suit design requirements of power transformers. In order to ensure a correctly functioning transformer, it is therefore very important to characterize and to understand the mechanical properties of pressboard under different operating conditions. Pressboard is composed of natural polymeric chains, whose mechanical properties are affected by moisture and temperature. Moreover, temperature and moisture conditions in power transformers vary throughout manufacturing process and service/operation life-time. An accurate definition of the mechanical properties is, therefore, necessary. The present article focuses on the effect of different combinations of temperature/moisture and mechanical load on the deformation behavior of HD pressboard samples. A mechanical constitutive model is developed for finite element (FEM) simulation based on a viscoelastic–viscoplastic material description. Special attention is given on the complex through-thickness deformation behavior of HD pressboard. Thorough analyses are performed based on the comparisons between the results of experimental characterization and FEM modeling and simulations. The good agreement between experimental and modeling results shows a great potential for application in mechanical design of transformer insulation., QC 20150413

While solid mechanics codes are now conventional tools both in industry and research, the increasingly more exigent requirements of both sectors are fuelling the need for more computational power and more advanced algorithms. For obvious reasons, commercial codes are lagging behind academic codes often dedicated either to the implementation of one new technique, or the upscaling of current conventional codes to tackle massively large scale computational problems. Only in a few cases, both approaches have been followed simultaneously. In this article, a solid mechanics simulation strategy for parallel supercomputers based on a hybrid approach is presented. Hybrid parallelization exploits the thread-level parallelism of multicore architectures, combining MPI tasks with OpenMP threads. This paper describes the proposed strategy, programmed in Alya, a parallel multi-physics code. Hybrid parallelization is specially well suited for the current trend of supercomputers, namely large clusters of multicores. The strategy is assessed through transient non-linear solid mechanics problems, both for explicit and implicit schemes, running on thousands of cores. In order to demonstrate the flexibility of the proposed strategy under advance algorithmic evolution of computational mechanics, a non-local parallel overset meshes method (Chimera-like) is implemented and the conservation of the scalability is demonstrated., D.D.T and A.J acknowledge funding through SIMUCOMP and ERA-NET MATERA+ project financed by the Consejería de Educación y Empleo of the Comunidad de Madrid and by the European Union’s Seventh Framework Programme (FP7/2007-2013). Thisworkwas partially supported by the grant SEV-2011-00067, Severo Ochoa Program, awarded by the Spanish Government. The authors’ would like to acknowledge PRACE infrastructure support., Peer Reviewed, Postprint (published version)

Cellulose-based components constitute the bulk of the current insulation for transformers. Cellulose is an organic polymer material which combines excellent electrical properties and good mechanical performance. As a polymeric material, cellulose is very sensitive to moisture and temperature. These factors can influence the electrical and mechanical performance of a transformer throughout its lifetime. In order to ensure the quality of the product during transformer manufacturing, as well as during transformer life-time services, adequate models for predicting the physical properties of its constituents are needed. The present investigation tackles the mechanical description of pressboard. For this purpose, a three dimensional mechanical model is developed for simulating the in-plane and out-of-plane behavior of the pressboard material. The model is based on an anisotropic viscoelastic-viscoplastic constitutive law, which includes features that are particular for cellulose-based materials, e.g. the peculiar double nature of fiber-network-based and porous material. The material is orthotropic by nature, i.e. the in-plane mechanical properties markedly differ from the out-of-plane ones. Particular regard is taken when considering the effect of out-of-plane stresses which both cause viscous deformation and permanent compaction. The analyses on the mechanical behavior of pressboard are performed by comparing the experimental data on pressboard and the results of model simulations., QC 20131218

The recently developed hybrid discontinuous Galerkin/extrinsic cohesive law framework is extended to the study of intra-laminar fracture of composite materials. Toward this end, micro-volumes of different sizes are studied. The method captures the debonding process, which is herein proposed to be assimilated to a damaging process, before the strain softening onset, and the density of dissipated energy resulting from the damage (debonding) remains the same for the different studied cell sizes. Finally, during the strain softening phase a micro-crack initiates and propagates in agreement with experimental observations. We thus extract a resulting mesoscale cohesive law, which is independent on the cell sizes, using literature methods., QC 20131125