On the analytical solutions of the Hindmarsh-Rose neuronal model

In this article we analytically solve the Hindmarsh-Rose model (Proc R Soc Lond B221:87-102, 1984) by means of a technique developed for strongly nonlinear problems-the step homotopy analysis method. This analytical algorithm, based on a modification of the standard homotopy analysis method, allows us to obtain a one-parameter family of explicit series solutions for the studied neuronal model. The Hindmarsh-Rose system represents a paradigmatic example of models developed to qualitatively reproduce the electrical activity of cell membranes. By using the homotopy solutions, we investigate the dynamical effect of two chosen biologically meaningful bifurcation parameters: the injected current I and the parameter r, representing the ratio of time scales between spiking (fast dynamics) and resting (slow dynamics). The auxiliary parameter involved in the analytical method provides us with an elegant way to ensure convergent series solutions of the neuronal model. Our analytical results are found to be in excellent agreement with the numerical simulations.

From weak Allee effect to no Allee effect in Richards’ growth models

Population dynamics have been attracting interest since many years. Among the considered models, the Richards’ equations remain one of the most popular to describe biological growth processes. On the other hand, Allee effect is currently a major focus of ecological research, which occurs when positive density dependence dominates at low densities. In this chapter, we propose the dynamical study of classes of functions based on Richards’ models describing the existence or not of Allee effect. We investigate bifurcation structures in generalized Richards’ functions and we look for the conditions in the (β, r) parameter plane for the existence of a weak Allee effect region. We show that the existence of this region is related with the existence of a dovetail structure. When the Allee limit varies, the weak Allee effect region disappears when the dovetail structure also disappears. Consequently, we deduce the transition from the weak Allee effect to no Allee effect to this family of functions. To support our analysis, we present fold and flip bifurcation curves and numerical simulations of several bifurcation diagrams.

Projeto de um centro de roda para veículo de competição

A presente tese surgiu da cooperação entre a empresa Semog Racing e o Instituto Superior de Engenharia do Porto, no âmbito da unidade curricular de DPEST, englobado no 2º ano do Mestrado de Engenharia Mecânica, ramo de Construções Mecânicas. Esta dissertação tem como objetivo principal, o projeto de um novo centro de roda para substituir o existente num veículo de competição da Semog Racing. Este centro deverá ser dimensionado para suportar o peso do veículo e os esforços transmitidos à roda. Pretende-se igualmente que o novo modelo apresente um design apelativo e um baixo custo de produção para poder vir ser a comercializado. Na primeira fase, será estudado o centro de roda original fornecido pela Semog Racing. Este estudo engloba a modelação do componente através do software SolidWorks® e uma fase de simulações para diferentes condições de carregamento. O estudo será complementado com a realização de ensaios experimentais para validação do modelo numérico. A segunda fase é dedicada ao desenvolvimento do novo modelo de centro de roda, focando as características mecânicas e o design. Este tem como base o estudo numérico realizado para a roda original, tendo sempre como objetivo final garantir que o novo centro de roda cumpra todos os requisitos. O caminho seguido no processo de otimização é suportado em simulações numéricas pelo método dos elementos finitos, o qual permite aferir quais os pontos críticos a corrigir. No final, será apresentado um novo modelo de centro de roda capaz de suportar as cargas previstas de serviço, que apresente um baixo custo de fabrico e design apelativo.

Modelação e reforço de vigas de betão armado com laminados de fibra de carbono

Com o atual estado da construção em Portugal, a reabilitação urbana é uma realidade. Com muitos dos edifícios a necessitarem de reforço, procurou-se abordar o comportamento real das estruturas, indo além da típica análise linear elástica. Desta forma, pretendeu-se aumentar o conhecimento acerca da modelação numérica não-linear de estruturas de betão armado, expondo modelos de cálculo relativamente simples e de fácil compreensão, com o objetivo de servir de base a uma avaliação da capacidade de carga de um elemento estrutural. O modelo de cálculo foi validado com recurso ao trabalho experimental de Bresler e Scordelis (1963). Analisou-se o comportamento até à rotura de três vigas ensaiadas à flexão. Posteriormente, foi realizado um estudo paramétrico de algumas propriedades do betão com vista à discussão do melhor de ajuste. Em seguida, já no campo do reforço estrutural, simulou-se numericamente vigas reforçadas com CFRP, com recurso à técnica EBR e NSM. Comparam-se os resultados numéricos com os ensaios experimentais de Cruz et al. (2011a). Avaliou-se ainda o desempenho de soluções alternativas com variações na área e comprimento dos laminados. Para finalizar, foi desenvolvida uma campanha experimental com diferentes áreas de reforço. Conceberam-se e executaram-se três vigas de betão armado sobre as quais se instalaram laminados de CFRP. Os resultados experimentais são apresentados e discutidos à luz dos resultados do respetivo modelo numérico. No cômputo geral, o presente trabalho permitiu aferir a validade de modelos não-lineares na previsão do comportamento efetivo das estruturas até à rotura. Assinala-se a concordância em vários resultados experimentais analisados. Ficaram também patentes os principais fenómenos ligados ao reforço de vigas com CFRP, focados nos respetivos modelos de cálculo e nos resultados experimentais apresentados.

Propagação de vírus informáticos baseada em modelos biológicos

A evolução digital proporcionou às sociedades uma facilidade extraordinária de comunicação. Com o número crescente de computadores e o aumento de acessos à internet, surgiu uma nova forma de criminologia, que cresceu em paralelo com o número de utilizadores. Desta forma, tornou-se comum a criação e difusão de vírus informáticos pelos chamados hackers. Neste trabalho estudam-se modelos de transmissão de vírus informáticos, usando modelos epidemiológicos. Começa-se por fazer uma revisão dos modelos existentes na literatura, de seguida sugerem-se alterações a esses modelos de forma a conseguir uma melhor aproximação à dinâmica real de transmissão de vírus informáticos. As simulações numéricas dos modelos permitem-nos inferir de que uma forma de controlar a transmissão de vírus informáticos é a diminuição da taxa de infeção, isto é, da taxa de transmissão do vírus. No último capítulo enumeram-se as conclusões do trabalho efetuado e indicam-se direções de trabalho futuro.

Avaliação do desempenho de ligações tubulares em estruturas de veículos pesados de passageiros por simulação numérica

Nesta Dissertação ir-se-á avaliar o desempenho de algumas ligações existentes nos veículos pesados de passageiros, através do Regulamento Eurocódigo 3. Nos últimos anos ocorreram diversos acidentes envolvendo este tipos de veículos, em que os mesmos causaram vítimas mortais e feridos graves. Serão testadas por simulação numérica algumas ligações pertencentes a elementos constituintes da superestrutura, em que esta é normalmente afectada com a ocorrência de acidentes. Assim sendo, o estudo de nós de ligação tem uma importância fulcral para que uma superestrutura suporte situações extremas e que resista a solicitações externas aplicadas. Iniciou-se esta Dissertação com o estudo da sinistralidade e de acidentes que envolvem veículos pesados de passageiros. No capítulo 2 abordou-se um programa que promove simulações numéricas de acidentes e estudo do comportamento dos passageiros em caso de acidente, sendo referido o Regulamento que homologa os veículos pesados de passageiros e os seus principais métodos. Abordaram-se os principais constituintes da estrutura de um veículo pesado de passageiros. No capítulo 3, é referido o Eurocódigo 3 em termos do estudo das ligações tubulares usadas neste tipo de veículos. No capítulo 4, fez-se o estudo e selecção de elementos a utilizar para a simulação numérica de casos preconizados pelo Eurocódigo 3 e estudaram-se três tipos de ligações que são usadas na construção da superestrutura deste tipo de veículos, tendo-se retirado conclusões deste estudo.

Effects of treatment, awareness and condom use in a coinfection model for HIV and HCV in MSM

We develop a new a coinfection model for hepatitis C virus (HCV) and the human immunodeficiency virus (HIV). We consider treatment for both diseases, screening, unawareness and awareness of HIV infection, and the use of condoms. We study the local stability of the disease-free equilibria for the full model and for the two submodels (HCV only and HIV only submodels). We sketch bifurcation diagrams for different parameters, such as the probabilities that a contact will result in a HIV or an HCV infection. We present numerical simulations of the full model where the HIV, HCV and double endemic equilibria can be observed. We also show numerically the qualitative changes of the dynamical behavior of the full model for variation of relevant parameters. We extrapolate the results from the model for actual measures that could be implemented in order to reduce the number of infected individuals.

Strange Dynamics in a Fractional Derivative of Complex-Order Network of Chaotic Oscillators

We study the peculiar dynamical features of a fractional derivative of complex-order network. The network is composed of two unidirectional rings of cells, coupled through a "buffer" cell. The network has a Z3 × Z5 cyclic symmetry group. The complex derivative Dα±jβ, with α, β ∈ R+ is a generalization of the concept of integer order derivative, where α = 1, β = 0. Each cell is modeled by the Chen oscillator. Numerical simulations of the coupled cell system associated with the network expose patterns such as equilibria, periodic orbits, relaxation oscillations, quasiperiodic motion, and chaos, in one or in two rings of cells. In addition, fixing β = 0.8, we perceive differences in the qualitative behavior of the system, as the parameter c ∈ [13, 24] of the Chen oscillator and/or the real part of the fractional derivative, α ∈ {0.5, 0.6, 0.7, 0.8, 0.9, 1.0}, are varied. Some patterns produced by the coupled system are constrained by the network architecture, but other features are only understood in the light of the internal dynamics of each cell, in this case, the Chen oscillator. What is more important, architecture and/or internal dynamics?

Matrix fractional systems

This paper addresses the matrix representation of dynamical systems in the perspective of fractional calculus. Fractional elements and fractional systems are interpreted under the light of the classical Cole–Cole, Davidson–Cole, and Havriliak–Negami heuristic models. Numerical simulations for an electrical circuit enlighten the results for matrix based models and high fractional orders. The conclusions clarify the distinction between fractional elements and fractional systems.

A simplified four-branch model for the analytical study of the out-of-plane performance of regular stone URM walls

During the last years, several studies have been made aiming to assess the out-of-plane seismic response of unreinforced stone masonry structures. This fact led to the development of a wide variety of models and approaches, ranging from simple kinematic based analytical models up to complex numerical simulations. Nevertheless, for the sake of simplicity, the out-of-plane seismic response of a masonry wall pier may be obtained by means of a simple single-degree-of-freedom system while still providing good results. In fact, despite the assumptions associated with such a simple formulation, it is also true that the epistemic uncertainty inherent with the selection of appropriate input parameters in more complex models may render them truly ineffective. In this framework, this paper focuses on the study of the out-of-plane bending of unreinforced stone masonry walls (cantilevers) by proposing a simplified analytical approach based on the construction of a linearized four-branch model, which is used to characterize the linear and nonlinear response of such structural elements through an overturning moment-rotation relationship. The formulation of the four-branch model is presented and described in detail and the meaningful parameters used for its construction are obtained from a set of experimental laboratory tests performed on six full-scale unreinforced regular sacco stone masonry specimens. Moreover, a parametric analysis aiming to evaluate the effect of these parameters’ variation on the final configuration of the model is presented and critically discussed. Finally, the results obtained from the application of the developed four-branch model on real unreinforced regular sacco stone masonry walls are thoroughly analysed and the main conclusions obtained from its application are summarized.

Long-term behaviour of railway transitions under dynamic loading application to soft soil sites

Dissertação para obtenção do Grau de Doutor em Engenharia Civil

Urban deforestation and urban development

This paper has developed a model of a single forest owner operating with perfect foresight in a dynamic open-city environment that allows for switching between alternative competing land uses (forest and urban use) at some point in the future. The model also incorporates external values of an even-aged standing forest in addition to the value of timber when it is harvested. Timber is exploited based on a multiple rotation model a la Faustmann with clear-cut harvesting. In contrast to previous models, our alternative land use to forest land is endogenous. Within this framework, we study the problem of the private owner as well as that of the social planner, when choosing the time to harvest, the time to convert land and the intensity of development. We also examine the extent to which the two-way linkage between urban development and forest management practices (timber production and provision of forest amenities) contributes to economic efficiency and improvements in non-market forest benefits. Finally, we consider policy options available to a regulator seeking to achieve improvements in efficiency including anti-sprawl policies (impact fees and density controls) and forest policies such a yield tax. Numerical simulations illustrate our analytical results.

Optimal choice between even- and uneven-aged forestry: the case of non-industrial private forest owners

An infinite-horizon discrete time model with multiple size-class structures using a transition matrix is built to assess optimal harvesting schedules in the context of Non-Industrial Private Forest (NIPF) owners. Three model specifications accounting for forest income, financial return on an asset and amenity valuations are considered. Numerical simulations suggest uneven-aged forest management where a rational forest owner adapts her or his forest policy by influencing the regeneration of trees or adjusting consumption dynamics depending on subjective time preference and market return rate dynamics on the financial asset. Moreover she or he does not value significantly non-market benefits captured by amenity valuations relatively to forest income.

Polyhydroxyalkanoates production by photosynthetic mixed cultures

Polyhydroxyalkanoates (PHAs) are natural biologically synthesized polymers that have been the subject of much interest in the last decades due to their biodegradability. Thus far, its microbial production is associated with high operational costs, which increases PHA prices and limits its marketability. To address this situation, this thesis’ work proposes the utilization of photosynthetic mixed cultures (PMC) as a new PHA production system that may lead to a reduction in operational costs. In fact, the operational strategies developed in this work led to the selection of PHA accumulating PMCs that, unlike the traditional mixed microbial cultures, do not require aeration, thus permitting savings in this significant operational cost. In particular, the first PHA accumulating PMC tested in this work was selected under non-aerated illuminated conditions in a feast and famine regime, being obtained a consortium of bacteria and algae, where photosynthetic bacteria accumulated PHA during the feast phase and consumed it for growth during the famine phase, using the oxygen produced by algae. In this symbiotic system, a maximum PHA content of 20% cell dry weight (cdw) was reached, proving for the first time, the capacity of a PMC to accumulate PHA. During adaptation to dark/light alternating conditions, the culture decreased its algae content but maintained its viability, achieving a PHA content of 30% cdw. Also, the PMC was found to be able to utilize different volatile fatty acids for PHA production, accumulating up to 20% cdw of a PHA co-polymer composed of 3-hydroxybutyrate (3HB) and 3-hydroxyvalerate (HV) monomers. Finally, a new selective approach for the enrichment of PMCs in PHA accumulating bacteria was tested. Instead of imposing a feast and famine regime, a permanent feast regime was used, thus selecting a PMC that was capable of simultaneously growing and accumulating PHA, being attained a maximum PHA content of 60% cdw, the highest value reported for a PMC thus far. The results presented in this thesis prospect the utilization of cheap, VFA-rich fermented wastes as substrates for PHA production, which combined with this new photosynthetic technology opens up the possibility for direct sunlight illumination, leading to a more cost-effective and environmentally sustainable PHA production process.

Shear strengthening of reinforced concrete beams with SHCC-FRP panels

Tese de Doutoramento em Engenharia Civil