Esse trabalho consiste na determinação das velocidades críticas do vento e das amplitudes das vibrações numa estrutura composta por uma viga engastada suspensa por um estai (cabo), submetida aos efeitos de vento e chuva. Foi considerada a deformação no cabo devido ao carregamento do peso próprio e o acoplamento não-linear das vibrações do cabo e da viga. Três modos de vibração são de especial interesse, chamados de primeiro modo global (flexão da viga e vibração no cabo), primeiro modo local (vibração no cabo, com flexão na viga desprezável) e primeiro modo à torção. O modelo foi reduzido a três graus de liberdade. A modelagem dos carregamentos aerodinâmicos aplicados na viga foi feita segundo procedimentos tradicionais. O carregamento aerodinâmico aplicado ao cabo sob efeito de chuva e vento também foi levado em consideração. Para a redução do modelo matemático, os coeficientes de rigidez e de amortecimento equivalente são definidos e dependem parametricamente da velocidade do vento. Os termos não-lineares são devidos ao acoplamento das vibrações do cabo e da viga à flexão (no plano do cabo) e também aos efeitos aeroelásticos no cabo. Os seguintes regimes instáveis são avaliados: o efeito de galope (galloping) no cabo, o drapejamento (flutter) unimodal na torção e o drapejamento (flutter) do modo de flexão da viga em conjunto com vibrações transversais do cabo. / This paper is concerned with determining wind critical velocities and post-critical vibration amplitudes in a cable-stayed beam, under wind-rain condition. It is considered the cable sag due to the dead load plus the non-linear coupling between the vibration of both the cable and the beam. Three modes are of special interest, namely the first global mode (beam bending & cable vibration), the first local mode (cable vibration & negligible beam bending) and the first torsion mode. A reduced mathematical model, with three degrees of freedom, is also developed. With regard to the modelling of the aerodynamic loads applied to the beam, it can be performed after extension of classical guidelines. The aerodynamic loads applied to the cable under rain are also taken into account. For the reduced mathematical model, equivalent damping and stiffness coefficients will be defined, which depend parametrically on the wind velocity. Non-linear terms appear due to the coupling between the cable and the beam bending vibrations, and also to the aero-elastic non-linear effects on the cable. Different unstable regimes are surveyed such as the cable galloping, the unimodal flutter in torsion and the unimodal flutter with beam bending and cable vibrations coupled.

Model order reduction is necessary even in a time where the parallel processing is usual in almost any personal computer. The recent Model Reduction Methods are useful tools nowadays on reducing the problem processing. This work intends to describe a combination between POD (Proper Orthogonal Decomposition) and Ritz vectors that achieve an efficient Galerkin projection that changes during the processing, comparing the development of the error and the convergence rate between the full space and the projection space, in addition to check the stability of the projection space, leading to an adaptive model order reduction for nonlinear dynamical problems more efficient. This model reduction is supported by a secant formulation, which is updated by BFGS (Broyden - Fletcher - Goldfarb - Shanno) method to accelerate convergence of the model, and a tangent formulation to correct the projection space. Furthermore, this research shows that this method permits a correction of the reduced model at low cost, especially when the classical POD is no more efficient to represent accurately the solution. / A Redução de ordem de modelo é necessária, mesmo em uma época onde o processamento paralelo é usado em praticamente qualquer computador pessoal. Os recentes métodos de redução de modelo são ferramentas úteis nos dias de hoje para a redução de processamento de um problema. Este trabalho pretende descrever uma combinação entre POD (Proper Orthogonal Decomposition) e vetores de Ritz para uma projecção de Galerkin eficiente que sofre alterações durante o processamento, comparando o desenvolvimento do erro e a taxa de convergência entre o espaço total e o espaço de projeção, além da verificação de estabilidade do espaço de projeção, levando a uma redução de ordem do modelo adaptativo mais eficiente para problemas dinâmicos não-lineares. Esta redução de modelo é assistida por uma formulação secante, que é atualizado pela formula de BFGS (Broyden - Fletcher- Goldfarb - Shanno) com o intuito de acelerar a convergência do modelo, e uma formulação tangente para a correção do espaço de projeção. Além disso, esta pesquisa mostra que este método permite a correção do modelo reduzido com baixo custo, especialmente quando o clássico POD não é mais eficiente para representar com precisão a solução.

Desenvolve-se uma formulação lagrangeana total do método dos elementos finitos para análise dinâmica de estruturas e mecanismos reticulados planos contendo ligações deslizantes sujeitas a grandes deslocamentos e rotações. Estas são introduzidas ao sistema mecânico na forma de juntas prismáticas e cilíndricas por meio do método dos multiplicadores de Lagrange, permitindo sua utilização na simulação de diversos tipos de estruturas e mecanismos. Também são consideradas rótulas entre as barras, estas introduzidas por meio da compatibilidade cinemática dos graus de liberdade dos nós comuns. A formulação do método dos elementos finitos adotada utiliza como parâmetros nodais as posições e os giros de modo desacoplado. Assim, pode-se utilizar a cinemática exata para barras de Reissner na análise de deslocamentos e giros finitos da estrutura. Adota-se o modelo constitutivo de Saint-Venant-Kirchhoff que relaciona a medida de deformação objetiva de Green-Lagrange com o tensor de tensões de Piola-Kirchhoff de segunda espécie. O equilíbrio dinâmico do sistema é obtido pelo princípio da energia total estacionária e a solução do sistema não linear de equações resultante é obtida pelo método de Newton-Raphson. A integração temporal é realizada pelo método de Newmark. São apresentados diversos exemplos para validação da formulação desenvolvida, os quais são comparados com soluções analíticas de modo a evidenciar as possibilidades de aplicação da formulação proposta. / A total lagrangian finite element method formulation is developed for the dynamic analysis of plane frame structures and mechanisms containing sliding joints that undergoes large displacements and rotations. Those connections are introduced in the mechanical system as prismatic and cylindrical joints by the method of Lagrange multipliers, allowing its use on the simulation of several types of structures and mechanisms. Hinges between bars are also considered by kinematic compatibility of the degrees of freedom on the common node. The adopted finite element formulation uses as nodal parameters uncoupled positions and angles. Therefore, Reissner exact kinematics for bars can be utilized for structural finite deformation. The Saint-Venant-Kirchhoff constitutive model, which relates the objective Green-Lagrange strain measure with the second Piola-Kirchhoff stress tensor, is adopted. The principle of stationary total energy is used to obtain the dynamic nonlinear equilibrium of the system and the solution of the resulting nonlinear system of equations is done by the Newton-Raphson method. The Newmark method is adopted for time integration. Several examples are presented for the validation of the developed formulation, and those are compared with analytical solutions in order to clarify the possibilities of application of the proposed formulation.

O trabalho trata da formulação e implementação numérica de um modelo de dano para o concreto. O objetivo é o de reproduzir laços de histerese observados experimentalmente em diagramas de tensão-deformação uniaxiais quando, a partir de um certo nível de solicitação, e dano correspondente, passam-se a realizar ciclos de descarregamento e recarregamento. O modelo propõe que os laços de histerese resultam da dissipação conjunta decorrente do dano evolutivo e do fenômeno de fricção interna nas faces das fissuras. O fenômeno de fricção interna está relacionado à suposição de que as fissuras não se abrem simplesmente por separação formando faces com superfícies lisas, mas tendem a combinar modos de separação e deslizamento segundo superfícies com uma certa rugosidade. O deslizamento entre as faces das fissuras dá origem a uma deformação por escorregamento, assumida como responsável por um comportamento plástico com encruamento cinemático não-linear. A abordagem termodinâmica considerada no sentido de contemplar a danificação e a fricção interna leva à dedução de que a tensão total num ponto do meio com dano resulta dividida em uma parcela dita tensão elasto-danificada e em uma outra denominada tensão de escorregamento, também dependente do nível de danificação da estrutura. Dois tipos básicos de aplicações são realizadas: as primeiras visam avaliar a influência das variáveis de estado e dos parâmetros do material na forma dos diagramas de tensão-deformação. As outras aplicações são dedicadas à utilização do modelo na análise do comportamento dinâmico de uma viga em concreto armado submetida à vibração forçada. Destaca-se o efeito dos laços de histerese sobre o amortecimento estrutural. Mostra-se que a utilização deste modelo de dano leva à geração de um amortecimento estrutural que com os modelos de dano usuais somente poderia ser obtido a partir da utilização do modelo de amortecimento de Rayleigh, por exemplo, na equação do movimento. / This work deals with numerical formulation and implementation of a damage model to concrete. The aim is to reproduce hysteresis loops observed experimentally in uniaxial stress-strain diagram when cycles of unloading and reloading are produced. The model proposes that the hysteresis loops are due to the coupling between the dissipation of damage process and the internal friction phenomenon between the faces of crack lips. The phenomenon of internal friction is related to the supposition that the crack doesn\'t open itself by separation mode creating plane lips, but there is a tendency to combine opening and sliding modes across irregular surfaces. The sliding between the crack surfaces generates a sliding strain, which is assumed as responsible by a plastic behaviour with non-linear kinematics hardening. The thermodynamic approach considering damage and internal friction leads to the deduction that the total stress at a point is composted by one component named elasto-damage stress and another named sliding stress, both dependent of the structure damage level. Two basic kinds of simulation are performed: the first one aims to evaluate the influence of the state variables and the material parameters over the stress-strain diagrams. The following simulations are dedicated to the use of the model in dynamic analysis of a reinforced concrete beam subjected to vibration. It is emphasized the effect of the structural damping from hysteresis loops. It is showed that the present damage model generates structural damping that could only be obtained with conventional damage model by using the Rayleigh damping approach in the equation of motion.

Orientador: Átila Madureira Bueno / Resumo: O som é uma onda mecânica e como tal transporta energia que age sobre partículas devido às forças de radiação acústica. O princípio para suspender corpos é aplicar uma força de tal forma a equilibrar seu peso. Na técnica de levitação acústica (AcLev) uma pequena esfera pode ser suspensa pela força de radiação acústica gerada por uma onda estacionária, sendo que o ponto de levitação está localizado na região em que o potencial acústico é mínimo, que é condição necessária para levitar uma esfera com raio muito menor que o comprimento de onda. Levitação acústica (AcLev) é uma ferramenta importante para manusear objetos sem contêineres. Nos anos recentes muitos dispositivos foram desenvolvidos com sucesso devido ao comportamento estável dos dispositivos AcLev. Como resultado, a maioria dos trabalhos sobre Aclev se concentram sobre simulações numéricas ou testes experimentais para estudar a geometria e arranjos dos emissores acústicos, ou a influência de vários tipos de perturbações, e a maioria desses modelos matemáticos considera somente o potencial acústico. Neste trabalho, a equação não linear de movimento para uma partícula levitada imersa em campo acústico de eixo único foi desenvolvida, considerando também forças dissipativas. O espaço parâmetro foi examinado buscando a existência de bifurcações, e faixas de projeto para os ganhos do dispositivo AcLev foram determinadas a partir da condição de existência de pontos de equilíbrio. Em adição, o comportamento dinâmico do dispos... (Resumo completo, clicar acesso eletrônico abaixo) / Abstract: Sound is a mechanical wave and aims to carry energy that acts on particles due to acoustic radiation forces, while the principle to suspend bodies is to apply a force in such a way as to balance their weight. In the acoustic levitation technique (Aclev) a small sphere can be suspended by the acoustic radiation force generated by a stationary wave and the levitation point is located in the region where the acoustic potential is minimal, which is a necessary condition for levitating a sphere with radius much smaller than the wavelength. AcLev is an important tool for handling objects without the use of containers. In recent years many devices have been successfully developed due to the stable behavior of AcLev. As a result, most works on Aclev focuses on numerical simulations and experimental tests to study the geometry and arrangement of acoustic emitters, or the influence of various kinds of perturbations, and most mathematical models consider only acoustic potential. In this work, the nonlinear equation of motion for a levitated particle immersed in an acoustic field with single axis was developed considering also dissipative forces. The parameter space was searched for the existence of bifurcations and the design range for AcLev device gains were determined from the condition of equilibrium points. In addition, the dynamic behavior of the AcLev device regarding gains has been studied, also considering the microgravity situation. Numerical simulations corroborated the analyt... (Complete abstract click electronic access below) / Mestre

Nesta tese, apresenta-se a análise numérica e experimental do comportamento dinâmico não linear de um pórtico excitado por uma fonte não ideal - um motor elétrico de corrente contínua desbalanceado - como exemplo dessa classe de problemas. Elabora-se um modelo matemático com quatro graus de liberdade: dois relacionados com o deslocamento horizontal e vertical do ponto de apoio da máquina e dois com os parâmetros de funcionamento do motor elétrico. Adota-se a formulação Lagrangeana para gerar as equações de movimento contendo termos não lineares até ordem cúbica. A solução numérica é obtida através do método de Runge-Kutta com passo adaptativo. Ensaios dinâmicos e estáticos foram realizados com o motor e com o sistema completo, com vistas a validar o modelo matemático. Um sistema de medição digital armazenou todos os dados obtidos em arquivos de fácil leitura pelos programas atualmente disponíveis. Filmagens de alta velocidade e fotografia estroboscópica registraram movimentos típicos relacionados com comportamentos não lineares. Os resultados numéricos e experimentais mostram boa correlação entre si, além de apresentarem alguns dos fenômenos associados ao se forçar uma estrutura a passar por uma de suas ressonâncias excitando-a com um dispositivo com potência limitada, como o efeito Sommerfeld. Outros fenômenos, devido ao comportamento geometricamente não linear da estrutura, são também detectados, tais como saturação modal e transferência de energia. / Numerical and experimental analysis of the nonlinear dynamics of a portal frame excited by a non-ideal source - an unbalanced direct current motor - is presented in this thesis as an example of this class of problems. A four degree of freedom model is elaborated: two of them related to the horizontal and vertical structural displacements and two others to the functioning parameters of the motor. A Lagrangian approach for deducing the equations of motion, up to cubic non-linear terms is followed. The numerical solution is obtained through Runge-Kutta algorithm with adaptive step. Static and dynamic tests were performed with the motor and with the complete system, in order to validate the mathematical model. A digital acquisition system recorded all data in computer files, ready to be read by available commercial programs. High speed filming and stroboscopic photography were used to register typical movements related to nonlinear behavior. Numerical and mathematical results show good correlation, as well as present some phenomena related to passage through resonance with limited power excitation such as the Sommerfeld effect. Some other phenomena, related to the nonlinear behavior of the structure are also detected, as modal saturation and energy transference.

O trabalho trata da formulação e implementação numérica de modelos matemáticos do comportamento de estruturas considerando-se as não-linearidades física e geométrica. O equilíbrio na posição deslocada é formulado via Princípio dos Trabalhos Virtuais, empregando-se o método dos elementos finitos para a discretização espacial das estruturas e busca de soluções aproximadas. Inicialmente destaca-se com base no caso de treliças planas o emprego de medidas de deformação e tensão conjugadas energeticamente. Particularizando-se a formulação geral do equilíbrio para os pórticos planos apresenta-se uma análise crítica das formulações lagrangiana total e atualizada. Em seguida, tendo-se em vista aplicações às estruturas em concreto armado, aborda-se o comportamento não-linear físico pela mecânica do dano em meios contínuos, empregando-se os modelos de dano para o concreto propostos por Mazars e La Borderie. Estendem-se os estudos do comportamento estrutural não-linear físico (dano) e geométrico incorporando-se a análise dinâmica. Utiliza-se para integração no domínio do tempo o método implícito de Newmark combinado com o procedimento incremental e iterativo de Newton-Raphson. O amortecimento é levado em conta por meio da regra de Rayleigh. Exemplos consistindo de análises não lineares estática e dinâmica de estruturas reticulares planas compõem numéricas. Os resultados obtidos ilustram o desempenho e as potencialidades das formulações empregadas. / The formulation and numerical implementation of mathematical models of the structural behavior of plane frames considering material and geometrical nonlinearities are treated in this work. The Principle of Virtual Work is presented in order to characterize the equilibrium in the displaced position. The structure is discretized by the finite element method. Plane trusses analysis is performed in order to show how important it is to take into account the strain and stress tensors energetically conjugated. Being particularized, the general equilibrium formulation for the plane frames a critical analysis of the total and updated lagrangian formulations is presented. The material nonlinear behavior is modeled by continuum damage mechanics by using both Mazars and La Boderie\'s damage models. The nonlinear analysis is extended including the dynamic response of the reinforced concrete plane frames. The usual iterative Newton-Raphson technique is used combined with implicit Newmark method in order to carry out the integration at time. The damping is introduced by means of the Rayleigh\'s rule. Numerical analysis by a suitable computer program show the theoretical results considering static and dynamic response of plane frames.

O trabalho se baseia em uma formulação numérica (Método dos Elementos Finitos Posicional) que é combinada a um modelo viscoelástico adequado (Kelvin-Voigt adaptado), o que direciona para o cumprimento do objetivo: a simulação de sólidos viscoelásticos em deformações finitas e de fluidos viscosos. A formulação desenvolvida é Lagrangeana total descrita para posições, permitindo aplicações em não linearidade dinâmica (com a utilização do método de Newton-Raphson para solução do sistema de equações não lineares e integração temporal via algoritmo implícito de Newmark) e sua combinação com um modelo viscoelástico coerente é deduzida neste trabalho. Inicialmente, são resolvidos problemas com elemento de chapa bidimensional, porém o elemento finito final utilizado é de sólido prismático de base triangular. Dois modelos são adotados para consideração do comportamento viscoelástico, (i) um modelo modificado de Kelvin-Voigt associado ao modelo constitutivo de Saint-Venant-Kirchhoff e (ii) um modelo visco-hiperelástico completo coerente para deformações finitas desenvolvido a partir da decomposição multiplicativa sobre o gradiente da função mudança de configuração em uma parcela volumétrica e duas isocóricas. Foram selecionados e comentados 15 exemplos em detalhe, abrangendo todas as etapas desta pesquisa, com problemas elásticos, dinâmicos, viscoelásticos em pequenas e grandes deformações, de flexão, de impacto e de fluidos viscosos. Os resultados obtidos para os exemplos de validação foram satisfatórios, coerentes com as referências, e o conjunto das análises conduzidas mostram a potencialidade da formulação alternativa desenvolvida neste trabalho. / The work is based on a numerical formulation (Positional Finite Element Method) combined with a suitable viscoelastic model (adapted Kelvin-Voigt), what directs to achieving its main goal: the simulation of viscoelastic solids in finite strain and of viscous fluids. The developed formulation is total Lagrangian described for positions, allowing applications in nonlinear dynamics (using the Newton-Raphson method for solution of the system of nonlinear equations, and performing time integration via an implicit Newmark algorithm); its combination with adequate viscoelastic model is shown step-by-step in this work. Initially, problems are solved using two-dimensional plate element, but the final finite element is a triangular-based prismatic solid. Two models are adopted in order to consider the viscoelastic behavior, (i) a modified Kelvin-Voigt model associated with the Saint-Venant-Kirchhoff constitutive model, and (ii) a coherent visco-hyperelastic model for finite deformations developed from the multiplicative decomposition over the deformation gradient in one volumetric and two isochoric parts. 15 examples were selected and commented in detail, comprehending all stages of this research, solving problems that are elastic, dynamic, viscoelastic under small and large strain, under flexural behavior, submitted to impact, and of viscous fluids problems. The results obtained for the validation examples were satisfactory, consistent with the references, and the whole of the conducted analysis shows the potentials of the alternative formulation developed in this work.

Introdução: as infecções são responsáveis por mais de um terço das mortes neonatais, sendo a sepse a de maior expressão. A sepse neonatal precoce (SNP) aumenta em cinco vezes a taxa de mortalidade entre recém-nascidos (RN) de termo e pré-termos tardios e em vinte e cinco vezes a taxa de mortalidade entre pré-termos extremos. A Variabilidade da Frequência Cardíaca (VFC) vem se mostrando capaz de antecipar o diagnóstico, reduzindo em 22% o risco relativo de mortalidade entre os pré-termos de muito baixo peso com sepse. Objetivo: verificar a VFC de RN com suspeita de SNP; identificar medidas significantes e verificar se, associadas a variáveis clínico-laboratoriais, aumentam a chance de identificar pacientes passíveis de desenvolver a doença. Método: utilizou-se um frequencímetro portátil para aferir as medidas da VFC nos domínios do tempo, da frequência e do Caos. Dentre elas, elegeu-se aquelas com curva ROC > 0,75 para, com seus valores de cut-off, aplicar regressão logística e identificar as que melhor se adequassem ao modelo. As variáveis clínico-laboratoriais significantes que restaram da análise univariada foram submetidas à análise multivariada a fim de identificar, também, as que melhor se adequassem ao modelo. Resultados: foram identificadas quatro variáveis clínico-laboratoriais (sexo masculino, idade gestacional < 34 semanas, Score de Rodwell e proteína C-reativa), e nove medidas da VFC (RRmean, SDNN, RR tri index, LF, HF, SD2 do plot de Poincaré, ShaEnt, CorrDim D2 e SamAsy). Cinco apresentaram área sob a curva ROC > 0,75: SDNN, RR tri index, LF, SD2 e CorrDim D2. Depois de aplicada a regressão logística, restou apenas a CorrDim D2 com cut-off = 0,1164. Identificadas, as variáveis foram submetidas ao cálculo do logit para, com isso, quantificar a chance do paciente suspeito de evoluir para a doença, na dependência das três variáveis eleitas: CorrDim D2, idade gestacional < 34 semanas e score de Rodwell. Conclusão: o estudo mostrou que tanto o modelo preditivo clinico como a CorrDim D2 são capazes, de forma independente, de estimar a chance que um futuro paciente com sepse suspeitada tem de efetivamente confirmar o diagnóstico de sepse. E que, além disso, se a CorrDim D2 for agregada ao modelo preditivo clínico, esta chance aumenta consideravelmente. / Introduction: infectious diseases are responsible for more than one-third of neonatal deaths, and sepsis is of great importance. Early Neonatal Sepsis (ENS) increases five times the mortality rate among term and late preterm newborns. Among extreme preterm babies the mortality rate increases by twenty-five times. Heart Rate Variability (HRV) has been showing to be capable of anticipating the diagnosis, thus reducing by 22 percent the relative risk of mortality of very low preterm babies with sepsis. Objective: to verify HRV of newborn babies suspects of having ENS; to identify significant measures of HRV and to verify if, associated to clinical-laboratorial variables, if they improve the chance of identifying patients exposed to developing ENS. Method: a portable frequencimeter was used to registering the measurements of HRV in the domains of time, of frequency and of chaos. Among them, only those with a ROC curve greater than 0.75 were chosen, with their cut-off values to enter a logistic regression with the aim of identifying those more suitable to the model. The clinical-laboratorial variables which were considered significant in a univariete analysis were also submitted to the multivariate analysis with the purpose of identifying those better fitted to the model. Results: four clinical-laboratorial variables (male sex, gestational age less than 34 weeks, the Rodwell Score and C-Reactive Protein); similarly, nine measures of HRV (RRmean, SDNN, RR tri index, LF, HF, SD2 of Poincaré plot, ShaEnt, CorrDimD2 and SamAsy). Five of them presented an area under the ROC curve greater than 0.75: SDNN, RR tri index, LF, SD2 and CorrDim D2. Once applied the logistic regression, only CorrDim D2 remained, with a cut-off value of 0.1164. The identified variables were, then, submitted to the calculation of the logit, with the purpose of quantifying the chance of a patient suspect of evolving to a sepsis, in the dependence of the three elected variables: CorrDim D2, gestational age under 34 weeks and the Score of Rodwell. Conclusion: the study showed that both, the predictive clinical model as CorrDim D2 are able, independently, to estimate the chance that a patient with suspected sepsis future has to effectively confirm the diagnosis of sepsis. And that, moreover, if the CorrDim D2 is aggregated to clinical predictive model, this chance increases considerably.

Introdução: as infecções são responsáveis por mais de um terço das mortes neonatais, sendo a sepse a de maior expressão. A sepse neonatal precoce (SNP) aumenta em cinco vezes a taxa de mortalidade entre recém-nascidos (RN) de termo e pré-termos tardios e em vinte e cinco vezes a taxa de mortalidade entre pré-termos extremos. A Variabilidade da Frequência Cardíaca (VFC) vem se mostrando capaz de antecipar o diagnóstico, reduzindo em 22% o risco relativo de mortalidade entre os pré-termos de muito baixo peso com sepse. Objetivo: verificar a VFC de RN com suspeita de SNP; identificar medidas significantes e verificar se, associadas a variáveis clínico-laboratoriais, aumentam a chance de identificar pacientes passíveis de desenvolver a doença. Método: utilizou-se um frequencímetro portátil para aferir as medidas da VFC nos domínios do tempo, da frequência e do Caos. Dentre elas, elegeu-se aquelas com curva ROC > 0,75 para, com seus valores de cut-off, aplicar regressão logística e identificar as que melhor se adequassem ao modelo. As variáveis clínico-laboratoriais significantes que restaram da análise univariada foram submetidas à análise multivariada a fim de identificar, também, as que melhor se adequassem ao modelo. Resultados: foram identificadas quatro variáveis clínico-laboratoriais (sexo masculino, idade gestacional < 34 semanas, Score de Rodwell e proteína C-reativa), e nove medidas da VFC (RRmean, SDNN, RR tri index, LF, HF, SD2 do plot de Poincaré, ShaEnt, CorrDim D2 e SamAsy). Cinco apresentaram área sob a curva ROC > 0,75: SDNN, RR tri index, LF, SD2 e CorrDim D2. Depois de aplicada a regressão logística, restou apenas a CorrDim D2 com cut-off = 0,1164. Identificadas, as variáveis foram submetidas ao cálculo do logit para, com isso, quantificar a chance do paciente suspeito de evoluir para a doença, na dependência das três variáveis eleitas: CorrDim D2, idade gestacional < 34 semanas e score de Rodwell. Conclusão: o estudo mostrou que tanto o modelo preditivo clinico como a CorrDim D2 são capazes, de forma independente, de estimar a chance que um futuro paciente com sepse suspeitada tem de efetivamente confirmar o diagnóstico de sepse. E que, além disso, se a CorrDim D2 for agregada ao modelo preditivo clínico, esta chance aumenta consideravelmente. / Introduction: infectious diseases are responsible for more than one-third of neonatal deaths, and sepsis is of great importance. Early Neonatal Sepsis (ENS) increases five times the mortality rate among term and late preterm newborns. Among extreme preterm babies the mortality rate increases by twenty-five times. Heart Rate Variability (HRV) has been showing to be capable of anticipating the diagnosis, thus reducing by 22 percent the relative risk of mortality of very low preterm babies with sepsis. Objective: to verify HRV of newborn babies suspects of having ENS; to identify significant measures of HRV and to verify if, associated to clinical-laboratorial variables, if they improve the chance of identifying patients exposed to developing ENS. Method: a portable frequencimeter was used to registering the measurements of HRV in the domains of time, of frequency and of chaos. Among them, only those with a ROC curve greater than 0.75 were chosen, with their cut-off values to enter a logistic regression with the aim of identifying those more suitable to the model. The clinical-laboratorial variables which were considered significant in a univariete analysis were also submitted to the multivariate analysis with the purpose of identifying those better fitted to the model. Results: four clinical-laboratorial variables (male sex, gestational age less than 34 weeks, the Rodwell Score and C-Reactive Protein); similarly, nine measures of HRV (RRmean, SDNN, RR tri index, LF, HF, SD2 of Poincaré plot, ShaEnt, CorrDimD2 and SamAsy). Five of them presented an area under the ROC curve greater than 0.75: SDNN, RR tri index, LF, SD2 and CorrDim D2. Once applied the logistic regression, only CorrDim D2 remained, with a cut-off value of 0.1164. The identified variables were, then, submitted to the calculation of the logit, with the purpose of quantifying the chance of a patient suspect of evolving to a sepsis, in the dependence of the three elected variables: CorrDim D2, gestational age under 34 weeks and the Score of Rodwell. Conclusion: the study showed that both, the predictive clinical model as CorrDim D2 are able, independently, to estimate the chance that a patient with suspected sepsis future has to effectively confirm the diagnosis of sepsis. And that, moreover, if the CorrDim D2 is aggregated to clinical predictive model, this chance increases considerably.