Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Numerical modeling of the non-isothermal liquid droplet impact on a hot solid substrate

The heat transfer from a solid phase to an impinging non-isothermal liquid droplet is studied numerically. A new approach based on an arbitrary Lagrangian-Eulerian (ALE) finite element method for solving the incompressible Navier Stokes equations in the liquid and the energy equation within the solid and the liquid is presented. The novelty of the method consists in using the ALE-formulation also in the solid phase to guarantee matching grids along the liquid solid interface. Moreover, a new technique is developed to compute the heat flux without differentiating the numerical solution. The free surface and the liquid solid interface of the droplet are represented by a moving mesh which can handle jumps in the material parameter and a temperature dependent surface tension. Further, the application of the Laplace-Beltrami operator technique for the curvature approximation allows a natural inclusion of the contact angle. Numerical simulation for varying Reynold, Weber, Peclet and Biot numbers are performed to demonstrate the capabilities of the new approach. (C) 2014 Elsevier Ltd. All rights reserved.

Design and Modeling of an Active Five-Axis Compliant Micromanipulator

This paper presents the design and modeling of an active five-axis compliant micromanipulator whose tip orientation can be independently controlled by large angles about two axes and the tip-position can be controlled in three dimensions. These features enable precise control of the contact point of the tip and the tip-sample interaction forces with three-dimensional nanoscale objects, including those features that are conventionally inaccessible. Control of the tip-motion is realized by means of electromagnetic actuation combined with a novel kinematic and structural design of the micromanipulator, which, in addition, also ensures compatibility with existing high-resolution motion-measurement systems. The design and analysis of the manipulator structure and those of the actuation system are first presented. Quasi-static and dynamic lumped-parameter (LP) models are then derived for the five-axis compliant micromanipulator. Finite element (FE) analysis is employed to validate these models, which are subsequently used to study the effects of tip orientation on the mechanical characteristics of the five-axis micromanipulator. Finally, a prototype of the designed five-axis manipulator is fabricated by means of focused ion-beam milling (FIB).

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.

Girsanov Transformation-Based Reliability Modeling and Testing of Actively Controlled Structures

The problem of estimation of the time-variant reliability of actively controlled structural dynamical systems under stochastic excitations is considered. Monte Carlo simulations, reinforced with Girsanov transformation-based sampling variance reduction, are used to tackle the problem. In this approach, the external excitations are biased by an additional artificial control force. The conflicting objectives of the two control forces-one designed to reduce structural responses and the other to promote limit-state violations (but to reduce sampling variance)-are noted. The control for variance reduction is fashioned after design-point oscillations based on a first-order reliability method. It is shown that for structures that are amenable to laboratory testing, the reliability can be estimated experimentally with reduced testing times by devising a procedure based on the ideas of the Girsanov transformation. Illustrative examples include studies on a building frame with a magnetorheologic damper-based isolation system subject to nonstationary random earthquake excitations. (C) 2014 American Society of Civil Engineers.

Weak conservation of structural features in the interfaces of homologous transient protein-protein complexes

Residue types at the interface of protein-protein complexes (PPCs) are known to be reasonably well conserved. However, we show, using a dataset of known 3-D structures of homologous transient PPCs, that the 3-D location of interfacial residues and their interaction patterns are only moderately and poorly conserved, respectively. Another surprising observation is that a residue at the interface that is conserved is not necessarily in the interface in the homolog. Such differences in homologous complexes are manifested by substitution of the residues that are spatially proximal to the conserved residue and structural differences at the interfaces as well as differences in spatial orientations of the interacting proteins. Conservation of interface location and the interaction pattern at the core of the interfaces is higher than at the periphery of the interface patch. Extents of variability of various structural features reported here for homologous transient PPCs are higher than the variation in homologous permanent homomers. Our findings suggest that straightforward extrapolation of interfacial nature and inter-residue interaction patterns from template to target could lead to serious errors in the modeled complex structure. Understanding the evolution of interfaces provides insights to improve comparative modeling of PPC structures.

Saturated duration for dynamic structural plastic response and fundamental periods of elastic vibration

The relationship is determined between saturated duration of rectangular pressure pulses applied to rigid, perfectly plastic structures and their fundamental periods of elastic vibration. It is shown that the ratio between the saturated duration and the fundamental period of elastic vibration of a structure is dependent upon two factors: the first one is the slenderness or thinness ratio of the structure; and the second one is the square root of ratio between the Young's elastic modulus and the yield stress of the structural material. Dimensional analysis shows that the aforementioned ratio is one of the basic similarity parameters for elastic-plastic modeling under dynamic loading.

A Discussion on Modeling Shape Memory Alloy Embedded in A Composite Laminate as Axial Force and Elastic Foundation

In this paper, the possible error sources of the composite natural frequencies due to modeling the shape memory alloy (SMA) wire as an axial force or an elastic foundation and anisotropy are discussed. The great benefit of modeling the SMA wire as an axial force and an elastic foundation is that the complex constitutive relation of SMA can be avoided. But as the SMA wire and graphite-epoxy are rigidly bonded together, such constraint causes the re-distribution of the stress in the composite. This, together with anisotropy, which also reduces the structural stiffness can cause the relatively large error between the experimental data and theoretical results.

Prediction of structural dynamic plastic shear failure by Johnson's damage number

It is proved that Johnson's damage number is the sole similarity parameter for dynamic plastic shear failure of structures loaded impulsively, therefore, dynamic plastic shear failure can be understood when damage number reaches a critical value. It is suggested that the damage number be generally used to predict the dynamic plastic shear failure of structures under various kinds of dynamic loads (impulsive loading, rectangular pressure pulse, exponential pressure pulse, etc.,). One of the advantages for using the damage number to predict such kind of failure is that it is conveniently used for dissimilar material modeling.

Centrifuge modeling of suction bucket foundations for platforms under ice-sheet-induced cyclic lateral loadings

This paper presents the results of a series of centrifuge model tests performed to study the behavior of suction bucket foundations for a tension leg platform in the Bohai Bay, China. The target lateral loadings were from ice-sheet-induced structural vibrations at a frequency of 0.8-1.0 Hz. The results indicate that excess pore water pressures reach the highest values within a depth of 1.0-1.5 in below the mud line. The pore pressures and the induced settlement and lateral displacement increase with the amplitude of the cyclic loading. Two failure modes were observed: liquefaction in early excitations and settlement-induced problems after long-term excitations. (c) 2006 Elsevier Ltd. All rights reserved.

Modeling Environmental Stress

The word stress when applied to ecosystems is ambiguous. Stress may be low-level, with accompanying near-linear strain, or it may be of finite magnitude, with nonlinear response and possible disintegration of the system. Since there are practically no widely accepted definitions of ecosystem strain, classification of models of stressed systems is tenuous. Despite appearances, most ecosystem models seem to fall into the low-level linear response category. Although they sometimes simulate systems behavior well, they do not provide necessary and sufficient information about sudden structural changes nor structure after transition. Dynamic models of finiteamplitude response to stress are rare because of analytical difficulties. Some idea as to future transition states can be obtained by regarding the behavior of unperturbed functions under limiting strain conditions. Preliminary work shows that, since community variables do respond in a coherent manner to stress, macroscopic analyses of stressed ecosystems offer possible alternatives to compartmental models.

Structural requirements for catalysis: studies on RTEM-1 β-lactamase

β-lactamases are a group of enzymes that confer resistance to penam and cephem antibiotics by hydrolysis of the β-lactam ring, thereby inactivating the antibiotic. Crystallographic and computer modeling studies of RTEM-1 β-lactamase have indicated that Asp 132, a strictly conserved residue among the class A β-lactamases, appears to be involved in substrate binding, catalysis, or both. To study the contribution of residue 132 to β-lactamase function, site saturation mutagenesis was used to generate mutants coding for all 20 amino acids at position 132. Phenotypic screening of all mutants indicated that position 132 is very sensitive to amino acid changes, with only N132C, N132D, N132E, and N132Q showing any appreciable activity. Kinetic analysis of three of these mutants showed increases in K_M, along with substantial decreases in k_(cat). Efforts to trap a stable acyl-enzyme intermediate were unsuccessfuL These results indicate that residue 132 is involved in substrate binding, as well as catalysis, and supports the involvement of this residue in acylation as suggested by Strynadka et al.

Crystallographic and computer modeling studies of RTEM-1 β-lactamase have indicated that Lys 73 and Glu 166, two strictly conserved residues among the class A β-lactamases, appear to be involved in substrate binding, catalysis, or both. To study the contribution of these residues to β-lactamase function, site saturation mutagenesis was used to generate mutants coding for all 20 amino acids at positions 73 and 166. Then all 400 possible combinations of mutants were created by combinatorial mutagenesis. The colonies harboring the mutants were screened for growth in the presence of ampicillin. The competent colonys' DNA were sequenced, and kinetic parameters investigated. It was found that lysine is essential at position 73, and that position 166 only tolerated fairly conservative changes (Aspartic acid, Histidine, and Tyrosine). These functional mutants exhibited decreased kcat's, but K_M was close to wild-type levels. The results of the combinatorial mutagenesis experiments indicate that Lysis absolutely required for activity at position 73; no mutation at residue 166 can compensate for loss of the long side chain amine. The active mutants found--K73K/E166D, K73KIE166H, and K73KIE166Y were studied by kinetic analysis. These results reaffirmed the function of residue 166 as important in catalysis, specifically deacylation.

The identity of the residue responsible for enhancing the active site serine (Ser 70) in RTEM-1 β-lactamase has been disputed for some time. Recently, analysis of a crystal structure of RTEM-1 β-lactamase with covalently bound intermediate was published, and it was suggested that Lys 73, a strictly conserved residue among the class A β-lactamases, was acting as a general base, activating Ser 70. For this to be possible, the pK_a of Lys 73 would have to be depressed significantly. In an attempt to assay the pK_a of Lys 73, the mutation K73C was made. This mutant protein can be reacted with 2-bromoethylamine, and activity is restored to near wild type levels. ^(15)N-2-bromoethylamine hydrobromide and ^(13)C-2-bromoethylamine hydrobromide were synthesized. Reacting these compounds with the K73C mutant gives stable isotopic enrichment at residue 73 in the form of aminoethylcysteine, a lysine homologue. The pK_a of an amine can be determined by NMR titration, following the change in chemical shift of either the ^(15)N-amine nuclei or adjacent Be nuclei as pH is changed. Unfortunately, low protein solubility, along with probable label scrambling in the Be experiment, did not permit direct observation of either the ^(15)N or ^(13)C signals. Indirect detection experiments were used to observe the protons bonded directly to the ^(13)C atoms. Two NMR signals were seen, and their chemical shift change with pH variation was noted. The peak which was determined to correspond to the aminoethylcysteine residue shifted from 3.2 ppm down to 2.8 ppm over a pH range of 6.6 to 12.5. The pK_a of the amine at position 73 was determined to be ~10. This indicates that residue 73 does not function as a general base in the acylation step of the reaction. However the experimental measurement takes place in the absence of substrate. Since the enzyme undergoes conformational changes upon substrate binding, the measured pK_a of the free enzyme may not correspond to the pK_a of the enzyme substrate complex.