Decoupled three-dimensional finite element computation of thermoelastic damping using Zener's approximation

We consider three dimensional finite element computations of thermoelastic damping ratios of arbitrary bodies using Zener's approach. In our small-damping formulation, unlike existing fully coupled formulations, the calculation is split into three smaller parts. Of these, the first sub-calculation involves routine undamped modal analysis using ANSYS. The second sub-calculation takes the mode shape, and solves on the same mesh a periodic heat conduction problem. Finally, the damping coefficient is a volume integral, evaluated elementwise. In the only other decoupled three dimensional computation of thermoelastic damping reported in the literature, the heat conduction problem is solved much less efficiently, using a modal expansion. We provide numerical examples using some beam-like geometries, for which Zener's and similar formulas are valid. Among these we examine tapered beams, including the limiting case of a sharp tip. The latter's higher-mode damping ratios dramatically exceed those of a comparable uniform beam.

Holographic optical element based single-mode hybrid fiber optic interferometer for realizing zero-order fringe

A holographic optical element (HOE) based single-mode hybrid fiber optic interferometer for realizing the zero-order fringe is described. The HOE proposed and used integrates the actions of a beam combiner and a lens, and endows the interferometer with high tolerance for repositioning errors. The proposed method is simple and offers advantages such as the elimination of in situ processing for the hologram.

Optimal pulse spacing for dynamical decoupling in the presence of a purely dephasing spin bath

Maintaining quantum coherence is a crucial requirement for quantum computation; hence protecting quantum systems against their irreversible corruption due to environmental noise is an important open problem. Dynamical decoupling (DD) is an effective method for reducing decoherence with a low control overhead. It also plays an important role in quantum metrology, where, for instance, it is employed in multiparameter estimation. While a sequence of equidistant control pulses the Carr-Purcell-Meiboom-Gill (CPMG) sequence] has been ubiquitously used for decoupling, Uhrig recently proposed that a nonequidistant pulse sequence the Uhrig dynamic decoupling (UDD) sequence] may enhance DD performance, especially for systems where the spectral density of the environment has a sharp frequency cutoff. On the other hand, equidistant sequences outperform UDD for soft cutoffs. The relative advantage provided by UDD for intermediate regimes is not clear. In this paper, we analyze the relative DD performance in this regime experimentally, using solid-state nuclear magnetic resonance. Our system qubits are C-13 nuclear spins and the environment consists of a H-1 nuclear spin bath whose spectral density is close to a normal (Gaussian) distribution. We find that in the presence of such a bath, the CPMG sequence outperforms the UDD sequence. An analogy between dynamical decoupling and interference effects in optics provides an intuitive explanation as to why the CPMG sequence performs better than any nonequidistant DD sequence in the presence of this kind of environmental noise.

Probing the inter-tube transport in aligned and random multiwall carbon nanotubes

The temperature and magnetic field dependence of conductivity has been used to probe the inter-tube transport in multiwall carbon nanotubes (MWNTs). The scanning electron microscopy images show highly aligned and random distribution of MWNTs. The conductivity in aligned carbon nanotube (ACNT) and random carbon nanotube (RCNT) samples at low temperature follows T-1/2 (at T < 8 K) and T-3/4 (at T > 8 K) dependence in accordance with the weak localization and electron-electron (e-e) interaction model. The values of diffusion coefficient in ACNT and RCNT are 0.25 x 10(-2) and 0.71 x 10(-2) cm(2) s(-1), respectively, indicating that larger number of inter-tube junctions in later enhances the bulk transport. The positive magnetoconductance (MC) data in both samples show that the weak localization contribution is dominant. However, the saturation of MC at higher fields and lower temperatures indicate that e-e interaction is quite significant in RCNT. The T-3/4 and T-1/2 dependence of inelastic scattering length (l(in)) in ACNT and RCNT samples show that the inelastic e-e scattering is more important in aligned tubes. (C) 2011 American Institute of Physics. doi:10.1063/1.3552911]

On the performance of strain smoothing for quadratic and enriched finite element approximations (XFEM/GFEM/PUFEM)

By using the strain smoothing technique proposed by Chen et al. (Comput. Mech. 2000; 25: 137-156) for meshless methods in the context of the finite element method (FEM), Liu et al. (Comput. Mech. 2007; 39(6): 859-877) developed the Smoothed FEM (SFEM). Although the SFEM is not yet well understood mathematically, numerical experiments point to potentially useful features of this particularly simple modification of the FEM. To date, the SFEM has only been investigated for bilinear and Wachspress approximations and is limited to linear reproducing conditions. The goal of this paper is to extend the strain smoothing to higher order elements and to investigate numerically in which condition strain smoothing is beneficial to accuracy and convergence of enriched finite element approximations. We focus on three widely used enrichment schemes, namely: (a) weak discontinuities; (b) strong discontinuities; (c) near-tip linear elastic fracture mechanics functions. The main conclusion is that strain smoothing in enriched approximation is only beneficial when the enrichment functions are polynomial (cases (a) and (b)), but that non-polynomial enrichment of type (c) lead to inferior methods compared to the standard enriched FEM (e.g. XFEM). Copyright (C) 2011 John Wiley & Sons, Ltd.

Asymptotic and finite element analyses of mode III dynamic crack growth at a ductile-brittle interface

In this work, dynamic crack growth along a ductile-brittle interface under anti-plane strain conditions is studied. The ductile solid is taken to obey the J(2) flow theory of plasticity with linear isotropic strain hardening, while the substrate is assumed to exhibit linear elastic behavior. Firstly, the asymptotic near-tip stress and velocity fields are derived. These fields are assumed to be variable-separable with a power singularity in the radial coordinate centered at the crack tip. The effects of crack speed, strain hardening of the ductile phase and mismatch in elastic moduli of the two phases on the singularity exponent and the angular functions are studied. Secondly, full-field finite element analyses of the problem under small-scale yielding conditions are performed. The validity of the asymptotic fields and their range of dominance are determined by comparing them with the results of the full-field finite element analyses. Finally, theoretical predictions are made of the variations of the dynamic fracture toughness with crack velocity. The influence of the bi-material parameters on the above variation is investigated.

Inter And Intramolecular Interactions Of Porphyrins With 5,5′-Dithiobis(2-Nitrobenzoic Acid)

meso-Tetraphenylporphyrin and its metal [zinc(II) and copper(II)] derivatives form both inter and intramolecular complexes with 5,5′-dithiobis(2-nitrobenzoic acid) (DTNB). The nature of interaction is predominantly charge transfer (CT) in origin, with the porphyrin functioning as a II-donor and DTNB as an acceptor. Among the covalently linked intramolecular systems, the magnitude of CT interaction varies with the position (of one of the aryl groups of the porphyrin) to which DTNB is attached as ortho meta > para. Steady-state and time-resolved fluorescence studies revealed electron transfer to be the dominant pathway for the fluorescence quenching in these systems. Steady-state photolysis experiments probed using EPR and optical absorption studies have shown that electron transfer (from the excited singlet state of the porphyrin) to DTNB results in the formation of thiyl radical and production of free thiolate anion. It is found that the products of electrochemical reduction of covalently linked porphyrin-DTNB systems are different from those observed for the photochemical studies.

Dynamic stiffness matrix of a general cable element

A computational scheme for determining the dynamic stiffness coefficients of a linear, inclined, translating and viscously/hysteretically damped cable element is outlined. Also taken into account is the coupling between inplane transverse and longitudinal forms of cable vibration. The scheme is based on conversion of the governing set of quasistatic boundary value problems into a larger equivalent set of initial value problems, which are subsequently numerically integrated in a spatial domain using marching algorithms. Numerical results which bring out the nature of the dynamic stiffness coefficients are presented. A specific example of random vibration analysis of a long span cable subjected to earthquake support motions modeled as vector gaussian random processes is also discussed. The approach presented is versatile and capable of handling many complicating effects in cable dynamics in a unified manner.

Characterization of a Negative cis-Acting DNA Element Regulating the Transcription of CYP2B1/B2 Gene in Rat Liver

The region -160 to -127 nt of the upstream of CYP-2B1/B2 gene has been found to function as a negative cis-acting element on the basis of DNase-I footprint and gel mobility shift assays as well as cell-free transcriptional assays using Bal-31 mutants. A reciprocal relationship in the interaction of the negative and the recently characterized positive elements with their respective protein factors has been found under repressed and induced conditions of the gene. The negative element also harbors the core glucocorticoid responsive sequence, TGTCCT. It is concluded that the negative element mediates the repressed state of the gene under the uninduced condition and also mediates the repressive effect of dexamethasone, when given along with the inducer phenobarbitone in rats. Dexamethasone is able to antagonize the effects of phenobarbitone at as low a concentration as 100 mu g/kg body wt in these animals. (C) 1995 Academic Press,Inc.

Engineered dimer interface mutants of triosephosphate isomerase: the role of inter-subunit interactions in enzyme function and stability

The role of inter-subunit interactions in maintaining optimal catalytic activity in triosephosphate isomerase (TIM) has been probed, using the Plasmodium falciparum enzyme as a model. Examination of subunit interface contacts in the crystal structures suggests that residue 75 (Thr, conserved) and residue 13 (Cys, variable) make the largest number of inter-subunit contacts. The mutants Cys13Asp (C13D) and Cys13Glu (C13E) have been constructed and display significant reduction in catalytic activity when compared with wild-type (WT) enzyme (similar to 7.4-fold decrease in k(cat) for the C13D and similar to 3.3-fold for the C13E mutants). Analytical gel filtration demonstrates that the C13D mutant dissociates at concentrations < 1.25 mu M, whereas the WT and the C13E enzymes retain the dimeric structure. The order of stability of the mutants in the presence of chemical denaturants, like urea and guanidium chloride, is WT > Cys13Glu > Cys13Asp. Irreversible thermal precipitation temperatures follow the same order as well. Modeling studies establish that the Cys13Asp mutation is likely to cause a significantly greater structural perturbation than Cys13Glu. Analysis of sequence and structural data for TIMs from diverse sources suggests that residues 13 and 82 form a pair of proximal sites, in which a limited number of residue pairs may be accommodated.

Estimation of low-velocity impact damage in laminated composite circular plates using nonlinear finite element analysis

Nonlinear finite element analysis is used for the estimation of damage due to low-velocity impact loading of laminated composite circular plates. The impact loading is treated as an equivalent static loading by assuming the impactor to be spherical and the contact to obey Hertzian law. The stresses in the laminate are calculated using a 48 d.o.f. laminated composite sector element. Subsequently, the Tsai-Wu criterion is used to detect the zones of failure and the maximum stress criterion is used to identify the mode of failure. Then the material properties of the laminate are degraded in the failed regions. The stress analysis is performed again using the degraded properties of the plies. The iterative process is repeated until no more failure is detected in the laminate. The problem of a typical T300/N5208 composite [45 degrees/0 degrees/-45 degrees/90 degrees](s) circular plate being impacted by a spherical impactor is solved and the results are compared with experimental and analytical results available in the literature. The method proposed and the computer code developed can handle symmetric, as well as unsymmetric, laminates. It can be easily extended to cover the impact of composite rectangular plates, shell panels and shells.

Finite element simulation of wedge indentation

A finite element simulation of frictionless wedge indentation of a copper strip has been carried out under plane strain conditions. The problem was first modelled using an one-pass contact algorithm. The difficulties associated with using this method to model wedge indentation problems are explained. An alternative procedure which alleviates some of the problems associated with the one-pass contact algorithm is proposed for modelling frictionless wedge indentation. Also, a re-meshing procedure which has to be carried out when the distortion of the elements around the indenter becomes significant, is discussed. A sample problem involving indentation of a 4 mm copper strip by a rigid wedge indenter has been modelled and the results are compared with experimental and theoretical results.

Study of Layering Procedures in Finite-Element Analysis of RC Flexural and Torsional Elements

A new postcracking formulation for concrete, along with both implicit and explicit layering procedures, is used in the analysis of reinforced-concrete (RC) flexural and torsional elements. The postcracking formulation accounts for tension stiffening in concrete along the rebar directions, compression softening in cracked concrete based on either stresses or strains, and aggregate interlock based on crack-confining normal stresses. Transverse shear stresses computed using the layering procedures are included in material model considerations that permit the development of inclined cracks through the RC cross section. Examples of a beam analyzed by both the layering techniques, a torsional element, and a column-slab connection region analyzed by the implicit layering procedure are presented here. The study highlights the primary advantages and disadvantages of each layering approach, identifying the class of problems where the application of either procedure is more suitable.

Generation of an eight node brick element using Papcovitch-Neuber functions

Some conventional finite elements suffer from drawbacks, such as shear locking, membrane locking, etc. To overcome them researchers have developed various techniques, termed as tricks by some and variational crimes by others. Many attempts have been made, but satisfactory explanations for why some of these techniques work have not been obtained, especially in the case of solid elements. This paper attempts a simple non-conforming solid element using assumed displacement fields which satisfy the Navier equation exactly. Its behaviour under simple loadings like bending, torsion and tension is examined and comparisons are made with existing elements.

Large deformation finite element analysis of laminated circular composite plates

A 48 d.o.f., four-noded quadrilateral laminated composite shell finite element is particularised to a sector finite element and is used for the large deformation analysis of circular composite laminated plates. The strain-displacement relationships for the sector element are obtained by reducing those of the quadrilateral shell finite element by substituting proper values for the geometric parameters. Subsequently, the linear and tangent stiffness matrices are formulated using conventional methods. The Newton-Raphson method is employed as the nonlinear solution technique. The computer code developed is validated by solving an isotropic case for which results are available in the literature. The method is then applied to solve problems of cylindrically orthotropic circular plates. Some of the results of cylindrically orthotropic case are compared with those available in the literature. Subsequently, application is made to the case of laminated composite circular plates having different lay-up schemes. The computer code can handle symmetric/unsymmetric lay-up schemes. The large displacement analysis is useful in estimating the damage in composite plates caused by low-velocity impact.