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.

A spectrum-driven damage Identification technique : application and validation through the numerical simulation of the Z24 Bridge

The present paper focuses on a damage identification method based on the use of the second order spectral properties of the nodal response processes. The explicit dependence on the frequency content of the outputs power spectral densities makes them suitable for damage detection and localization. The well-known case study of the Z24 Bridge in Switzerland is chosen to apply and further investigate this technique with the aim of validating its reliability. Numerical simulations of the dynamic response of the structure subjected to different types of excitation are carried out to assess the variability of the spectrum-driven method with respect to both type and position of the excitation sources. The simulated data obtained from random vibrations, impulse, ramp and shaking forces, allowed to build the power spectrum matrix from which the main eigenparameters of reference and damage scenarios are extracted. Afterwards, complex eigenvectors and real eigenvalues are properly weighed and combined and a damage index based on the difference between spectral modes is computed to pinpoint the damage. Finally, a group of vibration-based damage identification methods are selected from the literature to compare the results obtained and to evaluate the performance of the spectral index.

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.

Thermal conditions affecting heat transfer in FDM/FFE: a contribution towards the numerical modelling of the process

The performance of parts produced by Free Form Extrusion (FFE), an increasingly popular additive manufacturing technique, depends mainly on their dimensional accuracy, surface quality and mechanical performance. These attributes are strongly influenced by the evolution of the filament temperature and deformation during deposition and solidification. Consequently, the availability of adequate process modelling software would offer a powerful tool to support efficient process set-up and optimisation. This work examines the contribution to the overall heat transfer of various thermal phenomena developing during the manufacturing sequence, including convection and radiation with the environment, conduction with support and between adjacent filaments, radiation between adjacent filaments and convection with entrapped air. The magnitude of the mechanical deformation is also studied. Once this exercise is completed, it is possible to select the material properties, process variables and thermal phenomena that should be taken in for effective numerical modelling of FFE.

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

Slip flows of Newtonian and viscoelastic fluids in a 4:1 contraction

This work presents a numerical study of the 4:1 planar contraction flow of a viscoelastic fluid described by the simplified Phan-Thien–Tanner model under the influence of slip boundary conditions at the channel walls. The linear Navier slip law was considered with the dimensionless slip coefficient varying in the range ½0; 4500. The simulations were carried out for a small constant Reynolds number of 0.04 and Deborah numbers (De) varying between 0 and 5. Convergence could not be achieved for higher values of the Deborah number, especially for large values of the slip coefficient, due to the large stress gradients near the singularity of the reentrant corner. Increasing the slip coefficient leads to the formation of two vortices, a corner and a lip vortex. The lip vortex grows with increasing slip until it absorbs the corner vortex, creating a single large vortex that continues to increase in size and intensity. In the range De = 3–5 no lip vortex was formed. The flow is characterized in detail for De ¼ 1 as function of the slip coefficient, while for the remaining De only the main features are shown for specific values of the slip coefficient.

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.

Comparison of different numerical methods for the solution of the time-fractional reaction-diffusion equation with variable diffusion coefficient

In this work we perform a comparison of two different numerical schemes for the solution of the time-fractional diffusion equation with variable diffusion coefficient and a nonlinear source term. The two methods are the implicit numerical scheme presented in [M.L. Morgado, M. Rebelo, Numerical approximation of distributed order reaction- diffusion equations, Journal of Computational and Applied Mathematics 275 (2015) 216-227] that is adapted to our type of equation, and a colocation method where Chebyshev polynomials are used to reduce the fractional differential equation to a system of ordinary differential equations

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.

Metallic glass/PVDF magnetoelectric energy harvester working up to the radiofrequency range

Companies and researchers involved in developing miniaturized electronic devices face the basic problem of the needed batteries size, finite life of time and environmental pollution caused by their final deposition. The current trends to overcome this situation point towards Energy Harvesting technology. These harvesters (or scavengers) store the energy from sources present in the ambient (as wind, solar, electromagnetic, etc) and are costless for us. Piezoelectric devices are the ones that show a higher power density, and materials as ceramic PZT or polymeric PVDF have already demonstrated their ability to act as such energy harvester elements. Combinations between piezoelectric and electromagnetic mechanism have been also extensively investigated. Nevertheless, the power generated by these combinations is limited under the application of small magnetic fields, reducing the performance of the energy harvester [1]. In the last years the appearance of magnetoelectric (ME) devices, in which the piezoelectric deformation is driven by the magnetostrictive element, enables to extract the energy of very small electromagnetic signals through the generated magnetoelectric voltage at the piezoelectric element. However, very little work has been done testing PVDF polymer as piezoelectric constituent of the ME energy harvester device, and only to be proposed as a possibility of application [2]. Among the advantages of using piezopolymers for vibrational energy harvesting we can remember that they are ductile, resilient to shock, deformable and lightweight. In this work we demonstrate the feasibility of using magnetostrictive Fe-rich magnetic amorphous alloys/piezoelectric PVDF sandwich-type laminated ME devices as energy harvesters. A very simple experimental set-up will show how these laminates can extract energy, in amounts of μW, from an external AC field.

Anderson localization of light in disordered superlattices containing graphene layers

We theoretically investigate light propagation and Anderson localization in one-dimensional disordered superlattices composed of dielectric stacks with graphene sheets in between. Disorder is introduced either on graphene material parameters ({\it e.g.} Fermi energy) or on the widths of the dielectric stacks. We derive an analytic expression for the localization length $\xi$, and compare it to numerical simulations using transfer matrix technique; a very good agreement is found. We demonstrate that the presence of graphene may strongly attenuate the anomalously delocalised Breswter modes, and is at the origin of a periodic dependence of $\xi$ on frequency, in contrast to the usual asymptotic decay, $\xi \propto \omega^{-2}$. By unveiling the effects of graphene on Anderson localization of light, we pave the way for new applications of graphene-based, disordered photonic devices in the THz spectral range.

Numerical simulation of the reflow soldering process

Dissertação de mestrado integrado em Mechanical Engineering

A new combination in Peristethium (Loranthaceae) expands the genus' range into the Amazon-Cerrado ecotone

The genus Peristethium, characterized by determinate inflorescences protected by deciduous bracts, occurs in the northwest of South America, as well as Costa Rica and Panama. The main objective of this paper was to transfer one species to what we believe is its correct generic placement in Peristethium, that likewise implies in a shift of the genus' distribution beyond the Amazon. A new combination, Peristethium reticulatum, is proposed, based on Struthanthus reticulatus, described from Tocantins in 1980. The sexual dimorphism of the inflorescences of P. reticulatum (sessile male flowers and pedicellate female flowers) associated with male inflorescences that are inserted at leafless nodes are unique within the genus. The male flowers have dimorphic stamens, well-developed anthers and a pistiloid, whilst female flowers have robust styles and stigmas, and much reduced staminodes. Peristethium reticulatum and P. polystachyum occurs in the Amazon regions of Brazil, with the former recorded also in the ecotone with the Central Brazilian savannas (Cerrados).

Numerical and experimental analysis of full scale arches reinforced with GFRP materials

In this contribution, original limit analysis numerical results are presented dealing with some reinforced masonry arches tested at the University of Minho-UMinho, PT. Twelve in-scale circular masonry arches were considered, reinforced in various ways at the intrados or at the extrados. GFRP reinforcements were applied either on undamaged or on previously damaged elements, in order to assess the role of external reinforcements even in repairing interventions. The experimental results were critically discussed at the light of limit analysis predictions, based on a 3D FE heterogeneous upper bound approach. Satisfactory agreement was found between experimental evidences and the numerical results, in terms of failure mechanisms and peak load.