[pt] A análise de fluxo em solos parcialmente saturados exige o conhecimento das propriedades hidráulicas do solo: a curva característica, q(y), e a curva de permeabilidade, k(y). Entre os problemas geotécnicos que envolvem fluxo em meios parcialmente saturados pode-se citar a instabilização de encostas, muitas vezes associada a variações de sucção, e o transporte de contaminantes, pois, freqüentemente, a fonte contaminante está acima do nível d`água. O principal objetivo deste trabalho é a implementação numérica de um programa, baseado num código de elementos finitos, para a realização da retroanálise de parâmetros de fluxo em solos parcialmente saturados a partir de ensaios de campo, e que forneça também informações sobre a qualidade da estimativa obtida. O método de retroanálise utilizado permite a determinação simultânea da curva característica e da curva de permeabilidade a partir de dados de fluxo transiente. O programa desenvolvido foi utilizado para determinar os parâmetros de fluxo de um ensaio realizado no campo experimental II da PUC - Rio. Com base na análise dos resultados obtidos neste ensaio e na análise de sensibilidade foi proposto um ensaio de configuração mais simples que pode fornecer resultados satisfatórios na determinação das propriedades hidráulicas de solos parcialmente saturados. / [en] Flow analysis in partially saturated soils requires the determination of soil hydraulic properties: its characteristic curve, q(y), and its permeability curve k (y). The instabilization of slopes, often associated with suction variations, and the transport of contaminants, where the source of contaminants is often above water leve, are some of the geotechnical problems involving flows in partially saturated media. The main objective of this work is the numerical implementation of a program, based on finite element coding, for the back analysis of flow parameters in partially saturated soils, based on field experiments, capable of providing information on the quality of the estimates obtained. The back analysis method used here allows the simultaneous determination of the characteristic and permeability curves from transient flow data. The program was used to determine flow parameters of a field experiment conducted in Experimental Field II (PUC- Rio). A simpler configuration test is proposed based on the analysis of the results obtained in this field experiment and on a sensitivity analysis. This test may provide satisfactory results in the determination of hydraulic properties of partially saturated soils. / [es] El análisis de flujo en suelos parcialmente saturados exige el conocimiento de las propiedades hidráulicas del suelo: la curva característica, q(y), y la curva de permeabilidade, k(y). Entre los problemas geotécnicos que consideran flujos en medios parcialmente saturados se puede citar la inestabilidad de pendientes, muchas veces asociada a la variaciones de succión; y el transporte de contaminantes, ya que frecuentemente la fuente contaminante está encima del nível de agua. El principal objetivo de este trabajo es la implementación numérica de un programa, basado en un código de elementos finitos, para la realización del retroanálisis de parámetros de flujo en suelos parcialmente saturados a partir de ensayos de campo, y que brinde también informaciones sobre la calidad de la estimatición obtenida. El método de retroanálisis utilizado permite la determinación simultánea de la curva característica y de la curva de permeabilidad a partir de datos de flujo transitorio. El programa desarrollado fue utilizado para determinar los parámetros de flujo de un ensayo realizado en el campo experimental II de la PUC - Rio. Analizando los resultados obtenidos en este ensayo y el análisis de sensibilidad, se propone un ensayo de configuración más simple que puede dar resultados satisfactorios en la determinación de las propiedades hidráulicas de suelos parcialmente saturados.

Nesta tese são propostas duas distribuições discretas: COM-Poisson correlacionada (CPC) e COM-Poisson generalizada parcialmente correlacionada (CPGPC). Também foram propostos modelos de regressão para a distribuição Poisson generalizada parcialmente correlacionada (PGPC) (proposto por Luceño (1995)). Calculamos a função massa de probabilidade (fmp) para todas as distribuições com duas parametrizações. As distribuições foram construídas usando a mesma expansão feita por (Luceño, 1995) na construção da distribuição Poisson generalizada parcialmente correlacionada. A distribuição CPC(l;f;r) é a mesma expansão da distribuição COM-Poisson zero inflacionada ZICMP(m;f;r). Para a distribuição CPGPC(l;f;r;L;K) foi determinada a função característica, função geradora de probabilidade, momentos e a estimação pelo método de máxima verossimilhança para as duas parametrizações. Fizemos a fmp, quantil e gerador de números aleatórios das distribuições citadas no programa R. / In this thesis two discrete distributions are proposed: Correlated COM-Poisson (CPC) and Generalized partially correlated COM-Poisson (CPGPC). We have also proposed regression models for the Generalized partially correlated Poisson distribution (PGPC) (proposed by Luceño (1995)). We calculated the probability mass function for all distributions with two parametrizations. The distributions were constructed using the same expansion made by Luceño (1995) in the construction of the correlated generalized Poisson distribution. The CPC(l;f;r) Correlated COM-Poisson distribution is the same expansion of the zero-inflated COM-Poisson distribution ZICMP(m;f;r). For the CPGPC(l;f;r;L;K) Generalized partially correlated COM-Poisson distribution, the characteristic function, probability-generating function, moments, and the maximum likelihood estimation for the two parametrizations were determined. We performed the probability mass function, quantile and random number generator of the distributions quoted in program R.

[pt] O presente trabalho tem como finalidade avaliar os diferentes modelos de mistura para o cálculo da combustão de reagentes pré- misturados utilizando a abordagem de Reator Parcialmente Misturado (PaSR). Os modelos de mistura considerados neste trabalho foram os modelos IEM estendido, Langevin e Langevin estendido. Investiga-se aqui o grau de mistura previsto por tais modelos e sua influência sobre as propriedades termoquímicas em um processo de combustão. A primeira parte deste trabalho consiste na apresentação e avaliação destes modelos de mistura, considerando-se um campo escalar inerte em presença de um campo turbulento homogêneo e isotrópico. Uma vez que estes modelos de mistura envolvem formulações do tipo estocástico, sua implementação foi realizada utilizando o método de Monte Carlo, mediante a utilização de esquemas numéricos adequados à resolução de equações diferenciais estocásticas. Assim, estuda-se a evolução da Função Densidade de Probabilidade (PDF) e das principais propriedades do campo escalar para cada modelo implementado. Os resultados obtidos também são comparados com simulação numérica direta e com resultados analáticos disponsáveis. Um ótimo acordo em termos qualitativos e quantitativos é obtido. A segunda parte deste trabalho utiliza estes modelos para o estudo numérico de um PaSR no qual são modelados os processos difusivos e reativos presentes durante a combustão. O PaSR é usado para avaliar a influência dos modelos de mistura nas propriedades termoquímicas da mistura em uma situação de combustão de tipo pré-misturada, que é modelada utilizando-se uma variável de progresso de uma reação. Os resultados obtidos com os diferentes modelos de mistura são comparados para diferentes regimes de funcionamento do PaSR, mostrando que, em situações de mistura rápida e reação intensa, os diferentes modelos apresentam resultados similares. Porém, nos casos de mistura lenta e reação moderada, discrepancias importantes são observadas entre os resultados dos modelos; as quais atingem até 65% para o valor médio da variável de progresso da reação. / [en] The present work evaluates several mixing models for the prediction of premixed combustion in a Partially Stirred Reactor (PaSR). The models considered in this work were the extended IEM, Langevin and extended Langevin models. The degree of mixing and its influence on the termochemical properties in a combustion process are investigated here. The first part of this work consists on the presentation and the assesment of these mixing models in which a single scalar field was considered in presence of a homogeneous and isotropic turbulent field. Since these mixing models involve stochastic terms, their implementation is performed by the Monte Carlo method using numerical schemes which solve the corresponding Stochastic Differential Equations (SDE). The evolution of the Probability Density Function (PDF) and the main properties for a single scalar field are studied for each mixing model. The numerical results are compared with Direct Numerical Simulation and available analytical results. Excellent qualitative and quantitative agreements are obtained. In the second part of this work, mixing models are used for numerical simulation of a PaSR where the diffusive and reactive processes occur. The PaSR is used to assess the mixing model influence on the termochemical properties of the mixture in a premixed combustion process, which is modeled using a reaction progress variable. The results obtained with the different mixing models are compared in several operating regimes of the PaSR, showing that when mixing is fast and reaction is intense, the different models lead to similar results. However, when mixing is slow and reaction is weak, important discrepancies are observed between the model results, which reach 65%, as far as the averaged reaction progress variable is concerned

O objetivo principal deste trabalho é apresentar um modelo de análise probabilística do comportamento ao longo do tempo de estruturas de concreto. Para isso, é realizada a análise probabilística em conjunto com análise numérica. A análise numérica é realizada através de um programa computacional, baseado no método dos elementos finitos, que considera o comportamento não-linear e dependente do tempo dos materiais, assim como o processo evolutivo da construção. Para determinar o efeito das incertezas dos parâmetros é realizada uma análise probabilística, utilizando o método de amostragem por hipercubo latino. O modelo apresentado pode ser aplicado para a análise probabilística do comportamento ao longo do tempo das estruturas de concreto em geral. No entanto, foi aplicado na análise das flechas diferidas de lajes pré-moldadas formadas por vigotas com armação treliçada. Com esta análise, pôde-se propor um coeficiente multiplicador das flechas imediatas para a avaliação das flechas diferidas no tempo. O coeficiente multiplicador obtido para as lajes pré-moldadas formadas por vigotas com armação treliçada pode alcançar valores muito superiores ao fator 'alfa'f recomendado pela NBR 6118 (2003) para o caso de vigas de concreto armado. / The goal of this research is to present a probabilistic analysis model of the long-term behavior of concrete structures. For that, the probabilistic analysis is carried out together with numerical analysis. The numerical analysis is carried out using a software based on the finite element method that takes into account the nonlinear and time dependent behavior of the materials, as well as the evolutionary construction process. A probabilistic analysis is carried out in order to determine the effects of the uncertainties of the parameters, using latin hypercube sampling method. The presented model can be applied in the probabilistic analysis of the long-term behavior of concrete structures in general. However, it was applied in the analysis of the long-term deflections of precast slabs made by joist with lattice reinforcement. By this analysis, a multiplier coefficient of immediate deflections to evaluate long-term deflections could be proposed. The multiplier coefficient obtained for precast slabs made by joist with lattice reinforcement can reach values much higher than the factor 'alfa'f recommended by NBR 6118 (2003) for the case of reinforced concrete beams.

Neste trabalho estudamos relações entre expoente de Lyapunov e continuidade absoluta da folheação central para difeomorfismos parcialmente hiperbólicos conservativos de \'T POT. 3\'. Sobre tal tema, provamos que tipicamente (\'C POT. 1\' aberto e \'C POT. 2\' denso) os difeomorfismos parcialmente hiperbólicos, conservativos de classe \'C POT. 2\' , do toro \'T POT. 3\', apresentam folheação central não absolutamente contínua. Desta maneira, respondemos positivamente uma pergunta proposta em [20]. Também neste trabalho, estudamos entropia topológica para Sistema Iterado de Funções. Neste contexto, damos uma nova demonstração para uma conjectura proposta em [14] e provada primeiramente em [15]. Apresentamos um método geométrico que nos permite calcular entropia para transformações de \'S POT. 1\', como em [15]. Além de disso o método apresentado se verifica para casos mais gerais, como por exemplo: transformações não comutativas / In this work we study relations between Lyapunov exponents, absolute continuity of center foliation for conservative partially hyperbolic diffeomorphisms of \'T POT. 3\'. About this theme, (on a \'C POT. 1\' open and \'C POT. 2\'dense set) of conservative partially hyperbolic \'C POT. 2\' diffeomorphisms of the 3-torus presents non absolutely continuous center foliation. So, we answer positively a question proposed in [20]. Also in this work, we study topological entropy for Iterated Functions Systems. In this setting, we give a proof for a conjecture proposed in [14] and firstly proved in [15]. We present a geometrical method that allows us to calcule the entropy for transformations of \'S POT. 1\', like in [15]. Furthermore this method holds for more general cases, for example: non commutative transformations

A análise de texturas dinâmicas tem se apresentado como uma área de pesquisa crescente e em potencial nos últimos anos em visão computacional. As texturas dinâmicas são sequências de imagens de textura (i.e. vídeo) que representam objetos dinâmicos. Exemplos de texturas dinâmicas são: evolução de colônia de bactérias, crescimento de tecidos do corpo humano, escada rolante em movimento, cachoeiras, fumaça, processo de corrosão de metal, entre outros. Apesar de existirem pesquisas relacionadas com o tema e de resultados promissores, a maioria dos métodos da literatura possui limitações. Além disso, em muitos casos as texturas dinâmicas são resultado de fenômenos complexos, tornando a tarefa de caracterização um desafio ainda maior. Esse cenário requer o desenvolvimento de um paradigma de métodos baseados em complexidade. A complexidade pode ser compreendida como uma medida de irregularidade das texturas dinâmicas, permitindo medir a estrutura dos pixels e quantificar os aspectos espaciais e temporais. Neste contexto, o objetivo deste mestrado é estudar e desenvolver métodos para caracterização de texturas dinâmicas baseado em metodologias de complexidade advindas da área de sistemas complexos. Em particular, duas metodologias já utilizadas em problemas de visão computacional são consideradas: redes complexas e caminhada determinística parcialmente auto-repulsiva. A partir dessas metodologias, três métodos de caracterização de texturas dinâmicas foram desenvolvidos: (i) baseado em difusão em redes - (ii) baseado em caminhada determinística parcialmente auto-repulsiva - (iii) baseado em redes geradas por caminhada determinística parcialmente auto-repulsiva. Os métodos desenvolvidos foram aplicados em problemas de nanotecnologia e tráfego de veículos, apresentando resultados potenciais e contribuindo para o desenvolvimento de ambas áreas. / Dynamic texture analysis has been an area of research increasing and in potential in recent years in computer vision. Dynamic textures are sequences of texture images (i.e. video) that represent dynamic objects. Examples of dynamic textures are: evolution of the colony of bacteria, growth of body tissues, moving escalator, waterfalls, smoke, process of metal corrosion, among others. Although there are researches related to the topic and promising results, most literature methods have limitations. Moreover, in many cases the dynamic textures are the result of complex phenomena, making a characterization task even more challenging. This scenario requires the development of a paradigm of methods based on complexity. The complexity can be understood as a measure of irregularity of the dynamic textures, allowing to measure the structure of the pixels and to quantify the spatial and temporal aspects. In this context, this masters aims to study and develop methods for the characterization of dynamic textures based on methodologies of complexity from the area of complex systems. In particular, two methodologies already used in computer vision problems are considered: complex networks and deterministic walk partially self-repulsive. Based on these methodologies, three methods of characterization of dynamic textures were developed: (i) based on diffusion in networks - (ii) based on deterministic walk partially self-repulsive - (iii) based on networks generated by deterministic walk partially self-repulsive. The developed methods were applied in problems of nanotechnology and vehicle traffic, presenting potencial results and contribuing to the development of both areas.

Os pilares mistos são elementos estruturais compostos por um perfil metálico e concreto trabalhando em conjunto. Entre estes destaca-se o pilar parcialmente revestido que apresenta vantagens como a possibilidade de pré-fabricação, redução no uso de formas e eficiência estrutural, porém seu processo executivo envolve uma trabalhosa tarefa que é a ancoragem das armaduras ao perfil metálico. Neste contexto, o presente trabalho tem por objetivo estudar o comportamento estrutural de pilares mistos parcialmente revestidos submetidos a compressão simples e flexo-compressão e avaliar a possibilidade de substituição da armadura convencional por alternativas que possibilitem uma execução mais simples. Para isto foi desenvolvido um programa experimental envolvendo 23 exemplares de pilares mistos fabricados com o perfil W 150 x 22,5 e três configurações de armadura: armadura convencional, telas de aço soldadas entre as mesas do perfil e concreto com adição de fibras. Além disso, foram realizadas simulações numéricas e um estudo paramétrico utilizando o pacote computacional FX + DIANA a fim de extrapolar os resultados experimentais. Observou-se que o comportamento estrutural de pilares submetidos a forças excêntricas depende do eixo de flexão, de modo que a ruptura ocorre de modo mais gradual quando o pilar é submetido a flexão em torno do eixo de maior inércia. Além disso, respostas similares foram obtidas para as três configurações de armadura avaliadas, sendo o valor de força máxima obtido para os exemplares com concreto com fibras ligeiramente inferior aos demais. A substituição da armadura convencional por telas de aço ou concreto com fibras de aço não altera significativamente os valores de força máxima nem o comportamento pós-pico, mostrando a viabilidade da solução proposta. / Composite columns are structural elements made with steel profiles and concrete working together. Among these, there is the partially encased composite columns that presents advantages, such as, the possibility of pre-fabrication, reduction in the use of formwork and structural efficiency, but its fabrication process involves a laborious task that is the anchoring of the reinforcements to the steel profiles. In this context, the present work has the objective of studying the structural behavior of partially encased composite columns subjected to axial compression and eccentric compression and to evaluate the possibility of replacing conventional reinforcement with alternatives that allow a simpler execution. For this, an experimental program was developed involving 23 specimens of composite columns made with W 150 x 22.5 shape and three reinforcement configurations: welded steel mesh and fiber-added concrete. In addition, numerical simulations and a parametric study were performed using the FX + DIANA computational package in order to extrapolate the experimental results. The results show the structural behavior of columns subjected to eccentric forces depends on the axis of flexure, so that the rupture occurs more gradually when subjected to bending about the major inertia axis. In addition, similar responses were obtained for the three reinforcement configurations evaluated, and the peak load obtained for the specimens with fiber concrete was slightly lower than the others. Therefore, it can be observed that the replacement of conventional reinforcement by steel meshes or fiber concrete does not significantly alter the peak load or post-peak behavior, showing potential to be considered as a new alternative.

O pilar misto parcialmente revestido é fruto da associação de um perfil metálico e concreto trabalhando em conjunto. Nestes pilares o perfil metálico, normalmente perfil I ou H, não é totalmente envolvido por concreto, sendo apenas a região entre as mesas preenchida. No Brasil, a utilização de pilares mistos parcialmente revestidos é ainda incipiente. Contudo a introdução do dimensionamento de elementos mistos na ABNT NBR 8800:2008 e as vantagens construtivas inerentes aos elementos mistos tem viabilizado a utilização destes elementos em obras pelo país afora, especialmente naquelas de grande porte. Neste contexto, o objetivo principal do trabalho é analisar o comportamento estrutural destes pilares e avaliar a possibilidade de substituição das armaduras tradicionais compostas por vergalhões por alternativas como telas de aço soldadas ou adição de fibras de aço ao concreto (taxa de 1,5% de fibras de 25 mm de comprimento). Para isto, foram feitos ensaios experimentais de cinco modelos físicos de pilares mistos parcialmente revestidos solicitados por compressão centrada e de um modelo solicitado por compressão excêntrica. Os resultados demonstraram-se promissores no que tange à substituição da armadura, mesmo para a situação com excentricidade na força aplicada. Em complemento, foi desenvolvida uma simulação numérica dos modelos físicos ensaiados utilizando o pacote computacional DIANA®, fundamentado no método dos elementos finitos, com o pré e pós-processador FX+. Apesar da hipótese simplificadora de aderência perfeita entre o aço e o concreto, o modelo numérico representou adequadamente o comportamento dos pilares mistos. / Partially encased composite columns are the result of the association of a steel profile and concrete filling working together. The most commonly used steel profiles have I or H cross sections not completely covered by concrete, with filling only in the region between flanges. In Brazil, the use of partially encased composite columns is still incipient. However, the introduction of design procedures for these columns in Brazilian Code ABNT - NBR 8800:2008 and also their inherent constructive advantages have enabled the use of these elements in buildings across the country, especially for large constructions. In this context, the main objective of this work is to analyze the structural behavior of these columns and evaluate the possibility of replacing the traditional reinforcement using steel rebar for alternatives such as welded steel meshes or steel fibers added to the concrete (rate of 1.5% and 25 mm long). To achieve these goals, experimental tests on five models of partially encased composite columns submitted to axial loading and one model subjected to eccentric loading were performed. The results proved to be promising in terms of replacing the traditional reinforcement, even for the situation with eccentricity of the applied force. In addition, a numerical modeling of the studied columns was developed, using finite element based software DIANA® with FX +. Despite simplifying assumptions of perfect bond between steel and concrete, the numerical model adequately represented the behavior of the columns.

A Roteirização de Veículos desempenha papel fundamental nos processos modernos de distribuição de produtos e realização de serviços. A atual disseminação de recursos de tecnologia de informação e comunicação, de forma confiável e economicamente acessível, permite trabalhar com informações em tempo real e melhoram os padrões de nível de serviço associados. O presente trabalho apresenta uma solução para roteirização de veículos cujas equipes de bordo realizam serviços que justificam seu deslocamento, uma vez que as demandas estão geograficamente dispersas. Tais demandas são, em parte, conhecidas antes do despacho (permitem programação antecipada) dos veículos e suas equipes; outra parte surge durante a jornada de trabalho. Como exemplos podem-se citar os casos de serviços de montagem e manutenção de instalações, equipamentos, engenharia e inspeção de tráfego, policiamento etc. Trata-se da aplicação da roteirização parcialmente dinâmica, conforme Larsen (2000), cujas bases foram definidas por Psaraftis (1988,1995), Bertsimas et al (1993) no problema DTRP (Dynamic Travelling Repairman Problem). A função objetivo apresenta uma combinação de minimização dos custos de deslocamento, para os pedidos de serviços conhecidos antes da saída dos veículos e de minimização do tempo de resposta (chegada no local do cliente ou da ocorrência) para os casos de pedidos imediatos ou emergenciais. A solução do problema envolve um modelo computacional de testes e avaliação, heurística de Clarke e Wright (1964) para formação das rotas estáticas, no Método Húngaro (Kuhn, 1955) para designar o veículo que resulta no menor tempo de resposta no atendimento a um pedido emergencial e a heurística de Clarke e Wright modificada na otimização do restante dos pedidos quando o veículo voltar a sua rota original. O modelo computacional foi testado em uma empresa de manutenção de elevadores na cidade de São Paulo, Brasil, onde demonstrou resultados comparativamente melhores em relação ao sistema de roteirização utilizado atualmente pela empresa. / The Vehicle Routing Problem plays a critical role on modern processes related to physical distribution of goods and services. The present expansion of information and communication technology in a reliable, economic and accessible way allows real time information and requires the utilization of appropriate tools for real time decisions resulting in significant improvements in quality and service level related to dynamic vehicle routing. A dynamic routing problem is presented, in which vehicles serve geographic dispersed service demands that justify their movement in a fixed area. Such service demands are partially known before vehicles dispatching (allowing prior programming) whilst others are known during the work journey. As examples, one can mention cases concerning installation and maintenance of utilities, equipment, engineering and surveillance services that refer to applications of Partially Dynamic Routing according to Larsen (2000), the groundings of which were defined by Psaraftis (1988,1995), Bertsimas et al (1993) in the Dynamic Travelling Repairman Problem (DTRP). The objective function is a combination of the minimization of movement costs to serve the prior demands and the minimization of time to reach (time to response) Dynamic-or-emergency-demand sites. The proposed solution involves a computational model for testing and evaluating a set of heuristics and methods comprising the Clarke and Wright (1964) Heuristic to compose the static routes, the Hungarian Method (Kuhn, 1955) to assign vehicles to the dynamic demands that produces the lowest response time and, finally, a Clarke and Wright Modified Heuristic used to optimize the remainder of the route when each diverted vehicle returns to its static route. The Computational Model was applied to a lift maintenance company located in the city of São Paulo (Brazil) demonstrating better results as compared to the present routing system.

[pt] O presente trabalho tem por objetivo estudar o comportamento de colunas esbeltas parcialmente enterradas, quando submetidas a um carregamento axial de compressão. A fundação é representada, ora por um modelo linear, o qual considera que a reação exercida pelo solo é proporcional às deflexões da coluna, ora por um modelo não-linear no qual esta relação de proporcionalidade não é mais verificada. Para a modelagem da coluna, é usada a teoria inextensional de barras esbeltas. Inicialmente, mostra-se como são deduzidas as equações diferenciais do problema a partir dos funcionais de energia da coluna. No problema linear, buscam-se obter as cargas críticas e modos críticos da coluna. Neste caso, sua solução analítica é encontrada a partir da resolução do problema de valor de contorno usando-se um programa de álgebra simbólica. Também é obtida uma solução aproximada através do método de Ritz. Um estudo paramétrico detalhado analisa a influência das condições de apoio da coluna e altura e rigidez da fundação na carga e modo críticos. Entretanto, no caso nãolinear, as equações diferenciais são mais complexas, não permitindo a obtenção de uma solução analítica. É utilizado, então, o método de Ritz, no qual as soluções analíticas obtidas para o problema linear (autofunções) são usadas como funções de interpolação. Em seguida, chega-se à uma equação não- linear de equilíbrio, da qual se obtém o caminho pós-crítico da coluna. Os resultados do problema nãolinear são comparados com os obtidos através do método dos elementos finitos. / [en] In this thesis the behavior of slender, partially embedded columns under axial compressive forces is studied. The foundation is either represented by a linear model, which considers that the soil reaction is proportional to the column’s deflections or by a non-linear model in which this proportionality relation is not observed. The inextensional slender beam theory is used to model the column. Initially, the governing differential equations are deduced from the energy functional of the column-foundation system. In the linear problem, the critical loads and corresponding critical modes are looked for. In this case, an analytic solution is obtained by the solution of the associated boundary value problem, using a symbolic algebra software. An approximate solution is also found by Ritz’s method. A parametric study is conducted to study the influence of the column boundary conditions and foundation’s height and stiffness on critical loads and modes. However, in the non-linear case, differential equations are much more complex and an analytical solution is not possible. So, the Ritz’s method is used once again, in which the analytic solutions of the linear problem (eigenfunctions) are used as interpolation functions. After that, a non-linear equilibrium equation is obtained together with the column post-buckling path. These results are compared with the ones obtained using the finite element method.