A Viscoelastic Constitutive Model with Nonlinear Evolutionary Internal Variables

Based on the internal variable theory, a viscoelastic constitutive model of a highly deformable continuous medium is proposed. A set of second rank tensorial internal state variables corresponding to Biot's strain is introduced, and a nonlinear evolution

Discrete element methods for asphalt concrete : development and application of user-defined microstructural models and a viscoelastic micromechanical model

As an important Civil Engineering material, asphalt concrete (AC) is commonly used to build road surfaces, airports, and parking lots. With traditional laboratory tests and theoretical equations, it is a challenge to fully understand such a random composite material. Based on the discrete element method (DEM), this research seeks to develop and implement computer models as research approaches for improving understandings of AC microstructure-based mechanics. In this research, three categories of approaches were developed or employed to simulate microstructures of AC materials, namely the randomly-generated models, the idealized models, and image-based models. The image-based models were recommended for accurately predicting AC performance, while the other models were recommended as research tools to obtain deep insight into the AC microstructure-based mechanics. A viscoelastic micromechanical model was developed to capture viscoelastic interactions within the AC microstructure. Four types of constitutive models were built to address the four categories of interactions within an AC specimen. Each of the constitutive models consists of three parts which represent three different interaction behaviors: a stiffness model (force-displace relation), a bonding model (shear and tensile strengths), and a slip model (frictional property). Three techniques were developed to reduce the computational time for AC viscoelastic simulations. It was found that the computational time was significantly reduced to days or hours from years or months for typical three-dimensional models. Dynamic modulus and creep stiffness tests were simulated and methodologies were developed to determine the viscoelastic parameters. It was found that the DE models could successfully predict dynamic modulus, phase angles, and creep stiffness in a wide range of frequencies, temperatures, and time spans. Mineral aggregate morphology characteristics (sphericity, orientation, and angularity) were studied to investigate their impacts on AC creep stiffness. It was found that aggregate characteristics significantly impact creep stiffness. Pavement responses and pavement-vehicle interactions were investigated by simulating pavement sections under a rolling wheel. It was found that wheel acceleration, steadily moving, and deceleration significantly impact contact forces. Additionally, summary and recommendations were provided in the last chapter and part of computer programming codes wree provided in the appendixes.

Numerical simulation of abrupt contraction flows using the Double Convected Pom-Pom model

The Double Convected Pom-Pom model was recently introduced to circumvent some numerical and theological defects found in other formulations of the Pom-Pom concept. It is used here for the simulation of a benchmark problem: the flow in an abrupt planar contraction. The predictions are compared with birefringence measurements and show reasonable quantitative agreement with experimental data. A parametric study is also carried out with the aim of analysing the effect of the branching parameter on vortex dynamics and extrudate swell. The results show that the Double Convected Pom-Pom model (DCPP) model is able to discriminate between branched and linear macromolecular structures in accordance with experimental observations. In that respect, the role of the extensional properties in determining complex flow behaviour is stressed. Also, the ratio of the first normal stress difference to the shear stress appears to play a major role in die swell observation. For the time being, the role of the second normal stress difference appears to be less obvious to evaluate in this complex flow. (C) 2004 Elsevier B.V. All rights reserved.

Measuring the compressive viscoelastic mechanical properties of human cervical tissue using indentation.

The human cervix is an important mechanical barrier in pregnancy which must withstand the compressive and tensile forces generated from the growing fetus. Premature cervical shortening resulting from premature cervical remodeling and alterations of cervical material properties are known to increase a woman׳s risk of preterm birth (PTB). To understand the mechanical role of the cervix during pregnancy and to potentially develop indentation techniques for in vivo diagnostics to identify women who are at risk for premature cervical remodeling and thus preterm birth, we developed a spherical indentation technique to measure the time-dependent material properties of human cervical tissue taken from patients undergoing hysterectomy. In this study we present an inverse finite element analysis (IFEA) that optimizes material parameters of a viscoelastic material model to fit the stress-relaxation response of excised tissue slices to spherical indentation. Here we detail our IFEA methodology, report compressive viscoelastic material parameters for cervical tissue slices from nonpregnant (NP) and pregnant (PG) hysterectomy patients, and report slice-by-slice data for whole cervical tissue specimens. The material parameters reported here for human cervical tissue can be used to model the compressive time-dependent behavior of the tissue within a small strain regime of 25%.

Técnica para dedução do modelo algébrico PTT, aplicações e análises numéricas em escoamentos bidimensionais

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

Mechanical Alterations in airway smooth muscle after interaction with indoor pollutant – A mathematical model.

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully understood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholine (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle ¯bers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, incorporates the phenomenology of both MCh and FA and reproduces experimental results observed with in vitro exposure of smooth muscle to FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells in a tissue level model. The model can also be used in different biological scales.

Mechanical alterations in airway smooth muscle after interaction with indoor pollutant - A mathematical model

The viscoelasticity of mammalian lung is determined by the mechanical properties and structural regulation of the airway smooth muscle (ASM). The exposure to polluted air may deteriorate these properties with harmful consequences to individual health. Formaldehyde (FA) is an important indoor pollutant found among volatile organic compounds. This pollutant permeates through the smooth muscle tissue forming covalent bonds between proteins in the extracellular matrix and intracellular protein structure changing mechanical properties of ASM and inducing asthma symptoms, such as airway hyperresponsiveness, even at low concentrations. In the experimental scenario, the mechanical effect of FA is the stiffening of the tissue, but the mechanism behind this effect is not fully w1derstood. Thus, the aim of this study is to reproduce the mechanical behavior of the ASM, such as contraction and stretching, under FA action or not. For this, it was created a two-dimensional viscoelastic network model based on Voronoi tessellation solved using Runge-Kutta method of fourth order. The equilibrium configuration was reached when the forces in different parts of the network were equal. This model simulates the mechanical behavior of ASM through of a network of dashpots and springs. This dashpot-spring mechanical coupling mimics the composition of the actomyosin machinery of ASM through the contraction of springs to a minimum length. We hypothesized that formation of covalent bonds, due to the FA action, can be represented in the model by a simple change in the elastic constant of the springs, while the action of methacholinc (MCh) reduce the equilibrium length of the spring. A sigmoid curve of tension as a function of MCh doses was obtained, showing increased tension when the muscle strip was exposed to FA. Our simulations suggest that FA, at a concentration of 0.1 ppm, can affect the elastic properties of the smooth muscle fibers by a factor of 120%. We also analyze the dynamic mechanical properties, observing the viscous and elastic behavior of the network. Finally, the proposed model, although simple, ir1corporates the phenomenology of both MCh and FA and reproduces experirnental results observed with ir1 vitro exposure of smooth muscle to .FA. Thus, this new mechanical approach incorporates several well know features of the contractile system of the cells ir1 a tissue level model. The model can also be used in different biological scales.

Design and development of new pressure sensors for aerodynamic applications

This artwork reports on two different projects that were carried out during the three years of Doctor of the Philosophy course. In the first years a project regarding Capacitive Pressure Sensors Array for Aerodynamic Applications was developed in the Applied Aerodynamic research team of the Second Faculty of Engineering, University of Bologna, Forlì, Italy, and in collaboration with the ARCES laboratories of the same university. Capacitive pressure sensors were designed and fabricated, investigating theoretically and experimentally the sensor’s mechanical and electrical behaviours by means of finite elements method simulations and by means of wind tunnel tests. During the design phase, the sensor figures of merit are considered and evaluated for specific aerodynamic applications. The aim of this work is the production of low cost MEMS-alternative devices suitable for a sensor network to be implemented in air data system. The last two year was dedicated to a project regarding Wireless Pressure Sensor Network for Nautical Applications. Aim of the developed sensor network is to sense the weak pressure field acting on the sail plan of a full batten sail by means of instrumented battens, providing a real time differential pressure map over the entire sail surface. The wireless sensor network and the sensing unit were designed, fabricated and tested in the faculty laboratories. A static non-linear coupled mechanical-electrostatic simulation, has been developed to predict the pressure versus capacitance static characteristic suitable for the transduction process and to tune the geometry of the transducer to reach the required resolution, sensitivity and time response in the appropriate full scale pressure input A time dependent viscoelastic error model has been inferred and developed by means of experimental data in order to model, predict and reduce the inaccuracy bound due to the viscolelastic phenomena affecting the Mylar® polyester film used for the sensor diaphragm. The development of the two above mentioned subjects are strictly related but presently separately in this artwork.

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.

Andaman Postseismic Deformation Observations: Still Slipping after All These Years?

More than six years after the great (M-w 9.2) Sumatra-Andaman earthquake, postevent processes responsible for relaxation of the coseismic stress change remain controversial. Modeling of Andaman Islands Global Positioning System (GPS) displacements indicated early near-field motions were dominated by slip down-dip of the rupture, but various researchers ascribe elements of relaxation to dominantly poroelastic, dominantly viscoelastic, and dominantly fault slip processes, depending primarily on their measurement sampling and modeling tools used. After subtracting a pre-2004 interseismic velocity, significant transient motion during the 2008.5-2010.5 epoch confirms that postseismic relaxation processes continue in Andaman. Modeling three-component velocities as viscoelastic flow yields a weighted root-mean-square (wrms) misfit that always exceeds the wrms of the measured signal (26.3 mm/yr). The best-fitting models are those that yield negligible deformation, indicating the model parameters have no real physical meaning. GPS velocities are well fit (wrms 4.0 mm/yr) by combining a viscoelastic flow model that best fits the horizontal velocities with similar to 50 cm/yr thrust slip down-dip of the coseismic rupture. Both deep slip and flow respond to stress changes, and each can significantly change stress in the realm of the other; it therefore is reasonable to expect that both transient deep slip and viscoelastic flow will influence surface deformation long after a great earthquake.

心肌细胞大变形黏弹性模型及在实验中的应用

在衡量单个细胞力学行为的研究中,越来越多地采用结合实验的数值模拟方法.在连续介质力学框架下,发展了一种新的心肌细胞本构模型,并与微管吮吸实验结合,探讨了心肌细胞的力学特性.本构模型是对普遍使用的仅能用于小变形分析的标准线性固体模型的一种扩展,它将超弹性性能引入到黏弹性模型中,用以描述细胞的大变形黏弹性效应.基于改进的本构模型,对心肌细胞微管吮吸实验过程进行了有限元模拟,并将计算结果与实验结果以及经典理论解进行了对比.结果显示发展的本构模型适合细胞大变形问题的有限元数值模拟.

Análise de barragens de gravidade de concreto considerando a fase construtiva e a interação dinâmica barragem-reservatório-fundação

o objetivo deste trabalho é a análise de barragens de gravidade de concreto desde a faseda sua construção até sua completa entrada em serviço. Inicialmente é feita a análise da fase construtiva, onde o problema fundamental é devido às tensões térmicas decorrentes do calor de hidratação. O método dos elementos finitos é empregado para a solução dos problemasde transferência de calor e de tensões. A influência da construção em camadas é introduzidaatravés da redefinição da malha de elementos finitos, logo após o lançamento de cadacamada de concreto. Uma atenção especial é dada ao problema de fissuração em estruturas de concreto simples.Algunsmodelos usuais são apresentados, discutindo-se a eficiência dos mesmos. Os modelosde fissuração distribuída têm sido preferidos, em virtude dos vários inconvenientes apresentados pelas formulações discretas. Esses modelos, entretanto, fornecem resultados dependentesda malha de elementos finitos e alguma consideração adicional deve ser feita para corrigiressas distorções. Normalmente, tenta-se corrigir esse problema através da adoção de umaresistênciaà tração minorada que é definida em função da energia de fratura do material. Neste trabalho, é demonstrado que esse procedimento não é satisfatório e é proposta uma novaformulaçãopara a análise de grandes estruturas de concreto. A análise das tensões na etapa de construção da barragem é feita com o emprego de um modelo constitutivo viscoelástico com envelhecimento para o concreto. Em virtude do envelhecimento,a matriz de rigidez da estrutura é variável no tempo, devendo ser redefinida e triangularizadaem cada instante. Isto leva a um grande esforço computacional, sobretudo, quandoa barragem é construída em muitas camadas. Para evitar esse inconveniente, adota-se um procedimento iterativo que permite que a matriz de rigidez seja redefinida em poucas idadesde referência. Numa segunda etapa da análise, a barragem é submetida à pressão hidrostática e a uma excitação sísmica. A análise dinâmica é realizada considerando-se o movimento do sistema acoplado barragem-reservatório-fundação. O sismo é considerado um processo estocásticonão estacionário e a segurança da estrutura é determinada em relação aos principais modos de falha