Fabrication and Nonlinear Thermomechanical Analysis of SU8 Thermal Actuator

SU8-based micromechanical structures are widely used as thermal actuators in the development of compliant micromanipulation tools. This paper reports the design, nonlinear thermomechanical analysis, fabrication, and thermal actuation of SU8 actuators. The thermomechanical analysis of the actuator incorporates nonlinear temperature-dependent properties of SU8 polymer to accurately model its thermal response during actuation. The designed SU8 thermal actuators are fabricated using surface micromachining techniques and the electrical interconnects are made to them using flip-chip bonding. The issues due to thermal stress during fabrication are discussed and a novel strategy is proposed to release the thermal stress in the fabricated actuators. Subsequent characterization of the actuator using an optical profilometer reveals excellent thermal response, good repeatability, and low hysteresis. The average deflection is similar to 8.5 mu m for an actuation current of similar to 5 mA. The experimentally obtained deflection profile and the tip deflection at different currents are both shown to be in good agreement with the predictions of the nonlinear thermomechanical model. This underscores the need to consider nonlinearities when modeling the response of SU8 thermal actuators. 2015-0087]

Optimal Guidance for Accurate Lunar Soft Landing with Minimum Fuel Consumption using Model Predictive Static Programming

In this paper the soft lunar landing with minimum fuel expenditure is formulated as a nonlinear optimal guidance problem. The realization of pinpoint soft landing with terminal velocity and position constraints is achieved using Model Predictive Static Programming (MPSP). The high accuracy of the terminal conditions is ensured as the formulation of the MPSP inherently poses final conditions as a set of hard constraints. The computational efficiency and fast convergence make the MPSP preferable for fixed final time onboard optimal guidance algorithm. It has also been observed that the minimum fuel requirement strongly depends on the choice of the final time (a critical point that is not given due importance in many literature). Hence, to optimally select the final time, a neural network is used to learn the mapping between various initial conditions in the domain of interest and the corresponding optimal flight time. To generate the training data set, the optimal final time is computed offline using a gradient based optimization technique. The effectiveness of the proposed method is demonstrated with rigorous simulation results.

Traveling waves in trimer granular lattice I: Bifurcation structure of traveling waves in the unit-cell model

Present paper is the first one in the series devoted to the dynamics of traveling waves emerging in the uncompressed, tri-atomic granular crystals. This work is primarily concerned with the dynamics of one-dimensional periodic granular trimer (tri-atomic) chains in the state of acoustic vacuum. Each unit cell consists of three spherical particles of different masses subject to periodic boundary conditions. Hertzian interaction law governs the mutual interaction of these particles. Under the assumption of zero pre-compression, this interaction is modeled as purely nonlinear, which means the absence of linear force component. The dynamics of such chains is governed by the two system parameters that scale the mass ratios between the particles of the unit cell. Such a system supports two different classes of periodic solutions namely the traveling and standing waves. The primary objective of the present study is the numerical analysis of the bifurcation structure of these solutions with emphasis on the dynamics of traveling waves. In fact, understanding of the bifurcation structure of the traveling wave solutions emerging in the unit-cell granular trimer is rather important and can shed light on the more complex nonlinear wave phenomena emerging in semi-infinite trimer chains. (c) 2016 Elsevier B.V. All rights reserved.

Identification of Nonlinear Systems Through Time-Frequency Filtering Technique

In the previous paper, a class of nonlinear system is mapped to a so-called skeleton linear model (SLM) based on the joint time-frequency analysis method. Behavior of the nonlinear system may be indicated quantitatively by the variance of the coefficients of SLM versus its response. Using this model we propose an identification method for nonlinear systems based on nonstationary vibration data in this paper. The key technique in the identification procedure is a time-frequency filtering method by which solution of the SLM is extracted from the response data of the corresponding nonlinear system. Two time-frequency filtering methods are discussed here. One is based on the quadratic time-frequency distribution and its inverse transform, the other is based on the quadratic time-frequency distribution and the wavelet transform. Both numerical examples and an experimental application are given to illustrate the validity of the technique.

Modeling Nonlinear Peeling of Ductile Thin Films-Critical Assessment of Analytical Bending Models Using FE Simulations

Three analytical double-parameter criteria based on a bending model and a two-dimensional finite element analysis model are presented for the modeling of ductile thin film undergoing a nonlinear peeling process. The bending model is based on different governing parameters: (1) the interfacial fracture toughness and the separation strength, (2) the interfacial fracture toughness and the crack tip slope angle, and (3) the interfacial fracture toughness and the critical Mises effective strain of the delaminated thin film at the crack tip. Thin film nonlinear peeling under steady-state condition is solved with the different governing parameters. In addition, the peeling test problem is simulated by using the elastic-plastic finite element analysis model. A critical assessment of the three analytical bending models is made by comparison of the bending model solutions with the finite element analysis model solutions. Furthermore, through analyses and comparisons for solutions based on both the bending model and the finite element analysis model, some connections between the bending model and the finite element analysis model are developed. Moreover, in the present research, the effect of different selections for cohesive zone shape on the ductile film peeling solutions is discussed.

Stochastic Structural Model of Rock and Soil Aggregates by Continuum-Based Discrete Element Method

This paper first presents a stochastic structural model to describe the random geometrical features of rock and soil aggregates. The stochastic structural model uses mixture ratio, rock size and rock shape to construct the microstructures of aggregates,and introduces two types of structural elements (block element and jointed element) and three types of material elements (rock element, soil element, and weaker jointed element)for this microstructure. Then, continuum-based discrete element method is used to study the deformation and failure mechanism of rock and soil aggregate through a series of loading tests. It is found that the stress-strain curve of rock and soil aggregates is nonlinear, and the failure is usually initialized from weaker jointed elements. Finally, some factors such as mixture ratio, rock size and rock shape are studied in detail. The numerical results are in good agreement with in situ test. Therefore, current model is effective for simulating the mechanical behaviors of rock and soil aggregates.

High Aspect Ratio (L/D) Riser Viv Prediction Using Wake Oscillator Model

A two-dimensional (2-D) vortex-induced vibration (VIV) prediction model for high aspect ratio (LID) riser subjected to uniform and sheared flow is studied in this paper. The nonlinear structure equations are considered. The near wake dynamics describing the fluctuating nature of vortex shedding is modeled using classical van der Pol equation. A new approach was applied to calibrate the empirical parameters in the wake oscillator model. Compared the predicted results with the experimental data and computational fluid dynamic (CFD) results. Good agreements are observed. It can be concluded that the present model can be used as simple computational tool in predicting some aspects of VIV of long flexible structures. (C) 2008 Elsevier Ltd. All rights reserved.

Nonlinear Constitutive Behavior Of Ferroelectric Materials

The ferroelectric specimen is considered as an aggregation of many randomly oriented domains. According to this mechanism, a multi-domain mechanical model is developed in this paper. Each domain is represented by one element. The applied stress and electric field are taken to be the stress and electric field in the formula of the driving force of domain switching for each element in the specimen. It means that the macroscopic switching criterion is used for calculating the volume fraction of domain switching for each element. By using the hardening relation between the driving force of domain switching and the volume fraction of domain switching calibrated, the volume fraction of domain switching for each element is calculated. Substituting the stress and electric field and the volume fraction of domain switching into the constitutive equation of ferroelectric material, one can easily get the strain and electric displacement for each element. The macroscopic behavior of the ferroelectric specimen is then directly calculated by volume averaging. Meanwhile, the nonlinear finite element analysis for the ferroelectric specimen is carried out. In the finite element simulation, the volume fraction of domain switching for each element is calculated by using the same method mentioned above. The interaction between different elements is taken into account in the finite element simulation and the local stress and electric field for each element is obtained. The macroscopic behavior of the specimen is then calculated by volume averaging. The computation results involve the electric butterfly shaped curves of axial strain versus the axial electric field and the hysteresis loops of electric displacement versus the electric field for ferroelectric specimens under the uniaxial coupled stress and electric field loading. The present theoretical prediction agrees reasonably with the experimental results.

Energy Principle And Nonlinear Electric-Mechanical Behavior Of Ferroelectric Ceramics

Many experimental observations have shown that a single domain in a ferroelectric material switches by progressive movement of domain walls, driven by a combination of electric field and stress. The mechanism of the domain switch involves the following steps: initially, the domain has a uniform spontaneous polarization; new domains with the reverse polarization direction nucleate, mainly at the surface, and grow though the crystal thickness; the new domain expands sideways as a new domain continues to form; finally, the domain switch coalesces to complete the polarization reversal. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of the ferroelectric material and used to study the nonlinear constitutive behavior of a ferroelectric body in this paper. The principle of stationary total potential energy is put forward in which the basic unknown quantities are the displacement u(i), electric displacement D-i and volume fraction rho(I) of the domain switching for the variant I. The mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total potential energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion established by Hwang et al. [ 1]. Based on the domain switching criterion, a set of linear algebraic equations for determining the volume fraction rho(I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. If the volume fraction rho(I) of domain switching for each domain is prescribed, the unknown displacement and electric potential can be obtained based on the conventional finite element procedure. It is assumed that a domain switches if the reduction in potential energy exceeds a critical energy barrier. According to the experimental results, the energy barrier will strengthen when the volume fraction of the domain switching increases. The external mechanical and electric loads are increased step by step. The volume fraction rho(I) of domain switching for each element obtained from the last loading step is used as input to the constitutive equations. Then the strain and electric fields are calculated based on the conventional finite element procedure. The finite element analysis is carried out on the specimens subjected to uniaxial coupling stress and electric field. Numerical results and available experimental data are compared and discussed. The present theoretic prediction agrees reasonably with the experimental results.

An Efficient Front-Tracking Method For Fully Nonlinear Interfacial Waves

A fully nonlinear and dispersive model within the framework of potential theory is developed for interfacial (2-layer) waves. To circumvent the difficulties arisen from the moving boundary problem a viable technique based on the mixed Eulerian and Lagrangian concept is proposed: the computing area is partitioned by a moving mesh system which adjusts its location vertically to conform to the shape of the moving boundaries but keeps frozen in the horizontal direction. Accordingly, a modified dynamic condition is required to properly compute the boundary potentials. To demonstrate the effectiveness of the current method, two important problems for the interfacial wave dynamics, the generation and evolution processes, are investigated. Firstly, analytical solutions for the interfacial wave generations by the interaction between the barotropic tide and topography are derived and compared favorably with the numerical results. Furthermore simulations are performed for the nonlinear interfacial wave evolutions at various water depth ratios and satisfactory agreement is achieved with the existing asymptotical theories. (c) 2008 Elsevier Inc. All rights reserved.

Pull-In Instability Of Micro-Switch Actuators: Model Review

The existing three widely used pull-in theoretical models (i.e., one-dimensional lumped model, linear supposition model and planar model) are compared with the nonlinear beam mode in this paper by considering both cantilever and fixed-fixed type micro and nano-switches. It is found that the error of the pull-in parameters between one-dimensional lumped model and the nonlinear beam model is large because the denominator of the electrostatic force is minimal when the electrostatic force is computed at the maximum deflection along the beam. Since both the linear superposition model and the slender planar model consider the variation of electrostatic force with the beam's deflection, these two models not only are of the same type but also own little error of the pull-in parameters with the nonlinear beam model, the error brought by these two models attributes to that the boundary conditions are not completely satisfied when computing the numerical integration of the deflection.

Energy principle of ferroelectric ceramics and single domain mechanical model

Many physical experiments have shown that the domain switching in a ferroelectric material is a complicated evolution process of the domain wall with the variation of stress and electric field. According to this mechanism, the volume fraction of the domain switching is introduced in the constitutive law of ferroelectric ceramic and used to study the nonlinear constitutive behavior of ferroelectric body in this paper. The principle of stationary total energy is put forward in which the basic unknown quantities are the displacement u (i) , electric displacement D (i) and volume fraction rho (I) of the domain switching for the variant I. Mechanical field equation and a new domain switching criterion are obtained from the principle of stationary total energy. The domain switching criterion proposed in this paper is an expansion and development of the energy criterion. On the basis of the domain switching criterion, a set of linear algebraic equations for the volume fraction rho (I) of domain switching is obtained, in which the coefficients of the linear algebraic equations only contain the unknown strain and electric fields. Then a single domain mechanical model is proposed in this paper. The poled ferroelectric specimen is considered as a transversely isotropic single domain. By using the partial experimental results, the hardening relation between the driving force of domain switching and the volume fraction of domain switching can be calibrated. Then the electromechanical response can be calculated on the basis of the calibrated hardening relation. The results involve the electric butterfly shaped curves of axial strain versus axial electric field, the hysteresis loops of electric displacement versus electric filed and the evolution process of the domain switching in the ferroelectric specimens under uniaxial coupled stress and electric field loading. The present theoretic prediction agrees reasonably with the experimental results given by Lynch.

Nonlinear dynamics of atomic force microscopy with intermittent contact

When the atomic force microscopy (AFM) in tapping mode is in intermittent contact with a soft substrate, the contact time can be a significant portion of a cycle, resulting in invalidity of the impact oscillator model, where the contact time is assumed to be infinitely small. Furthermore, we demonstrate that the AFM intermittent contact with soft substrate can induce the motion of higher modes in the AFM dynamic response. Traditional ways of modeling AFM (one degree of freedom (DOF) system or single mode analysis) are shown to have serious mistakes when applied to this kind of problem. A more reasonable displacement criterion on contact is proposed, where the contact time is a function of the mechanical properties of AFM and substrate, driving frequencies/amplitude, initial conditions, etc. Multi-modal analysis is presented and mode coupling is also shown. (c) 2006 Published by Elsevier Ltd.

Nonlinear analysis of pile foundation jacket platform

In this paper, the nonlinear collapse of the BOHAI-8 pile foundation jacket platform has been analyzed. The ultimate load and collapse process of two computational models of the structure are given. One model is of fixed support whose length is eight times the pile leg diameter and the other considers the nonlinearity of the soil-pile interaction.

A unified model for dislocation nucleation, dislocation emission and dislocation free zone

In this paper, a unified model for dislocation nucleation, emission and dislocation free zone is proposed based on the Peierls framework. Three regions are identified ahead of the crack tip. The emitted dislocations, located away from the crack tip in the form of an inverse pileup, define the plastic zone. Between that zone and the cohesive zone immediately ahead of the crack tip, there is a dislocation free zone. With the stress field and the dislocation density field in the cohesive zone and plastic zone being, respectively, expressed in the first and second Chebyshev polynomial series, and the opening and slip displacements in trigonometric series, a set of nonlinear algebraic equations can be obtained and solved with the Newton-Raphson Method. The results of calculations for pure shearing and combined tension and shear loading after dislocation emission are given in detail. An approximate treatment of the dynamic effects of the dislocation emission is also developed in this paper, and the calculation results are in good agreement with those of molecular dynamics simulations.