On the Equivalence of the Riemann-Liouville and the Caputo Fractional Order Derivatives in Modeling of Linear Viscoelastic Materials

Mathematics Subject Classification: 26A33

Modelling of non-linear viscoelastic tissues with bar elements

In this work we develop a viscoelastic bar element that can handle multiple rheo- logical laws with non-linear elastic and non-linear viscous material models. The bar element is built by joining in series an elastic and viscous bar, constraining the middle node position to the bar axis with a reduction method, and stati- cally condensing the internal degrees of freedom. We apply the methodology to the modelling of reversible softening with sti ness recovery both in 2D and 3D, a phenomenology also experimentally observed during stretching cycles on epithelial lung cell monolayers.

Polydispersity index from linear viscoelastic data: Unimodal and bimodal linear polymer melts

This article describes a method for determining the polydispersity index Ip2=Mz/Mw of the molecular weight distribution (MWD) of linear polymeric materials from linear viscoelastic data. The method uses the Mellin transform of the relaxation modulus of a simple molecular rheological model. One of the main features of this technique is that it enables interesting MWD information to be obtained directly from dynamic shear experiments. It is not necessary to achieve the relaxation spectrum, so the ill-posed problem is avoided. Furthermore, a determinate shape of the continuous MWD does not have to be assumed in order to obtain the polydispersity index. The technique has been developed to deal with entangled linear polymers, whatever the form of the MWD is. The rheological information required to obtain the polydispersity index is the storage G′(ω) and loss G″(ω) moduli, extending from the terminal zone to the plateau region. The method provides a good agreement between the proposed theoretical approach and the experimental polydispersity indices of several linear polymers for a wide range of average molecular weights and polydispersity indices. It is also applicable to binary blends.

Experimental analysis of active-passive vibration control using viscoelastic materials and extension and shear piezoelectric actuators

Hybrid active-passive damping treatments combine the reliability, low cost and robustness of viscoelastic damping treatments and the high-performance, modal selective and adaptive piezoelectric active control. Numerous hybrid damping treatments have been reported in the literature. They differ mainly by the relative positions of viscoelastic treatments, sensors and piezoelectric actuators. In this work we present an experimental analysis of three active-passive damping design configurations applied to a cantilever beam. In particular, two design configurations based on the extension mode of piezoelectric actuators combined with viscoelastic constrained layer damping treatments and one design configuration with shear piezoelectric actuators embedded in a sandwich beam with viscoelastic core are analyzed. For comparison purposes, a purely active design configuration with an extension piezoelectric actuator bonded to an elastic beam is also analyzed. The active-passive damping performance of the four design configurations is compared. Results show that active-passive design configurations provide more reliable and wider-range damping performance than the purely active configuration.

Superensembles of linear viscoelastic models of polymer melts

The idea of incorporating multiple models of linear rheology into a superensemble, to forge a consensus forecast from the individual model predictions, is investigated. The relative importance of the individual models in the so-called multimodel superensemble (MMSE) was inferred by evaluating their performance on a set of experimental training data, via nonlinear regression. The predictive ability of the MMSE model was tested by comparing its predictions on test data that were similar (in-sample) and dissimilar (out-of-sample) to the training data used in the calibration. For the in-sample forecasts, we found that the MMSE model easily outperformed the best constituent model. The presence of good individual models greatly enhanced the MMSE forecast, while the presence of some bad models in the superensemble also improved the MMSE forecast modestly. While the performance of the MMSE model on the out-of-sample training data was not as spectacular, it demonstrated the robustness of this approach.

Genetic Algorithm for the Determination of Linear Viscoelastic Relaxation Spectrum from Experimental Data

Conventional procedures employed in the modeling of viscoelastic properties of polymer rely on the determination of the polymer`s discrete relaxation spectrum from experimentally obtained data. In the past decades, several analytical regression techniques have been proposed to determine an explicit equation which describes the measured spectra. With a diverse approach, the procedure herein introduced constitutes a simulation-based computational optimization technique based on non-deterministic search method arisen from the field of evolutionary computation. Instead of comparing numerical results, this purpose of this paper is to highlight some Subtle differences between both strategies and focus on what properties of the exploited technique emerge as new possibilities for the field, In oder to illustrate this, essayed cases show how the employed technique can outperform conventional approaches in terms of fitting quality. Moreover, in some instances, it produces equivalent results With much fewer fitting parameters, which is convenient for computational simulation applications. I-lie problem formulation and the rationale of the highlighted method are herein discussed and constitute the main intended contribution. (C) 2009 Wiley Periodicals, Inc. J Appl Polym Sci 113: 122-135, 2009

Structure characterization of molecular complexes for non-linear optical materials. II. 1 : 1 complexes of 4-methyl-morpholine-N-oxide (1) and 3-picoline-N-oxide (2) with 2,4,6-trinitrophenol, studied by X-ray, AM1 and DFT calculations

(1) C6H2N3O7- center dot C5H12NO2+, Mr = 346.26, P2(1)/c, a = 7.2356(6), b = 10.5765(9), c = 19.593(2) angstrom, 3 beta=95.101(6)degrees, V = 1493.5(2) angstrom(3), Z = 4, R-1 = 0.0414; (2) C6H2N3O7- center dot C6H8NO+, Mr = 38.24, P2(1)/n, a = 7.8713(5), b = 6.1979(7), c = 28.697(3) angstrom, beta = 90.028(7)degrees, V = 1400.0(2) angstrom(3), Z = 4, R-1 = 0.0416. The packing units in both compounds consist of hydrogen bonded cation-anion pairs. The (hyper)polarizabilities have been calculated for the crystallographic and optimized molecules, by AM1 and at the DFT/B3LYP(6-31G**) level.

NUMERICAL AND ANALYTICAL VERIFICATION OF A MULTISCALE COMPUTATIONAL MODEL FOR IMPACT PROBLEMS IN HETEROGENEOUS VISCOELASTIC MATERIALS WITH EVOLVING DAMAGE

Composites are engineered materials that take advantage of the particular properties of each of its two or more constituents. They are designed to be stronger, lighter and to last longer which can lead to the creation of safer protection gear, more fuel efficient transportation methods and more affordable materials, among other examples. This thesis proposes a numerical and analytical verification of an in-house developed multiscale model for predicting the mechanical behavior of composite materials with various configurations subjected to impact loading. This verification is done by comparing the results obtained with analytical and numerical solutions with the results found when using the model. The model takes into account the heterogeneity of the materials that can only be noticed at smaller length scales, based on the fundamental structural properties of each of the composite’s constituents. This model can potentially reduce or eliminate the need of costly and time consuming experiments that are necessary for material characterization since it relies strictly upon the fundamental structural properties of each of the composite’s constituents. The results from simulations using the multiscale model were compared against results from direct simulations using over-killed meshes, which considered all heterogeneities explicitly in the global scale, indicating that the model is an accurate and fast tool to model composites under impact loads. Advisor: David H. Allen

Operational Modal Analysis on laminated glass beams

Laminated glass is a sandwich element consisting of two or more glass sheets, with one or more interlayers of polyvinyl butyral (PVB). The dynamic response of laminated glass beams and plates can be predicted using analytical or numerical models in which the glass and the PVB are usually modelled as linear-elastic and linear viscoelastic materials, respectively. In this work the dynamic behavior of laminated glass beams are predicted using a finite element model and the analytical model of Ross-Kerwin-Ungar. The numerical and analytical results are compared with those obtained by operational modal analysis performed at different temperatures.

Re-wetting potential of limestone filter cakes under mechanical decompression

During pressure filtration, the filter cake can experience deformation. When the filter cake exhibits elastic rheological behaviour, it expands while the pressure is released and, if adjacent liquid is present, re-wetting may occur. Such an expanding filter cake can pick up liquid already removed, worsening solid-liquid separation performances. Undesired phenomena such as filter cake re-wetting run contrary to solid-liquid separation performances and, specifically, to the separation target of obtaining a high dry solid content. At the industrial level, even a small quantity of liquid that can be additionally removed is beneficial. Re-wetting phenomenon is investigated using two different limestone materials and different filters and felts, respectively. Water storage capacity of filters and felts and elastic properties of filters, felts, and filter cakes are investigated. The elastic tests performed show that the filters and the felts are non-linear viscoelastic materials and can have a potential for re-wetting, while the limestone filter cakes are not showing measurable elastic deformation under decompression. However, during a filtration test an additional felt layer placed under the filter seems to result in a slightly higher cake dryness, i.e., an increase of the cake dryness by around 1% only. This difference may not be attributed to the re-wetting effect.

Viscoelastic properties of high pressure and heat induced tofu gels

Tofu gels were rheologically examined to determine their storage or elastic (G′) and loss or viscous (G″) moduli as a function of frequency within their linear viscoelastic limits. The tofu gels were made using either glucono-δ-lactone (GDL) or calcium sulphate (CaSO4·2H2O), followed by either heat treatment (heated soymilk at 97 °C prior to coagulation and subsequently held at 70 °C for 60 min, HT) or high pressure treatment (400 MPa at 20 °C for 10 min, HP). The overall moduli values of the GDL gels and CaSO4·2H2O gels of both physical treatments were similar, each gave frequency profiles expected for weak viscoelastic materials. However, although both temperature and high pressure treatments could be used to produce tofu gels, the final products were not the same. Pressure formed gels, despite having a higher overall “consistency” (increasing values of their moduli), had a proportionately higher contribution from the loss modulus (increased tan δ). Differences could also be observed using confocal scanning laser microscopy. While such treatment may give rise to differing systems/structures, with new or modified organoleptic properties, the more “open” structures obtained by pressure treatment may well cause processing difficulties if subsequent reworking or moulding is required.

Viscoelastic properties of high pressure and heat induced tofu gels

Tofu gels were rheologically examined to determine their storage or elastic (G') and loss or viscous (G '') moduli as a function of frequency within their linear viscoelastic limits. The tofu gels were made using either glucono-delta-lactone (GDL) or calcium sulphate (CaSO4 center dot 2H(2)O), followed by either heat treatment (heated soymilk at >= 97 degrees C prior to coagulation and subsequently held at 70 degrees C for 60 min, HT) or high pressure treatment (400 MPa at 20 degrees C for 10 min, HP). The overall moduli values of the GDL gels and CaSO4 center dot 2H(2)O gels of both physical treatments were similar, each gave frequency profiles expected for weak viscoelastic materials. However, although both temperature and high pressure treatments could be used to produce tofu gels, the final products were not the same. Pressure formed gels, despite having a higher overall "consistency" (increasing values of their moduli), had a proportionately higher contribution from the loss modulus (increased tan delta). Differences could also be observed using confocal scanning laser microscopy. While such treatment may give rise to differing systems/structures, with new or modified organoleptic properties, the more "open" structures obtained by pressure treatment may well cause processing difficulties if subsequent reworking or moulding is required. (c) 2007 Elsevier Ltd. All rights reserved.

Dynamic analysis for different finite element models of beams with viscoelastic damping layer

Many new viscoelastic materials have been developed recently to help improve noise and vibration levels in mechanical structures for applications in automobile and aeronautical industry. The viscoelastic layer treatment applied to solid metal structures modifies two main properties which are related to the mass distribution and the damping mechanism. The other property controlling the dynamics of a mechanical system is the stiffness that does not change much with the viscoelastic material. The model of such system is usually complex, because the viscoelastic material can exhibit nonlinear behavior, in contrast with the many available tools for linear dynamics. In this work, the dynamic behavior of sandwich beam is modeled by finite element method using different element types which are then compared with experimental results developed in the laboratory for various beams with different viscoelastic layer materials. The finite element model is them updated to help understand the effects in the damping for various natural frequencies and the trade-off between attenuation and the mass add to the structure.

Materials properties from electron density distributions: Molecular magnetism and linear optical properties

The general goal of this thesis is correlating observable properties of organic and metal-organic materials with their ground-state electron density distribution. In a long-term view, we expect to develop empirical or semi-empirical approaches to predict materials properties from the electron density of their building blocks, thus allowing to rationally engineering molecular materials from their constituent subunits, such as their functional groups. In particular, we have focused on linear optical properties of naturally occurring amino acids and their organic and metal-organic derivatives, and on magnetic properties of metal-organic frameworks. For analysing the optical properties and the magnetic behaviour of the molecular or sub-molecular building blocks in materials, we mostly used the more traditional QTAIM partitioning scheme of the molecular or crystalline electron densities, however, we have also investigated a new approach, namely, X-ray Constrained Extremely Localized Molecular Orbitals (XC-ELMO), that can be used in future to extracted the electron densities of crystal subunits. With the purpose of rationally engineering linear optical materials, we have calculated atomic and functional group polarizabilities of amino acid molecules, their hydrogen-bonded aggregates and their metal-organic frameworks. This has enabled the identification of the most efficient functional groups, able to build-up larger electric susceptibilities in crystals, as well as the quantification of the role played by intermolecular interactions and coordinative bonds on modifying the polarizability of the isolated building blocks. Furthermore, we analysed the dependence of the polarizabilities on the one-electron basis set and the many-electron Hamiltonian. This is useful for selecting the most efficient level of theory to estimate susceptibilities of molecular-based materials. With the purpose of rationally design molecular magnetic materials, we have investigated the electron density distributions and the magnetism of two copper(II) pyrazine nitrate metal-organic polymers. High-resolution X-ray diffraction and DFT calculations were used to characterize the magnetic exchange pathways and to establish relationships between the electron densities and the exchange-coupling constants. Moreover, molecular orbital and spin-density analyses were employed to understand the role of different magnetic exchange mechanisms in determining the bulk magnetic behaviour of these materials. As anticipated, we have finally investigated a modified version of the X-ray constrained wavefunction technique, XC-ELMOs, that is not only a useful tool for determination and analysis of experimental electron densities, but also enables one to derive transferable molecular orbitals strictly localized on atoms, bonds or functional groups. In future, we expect to use XC-ELMOs to predict materials properties of large systems, currently challenging to calculate from first-principles, such as macromolecules or polymers. Here, we point out advantages, needs and pitfalls of the technique. This work fulfils, at least partially, the prerequisites to understand materials properties of organic and metal-organic materials from the perspective of the electron density distribution of their building blocks. Empirical or semi-empirical evaluation of optical or magnetic properties from a preconceived assembling of building blocks could be extremely important for rationally design new materials, a field where accurate but expensive first-principles calculations are generally not used. This research could impact the community in the fields of crystal engineering, supramolecular chemistry and, of course, electron density analysis.