Matrix Shear-Lag Parameter in a Shape Memory Alloy-Actuator-Reinforced Silicon Elastomer
Michael. O. Obaji
Department of Precision Machinery and Instrumentation, University of Science and Technology of China, Hefei, China
To cite this article:
Michael. O. Obaji. Matrix Shear-Lag Parameter in a Shape Memory Alloy-Actuator-Reinforced Silicon Elastomer. International Journal of Materials Science and Applications. Vol. 5, No. 5, 2016, pp. 194-201. doi: 10.11648/j.ijmsa.20160505.13
Received: August 11, 2016; Accepted: August 22, 2016; Published: September 9, 2016
Abstract: The excellent qualities possessed by silicon elastomer when used in designing flexible parts of mechanical systems, make it imperative that we analyze its deformations and distributed axial forces in a matrix of shape memory alloy (SMA) fibers designed as possible appendages for gripping robots. The essence of these analyses is to determine the shear-lag parameter which has influence on the axial distributed forces proposed as a gauge for testing the structure in a high yield, high force and high strain mechanical environment. The insertion of SMA fibers in flexible rods cast using silicon elastomer results in the deformation of the host medium once shape recovery of the fiber occurs. This paper aims to analyze the mechanics of the said shape recovery in a modeled silicon elastomer rod with a single off-axis reinforced SMA actuator. The compressive force distribution mechanism and the bending moment caused by phase transformation in the design are determined using an approximate analytical model. The deformations on the structure proposed as an appendage on gripping robots were further analyzed by determining their equations of equilibrium, force factors and their comparative shear-lag models to be able to estimate the force distribution on the structure.
Keywords: SMA Actuator, Mechanics, Silicon Elastomer, Deformation, Soft Robots
The phase changes often associated with SMA fibers via heating and cooling cycles is usually as a result of electrical heating from martensite to austenite. This is what causes the bending of the SMA fibers . And because resistive heating generated by an electric current is the major driving force for their actuation, it thus becomes somewhat easier to drive the fibers, especially when embedded in a rod . SMA fibers are known to exert high force and allow for high strain, thus making their insertion in silicon elastomer rods an interesting area of research, especially in designing mechanical systems.
Silicone elastomer is a polymeric material with visco-elasticity made of poly-siloxanes obtained from ketones and monomers such as carbon, hydrogen and oxygen. They have low Young’s modulus and high yield strain compared with any other materials . When cast at ambient temperatures, silicon elastomers are relatively soft and deformable. They are mainly used for molding flexible parts in aircrafts, seals in building construction & maintenance and adhesives in cars. They are also widely used in producing medical grade materials for treating patients in hospitals. Silicon elastomers are thermosets by nature, which means they are vulcanizable but could also be thermoplastic. Their molecular structure look like the spaghetti and meatball cooked as food in our homes. The meatball serves as a cross-link (or chain) which gives it the elastic property. This elastic property comes from the ability of the long chains to reconfigure themselves and distribute applied stress [Keller et al]. The covalent chains (or cross-linkages) ensure that the elastomer returns to its original configuration when the stress is removed. Thus, making its shape-recovery features similar to that of the SMA fiber. Although, the insertion of SMA fibers in silicon elastomer for the purpose of developing appendages for mechanical systems such as robots is still a novel idea being researched on, there is need to emphasize that some breakthroughs similar to this had been made in designing appendages for robots, although not with silicon elastomer, as Hitachi produced a 4-fingered robotic hand using 12 groups of 0.2mm SMA fibers that closed their hands upon activation. Similarly, there was an SMA-Based Six-Legged Biomimetic Robot called SMABOT built by the National Taiwan University. There was also a 4-Degree of Freedom Origami-enabled, SMA-actuated Robotic End-effectors for Minimally Invasive Surgery designed in King’s College London. All these feats so far achieved encourage more researchers to exploit the possibility of using elastomeric materials reinforced with SMA actuators in building components for mechanical systems such as robots. An accurate analysis of the force distribution and shear-lag parameter on the modeled structure which could help determine its effectiveness when designing mechanical systems is a case for concern, and this has to be investigated. This paper will also investigate the nagging issue of how to determine SMA fiber placement in a rod of a given shape.
2. Design and Deformation
The embedding of SMA fiber in a flexible rod is complex due to force transfer mechanism from the actuator to the rod [2,11]. However, in order to deal with this complexity, we have to consider a single SMA fiber embedded in a flexible rod as shown in figure 1(a).
Further to this, a suitable model for the force transfer from the SMA fiber to the cylindrical host medium is the cylindrical shear-lag model . According to this model, the axial strain of the cylinder and the partial shape recovery of SMA occur by the deformation of a shear layer around the SMA fiber.
In this paper, we shall establish the equilibrium equations for the shear layer, the host medium, and the SMA fiber.
The cylinder cast using silicon elastomer is intended to imitate the appendage of an octopus or that of an octopus robot.
3. Mechanism and Mathematical Analysis
During force transfer from SMA fiber to the silicon elastomer rod, there is usually an existence of initial permanent strain which normally occurs when they are well-bonded to the rod during martensite phase . Under this condition, we also assume that there exists a stress-free initial combined structure. Thus, if the temperature is raised above austenitic finish for the combined structure, there will be an inducement of strain around the host medium, as a result of changes in length of the fiber connected with the martensite-austenite phase changes. Further to this, we assume that the fiber is in the martensite state, with an initial permanent axial strain represented by . This symbol also represents the strain associated with the stress-free shape recovery when the heating temperature rises above the austenitic finish. The strain would be recovered as long as the SMA fiber could undergo the martensitic-austenitic phase changes.
As a result of the constraint introduced by the elastomeric rod hosting the actuator, a self-imposed equilibrium state of internal stress will instead develop. And the pre-strain will be recovered partially.
Thus, the equations of equilibrium state of internal stress created by the cylindrical host medium, the shear layer and the SMA fiber are related as follows:
The equilibrium force in the axial direction of rod and the angular change, of the SMA fiber have the following relationship:
Where is the radius of the fiber; is the rod radius; is deduced from (refer to figures 5, 6 and 7).
Furthermore, represents the distance between SMA fiber axis and the central axis of the rod, as seen in figure 5.
is the shear-lag parameter and perhaps, a function of ; is the normal stress in the actuator in the axial direction and is the normal stress in the axial direction in the fiber; represents the shear stress in the shear layer which transfers the force from the actuating fiber to the cylindrical rod.
The actuating fiber and the cylindrical rod are related via a constitutive equation given as follows:
where is the Young’s modulus of the actuating fiber in the austenite phase; is the Young’s modulus of the elastomeric rod and is the shear modulus of the rod.
The fiber, the host medium and the shear layer have axial displacements represented as; and, respectively. Since and depend on , the axial displacements will also depend on .
The above equations have been solved by integration, using a general concept of an off-axis inserted SMA actuator. The concept is by first of all integrating equation (2), and getting its outcome embedded in (1), while (4) is applied, followed by the constitutive equations (5-7). Thus:
where is a typical length of the problem described by
and represents the stress in the actuating fiber away from the ends, summarized by:
By applying the boundary conditions, = 0, (L) = 0; then equation (8) is resolved into:
The shear-lag parameter in equation (11) remains the only undetermined value. But in an off-axis actuator inserted in a circular rod, is assumed to be equivalent to .
The interface shear stress between the actuating fiber and the rod is determined from (1) when (11) is substituted thereat, and then:
Because of the interface stress created, a distributed bending moment and a distributed axial force is applied on the cylindrical rod by the SMA fiber (see figure 8).
The distributed axial force which was deduced from the integration of the interface shear upon the boundary of the actuating fiber through the axis of the rod becomes:
The formula above, which is approximate for is obtained, if replaces in the solutions of and .
Figure 10 below shows a comparison between the approximate formula for and the exact.
Equation (13) above was numerically integrated via Mathematica (Wolfram 1991).
And because of the smallness of the difference between the exact and the approximate, it is justifiable to use the simpler approximate formula for as an estimate for force distribution on the structure. Figure 10 above, further explains the use of shear-lag parameter , where generates a graph from the following cases: and . The sensitivity of , which is small to , is explained by the solution which depends on logarithmically.
With reference to the coordinate system chosen in figure 1, the distributed bending moments will have the following components (see figures 5, 6, 7, and 8):
SMA fiber was inserted off-axially in a cylindrical rod cast with silicon elastomer and modeled as an appendage for a gripping robot. It was driven electrically to test its bending moment and gripping capability. Their deformations were measured as seen in figures 11 & 12, with readings taken there from.
The equilibrium equations of a rigid body that has distributed bending moments and applied distributed forces and a quadratic strain energy function are expressed as follows.
In the light of the above equations, and are the shear forces in the robot arm while is the axial force in the arm .
, and are the curvatures and twist of the deformed shape of the structure, respectively (see figure 11).
, and are the cross-sectional moments of inertia of the structure with respect to the principal axes, and the polar moment of inertia, respectively.
In the same vein, ,, represent the external forces distributed on the structure while ,, refer to couples on the structure.
Furthermore, G2 and E2 refer to the shear Modulus and Young’s Modulus of the structure respectively.
Considering an off-axis SMA actuator inserted in a straight elastomeric rod, with a cylindrical configuration, the coordinate system can be chosen such that , , while represents the arc length of the central axis connecting the the structure (see figure 5). Thus, equations (13) and (14ii) will beget and respectively, if replaces
The resultant equations of equilibrium of couple of forces on the structure, under this condition would resolve into the following:
Equations (21), (22) & (23) are non-linear Ordinary Differential
Equations. They are supported by the following kinematic equations:
where is the cross sectional rotation of the structure about axis.
An appendage to imitate the bending and gripping tendency of a soft robot such as octopus or elephant tusk was designed and modeled. Their deformations, bending moment and compressive force distribution mechanism were analyzed using an approximate analytical model to determine their comparative shear-lag parameter, with a view to determining their fitness in a mechanical environment where there is need for high yield, high force and high strain.
By considering a finite length of the structure designed, equations (21) - (26), which are Ordinary Differential Equations, were numerically integrated by applying a 4th order Runge-Kutta scheme.
All calculations were symbolically done via Mathematica .
The integration of equations (21) - (26) was done from S = 0 to S = L via an initial presumption for , and; but , and had an error which was used to fill-in , and by adding some increments to them gotten from the solution to a system of linear equations shown as below:
*This work was supported by the National Natural Science Foundation of China (No. 50975270) and the Fundamental Research Funds for the Central Universities of China (No. WK2090090002).
Michael. O. Obaji holds a PhD in Instrumentation Engineering from the University of Science & Technology of China. He was invited in 2014 by the Ebonyi State Government of Nigeria to install and test-run four production lines of HDPE Pipes Production Machines imported from China. He was later invited again to be the General Manager of the facility now called Ebonyi Pipes Production Limited. He is also with the Department of Mechatronics Engineering of the Federal University, Oye-Ekiti, Nigeria. His research interests include smart materials and their applications in bio-robots; intelligent mechanics; 3D-printing & Mechatronics systems.