A three-dimensional hybrid smoothed finite element method (H-SFEM) for nonlinear solid mechanics problems

This paper presents a novel three-dimensional hybrid smoothed finite element method (H-SFEM) for solid mechanics problems. In 3D H-SFEM, the strain field is assumed to be the weighted average between compatible strains from the finite element method (FEM) and smoothed strains from the node-based smoothed FEM with a parameter α equipped into H-SFEM. By adjusting α, the upper and lower bound solutions in the strain energy norm and eigenfrequencies can always be obtained. The optimized α value in 3D H-SFEM using a tetrahedron mesh possesses a close-to-exact stiffness of the continuous system, and produces ultra-accurate solutions in terms of displacement, strain energy and eigenfrequencies in the linear and nonlinear problems. The novel domain-based selective scheme is proposed leading to a combined selective H-SFEM model that is immune from volumetric locking and hence works well for nearly incompressible materials. The proposed 3D H-SFEM is an innovative and unique numerical method with its distinct features, which has great potential in the successful application for solid mechanics problems.

Spatial Variability of the Depth of Weathered and Engineering Bedrock using Multichannel Analysis of Surface Wave Method

In this paper an attempt has been made to evaluate the spatial variability of the depth of weathered and engineering bedrock in Bangalore, south India using Multichannel Analysis of Surface Wave (MASW) survey. One-dimensional MASW survey has been carried out at 58 locations and shear-wave velocities are measured. Using velocity profiles, the depth of weathered rock and engineering rock surface levels has been determined. Based on the literature, shear-wave velocity of 330 ± 30 m/s for weathered rock or soft rock and 760 ± 60 m/s for engineering rock or hard rock has been considered. Depths corresponding to these velocity ranges are evaluated with respect to ground contour levels and top surface levels have been mapped with an interpolation technique using natural neighborhood. The depth of weathered rock varies from 1 m to about 21 m. In 58 testing locations, only 42 locations reached the depths which have a shear-wave velocity of more than 760 ± 60 m/s. The depth of engineering rock is evaluated from these data and it varies from 1 m to about 50 m. Further, these rock depths have been compared with a subsurface profile obtained from a two-dimensional (2-D) MASW survey at 20 locations and a few selected available bore logs from the deep geotechnical boreholes.

Denoising neural data with state-space smoothing: Method and application

Neural data are inevitably contaminated by noise. When such noisy data are subjected to statistical analysis, misleading conclusions can be reached. Here we attempt to address this problem by applying a state-space smoothing method, based on the combined use of the Kalman filter theory and the Expectation–Maximization algorithm, to denoise two datasets of local field potentials recorded from monkeys performing a visuomotor task. For the first dataset, it was found that the analysis of the high gamma band (60–90 Hz) neural activity in the prefrontal cortex is highly susceptible to the effect of noise, and denoising leads to markedly improved results that were physiologically interpretable. For the second dataset, Granger causality between primary motor and primary somatosensory cortices was not consistent across two monkeys and the effect of noise was suspected. After denoising, the discrepancy between the two subjects was significantly reduced.

A new method for the in vivo identification of mechanical properties in arteries from cine MRI images: Theoretical framework and validation

Quantifying the stiffness properties of soft tissues is essential for the diagnosis of many cardiovascular diseases such as atherosclerosis. In these pathologies it is widely agreed that the arterial wall stiffness is an indicator of vulnerability. The present paper focuses on the carotid artery and proposes a new inversion methodology for deriving the stiffness properties of the wall from cine-MRI (magnetic resonance imaging) data. We address this problem by setting-up a cost function defined as the distance between the modeled pixel signals and the measured ones. Minimizing this cost function yields the unknown stiffness properties of both the arterial wall and the surrounding tissues. The sensitivity of the identified properties to various sources of uncertainty is studied. Validation of the method is performed on a rubber phantom. The elastic modulus identified using the developed methodology lies within a mean error of 9.6%. It is then applied to two young healthy subjects as a proof of practical feasibility, with identified values of 625 kPa and 587 kPa for one of the carotid of each subject.

Stress analysis of carotid atheroma in a transient ischaemic attack patient using the MRI-based fluid-structure interaction method

Rupture of atherosclerotic plaque is a major cause of mortality. Plaque stress analysis, based on patient-specific multisequence in vivo MRI, can provide critical information for the understanding of plaque rupture and could eventually lead to plaque rupture prediction. However, the direct link between stress and plaque rupture is not fully understood. In the present study, the plaque from a patient who recently experienced a transient ischaemic attack (TIA) was studied using a fluid-structure interaction method to quantify stress distribution in the plaque region based on in vivo MR images. The results showed that wall shear stress is generally low in the artery with a slight increase at the plaque throat owing to minor luminal narrowing. The oscillatory shear index is much higher in the proximal part of the plaque. Both local wall stress concentrations and the relative stress variation distribution during a cardiac cycle indicate that the actual plaque rupture site is collocated with the highest rupture risk region in the studied patient.

In vivo non-invasive high resolution MR-based method for the determination of the elastic modulus of arterial vessels

The mechanical properties of arterial walls have long been recognized to play an essential role in the development and progression of cardiovascular disease (CVD). Early detection of variations in the elastic modulus of arteries would help in monitoring patients at high cardiovascular risk stratifying them according to risk. An in vivo, non-invasive, high resolution MR-phase-contrast based method for the estimation of the time-dependent elastic modulus of healthy arteries was developed, validated in vitro by means of a thin walled silicon rubber tube integrated into an existing MR-compatible flow simulator and used on healthy volunteers. A comparison of the elastic modulus of the silicon tube measured from the MRI-based technique with direct measurements confirmed the method's capability. The repeatability of the method was assessed. Viscoelastic and inertial effects characterizing the dynamic response of arteries in vivo emerged from the comparison of the pressure waveform and the area variation curve over a period. For all the volunteers who took part in the study the elastic modulus was found to be in the range 50-250 kPa, to increase during the rising part of the cycle, and to decrease with decreasing pressure during the downstroke of systole and subsequent diastole.

Direct precipitation of lead zirconate titanate by the hydrothermal method

Direct precipitation of fine powders of lead zirconate titanate (PZT) in the complete range of solid solution, is investigated under hydrothermal conditions, starting from lead oxide and titania/zirconia mixed gels. The perovskite phase is formed in the temperature range of 165 – 340°C. Sequence of the hydrothermal reactions is studied by identifying the intermediate phases. The initial formation of PbO: TiO2 solid solution is followed by the reaction of the same with the remaining mixed gels giving rise to X-ray amorphous PZT phase. Further, through crystallite growth, the X-ray crystalline PZT is formed. This method can be extended for the preparation of PLZT powder as well. The resulting powders are sinterable to high density ceramics.

Comparison of active-set method deconvolution and matched-filtering for derivation of an ultrasound transit time spectrum

The quality of ultrasound computed tomography imaging is primarily determined by the accuracy of ultrasound transit time measurement. A major problem in analysis is the overlap of signals making it difficult to detect the correct transit time. The current standard is to apply a matched-filtering approach to the input and output signals. This study compares the matched-filtering technique with active set deconvolution to derive a transit time spectrum from a coded excitation chirp signal and the measured output signal. The ultrasound wave travels in a direct and a reflected path to the receiver, resulting in an overlap in the recorded output signal. The matched-filtering and deconvolution techniques were applied to determine the transit times associated with the two signal paths. Both techniques were able to detect the two different transit times; while matched-filtering has a better accuracy (0.13 μs vs. 0.18 μs standard deviation), deconvolution has a 3.5 times improved side-lobe to main-lobe ratio. A higher side-lobe suppression is important to further improve image fidelity. These results suggest that a future combination of both techniques would provide improved signal detection and hence improved image fidelity.

A Method for the Deposition of Transparent Conducting Thin-Films of Tin Oxide

A modified method has been developed for the deposition of transparent semiconducting thin films of tin oxide, involving the chemical vapour phase oxidation of tin iodide. These films show sheet resistances greater than 100 Ω/□ and an average optical transmission in the visible range exceeding 80%. The method avoids uncontrolled contamination, resulting in better reproducibility of the films. The films showed direct and indirect transitions and the possibility of an indirect forbidden transition. X-ray diffraction studies reveal that the films are polycrystalline. The low mobility values of the films have been attributed to the grain boundary scattering effect.

Efficient parameter estimation in longitudinal data analysis using a hybrid GEE method

The method of generalized estimating equations (GEEs) provides consistent estimates of the regression parameters in a marginal regression model for longitudinal data, even when the working correlation model is misspecified (Liang and Zeger, 1986). However, the efficiency of a GEE estimate can be seriously affected by the choice of the working correlation model. This study addresses this problem by proposing a hybrid method that combines multiple GEEs based on different working correlation models, using the empirical likelihood method (Qin and Lawless, 1994). Analyses show that this hybrid method is more efficient than a GEE using a misspecified working correlation model. Furthermore, if one of the working correlation structures correctly models the within-subject correlations, then this hybrid method provides the most efficient parameter estimates. In simulations, the hybrid method's finite-sample performance is superior to a GEE under any of the commonly used working correlation models and is almost fully efficient in all scenarios studied. The hybrid method is illustrated using data from a longitudinal study of the respiratory infection rates in 275 Indonesian children.

An E-Convergence Method For The Evaluation Of Critical Parameters

A simple cconversence technique is applied to obtain accurate estimates of critical temperatures and critical it\ponmts of a few two- and threpdiniensional king models. When applied to the virial series for hard spheres and hard discs, this method predicts a divergence of the equation-of-state at the density of closest packing.

A simple method for estimating growth parameters from multiple length-frequency data in presence of continuous recruitment

The extended recruitment season for short-lived species such as prawns biases the estimation of growth parameters from length-frequency data when conventional methods are used. We propose a simple method for overcoming this bias given a time series of length-frequency data. The difficulties arising from extended recruitment are eliminated by predicting the growth of the succeeding samples and the length increments of the recruits in previous samples. This method requires that some maximum size at recruitment can be specified. The advantages of this multiple length-frequency method are: it is simple to use; it requires only three parameters; no specific distributions need to be assumed; and the actual seasonal recruitment pattern does not have to be specified. We illustrate the new method with length-frequency data on the tiger prawn Penaeus esculentus from the north-western Gulf of Carpentaria, Australia.

A quasi-likelihood method for fractal-dimension estimation

We propose a simple method of constructing quasi-likelihood functions for dependent data based on conditional-mean-variance relationships, and apply the method to estimating the fractal dimension from box-counting data. Simulation studies were carried out to compare this method with the traditional methods. We also applied this technique to real data from fishing grounds in the Gulf of Carpentaria, Australia

An extension of the continual reassessment method using decision theory

The primary goal of a phase I trial is to find the maximally tolerated dose (MTD) of a treatment. The MTD is usually defined in terms of a tolerable probability, q*, of toxicity. Our objective is to find the highest dose with toxicity risk that does not exceed q*, a criterion that is often desired in designing phase I trials. This criterion differs from that of finding the dose with toxicity risk closest to q*, that is used in methods such as the continual reassessment method. We use the theory of decision processes to find optimal sequential designs that maximize the expected number of patients within the trial allocated to the highest dose with toxicity not exceeding q*, among the doses under consideration. The proposed method is very general in the sense that criteria other than the one considered here can be optimized and that optimal dose assignment can be defined in terms of patients within or outside the trial. It includes as an important special case the continual reassessment method. Numerical study indicates the strategy compares favourably with other phase I designs.

Mesh-free approximations via the error reproducing kernel method and applications to nonlinear systems developing shocks

Error estimates for the error reproducing kernel method (ERKM) are provided. The ERKM is a mesh-free functional approximation scheme [A. Shaw, D. Roy, A NURBS-based error reproducing kernel method with applications in solid mechanics, Computational Mechanics (2006), to appear (available online)], wherein a targeted function and its derivatives are first approximated via non-uniform rational B-splines (NURBS) basis function. Errors in the NURBS approximation are then reproduced via a family of non-NURBS basis functions, constructed using a polynomial reproduction condition, and added to the NURBS approximation of the function obtained in the first step. In addition to the derivation of error estimates, convergence studies are undertaken for a couple of test boundary value problems with known exact solutions. The ERKM is next applied to a one-dimensional Burgers equation where, time evolution leads to a breakdown of the continuous solution and the appearance of a shock. Many available mesh-free schemes appear to be unable to capture this shock without numerical instability. However, given that any desired order of continuity is achievable through NURBS approximations, the ERKM can even accurately approximate functions with discontinuous derivatives. Moreover, due to the variation diminishing property of NURBS, it has advantages in representing sharp changes in gradients. This paper is focused on demonstrating this ability of ERKM via some numerical examples. Comparisons of some of the results with those via the standard form of the reproducing kernel particle method (RKPM) demonstrate the relative numerical advantages and accuracy of the ERKM.