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

Falk K. Wittel

Falk K. Wittel 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-scicond-mat.softphysics.comp-phcond-mat.stat-mechBiological Physicsphysics.flu-dynphysics.geo-phcond-mat.dis-nncond-mat.otherNumerical AnalysisComputational Engineering, Finance, and Sciencemath.NAnlin.AOphysics.med-ph

Catalog footprint

What is connected

35works
14topics
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

35 published item(s)

preprint2016arXiv

Evolution of Strength and Failure of SCC during Early Hydration

The early strength evolution of self-consolidating concrete (SCC) is studied by a set of non-standard mechanical tests for compressive, tensile, shear and bending failure. The results are applicable in an industrial environment for process control, e.g. of slip casting with adaptive molds in robotic fabrication. A procedure for collapsing data to a master evolution curve is presented that allows to distinguish two regimes in the evolution. In the first, the material is capable of undergoing large localized plastic deformation, as expected from thixotropic yield stress fluids. This is followed by a transition to cohesive frictional material behavior dominated by crack growth. The typical differences in tensile and compressive strength of hardened concrete are observed to originate at the transition. Finally, the evolution of a limit surface in principal stress space is constructed and discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Micro-mechanical Failure Analysis of Wet Granular Matter

We employ a novel fluid-particle model to study the shearing behavior of granular soils under different saturation levels, ranging from the dry material via the capillary bridge regime to higher saturation levels with percolating clusters. The full complexity of possible liquid morphologies is taken into account, implying the formation of isolated arbitrary-sized liquid clusters with individual Laplace pressures that evolve by liquid exchange via films on the grain surface. Liquid clusters can grow in size, shrink, merge and split, depending on local conditions, changes of accessible liquid and the pore space morphology determined by the granular phase. This phase is represented by a discrete particle model based on Contact Dynamics, where capillary forces exerted from a liquid phase add to the motion of spherical particles. We study the macroscopic response of the system due to an external compression force at various liquid contents with the help of triaxial shear tests. Additionally, the change in liquid cluster distributions during the compression due to the deformation of the pore space is evaluated close to the critical load.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Computational up-scaling of anisotropic swelling and mechanical behavior of hierarchical cellular material

The hygro-mechanical behavior of a hierarchical cellular material, i.e. growth rings of softwood is investigated using a two-scale micro-mechanics model based on a computational homogenization technique. The lower scale considers the individual wood cells of varying geometry and dimensions. Honeycomb unit cells with periodic boundary conditions are utilized to calculate the mechanical properties and swelling coefficients of wood cells. Using the cellular scale results, the anisotropy in mechanical and swelling behavior of a growth ring in transverse directions is investigated. Predicted results are found to be comparable to experimental data. It is found that the orthotropic swelling properties of the cell wall in thin-walled earlywood cells produce anisotropic swelling behavior while, in thick latewood cells, this anisotropy vanishes. The proposed approach provides the ability to consider the complex microstructure when predicting the effective mechanical and swelling properties of softwood.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Continuous Wire Reinforcement for Jammed Granular Architecture

The mechanical behavior of continuous fiber reinforced granular columns is simulated by means of a Discrete Element Model. Spherical particles are randomly deposited simultaneously with a wire, that is deployed following different patterns inside of a flexible cylinder for triaxial compression testing. We quantify the effect of three different fiber deployment patterns on the failure envelope, represented by Mohr-Coulomb cones, and derive suggestions for improved deployment strategies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Contraction gradient induced microcracking in hardened cement paste

Drying induced cracking of concrete surfaces and repair layers is a common problem. A principal cause for this type of cracking is the moisture and resulting contraction gradient that develops in the cement paste matrix upon drying. This phenomenon has been experimentally quantified in unconfined hardened cement paste samples using a fluorescent resin impregnation technique. The effects of sample thickness and drying method on surface crack density and crack penetration depth are reported and explained. Finite element modelling of moisture gradients indicate the important role of the film coefficient in desiccation cracking of unconfined samples. The critical thickness for samples to remain crack-free upon drying was in the range of 2-5 mm depending on drying method. In thicker samples a crack spacing doubling process was observed that is in agreement with theoretical predictions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Discrete Element Simulation of Transverse Cracking During the Pyrolysis of Carbon Fibre Reinforced Plastics to Carbon/Carbon Composites

The fracture behavior of fiber-ceramics like C/C-SiC strongly depends on the initial damage arising during the production process. We study the transverse cracking of the 90° ply in [0/90]S cross-ply laminates due to the thermochemical degradation of the matrix material during the carbonization process by means of a discrete element method. The crack morphology strongly depends on the fiber-matrix interface properties, the transverse ply thickness as well as on the carbonization process itself. To model the 90° ply a two-dimensional triangular lattice of springs is constructed where nodes of the lattice represent fibers. Springs with random breaking thresholds model the disordered matrix material and interfaces. The spring-lattice is coupled by interface springs to two rigid bars which capture the two 0° plies or adjacent sublaminates in the model. Molecular dynamics simulation is used to follow the time evolution of the model system. It was found that under gradual heating of the specimen, after some distributed cracking, segmentation cracks occur in the 90° ply which then develop into a saturated state where the ply cannot support additional load. The dependence of the micro-structure of damage on the ply thickness and on the disorder in spring properties is also studied. Crack density and porosity of the system are monitored as a function of the temperature and compared to an analytic approach and experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

From fracture to fragmentation: discrete element modeling -- Complexity of crackling noise and fragmentation phenomena revealed by discrete element simulations

Discrete element modelling (DEM) is one of the most efficient computational approaches to the fracture processes of heterogeneous materials on mesoscopic scales. From the dynamics of single crack propagation through the statistics of crack ensembles to the rapid fragmentation of materials DEM had a substantial contribution to our understanding over the past decades. Recently, the combination of DEM with other simulation techniques like Finite Element Modelling further extended the field of applicability. In this paper we briefly review the motivations and basic idea behind the DEM approach to cohesive particulate matter and then we give an overview of on-going developments and applications of the method focusing on two fields where recent success has been achieved. We discuss current challenges of this rapidly evolving field and outline possible future perspectives and debates.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Full-field moisture induced deformation in Norway spruce: intra-ring variation of transverse swelling

The transverse hygro-expansion of the Norway spruce wood is studied on the growth ring level using digital image correlation. This non-destructive technique offers the possibility to contactless study deformation fields of relatively large areas. The measured full-field strains are segmented into individual growth rings. Whereas radial strains closely follow the density progression with the maximum in the dense latewood, tangential and shear strain remain constant except for positions around the edges of the sample. A simple FE three phase growth ring model is in good agreement with the experimental values. The selective activation of individual phases like earlywood, transitionwood and latewood demonstrates that the radial hygro-expansion is dominated by the earlywood deformation, whereas tangential deformation is a complex interplay of expansion and compression that needs all tissues to fully develop.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Generic failure mechanisms in adhesive bonds

The failure of adhesive bondlines has been studied at the microscopic level via tensile tests. Stable crack propagation could be generated by means of samples with improved geometry, which made in-situ observations possible. The interaction of cracks with adhesive bondlines under various angles to the crack propagation was the focus of this study as well as the respective loading situations for the adhesives UF, PUR, and PVAc, which have distinctly different mechanical behaviors. It is shown how adhesive properties influence the occurrence of certain failure mechanisms and determine their appearance and order of magnitude. With the observed failure mechanisms, it becomes possible to predict the propagation path of a crack through the specimen.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Grain Scale Modeling of Arbitrary Fluid Saturation in Random Packings

We propose a model for increasing liquid saturation in a granular packing which can account for liquid redistribution at saturation levels beyond the well-studied capillary bridge regime. The model is capable of resolving and combining capillary bridges, menisci and fully saturated pores to form local liquid clusters of any shape. They can exchange volume due to the local Laplace pressure gradient via a liquid film on the surfaces of grains. Local instabilities like Haines jumps trigger the discontinuous evolution of the liquid front. The applicability of the model is demonstrated and compared to benchmark experiments on the level of individual liquid structures as well as on larger systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Hierarchical structures for a robustness-oriented capacity design

In this paper, we study the response of 2D framed structures made of rectangular cells, to the sudden removal of columns. We employ a simulation algorithm based on the Discrete Element Method, where the structural elements are represented by elasto-plastic Euler Bernoulli beams with elongation-rotation failure threshold. The effect of structural cell slenderness and of topological hierarchy on the dynamic residual strength after damage $\ROne$ is investigated. Topologically \textit{hierarchical} frames have a primary structure made of few massive elements, while \textit{homogeneous} frames are made of many thin elements. We also show how $\ROne$ depends on the activated collapse mechanisms, which are determined by the mechanical hierarchy between beams and columns, i.e. by their relative strength and stiffness. Finally, principles of robustness-oriented capacity design which seem to be in contrast to the conventional anti-seismic capacity design are addressed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Invasion Percolation with a Hardening Interface under Gravity

We propose a modified Invasion Percolation (IP) model to simulate the infiltration of glue into a porous medium under gravity in 2D. Initially, the medium is saturated with air and then invaded by a fluid that has a hardening effect taking place from the interface towards the interior by contact with the air. To take into account that interfacial hardening, we use an IP model where capillary pressures of the growth sites are increased with time. In our model, if a site stays for a certain time at interface, it becomes a hard site and cannot be invaded anymore. That represents the glue interface becoming hard due to exposition with the air. Buoyancy forces are included in this system through the Bond number which represents the competition between the hydrostatic and capillary forces. We then compare our results with results from literature of non-hardening fluids in each regime of Bond number. We see that the invasion patterns change strongly with hardening while the non-hardening behavior remains basically not affected.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Inverse determination of effective mechanical properties of adhesive bondlines

A new approach for determining effective mechanical bondline properties using a combined experimental-numerical modal analysis technique is proposed. After characterizing clear spruce wood boards, an adhesive layer is applied to the boards surfaces. The shift of the eigenfrequencies resulting from the adhesive layer together with information on the bondline geometry can then be used to inversely determine the mechanical properties of the adhesive layer using Finite Element Models. The calculated values for clear wood as well as for the adhesive layer lie within reasonable ranges, thus demonstrating the methods potential.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Lattice-Boltzmann Method for Geophysical Plastic Flows

We explore possible applications of the Lattice-Boltzmann Method for the simulation of geophysical flows. This fluid solver, while successful in other fields, is still rarely used for geotechnical applications. We show how the standard method can be modified to represent free-surface realization of mudflows, debris flows, and in general any plastic flow, through the implementation of a Bingham constitutive model. The chapter is completed by an example of a full-scale simulation of a plastic fluid flowing down an inclined channel and depositing on a flat surface. An application is given, where the fluid interacts with a vertical obstacle in the channel.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mechanical Properties of Norway Spruce: Intra-Ring Variation and Generic Behavior of Earlywood and Latewood until Failure

The alternating earlywood and latewood growth ring structure has a strong influence on the mechanical performance of Norway spruce. In the current study, tensile tests in the longitudinal and tangential directions were performed on a series of specimens representing one growth ring at varying relative humidities. All tested mechanical parameters, namely modulus of elasticity and ultimate tensile stress, followed the density distribution in the growth ring, with the minimum values in earlywood and the maximum values in latewood. The samples were conditioned at three the relative humidities 50%, 65% and 95%. With increasing relative humidity, the values of the mechanical parameters were found to decrease. However, due to the high local variability, this decrease was not statistically significant. The test in the tangential direction on a set of earlywood and latewood specimens at 65% relative humidity revealed a similar limit of linear elasticity for both early- and latewood. Where the strength of both tissues was equal, the strain at failure was significantly greater for earlywood. Furthermore, the portion of the non-linear stress/strain behavior for earlywood was significantly greater. A Weibull analysis on the ultimate tensile strength revealed a tissue-independent Weibull modulus, which indicates similar defect distributions. For both, the failure occurred in the middle lamella.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Mechanisms in Impact Fragmentation

The brittle fragmentation of spheres is studied numerically by a 3D Discrete Element Model. Large scale computer simulations are performed with models that consist of agglomerates of many spherical particles, interconnected by beam-truss elements. We focus on a detailed description of the fragmentation process and study several fragmentation mechanisms involved. The evolution of meridional cracks is studied in detail. These cracks are found to initiate in the inside of the specimen with quasi-periodic angular distribution and give a broad peak in the fragment mass distribution for large fragments that can be fitted by a two-parameter Weibull distribution. The results prove to be independent of the degree of disorder in the model, but mean fragment sizes scale with velocity. Our results reproduce many experimental observations of fragment shapes, impact energy dependence or mass distribution, and significantly improve the understanding of the fragmentation process for impact fracture since we have full access to the failure conditions and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Moisture-induced Damage Evolution in Laminated Beech

Structural elements made of laminated hardwood are increasingly used in timber engineering. In this combined numerical and experimental approach, damage onset and propagation in uni-directional and cross-laminated samples out of European beech due to climatic changes are studied. The inter- and intra-laminar damage evolution is characterized for various configurations adhesively bonded by three structural adhesive systems. Typical situations are simulated by means of a comprehensive moisture-dependent non-linear rheological finite element model for wood with the capability to capture delaminations. The simulations give insight into the role of different strain components such as visco-elastic, mechano-sorptive, plastic, and hygro-elastic deformations under changing moisture content in progressive damage and delamination. We show the stress buildup under cyclic hygric loading resulting in hygro-fatigue and modify an analytical micro-mechanics of damage model, originally developed for cross-ply laminates, to describe the problem of moisture-induced damage in beech lamellae.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Multiscale Modeling of Polymer Gels-Chemo-Electric Model versus Discrete Element Model

Polyelectrolyte gels are a very attractive class of actuation materials with remarkable electronic and mechanical properties with a great similarity to biological contractile tissues. They consist of a polymer network with ionizable groups and a liquid phase with mobile ions. Absorption and delivery of solvent lead to a large change of volume. This mechanism can be triggered by chemical (change of salt concentration or pH of solution surrounding the gel), electrical, thermal or optical stimuli. Due to this capability, these gels can be used as actuators for technical applications, where large swelling and shrinkage is desired. In the present work chemically stimulated polymer gels in a solution bath are investigated. To adequately describe the different complicated phenomena occurring in these gels, they can be modeled on different scales. Therefore, models based on the statistical theory and porous media theory, as well as a coupled multi-field model and a discrete element formulation are derived and employed. In this paper, the coupled multi-field model and the discrete element model for chemical stimulation of a polymer gel film with and without domain deformation are employed. Based on these results, the presented formulations are compared and conclusions on their applicability in engineering practice are finally drawn.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Orthotropic rotation-free thin shell elements

A method to simulate orthotropic behaviour in thin shell finite elements is proposed. The approach is based on the transformation of shape function derivatives, resulting in a new orthogonal basis aligned to a specified preferred direction for all elements. This transformation is carried out solely in the undeformed state leaving minimal additional impact on the computational effort expended to simulate orthotropic materials compared to isotropic, resulting in a straightforward and highly efficient implementation. This method is implemented for rotation-free triangular shells using the finite element framework built on the Kirchhoff--Love theory employing subdivision surfaces. The accuracy of this approach is demonstrated using the deformation of a pinched hemispherical shell (with a 18° hole) standard benchmark. To showcase the efficiency of this implementation, the wrinkling of orthotropic sheets under shear displacement is analyzed. It is found that orthotropic subdivision shells are able to capture the wrinkling behavior of sheets accurately for coarse meshes without the use of an additional wrinkling model.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Packing of elastic wires in flexible shells

The packing problem of long thin filaments that are injected into confined spaces is of fundamental interest for physicists and biologists alike. How linear threads pack and coil is well known only for the ideal case of rigid containers, though. Here, we force long elastic rods into flexible spatial confinement borne by an elastic shell to examine under which conditions recently acquired knowledge on wire packing in rigid spheres breaks down. We find that unlike in rigid cavities, friction plays a key role by giving rise to the emergence of two distinct packing patterns. At low friction, the wire densely coils into an ordered toroidal bundle with semi-ellipsoidal cross-section, while at high friction, it packs into a highly disordered, hierarchic structure. These two morphologies are shown to be separated by a continuous phase transition. Our findings demonstrate the dramatic impact of friction and confinement elasticity on filamentous packing and might drive future research on such systems in physics, biology and even medical technology toward including these mutually interacting effects.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pore space analysis of beech wood: the vessel network

Water transport in wood is vital for the survival of trees. With synchrotron radiation X-ray tomographic microscopy (SRXTM), it becomes possible to characterize and quantify the 3D network formed by vessels that are responsible for longitudinal transport. In the present paper, the spatial size dependence of vessels and the organization inside single growth rings in terms of vessel induced porosity was studied by SRXTM. Network characteristics, such as connectivity, were deduced by digital image analysis from the processed tomographic data and related to known complex network topologies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Progressive Collapse Mechanisms of Brittle and Ductile Framed Structures

In this paper, we study the progressive collapse of 3D framed structures made of reinforced concrete after the sudden loss of a column. The structures are represented by elasto-plastic Euler Bernoulli beams with elongation-rotation failure threshold. We performed simulations using the Discrete Element Method considering inelastic collisions between the structural elements. The results show what collapse initiation and impact-driven propagation mechanisms are activated in structures with different geometric and mechanical features. Namely, we investigate the influence of the cross sectional size and reinforcement $α$ and of the plastic capacity $β$ of the structural elements. We also study the final collapse extent and the fragment size distribution and their relation to $α$, $β$ and to the observed collapse mechanisms. Finally, we compare the damage response of structures with symmetric and asymmetric reinforcement in the beams.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Simulating Thin Sheets: Buckling, Wrinkling, Folding and Growth

Numerical simulations of thin sheets undergoing large deformations are computationally challenging. Depending on the scenario, they may spontaneously buckle, wrinkle, fold, or crumple. Nature's thin tissues often experience significant anisotropic growth, which can act as the driving force for such instabilities. We use a recently developed finite element model to simulate the rich variety of nonlinear responses of Kirchhoff-Love sheets. The model uses subdivision surface shape functions in order to guarantee convergence of the method, and to allow a finite element description of anisotropically growing sheets in the classical Rayleigh-Ritz formalism. We illustrate the great potential in this approach by simulating the inflation of airbags, the buckling of a stretched cylinder, as well as the formation and scaling of wrinkles at free boundaries of growing sheets. Finally, we compare the folding of spatially confined sheets subject to growth and shrinking confinement to find that the two processes are equivalent.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Single Particle Fragmentation in Ultrasound Assisted Impact Comminution

Impact fragmentation is the underlying principle of comminution milling of dry, bulk solids. Unfortunately the outcome of the fragmentation process is more or less determined by the dimensionality of the impactor and its impact velocity. Since fragmentation is dominated by interfering shock waves, manipulating traveling shock waves and adding energy to the system during its fragmentation could be a promising approach to manipulate fragment mass distributions and energy input. In a former study we explored mechanisms in impact fragmentation of spheres, using a three-dimensional Discrete Element Model (DEM). This work is focused on studying how single spheres fragment when impacted on a planar vibrating target.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Study on the fragmentation of shells

Fragmentation can be observed in nature and in everyday life on a wide range of length scales and for all kinds of technical applications. Most studies on dynamic failure focus on the behaviour of bulk systems in one, two and three dimensions under impact and explosive loading, showing universal power law behaviour of fragment size distribution. However, hardly any studies have been devoted to fragmentation of shells. We present a detailed experimental and theoretical study on the fragmentation of closed thin shells of various materials, due to an excess load inside the system and impact with a hard wall. Characteristic fragmentation mechanisms are identified by means of a high speed camera and fragment shapes and mass distributions are evaluated. Theoretical rationalisation is given by means of stochastic break-up models and large-scale discrete element simulations with spherical shell systems under different extreme loading situations. By this we ex-plain fragment shapes and distributions and prove a power law for the fragment mass distribution. Satisfactory agreement between experimental findings and nu-merical predictions of the exponents of the power laws for the fragment shapes is obtained.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The role of water in the behavior of wood

Wood, due to its biological origin, has the capacity to interact with water. Sorption/desorption of moisture is accompanied with swelling/shrinkage and softening/hardening of its stiffness. The correct prediction of the behavior of wood components undergoing environmental loading requires that the moisture behavior and mechanical behavior of wood are considered in a coupled manner. We propose a comprehensive framework using a fully coupled poromechanical approach, where its multiscale implementation provides the capacity to take into account, directly, the exact geometry of the wood cellular structure, using computational homogenization. A hierarchical model is used to take into account the subcellular composite-like organization of the material. Such advanced modeling requires high resolution experimental data for the appropriate determination of inputs and for its validation. High-resolution x-ray tomography, digital image correlation, and neutron imaging are presented as valuable methods to provide the required information.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fragmentation processes in two-phase materials

We investigate the fragmentation process of solid materials with crystalline and amorphous phases using the discrete element method. Damage initiates inside spherical samples above the contact zone in a region where the circumferential stress field is tensile. Cracks initiated in this region grow to form meridional planes. If the collision energy exceeds a critical value which depends on the material's internal structure, cracks reach the sample surface resulting in fragmentation. We show that this primary fragmentation mechanism is very robust with respect to the internal structure of the material. For all configurations, a sharp transition from the damage to the fragmentation regime is observed, with smaller critical collision energies for crystalline samples. The mass distribution of the fragments follows a power law for small fragments with an exponent that is characteristic for the branching merging process of unstable cracks. Moreover this exponent depends only on the dimensionally of the system and not on the micro structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Morphogenesis of filaments growing in flexible confinements

Space-saving design is a requirement that is encountered in biological systems and the development of modern technological devices alike. Many living organisms dynamically pack their polymer chains, filaments or membranes inside of deformable vesicles or soft tissue like cell walls, chorions, and buds. Surprisingly little is known about morphogenesis due to growth in flexible confinements - perhaps owing to the daunting complexity lying in the nonlinear feedback between packed material and expandable cavity. Here we show by experiments and simulations how geometric and material properties lead to a plethora of morphologies when elastic filaments are growing far beyond the equilibrium size of a flexible thin sheet they are confined in. Depending on friction, sheet flexibility and thickness, we identify four distinct morphological phases emerging from bifurcation and present the corresponding phase diagram. Four order parameters quantifying the transitions between these phases are proposed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multi-scale approach to invasion percolation of rock fracture networks

A multi-scale scheme for the invasion percolation of rock fracture networks with heterogeneous fracture aperture fields is proposed. Inside fractures, fluid transport is calculated on the finest scale and found to be localized in channels as a consequence of the aperture field. The channel network is characterized and reduced to a vectorized artificial channel network (ACN). Different realizations of ACNs are used to systematically calculate efficient apertures for fluid transport inside differently sized fractures as well as fracture intersection and entry properties. Typical situations in fracture networks are parameterized by fracture inclination, flow path length along the fracture and intersection lengths in the entrance and outlet zones of fractures. Using these scaling relations obtained from the finer scales, we simulate the invasion process of immiscible fluids into saturated discrete fracture networks, which were studied in previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Particle-fluid-structure interaction for debris flow impact on flexible barriers

Flexible barriers are increasingly used for the protection from debris flow in mountainous terrain due to their low cost and environmental impact. However, a numerical tool for rational design of such structures is still missing. In this work, a hybrid computational framework is presented, using a total Lagrangian formulation of the Finite Element Method (FEM) to represent a flexible barrier. The actions exerted on the structure by a debris flow are obtained from simultaneous simulations of the flow of a fluid-grain mixture, using two conveniently coupled solvers: the Discrete Element Method (DEM) governs the motion of the grains, while the free-surface non-Newtonian fluid phase is solved using the Lattice-Boltzmann Method (LBM). Simulations on realistic geometries show the dependence of the momentum transfer on the barrier on the composition of the debris flow, challenging typical assumptions made during the design process today. In particular, we demonstrate that both grains and fluid contribute in a non-negligible way to the momentum transfer. Moreover, we show how the flexibility of the barrier reduces its vulnerability to structural collapse, and how the stress is distributed on its fabric, highlighting potential weak points.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rheological Model for Wood

Wood as the most important natural and renewable building material plays an important role in the construction sector. Nevertheless, its hygroscopic character basically affects all related mechanical properties leading to degradation of material stiffness and strength over the service life. Accordingly, to attain reliable design of the timber structures, the influence of moisture evolution and the role of time- and moisture-dependent behaviors have to be taken into account. For this purpose, in the current study a 3D orthotropic elasto-plastic, visco-elastic, mechano-sorptive constitutive model for wood, with all material constants being defined as a function of moisture content, is presented. The corresponding numerical integration approach, with additive decomposition of the total strain is developed and implemented within the framework of the finite element method (FEM). Moreover to preserve a quadratic rate of asymptotic convergence the consistent tangent operator for the whole model is derived. Functionality and capability of the presented material model are evaluated by performing several numerical verification simulations of wood components under different combinations of mechanical loading and moisture variation. Additionally, the flexibility and universality of the introduced model to predict the mechanical behavior of different species are demonstrated by the analysis of a hybrid wood element. Furthermore, the proposed numerical approach is validated by comparisons of computational evaluations with experimental results.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Subdivision Shell Elements with Anisotropic Growth

A thin shell finite element approach based on Loop's subdivision surfaces is proposed, capable of dealing with large deformations and anisotropic growth. To this end, the Kirchhoff-Love theory of thin shells is derived and extended to allow for arbitrary in-plane growth. The simplicity and computational efficiency of the subdivision thin shell elements is outstanding, which is demonstrated on a few standard loading benchmarks. With this powerful tool at hand, we demonstrate the broad range of possible applications by numerical solution of several growth scenarios, ranging from the uniform growth of a sphere, to boundary instabilities induced by large anisotropic growth. Finally, it is shown that the problem of a slowly and uniformly growing sheet confined in a fixed hollow sphere is equivalent to the inverse process where a sheet of fixed size is slowly crumpled in a shrinking hollow sphere in the frictionless, quasi-static, elastic limit.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Finite Element Simulation of Dense Wire Packings

A finite element program is presented to simulate the process of packing and coiling elastic wires in two- and three-dimensional confining cavities. The wire is represented by third order beam elements and embedded into a corotational formulation to capture the geometric nonlinearity resulting from large rotations and deformations. The hyperbolic equations of motion are integrated in time using two different integration methods from the Newmark family: an implicit iterative Newton-Raphson line search solver, and an explicit predictor-corrector scheme, both with adaptive time stepping. These two approaches reveal fundamentally different suitability for the problem of strongly self-interacting bodies found in densely packed cavities. Generalizing the spherical confinement symmetry investigated in recent studies, the packing of a wire in hard ellipsoidal cavities is simulated in the frictionless elastic limit. Evidence is given that packings in oblate spheroids and scalene ellipsoids are energetically preferred to spheres.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

From invasion percolation to flow in rock fracture networks

The main purpose of this work is to simulate two-phase flow in the form of immiscible displacement through anisotropic, three-dimensional (3D) discrete fracture networks (DFN). The considered DFNs are artificially generated, based on a general distribution function or are conditioned on measured data from deep geological investigations. We introduce several modifications to the invasion percolation (MIP) to incorporate fracture inclinations, intersection lines, as well as the hydraulic path length inside the fractures. Additionally a trapping algorithm is implemented that forbids any advance of the invading fluid into a region, where the defending fluid is completely encircled by the invader and has no escape route. We study invasion, saturation, and flow through artificial fracture networks, with varying anisotropy and size and finally compare our findings to well studied, conditioned fracture networks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Self-contact and instabilities in the anisotropic growth of elastic membranes

We investigate the morphology of thin discs and rings growing in the circumferential direction. Recent analytical results suggest that this growth produces symmetric excess cones (e-cones). We study the stability of such solutions considering self-contact and bending stress. We show that, contrary to what was assumed in previous analytical solutions, beyond a critical growth factor, no symmetric \textit{e}-cone solution is energetically minimal any more. Instead, we obtain skewed e-cone solutions having lower energy, characterized by a skewness angle and repetitive spiral winding with increasing growth. These results are generalized to discs with varying thickness and rings with holes of different radii.

0 citations0 reviewsNo code linkedOpen access