Comparison between analytical and 3D finite element solutions for borehole stresses in anisotropic elastic rock

We present a rigorous validation of the analytical Amadei solution for the stress concentration around an arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients b11 and b55 are not equal. It is shown from theoretical considerations and published experimental data that the b11 and b55 are not equal for realistic rocks. Second, we develop a 3D finite element elastic model within a hybrid analytical–numerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic, transverse isotropic and orthorhombic symmetries. It is concluded that the analytical Amadei solution is valid with no restriction on the borehole orientation or the symmetry of the elastic anisotropy.

Higher-order non-linear analysis of steel structures. Part I : elastic second-order formulation

This paper presents a higher-order beam-column formulation that can capture the geometrically non-linear behaviour of steel framed structures which contain a multiplicity of slender members. Despite advances in computational frame software, analyses of large frames can still be problematic from a numerical standpoint and so the intent of the paper is to fulfil a need for versatile, reliable and efficient non-linear analysis of general steel framed structures with very many members. Following a comprehensive review of numerical frame analysis techniques, a fourth-order element is derived and implemented in an updated Lagrangian formulation, and it is able to predict flexural buckling, snap-through buckling and large displacement post-buckling behaviour of typical structures whose responses have been reported by independent researchers. The solutions are shown to be efficacious in terms of a balance of accuracy and computational expediency. The higher-order element forms a basis for augmenting the geometrically non-linear approach with material non-linearity through the refined plastic hinge methodology described in the companion paper.

A comparative study of the semi-elastic and fully-elastic MapReduce models

MapReduce is a computation model for processing large data sets in parallel on large clusters of machines, in a reliable, fault-tolerant manner. A MapReduce computation is broken down into a number of map tasks and reduce tasks, which are performed by so called mappers and reducers, respectively. The placement of the mappers and reducers on the machines directly affects the performance and cost of the MapReduce computation in cloud computing. From the computational point of view, the mappers/reducers placement problem is a generation of the classical bin packing problem, which is NP-complete. Thus, in this paper we propose a new heuristic algorithm for the mappers/reducers placement problem in cloud computing and evaluate it by comparing with other several heuristics on solution quality and computation time by solving a set of test problems with various characteristics. The computational results show that our heuristic algorithm is much more efficient than the other heuristics and it can obtain a better solution in a reasonable time. Furthermore, we verify the effectiveness of our heuristic algorithm by comparing the mapper/reducer placement for a benchmark problem generated by our heuristic algorithm with a conventional mapper/reducer placement which puts a fixed number of mapper/reducer on each machine. The comparison results show that the computation using our mapper/reducer placement is much cheaper than the computation using the conventional placement while still satisfying the computation deadline.

Comparison of mechanical and ultrasound elastic modulus of ovine tibial cortical bone

Finite element models of bones can be created by deriving geometry from anx-ray CT scan. Material properties such as the elastic modulus can then be applied using either a single or set of homogeneous values, or individual elements can have local values mapped onto them. Values for the elastic modulus can be derived from the CT density values using an elasticityversus density relationship. Many elasticity–density relationships have been reported in the literature for human bone. However, while ovine in vivo models are common in orthopaedic research, no work has been done to date on creating FE models of ovine bones. To create these models and apply relevant material properties, an ovine elasticity-density relationship needs to be determined. Using fresh frozen ovine tibias the apparent density of regions of interest was determined from a clinical CT scan. The bones were the sectioned into cuboid samples of cortical bone from the regions of interest. Ultrasound was used to determine the elastic modulus in each of three directions – longitudinally, radially and tangentially. Samples then underwent traditional compression testing in each direction. The relationships between apparent density and both ultrasound, and compression modulus in each directionwere determined. Ultrasound testing was found to be a highly repeatable non-destructive method of calculating the elastic modulus, particularly suited to samples of this size. The elasticity-density relationships determined in the longitudinal direction were very similar between the compression and ultrasound data over the density range examined.A clear difference was seen in the elastic modulus between the longitudinal and transverse directions of the bone samples, and a transverse elasticity-density relationship is also reported.

On helical shear of a compressible elastic circular cylindrical tube

In this paper we examine the combined azimuthal and axial shear of a compressible isotropic elastic circular cylindrical tube of finite extent, otherwise referred to as helical shear (which is an isochoric deformation). The equilibrium equations are formulated in terms of the principal stretches, and explicit necessary and sufficient conditions on the strain-energy function for the material to support this deformation are obtained and compared with those obtained previously for this problem. Several classes of strain-energy functions are derived and in some general cases complete solutions of the equilibrium equations are obtained. Existing results are recovered as special cases and some new results for the strain-energy functions derived are determined and discussed.

On extension and torsion of a compressible elastic circular cylinder

In this paper we examine the combined extension and torsion of a compressible isotropic elastic cylinder of finite extent. The equilibrium equations are formulated in terms of the principal stretches and then applied to the special case of pure torsion superimposed on a uniform extension (an isochoric deformation). Explicit necessary and sufficient conditions on the strain-energy function for the material to support this deformation with vanishing traction on the lateral surfaces of the cylinder are obtained. Some strain-energy functions satisfying these conditions are considered, existing results are recovered as special cases and new results are obtained. We also point out how the strain-energy functions generated from the considered isochoric deformation considered (of a compressible material) can be used to generate energy functions and corresponding solutions for the incompressible theory.

Revisiting borehole stresses in anisotropic elastic media : comparison of analytical versus numerical solutions

We present a rigorous validation of the analyticalAmadei solution for the stress concentration around arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients β11 and β55 are not equal. It is shown from theoretical considerations and published experimental data that the β11 and β55 are not equal for realistic rocks. Second, we develop a 3D finite-element elastic model within a hybrid analyticalnumerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic and transverse isotropic symmetries. It is concluded that the analytical Amadei solution is valid with no restrictions on the borehole orientation or elastic anisotropy symmetry.

Effect of constitutive material models on seismic response of two-story reinforced concrete frame

This paper focuses on the finite element (FE) response sensitivity and reliability analyses considering smooth constitutive material models. A reinforced concrete frame is modeled for FE sensitivity analysis followed by direct differentiation method under both static and dynamic load cases. Later, the reliability analysis is performed to predict the seismic behavior of the frame. Displacement sensitivity discontinuities are observed along the pseudo-time axis using non-smooth concrete and reinforcing steel model under quasi-static loading. However, the smooth materials show continuity in response sensitivity at elastic to plastic transition points. The normalized sensitivity results are also used to measure the relative importance of the material parameters on the structural responses. In FE reliability analysis, the influence of smoothness behavior of reinforcing steel is carefully noticed. More efficient and reasonable reliability estimation can be achieved by using smooth material model compare with bilinear material constitutive model.

Prediction of shrinkage and porosity during drying: Considering both material properties and process conditions

Drying of food materials offers a significant increase in the shelf life of food materials, along with the modification of quality attributes due to simultaneous heat and mass transfer. Shrinkage and variations in porosity are the common micro and microstructural changes that take place during the drying of mostly the food materials. Although extensive research has been carried out on the prediction of shrinkage and porosity over the time of drying, no single model exists which consider both material properties and process condition in the same model. In this study, an attempt has been made to develop and validate shrinkage and porosity models of food materials during drying considering both process parameters and sample properties. The stored energy within the sample, elastic potential energy, glass transition temperature and physical properties of the sample such as initial porosity, particle density, bulk density and moisture content have been taken into consideration. Physical properties and validation have been made by using a universal testing machine ( Instron 2kN), a profilometer (Nanovea) and a pycnometer. Apart from these, COMSOL Multiphysics 4.4 has been used to solve heat and mass transfer physics. Results obtained from models of shrinkage and porosity is quite consistent with the experimental data. Successful implementation of these models would ensure the use of optimum energy in the course of drying and better quality retention of dried foods.

Diffraction Of Elastic Waves By Two Parallel Rigid Strips Embedded In An Infinite Orthotropic Medium

The elastodynamic response of a pair of parallel rigid strips embedded in an infinite orthotropic medium due to elastic waves incident normally on the strips has been investigated. The mixed boundary value problem has been solved by the Integral Equation method. The normal stress and the vertical displacement have been derived in closed form. Numerical values of stress intensity factors at inner and outer edges of the strips and vertical displacement at points in the plane of the strips for several orthotropic materials have been calculated and plotted graphically to show the effect of material orthotropy.

On the sensitivity of elastic waves due to structural damages:Time-frequency based indexing method

Time-frequency analysis of various simulated and experimental signals due to elastic wave scattering from damage are performed using wavelet transform (WT) and Hilbert-Huang transform (HHT) and their performances are compared in context of quantifying the damages. Spectral finite element method is employed for numerical simulation of wave scattering. An analytical study is carried out to study the effects of higher-order damage parameters on the reflected wave from a damage. Based on this study, error bounds are computed for the signals in the spectral and also on the time-frequency domains. It is shown how such an error bound can provide all estimate of error in the modelling of wave propagation in structure with damage. Measures of damage based on WT and HHT is derived to quantify the damage information hidden in the signal. The aim of this study is to obtain detailed insights into the problem of (1) identifying localised damages (2) dispersion of multifrequency non-stationary signals after they interact with various types of damage and (3) quantifying the damages. Sensitivity analysis of the signal due to scattered wave based on time-frequency representation helps to correlate the variation of damage index measures with respect to the damage parameters like damage size and material degradation factors.

Elastic analysis of pin joints in plates under some combined pin and plate loads

Joints are primary sources of weakness in structures. Pin joints are very common and are used where periodic disassembly of components is needed. A circular pin in a circular hole in an infinitely large plate is an abstraction of such a pin joint. A two-dimensional plane-stress analysis of such a configuration is carried out, here, subjected to pin-bearing and/or biaxial-plate loading. The pin is assumed to be rigid compared to the plate material. For pin load the reactive stresses at the edges of the infinite plate tend to zero though their integral over the external boundary equals to the pin load. The pin-hole interface is unbonded and so beyond some load levels the plate separates from the pin and the extent of separation is a non-linear function of load level. The problem is solved by inverse technique where the extent of contact is specified and the causative loads are evaluated directly. In the situations where combined load is acting the separation-contact zone specification generally needs two parameters (angles) to be specified. The present report deals with analysing such a situation in metallic (or isotropic) plates. Numerical results are provided for parametric representation and the methodology is demonstrated.