A multiscale multisurface constitutive model for the thermo-plastic behavior of polyethylene

We present a multiscale model bridging length and time scales from molecular to continuum levels with the objective of predicting the yield behavior of amorphous glassy polyethylene (PE). Constitutive pa- rameters are obtained from molecular dynamics (MD) simulations, decreasing the requirement for ad- hoc experiments. Consequently, we achieve: (1) the identification of multisurface yield functions; (2) the high strain rate involved in MD simulations is upscaled to continuum via quasi-static simulations. Validation demonstrates that the entire multisurface yield functions can be scaled to quasi-static rates where the yield stresses are possibly predicted by a proposed scaling law; (3) a hierarchical multiscale model is constructed to predict temperature and strain rate dependent yield strength of the PE.

Effects of dynamic strain ageing on fatigue crack propagation and threshold behaviour in a Ni-base superalloy

In the temperature range 200-400 degree C the Ni-base superalloy, N901, develops marked dynamic strain ageing effects in its tensile behavior. These include inverse strain rate sensitivity, especially in UTS values, strongly serrated stress-strain curves and a heavily sheared failure mode at the higher test-temperatures. As for steels these effects seem to be due to interactions between the dislocations and the interstitial carbon atoms present. The results of tensile and fatigue threshold tests carried out between 20 degree C and 420 degree C are reported and the fatigue behavior is discussed in terms of the effects of surface roughness induced closure, temperature and strain aging interactions.

Studies on the fracture behaviour of polymer blends with special reference to PP/HDPE and PS/HIPS blends

Investigations on the fracture behaviour of polymer blends is the topic of this thesis. The blends selected are PP/HDPE and PS/HIPS. PP/HDPE blend is chosen due to its commercial importance and PS/HIPS blend is selected to study the transition from brittle fracture to ductile fracture.PP/HDPE blends were prepared at different compositions by melt blending at 180°C and fracture failure process was investigated by conducting notch sensitivity test and tensile test at different strain rates. The effects of two types of modifiers (particulate and elastomer) on the fracture behaviour and notch sensitivity of PP/HDPE blends were studied. The modifiers used are calcium carbonate, a hard particulate filler commonly used in plastics and Ethylene Propylene Diene Monomer (EPDM). They were added in 2%, 4% and 6% by weight of the blends.The study shows that the mechanical properties of PP/HDPE blends can be optimized by selecting proper blend compositions. The selected modifiers are found to alter and improve the fracture behaviour and notch sensitivity of the blends. Particulate fillers like calcium carbonate can be used for making the mechanical behaviour more stable at the various blend compositions. The resistance to notch sensitivity of the blends is found to be marginally lower in the presence of calcium carbonate. The elastomeric modifier EPDM produces a better stability of the mechanical behaviour. A low concentration of EPDM is sufficient to effect such a change. EPDM significantly improves the resistance to notch sensitivity of the blends. The study shows that judicious selection of modifiers can improve the fracture behaviour and notch sensitivity of PP/HDPE blends and help these materials to be used for critical applications.For investigating the transition in fracture behaviour and failure modes, PS/HIPS blends were selected. The blends were prepared by melt mixing followed by injection moulding to prepare the specimens for conducting tensile, impact and flexure tests. These tests were used to simulate the various conditions which promote failure.The tensile behaviour of unnotched and notched PS/HIPS blend samples were evaluated at slow speeds. Tensile strengths and moduli were found to increase at the higher testing speed for all the blend combinations whereas maximum strain at break was found to decrease. For a particular speed of testing, the tensile strength and modulus show only a very slight decrease as HIPS content is increased up to about 40%. However, there is a drastic decrease on increasing the HIPS content thereafter.The maximum strain at break shows only a very slight change up to about 40% HIPS content and thereafter shows a remarkable increase. The notched specimens also follow a comparable trend even though the notch sensitivity is seen high for PS rich blends containing up to 40% HIPS. The notch sensitivity marginally decreases with increase in HIPS content. At the same time, it is found to increase with the increase in strain rate. It is observed that blends containing more than 40% HIPS fail in ductile mode.The impact characteristics of PSIHIPS blends studied were impact strength, the energy absorbed by the test specimen and impact toughness. Remarkable increase in impact strength is observed as HIPS content in the blend exceeds 40%. The energy absorbed by the test specimens and the impact toughness also show a comparable trend.Flexural testing which helps to characterize the load bearing capacity was conducted on PS/HIPS blend samples at the two different testing speeds of 5mmlmin and 10 mm/min. The flexural strength increases with increase in testing speed for all the blend compositions. At both the speeds, remarkable reduction in flexural strength is observed as HIPS content in the blend exceeds 40%. The flexural strain and flexural energy absorbed by the specimens are found to increase with increase in HIPS content. At both the testing speeds, brittle fracture is observed for PS rich blends whereas HIPS rich blends show ductile mode of failure.Photoelastic investigations were conducted on PS/HIPS blend samples to analyze their failure modes. A plane polariscope with a broad source of light was utilized for the study. The coloured isochromatic fringes formed indicate the presence of residual stress concentration in the blend samples. The coverage made by the fringes on the test specimens varies with the blend composition and it shows a reducing trend with the increase in HIPS content. This indicates that the presence of residual stress is a contributing factor leading to brittle fracture in PS rich blends and this tendency gradually falls with increase in HIPS content and leads to their ductile mode of failure.

Life-cycle simulations of shallow frontal waves and the impact of deformation strain

The life-cycle of shallow frontal waves and the impact of deformation strain on their development is investigated using the idealised version of the Met Office non-hydrostatic Unified Model which includes the same physics and dynamics as the operational forecast model. Frontal wave development occurs in two stages; first, a deformation strain is applied to a front and a positive potential vorticity (PV) strip forms, generated by latent heat release in the frontal updraft; second, as the deformation strain is reduced the PV strip breaks up into individual anomalies. The circulations associated with the PV anomalies cause shallow frontal waves to form. The structure of the simulated frontal waves is consistent with the conceptual model of a frontal cyclone. Deeper frontal waves are simulated if the stability of the atmosphere is reduced. Deformation strain rates of different strengths are applied to the PV strip to determine whether a deformation strain threshold exists above which frontal wave development is suppressed. An objective method of frontal wave activity is defined and frontal wave development was found to be suppressed by deformation strain rates $\ge 0.4\times10^{-5}\mbox{s}^{-1}$. This value compares well with observed deformation strain rate thresholds and the analytical solution for the minimum deformation strain rate needed to suppress barotropic frontal wave development. The deformation strain rate threshold is dependent on the strength of the PV strip with strong PV strips able to overcome stronger deformation strain rates (leading to frontal wave development) than weaker PV strips.

Instability of surface-temperature filaments in strain and shear

The effects of uniform straining and shearing on the stability of a surface quasi-geostrophic temperature filament are investigated. Straining is shown to stabilize perturbations for wide filaments but only for a finite time until the filament thins to a critical width, after which some perturbations can grow. No filament can be stabilized in practice, since there are perturbations that can grow large for any strain rate. The optimally growing perturbations, defined as solutions that reach a certain threshold amplitude first, are found numerically for a wide range of parameter values. The radii of the vortices formed through nonlinear roll-up are found to be proportional to θ/s, where θ is the temperature anomaly of the filament and s the strain rate, and are not dependent on the initial size of the filament. Shearing is shown to reduce the normal-mode growth rates, but it cannot stabilize them completely when there are temperature discontinuities in the basic state; smooth filaments can be stabilized completely by shearing and a simple scaling argument provides the shear rate required. Copyright © 2010 Royal Meteorological Society

The stability of a two-dimensional vorticity filament under uniform strain

The quantitative effects of uniform strain and background rotation on the stability of a strip of constant vorticity (a simple shear layer) are examined. The thickness of the strip decreases in time under the strain, so it is necessary to formulate the linear stability analysis for a time-dependent basic flow. The results show that even a strain rate γ (scaled with the vorticity of the strip) as small as 0.25 suppresses the conventional Rayleigh shear instability mechanism, in the sense that the r.m.s. wave steepness cannot amplify by more than a certain factor, and must eventually decay. For γ < 0.25 the amplification factor increases as γ decreases; however, it is only 3 when γ e 0.065. Numerical simulations confirm the predictions of linear theory at small steepness and predict a threshold value necessary for the formation of coherent vortices. The results help to explain the impression from numerous simulations of two-dimensional turbulence reported in the literature that filaments of vorticity infrequently roll up into vortices. The stabilization effect may be expected to extend to two- and three-dimensional quasi-geostrophic flows.

Numerical prediction of three-dimensional time-dependent viscoelastic extrudate swell using differential and algebraic models

This study investigates the numerical simulation of three-dimensional time-dependent viscoelastic free surface flows using the Upper-Convected Maxwell (UCM) constitutive equation and an algebraic explicit model. This investigation was carried out to develop a simplified approach that can be applied to the extrudate swell problem. The relevant physics of this flow phenomenon is discussed in the paper and an algebraic model to predict the extrudate swell problem is presented. It is based on an explicit algebraic representation of the non-Newtonian extra-stress through a kinematic tensor formed with the scaled dyadic product of the velocity field. The elasticity of the fluid is governed by a single transport equation for a scalar quantity which has dimension of strain rate. Mass and momentum conservations, and the constitutive equation (UCM and algebraic model) were solved by a three-dimensional time-dependent finite difference method. The free surface of the fluid was modeled using a marker-and-cell approach. The algebraic model was validated by comparing the numerical predictions with analytic solutions for pipe flow. In comparison with the classical UCM model, one advantage of this approach is that computational workload is substantially reduced: the UCM model employs six differential equations while the algebraic model uses only one. The results showed stable flows with very large extrudate growths beyond those usually obtained with standard differential viscoelastic models. (C) 2010 Elsevier Ltd. All rights reserved.

Predicting the behaviour of structures under impact loads using geometrically distorted scaled models

When a scaled structure (model or replica) is used to predict the response of a full-size compound (prototype), the model geometric dimensions should relate to the corresponding prototype dimensions by a single scaling factor. However, owing to manufacturing technical restrictions, this condition cannot be accomplished for some of the dimensions in real structures. Accordingly, the distorted geometry will not comply with the overall geometric scaling factor, infringing the Pi theorem requirements for complete dynamic similarity. In the present study, a method which takes geometrical distortions into account is introduced, leading to a model similar to the prototype. As a means to infer the performance of this method, three analytical problems of structures subjected to dynamic loads are analysed. It is shown that the replica developed applying this technique is able to accurately predict the full-size structure behaviour even when the studied models have some of their dimensions severely distorted. (C) 2012 Elsevier Ltd. All rights reserved.

Mechanical and thermal behaviour of isotactic polypropylene reinforced with inorganic fullerene-like WS2 nanoparticles: Effect of filler loading and temperature

The thermal and mechanical behaviour of isotactic polypropylene (iPP) nanocomposites reinforced with different loadings of inorganic fullerene-like tungsten disulfide (IF-WS2) nanoparticles was investigated. The IF-WS2 noticeably enhanced the polymer stiffness and strength, ascribed to their uniform dispersion, the formation of a large nanoparticle?matrix interface combined with a nucleating effect on iPP crystallization. Their reinforcement effect was more pronounced at high temperatures. However, a drop in ductility and toughness was found at higher IF-WS2 concentrations. The tensile behaviour of the nanocomposites was extremely sensitive to the strain rate and temperature, and their yield strength was properly described by the Eyring s equation. The activation energy increased while the activation volume decreased with increasing nanoparticle loading, indicating a reduction in polymer chain motion. The nanoparticles improved the thermomechanical properties of iPP: raised the glass transition and heat deflection temperatures while decreased the coefficient of thermal expansion. The nanocomposites also displayed superior flame retardancy with longer ignition time and reduced peak heat release rate. Further, a gradual rise in thermal conductivity was found with increasing IF-WS2 loading both in the glassy and rubbery states. The results presented herein highlight the benefits and high potential of using IF-nanoparticles for enhancing the thermomechanical properties of thermoplastic polymers compared to other nanoscale fillers.

Multiscale modeling of the plastic behaviour in single crystal tungsten: from atomistic to crystal plasticity simulations

En una planta de fusión, los materiales en contacto con el plasma así como los materiales de primera pared experimentan condiciones particularmente hostiles al estar expuestos a altos flujos de partículas, neutrones y grandes cargas térmicas. Como consecuencia de estas diferentes y complejas condiciones de trabajo, el estudio, desarrollo y diseño de estos materiales es uno de los más importantes retos que ha surgido en los últimos años para la comunidad científica en el campo de los materiales y la energía. Debido a su baja tasa de erosión, alta resistencia al sputtering, alta conductividad térmica, muy alto punto de fusión y baja retención de tritio, el tungsteno (wolframio) es un importante candidato como material de primera pared y como posible material estructural avanzado en fusión por confinamiento magnético e inercial. Sin embargo, el tiempo de vida del tungsteno viene controlado por diversos factores como son su respuesta termo-mecánica en la superficie, la posibilidad de fusión y el fallo por acumulación de helio. Es por ello que el tiempo de vida limitado por la respuesta mecánica del tungsteno (W), y en particular su fragilidad, sean dos importantes aspectos que tienes que ser investigados. El comportamiento plástico en materiales refractarios con estructura cristalina cúbica centrada en las caras (bcc) como el tungsteno está gobernado por las dislocaciones de tipo tornillo a escala atómica y por conjuntos e interacciones de dislocaciones a escalas más grandes. El modelado de este complejo comportamiento requiere la aplicación de métodos capaces de resolver de forma rigurosa cada una de las escalas. El trabajo que se presenta en esta tesis propone un modelado multiescala que es capaz de dar respuestas ingenieriles a las solicitudes técnicas del tungsteno, y que a su vez está apoyado por la rigurosa física subyacente a extensas simulaciones atomísticas. En primer lugar, las propiedades estáticas y dinámicas de las dislocaciones de tipo tornillo en cinco potenciales interatómicos de tungsteno son comparadas, determinando cuáles de ellos garantizan una mayor fidelidad física y eficiencia computacional. Las grandes tasas de deformación asociadas a las técnicas de dinámica molecular hacen que las funciones de movilidad de las dislocaciones obtenidas no puedan ser utilizadas en los siguientes pasos del modelado multiescala. En este trabajo, proponemos dos métodos alternativos para obtener las funciones de movilidad de las dislocaciones: un modelo Monte Cario cinético y expresiones analíticas. El conjunto de parámetros necesarios para formular el modelo de Monte Cario cinético y la ley de movilidad analítica son calculados atomísticamente. Estos parámetros incluyen, pero no se limitan a: la determinación de las entalpias y energías de formación de las parejas de escalones que forman las dislocaciones, la parametrización de los efectos de no Schmid característicos en materiales bcc,etc. Conociendo la ley de movilidad de las dislocaciones en función del esfuerzo aplicado y la temperatura, se introduce esta relación como ecuación de flujo dentro de un modelo de plasticidad cristalina. La predicción del modelo sobre la dependencia del límite de fluencia con la temperatura es validada experimentalmente con ensayos uniaxiales en tungsteno monocristalino. A continuación, se calcula el límite de fluencia al aplicar ensayos uniaxiales de tensión para un conjunto de orientaciones cristalográticas dentro del triángulo estándar variando la tasa de deformación y la temperatura de los ensayos. Finalmente, y con el objetivo de ser capaces de predecir una respuesta más dúctil del tungsteno para una variedad de estados de carga, se realizan ensayos biaxiales de tensión sobre algunas de las orientaciones cristalográficas ya estudiadas en función de la temperatura.-------------------------------------------------------------------------ABSTRACT ----------------------------------------------------------Tungsten and tungsten alloys are being considered as leading candidates for structural and functional materials in future fusion energy devices. The most attractive properties of tungsten for the design of magnetic and inertial fusion energy reactors are its high melting point, high thermal conductivity, low sputtering yield and low longterm disposal radioactive footprint. However, tungsten also presents a very low fracture toughness, mostly associated with inter-granular failure and bulk plasticity, that limits its applications. As a result of these various and complex conditions of work, the study, development and design of these materials is one of the most important challenges that have emerged in recent years to the scientific community in the field of materials for energy applications. The plastic behavior of body-centered cubic (bcc) refractory metals like tungsten is governed by the kink-pair mediated thermally activated motion of h¿ (\1 11)i screw dislocations on the atomistic scale and by ensembles and interactions of dislocations at larger scales. Modeling this complex behavior requires the application of methods capable of resolving rigorously each relevant scale. The work presented in this thesis proposes a multiscale model approach that gives engineering-level responses to the technical specifications required for the use of tungsten in fusion energy reactors, and it is also supported by the rigorous underlying physics of extensive atomistic simulations. First, the static and dynamic properties of screw dislocations in five interatomic potentials for tungsten are compared, determining which of these ensure greater physical fidelity and computational efficiency. The large strain rates associated with molecular dynamics techniques make the dislocation mobility functions obtained not suitable to be used in the next steps of the multiscale model. Therefore, it is necessary to employ mobility laws obtained from a different method. In this work, we suggest two alternative methods to get the dislocation mobility functions: a kinetic Monte Carlo model and analytical expressions. The set of parameters needed to formulate the kinetic Monte Carlo model and the analytical mobility law are calculated atomistically. These parameters include, but are not limited to: enthalpy and energy barriers of kink-pairs as a function of the stress, width of the kink-pairs, non-Schmid effects ( both twinning-antitwinning asymmetry and non-glide stresses), etc. The function relating dislocation velocity with applied stress and temperature is used as the main source of constitutive information into a dislocation-based crystal plasticity framework. We validate the dependence of the yield strength with the temperature predicted by the model against existing experimental data of tensile tests in singlecrystal tungsten, with excellent agreement between the simulations and the measured data. We then extend the model to a number of crystallographic orientations uniformly distributed in the standard triangle and study the effects of temperature and strain rate. Finally, we perform biaxial tensile tests and provide the yield surface as a function of the temperature for some of the crystallographic orientations explored in the uniaxial tensile tests.

Yield behaviour of a melt-compounded polyethylene-intercalated montmorillonite nanocomposite

The yield behaviour of a series of melt-mixed polyethylene-modified montmorillonite nanocomposites has been studied as a function of temperature and strain rate and compared to the behaviour of the base polymer. The processing conditions used gave an intercalated structure as assessed by X-ray diffraction. Although there was a modest improvement in stiffness with clay content, the yield behaviour was insensitive to the addition of the clay. Both the base polymer and the nanocomposites showed double yield points. These were analysed as activated rate processes, with the activation energies consistent with the low strain yield point being associated with the alpha(2) molecular relaxation and the higher strain yield point with W axis slip. (C) 2003 Society of Chemical Industry.

The behaviour of granular material in pure shear, direct shear and simple shear

In biaxial compression tests, the stress calculations based on boundary information underestimate the principal stresses leading to a significant overestimation of the shear strength. In direct shear tests, the shear strain becomes highly concentrated in the mid-plane of the sample during the test. Although the stress distribution within the specimen is heterogeneous, the evolution of the stress ratio inside the shear band is similar to that inferred from the boundary force calculations. It is also demonstrated that the dilatancy in the shear band significantly exceeds that implied from the boundary displacements. In simple shear tests, the stresses acting on the wall boundaries do not reflect the internal state of stress but merely provide information about the average mobilised wall friction. It is demonstrated that the results are sensitive to the initial stress state defined by K0 = sh/sv. For all cases, non-coaxiality of the principal stress and strain-rate directions is examined and the corresponding flow rule is identified. Periodic cell simulations have been used to examine biaxial compression for a wide range of initial packing densities. Both constant volume and constant mean stress tests have been simulated. The characteristic behaviour at both the macroscopic and microscopic scales is determined by whether or not the system percolates (enduring connectivity is established in all directions). The transition from non-percolating to percolating systems is characterised by transitional behaviour of internal variables and corresponds to an elastic percolation threshold, which correlates well with the establishment of a mechanical coordination number of ca. 3.0. Strong correlations are found between macroscopic and internal variables at the critical state.

A study of the micro-mechanics of granular media

The development of more realistic constitutive models for granular media, such as sand, requires ingredients which take into account the internal micro-mechanical response to deformation. Unfortunately, at present, very little is known about these mechanisms and therefore it is instructive to find out more about the internal nature of granular samples by conducting suitable tests. In contrast to physical testing the method of investigation used in this study employs the Distinct Element Method. This is a computer based, iterative, time-dependent technique that allows the deformation of granular assemblies to be numerically simulated. By making assumptions regarding contact stiffnesses each individual contact force can be measured and by resolution particle centroid forces can be calculated. Then by dividing particle forces by their respective mass, particle centroid velocities and displacements are obtained by numerical integration. The Distinct Element Method is incorporated into a computer program 'Ball'. This program is effectively a numerical apparatus which forms a logical housing for this method and allows data input and output, and also provides testing control. By using this numerical apparatus tests have been carried out on disc assemblies and many new interesting observations regarding the micromechanical behaviour are revealed. In order to relate the observed microscopic mechanisms of deformation to the flow of the granular system two separate approaches have been used. Firstly a constitutive model has been developed which describes the yield function, flow rule and translation rule for regular assemblies of spheres and discs when subjected to coaxial deformation. Secondly statistical analyses have been carried out using data which was extracted from the simulation tests. These analyses define and quantify granular structure and then show how the force and velocity distributions use the structure to produce the corresponding stress and strain-rate tensors.

Fracture behaviour of an aluminide coating on a single crystal nickel base superalloy

Measurement of the coating fracture strain of an aluminide coating on a single crystal nickel base superalloy has been performed both in three-point bending and using variable wall thickness testpieces. As-aged specimens, 28 to 33 μm in thickness, were tested at room temperature, 600, 700 and 750 °C; specimens pre-exposed for 140 h at 850 and 1100 °C in air and vacuum were tested at room temperature. Fracture strains varied from 0.52 to 0.70% for as-aged specimens tested at temperatures up to 700 °C and specimens exposed at 850 °C and tested at room temperature. The crack path for these conditions was intergranular or transgranular in the main coating, along carbide-matrix interfaces in the coating transition zone, and at an angle of 30-45° to the original crack path in the substrate. The as-aged coating tested at 750 °C was ductile; a ductile-brittle transition occurs between 700 and 750 °C for the strain rate used (1 × 10-5 s-1). Following 1100 °C pre-exposure, specimens were ductile at room temperature with fractures strains of several percent. In this condition the crack morphology changed to one of subsurface nucleation in β grains and at β-γ′ interfaces. © 1993.