The FEM Model of the Pump Made of Dielectric Electroactive Polymer Membrane

Dielectric electroactive polymers (DEAPs) undergo large deformations when subject to an electric field, which make them an attractive material for use in novel actuator systems. This article presents the possibility of using DEAPs to model an innovative pumping actuator structure. The model was used to map important object parameters at individual operating points of the modeled pump. The experimental work involved designing the membrane and testing its changes in elasticity under the influence of varying forces and voltage supplies. In the further part of the work, a finite element model (FEM) of a pumping device was implemented. In the new construction of the pump, pressure is generated by membrane deformation. This is due to electrostatic compressive force between two electrodes applied to the polymer surface and forces generated by permanent magnets. The results are presented graphically, confirming the compliance of the model with the measurements.


Introduction
Smart materials are ones that can change their properties in a controlled manner in response to environmental stimuli. Such materials can be used both as sensors and actuators [1,2]. Electroactive materials like DEAPs (dielectric electroactive polymers) are included in the group of smart materials and are increasingly used in many automation and robotics systems [3]. This work presents the concept of a device using an electroactive polymer that changes its mechanical properties under the influence of voltage excitation.
Finite element model (FEM) modeling of the behavior of DEAPs is a useful tool to understand such systems better and help in the optimal design of prototypes. The modelling and simulating of DEAP actuators is a cost-effective way of providing a better understanding of the devices built with their usage. DEAPs offer excellent performance, are lightweight, flexible and inexpensive. Therefore, dielectric elastomers, provide many potential applications as micro-actuators. A technique to accurately model this material-taking into account its nonlinearities as well as its large deformations-is being developed in this study. This work presents a model of the innovative construction of the pump using DEAP membrane deformation.
There are many works that describe different modeling methods and analytical models [4]. One of the works presents a technique for accurate modeling of electroactive polymer taking into account nonlinearities and large deformations [5]. The authors simulated an algorithm in the ANSYS environment that considers the electromechanical coupling of the EAP (electroactive polymers) model. Simulation of charge distribution and electrostatic force is the subject of another work which uses COMSOL software [6]. Another work focuses on studying the kinematics of forming and loading, along with thermodynamics [7]. While much theoretical and computational modeling effort has gone into describing the ideal, time-independent behavior of these materials, viscoelasticity is a crucial component of the observed mechanical response, and hence has a significant effect on electromechanical actuation. In [8], the authors show that viscoelasticity provides stabilization that delays the onset of instability under monotonic loading and may fully suppress instabilities under sufficiently fast cyclic loading, which may be desirable for many applications. Although all of these models provide advantageous insights into the behavior of the elastomer material, they lack the ability to accurately predict the complete DEAP actuator structure. They focus on several aspects of modeling but are not able to reflect the full complexity.
In this work, an FEM model was made that was verified faithfully with measurements. First, the relationship between static deformation and applied force was presented. Next, three series of displacement measurements were compared for a given force during alternating voltage applied to the membrane. The results prepared in this way allowed us to identify the values of the virtual Young's modulus corresponding to individual points of states by collecting a sufficient number of real object measurements. Using the work state mapping method, it was possible to model devices with sufficient precision having the knowledge how the membrane will deform; it was not necessary to calculate the full complexity of the model each time.

Principle of Operation of DEAP Actuator
A DEAP membrane is made of a silicone membrane that is pre-stretched during the production process. Then, the carbon grease electrodes are printed. In principle, DEAPs are similar to flexible capacitors. As can be seen in the Figure 1, they consist of an elastomeric film sandwiched between two electrodes. Applying voltage to them leads to compression of the elastomer, which causes a change in capacity [9][10][11][12][13].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 2 of 13 crucial component of the observed mechanical response, and hence has a significant effect on electromechanical actuation. In [8], the authors show that viscoelasticity provides stabilization that delays the onset of instability under monotonic loading and may fully suppress instabilities under sufficiently fast cyclic loading, which may be desirable for many applications. Although all of these models provide advantageous insights into the behavior of the elastomer material, they lack the ability to accurately predict the complete DEAP actuator structure. They focus on several aspects of modeling but are not able to reflect the full complexity. In this work, an FEM model was made that was verified faithfully with measurements. First, the relationship between static deformation and applied force was presented. Next, three series of displacement measurements were compared for a given force during alternating voltage applied to the membrane. The results prepared in this way allowed us to identify the values of the virtual Young's modulus corresponding to individual points of states by collecting a sufficient number of real object measurements. Using the work state mapping method, it was possible to model devices with sufficient precision having the knowledge how the membrane will deform; it was not necessary to calculate the full complexity of the model each time.

Principle of Operation of DEAP Actuator
A DEAP membrane is made of a silicone membrane that is pre-stretched during the production process. Then, the carbon grease electrodes are printed. In principle, DEAPs are similar to flexible capacitors. As can be seen in the Figure 1, they consist of an elastomeric film sandwiched between two electrodes. Applying voltage to them leads to compression of the elastomer, which causes a change in capacity [9][10][11][12][13]. For this circular geometry, the deflection of the DEA membrane is described by the radial stretch λr, circumferential stretch λc and thickness stretch λz [14].
I propose a dynamic model composed of a set of nonlinear, time-invariant differential equations describing the dynamic relationship between the input voltage u and the output actuator displacement d. I present a model which accurately describes the behavior of the actuator, both in transient and steady state, and for different mass loads. The vertical force equilibrium causing membrane deformation can be defined as: For this circular geometry, the deflection of the DEA membrane is described by the radial stretch λ r , circumferential stretch λ c and thickness stretch λ z [14].
I propose a dynamic model composed of a set of nonlinear, time-invariant differential equations describing the dynamic relationship between the input voltage u and the output actuator displacement d. I present a model which accurately describes the behavior of the actuator, both in transient and steady state, and for different mass loads. The vertical force equilibrium causing membrane deformation can be defined as: m where d is a distance, θ is an angle, m is the mass of the initial load, g means standard gravity and F L defines the external load. The electromechanical coupling is produced by Maxwell Stress which compresses the membrane when a voltage u is applied to the DEAP. The actuation is in the thickness direction and due to constant volume produces the change of radial stretch y r . The value of F M is equal to: where ε 0 is vacuum permittivity and ε r is the relative permittivity of the polymer material. The parameters z 0 , l 0 and r describe the geometry of DEAP actuator and are defined in Table 1.
Force F u considers dynamic viscoelastic process according to: . There are many hyperelastic strain energy density functions, such as Yeoh, Ogden, Arruda-Boyce, etc. Similar to [11], I use the Ogden model and define two parameters β i and γ i to obtain to the expression for the hyperelastic stress σ e : where α i is 2, 4 and 6. The main difficulty encountered with the identification of the DEAP parameters in the dynamic case is to obtain the damping coefficients b and the coefficients of the viscoelastic k e and η e . In this case the optimization procedure is based on both a step response and a frequency response (the results of this work were published in the [15,16]). Viscous friction is defined as: where b is damping coefficient.

Static Characteristics
In order to develop a faithful mathematical model, it was necessary to obtain real device parameters by measuring and identification process. Figure 2 presents the laboratory set used to collect displacement measurements depending on the static force acting on the DEAP membrane. The laboratory set contained a high voltage amplifier TREK MODEL 10/10B-HS, laser distance sensor Micro-Epsilon optoNCDT ILD1320-10 with 1 µm accuracy and an Inteco RT-DAC/USB data acquisition card. The actuator presented in our research was made of 3M VHB tape which was pre-stretched during the production process (1 mm to 200 µm thickness). The dimensions of the actuator created by the author are presented in detail in Table 2.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 4 of 13 In order to develop a faithful mathematical model, it was necessary to obtain real device parameters by measuring and identification process. Figure 2 presents the laboratory set used to collect displacement measurements depending on the static force acting on the DEAP membrane. The laboratory set contained a high voltage amplifier TREK MODEL 10/10B-HS, laser distance sensor Micro-Epsilon optoNCDT ILD1320-10 with 1 μm accuracy and an Inteco RT-DAC/USB data acquisition card. The actuator presented in our research was made of 3M VHB tape which was prestretched during the production process (1 mm to 200 μm thickness). The dimensions of the actuator created by the author are presented in detail in Table 2.     The obtained data were used to identify the value of Young's modulus in the developed FEM model. Figure 3 shows the membrane's FEM model loaded with force 1.42 [N]. Table 3 and Figure 4 provide a faithful comparison of simulations with measurements in the full load range of the actuator.

The identification process
The operation of the device in an actuator mode using the DEAP membrane requires a power supply. To collect reference data, a number of experiments on the deformation of the actuator

The Identification Process
The operation of the device in an actuator mode using the DEAP membrane requires a power supply. To collect reference data, a number of experiments on the deformation of the actuator membrane were carried out under the influence of various supply voltages. The relatively narrow range of membrane supply was sampled at nine points from 5.0 [kV] to 7.0 [kV] with a 0.25 [kV] step. This approach allowed the determination of the virtual value of Young's modulus compensating for the nonlinear structure of the DEAP model.
To verify the correctness of the actuator modelling methodology, measurements were made for three different forces acting on the actuator. Tables 4-6

Voltage [kV]
measurement "model"      The results of these experiments show that the parameters obtained from the voltage test sufficiently accurately describe the behavior of the membrane and can be used in the model.

Pump FEM model with DEAP membrane
Based on the obtained values of the virtual Young's module, a pump model was developed. In order to increase the force and thus the scope of the pump operation, two permanent magnets were implemented (see Figure 8)  The results of these experiments show that the parameters obtained from the voltage test sufficiently accurately describe the behavior of the membrane and can be used in the model.

Pump FEM Model with DEAP Membrane
Based on the obtained values of the virtual Young's module, a pump model was developed. In order to increase the force and thus the scope of the pump operation, two permanent magnets were implemented (see Figure 8) [T]. The force between two identical cylindrical magnets was approximately calculated according to Gilbert's model [17,18]. This force in the start position without the pump supply was 0.8 [N], which caused a static deformation of membranes equal to 8 [mm].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 13 A pump chamber was proposed (see Figure 9), which is a segment of a sphere with a radius of 9 [cm]. A pump chamber was proposed (see Figure 9), which is a segment of a sphere with a radius of 9 [cm].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 13 Figure 8. A pump concept built using a DEAP membrane and permanent magnets.
A pump chamber was proposed (see Figure 9), which is a segment of a sphere with a radius of 9 [cm].  The boundary conditions applied in ANSYS software closely match the boundary conditions during the experiments. In the experiments the frame of the DEAP was clamped rigidly and the force was applied at the central part of the actuator. The actuator was designed to move out of plane by adding a fixed boundary constraint to the perimeter of the material. The scheme of the applied forces and the mechanical boundary conditions applied in ANSYS are as shown in the Figure 10. Since, the outer plastic frame was rigid as compared to the elastomer, the outer end of the elastomer could be considered fixed to a rigid frame. Similarly, the circular part of the frame at the center of the EAP was rigid and hence no radial displacement of the polymer occurs at the inner end. The force was applied at the inner boundary in the vertically downward direction. For electromechanical simulations, the area surrounding the EAP was specified as air. The directional distribution of forces attracting magnets is shown in the Figure 11. The boundary conditions applied in ANSYS software closely match the boundary conditions during the experiments. In the experiments the frame of the DEAP was clamped rigidly and the force was applied at the central part of the actuator. The actuator was designed to move out of plane by adding a fixed boundary constraint to the perimeter of the material. The scheme of the applied forces and the mechanical boundary conditions applied in ANSYS are as shown in the Figure 10. Since, the outer plastic frame was rigid as compared to the elastomer, the outer end of the elastomer could be considered fixed to a rigid frame. Similarly, the circular part of the frame at the center of the EAP was rigid and hence no radial displacement of the polymer occurs at the inner end. The force was applied at the inner boundary in the vertically downward direction. For electromechanical simulations, the area surrounding the EAP was specified as air. The directional distribution of forces attracting magnets is shown in the Figure 11. Since, the outer plastic frame was rigid as compared to the elastomer, the outer end of the elastomer could be considered fixed to a rigid frame. Similarly, the circular part of the frame at the center of the EAP was rigid and hence no radial displacement of the polymer occurs at the inner end. The force was applied at the inner boundary in the vertically downward direction. For electromechanical simulations, the area surrounding the EAP was specified as air. The directional distribution of forces attracting magnets is shown in the Figure 11. In the suction cycle, the inlet slot (valve) opens and the outlet slot (valve) closes. In order to pump out the medium, the valves change their state and a voltage of 7 [kV] is applied to the membranes. The change in volume of one pump chamber ranges from 135.64 cm 3 to 379.81 cm 3 and its graph over time is shown in Figure 12. In the suction cycle, the inlet valve opens, and the outlet valve closes. In order to pump out the liquid, the valves change their state and the voltage on the membrane electrodes is switched on. In the suction cycle, the inlet slot (valve) opens and the outlet slot (valve) closes. In order to pump out the medium, the valves change their state and a voltage of 7 [kV] is applied to the membranes. The change in volume of one pump chamber ranges from 135.64 cm 3 to 379.81 cm 3 and its graph over time is shown in Figure 12. In the suction cycle, the inlet valve opens, and the outlet valve closes. In order to pump out the liquid, the valves change their state and the voltage on the membrane electrodes is switched on.

Conclusions
DEAPs are a subclass of electroactive polymers in which actuation is produced by an elastic deformation resulting from the compressive electrostatic forces. Their most important advantages like energy efficiency, lightweight and scalability are key features for actuators in applications such as pumps. Dielectric electroactive polymer technology is able to fulfill requirements as well as commonly used technology e.g., solenoids, but its limitations concern relatively low force [19,20]. One way to increase force generated by these devices can be the application of neodymium permanent magnets. The concept for pump presented in this paper consists of a stack of DEAP membranes combined with a permanent magnets. These two components are combined in a novel pump construction, which allows a compact design by integrating the biasing mechanism with the DEAP membranes. Subsequently, the single components are manufactured, tested, and their forcedisplacement characteristic is documented. Identifying the parameters experimentally allowed to develop a faithful FEM model of the pump device and make its simulations for different voltage. FEM modeling of the DEAP membrane behavior is a useful tool to understand such systems better and helps to create the optimal design of prototypes. These modeling efforts must account for the electromechanical coupling in order to accurately predict their response to multiple loading conditions expected during real operating conditions. Because of the complexity of the nonlinear processes, the finite element analysis is reasonable to be used and study the behaviour of dielectric elastomer pump. The experiments under static and transient considerations have been performed to validate the design and implementation. It is shown, that the model is in close agreement with the measured responses of the real DEAP actuator. The collected data allow for the analysis of various other device designs using DEAP membrane in the future.

Conclusions
DEAPs are a subclass of electroactive polymers in which actuation is produced by an elastic deformation resulting from the compressive electrostatic forces. Their most important advantages like energy efficiency, lightweight and scalability are key features for actuators in applications such as pumps. Dielectric electroactive polymer technology is able to fulfill requirements as well as commonly used technology e.g., solenoids, but its limitations concern relatively low force [19,20]. One way to increase force generated by these devices can be the application of neodymium permanent magnets. The concept for pump presented in this paper consists of a stack of DEAP membranes combined with a permanent magnets. These two components are combined in a novel pump construction, which allows a compact design by integrating the biasing mechanism with the DEAP membranes. Subsequently, the single components are manufactured, tested, and their force-displacement characteristic is documented. Identifying the parameters experimentally allowed to develop a faithful FEM model of the pump device and make its simulations for different voltage. FEM modeling of the DEAP membrane behavior is a useful tool to understand such systems better and helps to create the optimal design of prototypes. These modeling efforts must account for the electromechanical coupling in order to accurately predict their response to multiple loading conditions expected during real operating conditions. Because of the complexity of the nonlinear processes, the finite element analysis is reasonable to be used and study the behaviour of dielectric elastomer pump. The experiments under static and transient considerations have been performed to validate the design and implementation. It is shown, that the model is in close agreement with the measured responses of the real DEAP actuator. The collected data allow for the analysis of various other device designs using DEAP membrane in the future.

Conflicts of Interest:
The author declare no conflict of interest.