Numerical solution of differential equations having multitude of scales in the solution field is one of the most challenging research areas, but highly demanded in scientific and industrial applications. One of the natural approaches for handling such problems is to separate the scales and approximate the solution of the segregated scales with appropriate numerical method. Variational multiscale method (VMS) is a predominant method in the paradigm of scale separation schemes. In our work we have used the VMS technique to develop a numerical scheme for computations of turbulent flows in time-dependent domains. VMS allows separation of the entire range of scales in the flow field into two or three groups, thereby enabling a different numerical treatment for the different groups. In the context of computational fluid dynamics(CFD), VMS is a significant new improvement over the classical large eddy simulation (LES). VMS does away with the commutation errors arising due to filtering in LES. Further, in a three-scale VMS approach the model for the subgrid scale can be contained to only a part of the resolved scales instead of effecting the entire range of resolved scales. The projection based VMS scheme that we have developed gives a robust and efficient method for solving problems of turbulent fluid flows in deforming domains, governed by incompressible Navier {Stokes equations. In addition to the existing challenges due to turbulence, the computational complexity of the problem increases further when the considered domain is time-dependent. In this work, we have used an arbitrary Lagrangian-Eulerian (ALE) based VMS scheme to account for the domain deformation. In the proposed scheme, the large scales are represented by an additional tensor valued space. The resolved large and small scales are computed in a single unified equation, and the effect of unresolved scales is confined only to the resolved small scales, by using a projection operator. The popular Smagorinsky eddy viscosity model is used to approximate the effects of unresolved scales. The used ALE approach consists of an elastic mesh update technique. Moreover, a computationally efficient scheme is obtained by the choice of orthogonal finite element basis function for the resolved large scales, which allows to reformulate the ALE-VMS system matrix into the standard form of the NSE system matrix. Thus, any existing Navier{Stokes solver can be utilized for this scheme, with modifications. Further, the stability and error estimates of the scheme using a linear model of the NSE are also derived. Finally, the proposed scheme has been validated by a number of numerical examples over a wide range of problems.

Este trabalho trata do estudo de robôs de arquitetura paralela, focando na modelagem e otimização dos mesmos. Não foi construído nenhum tipo de protótipo físico, contudo os modelos virtuais poderão, no futuro, habilitar tal façanha. Após uma busca por uma aplicação que se beneficie do uso de um robô de arquitetura paralela, fez-se uma pesquisa por arquiteturas viáveis já existentes ou relatadas na literatura. Escolheu-se a mais apta e prosseguiu-se com os estudos e modelagem cinemática e dinâmica, dando uma maior ênfase na cinemática e dinâmica inversa, esta última utilizando a formulação de Newton - Euler. Foi construído um simulador virtual em ambiente MATLAB 6.5, dotado de várias capacidades como interpolação linear e circular, avanço e uso de múltiplos eixos coordenados. Seu propósito principal é o de demonstrar a funcionalidade e eficácia dos métodos utilizados. Depois foi incorporado ao simulador um algoritmo de cálculo do volume de trabalho da máquina que utiliza alguns dados do usuário para calcular o volume, que pode ser aquele atrelado a uma postura em particular ou o volume de trabalho de orientação total. Algoritmos para medir o desempenho da máquina quanto à uniformidade e utilização da força dos atuadores foram construídos e também incorporados ao simulador, que consegue mostrar o elipsóide de forças ao longo de quaisquer movimentos executados pela plataforma móvel. Quanto à otimização, parte do ferramental previamente construído foi utilizado para que se pudesse chegar a um modelo de uma máquina que respeitasse restrições mínimas quanto ao tamanho e forma de seu volume de trabalho, mas ainda mantendo o melhor desempenho possível dentro deste volume. / This work is about the study of parallel architecture robots, focusing in modeling and optimization. No physical prototypes were built, although the virtual models can help those willing to do so. After searching for an application that could benefit from the use of a parallel robot, another search was made, this time for the right architecture type. After selecting the architecture, the next step was the kinematics and dynamics analysis. The dynamics model is developed using the Newton ? Euler method. A virtual simulator was also developed in MATLAB 6.5 environment. The simulator?s main purpose was to demonstrate that the methods applied were correct and efficient, so it has several features such as linear and circular interpolations, capacity to use multiple coordinate systems and others. After finishing the simulator, an algorithm to calculate the machine workspace was added. The algorithm receives as input some desired requirements regarding the manipulator pose and then calculates the workspace, taking into consideration imposed constraints. Lastly, algorithms capable to measure the manipulator?s performance regarding to its actuator and end-effector force relationship were also incorporated into the simulator that calculates the machine?s force ellipsoid during any movement, for each desired workspace point. For the optimization procedures, some previously developed tools were used, so that the resulting model was capable to respect some workspace constraints regarding size and shape, but also maintaining the best performance possible inside this volume.

Este trabalho propõe estudar o projeto de estruturas submetidas a grandes deslocamentos utilizando o Método de Otimização Topológica (MOT). O MOT é um método numérico capaz de fornecer de forma sistemática a distribuição ótima de material no domínio de uma estrutura de forma a atender a um dado requisito de projeto, por exemplo, o valor de flexibilidade máxima permitida em uma estrutura. Desde sua introdução, há quase três décadas, o MOT ganhou popularidade na área acadêmica e na indústria. Até o presente momento (2011), a maioria dos trabalhos relacionados com o método tem se preocupado com a otimização de estruturas com o comportamento linear, ou seja, pequenos deslocamentos. Um pequeno número de artigos e trabalhos tem sido relacionado com a modelagem e otimização topológica de estruturas submetidas a efeitos não-lineares. Este trabalho propõe compilar as formulações descritas na literatura e agregar novas técnicas na implementação da OT de forma a melhorar a robustez na obtenção de resultados sob não-linearidade geométrica. O MOT para o comportamento não-linear geométrico neste trabalho foi implementado utilizando o modelo de material SIMP. O comportamento não-linear geométrico é representado utilizando a formulação Lagrangiana para as leis de material de Kirchhoff-Saint Venant e neo-Hookiana. Ambas as leis de material foram implementadas utilizando o método de elementos finitos (MEF) e o equilíbrio estático da estrutura é obtido através de uma rotina incremental e iterativa de Newton incluindo todos os elementos (inclusive os de baixa densidade) dentro do domínio de projeto. A sensibilidade da função objetivo é deduzida utilizando o método adjunto e o problema de otimização é resolvido utilizando o Método das Assíntotas Móveis (MAM) em conjunto com uma função de Relaxação proposta para estabilizar a solução de OT não-linear. A função de projeção não-linear em conjunto com o Método da Continuação é utilizada para eliminar o problema de tabuleiro e independência de malha, melhorando a convergência dos resultados. A função objetivo para minimização da flexibilidade no ponto de aplicação do carregamento é testada, considerando um carregamento fixo. Neste trabalho, os exemplos mostram que as diferenças na rigidez das estruturas otimizadas utilizando modelagem linear e não-linear são geralmente pequenas para pequenos carregamentos, mas elas podem ser grandes em certos casos envolvendo grandes cargas, acarretando em instabilidades na estrutura, o que pode degenerar a solução obtida. / This work proposes studying the design of structures undergoing large displacement using Topology Optimization Method (TOM). The TOM is a numerical method capable of synthesizing the basic layout of a mechanical structure accomplishing to a given design requirement, for example the maximum strain energy allowed in the structure. Since its introduction nearly three decades, TOM has gained widespread popularity in academia and industry. So far, most papers dealing with the method have been concerned with the optimization of structures with linear geometric and material behavior. Even now a small number of works and articles have been concerned with the modeling and topology optimization of structures undergoing nonlinear effects. This work proposes to compile the formulations described in the literature and adding new techniques to improve the robustness for obtaining results of OT under geometric nonlinearity. The TOM for geometric nonlinear behavior in this work is implemented with Solid Isotropic Microstructure with Penalization (SIMP) material model. The geometrically nonlinear behavior of the structures is modeled using a Lagrangean description for hyperelastic constitutive models for Saint Venant-Kirchhoff and neo-Hookean. Both constitutive models are implemented using the Finite Element Method (FEM) and the static equilibrium of the structure is obtained using an incremental and iterative Full-Newton Method considering all elements and internal force of the design domain (elements called \"voids\"). The sensitivity of the objective function is derived using the adjoint method and the optimization problem is solved using the Optimality Criteria (OC) method and Method of Moving Asymptotes (MMA) together with a Relaxation Function proposed to stabilize the TO nonlinear solution. The nonlinear projection function in conjunction with the Continuation Method is used to obtain checkerboard-free and mesh-independent designs and to improve the convergence results. The objective function of end-compliance is tested, by minimizing it for a fixed load. In this work, some examples show that differences in stiffness of optimized structures using linear and nonlinear modeling are generally small, however they can be large in certain cases involving buckling or bifurcation point, that degenerate the solution obtained.

O trabalho trata da utilização de formulações não-convencionais de elementos finitos na obtenção de fatores de intensidade de tensão associados a múltiplas fissuras distribuídas num domínio bidimensional. A formulação do problema de múltiplas fissuras baseia-se numa abordagem de sobreposição proposta pelo Método da Partição (\"Splitting Method\"). Segundo essa abordagem a solução do problema pode ser encontrada a partir da sobreposição de três subproblemas combinados de tal forma que o fluxo de tensão resultante nas faces das fissuras seja nulo. O uso do Método dos Elementos Finitos (MEF) em sua forma convencional pode requerer um refinamento excessivo da rede nesse tipo de problema, aumentando o custo computacional da análise. Objetivando reduzir este custo, empregam-se duas formulações não-convencionais, de forma independente, num dos subproblemas, dito local: a formulação híbrido-Trefftz de tensão e o Método dos Elementos Finitos Generalizados (MEFG). Na formulação híbrido-Trefftz é adotado o recurso do enriquecimento seletivo mediante o refrno- p na aproximação dos campos de deslocamento no contorno do elemento. Já com relação ao MEFG, empregam-se funções polinomiais e a solução analítica da mecânica da fratura como funções enriquecedoras. Exemplos de simulação numérica são apresentados no sentido de comprovar que a utilização dessas formulações não-convencionais juntamente com o Método da Partição viabiliza a obtenção de resultados com boa aproximação com recurso a redes pouco refinadas, reduzindo significativamente o custo computacional de toda a análise. / This paper treats with the use of non-conventional finite element formulations to obtain the stress intensity factor of multiple cracks located in a two-dimensional domain. The formulation of the multiple cracks problem is based on an overlapping approach suggested by the Splitting Method. Accordingly, the solution of the problem can be achieved by dividing the problem in three steps, combined so that the resulting stress flux is zero on the cracks face. The use of the Finite Element Method (FEM) in its conventional formulation requires a mesh refinement in this kind of problem, then increasing the computational cost. Aiming to reduce this cost, two non-conventional formulations are used independently to solve the local problem: the Hybrid-Trefftz stress formulation and the Generalized Finite Elements Method (GFEM). The Hybrid-Trefftz formulation is applied with selective enrichment using p-refinement in the displacements field on the element boundaries. The GFEM employs polynomial functions and analytical solutions of the fracture mechanics as enrichment functions. Examples of numerical simulations are presented in order to show that non- conventional formulations and the Splitting Method can provide accurate results with coarse mesh, thus reducing the computational cost.

O objetivo central deste trabalho é implementar modelos constitutivos hiperelásticos não lineares em um código computacional que faz análise não linear geométrica de cascas. São necessários, para este propósito, conceitos sobre álgebras linear e tensorial, cinemática, deformação, tensão, balanços, princípios variacionais, métodos numéricos e hiperelasticidade. Tal programa usa a formulação Lagrangiana posicional, o método dos elementos finitos, o princípio dos trabalhos virtuais e o método iterativo de Newton-Raphson para solução das equações não lineares. O elemento finito de casca possui dez nós, sete parâmetros por nó e variação linear da deformação ao longo da espessura. Para dedução dos novos modelos usou-se a decomposição multiplicativa do gradiente da função mudança de configuração, o tensor deformação de Green-Lagrange e o tensor da tensão de Piola-Kirchhoff de segunda espécie. O código desenvolvido foi usado em simulações de diversos exemplos e apresentou boa precisão na análise mecânica de polímeros naturais altamente deformáveis. A ocorrência do fenômeno travamento não se manifestou nas análises realizadas. A presente pesquisa confirmou outros trabalhos, reforçou a necessidade de se usar modelos hiperelásticos não lineares para simular o comportamento mecânico de polímeros naturais e apresentou resultados condizentes com dados experimentais existentes na literatura científica e às respectivas soluções analíticas. / The main objective of this work is to implement nonlinear hyperelastic constitutive models in a computational code of geometrically nonlinear analysis of shells. For this purpose, concepts of linear and tensor algebras, kinematics, strain, stress, balances, variational principles, numerical methods and hyperelasticity are necessary. Such program uses the positional Lagrangian formulation, the finite element method, the principle of virtual work and the iterative method of Newton-Raphson for the solution of the nonlinear equations. The shell finite element has ten nodes, seven parameters per node and presents linear variation of the strain along the thickness. To achieve the new constitutive models the multiplicative decomposition of the deformation gradient, the Green-Lagrange strain tensor and the second Piola-Kirchhoff stress tensor are used. The developed code is tested for simulations of various examples and presents good accuracy in the mechanical analysis of highly deformable natural rubber. The locking phenomena didn\'t appear in the proposed analysis. The present research confirms other works, corroborates the need of using nonlinear hyperelastic models to simulate the mechanical behavior of natural rubber and presents suitable results when compared to existent experimental data of the scientific literature and to the respective analytical solutions.

Este trabalho trata da análise da propagação de fissuras, independente do tempo, em domínios bidimensionais utilizando uma formulação alternativa do método dos elementos de contorno (MEC). O MEC vem sendo utilizado com sucesso na análise de diversos problemas de engenharia. Considerando problemas de mecânica da fratura, o MEC é especialmente eficiente devido à redução da dimensionalidade de sua malha, o que permite a simulação do crescimento das fissuras sem as dificuldades do processo de remalhamento. Nesta pesquisa, desenvolvem-se formulações não lineares do MEC para a análise da propagação de fissuras em materiais quase frágeis. Nesses materiais, a zona de processo à frente da ponta da fissura introduz efeitos fisicamente não lineares no comportamento estrutural. Assim, para a simulação da presença da zona de processo, modelos não lineares são necessários. Classicamente a formulação dual do MEC é utilizada para modelar propagação de fissuras na quais equações singulares e hipersingulares são escritas para elementos definidos ao longo das faces das fissuras. O presente trabalho propõe uma segunda formulação utilizando um campo de tensões iniciais para a representação da zona coesiva. Nesta formulação, o termo de domínio da equação integral clássica do MEC é degenerado, de forma a atuar somente ao longo do caminho de crescimento das fissuras, sendo que esse procedimento dá origem a uma nova variável denominada dipolo, responsável por garantir o atendimento das condições de contorno. Em conjunto com essa nova formulação, se propõe o uso do operador tangente (OT), que é deduzido no trabalho, a fim de acelerar o processo de convergência da solução. Os resultados obtidos, por meio da formulação alternativa, são comparados tanto com dados experimentais quanto com o MEC dual, ambos disponíveis na literatura. As respostas encontradas foram satisfatórias no sentido de conseguir reproduzir o comportamento real da estrutura explorando as vantagens computacionais proporcionadas pelo OT. / This work presents a time-independent crack propagation analysis, in two-dimensional domains, using an alternative boundary element method (BEM) formulation. BEM has been used successfully to analyze several engineering problems. Considering fracture mechanics problems, BEM is especially efficient due to its mesh reduction aspects, which allows the simulation of crack growth without remeshing difficulties. In this research, nonlinear BEM formulations are develop in order to analyze crack propagation in quasi-brittle materials. Considering these materials, the process zone ahead of the crack tip leads to nonlinear effects related to structural behavior. Thus, nonlinear models are required for simulating the presence of the process zone. Classically, the dual BEM is used for modeling the crack propagation, in which singular and hyper-singular equations are written for elements defined along the crack faces. This work proposes an alternative formulation using the initial stress field to represent the cohesive zone. In this formulation, the classic domain integral term is degenerated in order to be non-null only at the crack growth path. This procedure leads the creation of new variable called dipole, which is responsible for ensuring the compliance of the boundary conditions. In addition to this new formulation, it is proposed the use of the tangent operator (TO), which is derived in this work, in order to accelerate the convergence. The results obtained using the new formulation, are compared with experimental data and dual BEM results available in the literature. The responses were found satisfactory in reproducing the behavior of real structures exploiting the computational advantages provided by the TO.

Este trabalho trata da combinação entre a formulação híbrida-mista de tensão (FHMT) (Freitas et al. (1996)), para a elasticidade plana, com o método dos elementos finitos generalizados (MEFG), Duarte et al. (2000). O MEFG se caracteriza como uma forma não-convencional do método dos elementos finitos (MEF) que resulta da incorporação a este de conceitos e técnicas dos métodos sem malha, como o enriquecimento nodal proposto do método das nuvens “hp”. Como na FHMT são aproximados dois campos no domínio (tensão e deslocamento) e um no contorno (deslocamento), diferentes possibilidades de enriquecimento nodal são exploradas. Para a discretização do modelo híbrido-misto empregam-se elementos finitos quadrilaterais com funções de forma bilineares para o domínio e elementos lineares para o contorno. Essas funções são enriquecidas por funções polinomiais, trigonométricas, polinômios que proporcionam distribuição de tensões auto-equilibradas ou mesmo funções especiais relacionadas às soluções dos problemas de fratura. Uma extensão do teste numérico abordado em Zienkiewicz et al. (1986) é proposta como investigação inicial das condições necessárias para garantia de estabilidade da resposta numérica. O estudo da estabilidade é completado com a análise da condição de Babuška-Brezzi (inf-sup). Esta condição é aplicada nos elementos finitos quadrilaterais híbridos-mistos enriquecidos por meio de um teste numérico, denominado de inf-sup teste, desenvolvido com base no trabalho de Chapelle e Bathe (1993). Exemplos numéricos revelam que a FHMT é uma interessante alternativa para obtenção de boas estimativas para os campos de tensões e deslocamentos, usando-se enriquecimento sobre alguns nós de malhas pouco refinadas / This work presents a combination of hybrid-mixed stress model formulation (HMSMF) (Freitas et al. (1996)), to treat plane elasticity problems, with generalized finite element method (GFEM), (Duarte et al. (2000)). GFEM is characterized as a nonconventional formulation of the finite element method (FEM). GFEM is the result of the incorporation of concepts and techniques from meshless methods. One example of these techniques is the nodal enrichment that was formulated in the “hp” clouds method. Since two fields in domain (stress and displacement) and one in boundary (displacement) are approximated in the HMSMF, different possibilities of nodal enrichment are tested. For the discretization of the hybrid-mixed model quadrilateral finite elements with bilinear shape functions for the domain and linear elements for the boundary were employed. These functions are enriched with polynomial functions, trigonometric functions, polynomials that generate self-equilibrated stress distribution, or, even special functions connected with solutions of fracture problems. An extension of the numerical test cited in Zienkiewicz et al. (1986) is proposed as initial investigation of necessary conditions to assure the stability of the numerical answer. The stability study is completed with the analysis of the Babuška-Brezzi (inf-sup) condition. This last condition is applied to hybrid-mixed enrichment quadrilaterals finite elements by means of a numerical test, denominated inf-sup test, which was developed based on paper of Chapelle and Bathe (1993). Numerical examples reveal that HMSMF is an interesting alternative to obtain good estimates of the stress and displacement fields, using enrichment over some nodes of poor meshes

O uso de antioxidantes de origem natural apresenta uma perspectiva promissora na prevenção dos efeitos adversos causados pela radiação ultravioleta, incluindo o fotoenvelhecimento e a fotocarcinogênese. As atividades antioxidante, fotoprotetora e antienvelhecimento cutâneo foram demonstradas para o extrato de Pothomorphe umbellata. No entanto, seu potencial protetor da fotocarcinogênese ainda não foi estabelecido. O objetivo deste estudo é avaliar os efeitos da aplicação tópica do gel de extrato de raiz de P. umbelallta nas alterações promovidas pelas exposições aguda (dose única - 0,20 KJ/m2) e crônica (4 semanas, dose total de 8,6 KJ/m2, 5 semanas, 13,17KJ/m2, 15 semanas, 55,51 KJ/m2) à radiação ultravioleta B na pele de camundongos sem pelo. Para avaliação dos antioxidantes não enzimáticos, ácido ascórbico e α-tocoferol, foram utilizados métodos de cromatografia líquida de alta eficiência e para avaliação dos enzimáticos, catalase, superóxido dismutase, glutationa peroxidase e glutationa redutase, métodos espectrofotométricos. As modulações das proteínas p53 e PCNA foram avaliadas por imunohistoquímica. Os resultados demonstram que o tratamento com gel de extrato de raiz de P. umbellata é capaz de modular o sistema antioxidante na pele de camundongos sem pêlo, proporcionando maior proteção frente à exposição à radiação UVB. Além disso, os resultados sugerem que a aplicação tópica de extrato de P. umbellata protege contra os danos causados pela radiação modulando vias de reparo e/ou morte celular. Estes resultados contribuem para elucidação do mecanismo de ação do extrato de P.umbellata na pele exposta à radiação UVB, indicando sua capacidade em preservar a pele sadia após a irradiação. / Topical application of botanical antioxidants has been appreciated as a viable approach to reduce the occurrence of deleterious damage cased by sunlight exposure, such as photoaging and photocarcinogenesis. The plant Pothomorphe umbellata has been shown to exert a potent antioxidant activity on the skin and to delay the onset and reduce the incidence of UVB-induced skin damage and photoaging. The purpose of this study was to evaluate the effects of topical P. umbellata treatment on acute (single dose 0.20 KJ/m2) and chronic (4 weeks, total dose 8.6 KJ/m2, 5 weeks, 13.17KJ/m2, 15 weeks, 55.51 KJ/m2) UVB-induced damage hairless mouse epidermis. Nonenzymic antioxidants, ascorbic acid and α-tocopherol, were assayed by high performance liquid chromatography and enzymic antioxidants, catalase, superoxide dismutase, glutathione peroxidase and glutathione reductase, spectrophotometrically. UVB-induced p53 and PCNA were deteted by immunohistochemical technique. The data collected here suggest that topical application of P. umbellata extract affords strong protection against UVB-induced damages in skin epidermis by either preventing DNA damage or enhancing repair and inhibiting UVB-caused apoptosis and sunburn cell formation. Moreover, the results also suggest that P. umbellata extract modulates the cutaneous antioxidant network resulting in reduced susceptibility to UVB-induced oxidative damage. Together these results contribute for the elucidation of the mechanism of action of P.umbellata on UVB exposed skin and indicate the capacity of this extract to preserve skin health after irradiation.

The use of computational systems to solve engineering problems is increasingly present practical in the design offices. The Finite Element Method is one of the numerical techniques used in many of these systems, like the DOOLINES, a framework developed in objectoriented language, that allows a time domain dynamic analysis of the mooring lines and risers. In DOOLINES, until this work, the mooring lines and risers were discretized in onedimensional finite elements using a simplified formulation with six degrees of freedom (translational) element, with the consideration of only the axial and bending stiffness, when desirable, being this indirectly. This work incorporates into the DOOLINES class structure a new complete formulation element type that uses the co-rotational theory, with twelve degrees of freedom (translational and rotational) element, able to represent the axial, bending and torsional stiffness. The lumped-mass matrix is used in both formulations. This paper makes a comparative study of these elements about the quality of the response and the time required for simulation. / O uso de sistemas computacionais para resolução de problemas de engenharia é prática cada vez mais presente nos escritórios de projeto. O Método dos Elementos Finitos é uma das técnicas numéricas utilizadas em muitos desses sistemas, a exemplo do DOOLINES, um framework desenvolvido em linguagem orientada a objetos, que possibilita a análise dinâmica, no domínio do tempo, de linhas de ancoragem e risers. No DOOLINES, até então, as linhas de ancoragem e os risers eram discretizados em elementos finitos unidimensionais que utilizam uma formulação simplificada com seis graus de liberdade (translacionais) por elemento, com a consideração apenas das rigidezes axial e à flexão, quando desejável, sendo essa de maneira indireta. Neste trabalho, incorpora-se à estrutura de classes do DOOLINES um novo tipo de elemento com formulação completa de viga que utiliza a teoria corrotacional, com doze graus de liberdade (translacionais e rotacionais) por elemento, capaz de representar diretamente as rigidezes axial, à flexão e à torção. Em ambas as formulações considera-se a matriz de massa concentrada, onde se faz um estudo comparativo desses elementos no que se referem à qualidade da resposta e ao tempo de processamento necessário para simulação.

O cloridrato de tramadol, analgésico sintético de ação central, possui dois centros quirais: o isômero cis é ativo e o isômero trans é uma impureza de processo. Ambos os enantiômeros do cloridrato de cis-tramadol contribuem para o efeito analgésico, mas através de mecanismos diferentes, complementares e interativos farmacologicamente. Os dois isômeros do cis-tramadol apresentam efeitos terapêuticos, e a presença de impurezas, incluindo os isômeros trans - decorrentes do processo de síntese ou devido à decomposição - podem comprometer a qualidade do produto comercializado. Assim, este trabalho teve como objetivo desenvolver e validar metodologia estereosseletiva para análise do cloridrato de cis-tramadol e das possíveis impurezas quirais ou não na matéria-prima e formulações farmacêuticas. Para a separação e quantificação dos enantiômeros do cloridrato de cis-tramadol e das impurezas trans-tramadol, 1,2-olefina e 1,6-olefina, foi utilizada a coluna Chiralcel® OD-H, fase móvel constituída por hexano (60% e 100% de n-hexano, 1:1, v/v):isopropanol:dietilamina:ácido trifluoracético (99,5:0,5:0,3:0,1, v/v/v/v), na vazão de 0,7 mL min-1 e detecção em 274 nm. A coluna Chiralpak® AD fase móvel constituída por hexano (60% de n-hexano):etanol absoluto:dietilamina (95:5:0,1, v/v/v), na vazão de 1,0 mL min-1 e o comprimento de onda para detecção dos compostos foi de 228 nm foi utilizada para a separação e quantificação das impurezas O-desmetiltramadol, N-desmetiltramadol e tramadol N-óxido. Os métodos desenvolvidos foram devidamente validados através dos parâmetros seletividade, linearidade, precisão, exatidão, intervalo, limite de detecção e limite de quantificação. Os resultados obtidos na validação mostraram que os métodos são adequados para a determinação do cis-tramadol e de suas impurezas na matéria prima e na formulação farmacêutica. / Tramadol hydrochloride is a centrally acting analgesic with two chiral centers; the cis isomer is the active drug and the trans isomer is a process impurity. Both enantiomers of cis-tramadol hydrochloride contribute to the analgesic effect through different, but complementary and interactive pharmacological mechanisms. Although both isomers of cis-tramadol hydrochloride show therapeutic effects, the presence of impurities, originated from the synthesis process or due to degradation, can compromise the quality of the marketed product. The aim of this present work was the development and validation of a stereosselective methodology for the analysis of the drug cis-tramadol hydrochloride and the possible chiral or non-chiral impurities in raw materials and pharmaceutical formulations. The separation and quantitation of cis-tramadol enantiomers and the impurities trans-tramadol, 1,2-olefin and 1,6-olefin were carried out using a Chiralcel® OD-H column, mobile phase of hexane (60% and 100% of n-hexane, 1:1, v/v):2-propanol:diethylamine:trifluoroacetic acid (99,5:0,5:0,3:0,1, v/v/v/v) at a flow rate of 0,7 mL min-1 and detection at 274 nm. For the separation and quantitation of the impurities O-desmethyltramadol, N-desmethyltramadol and tramadol N-oxide, a Chiralpak® AD column was used with a mobile phase of hexane (60% of n-hexane):ethanol absolute: diethylamine (95:5:0,1, v/v/v) at a flow rate of 1,0 mL min-1 and detection at 228 nm. The methods were validated using the parameters selectivity, linearity, precision, accuracy, range, detection limit and quantitation limit. The results obtained show that the methods are suitable for the analysis of cis-tramadol and its impurities in raw material and pharmaceutical formulation.