BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jan Zeman

Jan Zeman appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.mtrl-scimath.NAComputational Engineering, Finance, and Sciencemath.APNumerical Analysisphysics.comp-phComputational Geometrycond-mat.softcond-mat.stat-mechmath.OCphysics.flu-dyn

Catalog footprint

What is connected

33works
11topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

Modular-topology optimization of structures and mechanisms with free material design and clustering

Topology optimization of modular structures and mechanisms enables balancing the performance of automatically-generated individualized designs, as required by Industry 4.0, with enhanced sustainability by means of component reuse. For optimal modular design, two key questions must be answered: (i) what should the topology of individual modules be like and (ii) how should modules be arranged at the product scale? We address these challenges by proposing a bi-level sequential strategy that combines free material design, clustering techniques, and topology optimization. First, using free material optimization enhanced with post-processing for checkerboard suppression, we determine the distribution of elasticity tensors at the product scale. To extract the sought-after modular arrangement, we partition the obtained elasticity tensors with a novel deterministic clustering algorithm and interpret its outputs within Wang tiling formalism. Finally, we design interiors of individual modules by solving a single-scale topology optimization problem with the design space reduced by modular mapping, conveniently starting from an initial guess provided by free material optimization. We illustrate these developments with three benchmarks first, covering compliance minimization of modular structures, and, for the first time, the design of non-periodic compliant modular mechanisms. Furthermore, we design a set of modules reusable in an inverter and in gripper mechanisms, which ultimately pave the way towards the rational design of modular architectured (meta)materials.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Optimal FFT-accelerated Finite Element Solver for Homogenization

We propose a matrix-free finite element (FE) homogenization scheme that is considerably more efficient than generic FE implementations. The efficiency of our scheme follows from a preconditioned well-scaled reformulation allowing for the use of the conjugate gradient or similar iterative solvers. The geometrically-optimal preconditioner -- a discretized Green's function of a periodic homogeneous reference problem -- has a block-diagonal structure in the Fourier space which permits its efficient inversion using the fast Fourier transform (FFT) techniques for generic regular meshes. This implies that the scheme scales as $\mathcal{O}(n \log(n))$ like FFT, rendering it equivalent to spectral solvers in terms of computational efficiency. However, in contrast to classical spectral solvers, the proposed scheme works with FE shape functions with local supports and is free of the Fourier ringing phenomenon. We showcase that the scheme achieves the number of iterations that are almost independent of spatial discretisation and scales mildly with the phase contrast. Additionally, we discuss the equivalence between our displacement-based scheme and the recently proposed strain-based homogenization technique with finite-element projection.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Elimination of ringing artifacts by finite-element projection in FFT-based homogenization

Micromechanical homogenization is often carried out with Fourier-accelerated methods that are prone to ringing artifacts. We here generalize the compatibility projection introduced by Vondřejc, Zeman & Marek [Comput. Math. Appl. 68, 156 (2014)] beyond the Fourier basis. In particular, we formulate the compatibility projection for linear finite elements while maintaining Fourier-acceleration and the fast convergence properties of the original method. We demonstrate that this eliminates ringing artifacts and yields an efficient computational homogenization scheme that is equivalent to canonical finite-element formulations on fully structured grids.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global optimality in minimum compliance topology optimization of frames and shells by moment-sum-of-squares hierarchy

The design of minimum-compliance bending-resistant structures with continuous cross-section parameters is a challenging task because of its inherent non-convexity. Our contribution develops a strategy that facilitates computing all guaranteed globally optimal solutions for frame and shell structures under multiple load cases and self-weight. To this purpose, we exploit the fact that the stiffness matrix is usually a polynomial function of design variables, allowing us to build an equivalent non-linear semidefinite programming formulation over a semi-algebraic feasible set. This formulation is subsequently solved using the Lasserre moment-sum-of-squares hierarchy, generating a sequence of outer convex approximations that monotonically converges from below to the optimum of the original problem. Globally optimal solutions can subsequently be extracted using the Curto-Fialkow flat extension theorem. Furthermore, we show that a simple correction to the solutions of the relaxed problems establishes a feasible upper bound, thereby deriving a simple sufficient condition of global $\varepsilon$-optimality. When the original problem possesses a unique minimum, we show that this solution is found with a zero optimality gap in the limit. These theoretical findings are illustrated on several examples of topology optimization of frames and shells, for which we observe that the hierarchy converges in a finite (rather small) number of steps.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Guaranteed two-sided bounds on all eigenvalues of preconditioned diffusion and elasticity problems solved by the finite element method

A method of estimating all eigenvalues of a preconditioned discretized scalar diffusion operator with Dirichlet boundary conditions has been recently introduced in T. Gergelits, K.A. Mardal, B.F. Nielsen, Z. Strakoš: Laplacian preconditioning of elliptic PDEs: Localization of the eigenvalues of the discretized operator, SIAM Journal on Numerical Analysis 57(3) (2019), 1369-1394. Motivated by this paper, we offer a slightly different approach that extends the previous results in some directions. Namely, we provide bounds on all increasingly ordered eigenvalues of a general diffusion or elasticity operator with tensor data, discretized with the conforming finite element method, preconditioned by the inverse of a matrix of the same operator with different data. Our results hold for mixed Dirichlet and Robin or periodic boundary conditions applied to the original and preconditioning problems. The bounds are two-sided, guaranteed, easily accessible, and depend solely on the material data.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Level-set based design of Wang tiles for modelling complex microstructures

Microstructural geometry plays a critical role in the response of heterogeneous materials. Consequently, methods for generating microstructural samples are increasingly crucial to advanced numerical analyses. We extend Sonon et al.'s unified framework, developed originally for generating particulate and foam-like microstructural geometries of Periodic Unit Cells, to non-periodic microstructural representations based on the formalism of Wang tiles. This formalism has been recently proposed in order to generalize the Periodic Unit Cell approach, enabling a fast synthesis of arbitrarily large, stochastic microstructural samples from a handful of domains with predefined compatibility constraints. However, a robust procedure capable of designing complex, three-dimensional, foam-like and cellular morphologies of Wang tiles has not yet been proposed. This contribution fills the gap by significantly broadening the applicability of the tiling concept. Since the original Sonon et al.'s framework builds on a random sequential addition of particles enhanced with an implicit representation of particle boundaries by the level-set field, we first devise an analysis based on a connectivity graph of a tile set, resolving the question where a particle should be copied when it intersects a tile boundary. Next, we introduce several modifications to the original algorithm that are necessary to ensure microstructural compatibility in the generalized periodicity setting of Wang tiles. Having established a universal procedure for generating tile morphologies, we compare strictly aperiodic and stochastic sets with the same cardinality in terms of reducing the artificial periodicity in reconstructed microstructural samples. We demonstrate the superiority of the vertex-defined tile sets for two-dimensional problems and illustrate the capabilities of the algorithm with two- and three-dimensional examples.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Reduced integration schemes in micromorphic computational homogenization of elastomeric mechanical metamaterials

Exotic behaviour of mechanical metamaterials often relies on an internal transformation of the underlying microstructure triggered by its local instabilities, rearrangements, and rotations. Depending on the presence and magnitude of such a transformation, effective properties of a metamaterial may change significantly. To capture this phenomenon accurately and efficiently, homogenization schemes are required that reflect microstructural as well as macro-structural instabilities, large deformations, and non-local effects. To this end, a micromorphic computational homogenization scheme has recently been developed, which employs the particular microstructural transformation as a non-local mechanism, magnitude of which is governed by an additional coupled partial differential equation. Upon discretizing the resulting problem it turns out that the macroscopic stiffness matrix requires integration of macro-element basis functions as well as their derivatives, thus calling for a higher-order integration rules. Because evaluation of constitutive law in multiscale schemes involves an expensive solution of a non-linear boundary value problem, computational efficiency can be improved by reducing the number of integration points. Therefore, the goal of this paper is to investigate reduced-order schemes in computational homogenization, with emphasis on the stability of the resulting elements. In particular, arguments for lowering the order of integration from the expensive mass-matrix to a cheaper stiffness-matrix equivalent are first outlined. An efficient one-point integration quadrilateral element is then introduced and proper hourglass stabilization discussed. Performance of the resulting set of elements is finally tested on a benchmark bending example, showing that we achieve accuracy comparable to the full quadrature rules.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comparative study on low-memory iterative solvers for FFT-based homogenization of periodic media

In this paper, we assess the performance of four iterative algorithms for solving non-symmetric rank-deficient linear systems arising in the FFT-based homogenization of heterogeneous materials defined by digital images. Our framework is based on the Fourier-Galerkin method with exact and approximate integrations that has recently been shown to generalize the Lippmann-Schwinger setting of the original work by Moulinec and Suquet from 1994. It follows from this variational format that the ensuing system of linear equations can be solved by general-purpose iterative algorithms for symmetric positive-definite systems, such as the Richardson, the Conjugate gradient, and the Chebyshev algorithms, that are compared here to the Eyre-Milton scheme - the most efficient specialized method currently available. Our numerical experiments, carried out for two-dimensional elliptic problems, reveal that the Conjugate gradient algorithm is the most efficient option, while the Eyre-Milton method performs comparably to the Chebyshev semi-iteration. The Richardson algorithm, equivalent to the still widely used original Moulinec-Suquet solver, exhibits the slowest convergence. Besides this, we hope that our study highlights the potential of the well-established techniques of numerical linear algebra to further increase the efficiency of FFT-based homogenization methods.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A Variational Formulation of Dissipative Quasicontinuum Methods

Lattice systems and discrete networks with dissipative interactions are successfully employed as meso-scale models of heterogeneous solids. As the application scale generally is much larger than that of the discrete links, physically relevant simulations are computationally expensive. The QuasiContinuum (QC) method is a multiscale approach that reduces the computational cost of direct numerical simulations by fully resolving complex phenomena only in regions of interest while coarsening elsewhere. In previous work (Beex et al., J. Mech. Phys. Solids 64, 154-169, 2014), the originally conservative QC methodology was generalized to a virtual-power-based QC approach that includes local dissipative mechanisms. In this contribution, the virtual-power-based QC method is reformulated from a variational point of view, by employing the energy-based variational framework for rate-independent processes (Mielke and Roubíček, Rate-Independent Systems: Theory and Application, Springer-Verlag, 2015). By construction it is shown that the QC method with dissipative interactions can be expressed as a minimization problem of a properly built energy potential, providing solutions equivalent to those of the virtual-power-based QC formulation. The theoretical considerations are demonstrated on three simple examples. For them we verify energy consistency, quantify relative errors in energies, and discuss errors in internal variables obtained for different meshes and two summation rules.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Integrated Experimental-Numerical Study of the Performance of Lime-Based Mortars in Masonry Piers Under Eccentric Loading

Architectural conservation and repair are becoming increasingly important issues in many countries due to numerous prior improper interventions, including the use of inappropriate repair materials over time. As a result, the composition of repair masonry mortars is now being more frequently addressed in mortar research. Just recently, for example, it has become apparent that Portland cement mortars, extensively exploited as repair mortars over the past few decades, are not suitable for repair because of their chemical, physical, mechanical, and aesthetic incompatibilities with original materials. This paper focuses on the performance of various lime-based alternative materials intended for application in repairing historic structures when subjected to mechanical loading. Results of basic material tests indicate that the use of metakaolin as a pozzolanic additive produces mortars with superior strength and sufficiently low shrinkage. Moreover, mortar strength can be further enhanced by the addition of crushed brick fragments, which explains the longevity of Roman concretes rich in pozzolans and aggregates from crushed clay products such as tiles, pottery, or bricks. An integrated experimental-numerical approach was used to identify key mortar parameters influencing the load-bearing capacity of masonry piers subjected to a combination of compression and bending. The simulations indicate increased load-bearing capacities for masonry piers containing metakaolin-rich mortars with crushed brick fragments, as a result of their superior compressive strength.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational Homogenization of Fresh Concrete Flow Around Reinforcing Bars

Motivated by casting of fresh concrete in reinforced concrete structures, we introduce a numerical model of a steady-state non-Newtonian fluid flow through a porous domain. Our approach combines homogenization techniques to represent the reinforced domain by the Darcy law with an interface coupling of the Stokes and Darcy flows through the Beavers-Joseph-Saffman conditions. The ensuing two-scale problem is solved by the Finite Element Method with consistent linearization and the results obtained from the homogenization approach are verified against fully resolved direct numerical simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Guaranteed upper-lower bounds on homogenized properties by FFT-based Galerkin method

Guaranteed upper-lower bounds on homogenized coefficients, arising from the periodic cell problem, are calculated in a scalar elliptic setting. Our approach builds on the recent variational reformulation of the Moulinec-Suquet (1994) Fast Fourier Transform (FFT) homogenization scheme by Vondřejc et al. (2014), which is based on the conforming Galerkin approximation with trigonometric polynomials. Upper-lower bounds are obtained by adjusting the primal-dual finite element framework developed independently by Dvořák (1993) and Wieckowski (1995) to the FFT-based Galerkin setting. We show that the discretization procedure differs for odd and non-odd number of grid points. Thanks to the Helmholtz decomposition inherited from the continuous formulation, the duality structure is fully preserved for the odd discretizations. In the latter case, a more complex primal-dual structure is observed due to presence of the trigonometric polynomials associated with the Nyquist frequencies. These theoretical findings are confirmed with numerical examples. To conclude, the main advantage of the FFT-based approach over conventional finite-element schemes is that the primal and the dual problems are treated on the same basis, and this property can be extended beyond the scalar elliptic setting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Localization Study of a Regularized Variational Damage Model

The paper presents a detailed analysis and extended formulation of a rate-independent regularized damage model proposed by Mielke and Roubíček (2006). Localization properties are studied in the context of a simple one-dimensional problem, but the results reveal the fundamental features of the basic model and of its modified versions. The initial bifurcation from a uniform solution is described analytically while the complete failure process is studied numerically. Modifications of the regularizing term and of the dissipation distance are introduced and their effect on the global response is investigated. It is shown that, with a proper combination of model parameters, a realistic shape of the load-displacement diagram can be achieved and pathological effects such as extremely brittle response or expansion of the damage zone accompanied by stress locking can be eliminated.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Micromechanics-Based Simulations of Compressive and Tensile Testing on Lime-Based Mortars

The purpose of this paper is to propose a continuum micromechanics model for the simulation of uniaxial compressive and tensile tests on lime-based mortars, in order to predict their stiffness, compressive and tensile strengths, and tensile fracture energy. In tension, we adopt an incremental strain-controlled form of the Mori-Tanaka scheme with a damageable matrix phase, while a simple $J_2$ yield criterion is employed in compression. To reproduce the behavior of lime-based mortars correctly, the scheme must take into account shrinkage cracking among aggregates. This phenomenon is introduced into the model via penny-shaped cracks, whose density is estimated on the basis of a particle size distribution combined with the results of finite element analyses of a single crack formation between two spherical inclusions. Our predictions show a good agreement with experimental data and explain the advantages of compliant crushed brick fragments, often encountered in ancient mortars, over stiff sand particles. The validated model provides a reliable tool for optimizing the composition of modern lime-based mortars with applications in conservation and restoration of architectural heritage.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An FFT-based Galerkin Method for Homogenization of Periodic Media

In 1994, Moulinec and Suquet introduced an efficient technique for the numerical resolution of the cell problem arising in homogenization of periodic media. The scheme is based on a fixed-point iterative solution to an integral equation of the Lippmann-Schwinger type, with action of its kernel efficiently evaluated by the Fast Fourier Transform techniques. The aim of this work is to demonstrate that the Moulinec-Suquet setting is actually equivalent to a Galerkin discretization of the cell problem, based on approximation spaces spanned by trigonometric polynomials and a suitable numerical integration scheme. For the latter framework and scalar elliptic setting, we prove convergence of the approximate solution to the weak solution, including a-priori estimates for the rate of convergence for sufficiently regular data and the effects of numerical integration. Moreover, we also show that the variational structure implies that the resulting non-symmetric system of linear equations can be solved by the conjugate gradient method. Apart from providing a theoretical support to Fast Fourier Transform-based methods for numerical homogenization, these findings significantly improve on the performance of the original solver and pave the way to similar developments for its many generalizations proposed in the literature.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Aperiodic compression and reconstruction of real world material systems based on Wang tiles

The paper presents a concept/technique to compress and synthesize complex material morphologies that is based on Wang tilings. Specifically, a microstructure is stored in a set of Wang tiles and its reconstruction is performed by means of a stochastic tiling algorithm. A substantial part of the study is devoted to the setup of optimal parameters of the automatic tile design by means of parametric studies with statistical descriptors at heart. The performance of the method is demonstrated on four two-dimensional two-phase target systems, monodisperse media with hard and soft discs, sandstone, and high porosity metallic foam.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Finite element model based on refined plate theories for laminated glass units

Laminated glass units exhibit complex response as a result of different mechanical behavior and properties of glass and polymer foil. We aim to develop a finite element model for elastic laminated glass plates based on the refined plate theory by Mau. For a geometrically nonlinear description of the behavior of units, each layer behaves according to the Reissner-Mindlin kinematics, complemented with membrane effects and the von Kármán assumptions. Nodal Lagrange multipliers enforce the compatibility of independent layers in this approach. We have derived the discretized model by the energy-minimization arguments, assuming that the unknown fields are approximated by bi-linear functions at the element level, and solved the resulting system by the Newton method with consistent linearization. We have demonstrated through verification and validation examples that the proposed formulation is reliable and accurately reproduces the behavior of laminated glass units. This study represents a first step to the development of a comprehensive, mechanics-based model for laminated glass systems that is suitable for implementation in common engineering finite element solvers.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A framework for integrated design of algorithmic architectural forms

This paper presents a methodology and software tools for parametric design of complex architectural objects, called digital or algorithmic forms. In order to provide a flexible tool, the proposed design philosophy involves two open source utilities Donkey and MIDAS written in Grasshopper algorithm editor and C++, respectively, that are to be linked with a scripting-based architectural modellers Rhinoceros, IntelliCAD and the open source Finite Element solver OOFEM. The emphasis is put on the mechanical response in order to provide architects with a consistent learning framework and an insight into structural behaviour of designed objects. As demonstrated on three case studies, the proposed modular solution is capable of handling objects of considerable structural complexity, thereby accelerating the process of finding procedural design parameters from orders of weeks to days or hours.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Existence results for incompressible magnetoelasticity

We investigate a variational theory for magnetoelastic solids under the incompressibility constraint. The state of the system is described by deformation and magnetization. While the former is classically related to the reference configuration, magnetization is defined in the deformed configuration instead. We discuss the existence of energy minimizers without relying on higher-order deformation gradient terms. Then, by introducing a suitable positively $1$-homogeneous dissipation, a quasistatic evolution model is proposed and analyzed within the frame of energetic solvability.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical model of elastic laminated glass beams under finite strain

Laminated glass structures are formed by stiff layers of glass connected with a compliant plastic interlayer. Due to their slenderness and heterogeneity, they exhibit a complex mechanical response that is difficult to capture by single-layer models even in the elastic range. The purpose of this paper is to introduce an efficient and reliable finite element approach to the simulation of the immediate response of laminated glass beams. It proceeds from a refined plate theory due to Mau (1973), as we treat each layer independently and enforce the compatibility by the Lagrange multipliers. At the layer level, we adopt the finite-strain shear deformable formulation of Reissner (1972) and the numerical framework by Ibrahimbegović and Frey (1993). The resulting system is solved by the Newton method with consistent linearization. By comparing the model predictions against available experimental data, analytical methods and two-dimensional finite element simulations, we demonstrate that the proposed formulation is reliable and provides accuracy comparable to the detailed two-dimensional finite element analyzes. As such, it offers a convenient basis to incorporate more refined constitutive description of the interlayer.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A micromechanics-based model for stiffness and strength estimation of cocciopesto mortars

The purpose of this paper is to propose an inexpensive micromechanics-based scheme for stiffness homogenization and strength estimation of mortars containing crushed bricks, known as cocciopesto. The model utilizes the Mori-Tanaka method to determine the effective stiffness, combined with estimates of quadratic invariants of the deviatoric stresses inside phases to predict the compressive strength. Special attention is paid to the representation of C-S-H gel layer around bricks and interfacial transition zone around sand aggregates, which renders the predictions sensitive to particle sizes. Several parametric studies are performed to demonstrate that the method correctly reproduces data and trends reported in available literature. Moreover, the model is based exclusively on parameters with clear physical or geometrical meaning and as such it provides a convenient framework for its further experimental validation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Compressing Random Microstructures via Stochastic Wang Tilings

This paper presents a stochastic Wang tiling based technique to compress or reconstruct disordered microstructures on the basis of given spatial statistics. Unlike the existing approaches based on a single unit cell, it utilizes a finite set of tiles assembled by a stochastic tiling algorithm, thereby allowing to accurately reproduce long-range orientation orders in a computationally efficient manner. Although the basic features of the method are demonstrated for a two-dimensional particulate suspension, the present framework is fully extensible to generic multi-dimensional media.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Homogenization of plain weave composites with imperfect microstructure: Part II--Analysis of real-world materials

A two-layer statistically equivalent periodic unit cell is offered to predict a macroscopic response of plain weave multilayer carbon-carbon textile composites. Falling-short in describing the most typical geometrical imperfections of these material systems the original formulation presented in (Zeman and Šejnoha, International Journal of Solids and Structures, 41 (2004), pp. 6549--6571) is substantially modified, now allowing for nesting and mutual shift of individual layers of textile fabric in all three directions. Yet, the most valuable asset of the present formulation is seen in the possibility of reflecting the influence of negligible meso-scale porosity through a system of oblate spheroidal voids introduced in between the two layers of the unit cell. Numerical predictions of both the effective thermal conductivities and elastic stiffnesses and their comparison with available laboratory data and the results derived using the Mori-Tanaka averaging scheme support credibility of the present approach, about as much as the reliability of local mechanical properties found from nanoindentation tests performed directly on the analyzed composite samples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Localization analysis of variationally based gradient plasticity model

The paper presents analytical or semi-analytical solutions for the formation and evolution of localized plastic zone in a uniaxially loaded bar with variable cross-sectional area. A variationally based formulation of explicit gradient plasticity with linear softening is used, and the ensuing jump conditions and boundary conditions are discussed. Three cases with different regularity of the stress distribution are considered, and the problem is converted to a dimensionless form. Relations linking the load level, size of the plastic zone, distribution of plastic strain and plastic elongation of the bar are derived and compared to another, previously analyzed gradient formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Selected topics in homogenization of transport processes in historical masonry structures

The paper reviews several topics associated with the homogenization of transport processed in historical masonry structures. Since these often experience an irregular or random pattern, we open the subject by summarizing essential steps in the formulation of a suitable computational model in the form of Statistically Equivalent Periodic Unit Cell (SEPUC). Accepting SEPUC as a reliable representative volume element is supported by application of the Fast Fourier Transform to both the SEPUC and large binary sample of real masonry in search for effective thermal conductivities limited here to a steady state heat conduction problem. Fully coupled non-stationary heat and moisture transport is addressed next in the framework of two-scale first-order homogenization approach with emphases on the application of boundary and initial conditions on the meso-scale.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Simple Numerical Model of Laminated Glass Beams

This contribution presents a simple Finite Element model aimed at efficient simulation of layered glass units. The adopted approach is based on considering independent kinematics of each layer, tied together via Lagrange multipliers. Validation and verification of the resulting model against independent data demonstrate its accuracy, showing its potential for generalization towards more complex problems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Analysis of coupled transport phenomena in concrete at elevated temperatures

In this paper, we study a non-linear numerical scheme arising from the implicit time discretization of the Bažant-Thonguthai model for hygro-thermal behavior of concrete at high temperatures. Existence and uniqueness of the time-discrete solution in two dimensions is established using the theory of pseudomonotone operators in Banach spaces. Next, the spatial discretization is accomplished by the conforming finite element method. An illustrative numerical example shows that the numerical model reproduces well the rapid increase of pore pressure in wet concrete due to extreme heating. Such phenomenon is of particular interest for the safety assessment of concrete structures prone to thermally-induced spalling.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mori-Tanaka Based Estimates of Effective Thermal Conductivity of Various Engineering Materials

The purpose of this paper is to present a simple micromechanics-based model to estimate the effective thermal conductivity of real-world macroscopically isotropic materials of matrix-inclusion type. The methodology is based on the well-established Mori-Tanaka method for composite media reinforced with ellipsoidal inclusions, extended to account for imperfect thermal contact at the matrix-inclusion interface, random orientation of particles and particle size distribution. Using simple ensemble averaging arguments, we show that the original Mori-Tanaka relations are still applicable for these complex systems, provided that the inclusion conductivity is appropriately modified. Such conclusion is supported by the verification of the model against a detailed finite-element study as well as its validation against experimental data for a wide range of engineering material systems.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Some properties of strong solutions to nonlinear heat and moisture transport in multi-layer porous structures

The present paper deals with mathematical models of heat and moisture transport in layered building envelopes. The study of such processes generates a system of two doubly nonlinear evolution partial differential equations with appropriate initial and boundary conditions. The existence of the strong solution in two dimensions on a (short) time interval is proven. The proof rests on regularity results for elliptic transmission problem for isotropic composite-like materials.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-local energetics of random heterogeneous lattices

In this paper, we study the mechanics of statistically non-uniform two-phase elastic discrete structures. In particular, following the methodology proposed in (Luciano and Willis, Journal of the Mechanics and Physics of Solids 53, 1505-1522, 2005), energetic bounds and estimates of the Hashin-Shtrikman-Willis type are developed for discrete systems with a heterogeneity distribution quantified by second-order spatial statistics. As illustrated by three numerical case studies, the resulting expressions for the ensemble average of the potential energy are fully explicit, computationally feasible and free of adjustable parameters. Moreover, the comparison with reference Monte-Carlo simulations confirms a notable improvement in accuracy with respect to approaches based solely on the first-order statistics.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Macroscopic constitutive law for Mastic Asphalt Mixtures from multiscale modeling

A well established framework of an uncoupled hierarchical modeling approach is adopted here for the prediction of macroscopic material parameters of the Generalized Leonov (GL) constitutive model intended for the analysis of flexible pavements at both moderate and elevated temperature regimes. To that end, a recently introduced concept of a statistically equivalent periodic unit cell (SEPUC) is addressed to reflect a real microstructure of Mastic Asphalt mixtures (MAm). While mastic properties are derived from an extensive experimental program, the macroscopic properties of MAm are fitted to virtual numerical experiments performed on the basis of first order homogenization scheme. To enhance feasibility of the solution of the underlying nonlinear problem a two-step homogenization procedure is proposed. Here, the effective material properties are first found for a mortar phase, a composite consisting of a mastic matrix and a fraction of small aggregates. These properties are then introduced in place of the matrix in actual unit cells to give estimates of the model parameters on macroscale. Comparison with the Mori-Tanaka predictions is also provided suggesting limitations of classical micromechanical models.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Softening Gradient Plasticity: Analytical Study of Localization under Nonuniform Stress

Localization of plastic strain induced by softening can be objectively described by a regularized plasticity model that postulates a dependence of the current yield stress on a nonlocal softening variable defined by a differential (gradient) expression. This paper presents analytical solutions of the one-dimensional localization problem under certain special nonuniform stress distributions. The one-dimensional problem can be interpreted as describing either a tensile bar with variable cross section, or a beam subjected to a nonuniform bending moment. Explicit as well as implicit gradient formulations are considered. The evolution of the plastic strain profile and the shape of the load-displacement diagram are investigated. It is shown that even if the local constitutive law exhibits softening right from the onset of yielding, the global load-displacement diagram has a hardening part. The interplay between the internal length scales characterizing the material and the geometry is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Microstructure-based modeling of elastic functionally graded materials: One dimensional case

Functionally graded materials (FGMs) are two-phase composites with continuously changing microstructure adapted to performance requirements. Traditionally, the overall behavior of FGMs has been determined using local averaging techniques or a given smooth variation of material properties. Although these models are computationally efficient, their validity and accuracy remain questionable, since a link with the underlying microstructure (including its randomness) is not clear. In this paper, we propose a modeling strategy for the linear elastic analysis of FGMs systematically based on a realistic microstructural model. The overall response of FGMs is addressed in the framework of stochastic Hashin-Shtrikman variational principles. To allow for the analysis of finite bodies, recently introduced discretization schemes based on the Finite Element Method and the Boundary Element Method are employed to obtain statistics of local fields. Representative numerical examples are presented to compare the performance and accuracy of both schemes. To gain insight into similarities and differences between these methods and to minimize technicalities, the analysis is performed in the one-dimensional setting.

0 citations0 reviewsNo code linkedOpen access