Complex fiber bundle model for optimization of heterogeneous materials

We propose a complex fiber bundle model for the optimization of heterogeneous materials, which consists of many simple bundles. We also present an exact and compact recursion relation for the failure probability of a simple fiber bundle model with local load sharing, which is more efficient than the ones reported previously. Using a ''renormalization method'' and the recursion relation developed for the simple bundle, we calculate the failure probabilities of the complex fiber bundle. When the total number of fibers is given, we find that there exists an optimum way to organize the complex bundle, in which one gets a stronger bundle than in other ways.

Modeling mechanical response of heterogeneous materials

Heterogeneous materials are ubiquitous in nature and as synthetic materials. These materials provide unique combination of desirable mechanical properties emerging from its heterogeneities at different length scales. Future structural and technological applications will require the development of advanced light weight materials with superior strength and toughness. Cost effective design of the advanced high performance synthetic materials by tailoring their microstructure is the challenge facing the materials design community. Prior knowledge of structure-property relationships for these materials is imperative for optimal design. Thus, understanding such relationships for heterogeneous materials is of primary interest. Furthermore, computational burden is becoming critical concern in several areas of heterogeneous materials design. Therefore, computationally efficient and accurate predictive tools are highly essential. In the present study, we mainly focus on mechanical behavior of soft cellular materials and tough biological material such as mussel byssus thread. Cellular materials exhibit microstructural heterogeneity by interconnected network of same material phase. However, mussel byssus thread comprises of two distinct material phases. A robust numerical framework is developed to investigate the micromechanisms behind the macroscopic response of both of these materials. Using this framework, effect of microstuctural parameters has been addressed on the stress state of cellular specimens during split Hopkinson pressure bar test. A voronoi tessellation based algorithm has been developed to simulate the cellular microstructure. Micromechanisms (microinertia, microbuckling and microbending) governing macroscopic behavior of cellular solids are investigated thoroughly with respect to various microstructural and loading parameters. To understand the origin of high toughness of mussel byssus thread, a Genetic Algorithm (GA) based optimization framework has been developed. It is found that two different material phases (collagens) of mussel byssus thread are optimally distributed along the thread. These applications demonstrate that the presence of heterogeneity in the system demands high computational resources for simulation and modeling. Thus, Higher Dimensional Model Representation (HDMR) based surrogate modeling concept has been proposed to reduce computational complexity. The applicability of such methodology has been demonstrated in failure envelope construction and in multiscale finite element techniques. It is observed that surrogate based model can capture the behavior of complex material systems with sufficient accuracy. The computational algorithms presented in this thesis will further pave the way for accurate prediction of macroscopic deformation behavior of various class of advanced materials from their measurable microstructural features at a reasonable computational cost.

A double inclusion homogenization scheme for polycrystals with hierarchal topologies: Application to twinning in Mg alloys

The present work introduces a double inclusion elasto-plastic self-consistent (DI-EPSC) scheme for topologies in which crystals can contain subdomains (i.e. twins, etc.). The approach yields a direct coupling between the mechanical response of grains and their subdomains via a concentration relationship on mean fields derived from both the Eshelby and the Tanaka-Mori properties. The latent effect caused by twinning on the mechanical response is observed on both initially extruded and non-textured Mg alloys. For twinned grains, it is shown that deformation system activities and plastic strain distributions within twins drastically depend on the interaction with parent domains. Moreover, a quantitative study on the coupled influence of secondary slip activities on the material response is proposed. © 2014 Published by Elsevier Ltd.

Effects of Microstructural Heterogeneity on the Spallation Behavior of Materials

It is of utmost importance to understand the spallation behaviour of heterogeneous materials. In this paper, a driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. The spallation behavior of heterogeneity material is analyzed with this model. The heterogeniety of mesoscopic units is characterized in terms of Weibull modulus m of strength distibution and stress fluctuation parameter k. At high stress, the maximum damage increases with m; while at low stress, the maximum damage decreases. In addition, for low stress, severe stress fluctuation causes higher damage; while for high stress, causes lower damage.

Fracture statistics of brittle materials: Weibull or normal distribution

The fit of fracture strength data of brittle materials (Si3N4, SiC, and ZnO) to the Weibull and normal distributions is compared in terms of the Akaike information criterion. For Si3N4, the Weibull distribution fits the data better than the normal distribution, but for ZnO the result is just the opposite. In the case of SiC, the difference is not large enough to make a clear distinction between the two distributions. There is not sufficient evidence to show that the Weibull distribution is always preferred to other distributions, and the uncritical use of the Weibull distribution for strength data is questioned.

Effects of Microstructural Heterogeneity on The Spallation Behavior of Materials

It is of utmost importance to understand the spallation behaviour of heterogeneous materials. In this paper, a driven nonlinear threshold model with stress fluctuation is presented to study the effects of microstructural heterogeneity on continuum damage evolution. The spallation behavior of heterogeneity material is analyzed with this model. The heterogeniety of mesoscopic units is characterized in terms of Weibull modulus m of strength distibution and stress fluctuation parameter k. At high stress, the maximum damage increases with m; while at low stress, the maximum damage decreases. In addition, for low stress, severe stress fluctuation causes higher damage; while for high stress, causes lower damage.

Evolution of Localized Damage Zone in Heterogeneous Media

Evolution of localized damage zone is a key to catastrophic rupture in heterogeneous materials. In the present article, the evolutions of strain fields of rock specimens are investigated experimentally. The observed evolution of fluctuations and autocorrelations of strain fields under uniaxial compression demonstrates that the localization of deformation always appears ahead of catastrophic rupture. In particular, the localization evolves pronouncedly with increasing deformation in the rock experiments. By means of the definition of the zone with high strain rate and likely damage localization, it is found that the size of the localized zone decreases from the sample size at peak load to an eventual value. Actually, the deformation field beyond peak load is bound to suffer bifurcation, namely an elastic unloading part and a continuing but localized damage part will co-exist in series in a specimen. To describe this continuous bifurcation and localization process observed in experiments, a model on continuum mechanics is developed. The model can explain why the decreasing width of localized zone can lead stable deformation to unstable, but it still has not provided the complete equations governing the evolution of the localized zone.

An energy based nanomechnical properties evaluation method for cementitious materials

Advances in multiscale material modeling of structural concrete have created an upsurge of interest in the accurate evaluation of mechanical properties and volume fractions of its nano constituents. The task is accomplished by analyzing the response of a material to indentation, obtained as an outcome of a nanoindentation experiment, using a procedure called the Oliver and Pharr (OP) method. Despite its widespread use, the accuracy of this method is often questioned when it is applied to the data from heterogeneous materials or from the materials that show pile-up and sink-in during indentation, which necessitates the development of an alternative method. ^ In this study, a model is developed within the framework defined by contact mechanics to compute the nanomechanical properties of a material from its indentation response. Unlike the OP method, indentation energies are employed in the form of dimensionless constants to evaluate model parameters. Analysis of the load-displacement data pertaining to a wide range of materials revealed that the energy constants may be used to determine the indenter tip bluntness, hardness and initial unloading stiffness of the material. The proposed model has two main advantages: (1) it does not require the computation of the contact area, a source of error in the existing method; and (2) it incorporates the effect of peak indentation load, dwelling period and indenter tip bluntness on the measured mechanical properties explicitly. ^ Indentation tests are also carried out on samples from cement paste to validate the energy based model developed herein by determining the elastic modulus and hardness of different phases of the paste. As a consequence, it has been found that the model computes the mechanical properties in close agreement with that obtained by the OP method; a discrepancy, though insignificant, is observed more in the case of C-S-H than in the anhydrous phase. Nevertheless, the proposed method is computationally efficient, and thus it is highly suitable when the grid indentation technique is required to be performed. In addition, several empirical relations are developed that are found to be crucial in understanding the nanomechanical behavior of cementitious materials.^

Modified lattice model for mode-I fracture analysis of notched plain concrete beam using probabilistic approach

A modified lattice model using finite element method has been developed to study the mode-I fracture analysis of heterogeneous materials like concrete. In this model, the truss members always join at points where aggregates are located which are modeled as plane stress triangular elements. The truss members are given the properties of cement mortar matrix randomly, so as to represent the randomness of strength in concrete. It is widely accepted that the fracture of concrete structures should not be based on strength criterion alone, but should be coupled with energy criterion. Here, by incorporating the strain softening through a parameter ‘α’, the energy concept is introduced. The softening branch of load-displacement curves was successfully obtained. From the sensitivity study, it was observed that the maximum load of a beam is most sensitive to the tensile strength of mortar. It is seen that by varying the values of properties of mortar according to a normal random distribution, better results can be obtained for load-displacement diagram.

A Nonlinear Threshold Model Applied to Spallation Analysis

Spallation in heterogeneous media is a complex, dynamic process. Generally speaking, the spallation process is relevant to multiple scales and the diversity and coupling of physics at different scales present two fundamental difficulties for spallation modeling and simulation. More importantly, these difficulties can be greatly enhanced by the disordered heterogeneity on multi-scales. In this paper, a driven nonlinear threshold model for damage evolution in heterogeneous materials is presented and a trans-scale formulation of damage evolution is obtained. The damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that some dimensionless numbers govern the whole process of deformation and damage evolution. The effects of heterogeneity in terms of Weibull modulus on damage evolution in spallation process are also investigated.

A rigorous analytical method for doubly periodic cylindrical inclusions under longitudinal shear and its application

An infinite elastic solid containing a doubly periodic parallelogrammic array of cylindrical inclusions under longitudinal shear is studied. A rigorous and effective analytical method for exact solution is developed by using Eshelby's equivalent inclusion concept integrated with the new results from the doubly quasi-periodic Riemann boundary value problems. Numerical results show the dependence of the stress concentrations in such heterogeneous materials on the periodic microstructure parameters. The overall longitudinal shear modulus of composites with periodic distributed fibers is also studied. Several problems of practical importance, such as those of doubly periodic holes or rigid inclusions, singly periodic inclusions and single inclusion, are solved or resolved as special cases. The present method can provide benchmark results for other numerical and approximate methods. (C) 2003 Elsevier Ltd. All rights reserved.

非均质材料损伤破坏的数值模拟软件开发及强韧化机理研究

本文针对平面链网模型与空间链网模型，发展了模拟非均质材料的数值软件，旨在分析材料的损伤破坏特征，寻找材料细观结构同宏观力学行为之间的相互联系，探讨材料强韧化机理以及评估材料的力学性能。采用可视化技术设计的具有交互式功能的图象用户界面，能根据具体材料对象构造特定材料模型来分析其力学行为与服役能力。初步结果表明，无论是材料弹性阶段的性质还是损伤演化阶段材料中裂纹的萌生、扩展、汇集和材料的最终失效等现象都与已有的理论结果及试验结果符合得很好。另外，本文还对材料的强韧化机理作了探索。