3.1. Feedstock Preparation and Magnetorheological Properties
Figure. 3 shows the mixing torque of CIP/SE-130AB mixtures across a volumetric ratio of CIP loading, ranging from 4 to 43 vol% of the total volume of powder (CIP) and polymer (SE-130AB) mixtures The mixing torque gradually increased as the CIP volume fraction increased, with the critical solids loading point occurring at the maximum rate of change in mixing torque, measured at 32.6 vol%. The critical solids loading means the closest packing composition that powder and polymer are mixed without formation of void. Therefore, it serves as a threshold point for the injection molding process, beyond which viscosity begins to increase infinitely due to the severe particle frictions from excess powder fraction [
26]. Above the critical point, the flow index of the composition deteriorates significantly, rendering injection molding unfeasible. The critical point for the CIP/SE-130AB mixtures was measured to be approximately 32.6 vol%, and thus we set the optimal solids loading at 30 vol%, about 3 vol% lower than the critical point, based on the recommended optimal loading condition (2 to 5 vol% lower than the critical point) for powder/binder mixture injection molding [
24]. In summary, the feedstock with 30 vol% CIP addition was designed to contain maximum powder fraction, allowing the developed micropillar arrays to achieve maximum magnetic responsiveness during magnetic actuation.
The rheological behaviors of the feedstock were characterized by the effect of temperature, shear rate, and magnetic field. Firstly, the off-state (0 Tesla) shear viscosity was measured over a high shear rate range from 100 to 10,000 s-1 at temperatures 160, 180, 200 ˚C using a capillary rheometer, as shown in
Figure 4a. The recorded shear viscosity remained below 300 Pa·s, which is within the desirable viscosity range (below 1,000 Pa·s) for shear rates between 100 to 10,000 s-1 during molding [
24]. The off-state viscosity exhibited pseudo-plastic behavior, where the molten feedstock’s shear viscosity decreased as the shear rate increased. The neat polymer, SE-130AB, also demonstrated typical pseudoplastic behavior at high shear rates due to the unraveling of polymer molecules with increasing shear rate. This pseudoplastic behavior can be described using the Power Law model. Additionally, the temperature sensitivity of the shear viscosity is characterized by the Arrhenius model. The combined effects of strain rate and temperature on rheology can be represented by these models, as given below:
Where η is the apparent viscosity,
is the shear rate, n is the Power law exponent, T is the temperature, R is the gas constant, and E is the activation energy for Arrhenius model reflecting temperature sensitivity on the shear viscosity. A power law exponent of n < 1 indicates non-Newtonian fluid behavior with pseudoplasticity. The Power-law exponent, n, varied from 0.215, 0.229, and 0.257 for 160, 180, and 200 ˚C, respectively. Since the Power-law exponent, n, is also dependent on the temperature, it increased as the temperature increased with the linear progression relationship [
27]. The Arrhenius activation energy, E, was determined by calculating the slope of
versus
as shown in
Figure 4b. The Arrhenius activation energy, E, was calculated as 11,020 J/mol, which is significantly lower than that of other thermoplastic polymers, such as polypropylene (PP) and PDMS [
28,
29]. This result indicates that the CIP/SE-130AB feedstock is relatively insensitive to the molding temperature range of 160 to 200°C.
In
Figure 5, the dynamic magnetorheological response was analyzed using a magnetic frequency sweep test that gather the internal structure change of polymer or feedstock system within the linear viscoelastic (LVE) region under an external magnetic field. The storage modulus (G’) and loss modulus (G’‘) of the feedstock at 160°C are plotted in
Figure 5a, with the external magnetic field increasing from 0 to 0.7 Tesla. The results indicate that the feedstock exhibited a higher G’ than G’‘ under both magnetic and non-magnetic field conditions, demonstrating that the CIP/SE-130AB feedstock behaves like a solid, with its elastic portion being predominant over the viscous portion. As the magnetic field increased, both storage and loss moduli increased, as shown in
Figure 5a. The strength of the dipole-dipole interactions among the magnetic powders increased with the higher magnetization, so that the magnetorheological stiffening effect increased with increasing external field, and consequently forming a chain-like structure among the particles from a macroscale viewpoint. The CIPs are magnetically aligned along the magnetic field direction and subsequently agglomerated and assembled into magnetic columnar structures due to the dipole-dipole magnetic interactions among the adjacent neighboring particles [
30,
31]. Within the LVE region, the CIP/SE-130AB feedstock exhibited a stable or slightly increasing G’ as a function of angular frequency, which results from the extrinsically formed magnetic-column structure of particle-particle interactions as well as intrinsically retained sturdy structure of elastomeric matrix. Moreover, G’‘ above 0.3 T external field shows a decreasing trend, even with an increase in angular frequency, suggesting the accumulative stiffening effect of dipole interactions.
Figure 5b illustrates the results of complex viscosity (η*) with increasing magnetic field, comparing the on-state (0.1 to 0.7 T) and off-state (0 T) dynamic shear viscosities. Complex viscosity was calculated as the complex conjugate of the elastic and viscous components [
32]. η* progressively increased as the external magnetic field intensified, with the increment rate from 0 to 0.4 T being significantly higher than from 0.4 to 0.7 T. To better understand the dynamic magnetorheological response relative to field intensity, the normalized complex viscosity (
), was calculated for each angular frequency from 1.58 to 100 rad/s (
Figure 5c). The normalized increment plots displayed a sigmoidal S-shape, and the first derivative showing a bell-shaped increment rate (
Figure 5d). The rate of increase in η* peaked in the 0.3 to 0.4 T range, then gradually decreased after this inflection point (shaded area in
Figure 5d). This behavior in complex viscosity is attributed to particle alignment and the formation of columnar structures. Initially, in the magnetic field region from 0 to 0.4 T, the CIPs become magnetized and align along the field direction, forming chain-like structures. In this region, the CIPs rapidly magnetize due to the intrinsic low residual magnetization (0.75 emu/g). However, the chain-like structure still has a weak interaction among the particles since the external field is below the saturation magnetization field of around 0.6-0.7 T, as seen in the B-H hysteresis loop (
Figure 2b). Beyond the inflection point, stronger dipole interactions are developed as magnetization increases, leading to the formation of numerous columnar structures with more CIPs. Eventually, these interactions and the strength of the column structures reach a plateau as the field approaches 0.7 T. Furthermore, the slope of the normalized complex viscosity decreased with increasing angular frequency. At low frequencies, magnetized particles form stable and robust chain-like structures in response to a strong external magnetic field, resulting in the maximum slope of the normalized increment ratio of complex viscosity. In contrast, as angular frequency increased, hydrodynamic interactions hindered the alignment of CIPs and the formation of column structure, leading to a decrease in the slope of the increment rate in complex viscosity [
33]. In summary, the variation in complex viscosity with magnetic field is directly related to the dynamic magnetorheological response in the feedstock and can be used as an indicator for designing the magnetic alignment of micropillar arrays in MREs.
3.2. Development of LIGA Process and Soft Magnetic Microarrays
The LIGA process for developing the micro-patterned PMMA insert mold consisted of four steps: 1) UV photolithography, 2) X-ray mask design, 3) X-ray lithography, and 4) development. The target design of the micro pattern is listed in
Table 2, with square pattern widths of 50, 70, 150, and 350 μm. Through X-ray lithography using high phonon flux energy, we achieved a high width-to-height ratio for the micropillar pattern, up to 10, with the pattern width of 50 μm.
Figure 6a (i) to (ii) show the intermediate photomask prior to the UV lithography process, which was made using a silicon wafer, polyimide film, Au/Cr seed layer, and SU-8 3010 (MicroChem, USA). A 4-inch silicon wafer was prepared as the base substrate for the photomask [
25]. A 450 μm photosensitive polyimide film was adhered on top of the silicon wafer, acting as a transmission layer with high transparency. Next, the seed layer, consisting of Au and Cr, was deposited on the polyimide film. Prior to UV lithography, a negative type UV photoresist, SU-8 3010, was coated onto the seed layer. The SU-8 3010 solution was spin-coated at speeds ranging from 300 to 1,000 rpm, resulting in an approximately 10 to 16 μm thick negative photoresist layer.
In the first step, the SU-8 3010 layer was exposed to a UV light source through the working mask, as shown in
Figure 6a (iii). Since SU-8 3010 is a negative photoresist, the exposed structure underwent cross-linking and strengthening during lithography. The SU-8 3010 layer was then baked at 65 ˚C for 1 min, followed by an successive heating at 95 ˚C for 5 mins. The non-exposed portions were dissolved by SU-8 liquid developers, and intermediate UV photomasks are shown in
Figure 6b. In the second step, a gold layer was electroplated onto the SU-8 3010 layer of the photomask to serve as an X-ray absorber as shown in
Figure 6a (iv). Au plating was performed on the SU-8 3010 pattern to develop the X-ray absorber. Au is widely used as an X-ray absorber material due to its extremely low transmittance for X-ray irradiation [
34]. The gold layer was plated to a thickness of around 10 to 15 μm.
Figure 6c shows the final X-ray gold mask for X-ray LIGA process. Next, the positive photoresist PMMA was prepared for the sacrificial insert mold. PMMA sheets with thickness of 500 and 1,000 μm were attached to a graphite substrate. Subsequently, X-ray lithography was performed through the X-ray gold mask onto the positive photoresist PMMA mold. The target X-ray exposure energy was set to 4.5 kJ/cm
3. After X-ray exposure, the exposed PMMA was etched using a GG developer in the fourth step. The GG developer consists of 5 vol% 2-aminoethanol, 15 vol% deionized water, 20 vol% tetrahydro-1, 4-oxazine, and 60 vol% 2-(2-butoxyethoxy) ethanol. The sacrificial PMMA insert molds, shown in
Figure 7, exhibited no visible defects. The micropatterns of the PMMA molds had pattern width ranging from 50 to 350 μm, with sheet thickness of 500 and 1,000 μm, achieving aspect ratios of up to 10.
Table 2 provides detailed information on the micropattern insert molds, including square pattern width, height, and interpillar spacing. Using the developed sacrificial PMMA molds in
Figure 7, the mold insert injection molding process was carried out under an external magnetic field aligned along the micropillar direction. The field intensity within the mold cavity was set to approximately 0.9 T, a strength sufficient to align the particles along the micropillars, surpassing the peak gradient of 0.3 to 0.4 T as well as the pleatue field intensity of 0.7 T as shown in
Figure 5c and d. After the field-induced injection molding into the insert mold, the PMMA insert molds were removed in the solvent debinding step. The molded samples were immersed in 99.5 % acetone solution for 3 hours at room temperature, allowing them to chemically dissolve and depolymerize.
Figure 8 organizes 3-D optical microscope image of the high-aspect-ratio soft magnetic microarrays.
An external magnetic field of 0.5 to 0.6 T was applied to induce the magneto-mechanical responses (referred to as magnetic actuation) of the developed magnetic microarrays. As shown in
Figure 9, the microarrays exhibited bending motion in response to the linear external magnetic field. The particle chain structures along the axial direction of the arrays generated magnetic dipole moments in the presence of the magnetic gradient field, producing magnetic bending torques along the linear direction of the external field. When a gradient magnetic field is applied to the magnetic microarrays, the resulting magnetic body force per unit volume generated by the soft-magnetic particles can be expressed as:
Where∇B is the gradient of the external magnetic field, and M is the magnetic dipole moments of the magnetic suspensions, respectively. Soft magnetic materials, such as CIP, are highly responsive to external magnetic fields due to their low coercivity (H
c), which is the force required to demagnetize the material. Additionally, their narrow hysteresis loop allows them to be easily magnetized and demagnetized in alignment with the applied magnetic field [
35,
36]. In addition to single bending actuation, we confirmed that by switching the external magnetic field on and off, the soft magnetic microarrays exhibited reversible actuation motion, demonstrating the possibility of reversible shape reconfiguration (see supplementary video S1). Our developed soft magnetic microarrays with square pattern width of 50, 70, and 150 μm, and aspect ratios of 10.00, 7.14, and 6.67, exhibited bending angles of 79, 65, and 81˚, respectively. This is attributed to the highly responsive magneto-mechanical bending moments inherent in the high-aspect ratio microarray structures.