Element-Free Galerkin Modelling of Crack Migration in Composite Laminates

The development of a virtual testing environment, as a cost-effective industrial design tool in the design and analysis of composite structures, requires the need to create models efficiently, as well as accelerate the analysis by reducing the number of degrees of freedom, while still satisfying the need for accurately tracking the evolution of a debond, delamination or crack front. The eventual aim is to simulate both damage initiation and propagation in components with realistic geometrical features, where crack propagation paths are not trivial. Meshless approaches, and the Element-Free Galerkin (EFG) method, are particularly suitable for problems involving changes in topology and have been successfully applied to simulate damage in homogeneous materials and concrete. In this work, the method is utilized to model initiation and mixed-mode propagation of cracks in composite laminates, and to simulate experimentally-observed crack migration which is difficult to model using standard finite element analysis. N

Improving composite damage modelling through automatic placement of cohesive elements

Numerous experimental studies of damage in composite laminates have shown that intralaminar (in-plane) matrix cracks lead to interlaminar delamination (out-of-plane) at ply interfaces. The smearing of in-plane cracks over a volume, as a consequence of the use of continuum damage mechanics, does not always effectively capture the full extent of the interaction between the two failure mechanisms. A more accurate representation is obtained by adopting a discrete crack approach via the use of cohesive elements, for both in-plane and out-of-plane damage. The difficulty with cohesive elements is that their location must be determined a priori in order to generate the model; while ideally the position of the crack migration, and more generally the propagation path, should be obtained as part of the problem’s solution. With the aim of enhancing current modelling capabilities with truly predictive capabilities, a concept of automatic insertion of interface elements is utilized. The consideration of a simple traction criterion in relation to material strength, evaluated at each node of the model (or of the regions of the model where it is estimated cracks might form), allows for the determination of initial crack location and subsequent propagation by the insertion of cohesive elements during the course of the analysis. Several experimental results are modelled using the commercial package ABAQUS/Standard with an automatic insertion subroutine developed in this work, and the results are presented to demonstrate the capabilities of this technique.

The Reliability of the Arc-Length Method in the Analysis of Mode-Jumping Problems

The Arc-Length Method is a solution procedure that enables a generic non-linear problem to pass limit points. Some examples are provided of mode-jumping problems solutions using a commercial nite element package, and other investigations are carried out on a simple structure of which the numerical solution can be compared with an analytical one. It is shown that Arc-Length Method is not reliable when bifurcations are present in the primary equilibrium path; also the presence of very sharp snap-backs or special boundary conditions may cause convergence diÆculty at limit points. An improvement to the predictor used in the incremental procedure is suggested, together with a reliable criteria for selecting either solution of the quadratic arc-length constraint. The gap that is sometimes observed between the experimantal load level of mode-jumping and its arc-length prediction is explained through an example.

Development and Validation of a Procedure for Numerical Vibration Analysis of an Oscillating Wave Surge Converter

During extreme sea states so called impact events can be observed on the wave energy converter Oyster. In small scale experimental tests these impact events cause high frequency signals in the measured load which decrease confidence in the data obtained. These loads depend on the structural dynamics of the model. Amplification of the loads can occur and is transferred through the structure from the point of impact to the load cell located in the foundation. Since the determination of design data and load cases for Wave Energy Converters originate from scale experiments, this lack of confidence has a direct effect on the development.

Numerical vibration analysis is a valuable tool in the research of the structural load response of Oyster to impact events, but must take into account the effect of the surrounding water. This can be done efficiently by adding an added mass distribution, computed with a linearised potential boundary element method. This paper presents the development and validation of a numerical procedure, which couples the OpenSource boundary element code NEMOH with the Finite Element Analysis tool CodeAster. Numerical results of the natural frequencies and mode shapes of the structure under the influence of added mass due to specific structural modes are compared with experimental results.

Molecular Dynamics of HIV1-Integrase in Complex with 93del - A Structural Perspective on the Mechanism of Inhibition

HIV1 integrase is an important target for the antiviral therapy. Guanine-rich quadruplex, such as 93del, have been shown to be potent inhibitors of this enzyme and thus representing a new class of antiviral agents. Although X-ray and NMR structures of HIV1 integrase and 93del have been reported, there is no structural information of the complex and the mechanism of inhibition still remains unexplored. A number of computational methods including automated protein-DNA docking and molecular dynamics simulation in explicit solvent were used to model the binding of 93del to HIV1 integrase. Analysis of the dynamic behaviour of the complex using principal components analysis and elastic network modelling techniques allow us to understand how the binding of 93del aptamer and its interactions with key residues affect the intrinsic motions of the catalytic loops by stabilising them in catalytically inactive conformations. Such insights into the structural mechanism of inhibition can aid in improving the design of anti-HIV aptamers.

Restricted rotation dynamics and far-infrared spectra of liquid water governed by elastic interactions

A semi-phenomenological model describing wideband dielectric and far-infrared spectra of liquid water was proposed recently by the same authors [J. Mol. Struct. 606 (2002) 9], where a small dipole-moment component changing harmonically with time determines a weak absorption band (termed here the R-band) centred at the wavenumber v similar to 200 cm(-1). In the present work, a rough molecular theory of the R-band based on the concept of elastic interactions is given. Stretching and bending of hydrogen bonds cause restricted rotation (RR) of a polar water molecule in terms of a dimer comprising the H- bonded molecules. Analytical expression for the RR frequency nu(str) is derived as a function of the RR amplitude, geometrical parameters and force constants. The density g(nu(str)) of frequency distribution is shown to be centred in the R-band. The spectrum of the dipolar auto-correlation function calculated for this structural-dynamical model is found. A composite model comprising two intermolecular potentials is proposed, which yields for water a good description of the experimental wideband (from 0 to 1000 cm(- 1)) spectra of complex permittivity and of absorption coefficient. The presented interpretation of these spectra is based on a concept that water presents a two-component solution, with components differing by the types of molecular rotation. (C) 2003 Elsevier B.V. All rights reserved.

Molecular dynamics studies of trinucleotide repeat DNA involved in neurodegenerative disorders.

Expansion of trinucleotide repeat DNA of the classes CAG�·CTG, CGG�·CCG and GAA�·TTC are found to be associated with several neurodegenerative disorders. Different mechanisms have been attributed to the expansion of triplets, mainly involving the formation of alternate secondary structures by such repeats. This paper reports the molecular dynamics simulation of triplet repeat DNA sequences to study the basic structural features of DNA that are responsible for the formation of structures such as hairpins and slip-strand DNA leading to expansion. All the triplet repeat sequences studied were found to be more flexible compared to the control sequence unassociated with disease. Moreover, flexibility was found to be in the order CAG�·CTG > CGG�·CCG = GAA�·TTC, the highly flexible CAG�·CTG repeat being the most common cause of neurodegenerative disorders. In another simulation, a single G�·C to T�·A mutation at the 9th position of the CAG�·CTG repeat exhibited a reduction in bending compared to the pure 15-mer CAGâ�¢CTG repeat. EPM1 dodecamer repeat associated with the pathogenesis of progressive myoclonus epilepsy was also simulated and showed flexible nature suggesting a similar expansion mechanism.

Relative energetics and structural properties of zirconia using a self-consistent tight-binding model

We describe an empirical, self-consistent, orthogonal tight-binding model for zirconia, which allows for the polarizability of the anions at dipole and quadrupole levels and for crystal field splitting of the cation d orbitals, This is achieved by mixing the orbitals of different symmetry on a site with coupling coefficients driven by the Coulomb potentials up to octapole level. The additional forces on atoms due to the self-consistency and polarizabilities are exactly obtained by straightforward electrostatics, by analogy with the Hellmann-Feynman theorem as applied in first-principles calculations. The model correctly orders the zero temperature energies of all zirconia polymorphs. The Zr-O matrix elements of the Hamiltonian, which measure covalency, make a greater contribution than the polarizability to the energy differences between phases. Results for elastic constants of the cubic and tetragonal phases and phonon frequencies of the cubic phase are also presented and compared with some experimental data and first-principles calculations. We suggest that the model will be useful for studying finite temperature effects by means of molecular dynamics.

Structure and dynamics of a confined ionic liquid. topics of relevance to dye-sensitized solar cells

The behavior of a model ionic liquid (IL) confined between two flat parallel walls was studied at various interwall distances using computer simulations. The results focus both on structural and dynamical properties. Mass and charge density along the confinement axis reveal a structure of layers parallel to the walls that leads to an oscillatory profile in the electrostatic potential. Orientational correlation functions indicate that cations at the interface orient tilted with respect to the surface and that any other orientational order is lost thereafter. The diffusion coefficients of the ions exhibit a maximum as a function of the confinement distance, a behavior that results from a combination of the structure of the liquid as a whole and a faster molecular motion in the vicinity of the walls. We discuss the relevance of the present results and elaborate on topics that need further attention regarding the effects of ILs in the functioning of IL-based dye-sensitized solar cells.

Thermodynamics, structure, and dynamics in room temperature ionic liquids: The case of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6])

A detailed investigation of the phase diagram of 1-butyl-3-methyl imidazolium hexafluorophosphate ([bmim][PF6]) is presented on the basis of a wide set of experimental data accessing thermodynamic, structural, and dynamical properties of this important room temperature ionic liquid (RTIL). The combination of quasi adiabatic, continuous calorimetry, wide angle neutron and X-ray diffraction, and quasi elastic neutron scattering allows the exploration of many novel features of this material. Thermodynamic and microscopic structural information is derived on both glassy and crystalline states and compared with results that recently appeared in the literature allowing direct information to be obtained on the existence of two crystalline phases that were not previously characterized and confirming the view that RTILs show a substantial degree of order (even in their amorphous states), which resembles the crystalline order. We highlight a strong connection between structure and dynamics, showing the existence of three temperature ranges in the glassy state across which both the spatial correlation and the dynamics change. The complex crystalline polymorphism in [bmim][PF6] also is investigated; we compare our findings with the corresponding findings for similar RTILs. These results provide a strong experimental basis for the exploration of the features of the phase diagram of RTILs and for the further study of longer alkyl chain salts.

Structural basis of substrate selectivity in the glycerol-3-phosphate: phosphate antiporter GlpT.

Major facilitators represent the largest superfamily of secondary active transporter proteins and catalyze the transport of an enormous variety of small solute molecules across biological membranes. However, individual superfamily members, although they may be architecturally similar, exhibit strict specificity toward the substrates they transport. The structural basis of this specificity is poorly understood. A member of the major facilitator superfamily is the glycerol-3-phosphate (G3P) transporter (GlpT) from the Escherichia coli inner membrane. GlpT is an antiporter that transports G3P into the cell in exchange for inorganic phosphate (Pi). By combining large-scale molecular-dynamics simulations, mutagenesis, substrate-binding affinity, and transport activity assays on GlpT, we were able to identify key amino acid residues that confer substrate specificity upon this protein. Our studies suggest that only a few amino acid residues that line the transporter lumen act as specificity determinants. Whereas R45, K80, H165, and, to a lesser extent Y38, Y42, and Y76 contribute to recognition of both free Pi and the phosphate moiety of G3P, the residues N162, Y266, and Y393 function in recognition of only the glycerol moiety of G3P. It is the latter interactions that give the transporter a higher affinity to G3P over Pi.

Salt-bridge dynamics control substrate-induced conformational change in the membrane transporter GlpT.

Active transport of substrates across cytoplasmic membranes is of great physiological, medical and pharmaceutical importance. The glycerol-3-phosphate (G3P) transporter (GlpT) of the E. coli inner membrane is a secondary active antiporter from the ubiquitous major facilitator superfamily that couples the import of G3P to the efflux of inorganic phosphate (Pi) down its concentration gradient. Integrating information from a novel combination of structural, molecular dynamics simulations and biochemical studies, we identify the residues involved directly in binding of substrate to the inward-facing conformation of GlpT, thus defining the structural basis for the substrate-specificity of this transporter. The substrate binding mechanism involves protonation of a histidine residue at the binding site. Furthermore, our data suggest that the formation and breaking of inter- and intradomain salt bridges control the conformational change of the transporter that accompanies substrate translocation across the membrane. The mechanism we propose may be a paradigm for organophosphate:phosphate antiporters.

Structural motifs of cholesterol nanoparticles

The growth sequence of gas-phase cholesterol clusters (Ch(N)) with up to N=36 molecules has been investigated by atomistic simulation based on an empirical force field model. The results of long annealings from high temperature show that the geometric motifs characterizing the structure of pure cholesterol crystals already appear in nanometric aggregates. In all clusters molecules tend to align along a common direction. For cluster sizes above the smallest ones, dispersion interactions among the hydrocarbon body and tails of cholesterol cooperate with hydrogen bonding to give rise to a bilayer structure. Analysis of snapshots from the annealing shows that the condensation of hydrogen bonds into a connected network of rings and chains is an important step in the self-organization of cholesterol clusters. The effect of solvation on the equilibrium properties of medium-size aggregates is investigated by short molecular dynamics simulations for the N=30 and N=40 clusters in water at near ambient conditions and in supercritical carbon dioxide at T=400 K.

Molecular Modeling on Inhibitor Complexes and Active-Site Dynamics of Cytochrome P450 C17, a Target for Prostate Cancer Therapy

A molecular model for the P450 enzyme cytochrome P450 C17 (CYP17) is presented based on sequence alignments of multiple template structures and homology modeling. This enzyme plays a central role in the biosynthesis of testosterone and is emerging as a major target in prostate cancer, with the recently developed inhibitor abiraterone currently in advanced clinical trials. The model is described in detail, together with its validation, by providing structural explanations to available site-directed mutagenesis data. The CYP17 molecule in this model is in the form of a triangular prism, with an edge of similar to 55 angstrom and a thickness of similar to 37 angstrom. It is predominantly helical, comprising 13 alpha helices interspersed by six 3(10) helices and 11 beta-sheets. Multinanosecond molecular dynamics simulations in explicit solvent have been carried out, and principal components analysis has been used to reveal the details of dynamics around the active site. Coarse-grained methods have also been used to verify low-frequency motions, which have been correlated with active-site gating. The work also describes the results of docking synthetic inhibitors, including the drug abiraterone and the natural substrate pregnenolone, in the CYP17 active site together with molecular dynamics simulations on the complexes. (C) 2010 Elsevier Ltd. All rights reserved.

Transonic Aeroelastic Stability Predictions Under the Influence of Structural Variability

Flutter prediction as currently practiced is almost always deterministic in nature, based on a single structural model that is assumed to represent a fleet of aircraft. However, it is also recognized that there can be significant structural variability, even for different flights of the same aircraft. The safety factor used for flutter clearance is in part meant to account for this variability. Simulation tools can, however, represent the consequences of structural variability in the flutter predictions, providing extra information that could be useful in planning physical tests and assessing risk. The main problem arising for this type of calculation when using high-fidelity tools based on computational fluid dynamics is the computational cost. The current paper uses an eigenvalue-based stability method together with Euler-level aerodynamics and different methods for propagating structural variability to stability predictions. The propagation methods are Monte Carlo, perturbation, and interval analysis. The feasibility of this type of analysis is demonstrated. Results are presented for the Goland wing and a generic fighter configuration.