Direct method-based statistical limit analysis of wc-co composites

In this paper, a direct method-based prediction of load-bearing capacity of nonperiodic WC-Co composites is presented. The main goal is to generalize the methodology of limit analysis on periodic heterogeneous media to materials with random microstructures. For such materials, the admissible macroscopic loading domains demonstrate remarkable scatter among RVE models of identical size and constituents but different morphologies. Limit analysis is performed on samples of a group of RVE models converted automatically from scanning electron microscopy (SEM) images. The corresponding admissible loading domains are numerically determined and statistically interpreted. The obtained results for plastic limit loads by direct method are compared with those from conventional incremental analysis

A hysteretic multiscale formulation for nonlinear dynamic analysis of composite materials

This article has been made available through the Brunel Open Access Publishing Fund.A new multiscale finite element formulation is presented for nonlinear dynamic analysis of heterogeneous structures. The proposed multiscale approach utilizes the hysteretic finite element method to model the microstructure. Using the proposed computational scheme, the micro-basis functions, that are used to map the microdisplacement components to the coarse mesh, are only evaluated once and remain constant throughout the analysis procedure. This is accomplished by treating inelasticity at the micro-elemental level through properly defined hysteretic evolution equations. Two types of imposed boundary conditions are considered for the derivation of the multiscale basis functions, namely the linear and periodic boundary conditions. The validity of the proposed formulation as well as its computational efficiency are verified through illustrative numerical experiments

Dynamic Response of Tunable Phononic Crystals and New Homogenization Approaches in Magnetoactive Composites

This research investigates dynamic response of tunable periodic structures and homogenization methods in magnetoelastic composites (MECs). The research on tunable periodic structures is focused on the design, modeling and understanding of wave propagation phenomena and the dynamic response of smart phononic crystals. High amplitude wrinkle formation is employed to study a one-dimensional phononic crystal slab consists of a thin film bonded to a thick compliant substrate. Buckling induced surface instability generates a wrinkly structure triggered by a compressive strain. It is demonstrated that surface periodic pattern and the corresponding large deformation can control elastic wave propagation in the low thickness composite slab. Simulation results show that the periodic wrinkly structure can be used as a smart phononic crystal which can switch band diagrams of the structure in a transformative manner. A magnetoactive phononic crystal is proposed which its dynamic properties are controlled by combined effects of large deformations and an applied magnetic field. Finite deformations and magnetic induction influence phononic characteristics of the periodic structure through geometrical pattern transformation and material properties. A magnetoelastic energy function is proposed to develop constitutive laws considering large deformations and magnetic induction in the periodic structure. Analytical and finite element methods are utilized to compute dispersion relation and band structure of the phononic crystal for different cases of deformation and magnetic loadings. It is demonstrated that magnetic induction not only controls the band diagram of the structure but also has a strong effect on preferential directions of wave propagation. Moreover, a thermally controlled phononic crystal is designed using ligaments of bi-materials in the structure.Comment: PhD mechanical engineering, University of Nevada, Reno (2015

Numerical modelling based on the multiscale homogenization theory. Application in composite materials and structures

A multi-domain homogenization method is proposed and developed in this thesis based on a two-scale technique. The method is capable of analyzing composite structures with several periodic distributions by partitioning the entire domain of the composite into substructures making use of the classical homogenization theory following a first-order standard continuum mechanics formulation. The need to develop the multi-domain homogenization method arose because current homogenization methods are based on the assumption that the entire domain of the composite is represented by one periodic or quasi-periodic distribution. However, in some cases the structure or composite may be formed by more than one type of periodic domain distribution, making the existing homogenization techniques not suitable to analyze this type of cases in which more than one recurrent configuration appears. The theoretical principles used in the multi-domain homogenization method were applied to assemble a computational tool based on two nested boundary value problems represented by a finite element code in two scales: a) one global scale, which treats the composite as an homogeneous material and deals with the boundary conditions, the loads applied and the different periodic (or quasi-periodic) subdomains that may exist in the composite; and b) one local scale, which obtains the homogenized response of the representative volume element or unit cell, that deals with the geometry distribution and with the material properties of the constituents. The method is based on the local periodicity hypothesis arising from the periodicity of the internal structure of the composite. The numerical implementation of the restrictions on the displacements and forces corresponding to the degrees of freedom of the domain's boundary derived from the periodicity was performed by means of the Lagrange multipliers method. The formulation included a method to compute the homogenized non-linear tangent constitutive tensor once the threshold of nonlinearity of any of the unit cells has been surpassed. The procedure is based in performing a numerical derivation applying a perturbation technique. The tangent constitutive tensor is computed for each load increment and for each iteration of the analysis once the structure has entered in the non-linear range. The perturbation method was applied at the global and local scales in order to analyze the performance of the method at both scales. A simple average method of the constitutive tensors of the elements of the cell was also explored for comparison purposes. A parallelization process was implemented on the multi-domain homogenization method in order to speed-up the computational process due to the huge computational cost that the nested incremental-iterative solution embraces. The effect of softening in two-scale homogenization was investigated following a smeared cracked approach. Mesh objectivity was discussed first within the classical one-scale FE formulation and then the concepts exposed were extrapolated into the two-scale homogenization framework. The importance of the element characteristic length in a multi-scale analysis was highlighted in the computation of the specific dissipated energy when strain-softening occurs. Various examples were presented to evaluate and explore the capabilities of the computational approach developed in this research. Several aspects were studied, such as analyzing different composite arrangements that include different types of materials, composites that present softening after the yield point is reached (e.g. damage and plasticity) and composites with zones that present high strain gradients. The examples were carried out in composites with one and with several periodic domains using different unit cell configurations. The examples are compared to benchmark solutions obtained with the classical one-scale FE method.En esta tesis se propone y desarrolla un método de homogeneización multi-dominio basado en una técnica en dos escalas. El método es capaz de analizar estructuras de materiales compuestos con varias distribuciones periódicas dentro de un mismo continuo mediante la partición de todo el dominio del material compuesto en subestructuras utilizando la teoría clásica de homogeneización a través de una formulación estándar de mecánica de medios continuos de primer orden. La necesidad de desarrollar este método multi-dominio surgió porque los métodos actuales de homogeneización se basan en el supuesto de que todo el dominio del material está representado por solo una distribución periódica o cuasi-periódica. Sin embargo, en algunos casos, la estructura puede estar formada por más de un tipo de distribución de dominio periódico. Los principios teóricos desarrollados en el método de homogeneización multi-dominio se aplicaron para ensamblar una herramienta computacional basada en dos problemas de valores de contorno anidados, los cuales son representados por un código de elementos finitos (FE) en dos escalas: a) una escala global, que trata el material compuesto como un material homogéneo. Esta escala se ocupa de las condiciones de contorno, las cargas aplicadas y los diferentes subdominios periódicos (o cuasi-periódicos) que puedan existir en el material compuesto; y b) una escala local, que obtiene la respuesta homogenizada de un volumen representativo o celda unitaria. Esta escala se ocupa de la geometría, y de la distribución espacial de los constituyentes del compuesto así como de sus propiedades constitutivas. El método se basa en la hipótesis de periodicidad local derivada de la periodicidad de la estructura interna del material. La implementación numérica de las restricciones de los desplazamientos y las fuerzas derivadas de la periodicidad se realizaron por medio del método de multiplicadores de Lagrange. La formulación incluye un método para calcular el tensor constitutivo tangente no-lineal homogeneizado una vez que el umbral de la no-linealidad de cualquiera de las celdas unitarias ha sido superado. El procedimiento se basa en llevar a cabo una derivación numérica aplicando una técnica de perturbación. El tensor constitutivo tangente se calcula para cada incremento de carga y para cada iteración del análisis una vez que la estructura ha entrado en el rango no-lineal. El método de perturbación se aplicó tanto en la escala global como en la local con el fin de analizar la efectividad del método en ambas escalas. Se lleva a cabo un proceso de paralelización en el método con el fin de acelerar el proceso de cómputo debido al enorme coste computacional que requiere la solución iterativa incremental anidada. Se investiga el efecto de ablandamiento por deformación en el material usando el método de homogeneización en dos escalas a través de un enfoque de fractura discreta. Se estudió la objetividad en el mallado dentro de la formulación clásica de FE en una escala y luego los conceptos expuestos se extrapolaron en el marco de la homogeneización de dos escalas. Se enfatiza la importancia de la longitud característica del elemento en un análisis multi-escala en el cálculo de la energía específica disipada cuando se produce el efecto de ablandamiento. Se presentan varios ejemplos para evaluar la propuesta computacional desarrollada en esta investigación. Se estudiaron diferentes configuraciones de compuestos que incluyen diferentes tipos de materiales, así como compuestos que presentan ablandamiento después de que el punto de fluencia del material se alcanza (usando daño y plasticidad) y compuestos con zonas que presentan altos gradientes de deformación. Los ejemplos se llevaron a cabo en materiales compuestos con uno y con varios dominios periódicos utilizando diferentes configuraciones de células unitarias. Los ejemplos se comparan con soluciones de referencia obtenidas con el método clásico de elementos finitos en una escala

A parallel elastic and inelastic heterogeneous multiscale method for rate-independent materials

In this thesis, we consider a heterogeneous multiscale method for elasticity and its extension to inelasticity within a two-scale energetic approach for rate-independent material models. For the efficiency, we propose a parallel realization of an FE² method which allows the computation of 3D microstructures with fine resolutions. The efficiency of this method is demonstrated for a simple damage model combined with elasto-plasticity describing a PBT matrix material with firmly embedded glass fiber inclusions

On the effective behavior of nonlinear inelastic composites: I. Incremental variational principles.

International audienceA new method for determining the overall behavior of composite materials comprised of nonlinear inelastic constituents is presented. Upon use of an implicit time-discretization scheme, the evolution equations describing the constitutive behavior of the phases can be reduced to the minimization of an incremental energy function. This minimization problem is rigorously equivalent to a nonlinear thermoelastic problem with a transformation strain which is a nonuniform field (not even uniform within the phases). In this first part of the study the variational technique of Ponte Castaneda is used to approximate the nonuniform eigenstrains by piecewise uniform eigenstrains and to linearize the nonlinear thermoelastic problem. The resulting problem is amenable to simpler calculations and analytical results for appropriate microstructures can be obtained. The accuracy of the proposed scheme is assessed by comparison of the method with exact results