Thermodynamic properties and atomistic structure of the dry amorphous silica surface from a reactive force field model

A force field model of the Keating type supplemented by rules to break, form, and interchange bonds is applied to investigate thermodynamic and structural properties of the amorphous SiO2 surface. A simulated quench from the liquid phase has been carried out for a silica sample made of 3888 silicon and 7776 oxygen atoms arranged on a slab similar to 40 angstrom thick, periodically repeated along two directions. The quench results into an amorphous sample, exposing two parallel square surfaces of similar to 42 nm(2) area each. Thermal averages computed during the quench allow us to determine the surface thermodynamic properties as a function of temperature. The surface tension turns out to be gamma=310 +/- 20 erg/cm(2) at room temperature and gamma=270 +/- 30 at T=2000 K, in fair agreement with available experimental estimates. The entropy contribution Ts-s to the surface tension is relatively low at all temperatures, representing at most similar to 20% of the surface energy. Almost without exceptions, Si atoms are fourfold coordinated and oxygen atoms are twofold coordinated. Twofold and threefold rings appear only at low concentration and are preferentially found in proximity of the surface. Above the glass temperature T-g=1660 +/- 50 K, the mobility of surface atoms is, as expected, slightly higher than that of bulk atoms. The computation of the height-height correlation function shows that the silica surface is rough in the equilibrium and undercooled liquid phase, becoming smooth below the glass temperature T-g.

Numerical study of the effects of reinforcement/matrix interphase on stress-strain behavior of YAl2 particle reinforced MgLiAl composites

For the potential influence produced by the reinforcement/matrix interphase in particle reinforced metal matrix composites (PMMCs), a unit cell model with transition interphase was proposed. Uniaxial tensile loading was simulated and the stress/strain behavior was predicted. The results show that a transition interphase with both appropriate strength and thickness could affect the failure mode, reduce the stress concentration, and enhance the maximum strain value of the composite.

Measurement of fibre fracture toughness using an alternative specimen geometry

This paper presents an experimental and numerical study focused on the tensile fibre fracture toughness characterisation of hybrid plain weave composite laminates using non-standardized Overheight Compact Tension (OCT) specimens. The position as well as the strain field ahead of the crack tip in the specimens was determined using a digital speckle photogrammetry system. The limitation on the applicability of standard data reduction schemes for the determination of the intralaminar fibre fracture toughness of composites is presented and discussed. A methodology based on the numerical evaluation of the strain energy release rate using the J-integral method is proposed to derive new geometric correction functions for the determination of stress intensity factor for alternative composite specimen geometries. A comparison between different methods currently available to compute the intralaminar fracture toughness in composites is also presented and discussed. Good agreement between numerical and experimental results using the proposed methodology was obtained.

Intralaminar toughness characterisation of unbalanced hybrid plain weave laminates

<p>A numerical and experimental investigation on the mode-I intralaminar toughness of a hybrid plain weave composite laminate manufactured using resin infusion under flexible tooling (RIFT) process is presented in this paper. The pre-cracked geometries consisted of overheight compact tension (OCT), double edge notch (DEN) and centrally cracked four-point-bending (4PBT) test specimens. The position as well as the strain field ahead of the crack tip during the loading stage was determined using a digital speckle photogrammetry system. The limitation on the applicability of the standard data reduction schemes for the determination of intralaminar toughness of composite materials is presented and discussed. A methodology based on the numerical evaluation of the strain energy release rate using the J-integral method is proposed to derive new geometric correction functions for the determination of the stress intensity factor for composites. The method accounts for material anisotropy and finite specimen dimension effects regardless of the geometry. The approach has been validated for alternative non-standard specimen geometries. A comparison between different methods currently available for computing the intralaminar fracture toughness in composite laminates is presented and a good agreement between numerical and experimental results using the proposed methodology was obtained. </p>

Fracture mechanics using a 3D composite element

<p>The formulation of a 3D composite element and its use in a mixed-mode fracture mechanics example is presented. This element, like a conventional 3D finite element, has three degrees of freedom per node although, like a plate element, the strains are defined in the local directions of the mid-plane surface. The stress-strain property matrix of this element was modified to decouple the stresses in the local mid-plane and the strains normal to this plane thus preventing the element from being too stiff in bending. A main advantage of this formulation is the ability to model a laminate with a single 3D element. The motivation behind this work was to improve the computational efficiency associated with the calculation of strain energy release rates in laminated structures. A comparison of mixed-mode results using different elements of an in-house finite element package are presented. Good agreement was achieved between the results obtained using the new element and coventional higher-order elements</p>

Crack propagation in non-homogenous materials: Evaluation of mixed-mode SIFs, T-stress and kinking angle using a variant EFG method

A new variant of the Element-Free Galerkin (EFG) method, that combines the diffraction method, to characterize the crack tip solution, and the Heaviside enrichment function for representing discontinuity due to a crack, has been used to model crack propagation through non-homogenous materials. In the case of interface crack propagation, the kink angle is predicted by applying the maximum tangential principal stress (MTPS) criterion in conjunction with consideration of the energy release rate (ERR). The MTPS criterion is applied to the crack tip stress field described by both the stress intensity factor (SIF) and the T-stress, which are extracted using the interaction integral method. The proposed EFG method has been developed and applied for 2D case studies involving a crack in an orthotropic material, crack along an interface and a crack terminating at a bi-material interface, under mechanical or thermal loading; this is done to demonstrate the advantages and efficiency of the proposed methodology. The computed SIFs, T-stress and the predicted interface crack kink angles are compared with existing results in the literature and are found to be in good agreement. An example of crack growth through a particle-reinforced composite materials, which may involve crack meandering around the particle, is reported.

Estimation of the small strain stiffness of glacial till using geophysical methods and barometric loading response

Stiffness values in geotechnical structures can range over many orders of magnitude for relatively small operational strains. The typical strain levels where soil stiffness changes most dramatically is in the range 0.01-0.1%, however soils do not exhibit linear stress-strain behaviour at small strains. Knowledge of the in situ stiffness at small strain is important in geotechnical numerical modelling and design. The stress-strain regime of cut slopes is complex, as we have different principle stress directions at different positions along the potential failure plane. For example, loading may be primarily in extension near the toe of the slope, while compressive loading is predominant at the crest of a slope. Cuttings in heavily overconsolidated clays are known to be susceptible to progressive failure and subsequent strain softening, in which progressive yielding propagates from the toe towards the crest of the slope over time. In order to gain a better understanding of the rate of softening it would be advantageous to measure changes in small strain stiffness in the field.

Towards An Acoustic Model-Based Poroelastic Imaging Method: I. Theoretical Foundation

The ultrasonic measurement and imaging of tissue elasticity is currently under wide investigation and development as a clinical tool for the assessment of a broad range of diseases, but little account in this field has yet been taken of the fact that soft tissue is porous and contains mobile fluid. The ability to squeeze fluid out of tissue may have implications for conventional elasticity imaging, and may present opportunities for new investigative tools. When a homogeneous, isotropic, fluid-saturated poroelastic material with a linearly elastic solid phase and incompressible solid and fluid constituents is subjected to stress, the behaviour of the induced internal strain field is influenced by three material constants: the Young's modulus (E(s)) and Poisson's ratio (nu(s)) of the solid matrix and the permeability (k) of the solid matrix to the pore fluid. New analytical expressions were derived and used to model the time-dependent behaviour of the strain field inside simulated homogeneous cylindrical samples of such a poroelastic material undergoing sustained unconfined compression. A model-based reconstruction technique was developed to produce images of parameters related to the poroelastic material constants (E(s), nu(s), k) from a comparison of the measured and predicted time-dependent spatially varying radial strain. Tests of the method using simulated noisy strain data showed that it is capable of producing three unique parametric images: an image of the Poisson's ratio of the solid matrix, an image of the axial strain (which was not time-dependent subsequent to the application of the compression) and an image representing the product of the aggregate modulus E(s)(1-nu(s))/(1+nu(s))(1-2nu(s)) of the solid matrix and the permeability of the solid matrix to the pore fluid. The analytical expressions were further used to numerically validate a finite element model and to clarify previous work on poroelastography.

Co-melt fluidised bed granulation of pharmaceutical powders: Improvements in drug bioavailability

This study investigates the use of co-melt fluidised bed granulation for the agglomeration of model pharmaceutical powders, namely, lactose mono-hydrate, PEG 10000, poly-vinyl pyrolidone and ibuprofen as a model drug. Granulation within the co-melt system was found to follow a nucleationâ??steady growthâ??coating regime profile. Using high molecular weight PEG binder, the granulation mechanism and thus the extent of granulation was found to be significantly influenced by binder viscosity. The compression properties of the granulate within the hot fluidised bed were correlated using a novel high temperature experimental procedure. It was found that the fracture stress and fractural modulus of the materials under hot processing conditions were orders of magnitude lower than those measured under ambient conditions. A range of particle velocities within the granulator were considered based on theoretical models. After an initial period of nucleation, the Stokes deformation number analysis indicated that only velocities within the high shear region of the fluidised bed were sufficient to promote significant granule deformation and therefore, coalescence. The data also indicated that larger granules de-fluidised preventing agglomeration by coalescence. Furthermore, experimental data indicated that dissipation of the viscous molten binder to the surface was the most important factor in the latter stages of the granulation process. From a pharmaceutical perspective the inclusion of the model drug, ibuprofen, combined with PVP in the co-melt process proved to be highly significant. It was found that using DSC analysis on the formulations that the decrease in the heat of fusion associated with the melting of ibuprofen within the FHMG systems may be attributed to interaction between PVP and ibuprofen through inter-molecular hydrogen bonding. This interaction decreases the crystallinity of ibuprofen and facilitates solubilisation and bioavailability within the solid matrix.

Effects of sintering temperature on morphology and mechanical characteristics of 3D printed porous titanium used as dental implant

Porous titanium samples were manufactured using the 3D printing and sintering method in order to determine the effects of final sintering temperature on morphology and mechanical properties. Cylindrical samples were printed and split into groups according to a final sintering temperature (FST). Irregular geometry samples were also printed and split into groups according to their FST. The cylindrical samples were used to determine part shrinkage, in compressive tests to provide stress-strain data, in microCT scans to provide internal morphology data and for optical microscopy to determine surface morphology. All of the samples were used in microhardness testing to establish the hardness. Below 1100 C FST, shrinkage was in the region of 20% but increased to approximately 30% by a FST of 1300 C. Porosity varied from a maximum of approximately 65% at the surface to the region of 30% internally. Between 97 and 99% of the internal porosity is interconnected. Average pore size varied between 24 µm at the surface and 19 µm internally. Sample hardness increased to in excess of 300 HV0.05 with increasing FST while samples with an FST of below 1250 C produced an elastic-brittle stress/strain curve and samples above this displayed elastic-plastic behaviour. Yield strength increased significantly through the range of sintering temperatures while the Young's modulus remained fairly consistent. © 2013 Elsevier B.V.

Review of constitutive models for describing the time dependent behaviour of soft clays

Different classes of constitutive models have been proposed to capture the time-dependent behaviour of soft soil (creep, stress relaxation, rate dependency). This paper critically reviews many of the models developed based on understanding of the time dependent stress-strain-stress rate-strain rate behaviour of soils and viscoplasticity in terms of their strengths and weaknesses. Some discussion is also made on the numerical implementation aspects of these models. Typical findings from numerical analyses of geotechnical structures constructed on soft soils are also discussed. The general elastic viscoplastic (EVP) models can roughly be divided into two categories: models based on the concept of overstress and models based on non-stationary flow surface theory. Although general in structure, both categories have their own strengths and shortcomings. This review indicates that EVP analysis is yet to be vastly used by the geotechnical engineers, apparently due to the mathematical complication involved in the formulation of the constitutive models, unconvincing benefit in terms of the accuracy of performance prediction, requirement of additional soil parameter(s), difficulties in determining them, and the necessity of excessive computing resources and time. © 2013 Taylor &amp; Francis.

Three dimensional morphology and compressive behaviour of sintered biodegradable composite scaffolds

Porous poly-L-lactide acid (PLA) scaffolds are prepared using polymer sintering and porogen leaching method. Different weight fractions of the Hydroxyapatite (HA) are added to the PLA to control the acidity and degradation rate. The three dimensional morphology and surface porosity are tested using micro CT, optical microscopy and scanning electron microscopy (SEM). Results indicate that the surface porosity does not change by addition of HA. The micro Ct examinations show slight decrease in the pore size and increase in wall thickness accompanied with reduced anisotropy for the scaffolds containing HA. SEM micrographs show detectable interconnected pores for the scaffold with pure PLA. Addition of the HA results in agglomeration of the HA which blocks some of the pores. Compression tests of the scaffold identify three stages in the stress-strain curve. The addition of HA adversely affects the modulus of the scaffold at the first stage, but this was reversed for the second and third stages of the compression. The results of these tests are compared with the cellular material model. The manufactured scaffold have acceptable properties for a scaffold, however improvement to the mixing of the phases of PLA and HA is required to achieve better integrity of the composite scaffolds.

Morphology and compression behaviour of biodegradable scaffolds produced by the sintering process

Porous poly(L-lactic acid) (PLA) scaffolds of 85 per cent and 90 per cent porosity are prepared using polymer sintering and porogen leaching method. Different weight fractions of 10 per cent, 30 per cent, and 50 per cent of hydroxyapatite (HA) are added to the PLA to control the acidity and degradation rate. The three-dimensional (3D) morphology and surface porosity are tested using micro-computer tomography (micro-CT), optical microscopy, and scanning electron microscopy (SEM). Results indicate that the surface porosity does not change on the addition of HA. The micro-CT examinations show a slight decrease in the pore size and increase in the wall thickness accompanied by reduced anisotropy for the scaffolds containing HA. Scanning electron micrographs show detectable interconnected pores for the scaffold with pure PLA. Addition of the HA results in agglomeration of the HA particles and reduced leaching of the porogen. Compression tests of the scaffold identify three stages in the stress-strain curve. The addition of HA results in a reduction in the modulus of the scaffold at the first stage of elastic bending of the wall, but this is reversed for the second and third stages of collapse of the wall and densification in the compression tests. In the scaffolds with 85 per cent porosity, the addition of a high percentage of HA could result in 70 per cent decrease in stiffness in the first stage, 200 per cent increase in stiffness in the second stage, and 20 per cent increase in stiffness in the third stage. The results of these tests are compared with the Gibson cellular material model that is proposed for prediction of the behaviour of cellular material under compression. The pH and molecular weight changes are tracked for the scaffolds within a period of 35 days. The addition of HA keeps the pH in the alkaline region, which results in higher rate of degradation at an early period of observation, followed by a reduced rate of degradation later in the process. The final molecular weight is higher for the scaffolds with HA than for scaffolds of pure PLA. The manufactured scaffolds offer acceptable properties in terms of the pore size range and interconnectivity of the pores and porosity for non-load-bearing bone graft substitute; however, improvement to the mixing of the phases of PLA and HA is required to achieve better integrity of the composite scaffolds. © 2008 IMechE.

Finite element analysis of interfacial friction and slip in composites with fully unbonded filler

A finite element model is developed to predict the stress-strain behaviour of particulate composites with fully unbonded filler particles. This condition can occur because of the lack of adhesion property of the filler surface. Whilst part of the filler particle is separated from the matrix, another section of filler keeps in contact with the matrix because of the lateral compressive displacement of the matrix. The slip boundary condition is imposed on the section of the interface that remains closed. The states of stress and displacement fields are obtained. The location of any further deformation through crazing or shear band formation is identified. A completely unbonded inclusion with partial slip at a section of the interface reduces the concentration of the stress at the interface significantly. Whereas this might lead to slightly higher strength, it decreases the load transfer efficiency and stiffness of this type of composite.