A vertex-based finite volume method applied to non-linear material problems in computational solid mechanics

A vertex-based finite volume (FV) method is presented for the computational solution of quasi-static solid mechanics problems involving material non-linearity and infinitesimal strains. The problems are analysed numerically with fully unstructured meshes that consist of a variety of two- and threedimensional element types. A detailed comparison between the vertex-based FV and the standard Galerkin FE methods is provided with regard to discretization, solution accuracy and computational efficiency. For some problem classes a direct equivalence of the two methods is demonstrated, both theoretically and numerically. However, for other problems some interesting advantages and disadvantages of the FV formulation over the Galerkin FE method are highlighted.

A Simplified Method for Ductility Calculation in RC Jacketed Columns

Reinforced concrete (RC) jacketing is a common method to retrofit existing columns with poor structural performance. It can be applied in two different ways: if the continuity of the jacket is ensured, the axial load of the column can be transferred to the jacket, which will be directly loaded; conversely, if no continuity is provided, the jacket induces only confinement action. In both cases the strength and ductility evaluation is rather complex, due to the different physical phenomena included, such as confinement, composite action core-jacket, preload, buckling of longitudinal bars.
Although different theoretical studies have been carried out to calculate the confinement effects, a practical approach to evaluate the flexural capacity and ductility is still missing. The calculation of these quantities is often related to the use of commercial computer programs, taking advantage of numerical methods such as fiber method or finite element method.
This paper presents a simplified approach to calculate the flexural strength and ductility of square RC jacketed sections subjected to axial load and bending moment. In particular the proposed approach is based on the calibration of the stress-block parameters including the confinement effect. Equilibrium equations are determined and buckling of longitudinal bars is modeled with a suitable stress-strain law. Moment-curvature curves are derived with simple calculations. Finally, comparisons are made with numerical analyses carried out with the code OpenSees and with experimental data available in the literature, showing good agreement.

Strength and ductility of RC jacketed columns: a simplified analytical method

Reinforced concrete (RC) jacketing is a common method for retrofitting existing columns with poor structural performance. It can be applied in two different ways: if the continuity of the jacket is ensured, the axial load of the column can be transferred to the jacket, which will be directly loaded; conversely, if no continuity is provided, the jacket will induce only confinement action. In both cases the strength and ductility evaluation is rather complex, due to the different physical phenomena included, such as confinement, core-jacket composite action, preload and buckling of longitudinal bars.
Although different theoretical studies have been carried out to calculate the confinement effects, a practical approach to evaluate the flexural capacity and ductility is still missing. The calculation of these quantities is often related to the use of commercial software, taking advantage of numerical methods such as fibre method or finite element method.
This paper presents a simplified approach to calculate the flexural strength and ductility of square RC jacketed sections subjected to axial load and bending moment. In particular the proposed approach is based on the calibration of the stress-block parameters including the confinement effect. Equilibrium equations are determined and buckling of longitudinal bars is modelled with a suitable stress-strain law. Moment-curvature curves are derived with simple calculations. Finally, comparisons are made with numerical analyses carried out with the code OpenSees and with experimental data available in the literature, showing good agreement.

Developing a numerical inverse-theory-based extraction of orientation-dependent relaxation rates from partially- relaxed spectra

Second-rank tensor interactions, such as quadrupolar interactions between the spin- 1 deuterium nuclei and the electric field gradients created by chemical bonds, are affected by rapid random molecular motions that modulate the orientation of the molecule with respect to the external magnetic field. In biological and model membrane systems, where a distribution of dynamically averaged anisotropies (quadrupolar splittings, chemical shift anisotropies, etc.) is present and where, in addition, various parts of the sample may undergo a partial magnetic alignment, the numerical analysis of the resulting Nuclear Magnetic Resonance (NMR) spectra is a mathematically ill-posed problem. However, numerical methods (de-Pakeing, Tikhonov regularization) exist that allow for a simultaneous determination of both the anisotropy and orientational distributions. An additional complication arises when relaxation is taken into account. This work presents a method of obtaining the orientation dependence of the relaxation rates that can be used for the analysis of the molecular motions on a broad range of time scales. An arbitrary set of exponential decay rates is described by a three-term truncated Legendre polynomial expansion in the orientation dependence, as appropriate for a second-rank tensor interaction, and a linear approximation to the individual decay rates is made. Thus a severe numerical instability caused by the presence of noise in the experimental data is avoided. At the same time, enough flexibility in the inversion algorithm is retained to achieve a meaningful mapping from raw experimental data to a set of intermediate, model-free

A cell by cell anisotropic adaptive mesh ALE method

A cell by cell anisotropic adaptive mesh Arbitrary Lagrangian Eulerian (ALE) method for the solution of the Euler equations is described. An efficient approach to equipotential mesh relaxation on anisotropically refined meshes is developed. Results for two test problems are presented.