Estabilidade e controle dinâmico de roll waves

Pós-graduação em Engenharia Elétrica - FEIS

Influência do teor de cinza volante na reologia de argamassa de assentamento: fabricado com reciclagem de resíduos

As argamassas podem ser consideradas como um conjunto de partículas sólidas em suspensão (agregado miúdo), aglomerantes e água (pasta de cimento). Sendo assim, seu comportamento reológico durante a cura pode ser influenciado por fatores como a granulometria, composição e teor de água incorporado. Este trabalho estudou as características reológicas de seis composições de argamassas, em função do teor de cinza volante, como material de atividade pozolânica, em substituição parcial do cimento. A argamassa estudada utilizou como base agregado o resíduo da construção civil (RCC), devidamente caracterizado por difratometria e fluorescência de raios-X. Formularam-se composições de argamassas para os ensaios reológicos com adição de resíduo de construção civil (RCC) nas proporções de 95, 90, 85, 80, 75 e 70%; cinzas volantes (CV) 0, 5, 10, 15, 20 e 25% e 5% de Cimento Portland Comum (CP II- E 32), sendo que em cada amostra foi incorporado 35% de água. Utilizou-se o viscosímetro modelo VT 550, com sensor tipo cilindros coaxiais SV₁ – Haake a temperatura de 28ºC. Para avaliar o tempo de cura do material programou-se uma taxa de cisalhamento constante 53,4 s^-1 em 1h e 45min, sendo avaliado o torque em intervalos de 15min. Para a elaboração das curvas de fluxo e curvas de histerese, utilizou-se a taxa de cisalhamento entre 0 e 600 s^-1 no intervalo de tempo de 120s. Os resultados mostraram que a incorporação de cinza volante foi o principal fator para aumento do torque no tempo de cura das argamassas e mudança do comportamento de reópetico para tixotrópico, nos traços analisados. Os dados experimentais das formulações propostas ajustaram-se ao modelo reológico de Herschel-Bulkley.

Control of instabilities in non-Newtonian free surface fluid flows

Free surface flows in inclined channels can develop periodic instabilities that are propagated downstream as shock waves with well-defined wavelengths and amplitudes. Such disturbances are called roll waves and are common in channels, torrential lava, landslides, and avalanches. The prediction and detection of such waves over certain types of structures and environments are useful for the prevention of natural risks. In this work, a mathematical model is established using a theoretical approach based on Cauchy's equations with the Herschel-Bulkley rheological model inserted into the viscous part of the stress tensor. This arrangement can adequately represent the behavior of muddy fluids, such as water-clay mixture. Then, taking into account the shallow water and the Rankine-Hugoniot's (shock wave) conditions, the equation of the roll wave and its properties, profile, and propagation velocity are determined. A linear stability analysis is performed with an emphasis on determining the condition that allows the generation of such instabilities, which depends on the minimum Froude number. A sensitivity analysis on the numerical parameters is performed, and numerical results including the influence of the Froude number, the index flow and dimensionless yield stress on the amplitude, the wavelength of roll waves and the propagation velocity of roll waves are shown. We show that our numerical results were in agreement with Coussot's experimental results (1994).

Determinação de parâmetros reológicos de argamassas e concreto através de técnicas alternativas

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Implementação de aparato experimental para medição de instabilidades tipo Roll Waves em fluidos não- newtonianos

Pós-graduação em Engenharia Mecânica - FEIS

Estudo da influência da cinza de casca de arroz e da cal nas propriedades do solo-cimento auto-adensável

Pós-graduação em Engenharia Civil - FEIS

In situ micropillar compression reveals superior strength and ductility but an absence of damage in lamellar bone

Ageing 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 hierarchical organization of bone. Here, we investigate the mechanical response under monotonic and cyclic compression of both single osteonal lamellae and macroscopic samples containing numerous osteons. Micropillar compression tests in a scanning electron microscope, microindentation and macroscopic compression tests were performed on dry ovine bone to identify the elastic modulus, yield stress, plastic deformation, damage accumulation and failure mechanisms. We found that isolated lamellae exhibit a plastic behaviour, with higher yield stress and ductility but no damage. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behaviour of bone at the microscale to a quasi-brittle response driven by the growth of cracks along interfaces or in the vicinity of pores at the macroscale.

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.

Modelling of landslides propagation with SPH : effects of rheology and pore water pressure

El estudio desarrollado en este trabajo de tesis se centra en la modelización numérica de la fase de propagación de los deslizamientos rápidos de ladera a través del método sin malla Smoothed Particle Hydrodynamics (SPH). Este método tiene la gran ventaja de permitir el análisis de problemas de grandes deformaciones evitando operaciones costosas de remallado como en el caso de métodos numéricos con mallas tal como el método de los Elementos Finitos. En esta tesis, particular atención viene dada al rol que la reología y la presión de poros desempeñan durante estos eventos. El modelo matemático utilizado se basa en la formulación de Biot-Zienkiewicz v - pw, que representa el comportamiento, expresado en términos de velocidad del esqueleto sólido y presiones de poros, de la mezcla de partículas sólidas en un medio saturado. Las ecuaciones que gobiernan el problema son: • la ecuación de balance de masa de la fase del fluido intersticial, • la ecuación de balance de momento de la fase del fluido intersticial y de la mezcla, • la ecuación constitutiva y • una ecuación cinemática. Debido a sus propiedades geométricas, los deslizamientos de ladera se caracterizan por tener una profundidad muy pequeña frente a su longitud y a su anchura, y, consecuentemente, el modelo matemático mencionado anteriormente se puede simplificar integrando en profundidad las ecuaciones, pasando de un modelo 3D a 2D, el cual presenta una combinación excelente de precisión, sencillez y costes computacionales. El modelo propuesto en este trabajo se diferencia de los modelos integrados en profundidad existentes por incorporar un ulterior modelo capaz de proveer información sobre la presión del fluido intersticial a cada paso computacional de la propagación del deslizamiento. En una manera muy eficaz, la evolución de los perfiles de la presión de poros está numéricamente resuelta a través de un esquema explicito de Diferencias Finitas a cada nodo SPH. Este nuevo enfoque es capaz de tener en cuenta la variación de presión de poros debida a cambios de altura, de consolidación vertical o de cambios en las tensiones totales. Con respecto al comportamiento constitutivo, uno de los problemas principales al modelizar numéricamente deslizamientos rápidos de ladera está en la dificultad de simular con la misma ley constitutiva o reológica la transición de la fase de iniciación, donde el material se comporta como un sólido, a la fase de propagación donde el material se comporta como un fluido. En este trabajo de tesis, se propone un nuevo modelo reológico basado en el modelo viscoplástico de Perzyna, pensando a la viscoplasticidad como a la llave para poder simular tanto la fase de iniciación como la de propagación con el mismo modelo constitutivo. Con el fin de validar el modelo matemático y numérico se reproducen tanto ejemplos de referencia con solución analítica como experimentos de laboratorio. Finalmente, el modelo se aplica a casos reales, con especial atención al caso del deslizamiento de 1966 en Aberfan, mostrando como los resultados obtenidos simulan con éxito estos tipos de riesgos naturales. The study developed in this thesis focuses on the modelling of landslides propagation with the Smoothed Particle Hydrodynamics (SPH) meshless method which has the great advantage of allowing to deal with large deformation problems by avoiding expensive remeshing operations as happens for mesh methods such as, for example, the Finite Element Method. In this thesis, special attention is given to the role played by rheology and pore water pressure during these natural hazards. The mathematical framework used is based on the v - pw Biot-Zienkiewicz formulation, which represents the behaviour, formulated in terms of soil skeleton velocity and pore water pressure, of the mixture of solid particles and pore water in a saturated media. The governing equations are: • the mass balance equation for the pore water phase, • the momentum balance equation for the pore water phase and the mixture, • the constitutive equation and • a kinematic equation. Landslides, due to their shape and geometrical properties, have small depths in comparison with their length or width, therefore, the mathematical model aforementioned can then be simplified by depth integrating the equations, switching from a 3D to a 2D model, which presents an excellent combination of accuracy, computational costs and simplicity. The proposed model differs from previous depth integrated models by including a sub-model able to provide information on pore water pressure profiles at each computational step of the landslide's propagation. In an effective way, the evolution of the pore water pressure profiles is numerically solved through a set of 1D Finite Differences explicit scheme at each SPH node. This new approach is able to take into account the variation of the pore water pressure due to changes of height, vertical consolidation or changes of total stress. Concerning the constitutive behaviour, one of the main issues when modelling fast landslides is the difficulty to simulate with the same constitutive or rheological model the transition from the triggering phase, where the landslide behaves like a solid, to the propagation phase, where the landslide behaves in a fluid-like manner. In this work thesis, a new rheological model is proposed, based on the Perzyna viscoplastic model, thinking of viscoplasticity as the key to close the gap between the triggering and the propagation phase. In order to validate the mathematical model and the numerical approach, benchmarks and laboratory experiments are reproduced and compared to analytical solutions when possible. Finally, applications to real cases are studied, with particular attention paid to the Aberfan flowslide of 1966, showing how the mathematical model accurately and successfully simulate these kind of natural hazards.

Estudio de la fluencia de una calcarenita: la Piedra de San Julián (Alicante)

La fluencia de las rocas tiene gran importancia en la evaluación del comportamiento a largo plazo de elementos construidos con estos materiales. En este trabajo, se ha caracterizado física y mecánicamente una calcarenita porosa bien conocida localmente como Piedra de San Julián. Se han realizado ensayos de compresión uniaxial de 96 h. a carga constante. Se ha utilizado un modelo de fluencia bien conocido, el Código-modelo CEB-FIP 2010, usado para modelizar otro material pétreo (hormigón). Además, se ha propuesto un modelo reológico. El objetivo principal de este trabajo es investigar la posibilidad de aprovechar la gran experiencia acumulada en el estudio del hormigón, con el fin de obtener un enfoque para el comportamiento de la roca, para tiempos de prueba muy largos difíciles de implementar en laboratorio. Se propone una función de fluencia adaptada a la roca estudiada dependiente sólo de sus características elásticas y mecánicas.

Comparative Analysis of Viscoelastic Models Involving Fractional Derivatives of Different Orders

Mathematics Subject Classification: 74D05, 26A33

Estudo da incorporação de resíduos de perfuração para aplicação em cimentação de poços de petróleo

During the drilling of oil and natural gas are generated solid waste, liquid and gaseous. These solid fragments, which are known as cuttings, are carried to the surface through the drilling fluid. Furthermore, this fluid serves to cool the bit, keeping the internal pressure of the well, and others. This solid residue is very polluting, because it has incorporated beyond the drilling fluid, which has several chemical additives harmful to the environment, some heavy metals that are harmful to the environment, such as lead. To minimize the residue generated, are currently being studied numerous techniques to mitigate the problems that such waste can cause to the environment, like addition of cuttings in the composition of soil cement brick masonry construction, addition of cuttings on the clay matrix for the manufacture of solid masonry bricks and ceramic blocks and coprocessing of the cuttings in cement. So, the main objective of this work is the incorporation of cuttings drilling of oil wells, the cement slurry used in the cementing operation of the well. This cuttings used in this study, arising from the formation Pendências, was milled and separated in a sieve of 100 mesh. After grinding had a mean particle sike in order of 86 mm and crystal structure containing phases of quartz and calcite type, characteristic of the Portland cement. Were formulated and prepared slurries of cement with density 13 lb / gal, containing different concentrations of gravel, and realized characterization tests API SPEC 10A and RP 10B. Free water tests showed values lower than 5.9% and the rheological model that best described the behavior of the mixtures was the power. The results of compressive strength (10.3 MPa) and stability (Dr <0.5 lb / gal) had values within the set of operational procedures. Thus, the gravel from the drilling operation, may be used as binders in addition to Portland cement oil wells, in order to reuse this waste and reduce the cost of the cement paste.

Suppression of vibrations in a nonlinear half-car model using a magneto-rheological damper

The present work investigates the nonlinear response of a half-car model. The disturbances of the road are assumed to be sinusoidal. After constructing the bifurcation diagram, we use the 0-1 test to identify chaotic motions. The main objective of this study is to eliminate chaotic behavior of the chassis and reduce its vibrations. To accomplish this, a semi-active vehicle suspension control system, using magneto-rheological dampers, is proposed. The proposed semi-active control strategy consists of two nonlinear control laws: a feedforward control, and a feedback control. They are obtained by considering the SDRE (State Dependent Riccati Equation) control, where the control parameter is the voltage applied to the coils of the magneto-rheological dampers. Numerical results show that the proposed control method is effective in significantly reducing of the chassis vibration, increasing, therefore, passenger comfort.

Hysteresis modeling of Magneto-Rheological (MR) fluid damper by self tuning fuzzy control

Magneto-rheological (MR) fluid damper is a semi-active control device that has recently received more attention by the vibration control community. But inherent nonlinear hysteresis character of magneto-rheological fluid dampers is one of the challenging aspects for utilizing this device to achieve high system performance. So the development of accurate model is necessary to take the advantage their unique characteristics. Research by others [3] has shown that a system of nonlinear differential equations can successfully be used to describe the hysteresis behavior of the MR damper. The focus of this paper is to develop an alternative method for modeling a damper in the form of centre average fuzzy interference system, where back propagation learning rules are used to adjust the weight of network. The inputs for the model are used from the experimental data. The resulting fuzzy interference system is satisfactorily represents the behavior of the MR fluid damper with reduced computational requirements. Use of the neuro-fuzzy model increases the feasibility of real time simulation.

Modeling of a magneto-rheological (MR) fluid damper using a self tuning fuzzy mechanism

A magneto-rheological (MR) fluid damper is a semi-active control device that has recently begun to receive more attention in the vibration control community. However, the inherent nonlinear nature of the MR fluid damper makes it challenging to use this device to achieve high damping control system performance. Therefore the development of an accurate modeling method for a MR fluid damper is necessary to take advantage of its unique characteristics. Our goal was to develop an alternative method for modeling a MR fluid damper by using a self tuning fuzzy (STF) method based on neural technique. The behavior of the researched damper is directly estimated through a fuzzy mapping system. In order to improve the accuracy of the STF model, a back propagation and a gradient descent method are used to train online the fuzzy parameters to minimize the model error function. A series of simulations had been done to validate the effectiveness of the suggested modeling method when compared with the data measured from experiments on a test rig with a researched MR fluid damper. Finally, modeling results show that the proposed STF interference system trained online by using neural technique could describe well the behavior of the MR fluid damper without need of calculation time for generating the model parameters.