Made available in DSpace on 2015-03-03T11:52:53Z (GMT). No. of bitstreams: 0 Previous issue date: 2014-07-25Bitstream added on 2015-03-03T12:06:50Z : No. of bitstreams: 1 000807287.pdf: 653649 bytes, checksum: 86f360131853a430d43838283923d6a1 (MD5) / Os objetivos deste estudo foram construir uma plataforma de força capaz de simular várias alturas de salto do calçado, analisar a distribuição de percentual de peso aplicado em antepés e retropés e a angulação da coluna lombar e báscula pélvica em três alturas de salto, sendo elas: 1,70 cm, 5,20 cm e 7,92 cm. Participaram deste estudo 20 voluntárias com média de idade de 20,4 (2,7), IMC 21,4 (2,1) e número do calçado 36,3 (1,2). Segundo os resultados observados, quanto maior a altura do salto utilizado, maior é também o percentual de peso aplicado sobre antepés e menor é o percentual de peso aplicado sobre os retropés. Na altura de salto de 1,70 cm, a média obtida do percentual de peso aplicado em antepés não diferiu da média encontrada por Cavanagh (1987). Concluiu-se, então, que com altura de salto de 1,70 cm não ocorre sobrecarga na região plantar de antepés. Na análise da angulação da coluna lombar, encontrou-se que o aumento da altura do salto está associado ao aumento da lordose lombar, de forma que os valores da angulação lombar mostraram-se diferentes e menores com o aumento da altura do salto. Na análise da angulação da báscula pélvica não foram encontradas alterações significativas / The objectives of this study were to build a force platform able to simulate various heel heights of shoes, to analyze the distribution of percentage of weight applied on forefeet and hindfeet and the angulation of the lumbar spine and pelvis in three heel heights, namely: 1.70 cm, 5.20 cm and 7.92 cm. The study included 20 volunteers with a mean age of 20.4 (2.7), BMI 21.4 (2.1) and shoe size 7.5 (1.2). According to the observed results, the higher the heel height used, the largest the percentage of weight applied on forefeet and the lower the percentage of weight applied on hindfeet. In heel height of 1.70 cm, the mean obtained from the percentage of weight applied on forefeet did not differ from the mean found by Cavanagh (1987). Then, it was concluded that there is no overload on the plantar region of forefeet with heel height of 1.70 cm. In the analysis of the angulation of the lumbar spine, it was found that the increase in height of the heel is associated with an increase in lumbar lordosis, so that the lumbar angle values were different and smaller by increasing the heel height. In the analysis of the angulation of the pelvis no significant changes were found

Orientador: Roger Augusto Rodrigues / Coorientador: Adriano Souza / Resumo: O comportamento dos solos não saturados é um tema amplamente estudado, tanto por meio de ensaios de laboratório quanto de campo. Fenômenos como o colapso dos solos são relativamente bem explicados pelos modelos elastoplásticos existentes, dos quais o Modelo Básico de Barcelona é um dos mais utilizados e difundidos no meio acadêmico. Apesar dos avanços obtidos, ainda existe certa escassez de estudos que comparem o comportamento previsto por tais modelos a resultados de situações reais de engenharia. Esta pesquisa teve como objetivo simular numericamente resultados de recalques devidos a carregamentos e inundação de provas de carga em placa e de protótipos de sapata construídas sobre solo natural e compactado em cava. Com os resultados buscou-se analisar a eficiência da técnica de compactação para mitigar o problema do colapso do solo. Para isso, ensaios de compressão triaxial e edométricos com sucção controlada, bem como ensaios de curva de retenção e condutividade hidráulica foram realizados para determinação dos parâmetros constitutivos hidráulicos e mecânicos do solo natural e compactado que foram utilizados no programa Code_Bright para simulação numérica dos resultados de campo. Os resultados demonstraram a capacidade do modelo em reproduzir o comportamento do solo no campo. Demonstraram também que o uso da compactação do solo que preconiza compactar em cava uma camada de espessura equivalente à largura da sapata é viável para redução da colapsibilidade, mas depende do nív... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: The behaviour of unsaturated soils is a widely studied theme around the world, both by utilising laboratory and field tests. Phenomena like soil collapsibility are relatively well explained by the existing elastoplastic models, among which the Barcelona Basic Model is the most widely used in the academic field. Despite the recent advances, there is still a scarcity of papers that compare the behaviour predicted by constitutive models and results of real engineering problems. This research aims to numerically simulate the results of settlement due to loading and wetting of load tests built both on undisturbed and compacted soil. The results were analysed in order to study the efficacy of the soil compaction technique in mitigating soil collapse. Suction controlled triaxial compression and oedometer tests, as well as soil water retention curve and hydraulic conductivity experiments were conducted to determine the hydraulic and mechanic parameters used in the software Code_Bright to numerically simulate field results. The results highlight the capability of the model in reproducing the analysed field behaviour. They also show that the technique that involves compaction of the soil equal to the footing's width is viable to reduce its collapsibility, but its efficacy will also depend on the tensions applied on the foundation and on the tolerable settlement values. Tensions that are propagated to the natural soil layer that are higher than its preconsolidation stress will cause col... (Complete abstract click electronic access below) / Mestre

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico. / Made available in DSpace on 2012-10-18T00:09:59Z (GMT). No. of bitstreams: 0Bitstream added on 2014-09-25T19:32:56Z : No. of bitstreams: 1 176596.pdf: 3279357 bytes, checksum: 6e094f8a64844fb30049288f5944b9d8 (MD5) / A simulação do processo de extração de óleo em um reservatório de petróleo é um problema complexo que envolve o escoamento multifásico em um meio poroso de geometria arbitrária, podendo conter falhas geológicas e outras irregularidades. A modelagem numérica implícita deste problema requer a solução de um complexo sistema de equações não-lineares que pode ser resolvido de diversas formas, mais ou menos eficientes. A forma mais comum de solução deste sistema emprega o método de linearização de Newton ou alguma de suas variantes. O foco deste trabalho está na análise de uma extensa classe de métodos de Newton conhecidos como métodos de Newton inexatos. Estes métodos empregam um segundo método iterativo na solução da equação linearizada resolvendo-a até uma certa tolerância; daí serem chamados de inexatos. Os métodos inexatos abordados neste trabalho concentram-se na escolha de processos iterativos baseados em métodos de projeção no espaço de Krylov, tais como: GMRES, BICGSTAB, TFQMR. Estes métodos, que formam uma subclasse dos métodos de Newton inexatos, são conhecidos como métodos de Krylov-Newton. São abordados também métodos quase-newtonianos baseados na aproximação da matriz Jacobiana considerando apenas os termos ortogonais da discretização em sistemas curvilíneos não-ortogonais. Vários problemas com diferentes graus de não-ortogonalidade da malha e não-linearidades são resolvidos e comparados usando os diferentes métodos implementados.

Dissertação (Mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico. / Made available in DSpace on 2012-10-17T07:54:44Z (GMT). No. of bitstreams: 0Bitstream added on 2016-01-08T23:05:33Z : No. of bitstreams: 1 145052.pdf: 49228177 bytes, checksum: 92aaa98b4beff1698dab154a25bf5ff5 (MD5)

Dissertação (mestrado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Mecânica, Florianópolis, 2015. / Made available in DSpace on 2015-11-17T03:05:26Z (GMT). No. of bitstreams: 1 336029.pdf: 5635838 bytes, checksum: d12d5cb0c2293c0248e67f7011d9d368 (MD5) Previous issue date: 2015 / Com o contínuo desenvolvimento e competitividade da indústria automotiva, tem-se a crescente necessidade de melhorias, tanto de durabilidade quanto de desempenho dos componentes e sistemas, sendo acompanhado da constante busca de menores prazos e custos possíveis, o que gera intenso desenvolvimento de modelos para testes. Neste trabalho são desenvolvidos modelos de dinâmica veicular através de modelagem dinâmica (1D) com o uso da metodologia fluxo de potência, no ambiente computacional AMESim. Essa metodologia é largamente utilizada para simulação de sistemas de várias áreas da engenharia, como hidráulica e pneumática, mas ainda é pouco utilizada na análise de veículos quando comparada com modelagens computacionais por Elementos Finitos (FEA) ou Multi Corpos, apesar de sua fundamentação teórica. No trabalho é realizada a construção de modelos de subsistemas de veículos da categoria universitária Baja SAE, utilizados como estudo de caso. Mais especificamente são estudados os sistemas de transmissão e suspensão, no intuito de representar situações reais de falhas e de desempenho adequado, ocorridas nos veículos da Equipe UFSC Baja SAE. A pesquisa também tem o objetivo de demonstrar os benefícios da modelagem 1D para simulação de sistemas automotivos, destacando essa alternativa como uma ferramenta pra aprimorar o projeto desde suas fases iniciais, fase de projeto preliminar do processo de desenvolvimento de produtos (PDP), diminuindo custos e tempo do projeto. Para contextualizar o trabalho, são apresentadas metodologias de desenvolvimento de produto, da área academia e do setor automotivo; a metodologia fluxo de potência é apresentada; seguida de uma descrição dos aspectos gerais da dinâmica veicular e de alguns dos subsistemas que compõe um automóvel; e por fim são apresentados os modelos propostos, juntamente com as hipóteses matemáticas, simplificações adotadas e resultados obtidos.<br> / Abstract : With increasing development and competitiveness of the automotive industry, there is a growing need for improvements both in durability and performance of components and systems, followed by a constant search for shorter terms and costs, which generates intense development of models for tests. In this work, dynamic models are developed by 1D modeling using power flow methodology in AMESim computing environment. This methodology is widely used for systems simulation in several areas of engineering, such as hydraulics and pneumatics, but despite its theoretical foundation, it has been less used in vehicles analysis compared to other computational modeling, such as finite elements (FEM) or Multibody Systems. This research presents the development of vehicle subsystems models, using Baja SAE competition vehicles as study case. Specifically transmission and suspension systems are studied in order to represent real situations of failures and adequate performance that occurred in vehicles of Team UFSC Baja SAE. The research also aims to demonstrate the benefits of 1D modeling for automotive systems simulation, highlighting it as one tool to enhance the design process since its early stages, preliminary design phase of the product development process (PDP), reducing costs and project time. To contextualize the work, product development methodologies are presented, both from academic area and automotive sector; the power flow methodology is presented followed by a description of the general aspects of vehicle dynamics and some of the subsystems that compose a car and finally the proposed models are presented, along with the mathematical assumptions, simplifications adopted and results.

Tese (doutorado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Química, Florianópolis, 2014. / Made available in DSpace on 2015-04-29T21:03:07Z (GMT). No. of bitstreams: 1 333228.pdf: 10004087 bytes, checksum: 81bb960774686a26cd25ae783ad8e01c (MD5) Previous issue date: 2014 / Nesta tese, através da aplicação do Método da Média no Volume, desenvolveu-se o modelo matemático de duas equações para descrever o processo de separação e isomerização do p-xileno, visando contribuir na predição de melhores condições de operação de coluna de leito fixo. Este modelo também foi aplicado na separação de compostos BTX. O modelo matemático permitiu carregar hierarquicamente as informações físicas que descrevem os mecanismos de transferência de massa entre as escalas de comprimento (microescala até a macroescala). O modelo matemático de duas equações foi obtido para coluna de leito fixo empacotada com pellet catalisador usando o mecanismo de reação triangular para isomerização do p-xileno, considerando condição de salto na interface da partícula. Também foi obtido o modelo matemático de duas equações para coluna de leito fixo empacotada com partículas de adsorvente utilizando a isoterma linear, a isoterma de Langmuir e a isoterma competitiva de Langmuir, considerando condição de salto na interface da partícula. Foi aplicado o Método de Volumes Finitos na discretização das equações dos Problemas de Fechamento da microescala e macroescala e equações governantes de transporte de massa na escala de Darcy. As soluções foram encontradas numericamente através de um código computacional e alguns pacotes do software livre OpenFOAM®, versão 2.2.x. Foi possível resolver numericamente os Problemas de Fechamento originados da modelagem da microescala do pellet catalisador e da partícula do adsorvente para encontrar o tensor de difusividade efetiva da microescala sobre arranjos bidimensionais de cilindros e tridimensionais de esferas. Foram determinados numericamente os tensores de difusividade efetiva, de dispersão total, convectivo de transferência de massa, e outros tensores cruzados, sobre os arranjos bidimensionais de cilindros em linha e tridimensionais de esferas em bloco hexagonal. Resolveram-se numericamente as equações governantes de transporte de massa para a escala de Darcy (escala de projeto) usando os tensores de transporte calculados numericamente pelos problemas de fechamento e corroborou-se com dados da literatura. Foi realizada uma análise da influência de cada tensor de transporte no modelo de duas equações e uma análise entre as correlações de Wakao e Funazkri (1978) e Wilson e Geankoplis (1966), aplicadas na separação dos compostos BTX. O modelo de duas equações apresentado neste trabalho é o modelo mais completo, considerando que contém todos os tensores de transportes calculados teoricamente.<br> / Abstract : In this thesis, by applying the Method of Volume Averaging, the twoequation mathematical model is developed to describe the process of separation and isomerization of p-xylene in order to contribute better to predict fixed bed column operating conditions. This model was also applied in the separation of BTX compounds. The two-equation mathematical model allows loanding the physical hierarchical information describing the mass transfer mechanisms between the length scales (microscale until macroscale). The two-equation mathematical models wereobtainedto fixed bed column packed with catalyst pellets using the triangular reaction mechanism for p-xylene isomerization considering jump condition in the particle interface. It was also obtained the mathematical model of two equations for fixed bed column packed with adsorbent particles using linear isotherm, the Langmuir isotherm and the competitive Langmuir isotherm, considering jump condition in the particle interface. The equations from closure problems and Darcy s scale of this work were discretized using the finite volumes method. The solutions are found numerically through of a computational code and some packages from the free software OpenFOAM®, version 2.2.x. It was possible to solve numerically the closure problems, originated from modeling of microscale catalyst pellet and adsorbent particle to find the effective diffusivity tensor of microscale on arrangements two-dimensional cylinders and threedimensional of spheres. It was determined numerically the effective diffusivity tensor, the total dispersion tensor, the convective mass transfer tensor, and other cross tensors on arrangements twodimensional cylinders in-line and three-dimensional in hexagonal block. It was solved numerically the mass transport governing equations for the Darcy scale (project scale) using the transport tensor calculated numerically by solving the closure problems and corroborated with literature data. An analysis of the influence of each transport tensor was performed in the two-equation mathematical model and analysis between the correlations Wakao and Funazkri (1978) and Wilson and Geankoplis (1966), applied in the separation of BTX compounds. The two-equation mathematical model presented in this work is the most complete model, considering that contains all transport tensors theoretically calculated.

Tese (doutorado) - Universidade Federal de Santa Catarina, Centro Tecnológico, Programa de Pós-Graduação em Engenharia Elétrica, Florianópolis, 2015. / Made available in DSpace on 2015-12-01T03:15:02Z (GMT). No. of bitstreams: 1 336271.pdf: 7722227 bytes, checksum: c33c8ca3e76b93528dfdd5114e97b122 (MD5) Previous issue date: 2015 / A qualidade dos modelos tem sido uma constante preocupação dos engenheiros de planejamento, desde que a tarefa atribuída é o de predizer o desempenho do sistema de potência por meio de simulações digitais. Baseados nestes estudos são estabelecidos limites de transferências e políticas operativas que buscam preservar a segurança do sistema. O nível de confiança das simulações, como resposta esperada às perturbações, é principalmente influenciado pela qualidade dos modelos das unidades geradoras. Neste trabalho, é proposta uma nova metodologia para a validação de modelos de unidades geradoras utilizando registros obtidos de unidades de medição sincrofasorial. Das três fases propostas para o procedimento de validação, contribuições importantes são feitas na verificação e calibração. Na verificação é proposta uma nova função de discrepância global GDI (Global Discrepancy Indicator) para quantificar objetivamente a qualidade do modelo. A identificação dos parâmetros causadores das discrepâncias, factíveis de serem ajustados, é fundamental para o sucesso da calibração, devido à grande quantidade de parâmetros e as poucas variáveis registradas. Neste trabalho é mostrada a necessidade de se estabelecer, primeiro, quais parâmetros são observáveis; segundo, quais são os problemáticos e, terceiro, quais são determináveis de maneira única. Para cada um destes passos são propostas estratégias novas, baseadas nas propriedades da função GDI e nas sensibilidades das trajetórias do sistema. Para melhorar a observabilidade e validar cada um dos modelos componentes da unidade geradora é proposto o projeto de um tipo especial de unidade de medição sincrofasorial, com capacidade de registrar medidas fasoriais e medidas não fasoriais sincronizadas. Assim, é desenvolvida uma extensão da metodologia proposta envolvendo a ampliação do escopo de aplicação da simulação dinâmica híbrida, bem como da função GDI e de suas propriedades. A metodologia é validada realizando-se experimentos computacionais com registros virtuais e registros reais de uma unidade geradora de 270 MVA. Os resultados obtidos evidenciam as vantagens das propostas feitas e são considerados satisfatórios, aplicáveis e reproduzíveis para qualquer outro subsistema.<br> / Abstract : The quality of models has been a constant concern of planning engineers, since their main goal is to predict the performance of electrical power system through digital simulations. Based on these studies both transfer limits and operational policies are established in order to preserve system security. The confidence level of the simulations, measured as the system's fit to expected disturbances, is mainly influenced by the quality of the models that represent the generating units. This initiative proposes a new methodology to validate the model of generating units by using digital records sampled in synchrophasor measurement units. From the three stages proposed for the validation procedure, important contributions are made in the verification and calibration phases. In the verification, a new discrepancy function GDI (Global Discrepancy Indicator) is proposed to quantify in an objective way the quality of the model. The identification of the causative parameters of feasible discrepancies to be adjusted, define the final success of calibration, due to the large set of parameters and the few recorded variables. This work shows the need to establish, first, what parameters are observable; second, which ones are problematic, and third, which ones are determinable in a unique way. For each of these steps new strategies are proposed, based on GDI characteristics and trajectory sensitivities. To improve the observability and validate each model of the generator unit is proposed the design a special type of synchrophasor measurement unit with the ability to record synchronized phasor measurements and not phasor measurements. Additionally, it is proposed to extend the methodology in conjunction with the extended dynamic hybrid simulation and extended GDI functions. The methodology is verified by performing computational experiments with virtual records and real records of a 270 MVA generating unit. The results show the advantages of the proposals made and are considered satisfactory, applicable and reproducible for any other subsystem. Furthermore, when the special synchrophasor measurement unit proposed is used, the calibration phase performance is higher, due to the increase of observability

Tese (doutorado) - Universidade Federal de Santa Catarina, Centro de Ciências Físicas e Matemáticas. Programa de Pós-Graduação em Física. / Made available in DSpace on 2012-10-22T04:44:54Z (GMT). No. of bitstreams: 1 204258.pdf: 1972812 bytes, checksum: f2133f2deada68c59f0a45cd8e721e5a (MD5) / Neste trabalho estudamos alguns modelos estatísticos para a água e para soluções de moléculas anfifílicas. Obtivemos as propriedades termodinâmicas do modelo da água quadrada em um campo elétrico externo a partir de cálculos analíticos exatos e de soluções de campo médio, bem como através de simulações de Monte Carlo no ensemble canônico. Observamos o aparecimento de uma transição de fases induzida pelo campo quando a temperatura do sistema tende a zero, e verificamos um máximo no número de pontes de hidrogênio para um valor não nulo do campo externo. Em um outro estudo realizamos cálculos de campo médio em um modelo geométrico para a água em duas dimensões, obtendo o seu diagrama de fases no plano pressão versus temperatura, e a dependência da densidade e do número de pontes de hidrogênio com a temperatura e a pressão do sistema. Também para a água, desenvolvemos um modelo de trímeros interagentes em uma rede triangular, e estimamos suas propriedades termodinâmicas a partir de simulações de Monte Carlo no ensemble multicanônico. Essas simulações forneceram medidas diretas da entropia configuracional do sistema, bem como facilitaram a identificação de transições de fases de primeira ordem entre estados de diferentes densidades. Investigamos também modelos de agregação com propriedades semelhantes àquelas observadas em sistemas micelares. Dois modelos do tipo gás de rede, um com n estados e interações de pares, e outro com interações de pares e tripletos, foram explorados usando-se aproximações de campo médio usual e simulações computacionais. Verificamos que, apesar de sua simplicidade, os modelos apresentam comportamento similar ao de soluções de anfifílicas: a presença de uma concentração micelar crítica, e os mínimo e máximo locais na curva de distribuição de agregados. Por fim, analisamos alguns aspectos dinâmicos da agregação micelar, através da determinação dos tempos característicos de entrada e saída de anfifílicas de agregados (t1), e de formação/dissolução micelar (t2).

Made available in DSpace on 2016-01-13T13:26:48Z (GMT). No. of bitstreams: 0 Previous issue date: 2015-09-02. Added 1 bitstream(s) on 2016-01-13T13:32:06Z : No. of bitstreams: 1 000855350.pdf: 5394081 bytes, checksum: 5f4d040bbbe3297c23997fc77559961e (MD5) / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Os fornos contínuos do tipo túnel, dada a sua maior flexibilidade no aquecimento e alta produtividade, tornaram-se a melhor opção para o processamento de produtos alimentícios industrializados. O combustível utilizado nos queimadores é normalmente o gás natural, que apesar das incertezas e oscilações dos custos, vem ganhando espaço na indústria nacional. No presente trabalho foram desenvolvidos diversos modelos físicos tridimensionais - classificados em modelos de testes, de análises paramétricas e de otimização - e utilizados no estudo, modelagem e otimização de fornos alimentícios contínuos do tipo túnel com aquecimento indireto. Utilizou-se o software livre de código aberto OpenFOAM® para a modelagem dos problemas propostos. A geometria e condições operacionais utilizadas foram obtidas junto a uma empresa alimentícia bauruense para o processo de pré-assamento de pizzas. Os efeitos da geometria, malha, condições iniciais e de contorno, algoritmo de acoplamento pressão-velocidade, forças de campo, trocas de calor por radiação através dos modelos fvDOM e view/Factor e modelagem RANS k- E para a turbulência foram avaliados sobre os parâmetros do escoamento. O algoritmo de acoplamento pressão-velocidade mais indicado ao problema proposto é o SIMPLE. A utiilzação do modelo de radição viewFactor mostrou-se mais adequado. As trocas de calor por convecção sobre a esteira mostraram-se pouco significantes em relação as trocas por radiação, demonstrando a importância da correta modelagem desse fenômeno. Testes para diversas geometrias de forno foram realizados. O forno'com menor razão de aspecto geométrico e consequentemente, de menor altura, mostrou-se mais eficiente, embora isto possa provocar, na prática, limitações na sua utilização / Because of the better flexibility in warning and high production, the continuous furnace tunnel is the better option to the processing of industrialized food products. The fuel used in the burners is usually natural gas, which despite the uncertainties and fluctuantions in costs, is becoming more popular in national industry. In this work, several three-dimensional physical models were developed - classified in testing models, parametric analysis and optimization - and used to the study, modeling and optimization of continuous food furnaces of tunnel type with indirect heating. It was chosen the OpenFOAM® opensource software for modeling of proposed problems. The geometry and the operating conditions were based in a bauruense food company to the process of pre-baking of pizzas. The effects of geometry, mesh, initial and boundary conditions, pressure-velocity coupling algorithm, field strengths, radiation heat exchanges using fvDOM and view/Factor, and modeling RANS k - E for turbulence to were evaluated about the flow parameters. The SIMPLE is the most indicated pressure-velocity coupling algorithm to this problem. The use of viewFactor radiation model proved more adequate. The heat exchanges by convection over the treadmill have proved little significant compared the heat exchanges by radiation, demonstrating the importance of correct modeling of this phenomenon. Differents geometries tests were made, and the furnace with the smaller ratio aspect geometric and, consequently, smaller height, showed to be most efficient, although this may result in practice, in limitations in its use

Made available in DSpace on 2014-12-02T11:16:54Z (GMT). No. of bitstreams: 0 Previous issue date: 2014-02-17Bitstream added on 2014-12-02T11:20:56Z : No. of bitstreams: 1 000792167.pdf: 2217927 bytes, checksum: dca6b4b2a633a970667631439071e874 (MD5) / Neste trabalho buscam-se formulações de modelos analíticos e d aelaboração de programas computacionais para realizar uma busca por órbitas estáveis em torno de satélites planetários que poderão contribuir no planejamento de missões espaciais a serem conduzidas para o esudo destes corpos celestes. O estudo leva em consideração órbitas de satélites artificiais em torno de satélites plenetários sob a influência da petubação de terceiro corpo (a atração gravitacional de Júpiter ) e das pertubações devidas à distribuição não uniforme de massa (J2 e J3) do corpo principal (central) ... (Resumo completo, clicar acesso eletrônico abaixo) / This work seeks to formulations of analytical models and the development of computer programs to perform a search for stable orbits around planetary satellites that maybe helpful in planning space missions to be conducted to study these celestial bodies. The study takes into account artificial satellite orbits around planetary satellites under the influence of the third body pertubation (gravitational attraction of Jupiter) and the pertubations due to non-uniform distribution of mass (J2 and J3) of the main (central) body ... (Complete abstract click eletronic access below)