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.

Structural optimization is a useful tool for engineers and designers in construction technology as well asvehicle and mechanical engineering. With structure optimization, a computer can, with the help of finiteelement analysis, calculate the smallest possible amount of material needed to meet the requirements onthe part to be produced.The purpose of this report is to use two different implementations for finite element calculations fortopology optimization of a beam. Results from the optimizations will then be 3D printed with differentsettings. The beam will be tested for displacement, stress and strain in a universal testing machine. Theresults from the experiment will be compared with computed simulations of the same beam.For the structural optimization, two methods are used and compared: Solid Isotropic Material withPenalization and Bidirectional Evolutionary Structural Optimization. A total of eight beams, four fromeach method, were printed with a 3D printer with two different positions on the printer bed and withdifferent degrees of infill ratios. These were tested with a machine that could register both pressure anddeformation and were filmed to be able to see the strain. The deformation of the beams was alsosimulated in a software computer program to see what deformation difference there was betweenexperiment and reality.It turned out that the beams that were printed behaved anisotropic even though solid plastic should beincluded among isotropic materials. The deformation of the model looked like the finite elementcalculation, but the actual deformation was significantly larger than what was calculated by the computersoftware.

Esta investigación se enfoca en el análisis de vibración libre de vigas Timoshenko utilizando el método de elementos finitos. Se desarrolla el modelo utilizando el principio de Hamilton y la teoría de vigas Timoshenko que incluye deformaciones por corte. Se asume interpolaciones de alto orden para la aproximación de las variables fundamentales. Los materiales para emplear son isotrópicos. Se implementa un programa para estos materiales en MATLAB. Se comparan resultados con otros obtenidos en la literatura para validar el modelo. Se realiza un estudio paramétrico con una misma longitud y diferentes esbelteces. Se verifica que la formulación sea bastante precisa con resultados muy satisfactorios. / This research focuses on the free vibration analysis of Timoshkenko beams using the finite element method. The model is developed using the Hamilton principle and the Timoshenko beam theory that includes shear deformations. high order interpolations are assumed for the approximation of the fundamental variables. The materials to be used are isotropic. A program for these materials is implemented in MATLAB. Results are compared with others obtained in the literature to validate the model. A parametric study is carried out with the same length and different slenderness. It is verified that the formulation is quite precise with satisfactory results to the investigation. / Trabajo de investigación