On the initiation of boudinage and folding instabilities in layered elasto-visco-plastic solids - Insights from natural microstructures and numerical simulations

The present understanding of the initiation of boudinage and folding structures is based on viscosity contrasts and stress exponents, considering an intrinsically unstable state of the layer. The criterion of localization is believed to be prescribed by geometry-material interactions, which are often encountered in natural structures. An alternative localization phenomenon has been established for ductile materials, in which instability emerges for critical material parameters and loading rates from homogeneous conditions. In this thesis, conditions are sought under which this type of instability prevails and whether localization in geological materials necessarily requires a trigger by geometric imperfections. The relevance of critical deformation conditions, material parameters and the spatial configuration of instabilities are discussed in a geological context. In order to analyze boudinage geometries, a numerical eigenmode analysis is introduced. This method allows determining natural frequencies and wavelengths of a structure and inducing perturbations on these frequencies. In the subsequent coupled thermo-mechanical simulations, using a grain size evolution and end-member flow laws, localization emerges when material softening through grain size sensitive viscous creep sets in. Pinch-and-swell structures evolve along slip lines through a positive feedback between the matrix response and material bifurcations inside the layer, independent from the mesh-discretization length scale. Since boudinage and folding are considered to express the same general instability, both structures should arise independently of the sign of the loading conditions and for identical material parameters. To this end, the link between material to energy instabilities is approached by means of bifurcation analyses of the field equations and finite element simulations of the coupled system of equations. Boudinage and folding structures develop at the same critical energy threshold, where dissipative work by temperature-sensitive creep overcomes the diffusive capacity of the layer. This finding provides basis for a unified theory for strain localization in layered ductile materials. The numerical simulations are compared to natural pinch-and-swell microstructures, tracing the adaption of grain sizes, textures and creep mechanisms in calcite veins. The switch from dislocation to diffusion creep relates to strain-rate weakening, which is induced by dissipated heat from grain size reduction, and marks the onset of continuous necking. The time-dependent sequence uncovers multiple steady states at different time intervals. Microstructurally and mechanically stable conditions are finally expressed in the pinch-and-swell end members. The major outcome of this study is that boudinage and folding can be described as the same coupled energy-mechanical bifurcation, or as one critical energy attractor. This finding allows the derivation of critical deformation conditions and fundamental material parameters directly from localized structures in the field.

Effects of elastic indenter deformation on spherical instrumented indentation tests: the reduced elastic modulus

Although the Hertz theory is not applicable in the analysis of the indentation of elastic-plastic materials, it is common practice to incorporate the concept of indenter/specimen combined modulus to consider indenter deformation. The appropriateness was assessed of the use of reduced modulus to incorporate the effect of indenter deformation in the analysis of the indentation with spherical indenters. The analysis based on finite element simulations considered four values of the ratio of the indented material elastic modulus to that of the diamond indenter, E/E(i) (0, 0.04, 0.19, 0.39), four values of the ratio of the elastic reduced modulus to the initial yield strength, E(r)/Y (0, 10, 20, 100), and two values of the ratio of the indenter radius to maximum total displacement, R/delta(max) (3, 10). Indenter deformation effects are better accounted for by the reduced modulus if the indented material behaves entirely elastically. In this case, identical load-displacement (P - delta) curves are obtained with rigid and elastic spherical indenters for the same elastic reduced modulus. Changes in the ratio E/E(i), from 0 to 0.39, resulted in variations lower than 5% for the load dimensionless functions, lower than 3% in the contact area, A(c), and lower than 5% in the ratio H/E(r). However, deformations of the elastic indenter made the actual radius of contact change, even in the indentation of elastic materials. Even though the load dimensionless functions showed only a little increase with the ratio E/E(i), the hardening coefficient and the yield strength could be slightly overestimated when algorithms based on rigid indenters are used. For the unloading curves, the ratio delta(e)/delta(max), where delta(e) is the point corresponding to zero load of a straight line with slope S from the point (P(max), delta(max)), varied less than 5% with the ratio E/E(i). Similarly, the relationship between reduced modulus and the unloading indentation curve, expressed by Sneddon`s equation, did not reveal the necessity of correction with the ratio E/E(i). The most affected parameter in the indentation curve, as a consequence of the indentation deformation, was the ratio between the residual indentation depth after complete unloading and the maximum indenter displacement, delta(r)/delta(max) (up to 26%), but this variation did not significantly decrease the capability to estimate hardness and elastic modulus based on the ratio of the residual indentation depth to maximum indentation depth, h(r)/h(max). In general, the results confirm the convenience of the use of the reduced modulus in the spherical instrumented indentation tests.

Ultimate and failure strengths of normal muds and serpentinite muds of ODP Leg 125 samples

Serpentinite seamounts in the Mariana forearc have been explained as diapirs rising from the Benioff zone. This hypothesis predicts that the serpentinites should have low strengths as well as low densities relative to the surrounding rocks. Drilling during Leg 125 showed that the materials forming Conical Seamount in the Mariana forearc and Torishima Forearc Seamount in the Izu-Bonin forearc are water-charged serpentinite muds of density <2 g/cm**3. Wykeham-Farrance torsion-vane tests showed that they are plastic solids with a rheology that bears many similarities to the idealized Cam clay soil model and is well described by critical-state soil mechanics. The serpentinite muds have ultimate strengths of 1.3 to 273.7 kPa and yield strengths of approximately 1.0 to 50 kPa. These muds thus are orders of magnitude weaker than salt and are, in fact, comparable in density and strength to common deep-sea clay muds. Such weak and low-density materials easily become diapiric. Serpentinite muds from the summit of Conical Seamount are weaker and more ductile than those on its flanks or on Torishima Forearc Seamount. Moreover, the summit muds do not contain the pronounced pure- and simple-shear fabrics that characterize those on the seamount flanks. The seamounts are morphologically similar to shield volcanoes, and anastomosing serpentinite debris flows descending from their summits are similar in map view to pahoehoe flows. These morphologic features, together with the physical properties of the muds and their similarities to other oceanic muds and the geochemistry of the entrained waters, suggest that many forearc serpentinite seamounts are gigantic (10-20 km wide, 1.5-2.0 km high) mud volcanoes that formed by the eruption of highly liquid serpentinite muds. Torishima Forearc Seamount, which is blanketed by more ìnormalî pelagic/volcaniclastic sediment, has probably been inactive since the Miocene. Conical Seamount, which seems to consist entirely of serpentinite mud and is venting fresh water of unusual chemistry from its summit, is presently active. Muds from the flanks of Conical Seamount are stronger and more brittle than those from the summit site, and muds from Torishima Forearc Seamount are stronger yet; this suggests that the serpentinite debris flows are compacted and dewatered as they mature. The shear fabrics probably result from downslope creep and flow, but may also be inherited.

GENERALIZED FINITE ELEMENT METHOD FOR NONLINEAR THREE-DIMENSIONAL ANALYSIS OF SOLIDS

The Generalized Finite Element Method (GFEM) is employed in this paper for the numerical analysis of three-dimensional solids tinder nonlinear behavior. A brief summary of the GFEM as well as a description of the formulation of the hexahedral element based oil the proposed enrichment strategy are initially presented. Next, in order to introduce the nonlinear analysis of solids, two constitutive models are briefly reviewed: Lemaitre`s model, in which damage and plasticity are coupled, and Mazars`s damage model suitable for concrete tinder increased loading. Both models are employed in the framework of a nonlocal approach to ensure solution objectivity. In the numerical analyses carried out, a selective enrichment of approximation at regions of concern in the domain (mainly those with high strain and damage gradients) is exploited. Such a possibility makes the three-dimensional analysis less expensive and practicable since re-meshing resources, characteristic of h-adaptivity, can be minimized. Moreover, a combination of three-dimensional analysis and the selective enrichment presents a valuable good tool for a better description of both damage and plastic strain scatterings.

Plastic-damage smeared crack model to simulate the behaviour of structures made by cement based materials

This work proposes a constitutive model to simulate nonlinear behaviour of cement based materials subjected to different loading paths. The model incorporates a multidirectional fixed smeared crack approach to simulate crack initiation and propagation, whereas the inelastic behaviour of material between cracks is treated by a numerical strategy that combines plasticity and damage theories. For capturing more realistically the shear stress transfer between the crack surfaces, a softening diagram is assumed for modelling the crack shear stress versus crack shear strain. The plastic damage model is based on the yield function, flow rule and evolution law for hardening variable, and includes an explicit isotropic damage law to simulate the stiffness degradation and the softening behaviour of cement based materials in compression. This model was implemented into the FEMIX computer program, and experimental tests at material scale were simulated to appraise the predictive performance of this constitutive model. The applicability of the model for simulating the behaviour of reinforced concrete shear wall panels submitted to biaxial loading conditions, and RC beams failing in shear is investigated.

An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.

Hygro-elasto-plastic behavior of planar orthotropic material

The mechanical and hygroscopic properties of paper and board are factors affecting the whole lifecycle of a product, including paper/board quality, production, converting, and material and energy savings. The progress of shrinkage profiles, loose edges of web, baggy web causing wrinkling and misregistration in printing are examples of factors affecting runnability and end product quality in the drying section and converting processes, where paper or board is treated as a moving web. The structural properties and internal stresses or plastic strain differences built up during production also cause the end-product defects related to distortion of the shape of the product such as sheet or box. The objective of this work was to construct a model capable of capturing the characteristic behavior of hygroscopic orthotropic material under moisture change, during different external in-plane stretch or stress conditions. Two independent experimental models were constructed: the elasto-plastic material model and the hygroexpansivity-shrinkage model. Both describe the structural properties of the sheet with a fiber orientation probability distribution, and both are functions of the dry solids content and fiber orientation anisotropy index. The anisotropy index, introduced in this work, simplifies the procedure of determining the constitutive parameters of the material model and the hygroexpansion coefficients in different in-plane directions of the orthotropic sheet. The mathematically consistent elasto-plastic material model and the dry solids content dependent hygroexpansivity have been constructed over the entire range from wet to dry. The presented elastoplastic and hygroexpansivity-shrinkage models can be used in an analytical approach to estimate the plastic strain and shrinkage in simple one-dimensional cases. For studies of the combined and more complicated effects of hygro-elasto-plastic behavior, both models were implemented in a finite element program for a numerical solution. The finite element approach also offered possibilities for studying different structural variations of orthotropic planar material, as well as local buckling behavior and internal stress situations of the sheet or web generated by local strain differences. A comparison of the simulation examples presented in this work to results published earlier confirms that the hygro-elasto-plastic model provides at least qualitatively reasonable estimates. The application potential of the hygro-elasto-plastic model is versatile, including several phenomena and defects appearing in the drying, converting and end-use conditions of the paper or board webs and products, or in other corresponding complex planar materials.

The ecology of percolating filters containing a plastic filter medium in relation to their efficiency in The treatment of domestic sewage

The suitability of a new plastic supporting medium for biofiltration was tested over a three year period. Tests were carried out on the stability, surface properties, mechanical strength, and dimensions of the medium. There was no evidence to suggest that the medium was deficient in any of these respects. The specific surface (320m2m-3) and the voidage (94%) of the new medium are unlike any other used in bio-filtration and a pilot plant containing two filters was built to observe its effects on ecology and performance. Performance was estimated by chemical analysis and ecology studied by film examination and fauna counts. A system of removable sampling baskets was designed to enable samples to be obtained from two intermediate depths of filter. One of the major operating problems of percolating filters is excessive accumulation of film. The amount of film is influenced by hydraulic and organic load and each filter was run at a different loading. One was operated at 1.2m3m-3day-1 (DOD load 0.24kgm-3day-1) judged at the time to be the lowest filtration rate to offer advantages over conventional media. The other filter was operated at more than twice this loading (2.4m3m-3day-lBOD load 0.55kgm-3day-1) giving a roughly 2.5x and 6x the conventional loadings recommended for a Royal Commission effluent. The amount of film in each filter was normally low (0.05-3kgm(3 as volatile solids) and did not affect efficiency. The evidence collected during the study indicated that the ecology of the filters was normal when compared with the data obtained from the literature relating to filters with mineral media. There were indications that full ecological stability was yet to be reached and this was affecting the efficiency of the filters. The lower rate filter produced an average 87% BOD removal giving a consistent Royal Commission effluent during the summer months. The higher rate filter produced a mean 83% BOD removal but at no stage a consistent Royal Commission effluent. From the data on ecology and performance the filters resembled conventional filters rather than high rate filters.

Influência da embalagem de polietileno de baixa densidade e da temperatura na conservação do repolho minimamente processado

This research studied the effect of low density polyethylene packaging and storage temperature on the preservation of fresh-cut (minimally processed) cabbage. The cabbages, previously cooled to a temperature of 10 ºC, were selected, washed, cut in four parts (with the central stalk removed), sanitized, cut in strips, rinsed, put in the centrifuge, weighed and stored in plastic packaging of low density polyethylene (70 µm), and then stored in cold chambers at temperatures of 1 and 10 ºC for 20 days. The following aspects were evaluated: carbon dioxide, oxygen and ethylene in the internal atmosphere of the package as well as, pH, titratable acidity, total soluble solids, vitamin C, loss of fresh mass and the total soluble solids/acidity in the fresh-cut cabbage ratio. The experimental design was entirely casual, with three repetitions. The analysis parameters, except for the vitamin C, loss of fresh mass and ethylene, presented significant variation between the temperatures and days of storage. The cabbage stored at a temperature of 1 ºC presented a shelf life of around 15 days, significantly higher than that stored at 10 ºC. At this temperature, on the 8th day of storage, the product was completely decayed, unfit for commercialization or consumption.

Analysis of the tip roundness effects on the micro- and macroindentation response of elastic-plastic materials

In this work, the effects of indenter tip roundness oil the load-depth indentation curves were analyzed using finite element modeling. The tip roundness level was Studied based on the ratio between tip radius and maximum penetration depth (R/h(max)), which varied from 0.02 to 1. The proportional Curvature constant (C), the exponent of depth during loading (alpha), the initial unloading slope (S), the correction factor (beta), the level of piling-up or sinking-in (h(c)/h(max)), and the ratio h(max)/h(f) are shown to be strongly influenced by the ratio R/h(max). The hardness (H) was found to be independent of R/h(max) in the range studied. The Oliver and Pharr method was successful in following the variation of h(c)/h(max) with the ratio R/h(max) through the variation of S with the ratio R/h(max). However, this work confirmed the differences between the hardness values calculated using the Oliver-Pharr method and those obtained directly from finite element calculations; differences which derive from the error in area calculation that Occurs when given combinations of indented material properties are present. The ratio of plastic work to total work (W(p)/W(t)) was found to be independent of the ratio R/h(max), which demonstrates that the methods for the Calculation of mechanical properties based on the *indentation energy are potentially not Susceptible to errors caused by tip roundness.

Analysis of the effects of conical indentation variables on the indentation response of elastic-plastic materials through factorial design of experiment

In this work, the effects of conical indentation variables on the load-depth indentation curves were analyzed using finite element modeling and dimensional analysis. A factorial design 2(6) was used with the aim of quantifying the effects of the mechanical properties of the indented material and of the indenter geometry. Analysis was based on the input variables Y/E, R/h(max), n, theta, E, and h(max). The dimensional variables E and h(max) were used such that each value of dimensionless Y/E was obtained with two different values of E and each value of dimensionless R/h(max) was obtained with two different h(max) values. A set of dimensionless functions was defined to analyze the effect of the input variables: Pi(1) = P(1)/Eh(2), Pi(2) = h(c)/h, Pi(3) = H/Y, Pi(4) = S/Eh(max), Pi(6) = h(max)/h(f) and Pi(7) = W(P)/W(T). These six functions were found to depend only on the dimensionless variables studied (Y/E, R/h(max), n, theta). Another dimension less function, Pi(5) = beta, was not well defined for most of the dimensionless variables and the only variable that provided a significant effect on beta was theta. However, beta showed a strong dependence on the fraction of the data selected to fit the unloading curve, which means that beta is especially Susceptible to the error in the Calculation of the initial unloading slope.

Evolution of magnetic properties and crystallographic texture in electrical steel with large plastic deformation

Deformation leads to a hardening of steel due to an increase in the density of dislocations and a reduction in their mobility, giving rise to a state of elevated residual stresses in the crystal lattice. In the microstructure, one observes an increase in the contribution of crystalline orientations which are unfavorable to the magnetization, as seen, for example, by a decrease in B(50), the magnetic flux density at a field of 50 A/cm. The present study was carried out with longitudinal strips of fully processed non-oriented (NO) electrical steel, with deformations up to 70% resulting from cold rolling in the longitudinal direction. With increasing plastic deformation, the value of B(50) gradually decreases until it reaches a minimum value, where it remains even for larger deformations. On the other hand, the coercive field H(c) continually increases. Magnetometry results and electron backscatter diffraction results are compared and discussed. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3560895]

Systematic investigation of a family of gradient-dependent functionals for solids

Eleven density functionals are compared with regard to their performance for the lattice constants of solids. We consider standard functionals, such as the local-density approximation and the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), as well as variations of PBE GGA, such as PBEsol and similar functionals, PBE-type functionals employing a tighter Lieb-Oxford bound, and combinations thereof. On a test set of 60 solids, we perform a system-by-system analysis for selected functionals and a full statistical analysis for all of them. The impact of restoring the gradient expansion and of tightening the Lieb-Oxford bound is discussed, and confronted with previous results obtained from other codes, functionals or test sets. No functional is uniformly good for all investigated systems, but surprisingly, and pleasingly, the simplest possible modifications to PBE turn out to have the most beneficial effect on its performance. The atomization energy of molecules was also considered and on a testing set of six molecules, we found that the PBE functional is clearly the best, the others leading to strong overbinding.

Gradient-dependent density functionals of the Perdew-Burke-Ernzerhof type for atoms, molecules, and solids

One of the standard generalized-gradient approximations (GGAs) in use in modern electronic-structure theory [Perdew-Burke-Ernzerhof (PBE) GGA] and a recently proposed modification designed specifically for solids (PBEsol) are identified as particular members of a family of functionals taking their parameters from different properties of homogeneous or inhomogeneous electron liquids. Three further members of this family are constructed and tested, together with the original PBE and PBEsol, for atoms, molecules, and solids. We find that PBE, in spite of its popularity in solid-state physics and quantum chemistry, is not always the best performing member of the family and that PBEsol, in spite of having been constructed specifically for solids, is not the best for solids. The performance of GGAs for finite systems is found to sensitively depend on the choice of constraints stemming from infinite systems. Guidelines both for users and for developers of density functionals emerge from this work.

High Solids Content, Soap-Free, Film-Forming Latexes Stabilized by Laponite Clay Platelets

High solids content film-forming poly[styrene-co-(n-butyl acrylate)] [poly(Sty-co-BuA)] latexes armored with Laponite clay platelets have been synthesized by soap-free emulsion copolymerization of styrene and n-butyl acrylate. The polymerizations were performed in batch in the presence of Laponite and a methyl ether acrylateterminated poly(ethylene glycol) macromonomer in order to promote polymer/clay association. The overall polymerization kinetics showed a pronounced effect of clay on nucleation and stabilization of the latex particles. Cryo-transmission electron microscopy observation confirmed the armored morphology and indicated that the majority of Laponite platelets were located at the particle surface. The resulting nanostructured films displayed enhanced mechanical properties.