Effect of environmental ageing on the numerical response of FRP-strengthened masonry walls

Recent durability studies have shown the susceptibility of bond in fiber-reinforced polymer (FRP) strengthened masonry components to hygrothermal exposures. However, it is not clear how this local material degradation affects the global behavior of FRP-strengthened masonry structures. This study addresses this issue by numerically investigating the nonlinear behavior of FRP-masonry walls after aging in two different environmental conditions. A numerical modeling strategy is adopted and validated with existing experimental tests on FRP-strengthened masonry panels. The model, once validated, is used for modeling of four hypothetical FRP-strengthened masonry walls with different boundary conditions, strengthening schemes, and reinforcement ratios. The nonlinear behavior of the walls is then simulated before and after aging in two different environmental conditions. The degradation data are taken from previous accelerated aging tests. The changes in the failure mode and nonlinear response of the walls after aging are presented and discussed.

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.

Numerical approaches for the analysis of timber frame walls

A numerical approach to simulate the behaviour of timber shear walls under both static and dynamic loading is proposed. Because the behaviour of timber shear walls hinges on the behaviour of the nail connections, the force-displacement behaviour of sheathing-to-framing nail connections are first determined and then used to define the hysteretic properties of finite elements representing these connections. The model nails are subsequently implemented into model walls. The model walls are verified using experimental results for both monotonic and cyclic loading. It is demonstrated that the complex hysteretic behaviour of timber shear walls can be reasonably represented using model shear walls in which nonlinear material failure is concentrated only at the sheathing-to-framing nail connections.

Analytical and numerical study of the electro-osmotic annular flow of viscoelastic fluids

In this work we present semi-analytical solutions for the electro-osmotic annular flow of viscoelastic fluids modeled by the Linear and Exponential PTT models. The viscoelastic fluid flows in the axial direction between two concentric cylinders under the combined influences of electrokinetic and pressure forcings. The analysis invokes the Debye-Hückel approximation and includes the limit case of pure electro-osmotic flow. The solution is valid for both no slip and slip velocity at the walls and the chosen slip boundary condition is the linear Navier slip velocity model. The combined effects of fluid rheology, electro-osmotic and pressure gradient forcings on the fluid velocity distribution are also discussed.

Annular flow of viscoelastic fluids: Analytical and numerical solutions

This work provides analytical and numerical solutions for the linear, quadratic and exponential Phan–Thien–Tanner (PTT) viscoelastic models, for axial and helical annular fully-developed flows under no slip and slip boundary conditions, the latter given by the linear and nonlinear Navier slip laws. The rheology of the three PTT model functions is discussed together with the influence of the slip velocity upon the flow velocity and stress fields. For the linear PTT model, full analytical solutions for the inverse problem (unknown velocity) are devised for the linear Navier slip law and two different slip exponents. For the linear PTT model with other values of the slip exponent and for the quadratic PTT model, the polynomial equation for the radial location (β) of the null shear stress must be solved numerically. For both models, the solution of the direct problem is given by an iterative procedure involving three nonlinear equations, one for β, other for the pressure gradient and another for the torque per unit length. For the exponential PTT model we devise a numerical procedure that can easily compute the numerical solution of the pure axial flow problem

Fractional pennes' bioheat equation: theoretical and numerical studies

In this work we provide a new mathematical model for the Pennes’ bioheat equation, assuming a fractional time derivative of single order. Alternative versions of the bioheat equation are studied and discussed, to take into account the temperature-dependent variability in the tissue perfusion, and both finite and infinite speed of heat propagation. The proposed bioheat model is solved numerically using an implicit finite difference scheme that we prove to be convergent and stable. The numerical method proposed can be applied to general reaction diffusion equations, with a variable diffusion coefficient. The results obtained with the single order fractional model, are compared with the original models that use classical derivatives.

Development of a new 3D OpenFOAM® solver to model the cooling stage in profile extrusion

[Excerpt] The advantages resulting from the use of numerical modelling tools to support the design of processing equipment are almost consensual. The design of calibration systems in profile extrusion is not an exception . H owever , the complex geome tries and heat exchange phenomena involved in this process require the use of numerical solvers able to model the heat exchange in more than one domain ( calibrator and polymer), the compatibilization of the heat transfer at the profile - calibrator interface and with the ability to deal with complex geometries. The combination of all these features is usually hard to find in commercial software. Moreover , the dimension of the meshes required to ob tain accurate results, result in computational times prohibitive for industrial application. (...)

Improvign mixing processes with the help of numerical tools

In several industrial applications, highly complex behaviour materials are used together with intricate mixing processes, which difficult the achievement of the desired properties for the produced materials. This is the case of the well-known dispersion of nano-sized fillers in a melt polymer matrix, used to improve the nanocomposite mechanical and/or electrical properties. This mixing is usually performed in twin-screw extruders, that promote complex flow patterns, and, since an in loco analysis of the material evolution and mixing is difficult to perform, numerical tools can be very useful to predict the evolution and behaviour of the material. This work presents a numerical based study to improve the understanding of mixing processes. Initial numerical studies were performed with generalized Newtonian fluids, but, due to the null relaxation time that characterize this type of fluids, the assumption of viscoelastic behavior was required. Therefore, the polymer melt was rheologically characterized, and, a six mode Phan-Thien-Tanner and Giesekus models were used to fit the rheological data. These viscoelastic rheological models were used to model the process. The conclusions obtained in this work provide additional and useful data to correlate the type and intensity of the deformation history promoted to the polymer nanocomposite and the quality of the mixing obtained.

Numerical modeling and assessment of the thermal environment in an industrial premise of automotive components

Dissertação de mestrado em Engenharia Industrial

Implementation and validation of a strain rate dependent anisotropic continuum model for masonry

A newly developed strain rate dependent anisotropic continuum model is proposed for impact and blast applications in masonry. The present model adopted the usual approach of considering different yield criteria in tension and compression. The analysis of unreinforced block work masonry walls subjected to impact is carried out to validate the capability of the model. Comparison of the numerical predictions and test data revealed good agreement. Next, a parametric study is conducted to evaluate the influence of the tensile strengths along the three orthogonal directions and of the wall thickness on the global behavior of masonry walls.

Dynamic interface model for masonry walls subjected to high strain rate out-of-plane loads

The present study proposes a dynamic constitutive material interface model that includes non-associated flow rule and high strain rate effects, implemented in the finite element code ABAQUS as a user subroutine. First, the model capability is validated with numerical simulations of unreinforced block work masonry walls subjected to low velocity impact. The results obtained are compared with field test data and good agreement is found. Subsequently, a comprehensive parametric analysis is accomplished with different joint tensile strengths and cohesion, and wall thickness to evaluate the effect of the parameter variations on the impact response of masonry walls.

Experimental and numerical studies of hybrid PCM embedded in plastering mortar for enhanced thermal behaviour of buildings

This paper proposes a methodology for improvement of energy efficiency in buildings through the innovative simultaneous incorporation of three distinct phase change materials (here termed as hybrid PCM) in plastering mortars for façade walls. The thermal performance of a hybrid PCM mortar was experimentally evaluated by comparing the behaviour of a prototype test cell (including hybrid PCM plastering mortar) subjected to realistic daily temperature profiles, with the behaviour of a similar prototype test cell, in which no PCM was added. A numerical simulation model was employed (using ANSYS-FLUENT) to validate the capacity of simulating temperature evolution within the prototype containing hybrid PCM, as well as to understand the contribution of hybrid PCM to energy efficiency. Incorporation of hybrid PCM into plastering mortars was found to have the potential to significantly reduce heating/cooling temperature demands for maintaining the interior temperature within comfort levels when compared to normal mortars (without PCM), or even mortars comprising a single type of PCM.

Analytical and numerical investigation of two-phase flows

Multiphase flows, hyperbolic model, Godunov method, nozzle flow, nonstrictly hyperbolic

Double-Exponential Model of a Photovoltaic Cell with Temperature Dependence

The purpose of the research is the creation of mathematical models in MATLAB based on the double exponential model of the photovoltaic cell. The developed model allows for different physical and environmental parameters. An equivalent circuit of the model includes a photocurrent source, two diodes, and a series and parallel resistance. The paper presents the simulation results for each parameter. The simulation data are displayed graphically and numerical results are saved in a file.

Convergence of the mass-transport steepest descent scheme for the sub-critical Patlak-Keller-Segel model

Variational steepest descent approximation schemes for the modified Patlak-Keller-Segel equation with a logarithmic interaction kernel in any dimension are considered. We prove the convergence of the suitably interpolated in time implicit Euler scheme, defined in terms of the Euclidean Wasserstein distance, associated to this equation for sub-critical masses. As a consequence, we recover the recent result about the global in time existence of weak-solutions to the modified Patlak-Keller-Segel equation for the logarithmic interaction kernel in any dimension in the sub-critical case. Moreover, we show how this method performs numerically in one dimension. In this particular case, this numerical scheme corresponds to a standard implicit Euler method for the pseudo-inverse of the cumulative distribution function. We demonstrate its capabilities to reproduce easily without the need of mesh-refinement the blow-up of solutions for super-critical masses.