Skip to main content

research

Micro cantilever pre-stress

Submitted by Iskandar Samad on
Choose a channel featured in the header of iMechanica

Dear all,

Im a PhD student in Cambridge Uni, UK working in the field of MEMS, and as part of my work, of late Ive been looking at deriving materials properties of MEMS thin film materials by means of resonant testing. The basic outline of the experiment is first creating free standing rectangular cantilevers of the material under test, evaporate with gold to increase reflectivity (when needed), then (under reasonable vacuum) applying a base excitation using a chirp signal into a piezo actuator and logging the cantilever tip response using a laser/photodetector setup. The frequency response is then calculated and the modal frequencies noted.

To determine the materials' properties, both an analytical model (with bending/torsion modes) and finite element model (using 2D mindlin) are created with similar geometry as the sample, and by minimising the squared relative error between the measured modes and those from the models, the value of Young's modulus (known density) and poissons ratio may be determined iteratively. These yield fairly consistent results.

To take the work further I now feel I should also include the effects of residual stress in the cantilevers. The method Ive been looking at is by using finite element (via COMSOL) - the beam geometry is created and loaded with the stress model ('surface stress' as a force tangential to the top boundary, and gradient stress as a 'tangential' force that varies from +F to -F from top to bottom boundary of the cantilever). The model is solved statically, and the deformed shape is then saved as the linearisation point for the next model, which then computes the eigenfrequencies. Btw I can only do this in 3D FE, which makes computation times quite long hence using iteration to quantify this stress highly unlikely.

In any case, is there an analytical model I can use to model the effect of this stress on multiple modes of a cantilever. Id like to verify whether the FE is giving me anything close to ball park numbers before I work out a means to compare them with experimental results. I was thinking of using the Rayleigh method by representing the effect of prestress as an additional term in the potential energy. The original mode shapes, with some modification will be used to evaluate the two energy integrals. The potential energy due to stress is worked out by measuring the static deflected shape using a zygo inteferometer - some rough model is used, with the beam curvature and peak deflection as input to work out the amount of this energy. Not having much experience in mechanics (i was an electronics undergrad!), Im not sure how good an estimate this would be, if at all its a useable or even possible one. Will the extra energy factor in to the torsion modes just the same?

Model Reduction of Large Proteins for Normal Mode Studies

Submitted by Kilho Eom on

Recently, I reported the model reduction method for large proteins for understanding large protein dynamics based on low-frequency normal modes. This work was pubslihed at Journal of Computational Chemistry (click here).

Coarse-Graining of protein structures for the normal mode studies

Abstracts 

Quantum Stability of Metallic Thin Films and Nanostructures

Submitted by Zhenyu Zhang on

When a metal system shrinks its dimension(s), the conduction electrons inside the metal feel the squeezing, and are forced into (discrete) quantum states. Such confined motion of the conduction electrons may influence the global or local stability of the low dimensional systems, and in the case of a thin film on a foreign substrate this "quantum energy" of electronic origin can easily overwhelm the strain effects in definging the film stability, thereby severely influencing the preferred growth mode (see, e.g., Suo and Zhang, Phys. Rev. B 58, 5116 (1998)).

Mechanics and deformation of the nucleus in micropipette aspiration experiment

Submitted by Ashkan Vaziri on

Robust biomechanical models are essential for studying the nuclear mechanics and can help shed light on the underlying mechanisms of stress transition in nuclear elements. Here, we develop a computational model for an isolated nucleus undergoing micropipette aspiration. Our model includes distinct components representing the nucleoplasm and the nuclear envelope. The nuclear envelope itself comprises three layers: inner and outer nuclear membranes and one thicker layer representing the nuclear lamina.

Atomistic simulations for the evolution of a U-shaped dislocation in fcc Al

Submitted by Xiaoyan Li on

We show, through MD simulations, a new evolution pattern of the U-shaped dislocation in fcc Al that would enrich the FR mechanism. Direct atomistic investigation indicates that a U-shaped dislocation may behave in different manners when it emits the first dislocation loop by bowing out of an extended dislocation. One manner is that the glissile dislocation segment always bows in the original glide plane, as the conventional FR mechanism. Another is that non-coplanar composite dislocations appear owing to conservative motion of polar dislocation segments, and then bow out along each slip plane, creating a closed helical loop. The motion of these segments involves a cross-slip mechanism by which a dislocation with screw component moves from one slip plane into another. Ultimately, such non-coplanar evolution results in the formation of a FR source.

Some numerical mechanics software

Submitted by Mogadalai Gururajan on

Recently, during one of my net searches, I came across this page of RPI, where I learnt about a couple of numerical mechanics software which might be of interest to some of you.

FMDB:

As for the effort toward the scalable engineering simulations on distributed environements, we addressed this challenge by developing a distributed mesh data management infrastructure that satisfies the needs of distributed domain of applications.

A structure-based sliding-rebinding mechanism for catch bonds

Submitted by Cheng Zhu on

This is a paper by Jizhong Lou and myself, which is in press in Biophysical Journal.

Abstract.  Catch bonds, whose lifetimes are prolonged by force, have been observed in selectin-ligand interactions and other systems. Several biophysical models have been proposed to explain this counter-intuitive phenomenon, but none was based on the structure of the interacting molecules and the noncovalent interactions at the binding interface. Here we used molecular dynamics simulations to study changes in structure and atomic-level interactions during forced unbinding of P-selectin from P-selectin glycoprotein ligand-1. A mechanistic model for catch bonds was developed based on these observations. In the model, "catch" results from forced opening of an interdomain hinge that tilts the binding interface to allow two sides of the contact to slide against each other. Sliding promotes formation of new interactions and even rebinding to the original state, thereby slowing dissociation and prolonging bond lifetimes. Properties of this sliding-rebinding mechanism were explored using a pseudo-atom representation and Monte Carlo simulations. The model has been supported by its ability to fit experimental data and can be related to previously proposed two-pathway models.

How can we obtain more information from protein structure?

Submitted by Cheng Zhu on

We know - or believe - protein function is determined by structure. Crystallographic and NMR studies can provide protein structures with atomic-level details at equilibrium. MD simulations can follow protein conformational changes in time with fs temporal resolution in the absence or presence of a bias mechanism, e.g., applied force, used to induce such changes.