Miniature Compliant Grippers With Vision-Based Force Sensing

This paper is concerned with grasping biological cells in aqueous medium with miniature grippers that can also help estimate forces using vision-based displacement measurement and computation. We present the design, fabrication, and testing of three single-piece, compliant miniature grippers with parallel and angular jaw motions. Two grippers were designed using experience and intuition, while the third one was designed using topology optimization with implicit manufacturing constraints. These grippers were fabricated using different manufacturing techniques using spring steel and polydimethylsiloxane ( PDMS). The grippers also serve the purpose of a force sensor. Toward this, we present a vision-based force-sensing technique by solving Cauchy's problem in elasticity using an improved algorithm. We validated this technique at the macroscale, where there was an independent method to estimate the force. In this study, the gripper was used to hold a yeast ball and a zebrafish egg cell of less than 1 mm in diameter. The forces involved were estimated to be about 30 and 10 mN for the yeast ball and the zebrafish egg cell, respectively.

Frame Soil Interaction And Winkler Model

The Winkler spring model is the most convenient representation of soil support in the domain of linear elasticity for framed structure-soil interaction analyses. The closeness of the analytical results obtained using this model with those corresponding to the elastic half-space continuum has been investigated in the past for foundation beams. The findings, however, are not applicable to framed structures founded on beam or strip footings. Moreover, the past investigations employ the concept of characteristic length which does not adequately account for the stiffness contribution of the superstructure. A framed structure on beam foundation can be described parametrically by the ratios of stiffnesses of superstructure and foundation beams to that of soil. For a practical range of soil allowable pressures, the ranges of these relative stiffness ratios have been established. The present study examines the variation between interactive analyses based on Winkler springs with those using the half-space continuum over these ranges of relative stiffness ratios. The findings enable the analyst to undertake a Winkler spring-based-interaction analysis with knowledge of the likely variation of values with those derived for the more computation-intensive half-space continuum.

Spectral element approach to wave propagation in uncertain beam structures

This paper presents a study on the uncertainty in material parameters of wave propagation responses in metallic beam structures. Special effort is made to quantify the effect of uncertainty in the wave propagation responses at high frequencies. Both the modulus of elasticity and the density are considered uncertain. The analysis is performed using a Monte Carlo simulation (MCS) under the spectral finite element method (SEM). The randomness in the material properties is characterized by three different distributions, the normal, Weibull and extreme value distributions. Their effect on wave propagation in beams is investigated. The numerical study shows that the CPU time taken for MCS under SEM is about 48 times less than for MCS under a conventional one-dimensional finite element environment for 50 kHz loading. The numerical results presented investigate effects of material uncertainties on high frequency modes. A study is performed on the usage of different beam theories and their uncertain responses due to dynamic impulse load. These studies show that even for a small coefficient of variation, significant changes in the above parameters are noticed. A number of interesting results are presented, showing the true effects of uncertainty response due to dynamic impulse load.

Income inequality in the process of economic development : An empirical approach

This thesis studies the effect of income inequality on economic growth. This is done by analyzing panel data from several countries with both short and long time dimensions of the data. Two of the chapters study the direct effect of inequality on growth, and one chapter also looks at the possible indirect effect of inequality on growth by assessing the effect of inequality on savings. In Chapter two, the effect of inequality on growth is studied by using a panel of 70 countries and a new EHII2008 inequality measure. Chapter contributes on two problems that panel econometric studies on the economic effect of inequality have recently encountered: the comparability problem associated with the commonly used Deininger and Squire s Gini index, and the problem relating to the estimation of group-related elasticities in panel data. In this study, a simple way to 'bypass' vagueness related to the use of parametric methods to estimate group-related parameters is presented. The idea is to estimate the group-related elasticities implicitly using a set of group-related instrumental variables. The estimation results with new data and method indicate that the relationship between income inequality and growth is likely to be non-linear. Chapter three incorporates the EHII2.1 inequality measure and a panel with annual time series observations from 38 countries to test the existence of long-run equilibrium relation(s) between inequality and the level of GDP. Panel unit root tests indicate that both the logarithmic EHII2.1 inequality measure and the logarithmic GDP per capita series are I(1) nonstationary processes. They are also found to be cointegrated of order one, which implies that there is a long-run equilibrium relation between them. The long-run growth elasticity of inequality is found to be negative in the middle-income and rich economies, but the results for poor economies are inconclusive. In the fourth Chapter, macroeconomic data on nine developed economies spanning across four decades starting from the year 1960 is used to study the effect of the changes in the top income share to national and private savings. The income share of the top 1 % of population is used as proxy for the distribution of income. The effect of inequality on private savings is found to be positive in the Nordic and Central-European countries, but for the Anglo-Saxon countries the direction of the effect (positive vs. negative) remains somewhat ambiguous. Inequality is found to have an effect national savings only in the Nordic countries, where it is positive.

Direct tests of the entropic model of quasicrystals

Assuming an entropic origin for phason elasticity in quasicrystals, we derive predictions for the temperature dependence of grain-boundary structure and free energy, the nature of the elastic instability in these systems, and the behavior of sound damping near the instability. We believe that these will provide decisive tests of the entropic model for quasicrystals.

Strong nonlocalization induced by small scale parameter on terahertz flexural wave dispersion characteristics of a monolayer graphene

This paper presents the strong nonlocal scale effect on the flexural wave propagation in a monolayer graphene sheet. The graphene is modeled as an isotropic plate of one atom thick. Nonlocal governing equation of motion is derived and wave propagation analysis is performed using spectral analysis. The present analysis shows that the flexural wave dispersion in graphene obtained by local and nonlocal elasticity theories is quite different. The nonlocal elasticity calculation shows that the wavenumber escapes to infinite at certain frequency and the corresponding wave velocity tends to zero at that frequency indicating localization and stationary behavior. This behavior is captured in the spectrum and dispersion curves. The cut-off frequency of flexural wave not only depend on the axial wavenumber but also on the nonlocal scaling parameter. The effect of axial wavenumber on the wave behavior in graphene is also discussed in the present manuscript. (C) 2010 Elsevier B.V. All rights reserved.

Instability and dewetting of ultrathin solid viscoelastic films on homogeneous and heterogeneous substrates

Instability and dewetting engendered by the van der Waals force in soft thin (<100 nm) linear viscoelastic solid (e. g., elastomeric gel) films on uniform and patterned surfaces are explored. Linear stability analysis shows that, although the elasticity of the film controls the onset of instability and the corresponding critical wavelength, the dominant length-scale remains invariant with the elastic modulus of the film. The unstable modes are found to be long-wave, for which a nonlinear long-wave analysis and simulations are performed to uncover the dynamics and morphology of dewetting. The stored elastic energy slows down the temporal growth of instability significantly. The simulations also show that a thermodynamically stable film with zero-frequency elasticity can be made unstable in the presence of physico-chemical defects on the substrate and can follow an entirely different pathway with far fewer holes as compared to the viscous films. Further, the elastic restoring force can retard the growth of a depression adjacent to the hole-rim and thus suppress the formation of satellite holes bordering the primary holes. These findings are in contrast to the dewetting of viscoelastic liquid films where nonzero frequency elasticity accelerates the film rupture and promotes the secondary instabilities. Thus, the zero-frequency elasticity can play a major role in imposing a better-defined long-range order to the dewetted structures by arresting the secondary instabilities. (C) 2011 American Institute of Physics. doi: 10.1063/1.3554748]

3-D Warping in Four-Bar Laminated Linkages

This paper deals with the evaluation of the component-laminate load-carrying capacity, i.e., to calculate the loads that cause the failure of the individual layers and the component-laminate as a whole in four-bar mechanism. The component-laminate load-carrying capacity is evaluated using the Tsai-Wu-Hahn failure criterion for various layups. The reserve factor of each ply in the component-laminate is calculated by using the maximum resultant force and the maximum resultant moment occurring at different time steps at the joints of the mechanism. Here, all component bars of the mechanism are made of fiber reinforced laminates and have thin rectangular cross-sections. They could, in general, be pre-twisted and/or possess initial curvature, either by design or by defect. They are linked to each other by means of revolute joints. We restrict ourselves to linear materials with small strains within each elastic body (beam). Each component of the mechanism is modeled as a beam based on geometrically nonlinear 3-D elasticity theory. The component problems are thus split into 2-D analyses of reference beam cross-sections and nonlinear 1-D analyses along the three beam reference curves. For the thin rectangular cross-sections considered here, the 2-D cross-sectional nonlinearity is also overwhelming. This can be perceived from the fact that such sections constitute a limiting case between thin-walled open and closed sections, thus inviting the nonlinear phenomena observed in both. The strong elastic couplings of anisotropic composite laminates complicate the model further. However, a powerful mathematical tool called the Variational Asymptotic Method (VAM) not only enables such a dimensional reduction, but also provides asymptotically correct analytical solutions to the nonlinear cross-sectional analysis. Such closed-form solutions are used here in conjunction with numerical techniques for the rest of the problem to predict more quickly and accurately than would otherwise be possible. Local 3-D stress, strain and displacement fields for representative sections in the component-bars are recovered, based on the stress resultants from the 1-D global beam analysis. A numerical example is presented which illustrates the failure of each component-laminate and the mechanism as a whole.

Analysis of a cracked lug loaded by an interference fit pin

In the present study, a lug joint fitted with an interference fit (oversized) pin is considered with radial through cracks situated at diametrically opposite points perpendicular to the loading direction. A finite element contact stress algorithm is developed with linear elastic assumptions to deal with varying partial contact/separation at the pin-plate interface using a marching solution. Stress Intensity Factor (SIF) at the crack tips is evaluated using the Modified Crack Closure Integral (MCCI) method. The effect of change in crack length and edge distance on the load-contact relation, SIFs and stress distributions are studied. A rigorous plane stress elasticity solution of the pin-plate interface at the crack mouth confirmed the existence of the stress concentration leading to a local peak in the radial stress at the crack mouth and provided a method of estimating it quantitatively. Copyright (C) 1996 Elsevier Science Ltd.

Analysis of a long thick orthotropic circular cylindrical shell panel

An exact three-dimensional elasticity solution has been obtained for an infinitely long, thick transversely isotropic circular cylindrical shell panel, simply supported along the longitudinal edges and subjected to a radial patch load. Using a set of three displacement functions, the boundary value problem is reduced to Bessel's differential equation. Numerical results are presented for different thickness to mean radius ratios and semicentral angles of the shell panel. Classical and first-order shear deformation orthotropic shell theories have been examined in comparison with the present elasticity solution.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

Stability of the viscous flow of a fluid through a flexible tube

The stability of Hagen-Poiseuille flow of a Newtonian fluid of viscosity eta in a tube of radius R surrounded by a viscoelastic medium of elasticity G and viscosity eta(s) occupying the annulus R < r < HR is determined using a linear stability analysis. The inertia of the fluid and the medium are neglected, and the mass and momentum conservation equations for the fluid and wall are linear. The only coupling between the mean flow and fluctuations enters via an additional term in the boundary condition for the tangential velocity at the interface, due to the discontinuity in the strain rate in the mean flow at the surface. This additional term is responsible for destabilizing the surface when the mean velocity increases beyond a transition value, and the physical mechanism driving the instability is the transfer of energy from the mean flow to the fluctuations due to the work done by the mean flow at the interface. The transition velocity Gamma(t) for the presence of surface instabilities depends on the wavenumber k and three dimensionless parameters: the ratio of the solid and fluid viscosities eta(r) = (eta(s)/eta), the capillary number Lambda = (T/GR) and the ratio of radii H, where T is the surface tension of the interface. For eta(r) = 0 and Lambda = 0, the transition velocity Gamma(t) diverges in the limits k much less than 1 and k much greater than 1, and has a minimum for finite k. The qualitative behaviour of the transition velocity is the same for Lambda > 0 and eta(r) = 0, though there is an increase in Gamma(t) in the limit k much greater than 1. When the viscosity of the surface is non-zero (eta(r) > 0), however, there is a qualitative change in the Gamma(t) vs. k curves. For eta(r) < 1, the transition velocity Gamma(t) is finite only when k is greater than a minimum value k(min), while perturbations with wavenumber k < k(min) are stable even for Gamma--> infinity. For eta(r) > 1, Gamma(t) is finite only for k(min) < k < k(max), while perturbations with wavenumber k < k(min) or k > k(max) are stable in the limit Gamma--> infinity. As H decreases or eta(r) increases, the difference k(max)- k(min) decreases. At minimum value H = H-min, which is a function of eta(r), the difference k(max)-k(min) = 0, and for H < H-min, perturbations of all wavenumbers are stable even in the limit Gamma--> infinity. The calculations indicate that H-min shows a strong divergence proportional to exp (0.0832 eta(r)(2)) for eta(r) much greater than 1.

Lamination-dependent shear deformation models for bending of angle-ply laminated plates

Lamination-dependent shear corrective terms in the analysis of bending of laminated plates are derived from a priori assumed linear thicknesswise distributions for gradients of transverse shear stresses by using CLPT inplane stresses in the two in-plane equilibrium equations of elasticity in each ply. In the development of a general model for angle-ply laminated plates, special cases like cylindrical bending of laminates in either direction, symmetric laminates, cross-ply laminates, antisymmetric angle-ply laminates, homogeneous plates are taken into consideration. Adding these corrective terms to the assumed displacements in (i) Classical Laminate Plate Theory (CLPT) and (ii) Classical Laminate Shear Deformation Theory (CLSDT), two new refined lamination-dependent shear deformation models are developed. Closed form solutions from these models are obtained for antisymmetric angle-ply laminates under sinusoidal load for a type of simply supported boundary conditions. Results obtained from the present models and also from Ren's model (1987) are compared with each other.

Lamination-dependent shear deformation models for cylindrical bending of angle-ply laminates

Lamination-dependent shear corrective terms in the analysis of flexure of laminates are derived from a priori assumed linear thicknesswise distributions for gradients of transverse shear stresses and using them in the two in-plane equilibrium equations of elasticity in each ply. Adding these corrective terms to (i) Classical Laminate Plate Theory (CLPT) displacements and (ii) Classical Laminate Shear Deformation Theory (CLSDT) displacements, four new refined lamination-dependent shear deformation models for angle-ply laminates are developed. Performance of these models is evaluated by comparing the results from these models with those from exact elasticity solutions for antisymmetric 2-ply laminates and for 4-ply [15/-15](s) laminates. In general, the model with shear corrective terms based on CLPT and added to CLSDT displacements is sufficient and predicts good estimates, both qualitatively and quantitatively, for all displacements and stresses.

The Equations of Equilibrium in Orthogonal Curvilinear Reference Coordinates

Analytical solutions to problems in finite elasticity are most often derived using the semi-inverse approach along with the spatial form of the equations of motion involving the Cauchy stress tensor. This procedure is somewhat indirect since the spatial equations involve derivatives with respect to spatial coordinates while the unknown functions are in terms of material coordinates, thus necessitating the use of the chain rule. In this classroom note, we derive compact expressions for the components of the divergence, with respect to orthogonal material coordinates, of the first Piola-Kirchhoff stress tensor. The spatial coordinate system is also assumed to be an orthogonal curvilinear one, although, not necessarily of the same type as the material coordinate system. We show by means of some example applications how analytical solutions can be derived more directly using the derived results.