The adhesive contact problems of a conical frustum punch indenting a transversely isotropic elastic half space and an orthotropic elastic half space are analytically studied in this thesis work. To solve the problem involving a transversely isotropic half space, the harmonic potential function method and the Hankel transform are employed, which lead to a general closed-form solution for the adhesive contact problem. For the case with an orthotropic half space, the problem of a point load applied on the half space is first solved by using the double Fourier transform method. The solution for the adhesive contact problem is then obtained through integrating the former solutions over the punch surface.

Shear Wave Elastography (SWE) is an ultrasound based technique which is able to measure tissue stiffness through the speed of induced shear waves. Tissue stiffness is often related to pathological conditions and detecting mechanical changes can help the recognition of potential diseases. The clinical use of SWE is limited to isotropic tissue due to the difficulty in assessing a theoretical model for more complex tissue and this project therefore aimed to evaluate the possibility of obtaining a full mechanical characterization of a transversely isotropic (TI) phantom with dual-probe SWE. A TI hydrogel phantom was developed and mechanical tests were performed to verify its anisotropy and determine the elastic moduli in both the perpendicular and longitudinal directions. Shear moduli were estimated using conventional and dual-probe SWE comparing the results to theoretical pure-transverse (PT) and quasi-transverse (QT) wave propagation modes. Both mechanical and SWE tests showed that the phantoms were transversely isotropic ET/EL=0.81. Moreover, multiple wave propagation modes calculated with dual-probe SWE showed a good agreement with the theoretical curves and indicated the possibility of measuring all the elasticity constants needed to fully characterize an incompressible, TI tissue with dual-probe SWE.

No description available.

Herein we consider material tensors that pertain to thin sheets or thin films, which we model as two-dimensional objects. We assume that the thin sheet in question carries a crystallographic texture characterized by an orientation distribution function defined on the rotation group SO(3), which is almost transversely-isotropic about the sheet normal so that mechanical and physical properties of the thin sheet have weak planar-anisotropy. We present a procedure by which a special orthonormal basis can be determined in each tensor subspace invariant under the action of the orthogonal group O(2). We call members of such special bases irreducible basis tensors under O(2). For the class of thin sheets in question, we derive a representation formula in which each tensor in any given tensor subspace Z is written as the sum of a transversely-isotropic term and a linear combination of orthonormal irreducible basis tensors in Z, where the coefficients are given explicitly in terms of texture coefficients and undetermined material parameters. In addition to the general representation formula, we present also the specialized form for subspaces of tensor products of second-order symmetric tensors, a type commonly found in mechanics of materials.

Orientador: Paulo Sollero / Tese (doutorado) - Universidade Estadual de Campinas. Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-06T09:16:06Z (GMT). No. of bitstreams: 1 Noritomi_PedroYoshito_D.pdf: 4114284 bytes, checksum: 806cfb9e9589c5cd848086ed60178957 (MD5) Previous issue date: 2005 / Resumo: Este trabalho mostra o desenvolvimento de uma metodologia para análise de problemas de bioengenharia, aplicando modelagem numérica elastostática de tensões e deformações, baseada no método dos elementos de contorno com formulação 3D para meios transversalmente isotrópicos lineares, incluindo a capacidade de simulação do comportamento de remodelagem óssea superficial. A implementação do núcleo transversalmente isotrópico sobre a estrutura básica de análise por elementos de contorno 3D utilizou a solução fundamental proposta por Pan & Chou e revisada por Loloi, tendo exigido o cálculo adicional das soluções fundamentais de força de superfície a partir da derivação das soluções fundamentais de deslocamento. O modelo de remodelagem óssea superficial baseou-se na hipótese de estímulo biológico por campo de deformação, partindo de um modelo 2D, adaptado para o espaço 3D com o uso de deformações principais como grandezas de referência. As implementações foram testadas através de análises numéricas de problemas com solução analítica e validações com resultados de aplicações comerciais baseadas em elementos finitos, para problemas padrão de engenharia, bem como comparações com resultados da literatura para problemas de bioengenharia. A análise dos resultados mostrará que, tanto a metodologia quanto as implementações são funcionais, oferecendo uma base sólida para desenvolvimento e teste de novas soluções de bioengenharia / Abstract: This work shows the development of a methodology to analyse bioengineering problems using elastostatic stress-strain numerical modeling based on a 3D transversely isotropic linear boundary element formulation including surface bone remodeling simulation capabilities. The transversely isotropic kernel implementation on the basic 3D boundary element analysis program used the fundamental solution purposed by Pan & Chou and revised by Loloi, with additional fundamental solutions for traction calculation made with the displacement fundamental solution derivatives. The surface bone remodeling model was based on a 2D strain field biological stimulus, extended to the 3D space by using the principal strain as reference values. The implementations were tested through numerical analysis of problems with analytical solution and validation with commercial finite elements applications for standard engineering problems, as well as comparison with literature data for bioengineering problems. The analysis of results will show that both, the methodology and the implementations are fully functional, offering a solid start for development and test of new bioengineering solutions / Doutorado / Mecanica dos Sólidos e Projeto Mecanico / Doutor em Engenharia Mecânica

Fiber-reinforced rubberlike materials (FRRM) commonly used in tires undergo large deformations, and exhibit different response in tension and compression along the fiber direction. Assuming that the response of a fiber-reinforced rubberlike material can be modeled as transversely isotropic with the fiber direction as the axis of transverse isotropy, we express the stored energy function, W, in terms of the five invariants of the right Cauchy-Green strain tensor and the fiber direction, and account for different response in tension and compression along the fiber direction. It has been shown in the literature that in shear-dominated deformations, the 5th invariant, I5, significantly contribution to the stress-strain curve. We have implemented the constitutive relation in the commercial software, LS-DYNA. The numerical solutions of several boundary value problems studied here agree with their analytical solutions derived by using Ericksen's inverse approach, in which a part of the solution is assumed and unknowns in the presumed solution are then found by analyzing the pertinent boundary value problem. However, computed results have not been compared with experimental findings. For W of the FRRMs an expression that is a complete quadratic function of the five invariants is also examined. Homogeneous deformations such as simple extension, simple shear, and biaxial loading problems are studied to delineate the mechanical behaviors of FRRMs. Consistency with the infinitesimal deformation theory requires that linear terms in the 4th and 5th invariants, I4 and I5, be included in the expression for W. Stability analysis of deformations reveals the qualitative changes triggered by the second order terms of the quadratic function. Analytical solutions for inflation, extension and twist deformations caused by internal pressure, end torque, and axial force for a pressurized cylindrical laminate are derived using Ericksen's inverse method. Effects of fiber orientations on the mechanical behaviors of a +/-α angle-ply cylindrical tube are investigated using the derived analytical solutions. The T-peel test, widely used for characterizing adhesion across a plethora of adhesives, adherends, and geometries, results in a range of responses that may complicate meaningful interpretation of the test data. This research effort, involving several specific specimen types, was undertaken to investigate concerns that commonly used configurations may not always result in plateaus in the force-displacement response. We experimentally and numerically study debonding of T-peel specimens having 75 mm bond length and 0.81 mm thick adherends made of either 6061 aluminum (Al) or one of the three steels (G70 70U hot dip galvanized, E60 elctrogalvanized (EGZ), 1010 cold-rolled steel (CRS) bonded with either LORD® 406 or Maxlok™ acrylic adhesive. For the EGZ and the Al adherends, specimens with a bond length of 250 mm and adherend thickness of 1.60 mm are also examined. Effects of adherend materials and thicknesses, bond lengths, and adhesives on test results are examined using three metrics to interpret the T-peel bond performance. We find a limited correlation between the commonly used "T-peel strength" and the energy dissipated per unit debond area. For those two metrics, the relative performances of the CRS and the Al specimens are quite different. Quasi-static plane strain deformations of the test specimens are analyzed by the finite element method (FEM) and a cohesive zone model using the commercial software, ABAQUS, to help interpret the test data. Numerical results provided energies required to elastically and plastically deform the adherends, and help determine the transition from non-self-similar to self-similar debonding. The FE simulations also facilitate determination of the fraction of the crosshead displacement at which self-similar debonding occurs. Results reported herein should help practitioners select appropriate specimen dimensions for extracting meaningful data for adhesive performance. / Ph. D. / Tire belts, seals, and impact absorbing cushions are usually made of fiber-reinforced rubberlike materials (FRRMs), but are difficult to analyze because their response to complex loading situations is strongly dependent on a variety of material properties. Many biological soft tissues, such as tendons, ligaments and arteries are also typically modeled as FRRMs. We assume that a fiber-reinforced rubberlike material can be modeled as nonlinear, incompressible and directionally dependent, with different response in tension and compression along the fiber direction. For such a material, the stored energy functions, W, depends upon five invariant metrics of the imposed strain state and the fiber direction. Explicit expressions for the stresses are derived for two polynomial functions of the five invariants for W. Homogeneous deformations such as simple extension, simple shear, and biaxial loading problems, nonhomogeneous deformations such as plane strain bending of a rectangle beam into a circular one, and inflation, twist and extension of a pressurized cylindrical laminate, are analyzed to reveal the mechanical behaviors descried by the developed material models. To enable the numerical solutions, the developed material models are incorporated in the commercial software, LS-DYNA, as user-defined subroutines. The implementations have been verified by ensuring that the computed solutions of several boundary value problems agree well with the derived analytical solutions or those available in the literature. The work provides theoretical guidelines for using quadratic polynomial functions for material models of FRRM, and delivers the software (user-defined material subroutines) capable of numerically analyzing large deformations of FRRM with different responses in tension and compressions. Large elasto-plastic deformations of T-peel joints have been analyzed using the commercial software, ABAQUS, to delineate conditions that result in self-similar debonding, enabling one to appropriately partition the energy involved in bending the adherends and propagating a debond. Using experimentally measured fracture energies from separate double cantilever beam (DCB) tests, implemented in a traction-separation law, accurate estimates of required peel force, crosshead displacements at break, and plastically deformed peel arm shapes are made. The demonstrated success of predicting load-displacement curves, deformed shape, and various energy metrics by using the traction-separation law in ABAQUS provides us with a framework to use in the future assessment of T-peel configurations being addressed in this study.

No description available.

No description available.

A numerical study of dilute suspensions based on a non-Newtonian matrix fluid and rigid spheroidal particles is performed. A Carreau fluid describes the non-Newtonian matrix. The special case of rigid spherical particles is considered. Here, a uniaxial elongational flow around a sphere is simulated and numerical homogenization is used to obtain the bulk viscosity of the dilute suspension for different applied rates of deformation and different thinning exponents. In the Newtonian regime the well-known Einstein result for the viscosity of a dilute suspension of rigid spherical particles is obtained. In the power-law regime it is found that the intrinsic viscosity depends only on the thinning exponent. Utilizing the simulation results a modification of the Carreau model for dilute suspensions with a non-Newtonian matrix fluid is proposed. To investigate the influence of the particle shape another numerical study is performed. In particular, different flows around spheroidal particles with different orientations are simulated and numerical homogenization is used to obtain the intrinsic viscosity of the suspension as function of applied rate of deformation, thinning exponent and aspect ratio. From the results it is possible to extract the rheological coefficients of the Lipscomb model. In the Newtonian regime the simulation results coincide with Lipscomb’s predictions. In the power-law regime the rheological coefficients depend strongly on the thinning exponent. Furthermore, simulation results indicate that the rheological coefficients additionally depend on the particle orientation in the non-linear regime. / Une étude numérique sur des suspensions diluées à base d’un fluide non newtonien et de particules sphéroïdales rigides est réalisée. Le comportement de la matrice est décrit par un fluide de type Carreau. De particules sphériques et rigides est considéré en premier. Un écoulement en élongation uniaxiale autour d’une sphère est simulée. Ensuite, l’homogénéisation numérique est utilisée pour déterminer la viscosité apparente de la suspension pour différents taux de déformation et d’indices pseudoplastiques. Dans le domaine newtonien, le résultat d’Einstein donnant la viscosité d’une suspension diluée de particules sphériques et rigides est obtenu. Dans le régime en loi de puissance on constate que la viscosité intrinsèque dépend uniquement de l’indice pseudoplastique. Une autre étude numérique est effectuée pour investiguer l’influence de la forme des particules. Plusieurs écoulements autour d’une particule sphéroïdale sont simulés pour différentes orientations. Une homogénéisation numérique est ensuite utilisée pour obtenir la viscosité intrinsèque de la suspension en fonction du taux de déformation appliqué, de l’indice d’écoulement et du rapport de forme de la particule. A partir de ces résultats, il est possible d’exprimer les coefficients rhéologiques du modèle de Lipscomb. Dans le régime newtonien, les résultats coïncident avec les prédictions de Lipscomb. Dans le domaine en loi de puissance, les coefficients rhéologiques deviennent fortement dépendent de l’indice pseudoplastique. En outre, les résultats des simulations montrent que ces coefficients rhéologiques dépendent également de l’orientation des particules dans le régime non linéaire. / Numerische Untersuchung zu verdünnten Suspensionen basierend auf einer nicht Newtonschen Matrixflüssigkeit und harten spheroidalen Partikeln wurde durchgeführt. Ein Carreau Fluid beschreibt die nicht Newtonsche Matrix. Zuerst wird der Spezialfall harter Kugeln betrachtet. Hierzu wird eine uniaxiale Dehnströmung um eine Kugel simuliert und numerische Homogenisierung wird verwendet um die effektive Viskosität der Suspension für verschieden aufgebrachte Deformationsgeschwindigkeiten und Verdünnungsexponenten zu bestimmen. Im Newtonschen Bereich wird die bekannte Lösung Einsteins für die Viskosität einer verdünnten Suspension harter Kugeln erhalten. Im power-law Bereich ist die intrinsische Viskosität einzig eine Funktion des Verdünnungsexponenten. Unter Nutzung der Simulationsergebnisse wird eine Modifikation des Carreau Modells vorgeschlagen. Um den Einfluss der Partikelform auf die nichtlinearen Eigenschaften zu untersuchen wird eine weitere numerische Simulationen durchgeführt. Dabei werden verschiedene Strömungen um spheroidale Partikel mit unterschiedlicher Orientierung simuliert und numerische Homogenisierung wird verwendet um die intrinsische Viskosität als Funktion der aufgebrachten Deformationsgeschwindigkeit, des Verdünnungsexponenten und des Partikelaspektverhältnisses zu bestimmen. Es ist möglich die rheologischen Parameter des Lipscomb Modells aus den Simulationsergebnissen zu bestimmen. Im Newtonschen Bereich stimmen die numerisch bestimmten Werte mit der Vorhersage Lipscomb‘s überein. Im power-law Bereich hängen die rheologischen Parameter stark vom Verdünnungsexponenten ab. Weiter kann man aus den Ergebnissen auf eine zusätzliche Abhängigkeit der rheologischen Parameter von der Partikelorientierung schließen.

A numerical study of dilute suspensions based on a non-Newtonian matrix fluid and rigid spheroidal particles is performed. A Carreau fluid describes the non-Newtonian matrix. The special case of rigid spherical particles is considered. Here, a uniaxial elongational flow around a sphere is simulated and numerical homogenization is used to obtain the bulk viscosity of the dilute suspension for different applied rates of deformation and different thinning exponents. In the Newtonian regime the well-known Einstein result for the viscosity of a dilute suspension of rigid spherical particles is obtained. In the power-law regime it is found that the intrinsic viscosity depends only on the thinning exponent. Utilizing the simulation results a modification of the Carreau model for dilute suspensions with a non-Newtonian matrix fluid is proposed. To investigate the influence of the particle shape another numerical study is performed. In particular, different flows around spheroidal particles with different orientations are simulated and numerical homogenization is used to obtain the intrinsic viscosity of the suspension as function of applied rate of deformation, thinning exponent and aspect ratio. From the results it is possible to extract the rheological coefficients of the Lipscomb model. In the Newtonian regime the simulation results coincide with Lipscomb’s predictions. In the power-law regime the rheological coefficients depend strongly on the thinning exponent. Furthermore, simulation results indicate that the rheological coefficients additionally depend on the particle orientation in the non-linear regime. / Une étude numérique sur des suspensions diluées à base d’un fluide non newtonien et de particules sphéroïdales rigides est réalisée. Le comportement de la matrice est décrit par un fluide de type Carreau. De particules sphériques et rigides est considéré en premier. Un écoulement en élongation uniaxiale autour d’une sphère est simulée. Ensuite, l’homogénéisation numérique est utilisée pour déterminer la viscosité apparente de la suspension pour différents taux de déformation et d’indices pseudoplastiques. Dans le domaine newtonien, le résultat d’Einstein donnant la viscosité d’une suspension diluée de particules sphériques et rigides est obtenu. Dans le régime en loi de puissance on constate que la viscosité intrinsèque dépend uniquement de l’indice pseudoplastique. Une autre étude numérique est effectuée pour investiguer l’influence de la forme des particules. Plusieurs écoulements autour d’une particule sphéroïdale sont simulés pour différentes orientations. Une homogénéisation numérique est ensuite utilisée pour obtenir la viscosité intrinsèque de la suspension en fonction du taux de déformation appliqué, de l’indice d’écoulement et du rapport de forme de la particule. A partir de ces résultats, il est possible d’exprimer les coefficients rhéologiques du modèle de Lipscomb. Dans le régime newtonien, les résultats coïncident avec les prédictions de Lipscomb. Dans le domaine en loi de puissance, les coefficients rhéologiques deviennent fortement dépendent de l’indice pseudoplastique. En outre, les résultats des simulations montrent que ces coefficients rhéologiques dépendent également de l’orientation des particules dans le régime non linéaire. / Numerische Untersuchung zu verdünnten Suspensionen basierend auf einer nicht Newtonschen Matrixflüssigkeit und harten spheroidalen Partikeln wurde durchgeführt. Ein Carreau Fluid beschreibt die nicht Newtonsche Matrix. Zuerst wird der Spezialfall harter Kugeln betrachtet. Hierzu wird eine uniaxiale Dehnströmung um eine Kugel simuliert und numerische Homogenisierung wird verwendet um die effektive Viskosität der Suspension für verschieden aufgebrachte Deformationsgeschwindigkeiten und Verdünnungsexponenten zu bestimmen. Im Newtonschen Bereich wird die bekannte Lösung Einsteins für die Viskosität einer verdünnten Suspension harter Kugeln erhalten. Im power-law Bereich ist die intrinsische Viskosität einzig eine Funktion des Verdünnungsexponenten. Unter Nutzung der Simulationsergebnisse wird eine Modifikation des Carreau Modells vorgeschlagen. Um den Einfluss der Partikelform auf die nichtlinearen Eigenschaften zu untersuchen wird eine weitere numerische Simulationen durchgeführt. Dabei werden verschiedene Strömungen um spheroidale Partikel mit unterschiedlicher Orientierung simuliert und numerische Homogenisierung wird verwendet um die intrinsische Viskosität als Funktion der aufgebrachten Deformationsgeschwindigkeit, des Verdünnungsexponenten und des Partikelaspektverhältnisses zu bestimmen. Es ist möglich die rheologischen Parameter des Lipscomb Modells aus den Simulationsergebnissen zu bestimmen. Im Newtonschen Bereich stimmen die numerisch bestimmten Werte mit der Vorhersage Lipscomb‘s überein. Im power-law Bereich hängen die rheologischen Parameter stark vom Verdünnungsexponenten ab. Weiter kann man aus den Ergebnissen auf eine zusätzliche Abhängigkeit der rheologischen Parameter von der Partikelorientierung schließen.