Improved hybrid elements for structural analysis

Hybrid elements, which are based on a two-field variational formulation with the displacements and stresses interpolated separately, are known to deliver very high accuracy, and to alleviate to a large extent problems of locking that plague standard displacement-based formulations. The choice of the stress interpolation functions is of course critical in ensuring the high accuracy and robustness of the method. Generally, an attempt is made to keep the stress interpolation to the minimum number of terms that will ensure that the stiffness matrix has no spurious zero-energy modes, since it is known that the stiffness increases with the increase in the number of terms. Although using such a strategy of keeping the number of interpolation terms to a minimum works very well in static problems, it results either in instabilities or fails to converge in transient problems. This is because choosing the stress interpolation functions merely on the basis of removing spurious energy modes can violate some basic principles that interpolation functions should obey. In this work, we address the issue of choosing the interpolation functions based on such basic principles of interpolation theory and mechanics. Although this procedure results in the use of more number of terms than the minimum (and hence in slightly increased stiffness) in many elements, we show that the performance continues to be far superior to displacement-based formulations, and, more importantly, that it also results in considerably increased robustness.

Elasto-Electrostatic Analysis of Circular Microplates Used in Capacitive Micromachined Ultrasonic Transducers

The active structural component of a capacitive micromachined ultrasonic transducer (CMUT) is the top plate which vibrates under the influence of a time-varying electrostatic force thereby producing ultrasound waves of the desired frequency in the surrounding medium. Analysis of MEMS devices which rely on electrostatic actuation is complicated due to the fact that the structural deformations alter the electrostatic forces, which redistribute and modify the applied loads. Hence, it becomes imperative to consider the electrostatics-structure coupling aspect in the design of these devices. This paper presents an approximate analytical solution for the static deflection of a thin, clamped circular plate caused by electrostatic forces which are inherently nonlinear. Traditionally, finite element simulations using some commercial software such as ANSYS are employed to determine the structural deflections caused by electrostatic forces. Since the structural deformation alters the electrostatic field, a coupled-field simulation is required wherein the electrostatic mesh is continuously updated to coincide with the deflection of the structure. Such simulations are extremely time consuming, in addition to being nontransparent and somewhat hard to implement. We employ the classical thin-plate theory which is adequate when the ratio of the diameter to thickness of the plate is very large, a situation commonly prevalent in many MEMS devices, especially the CMUTs. We solve the thin-plate electrostatic-elastic equation using the Galerkin-weighted residual technique, under the assumption that the deflections are small in comparison to the thickness of the plate. The evaluation of the electrostatic force between the two plates is simplified due to the fact that the electrostatic gap is much smaller than the lateral dimensions of the device. The results obtained are compared to those found from ANSYS simulations and an excellent agreement is observed between the two. The pull-in voltage predicted by our model is close to the value predicted by ANSYS simulations.

A real transformation for modal analysis of flexible damped spacecraft

Modal approach is widely used for the analysis of dynamics of flexible structures. However, space analysts yet lack an intimate modal analysis of current spacecraft which are rich with flexibility and possess both structural and discrete damping. Mathematical modeling of such spacecraft incapacitates the existing real transformation procedure, for it cannot include discrete damping, demands uncomputable inversion of a modal matrix inaccessible due to its overwhelming size and does not permit truncation. On the other hand, complex transformation techniques entail more computational time and cannot handle structural damping. This paper presents a real transformation strategy which averts inversion of the associated real transformation matrix, allows truncation and accommodates both forms of damping simultaneously. This is accomplished by establishing a key relation between the real transformation matrix and its adjoint. The relation permits truncation of the matrices and leads to uncoupled pairs of coupled first order equations which contain a number of adjoint eigenvectors. Finally these pairs are solved to obtain a literal modal response of forced gyroscopic damped flexibile systems at arbitrary initial conditions.

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Chemical toughening of epoxies. I. Structural modification of epoxy resins by reaction with hydroxy-terminated poly (butadiene-co-acrylonitrile)

In an attempt to toughen the epoxy resin matrix for fiber-reinforced composite applications, a chemical modification procedure of a commercially available bisphenol-A-based epoxy resin using reactive liquid rubber HTBN [hydroxy-terminated poly(butadiene-co-acrylonitrile)] and TDI (tolylene diisocyanate) is described. The progress of the reaction and the structural changes during modification process are studied using IR spectroscopy, viscosity data, and chemical analysis (epoxy value determination). The studies support the proposition that TDI acts as a coupling agent between the epoxy and HTBN, forming a urethane linkage with the former and an oxazolidone ring with the latter. The chemical reactions that possibly take place during the modification are discussed.

Chemical toughening of epoxies. I. Structural modification of epoxy resins by reaction with hydroxy-terminated poly(butadiene-co-acrylonitrile)

In an attempt to toughen the epoxy resin matrix for fiber-reinforced composite applications, a chemical modification procedure of a commercially available bisphenol-A-based epoxy resin using reactive liquid rubber HTBN [hydroxy-terminated poly(butadiene-co-acrylonitrile)] and TDI (tolylene diisocyanate) is described. The progress of the reaction and the structural changes during modification process are studied using IR spectroscopy, viscosity data, and chemical analysis (epoxy value determination). The studies support the proposition that TDI acts as a coupling agent between the epoxy and HTBN, forming a urethane linkage with the former and an oxazolidone ring with the latter. The chemical reactions that possibly take place during the modification are discussed.

Analysis of the Ce 3d-4d4d Auger spectrum with the use of synchrotron radiation

We report 3d-4d4d Auger spectra of Ce metal with the use of synchrotron radiation to excite the initial core hole. By sweeping the excitation energy through the 3d-->4f threshold, it has been possible to excite different initial states selectively, enabling us to analyze the complex spectrum in terms of different contributions arising from various deca channels.

Analysis of the Gd 4d-XY Auger spectrum using synchrotron radiation

We report the 4d-XY (X, Y = 5p, 4f, and the conduction band) Auger spectra of clean Gd using a monochromatic photon source with energies above and below the 3d threshold. The spectra with higher hv show the existence of intense spectator-hole Auger transitions. Comparison of these spectra with those obtained with a primary electron source allows detailed interpretation of the various features and explains the unusual spin polarization of the electron-induced spectrum reported earlier.

A macroeconomic analysis of adaptation to climate change impacts on forests in India

We examine the potential for adaptation to climate change in Indian forests, and derive the macroeconomic implications of forest impacts and adaptation in India. The study is conducted by integrating results from the dynamic global vegetation model IBIS and the computable general equilibrium model GRACE-IN, which estimates macroeconomic implications for six zones of India. By comparing a reference scenario without climate change with a climate impact scenario based on the IPCC A2-scenario, we find major variations in the pattern of change across zones. Biomass stock increases in all zones but the Central zone. The increase in biomass growth is smaller, and declines in one more zone, South zone, despite higher stock. In the four zones with increases in biomass growth, harvest increases by only approximately 1/3 of the change in biomass growth. This is due to two market effects of increased biomass growth. One is that an increase in biomass growth encourages more harvest given other things being equal. The other is that more harvest leads to higher supply of timber, which lowers market prices. As a result, also the rent on forested land decreases. The lower prices and rent discourage more harvest even though they may induce higher demand, which increases the pressure on harvest. In a less perfect world than the model describes these two effects may contribute to an increase in the risk of deforestation because of higher biomass growth. Furthermore, higher harvest demands more labor and capital input in the forestry sector. Given total supply of labor and capital, this increases the cost of production in all the other sectors, although very little indeed. Forestry dependent communities with declining biomass growth may, however, experience local unemployment as a result.

Analysis of the effects of amino-acid-sequence on the structure of proteins

The conformation of amino acid side chains as observed in well-determined structures of globular proteins has earlier been extensively investigated. In contrast, the structural features of the polypeptide backbone that result from the occurrence of specific amino acids along the polypeptide have not been analysed. In this article, we present the statistically significant features in the backbone geometry that appear to be a consequence of the occurrence of rotamers of different amino acid side chains by analysing 102 well-refined structures that form a random collection of proteins. It is found that the persistence of helical segments around each residue is influenced by the residue type. Several residues exert asymmetrical influence between the carboxyl and amino terminal polypeptide segments. The degree to which secondary structures depart from an average geometry also appears to depend on residue type. These departures are correlated to the corresponding Chou and Fasman parameters of amino acid residues. The frequency distribution of the side chain rotamers is influenced by polypeptide secondary structure. In turn, the rotamer conformation of side chain affects the extension of the secondary structure of the backbone. The strongest correlation is found between the occurrence of g+ conformation and helix propagation on the carboxyl side of many residues.

Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

Structural characterization of gel-grown neodymium copper oxalate single crystals

Sparingly soluble neodymium copper oxalate (NCO) single crystals were grown by gel method, by the diffusion of a mixture of neodymium nitrate and cupric nitrate into the set gel containing oxalic acid. Tabular crystal, revealing well-defined dissolution figures has been recorded. X-ray diffraction studies of the powdered sample reveal that NCO is crystalline. Infrared absorption spectrum confirmed the formation of oxalato complex with water of crystallization, while energy dispersive X-ray analysis established the presence of neodymium dominant over copper in the sample. X-ray photoelectron spectroscopic studies established the presence of Nd and Cu in oxide states besides (C2O4)(2-) oxalate group. The intensities of Nd (3d(5/2)) and Cu (2p(3/2)) peaks measured in terms of maximum photoelectron count rates also revealed the presence of Nd in predominance. The inductively coupled plasma analysis supports the EDAX and XPS data by the estimation of neodymium percentage by weight to that of copper present in the NCO sample. On the basis of these findings, an empirical structure for NCO has been proposed. The implications are discussed.

Studies on molecular evolution and structural features of double-headed inhibitors of ?-amylase and trypsin in plants

Plant seeds usually have high concentrations of proteinase and amylase inhibitors. These inhibitors exhibit a wide range of specificity, stability and oligomeric structure. In this communication, we report analysis of sequences that show statistically significant similarity to the double-headed alpha-amylase/trypsin inhibitor of ragi (Eleusine coracana). Our aim is to understand their evolutionary and structural features. The 14 sequences of this family that are available in the SWISSPROT database form three evolutionarily distinct branches. The branches relate to enzyme specificities and also probably to the oligomeric state of the proteins and not to the botanical class of the plant from which the enzymes are derived. This suggests that the enzyme specificities of the inhibitors evolved before the divergence of commercially cultivated cereals. The inhibitor sequences have three regions that display periodicity in hydrophobicity. It is likely that this feature reflects extended secondary structure in these segments. One of the most variable regions of the polypeptide corresponds to a loop, which is most probably exposed in the native structure of the inhibitors and is responsible for the inhibitory property.

The primary structure of sheep liver cytosolic serine hydroxymethyltransferase and an analysis of the evolutionary relationships among serine hydroxymethyltransferases

The complete amino-acid sequence of sheep liver cytosolic serine hydroxymethyltransferase was determined from an analysis of tryptic, chymotryptic, CNBr and hydroxylamine peptides. Each subunit of sheep liver serine hydroxymethyltransferase consisted of 483 amino-acid residues. A comparison of this sequence with 8 other serine hydroxymethyltransferases revealed that a possible gene duplication event could have occurred after the divergence of animals and fungi. This analysis also showed independent duplication of SHMT genes in Neurospora crassa. At the secondary structural level, all the serine hydroxymethyltransferases belong to the alpha/beta category of proteins. The predicted secondary structure of sheep liver serine hydroxymethyltransferase was similar to that of the observed structure of tryptophan synthase, another pyridoxal 5'-phosphate containing enzyme, suggesting that sheep liver serine hydroxymethyltransferase might have a similar pyridoxal 5'-phosphate binding domain. In addition, a conserved glycine rich region, G L Q G G P, was identified in all the serine hydroxymethyltransferases and could be important in pyridoxal 5'-phosphate binding. A comparison of the cytosolic serine hydroxymethyltransferases from rabbit and sheep liver with other proteins sequenced from both these sources showed that serine hydroxymethyltransferase was a highly conserved protein. It was slightly less conserved than cytochrome c but better conserved than myoglobin, both of which are well known evolutionary markers. C67 and C203 were specifically protected by pyridoxal 5'-phosphate against modification with [C-14]iodoacetic acid, while C247 and C261 were buried in the native serine hydroxymethyltransferase. However, the cysteines are not conserved among the various serine hydroxymethyltransferases. The exact role of the cysteines in the reaction catalyzed by serine hydroxymethyltransferase remains to be elucidated.