Fourier-transform rheology on anionically synthesised polymer melts and solutions of various topology

Aim of this thesis was to further extend the applicability of the Fourier-transform (FT) rheology technique especially for non-linear mechanical characterisation of polymeric materials on the one hand and to investigated the influence of the degree of branching on the linear and non-linear relaxation behaviour of polymeric materials on the other hand. The latter was achieved by employing in particular FT-rheology and other rheological techniques to variously branched polymer melts and solutions. For these purposes, narrowly distributed linear and star-shaped polystyrene and polybutadiene homo-polymers with varying molecular weights were anionically synthesised using both high-vacuum and inert atmosphere techniques. Furthermore, differently entangled solutions of linear and star-shaped polystyrenes in di-sec-octyl phthalate (DOP) were prepared. The several linear polystyrene solutions were measured under large amplitude oscillatory shear (LAOS) conditions and the non-linear torque response was analysed in the Fourier space. Experimental results were compared with numerical predictions performed by Dr. B. Debbaut using a multi-mode differential viscoelastic fluid model obeying the Giesekus constitutive equation. Apart from the analysis of the relative intensities of the harmonics, a detailed examination of the phase information content was developed. Further on, FT-rheology allowed to distinguish polystyrene melts and solutions due to their different topologies where other rheological measurements failed. Significant differences occurred under LAOS conditions as particularly reflected in the phase difference of the third harmonic, Ħ3, which could be related to shear thinning and shear thickening behaviour.

An anisotropic elastic-viscoplastic damage model for bone tissue

A new anisotropic elastic-viscoplastic damage constitutive model for bone is proposed using an eccentric elliptical yield criterion and nonlinear isotropic hardening. A micromechanics-based multiscale homogenization scheme proposed by Reisinger et al. is used to obtain the effective elastic properties of lamellar bone. The dissipative process in bone is modeled as viscoplastic deformation coupled to damage. The model is based on an orthotropic ecuntric elliptical criterion in stress space. In order to simplify material identification, an eccentric elliptical isotropic yield surface was defined in strain space, which is transformed to a stress-based criterion by means of the damaged compliance tensor. Viscoplasticity is implemented by means of the continuous Perzyna formulation. Damage is modeled by a scalar function of the accumulated plastic strain D(κ) , reducing all element s of the stiffness matrix. A polynomial flow rule is proposed in order to capture the rate-dependent post-yield behavior of lamellar bone. A numerical algorithm to perform the back projection on the rate-dependent yield surface has been developed and implemented in the commercial finite element solver Abaqus/Standard as a user subroutine UMAT. A consistent tangent operator has been derived and implemented in order to ensure quadratic convergence. Correct implementation of the algorithm, convergence, and accuracy of the tangent operator was tested by means of strain- and stress-based single element tests. A finite element simulation of nano- indentation in lamellar bone was finally performed in order to show the abilities of the newly developed constitutive model.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

Caputo Derivatives in Viscoelasticity: A Non-Linear Finite-Deformation Theory for Tissue

Mathematics Subject Classification: 26A33, 74B20, 74D10, 74L15

Polymer Process Partitioning: Extruding Plastic Pipe

When plastic pipe is solidified, it proceeds through a long cooling chamber. Inside this chamber, inside the hollow extrudate, the plastic is molten, and this inner surface solidifies last. Sag, the flow due to the self-weight of the molten plastic, then happens in this cooling chamber, and sometimes, thickened regions (called knuckles) arise in the lower quadrants, especially of large diameter thickwalled pipes. To compensate for sag, engineers normally shift the die centerpiece downward. This thesis focuses on the consequences of this decentering. Specifically, when the molten polymer is viscoelastic, as is normally the case, a downward lateral force is exerted on the mandrel. Die eccentricity also affects the downstream axial force on the mandrel. These forces govern how rigidly the mandrel must be attached (normally, on a spider die). We attack this flow problem in eccentric cylindrical coordinates, using the Oldroyd 8-constant constitutive model framework. Specifically, we revise the method of Jones (1964), called polymer process partitioning. We estimate both axial and lateral forces. We develop a corresponding map to help plastics engineers predict the extrudate shape, including extrudate knuckles. From the mass balance over the postdie region, we then predict the shape of the extrudate entering the cooling chamber. We further include expressions for the stresses in the extruded polymer melt. We include detailed dimensional worked examples to show process engineers how to use our results to design pipe dies, and especially to suppress extrudate knuckling.

Semi-implicit finite strain constitutive integration and mixed strain/stress control based on intermediate configurations

A new semi-implicit stress integration algorithm for finite strain plasticity (compatible with hyperelas- ticity) is introduced. Its most distinctive feature is the use of different parameterizations of equilibrium and reference configurations. Rotation terms (nonlinear trigonometric functions) are integrated explicitly and correspond to a change in the reference configuration. In contrast, relative Green–Lagrange strains (which are quadratic in terms of displacements) represent the equilibrium configuration implicitly. In addition, the adequacy of several objective stress rates in the semi-implicit context is studied. We para- metrize both reference and equilibrium configurations, in contrast with the so-called objective stress integration algorithms which use coinciding configurations. A single constitutive framework provides quantities needed by common discretization schemes. This is computationally convenient and robust, as all elements only need to provide pre-established quantities irrespectively of the constitutive model. In this work, mixed strain/stress control is used, as well as our smoothing algorithm for the complemen- tarity condition. Exceptional time-step robustness is achieved in elasto-plastic problems: often fewer than one-tenth of the typical number of time increments can be used with a quantifiable effect in accuracy. The proposed algorithm is general: all hyperelastic models and all classical elasto-plastic models can be employed. Plane-stress, Shell and 3D examples are used to illustrate the new algorithm. Both isotropic and anisotropic behavior is presented in elasto-plastic and hyperelastic examples.

Evaluation of the geotechnical properties of MSW in two Brazilian landfills

The characteristics of municipal solid waste (MSW) play a key role in many aspects of waste disposal facilities and landfills. Because most of a landfill is made up of MSW, the overall stability of the landfill slopes are governed by the strength parameters and physical properties of the MSW. These parameters are also important in interactions involving the waste body and the landfill structures: cover liner, leachate and gas collection systems. On the other hand, the composition of the waste, which affects the geotechnical behavior of the MSW, is dependent on a variety of factors such as climate, disposal technology, the culture and habits of the local community. It is therefore essential that the design and stability evaluations of landfills in each region be performed based on the local conditions and the geotechnical characteristic of the MSW. The Bandeirantes Landfill, BL, in Sao Paulo and the Metropolitan Center Landfill, MCL, in Salvador, are among the biggest landfills in Brazil. These two disposal facilities have been used for the development of research involving waste mechanics in recent years. Considerable work has been made in the laboratory and in the field to evaluate parameters such as water and organic contents, composition, permeability, and shear strength. This paper shows and analyzes the results of tests performed on these two landfills. The authors believe that these results could be a good reference for certain aspects and geotechnical properties of MSW materials in countries with similar conditions. (C) 2010 Elsevier Ltd. All rights reserved.

Texture alignment in simple shear

We illustrate the flow behaviour of fluids with isotropic and anisotropic microstructure (internal length, layering with bending stiffness) by means of numerical simulations of silo discharge and flow alignment in simple shear. The Cosserat theory is used to provide an internal length in the constitutive model through bending stiffness to describe isotropic microstructure and this theory is coupled to a director theory to add specific orientation of grains to describe anisotropic microstructure. The numerical solution is based on an implicit form of the Material Point Method developed by Moresi et al. [1].

Finite element studies of the mechanical behaviour of the diaphragm in normal and pathological cases

The diaphragm is a muscular membrane separating the abdominal and thoracic cavities, and its motion is directly linked to respiration. In this study, using data from a 59-year-old female cadaver obtained from the Visible Human Project, the diaphragm is reconstructed and, from the corresponding solid object, a shell finite element mesh is generated and used in several analyses performed with the ABAQUS 6.7 software. These analyses consider the direction of the muscle fibres and the incompressibility of the tissue. The constitutive model for the isotropic strain energy as well as the passive and active strain energy stored in the fibres is adapted from Humphrey's model for cardiac muscles. Furthermore, numerical results for the diaphragmatic floor under pressure and active contraction in normal and pathological cases are presented.

Predictive analytical and numerical modeling for orthogonal cutting

In this thesis, a predictive analytical and numerical modeling approach for the orthogonal cutting process is proposed to calculate temperature distributions and subsequently, forces and stress distributions. The models proposed include a constitutive model for the material being cut based on the work of Weber, a model for the shear plane based on Merchants model, a model describing the contribution of friction based on Zorev’s approach, a model for the effect of wear on the tool based on the work of Waldorf, and a thermal model based on the works of Komanduri and Hou, with a fraction heat partition for a non-uniform distribution of the heat in the interfaces, but extended to encompass a set of contributions to the global temperature rise of chip, tool and work piece. The models proposed in this work, try to avoid from experimental based values or expressions, and simplifying assumptions or suppositions, as much as possible. On a thermo-physical point of view, the results were affected not only by the mechanical or cutting parameters chosen, but also by their coupling effects, instead of the simplifying way of modeling which is to contemplate only the direct effect of the variation of a parameter. The implementation of these models was performed using the MATLAB environment. Since it was possible to find in the literature all the parameters for AISI 1045 and AISI O2, these materials were used to run the simulations in order to avoid arbitrary assumption.

Strengthening to seismic action of reinforced columns with rectangular cross-section

The first part of this research work regards the assessment of the mathematical modelling of reinforced concrete columns confined with carbon fibre (CFRP) sheets under axial loading. The purpose was to evaluate existing analytical models, contribute to possible improvements and choose the best model(s) to be part of a new model for the prediction of the behaviour of confined columns under bending and compression. For circular columns, a wide group of authors have proposed several models specific for FRP-confined concrete. The analysis of some of the existing models was carried out by comparing these with several tested columns. Although several models predict fairly the peak load only few can properly estimate the load-strain and dilation behaviour of the columns. Square columns confined with CFRP show a more complex interpretation of their behaviour. Accordingly, the analysis of two experimental programs was carried out to propose new modelling equations for the whole behaviour of columns. The modelling results show that the analytical curves are in general agreement with the presented experimental curves for a wide range of dimensions. An analysis similar to the one done for circular columns was this turn carried out for square columns. Few models can fairly estimate the whole behaviour of the columns and with less accuracy at all levels when compared with circular columns. The second part of this study includes seven experimental tests carried out on reinforced concrete rectangular columns with rounded corners, different damage condition and with confinement and longitudinal strengthening systems. It was concluded that the use of CFRP confinement is viable and of effective performance enhancement alone and combined with other techniques, maintaining a good ductile behaviour for established threshold displacements. As regards the use of external longitudinal strengthening combined with CFRP confinement, this system is effective for the performance enhancement and viable in terms of execution. The load capacity was increased significantly, preserving also in this case a good ductile behaviour for threshold displacements. As to the numerical nonlinear modelling of the tested columns, the results show a variation of the peak load of 1% to 10% compared with tests results. The good results are partly due to the inclusion of the concrete constitutive model by Mander et al. modified by Faustino, Chastre & Paula taking into account the confinement effect. Despite the reasonable approximation to tests results, the modelling results showed higher unloading, which leads to an overestimate dissipated energy and residualdisplacement.

Exploring the use of HPFRC and GFRP grids for the production of manhole covers

Corrosion of the steel reinforced concrete elements is one of the common pathologies that limits the long-term performance of urban infrastructures. This problem causes the loss of structural serviceability by decreasing the concrete-steel bond strength and reducing the cross section of the reinforcements. The present study introduces a new system for developing free-corrosion resistance prefabricated manhole covers for applications in the aggressive environments, i.e. wastewater collector systems, sewer systems, stormwater systems, etc. Fibre reinforced cement composites were applied in this system in order to suppress the corrodible steel mesh and maintain the structural ductility as well. Application of fibre reinforced polymer (FRP) system is adopted as the additional solution for increasing the load carrying capacity of these elements without concerns about corrosion. The effectiveness of the applied strategy in developing the manhole covers in terms of load carrying capacity and failure mode is evaluated in this research. Furthermore, this paper discusses a FEM-based simulation, aiming to address the possibility of calibrating the constitutive model parameters related to fracture modes I and II.

High performance fiber reinforced concrete for the shear reinforcement: experimental and numerical research

High performance fiber reinforced concrete (HPFRC) is developing rapidly to a modern structural material with unique rheological and mechanical characteristics. Despite applying several methodologies to achieve self15 compacting requirements, some doubts still remain regarding the most convenient strategy for developing a HPFRC. In the present study, an innovative mix design method is proposed for the development of high17 performance concrete reinforced with a relatively high dosage of steel fibers. The material properties of the developed concrete are assessed, and the concrete structural behavior is characterized under compressive, flexural and shear loading. This study better clarifies the significant contribution of fibers for shear resistance of concrete elements. This paper further discusses a FEM-based simulation, aiming to address the possibility of calibrating the constitutive model parameters related to fracture modes I and II.

Retroanálise de parâmetros geomecânicos baseada em algoritmos de otimização

Dissertação de mestrado integrado em Engenharia Civil

Zur Druck-Querzug-Festigkeit von Stahlbeton und stahlfaserverstärktem Stahlbeton in scheibenförmigen Bauteilen

Am Fachgebiet Massivbau (Institut für Konstruktiven Ingenieurbau – IKI) des Fachbereichs Bauingenieurwesen der Universität Kassel wurden Bauteilversuche an zweiaxial auf Druck-Zug belasteten, faserfreien und faserverstärkten Stahlbetonscheiben durchgeführt. Dabei wurden die Auswirkungen der Querzugbeanspruchung und der Rissbildung auf die Druckfestigkeit, auf die Stauchung bei Erreichen der Höchstlast sowie auf die Drucksteifigkeit des stabstahl- und faserbewehrten Betons an insgesamt 56 faserfreien und faserverstärkten Beton- und Stahlbetonscheiben untersucht. Auf der Grundlage der experimentell erhaltenen Ergebnisse wird ein Vorschlag zur Abminderung der Druckfestigkeit des gerissenen faserfreien und faserbewehrten Stahlbetons in Abhängigkeit der aufgebrachten Zugdehnung formuliert. Die Ergebnisse werden den in DIN 1045-1 [D4], Eurocode 2 [E3, E4], CEB-FIP Model Code 1990 [C1] und ACI Standard 318-05 [A1] angegebenen Bemessungsregeln für die Druckstrebenfestigkeit des gerissenen Stahlbetons gegenübergestellt und mit den Untersuchungen anderer Wissenschaftler verglichen. Die bekannten Widersprüche zwischen den Versuchsergebnissen, den vorgeschlagenen Modellen und den Regelwerken aus U.S.A., Kanada und Europa können dabei weitgehend aufgeklärt werden. Für nichtlineare Verfahren der Schnittgrößenermittlung und für Verformungsberechnungen wird ein Materialmodell des gerissenen faserfreien und faserbewehrten Stahlbetons abgeleitet. Hierzu wird die für einaxiale Beanspruchungszustände gültige Spannungs-Dehnungs-Linie nach Bild 22 der DIN 1045-1 auf den Fall der zweiaxialen Druck-Zug-Beanspruchung erweitert.