Monte Carlo simulation of complex cohesive fracture in random heterogeneous quasi-brittle materials

A numerical method is developed to simulate complex two-dimensional crack propagation in quasi-brittle materials considering random heterogeneous fracture properties. Potential cracks are represented by pre-inserted cohesive elements with tension and shear softening constitutive laws modelled by spatially varying Weibull random fields. Monte Carlo simulations of a concrete specimen under uni-axial tension were carried out with extensive investigation of the effects of important numerical algorithms and material properties on numerical efficiency and stability, crack propagation processes and load-carrying capacities. It was found that the homogeneous model led to incorrect crack patterns and load–displacement curves with strong mesh-dependence, whereas the heterogeneous model predicted realistic, complicated fracture processes and load-carrying capacity of little mesh-dependence. Increasing the variance of the tensile strength random fields with increased heterogeneity led to reduction in the mean peak load and increase in the standard deviation. The developed method provides a simple but effective tool for assessment of structural reliability and calculation of characteristic material strength for structural design.

Evidence and Spatial Planning for Sustainable Development: A Model for Spatially Allocating Material Flows

The aspiration the spatial planning should act as the main coordinating function for the transition to a sustainable society is grounded on the assumption that it is capable of incorporating both a strong evidence base of environmental accounting for policy, coupled with opportunities for open, deliberative decision-making. While there are a number of increasingly sophisticated methods (such as material flow analysis and ecological footprinting) that can be used to longitudinally determine the impact of policy, there are fewer that can provide a robust spatial assessment of sustainability policy. In this paper, we introduce the Spatial Allocation of Material Flow Analysis (SAMFA) model, which uses the concept of socio-economic metabolism to extrapolate the impact of local consumption patterns that may occur as a result of the local spatial planning process at multiple spatial levels. The initial application the SAMFA model is based on County Kildare in the Republic of Ireland, through spatial temporal simulation and visualisation of construction material flows and associated energy use in the housing sector. Thus, while we focus on an Ireland case study, the model is applicable to spatial planning and sustainability research more generally. Through the development and evaluation of alternative scenarios, the model appears to be successful in its prediction of the cumulative resource and energy impacts arising from consumption and development patterns. This leads to some important insights in relation to the differential spatial distribution of disaggregated allocation of material balance and energy use, for example that rural areas have greater resource accumulation (and are therefore in a sense “less sustainable”) than urban areas, confirming that rural housing in Ireland is both more material and energy intensive. This therefore has the potential to identify hotspots of higher material and energy use, which can be addressed through targeted planning initiatives or focussed community engagement. Furthermore, due to the ability of the model to allow manipulation of different policy criteria (increased density, urban conservation etc), it can also act as an effective basis for multi-stakeholder engagement.

Strong electronic correlation in the hydrogen chain: A variational Monte Carlo study

In this paper, we report a fully ab initio variational Monte Carlo study of the linear and periodic chain of hydrogen atoms, a prototype system providing the simplest example of strong electronic correlation in low dimensions. In particular, we prove that numerical accuracy comparable to that of benchmark density-matrix renormalization-group calculations can be achieved by using a highly correlated Jastrow-antisymmetrized geminal power variational wave function. Furthermore, by using the so-called "modern theory of polarization" and by studying the spin-spin and dimer-dimer correlations functions, we have characterized in detail the crossover between the weakly and strongly correlated regimes of this atomic chain. Our results show that variational Monte Carlo provides an accurate and flexible alternative to highly correlated methods of quantum chemistry which, at variance with these methods, can be also applied to a strongly correlated solid in low dimensions close to a crossover or a phase transition.

1-D constitutive model for evolution of stress-induced R-phase and localized Lüders-like stress-induced martensitic transformation of super-elastic NiTi wires

NiTi alloys have been widely used in the applications for micro-electro-mechanical-systems (MEMS), which often involve some precise and complex motion control. However, when using the NiTi alloys in MEMS application, the main problem to be considered is the degradation of functional property during cycling loading. This also stresses the importance of accurate prediction of the functional behavior of NiTi alloys. In the last two decades, a large number of constitutive models have been proposed to achieve the task. A portion of them focused on the deformation behavior of NiTi alloys under cyclic loading, which is a practical and non-negligible situation. Despite of the scale of modeling studies of the field in NiTi alloys, two experimental observations under uniaxial tension loading have not received proper attentions. First, a deviation from linearity well before the stress-induced martensitic transformation (SIMT) has not been modeled. Recent experiments confirmed that it is caused by the formation of stress-induced R phase. Second, the influence of the well-known localized Lüders-like SIMT on the macroscopic behavior of NiTi alloys, in particular the residual strain during cyclic loading, has not been addressed. In response, we develop a 1-D phenomenological constitutive model for NiTi alloys with two novel features: the formation of stress-induced R phase and the explicit modeling of the localized Lüders-like SIMT. The derived constitutive relations are simple and at the same time sufficient to describe the behavior of NiTi alloys. The accumulation of residual strain caused by R phase under different loading schemes is accurately described by the proposed model. Also, the residual strain caused by irreversible SIMT at different maximum loading strain under cyclic tension loading in individual samples can be explained by and fitted into a single equation in the proposed model. These results show that the proposed model successfully captures the behavior of R phase and the essence of localized SIMT.

Global quantum correlations in finite-size spin chains

We perform an extensive study of the properties of global quantum correlations in finite-size one-dimensional quantum spin models at finite temperature. By adopting a recently proposed measure for global quantum correlations (Rulli and Sarandy 2011 Phys. Rev. A 84 042109), called global discord, we show that critical points can be neatly detected even for many-body systems that are not in their ground state. We consider the transverse Ising model, the cluster-Ising model where three-body couplings compete with an Ising-like interaction, and the nearest-neighbor XX Hamiltonian in transverse magnetic field. These models embody our canonical examples showing the sensitivity of global quantum discord close to criticality. For the Ising model, we find a universal scaling of global discord with the critical exponents pertaining to the Ising universality class.

Finite Element Analysis of Drilling of Carbon Fibre Reinforced Composites

Despite the increased applications of the composite materials in aerospace due to their exceptional physical and mechanical properties, the machining of composites remains a challenge. Fibre reinforced laminated composites are prone to different damages during machining process such as delamination, fibre pull-out, microcracks, thermal damages. Optimization of the drilling process parameters can reduces the probability of these damages. In the current research, a 3D finite element (FE) model is developed of the process of drilling in the carbon fibre reinforced composite (CFC). The FE model is used to investigate the effects of cutting speed and feed rate on thrust force, torque and delamination in the drilling of carbon fiber reinforced laminated composite. A mesoscale FE model taking into account of the different oriented plies and interfaces has been proposed to predict different damage modes in the plies and delamination. For validation purposes, experimental drilling tests have been performed and compared to the results of the finite element analysis. Using Matlab a digital image analysis code has been developed to assess the delamination factor produced in CFC as a result of drilling. © Springer Science+Business Media B.V. 2011.

Strudel: A corpus-based semantic model based on properties and types

Computational models of meaning trained on naturally occurring text successfully model human performance on tasks involving simple similarity measures, but they characterize meaning in terms of undifferentiated bags of words or topical dimensions. This has led some to question their psychological plausibility (Murphy, 2002; Schunn, 1999). We present here a fully automatic method for extracting a structured and comprehensive set of concept descriptions directly from an English part-of-speech-tagged corpus. Concepts are characterized by weighted properties, enriched with concept-property types that approximate classical relations such as hypernymy and function. Our model outperforms comparable algorithms in cognitive tasks pertaining not only to concept-internal structures (discovering properties of concepts, grouping properties by property type) but also to inter-concept relations (clustering into superordinates), suggesting the empirical validity of the property-based approach. Copyright © 2009 Cognitive Science Society, Inc. All rights reserved.

Physical and numerical constraints in source modeling for finite difference simulation of room acoustics

In finite difference time domain simulation of room acoustics, source functions are subject to various constraints. These depend on the way sources are injected into the grid and on the chosen parameters of the numerical scheme being used. This paper addresses the issue of selecting and designing sources for finite difference simulation, by first reviewing associated aims and constraints, and evaluating existing source models against these criteria. The process of exciting a model is generalized by introducing a system of three cascaded filters, respectively, characterizing the driving pulse, the source mechanics, and the injection of the resulting source function into the grid. It is shown that hard, soft, and transparent sources can be seen as special cases within this unified approach. Starting from the mechanics of a small pulsating sphere, a parametric source model is formulated by specifying suitable filters. This physically constrained source model is numerically consistent, does not scatter incoming waves, and is free from zero- and low-frequency artifacts. Simulation results are employed for comparison with existing source formulations in terms of meeting the spectral and temporal requirements on the outward propagating wave.

Criticality, factorization, and long-range correlations in the anisotropic XY model

We study the long-range quantum correlations in the anisotropic XY model. By first examining the thermodynamic limit, we show that employing the quantum discord as a figure of merit allows one to capture the main features of the model at zero temperature. Furthermore, by considering suitably large site separations we find that these correlations obey a simple scaling behavior for finite temperatures, allowing for efficient estimation of the critical point. We also address ground-state factorization of this model by explicitly considering finite-size systems, showing its relation to the energy spectrum and explaining the persistence of the phenomenon at finite temperatures. Finally, we compute the fidelity between finite and infinite systems in order to show that remarkably small system sizes can closely approximate the thermodynamic limit.

Effect of the crown design and interface lute parameters on the stress-state of a machined crown-tooth system: A finite element analysis

The effect of preparation design and the physical properties of the interface lute on the restored machined ceramic crown-tooth complex are poorly understood. The aim of this work was to determine, by means of three-dimensional finite element analysis (3D FEA) the effect of the tooth preparation design and the elastic modulus of the cement on the stress state of the cemented machined ceramic crown-tooth complex. The three-dimensional structure of human premolar teeth, restored with adhesively cemented machined ceramic crowns, was digitized with a micro-CT scanner. An accurate, high resolution, digital replica model of a restored tooth was created. Two preparation designs, with different occlusal morphologies, were modeled with cements of 3 different elastic moduli. Interactive medical image processing software (mimics and professional CAD modeling software) was used to create sophisticated digital models that included the supporting structures; periodontal ligament and alveolar bone. The generated models were imported into an FEA software program (hypermesh version 10.0, Altair Engineering Inc.) with all degrees of freedom constrained at the outer surface of the supporting cortical bone of the crown-tooth complex. Five different elastic moduli values were given to the adhesive cement interface 1.8 GPa, 4 GPa, 8 GPa, 18.3 GPa and 40 GPa; the four lower values are representative of currently used cementing lutes and 40 GPa is set as an extreme high value. The stress distribution under simulated applied loads was determined. The preparation design demonstrated an effect on the stress state of the restored tooth system. The cement elastic modulus affected the stress state in the cement and dentin structures but not in the crown, the pulp, the periodontal ligament or the cancellous and cortical bone. The results of this study suggest that both the choice of the preparation design and the cement elastic modulus can affect the stress state within the restored crown-tooth complex.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 1:Particulate filled composites

Micro-mechanical analysis of polymeric composites provides a powerful means for the quantitative assessment of their bulk behavior. In this paper we describe a robust finite element model (FEM) for the micro-structural modeling of the behavior of particulate filled polymer composites under external loads. The developed model is applied to simulate stress distribution in polymer composites containing particulate fillers. Quantitative information about the magnitude and location of maximum stress concentrations obtained from these simulations is used to predict the dominant failure and crack growth mechanisms in these composites. The model predictions are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show the range of the validity of the developed model and its predictive potential.

Prediction of failure and fracture mechanisms of polymeric composites using finite element analysis. Part 2:Fiber reinforced composites

A robust finite element scheme for the micro-mechanical modeling of the behavior of fiber reinforced polymeric composites under external loads is developed. The developed model is used to simulate stress distribution throughout the composite domain and to identify the locations where maximum stress concentrations occur. This information is used as a guide to predict dominant failure and crack growth mechanisms in fiber reinforced composites. The differences between continuous fibers, which are susceptible to unidirectional transverse fracture, and short fibers have been demonstrated. To assess the validity and range of applicability of the developed scheme, numerical results obtained by the model are compared with the available experimental data and also with the values found using other methods reported in the literature. These comparisons show that the present finite element scheme can generate meaningful results in the analysis of fiber reinforced composites.

Finite element analysis of drilling of titanium alloy

Drilling is a highly demanding machining process due to complex tool geometry and the progressive material failure on the work piece. In this study, a 3D model is developed using commercial finite element software ABAQUS/Explicit. The proposed model simulates the drilling process by taking into account of the damage initiation and evolution of the work piece material, a contact model at the interface between drill bit and work piece and the process parameters. The results of the simulations demonstrate the effects of machining parameters on drilling. The results also confirm the capability and advantage of FE simulation of the drilling process. © 2011 Published by Elsevier Ltd.

Finite element analysis of ultrasonic insertion of SiC fibre in aluminium alloy 6061

In this work, a computational framework has been proposed to successfully simulate the fibre embedding using ultrasonic consolidation process. The main components of the proposed computational approach are a developed constitutive model and a friction model which are especially suitable for the condition of ultrasonic process. The effect of different process parameters, such as velocity of sonotrode, displacement amplitude of ultrasonic vibration and applied loads are studied. The presented work especially focuses on the quality of the developed weld and the fibre coverage due to the plastic flow around the fibre. The areas of maximum plastic flow predicted by the simulation are confirmed by the EBSD microstructural studies. © 2011 Inderscience Enterprises Ltd.

An application of bi-directional evolutionary structural optimisation for optimising energy absorbing structures using a material damage model

The Bi-directional Evolutionary Structural Optimisation (BESO) method is a numerical topology optimisation method developed for use in finite element analysis. This paper presents a particular application of the BESO method to optimise the energy absorbing capability of metallic structures. The optimisation objective is to evolve a structural geometry of minimum mass while ensuring that the kinetic energy of an impacting projectile is reduced to a level which prevents perforation. Individual elements in a finite element mesh are deleted when a prescribed damage criterion is exceeded. An energy absorbing structure subjected to projectile impact will fail once the level of damage results in a critical perforation size. It is therefore necessary to constrain an optimisation algorithm from producing such candidate solutions. An algorithm to detect perforation was implemented within a BESO framework which incorporated a ductile material damage model.