Sliding Modes for the Simulation of Mechanical and Electrical Systems Defined by Differential-Algebraic Equations

We offer a technique, motivated by feedback control and specifically sliding mode control, for the simulation of differential-algebraic equations (DAEs) that describe common engineering systems such as constrained multibody mechanical structures and electric networks. Our algorithm exploits the basic results from sliding mode control theory to establish a simulation environment that then requires only the most primitive of numerical solvers. We circumvent the most important requisite for the conventionalsimulation of DAEs: the calculation of a set of consistent initial conditions. Our algorithm, which relies on the enforcement and occurrence of sliding mode, will ensure that the algebraic equation is satisfied by the dynamic system even for inconsistent initial conditions and for all time thereafter. [DOI:10.1115/1.4001904]

Mechanical Properties of Al-Mica Particulate Composite Material

Cast aluminium alloy mica particle composites of varying mica content were tested in tension, compression, and impact. With 2.2 percent mica (size range 40µm – 120µm) the tensile and compression strengths of aluminium alloy decreased by 56 and 22 percent, respectively. The corresponding decreases in percent elongation and percent reduction are 49 and 39 percent. Previous work [2] shows that despite this decrease in strength the composite with 2.5 percent mica and having an UTS of 15 kg/mm2 and compression strength of 28 kg/mm2 performs well as a bearing material under severe running conditions. The differences in strength characteristics of cast aluminium-mica particle composites between tension and compression suggests that, as in cast iron, expansion of voids at the matrix particle interface may be the guiding mechanism of the deformation. SEM studies show that on the tensile fractured specimen surface, there are large voids at the particle matrix interface.

COMPOSITE AND MULTI-SCALE COMPLIANT MECHANISMS FOR MANIPULATION AND MECHANICAL CHARACTERIZATION

In this paper, we integrate two or more compliant mechanisms to get enhanced functionality for manipulating and mechanically characterizing the grasped objects of varied size (cm to sub-mm), stiffness (1e5 to 10 N/m), and materials (cement to biological cells). The concepts of spring-lever (SL) model, stiffness maps, and non-dimensional kinetoelastostatic maps are used to design composite and multi-scale compliant mechanisms. Composite compliant mechanisms comprise two or more different mechanisms within a single elastic continuum while multi-scale ones possess the additional feature of substantial difference in the sizes of the mechanisms that are combined into one. We present three applications: (i) a composite compliant device to measure the failure load of the cement samples; (ii) a composite multi-scale compliant gripper to measure the bulk stiffness of zebrafish embryos; and (iii) a compliant gripper combined with a negative-stiffness element to reduce the overall stiffness. The prototypes of all three devices are made and tested. The cement sample needed a breaking force of 22.5 N; the zebrafish embryo is found to have bulk stiffness of about 10 N/m; and the stiffness of a compliant gripper was reduced by 99.8 % to 0.2 N/m.

Effect Of Directionality Of Grinding Marks On Friction And Formation Of Transfer Layer

Surface texture influences friction and transfer layer formation during sliding. In the present investigation, basic studies were conducted using inclined pin-on-plate sliding tester to understand the effect of directionality of surface grinding marks of hard material on friction and transfer layer formation during sliding against soft materials. 080 M40 steel plates were ground to attain different surface roughness with unidirectional grinding marks. Then pins made of soft materials such as pure Al, pure Mg and Al-Mg alloy were slid against the prepared steel plates. Grinding angle (i.e., the angle between direction of sliding and grinding marks) was varied between 0 degrees and 90 degrees in the tests. Experiments were conducted under both dry and lubricated conditions on each plate in ambient environment. It was observed that the transfer layer formation and the coefficient of friction, which has two components adhesion and plowing - depend primarily on the directionality of grinding marks of the harder mating surface, and independent of surface roughness of the harder mating surface. For the case of pure Mg, stick-slip phenomenon was observed under dry condition for all grinding angles and it was absent upto 20 degrees grinding angles under lubricated condition. However, for the case of Al, it was observed only under lubricated conditions for angles exceeding 20 degrees. As regards the alloy, namely, Al-Mg alloy, it, was absent in both conditions. For the case of pure Mg and Al, it was observed that the amplitude of stick-slip motion primarily depends on plowing component of friction. The grinding angle effect on coefficient of friction was attributed to the variation of plowing component of friction with grinding angle.

Studies On Friction And Formation Of Transfer Layer In Hcp Metals

Surface texture plays an important role in the frictional behavior and transfer layer formation of contacting surfaces. In the present investigation, basic experiments were conducted using an inclined pin-on-plate sliding apparatus to better understand the role of surface texture on the coefficient of friction and the formation of a transfer layer. In the experiments, soft HCP materials such as pure Mg and pure Zn were used for the pins and a hardened 080 M40 steel was used for the plate. Two surface parameters of the steel plates—roughness and texture—were varied in tests that were conducted at a sliding speed of 2 mm/s in ambient conditions under both dry and lubricated conditions. The morphologies of the worn surfaces of the pins and the formation of the transfer layer on the counter surfaces were observed using a scanning electron microscope. In the experiments, the occurrence of stick-slip motion, the formation of a transfer layer, and the value of friction were recorded. With respect to the friction, both adhesion and plowing components were analyzed. Based on the experimental results, the effect of surface texture on the friction was attributed to differences in the amount of plowing. Both the plowing component of friction and the amplitude of stick-slip motion were determined to increase surface textures that promote plane strain conditions and decrease the textures that favor plane stress conditions.

Correlating The Features Of Topography To Friction By Sliding Experiments

Friction can influence the quality of the finished product to a large extent in certain manufacturing processes. Sheet metal forming is a particular case, where the friction between the hard-die and the relatively soft work-piece can be extremely important. Under such conditions, topography of the harder surface can influence the resistance to traction at the interface. This paper discusses about the correlation between certain features of the surface; topography and coefficient of friction based on experiments involving sliding of a few soft metal pins against a harder material. A brief description of the experimental procedure and the analysis are presented. A hybrid parameter which encapsulates both the amplitude features as well as the relative packing of peaks is shown to correlate well with the coefficient of friction.

Fatigue Laws Via Functional Equations

In routine industrial design, fatigue life estimation is largely based on S-N curves and ad hoc cycle counting algorithms used with Miner's rule for predicting life under complex loading. However, there are well known deficiencies of the conventional approach. Of the many cumulative damage rules that have been proposed, Manson's Double Linear Damage Rule (DLDR) has been the most successful. Here we follow up, through comparisons with experimental data from many sources, on a new approach to empirical fatigue life estimation (A Constructive Empirical Theory for Metal Fatigue Under Block Cyclic Loading', Proceedings of the Royal Society A, in press). The basic modeling approach is first described: it depends on enforcing mathematical consistency between predictions of simple empirical models that include indeterminate functional forms, and published fatigue data from handbooks. This consistency is enforced through setting up and (with luck) solving a functional equation with three independent variables and six unknown functions. The model, after eliminating or identifying various parameters, retains three fitted parameters; for the experimental data available, one of these may be set to zero. On comparison against data from several different sources, with two fitted parameters, we find that our model works about as well as the DLDR and much better than Miner's rule. We finally discuss some ways in which the model might be used, beyond the scope of the DLDR.

Axially Loaded Timoshenko Rotors From A Prestressed Continuum Approach

We consider an axially loaded Timoshenko rotor rotating at a constant speed and derive its governing equations from a continuum viewpoint. The primary aim of this paper is to understand the source and role of gyroscopic terms, when the rotor is viewed not as a Timoshenko beam but as a genuine 3D continuum. We offer the primary insight that macroscopically observed gyroscopic terms may also, quite equivalently, be viewed as external manifestations of internally existing spin-induced prestresses at the continuum level. To demonstrate this idea with an analytical example (the Timoshenko rotor), we have studied the reliable equations of Choi et al. (Journal of Vibration and Acoustics, 114, 1992, 249-259). Using a straightforward application of our insight in the framework of nonlinear elasticity, we obtain equations that exactly match Choi et al. for the case with no axial load. For the case of axial preload, our straightforward formulation leads to a slightly different set of equations that have negligible numerical consequence for solid rotors. However, we offer a macroscopic, intuitive, justification for modifying our formulation so as to obtain the exact equations of Choi et al. with the axial load included.

Infinite Dimensional Slow Modulations In A Delayed Model For Orthogonal Cutting Vibrations

We apply the method of multiple scales (MMS) to a well known model of regenerative cutting vibrations in the large delay regime. By ``large'' we mean the delay is much larger than the time scale of typical cutting tool oscillations. The MMS upto second order for such systems has been developed recently, and is applied here to study tool dynamics in the large delay regime. The second order analysis is found to be much more accurate than first order analysis. Numerical integration of the MMS slow flow is much faster than for the original equation, yet shows excellent accuracy. The main advantage of the present analysis is that infinite dimensional dynamics is retained in the slow flow, while the more usual center manifold reduction gives a planar phase space. Lower-dimensional dynamical features, such as Hopf bifurcations and families of periodic solutions, are also captured by the MMS. Finally, the strong sensitivity of the dynamics to small changes in parameter values is seen clearly.

A Global Constraint On Relative Rotation To Avoid Lumped Compliant Mechanisms In Topology Optimization

Some of the well known formulations for topology optimization of compliant mechanisms could lead to lumped compliant mechanisms. In lumped compliance, most of the elastic deformation in a mechanism occurs at few points, while rest of the mechanism remains more or less rigid. Such points are referred to as point-flexures. It has been noted in literature that high relative rotation is associated with point-flexures. In literature we also find a formulation of local constraint on relative rotations to avoid lumped compliance. However it is well known that a global constraint is easier to handle than a local constraint, by a numerical optimization algorithm. The current work presents a way of putting global constraint on relative rotations. This constraint is also simpler to implement since it uses linearized rotation at the center of finite-elements, to compute relative rotations. I show the results obtained by using this constraint oil the following benchmark problems - displacement inverter and gripper.

A Jacobian-Based Algorithm for Planning Attitude Maneuvers Using Forward and Reverse Rotations

Algorithms for planning quasistatic attitude maneuvers based on the Jacobian of the forward kinematic mapping of fully-reversed (FR) sequences of rotations are proposed in this paper. An FR sequence of rotations is a series of finite rotations that consists of initial rotations about the axes of a body-fixed coordinate frame and subsequent rotations that undo these initial rotations. Unlike the Jacobian of conventional systems such as a robot manipulator, the Jacobian of the system manipulated through FR rotations is a null matrix at the identity, which leads to a total breakdown of the traditional Jacobian formulation. Therefore, the Jacobian algorithm is reformulated and implemented so as to synthesize an FR sequence for a desired rotational displacement. The Jacobian-based algorithm presented in this paper identifies particular six-rotation FR sequences that synthesize desired orientations. We developed the single-step and the multiple-step Jacobian methods to accomplish a given task using six-rotation FR sequences. The single-step Jacobian method identifies a specific FR sequence for a given desired orientation and the multiple-step Jacobian algorithm synthesizes physically feasible FR rotations on an optimal path. A comparison with existing algorithms verifies the fast convergence ability of the Jacobian-based algorithm. Unlike closed-form solutions to the inverse kinematics problem, the Jacobian-based algorithm determines the most efficient FR sequence that yields a desired rotational displacement through a simple and inexpensive numerical calculation. The procedure presented here is useful for those motion planning problems wherein the Jacobian is singular or null.

Approximate closed-form solutions of realistic true proportional navigation guidance using the Adomian decomposition method

Approximate closed-form solutions of the non-linear relative equations of motion of an interceptor pursuing a target under the realistic true proportional navigation (RTPN) guidance law are derived using the Adomian decomposition method in this article. In the literature, no study has been reported on derivation of explicit time-series solutions in closed form of the nonlinear dynamic engagement equations under the RTPN guidance. The Adomian method provides an analytical approximation, requiring no linearization or direct integration of the non-linear terms. The complete derivation of the Adomian polynomials for the analysis of the dynamics of engagement under RTPN guidance is presented for deterministic ideal case, and non-ideal dynamics in the loop that comprises autopilot and actuator dynamics and target manoeuvre, as well as, for a stochastic case. Numerical results illustrate the applicability of the method.

Rotating Beams and Nonrotating Beams With Shared Eigenpair

In this paper, we look for rotating beams whose eigenpair (frequency and mode-shape) is the same as that of uniform nonrotating beams for a particular mode. It is found that, for any given mode, there exist flexural stiffness functions (FSFs) for which the jth mode eigenpair of a rotating beam, with uniform mass distribution, is identical to that of a corresponding nonrotating uniform beam with the same length and mass distribution. By putting the derived FSF in the finite element analysis of a rotating cantilever beam, the frequencies and mode-shapes of a nonrotating cantilever beam are obtained. For the first mode, a physically feasible equivalent rotating beam exists, but for higher modes, the flexural stiffness has internal singularities. Strategies for addressing the singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test-functions for rotating beam codes and for targeted destiffening of rotating beams.

Vibrations of a Beam in Variable Contact With a Flat Surface

We study small vibrations of cantilever beams contacting a rigid surface. We study two cases: the first is a beam that sags onto the ground due to gravity, and the second is a beam that sticks to the ground through reversible adhesion. In both cases, the noncontacting length varies dynamically. We first obtain the governing equations and boundary conditions, including a transversality condition involving an end moment, using Hamilton's principle. Rescaling the variable length to a constant value, we obtain partial differential equations with time varying coefficients, which, upon linearization, give the natural frequencies of vibration. The natural frequencies for the first case (gravity without adhesion) match that of a clamped-clamped beam of the same nominal length; frequencies for the second case, however, show no such match. We develop simple, if atypical, single degree of freedom approximations for the first modes of these two systems, which provide insights into the role of the static deflection profile, as well as the end moment condition, in determining the first natural frequencies of these systems. Finally, we consider small transverse sinusoidal forcing of the first case and find that the governing equation contains both parametric and external forcing terms. For forcing at resonance, w find that either the internal or the external forcing may dominate.

Active vibration suppression of non-linear beams using optimal dynamic inversion

Euler–Bernoulli beams are distributed parameter systems that are governed by a non-linear partial differential equation (PDE) of motion. This paper presents a vibration control approach for such beams that directly utilizes the non-linear PDE of motion, and hence, it is free from approximation errors (such as model reduction, linearization etc.). Two state feedback controllers are presented based on a newly developed optimal dynamic inversion technique which leads to closed-form solutions for the control variable. In one formulation a continuous controller structure is assumed in the spatial domain, whereas in the other approach it is assumed that the control force is applied through a finite number of discrete actuators located at predefined discrete locations in the spatial domain. An implicit finite difference technique with unconditional stability has been used to solve the PDE with control actions. Numerical simulation studies show that the beam vibration can effectively be decreased using either of the two formulations.