Active and passive behaviors of soft tissues: Pelvic floor muscles

: A new active-contraction visco-elastic numerical model of the pelvic floor (skeletal) muscle is presented. Our model includes all elements that represent the muscle constitutive behavior, contraction and relaxation. In contrast with the previous models, the activation function can be null. The complete equations are shown and exactly linearized. Small verification and validation tests are performed and the pelvis is modeled using the data from the intra-abdominal pressure tests

Derivative-free optimization and filter methods to solve nonlinear constrained problems

In real optimization problems, usually the analytical expression of the objective function is not known, nor its derivatives, or they are complex. In these cases it becomes essential to use optimization methods where the calculation of the derivatives, or the verification of their existence, is not necessary: the Direct Search Methods or Derivative-free Methods are one solution. When the problem has constraints, penalty functions are often used. Unfortunately the choice of the penalty parameters is, frequently, very difficult, because most strategies for choosing it are heuristics strategies. As an alternative to penalty function appeared the filter methods. A filter algorithm introduces a function that aggregates the constrained violations and constructs a biobjective problem. In this problem the step is accepted if it either reduces the objective function or the constrained violation. This implies that the filter methods are less parameter dependent than a penalty function. In this work, we present a new direct search method, based on simplex methods, for general constrained optimization that combines the features of the simplex method and filter methods. This method does not compute or approximate any derivatives, penalty constants or Lagrange multipliers. The basic idea of simplex filter algorithm is to construct an initial simplex and use the simplex to drive the search. We illustrate the behavior of our algorithm through some examples. The proposed methods were implemented in Java.

Análise estática não linear de um edifício em betão armado

Trabalho Final de Mestrado para obtenção do grau de Mestre em Engenharia Civil

Exciting dynamical behavior in a network of two coupled rings of Chen oscillators

We study exotic patterns appearing in a network of coupled Chen oscillators. Namely, we consider a network of two rings coupled through a “buffer” cell, with Z3×Z5 symmetry group. Numerical simulations of the network reveal steady states, rotating waves in one ring and quasiperiodic behavior in the other, and chaotic states in the two rings, to name a few. The different patterns seem to arise through a sequence of Hopf bifurcations, period-doubling, and halving-period bifurcations. The network architecture seems to explain certain observed features, such as equilibria and the rotating waves, whereas the properties of the chaotic oscillator may explain others, such as the quasiperiodic and chaotic states. We use XPPAUT and MATLAB to compute numerically the relevant states.

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.

Evaluation of seismic vulnerability assessment parameters for portuguese vernacular constructions with nonlinear numerical analysis

Considering that vernacular architecture may bear important lessons on hazard mitigation and that well-constructed examples showing traditional seismic resistant features can present far less vulnerability than expected, this study aims at understanding the resisting mechanisms and seismic behavior of vernacular buildings through detailed finite element modeling and nonlinear static (pushover) analysis. This paper focuses specifically on a type of vernacular rammed earth constructions found in the Portuguese region of Alentejo. Several rammed earth constructions found in the region were selected and studied in terms of dimensions, architectural layout, structural solutions, construction materials and detailing and, as a result, a reference model was built, which intends to be a simplified representative example of these constructions, gathering the most common characteristics. Different parameters that may affect the seismic response of this type of vernacular constructions have been identified and a numerical parametric study was defined aiming at evaluating and quantifying their influence in the seismic behavior of this type of vernacular buildings. This paper is part of an ongoing research which includes the development of a simplified methodology for assessing the seismic vulnerability of vernacular buildings, based on vulnerability index evaluation methods.

In vitro mechanical fatigue behavior of poly-ε-caprolactone macroporous scaffolds for cartilage tissue engineering: Influence of pore filling by a poly(vinyl alcohol) gel

Polymeric scaffolds used in regenerative therapies are implanted in the damaged tissue and subjected to repeated loading cycles. In the case of articular cartilage engineering, an implanted scaffold is typically subjected to long term dynamic compression. The evolution of the mechanical properties of the scaffold during bioresorption has been deeply studied in the past, but the possibility of failure due to mechanical fatigue has not been properly addressed. Nevertheless, the macroporous scaffold is susceptible to failure after repeated loading-unloading cycles. In this work fatigue studies of polycaprolactone scaffolds were carried by subjecting the scaffold to repeated compression cycles in conditions simulating the scaffold implanted in the articular cartilage. The behaviour of the polycaprolactone sponge with the pores filled with a poly(vinyl alcohol) gel simulating the new formed tissue within the pores was compared with that of the material immersed in water. Results were analyzed with Morrow’s criteria for failure and accurate fittings are obtained just up to 200 loading cycles. It is also shown that the presence of poly(vinyl alcohol) increases the elastic modulus of the scaffolds, the effect being more pronounced with increasing the number of freeze/thawing cycles.

Pure elastic contact force models

"Series: Solid mechanics and its applications, vol. 226"

Contact force models for multibody dynamics

"Series: Solid mechanics and its applications, vol. 226"

An objective anisotropic elastic plastic model and Algorithm applicable to bone mechanics

The implicit projection algorithm of isotropic plasticity is extended to an objective anisotropic elastic perfectly plastic model. The recursion formula developed to project the trial stress on the yield surface, is applicable to any non linear elastic law and any plastic yield function.A curvilinear transverse isotropic model based on a quadratic elastic potential and on Hill's quadratic yield criterion is then developed and implemented in a computer program for bone mechanics perspectives. The paper concludes with a numerical study of a schematic bone-prosthesis system to illustrate the potential of the model.

A mixed finite element method for nonlinear diffusion equations

We propose a mixed finite element method for a class of nonlinear diffusion equations, which is based on their interpretation as gradient flows in optimal transportation metrics. We introduce an appropriate linearization of the optimal transport problem, which leads to a mixed symmetric formulation. This formulation preserves the maximum principle in case of the semi-discrete scheme as well as the fully discrete scheme for a certain class of problems. In addition solutions of the mixed formulation maintain exponential convergence in the relative entropy towards the steady state in case of a nonlinear Fokker-Planck equation with uniformly convex potential. We demonstrate the behavior of the proposed scheme with 2D simulations of the porous medium equations and blow-up questions in the Patlak-Keller-Segel model.

A numerical solver for a nonlinear Fokker-Planck equation representation of neuronal network dynamics

To describe the collective behavior of large ensembles of neurons in neuronal network, a kinetic theory description was developed in [13, 12], where a macroscopic representation of the network dynamics was directly derived from the microscopic dynamics of individual neurons, which are modeled by conductance-based, linear, integrate-and-fire point neurons. A diffusion approximation then led to a nonlinear Fokker-Planck equation for the probability density function of neuronal membrane potentials and synaptic conductances. In this work, we propose a deterministic numerical scheme for a Fokker-Planck model of an excitatory-only network. Our numerical solver allows us to obtain the time evolution of probability distribution functions, and thus, the evolution of all possible macroscopic quantities that are given by suitable moments of the probability density function. We show that this deterministic scheme is capable of capturing the bistability of stationary states observed in Monte Carlo simulations. Moreover, the transient behavior of the firing rates computed from the Fokker-Planck equation is analyzed in this bistable situation, where a bifurcation scenario, of asynchronous convergence towards stationary states, periodic synchronous solutions or damped oscillatory convergence towards stationary states, can be uncovered by increasing the strength of the excitatory coupling. Finally, the computation of moments of the probability distribution allows us to validate the applicability of a moment closure assumption used in [13] to further simplify the kinetic theory.

The renin-angiotensin system and arterial wall behavior.

To assess the behavior of the arterial wall in hypertensive patients, we developed a noninvasive ultrasonic device. Simultaneous recordings of internal diameter and blood pressure over the whole cardiac cycle are used to establish compliance-pressure curves. Blood pressure, which is a co-determinant of compliance, is thus taken into account. This method allows one to compare arteries from patients with different blood pressures. Arterial compliance and distensibility were first investigated in healthy young volunteers administered either lisinopril (20 mg), atenolol (100 mg) or nitrendipine (20 mg) once a day. After 8 days of treatment, only lisinopril was found to increase arterial compliance. Subsequently, we compared arterial diameter- and distensibility-pressure curves from newly diagnosed and untreated hypertensive patients with those of matched normotensive control patients. Diameter-pressure curves did not differ significantly between the groups and distensibility was not reduced. Similar findings were later obtained in an animal model, when mechanical properties of carotid arteries were compared between spontaneously hypertensive rats and normotensive counterparts (Wistar-Kyoto rats). These results, although interesting by providing noninvasive information on the elastic response of the wall, call for further development of the technique to be able to measure arterial wall thickness. Stress-strain relationship could ultimately be established to thoroughly characterize physical properties of blood vessel walls.

Quantitative comparison of simulations of seismic wave propagation in heterogeneous poro-elastic media involving fluid-solid interfaces and in equivalent visco-elastic solids

There is increasing evidence to suggest that the presence of mesoscopic heterogeneities constitutes the predominant attenuation mechanism at seismic frequencies. As a consequence, centimeter-scale perturbations of the subsurface physical properties should be taken into account for seismic modeling whenever detailed and accurate responses of the target structures are desired. This is, however, computationally prohibitive since extremely small grid spacings would be necessary. A convenient way to circumvent this problem is to use an upscaling procedure to replace the heterogeneous porous media by equivalent visco-elastic solids. In this work, we solve Biot's equations of motion to perform numerical simulations of seismic wave propagation through porous media containing mesoscopic heterogeneities. We then use an upscaling procedure to replace the heterogeneous poro-elastic regions by homogeneous equivalent visco-elastic solids and repeat the simulations using visco-elastic equations of motion. We find that, despite the equivalent attenuation behavior of the heterogeneous poro-elastic medium and the equivalent visco-elastic solid, the seismograms may differ due to diverging boundary conditions at fluid-solid interfaces, where there exist additional options for the poro-elastic case. In particular, we observe that the seismograms agree for closed-pore boundary conditions, but differ significantly for open-pore boundary conditions. This is an interesting result, which has potentially important implications for wave-equation-based algorithms in exploration geophysics involving fluid-solid interfaces, such as, for example, wave field decomposition.

Shot noise anomalies in elastic nondegenerate diffusive conductors

We present a theoretical investigation of shot-noise properties in nondegenerate elastic diffusive conductors. Both Monte Carlo simulations and analytical approaches are used. Two interesting phenomena are found: (i) the display of enhanced shot noise for given energy dependences of the scattering time, and (ii) the recovery of full shot noise for asymptotic high applied bias. The first phenomenon is associated with the onset of negative differential conductivity in energy space that drives the system towards a dynamical electrical instability in excellent agreement with analytical predictions. The enhancement is found to be strongly amplified when the dimensionality in momentum space is lowered from three to two dimensions. The second phenomenon is due to the suppression of the effects of long-range Coulomb correlations that takes place when the transit time becomes the shortest time scale in the system, and is common to both elastic and inelastic nondegenerate diffusive conductors. These phenomena shed different light in the understanding of the anomalous behavior of shot noise in mesoscopic conductors, which is a signature of correlations among different current pulses.