The partition of unity finite element method for elastic wave propagation in Reissner-Mindlin plates

This paper reports a numerical method for modelling the elastic wave propagation in plates. The method is based on the partition of unity approach, in which the approximate spectral properties of the infinite dimensional system are embedded within the space of a conventional finite element method through a consistent technique of waveform enrichment. The technique is general, such that it can be applied to the Lagrangian family of finite elements with specific waveform enrichment schemes, depending on the dominant modes of wave propagation in the physical system. A four-noded element for the Reissner-indlin plate is derived in this paper, which is free of shear locking. Such a locking-free property is achieved by removing the transverse displacement degrees of freedom from the element nodal variables and by recovering the same through a line integral and a weak constraint in the frequency domain. As a result, the frequency-dependent stiffness matrix and the mass matrix are obtained, which capture the higher frequency response with even coarse meshes, accurately. The steps involved in the numerical implementation of such element are discussed in details. Numerical studies on the performance of the proposed element are reported by considering a number of cases, which show very good accuracy and low computational cost. Copyright (C)006 John Wiley & Sons, Ltd.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

A parametric study on the buckling of functionally graded material plates with internal discontinuities using the partition of unity method

In this paper, the effect of local defects, viz., cracks and cutouts on the buckling behaviour of functionally graded material plates subjected to mechanical and thermal load is numerically studied. The internal discontinuities, viz., cracks and cutouts are represented independent of the mesh within the framework of the extended finite element method and an enriched shear flexible 4-noded quadrilateral element is used for the spatial discretization. The properties are assumed to vary only in the thickness direction and the effective properties are estimated using the Mori-Tanaka homogenization scheme. The plate kinematics is based on the first order shear deformation theory. The influence of various parameters, viz., the crack length and its location, the cutout radius and its position, the plate aspect ratio and the plate thickness on the critical buckling load is studied. The effect of various boundary conditions is also studied. The numerical results obtained reveal that the critical buckling load decreases with increase in the crack length, the cutout radius and the material gradient index. This is attributed to the degradation in the stiffness either due to the presence of local defects or due to the change in the material composition. (C) 2013 Elsevier Masson SAS. All rights reserved.

Static analysis of infinite laminated cylindrical shell panel

This paper presents an approximate three-dimensional elasticity solution for an infinitely long, cross-ply laminated circular cylindrical shell panel with simply supported boundary conditions, subjected to an arbitrary discontinuous transverse loading. The solution is based on the principal assumption that the ratio of the thickness of the lamina to its middle surface radius is negligible compared to unity. The validity of this assumption and the range of application of this approximate solution have been established through a comparison with an exact solution. Results of classical and first-order shear deformation shell theories have been compared with the results of the present solution to bring out the accuracy of these theories. It is also shown that for very shallow shell panels the definition of a thin shell should be based on the ratio of thickness to chord width rather than the ratio of thickness to mean radius.

Thermodynamics of ligand (substrate end product) binding to endoxylanase from Chainia sp. (NCL-82-5-1): isothermal calorimetry and fluorescence titration studies

The binding of xylo-oligosaccharides to Chainia endoxylanase resulted in a decrease in fluorescence intensity of the enzyme with the formation of 1:1 complex. Equilibrium and thermodynamic parameters of ligand binding were determined by fluorescence titrations and titration calorimetry. The affinity of xylanase for the oligosaccharides increases in the order X-2 < X-3 < X-4 less than or equal to X-5. Contributions from the enthalpy towards the free energy change decreased with increasing chain length from X-2 to X-4, whereas an increase in entropy was observed, the change in enthalpy and entropy of binding being compensatory. The entropically driven binding process suggested that hydrophobic interactions as well as hydrogen bonds play a predominant role in ligand binding.

A Novel Operator-Splitting Technique for One-Dimensional Laminar Premixed Flames

This paper is concerned the calculation of flame structure of one-dimensional laminar premixed flames using the technique of operator-splitting. The technique utilizes an explicit method of solution with one step Euler for chemistry and a novel probabilistic scheme for diffusion. The relationship between diffusion phenomenon and Gauss-Markoff process is exploited to obtain an unconditionally stable explicit difference scheme for diffusion. The method has been applied to (a) a model problem, (b) hydrazine decomposition, (c) a hydrogen-oxygen system with 28 reactions with constant Dρ 2 approximation, and (d) a hydrogen-oxygen system (28 reactions) with trace diffusion approximation. Certain interesting aspects of behaviour of the solution with non-unity Lewis number are brought out in the case of hydrazine flame. The results of computation in the most complex case are shown to compare very favourably with those of Warnatz, both in terms of accuracy of results as well as computational time, thus showing that explicit methods can be effective in flame computations. Also computations using the Gear-Hindmarsh for chemistry and the present approach for diffusion have been carried out and comparison of the two methods is presented.

Bearing capacity factor N-c under phi=0 condition for piles in clays

Bearing capacity factor N-c for axially loaded piles in clays whose cohesion increases linearly with depth has been estimated numerically under undrained (phi=0) condition. The Study follows the lower bound limit analysis in conjunction With finite elements and linear programming. A new formulation is proposed for solving an axisymmetric geotechnical stability problem. The variation of N-c with embedment ratio is obtained for several rates of the increase of soil cohesion with depth; a special case is also examined when the pile base was placed on the stiff clay stratum overlaid by a soft clay layer. It was noticed that the magnitude of N-c reaches almost a constant value for embedment ratio greater than unity. The roughness of the pile base and shaft affects marginally the magnitudes of N-c. The results obtained from the present study are found to compare quite well with the different numerical solutions reported in the literature.

Three-dimensional elastic analysis of cracked thick plates under bending fields

A three-dimensional analysis is presented for the bending problem of finite thick plates with through-the-thickness cracks. A general solution is obtained for Navier's equations of the theory of elasticity. It is found that the in-plane stresses and the transverse normal stress at the crack front are singular with an inverse square root singularity, while the transverse shear stresses are of the order of unity. Results from a numerical study indicate that the stress intensity factor, which varies across the thickness, is influenced by the thickness ratio in a significant manner. Results from a parametric study and those from a comparative study with existing finite element values are presented.

Laser flash photolysis study of triplets of cyclobutanethiones

Five cyclobutanethiones with different chromophores at the 3-position were examined for triplet state behaviour in benzene using laser excitation into their low lying nπ*1 band systems. A weak transient absorption attributable to the triplet state is observed in all these cases. Results concerning triplet lifetimes, intersystem crossing yields (S1 → T1), self-quenching kinetics and kinetics of energy transfer to all-trans-1,6-diphenyl-1,3,5-hexatriene and oxygen and quenching by di-t-butyl nitroxide (DTBN) are presented. Intersystem crossing yields estimated with reference to p,p′-dimethoxythiobenzophenone are roughly unity in all five cases. Self-quenching rates are found to be less than diffusion limited and this is attributed to steric crowding at the α positions (dimethyl group). The rates of oxygen and DTBN quenching compare well with those reported for several other thiones in the literature. No transients other than the triplet were detected in the above cyclobutane-thiones.

Protein A: nature's universal anti-antibody

Staphylococcal protein A specifically interacts with immunogobulins. This fact is being used in various disciplines of biology and some of the unique properties of protein A and their applications are summarized in this review.

Crystallization and preliminary X-ray studies of the anti-T lectin from peanut (Arachis hypogaea)

The anti-T lectin from peanut (Arachis hypogaea) crystallizes in the orthorhombic space group P21212 with one tetrameric molecule (Mr 110,000) in the asymmetric unit in a cell of dimensions a = 129.3 Å, B = 126.9 Å and C = 76.9 Å. The crystals are suitable for high resolution work.

Natural frequencies of rectangular orthotropic plates with a pair of parallel edges simply supported

Solutions of the exact characteristic equations for the title problem derived earlier by an extension of Bolotin's asymptotic method are considered. These solutions, which correspond to flexural modes with frequency factor, R, greater than unity, are expressed in convenient forms for all combinations of clamped, simply supported and free conditions at the remaining pair of parallel edges. As in the case of uniform beams, the eigenvalues in the CC case are found to be equal to those of elastic modes in the FF case provided that the Kirchoff's shear condition at a free edge is replaced by the condition. The flexural modes with frequency factor less than unity are also investigated in detail by introducing a suitable modification in the procedure. When Poisson's ratios are not zero, it is shown that the frequency factor corresponding to the first symmetric mode in the free-free case is less than unity for all values of side ratio and rigidity ratios. In the case of one edge clamped and the other free it is found that modes with frequency factor less than unity exist for certain dimensions of the plate—a fact hitherto unrecognized in the literature.

Dilute Solution Properties of Methyl Methacrylate/Acrylonitrile Random Copolymers, 2a)

The Stockmayer-Fixman relation was used to evaluate the short range and long range interaction parameters for methyl methacrylate/acrylonitrile copolymers of 0,566 and 0,657 mole fraction of monomeric units of acrylonitrile in the solvents acetonitrile, 2-butanone, dimethyl formamide, and y-butyrolactone, at different temperatures (30, 45, and 60 “C). The values of KO were found to be lower than those of the parent homopolymers, and their values depend on both solvent and temperature. Even negative Ko-values were obtained, in cases in which the Mark Houwink exponent a is nearly unity. The values of the polymer-solvent interaction parameter, x, , are high and close to 0,5, indicating that these solvents are not good. The values of the excess interaction parameter, xAB, are negative and are not affected by temperature. The large extension of these copolymer chains, as exhibited by a and a;-values, can be understood in terms of unusual short range interactions only. Similar results were obtained for some cellulose derivatives.

Protein A: nature's universal anti-antibody

Staphylococcal protein A specifically interacts with immunogobulins. This fact is being used in various disciplines of biology and some of the unique properties of protein A and their applications are summarized in this review.

Zero-crossings in turbulent signals

A primary motivation for this work arises from the contradictory results obtained in some recent measurements of the zero-crossing frequency of turbulent fluctuations in shear flows. A systematic study of the various factors involved in zero-crossing measurements shows that the dynamic range of the signal, the discriminator characteristics, filter frequency and noise contamination have a strong bearing on the results obtained. These effects are analysed, and explicit corrections for noise contamination have been worked out. New measurements of the zero-crossing frequency N0 have been made for the longitudinal velocity fluctuation in boundary layers and a wake, for wall shear stress in a channel, and for temperature derivatives in a heated boundary layer. All these measurements show that a zero-crossing microscale, defined as Λ = (2πN0)−1, is always nearly equal to the well-known Taylor microscale λ (in time). These measurements, as well as a brief analysis, show that even strong departures from Gaussianity do not necessarily yield values appreciably different from unity for the ratio Λ/λ. Further, the variation of N0/N0 max across the boundary layer is found to correlate with the familiar wall and outer coordinates; the outer scaling for N0 max is totally inappropriate, and the inner scaling shows only a weak Reynolds-number dependence. It is also found that the distribution of the interval between successive zero-crossings can be approximated by a combination of a lognormal and an exponential, or (if the shortest intervals are ignored) even of two exponentials, one of which characterizes crossings whose duration is of the order of the wall-variable timescale ν/U2*, while the other characterizes crossings whose duration is of the order of the large-eddy timescale δ/U[infty infinity]. The significance of these results is discussed, and it is particularly argued that the pulse frequency of Rao, Narasimha & Badri Narayanan (1971) is appreciably less than the zero-crossing rate.