Forming of Dynamic Microstructure of Flexible Polymer

This study focuses on the forming of dynamic microstructure of flexible polymer. The dynamic tensile control microstructure of the polymer mold, along with the gasbag, was used to exert pressure to achieve forming. This study simulated the dynamic control of the flexible mold, and proposed four mechanical models of material viscoelastic response for modeling and evaluation. MATLAB software was used to calculate the imprint prediction calculation theory construction according to the imprint result of curved surface and asymmetric imprint forming. This study designed and developed a gasbag-assisted dynamic forming system, and tested the proposed system for verification. The test results showed that the mechanical stability, curved surface, and asymmetric imprint prediction calculation of the mechanical model of the viscoelastic response of flexible mold material, as proposed in this study, can display the geometric features of the imprinted microstructure. The dynamic mold microstructure control process can accurately transfer a bifacial microstructure and construct the confidence interval for transfer printing forming.


Introduction
With the rapid advances in science and technology, continuous mass production is the future of industrial development, and system engineering follows the trend of microminiaturization. Under the establishment and development of microsystem technology, microsystems have been rapidly, effectively, and continuously used in science and technology industries. The microelectromechanical system (MEMS) process technology [1][2][3][4] and semiconductor process technology [5][6][7] are critical in the mass production of microsystem components. They are both integrated with multiple technologies, and generally arranged according to multiple sequential processes, including the cleaning and surface preparation of wafer preprocesses, photoresist coating and baking, alignment exposure, development, etching (or ion implantation), etc., as the final process of creating micrographic components or micro molds. The technologies are continuously updated and developed. The modeling process of microsystem components is required before bulk production. In the process, the mold material is usually used for fabricating accurate sized molds, where the MEMS technology is often used, while Mask, Photoresist, and Lithography technologies are used to fabricate microstructure components. As the MEMS technology can achieve forming of stable microstructure components, it has been widely used in the front end of line and middle end of line in the bulk production of most semiconductors. Although the MEMS technology has high accuracy, the equipment cost is high. Moreover, the silicon substrate is the main material, but the hard-brittle (fragile) silicon substrate may increase the defects in the microsystem structure components during bulk production, thus, the fabrication process requires experts in the field. In addition, a precision mold is also time-consuming to fabricate. The model molding technology has developed rapidly in recent years, and the common microstructure component model molding technologies at present include plastics micro-injection molding [8,9], micro-thermoforming [10][11][12], After elastic strain, as the strain increases, the velocity gradually reduces, reaching the final strain (ε ) state, as expressed by Equation (4) = + (4) Based on the above equations, the creep deformation mode can be expressed by Equation (5) for prediction computation, as shown in Figure 1b.
In addition, as the tensile is maintained in the dynamic tensile process, the polymer may have stress relaxation in this process, which is the key point of evaluation. This study sets the tensile deformation as a constant ( 0 ε ), the Maxwell model is used as the evaluation method. The threeparameter model is used to discuss and spread the preliminary equations, as expressed by Equation (6) After further compilation: However, due to σ = , after a long time, the final stress relaxation value is expressed as Equation (7) σ (7) In this study, the dynamic tensile is to control microstructure deformation, and this action is repeated and compared with stress relaxation, thus, the polymer creep phenomenon is more important.
The significance of the polymer is defined, as follows: where, S: Total tensile stress; : Tensile stress on elastic modulus; : Tensile stress on viscous modulus; : Total tensile strain; : Tensile strain on elastic modulus; : Tensile strain on viscous modulus (a) (b) Therefore, after the Maxwell model is compiled (Figure 2a), it is expressed as Equation (8) In a similar manner, the Kelvin-Voight model can be deduced (Figure 2b), and expressed as Equation (9) Given the fixed tensile, the estimated tensile stress (σ) is fixed, and then dε /dt + ( E /η)ε = 1 /η E 0 + E E 0 σ 0 Let, (k = d /dt, j = E /η, m, C be constants) After mathematical analysis, as expressed by Equation (2) ε = m /j+Ce −jt = ( E 0 +E /E 0 E)σ 0 +Ce −jt (2) According to the elastic strain in the initial stress range, ε (t=0) = ε 0 = σ 0/E 0 and then After elastic strain, as the strain increases, the velocity gradually reduces, reaching the final strain (ε f inal ) state, as expressed by Equation (4) Based on the above equations, the creep deformation mode can be expressed by Equation (5) for prediction computation, as shown in Figure 1b.
In addition, as the tensile is maintained in the dynamic tensile process, the polymer may have stress relaxation in this process, which is the key point of evaluation. This study sets the tensile deformation as a constant (ε 0 ), the Maxwell model is used as the evaluation method. The three-parameter model is used to discuss and spread the preliminary equations, as expressed by Equation (6) After further compilation: However, due to σ 0 = E 0 ε 0 , after a long time, the final stress relaxation value is expressed as Equation (7) σ f inal= Eσ0 /E0 + E In this study, the dynamic tensile is to control microstructure deformation, and this action is repeated and compared with stress relaxation, thus, the polymer creep phenomenon is more important.
The significance of the polymer is defined, as follows: where, S: Total tensile stress; S e : Tensile stress on elastic modulus; S a : Tensile stress on viscous modulus; ε: Total tensile strain; ε e : Tensile strain on elastic modulus; ε a : Tensile strain on viscous modulus Therefore, after the Maxwell model is compiled (Figure 2a), it is expressed as Equation (8)

Constructing Model A and Model B
The second series Model A and the complex in-tandem models are built, as shown in Figure 3a. The complex in-tandem model, Model B, is built, as shown in Figure 3b. The mechanical properties of mold dynamic tensile are established and evaluated, in order to determine the optimum model.
: can be replaced by ; : can be replaced by ; : can be replaced by ; : can be replaced by Reduced to Equation (13): Let, E , E , E be replaced by E; η , η , η can be replaced by η, then  In a similar manner, the Kelvin-Voight model can be deduced (Figure 2b), and expressed as Equation (9)

Constructing Model A and Model B
The second series Model A and the complex in-tandem models are built, as shown in Figure 3a. The complex in-tandem model, Model B, is built, as shown in Figure 3b. The mechanical properties of mold dynamic tensile are established and evaluated, in order to determine the optimum model.

Asymmetric Imprint Molding Prediction Geometric Shape Construction
This study used the array rectangular microstructure to establish the computational theory of the asymmetric imprint molding prediction method. The square column arrays microstructure is geometrically calculated. MATLAB software is used to establish the asymmetric imprint prediction calculation theory. Considering the square column arrays, the appearance of the nanoimprint patterns of each square column exhibits a rhomb shape. Referring to Weng's study [40] of the asymmetric theorem, this study set a single rectangular microstructure (Figure 4), array rectangular microstructure, array cylindrical microstructure, and the array microlens hemisphere structure for asymmetric imprint calculation prediction using MATLAB. Model A: η 1 : the first viscous modulus; η 2 : the sec ond viscous modulus; η3 : the third viscous modulus E 1 : the first elastic modulus; E 2 : the second elastic modulus; E 3 : the third elastic modulus; S: Tensile stress where, S e : Tensile stress on the first elastic modulus; S 2 : Tensile stress on lower Kelvin-Voight S 3 : Tensile stress on upper Kelvin-Voight; S a : Tensile stress on the first viscous modulus; where, ε: Total tensile strain; ε e : Tensile strain on the first elastic modulus; ε 2 : Tensile strain on lower Kelvin-Voight; ε 3 : Tensile strain on upper Kelvin-Voight; ε a : Tensile strain on the first viscous modulus; ε e : can be replaced by S e E 1 ; dε 2 dt : can be replaced by dt : can be replaced by ; dε a dt : can be replaced by S a Reduced to Equation (13): Let, E 1 , E 2 , E 3 be replaced by E; η 1 , η 2 , η 3 can be replaced by η, then Let, ε 2 + ε 3 be replaced by ε 2 , then, the Model A mold creep equation can be obtained, and expressed as Equation (14) dε Let, ε e be replaced by S e where, ε 2 : Tensile strain on the left of parallel joint; ε 3 : Tensile strain on the right of parallel joint; S 2 : Tensile stress on the left of parallel joint; S 3 : Tensile stress on the right of parallel joint where, S 2e : Tensile stress on the second elastic modulus on the left of parallel joint; S 2a : Tensile stress on the second viscous modulus on the left of parallel joint, and can be replaced by η 2 where, S 3e : Tensile stress on the third elastic modulus on the right of parallel joint; S 3a : Tensile stress on the third viscous modulus on the right of parallel joint, and can be replaced by η 3 Let, ε 2e be replaced by Let, dε a dt be replaced by S a η 1 ; Let, E 1 be replaced by E; η 1 , η 2 , η 3 be replaced by η Then, the Model B mold creep equation can be obtained, and expressed as Equation (27) dε

Asymmetric Imprint Molding Prediction Geometric Shape Construction
This study used the array rectangular microstructure to establish the computational theory of the asymmetric imprint molding prediction method. The square column arrays microstructure is geometrically calculated. MATLAB software is used to establish the asymmetric imprint prediction calculation theory. Considering the square column arrays, the appearance of the nanoimprint patterns of each square column exhibits a rhomb shape. Referring to Weng's study [40] of the asymmetric theorem, this study set a single rectangular microstructure (Figure 4), array rectangular microstructure, array cylindrical microstructure, and the array microlens hemisphere structure for asymmetric imprint calculation prediction using MATLAB.

Design and Development of the Gasbag-Assisted Dynamic Forming System with Microstructure Mold Preparation
This study developed an innovative gasbag-assisted dynamic forming system, which can perform single-sided and simultaneous bifacial gasbag-assisted dynamic forming. The forming resists are exposed and cured by a side array UV-LED lamp, as shown in Figure 5. The imprinting system developed in this study is comprised of the following main components: upper and lower filling gasbag cavities (filling pressure gas), a changeable transparent imprinting platform in the middle, eight array UV exposure systems on the side, and 8 angle adjustable high accuracy microstructure mold dynamic fine regulators in the middle. The proposed system can perform single-sided operations or unilateral exposure, and is flexible in the forming process. The mold microstructure is made by the laser processing of aluminum sheets and casting, and the cross microstructure mold material is PDMS (Polydimethylsiloxane, Sylgard™ 184, Dow Corning).

Gasbag-Assisted Compressive Stress Distribution Simulation and Microstructure Change Analysis
This study used Abaqus and ANSYS simulation software to analyze the gasbag-assisted pressure on cross mold stress and strain uniformity, in order to analyze the uniformity of gasbag filling pressure. The uniformity of stress on the cross microstructure mold in the tensile and forming processes was also analyzed to compare the errors in the actual experiment. In addition, the MATLAB software was used to predict the curvature change of the microstructure under pressure in the curved imprinting process for the preliminary identification of the curved imprint. The MATLAB software also calculated the microstructure replication forming shape under the stress of the asymmetric imprinting platform. The results were compared with the actual forming to verify the feasibility and effectiveness of the proposed system and prediction method.

Gasbag-Assisted Dynamic Microstructure Forming Uniformity Test and Molding Steps
The proposed gasbag-assisted dynamic forming system regulates unequal gas gage pressure. The imprint uniformity was tested by pressure sensitive adhesive tape. The experimental results showed the uniformly distributed pressure color blocks. The system forming steps include: (a) the cross microstructure mold is set up and put in the imprinting plate; (b) the location of the cross microstructure mold is regulated dynamically and the gasbag is filled to generate the force of impression, which is uniformly applied to the mold; (c) the relative imprinted depth of the imprinting plate and cross microstructure mold is controlled and the LED lamp is turned on for exposure forming; (d) the LED lamp is turned off and the gasbag is decompressed to obtain the finished microstructure component after transfer printing.

Gasbag-Assisted Compressive Stress Distribution Simulation and Microstructure Change Analysis
This study used Abaqus and ANSYS simulation software to analyze the gasbag-assisted pressure on cross mold stress and strain uniformity, in order to analyze the uniformity of gasbag filling pressure. The uniformity of stress on the cross microstructure mold in the tensile and forming processes was also analyzed to compare the errors in the actual experiment. In addition, the MATLAB software was used to predict the curvature change of the microstructure under pressure in the curved imprinting process for the preliminary identification of the curved imprint. The MATLAB software also calculated the microstructure replication forming shape under the stress of the asymmetric imprinting platform. The results were compared with the actual forming to verify the feasibility and effectiveness of the proposed system and prediction method.

Gasbag-Assisted Dynamic Microstructure Forming Uniformity Test and Molding Steps
The proposed gasbag-assisted dynamic forming system regulates unequal gas gage pressure. The imprint uniformity was tested by pressure sensitive adhesive tape. The experimental results showed the uniformly distributed pressure color blocks. The system forming steps include: (a) the cross microstructure mold is set up and put in the imprinting plate; (b) the location of the cross microstructure mold is regulated dynamically and the gasbag is filled to generate the force of impression, which is uniformly applied to the mold; (c) the relative imprinted depth of the imprinting plate and cross microstructure mold is controlled and the LED lamp is turned on for exposure forming; (d) the LED lamp is turned off and the gasbag is decompressed to obtain the finished microstructure component after transfer printing.

Cross Mold Dynamic Control Simulation Analysis
This study used Abaqus software for the simulation analysis of the related factors of cross mold in the actual process. After the mechanical properties of material were tested, various material parameters of the system were established as simulation parameters, as shown in Table 1. The simulation results showed that after the gasbag is uniformly filled with inert gas, the stress surface contacts the cross mold, and the microstructure surface and cross section have uniform displacement and strain, as shown in Figure 6. By von Mises stress and Tresca stress, the cross mold is in contact with the stress surface, microstructure surface, and cross section, displaying uniform and symmetric stress distribution, as shown in Figure 7. After dynamic tensile control, before and after the cross mold receives the gasbag pressure, the stress surface, microstructure surface, and cross section have highly uniform force distribution, as shown in Figure 8. As indicated above, the cross sections in all directions of the cross mold have highly uniform and symmetrical stress and strain. with the stress surface, microstructure surface, and cross section, displaying uniform and symmetric stress distribution, as shown in Figure 7. After dynamic tensile control, before and after the cross mold receives the gasbag pressure, the stress surface, microstructure surface, and cross section have highly uniform force distribution, as shown in Figure 8. As indicated above, the cross sections in all directions of the cross mold have highly uniform and symmetrical stress and strain.

Construction of Mold Mechanical Model
This study used different mixing ratios of Agent A and Agent B (5:1, 10:1, 15:1) for the mold material (PDMS) to discuss the creep phenomenon after dynamic regulation of the four mechanical models. The mechanical properties test data of PDMS material under unequal test conditions ( Table  2) were substituted in the four mechanical models to obtain the η value in the mechanical model ( Table 3). The results were compared with the η value of the four mechanical model creep equations to obtain the time-dependent creep behavior of dependent variables in the dynamic tensile process, as shown in Figures 9-12. After the four models were built according to the different mixing ratios of

Construction of Mold Mechanical Model
This study used different mixing ratios of Agent A and Agent B (5:1, 10:1, 15:1) for the mold material (PDMS) to discuss the creep phenomenon after dynamic regulation of the four mechanical models. The mechanical properties test data of PDMS material under unequal test conditions (Table 2) were substituted in the four mechanical models to obtain the η value in the mechanical model ( Table 3). The results were compared with the η value of the four mechanical model creep equations to obtain the time-dependent creep behavior of dependent variables in the dynamic tensile process, as shown in Figures 9-12. After the four models were built according to the different mixing ratios of Agent A and Agent B, and according to the actual test, the Kelvin-Voight model cannot effectively correspond when the mixing ratio is 15:1 (less able to display the mold materials mechanical properties, as shown in Figure 9), while the precisions of the other three mold mechanical models of sequence are Model A, Model B, and the Maxwell model, that the material properties with respect to creep are more stable, as shown in Figures 10-12. This study used Model A as the model to consider creep in the experimental process.

Effect of Curved Imprinting on Mold Microstructure
This study designed the replaceable imprinting platform as the curved surface (convex, concave) to discuss the changes in the mold microstructure in the imprinting process of the convex surface of the imprinting platform, where the material parameters are given (modulus of elasticity: 1.70 Mpa, Poisson's ratio: 0.495, density: 965 kg/m 3 ). ANSYS software was used to perform the simulation,

Effect of Curved Imprinting on Mold Microstructure
This study designed the replaceable imprinting platform as the curved surface (convex, concave) to discuss the changes in the mold microstructure in the imprinting process of the convex surface of the imprinting platform, where the material parameters are given (modulus of elasticity: 1.70 Mpa, Poisson's ratio: 0.495, density: 965 kg/m 3 ). ANSYS software was used to perform the simulation, under imprinting pressure of 0.12 Mpa and surface fillet of R42mm without tensile. According to von Mises stress simulation distribution, the central region of the surface was subject to higher stress, as shown in Figure 13.

Effect of Curved Imprinting Under Unequal Curvature Radius on Mold Microstructure
When the imprinting platform was designed as a concave surface, MATLAB software was used to analyze the changes in the cross mold microstructure under unequal concave curvature radius surface (25 mm, 30 mm). This study used a columnar microstructure for analysis. According to the analysis result, with a small curvature radius (25 mm), the longitudinal deformation and strain of the microstructure of the curved surface are larger than those of the curved surface with a large curvature radius (30 mm); the lower the microstructure, the lower the central area of the mold, the larger the longitudinal deformation and strain, and the larger the longitudinal deformation and strain on the mold edge, as shown in Figures 14 and 15, respectively.

Effect of Curved Imprinting Under Unequal Curvature Radius on Mold Microstructure
When the imprinting platform was designed as a concave surface, MATLAB software was used to analyze the changes in the cross mold microstructure under unequal concave curvature radius surface (25 mm, 30 mm). This study used a columnar microstructure for analysis. According to the analysis result, with a small curvature radius (25 mm), the longitudinal deformation and strain of the microstructure of the curved surface are larger than those of the curved surface with a large curvature radius (30 mm); the lower the microstructure, the lower the central area of the mold, the larger the longitudinal deformation and strain, and the larger the longitudinal deformation and strain on the mold edge, as shown in Figures 14 and 15, respectively.

Curved Imprinting Dynamic Regulation Microstructure Replication Forming Analysis
This study built an array cylinder with a diameter of 250 um and a height of 75 um in a cross microstructure mold. In the actual curved imprinting process, the uniaxial extension of the cross mold microstructure was adjusted (0, 2, 4, 6, 8, 10 mm), pressurized to 0.15 Mpa, and the pressure was maintained. After UV exposure, it was measured by a white light microscope, and the height of the

Curved Imprinting Dynamic Regulation Microstructure Replication Forming Analysis
This study built an array cylinder with a diameter of 250 µm and a height of 75 µm in a cross microstructure mold. In the actual curved imprinting process, the uniaxial extension of the cross mold microstructure was adjusted (0, 2, 4, 6, 8, 10 mm), pressurized to 0.15 Mpa, and the pressure was maintained. After UV exposure, it was measured by a white light microscope, and the height of the microstructure adjacent to the lowest columnar microstructure was close to the simulation analysis, as shown in Table 4, Figures 16 and 17. microstructure adjacent to the lowest columnar microstructure was close to the simulation analysis, as shown in Table 4, Figures 16 and 17.

Asymmetric Imprint MATLAB Prediction Analysis and Experiment Accuracy
This study built a rectangular microstructure for prediction analysis, and the imprint forecast result of the unequal angle was calculated. The predicted position variation was calculated according to a single rectangular microstructure (1 × 1) at rotation angle π/9° and tilt angle π/6°, as shown in Figure 18. The calculation prediction of the array rectangular microstructure (4 × 4), where the results of array (4 × 4) at rotation angle 0°, tilt angle π/6° and rotation angle π/ 6°, tilt angle π/6°, is shown in Figure 19. The dynamic tensile asymmetric imprinting experiment was performed. The array microlens hemispherical structure, as derived from the actual dynamic tensile, was reproduced by rotation angle π/6° and tilt angle π/36°. This result is close to the calculation result of MATLAB software, which is enough to show the imprinting stability of this system and the feasibility of the prediction method, as shown in Figure 20.

Asymmetric Imprint MATLAB Prediction Analysis and Experiment Accuracy
This study built a rectangular microstructure for prediction analysis, and the imprint forecast result of the unequal angle was calculated. The predicted position variation was calculated according to a single rectangular microstructure (1 × 1) at rotation angle π/9 • and tilt angle π/6 • , as shown in Figure 18. The calculation prediction of the array rectangular microstructure (4 × 4), where the results of array (4 × 4) at rotation angle 0 • , tilt angle π/6 • and rotation angle π/ 6 • , tilt angle π/6 • , is shown in Figure 19. The dynamic tensile asymmetric imprinting experiment was performed. The array microlens hemispherical structure, as derived from the actual dynamic tensile, was reproduced by rotation angle π/6 • and tilt angle π/36 • . This result is close to the calculation result of MATLAB software, which is enough to show the imprinting stability of this system and the feasibility of the prediction method, as shown in Figure 20.

Experimental Performance of Bifacial Gasbag-Assisted Dynamic Tensile Forming
The proposed system has simultaneous upper and lower cavities, thus, the dynamic tensile and gasbag-assisted imprinting process can be performed simultaneously. Therefore, the simultaneous forming process of the bifacial microstructure can be performed. For the array circular microstructure diameter of 180 μm, microstructure height of 130 μm feature dimensions, the dynamic tensile of 1 mm, 3 mm, 5 mm, 7 mm, 9 mm, and 11 mm in various directions were regulated. The gasbag-assisted imprint was measured by a surface profiler. The result showed that the imprint transcription rate of this system has accurate transfer printing height (overall average transcription rate: 96.985%) and transfer printing diameter (overall average transcription rate: 97.34%), and the transfer printing formability 95% confidence interval is established, as shown in Figures 21 and 22.   100μm 100μm

Experimental Performance of Bifacial Gasbag-Assisted Dynamic Tensile Forming
The proposed system has simultaneous upper and lower cavities, thus, the dynamic tensile and gasbag-assisted imprinting process can be performed simultaneously. Therefore, the simultaneous forming process of the bifacial microstructure can be performed. For the array circular microstructure diameter of 180 µm, microstructure height of 130 µm feature dimensions, the dynamic tensile of 1 mm, 3 mm, 5 mm, 7 mm, 9 mm, and 11 mm in various directions were regulated. The gasbag-assisted imprint was measured by a surface profiler. The result showed that the imprint transcription rate of this system has accurate transfer printing height (overall average transcription rate: 96.985%) and transfer printing diameter (overall average transcription rate: 97.34%), and the transfer printing formability 95% confidence interval is established, as shown in Figures 21 and 22.

Forming Operation Window of Bifacial Gasbag-Assisted Dynamic Tensile Imprint Reproduction
Regarding the array circular microstructure of the cross microstructure mold in the upper and lower cavities (diameter of 180 m, microstructure height of 130 m feature dimension), the dynamic tensile was performed to adjust 1 mm, 3 mm, 5 mm, 7 mm, 9 mm, and 11 mm in various directions. The gasbag-assisted imprinting process reproduction forming operation window was established by multiple experiments. The experimental results showed that, under the axial operating conditions, the gas-assisted pressure within 0.08Mpa~0.15Mpa has good reproduction formability, as shown in Figure 22. Dynamic tensile and gasbag-assisted imprinting process microstructure diameter imprint transcription rate.

Forming Operation Window of Bifacial Gasbag-Assisted Dynamic Tensile Imprint Reproduction
Regarding the array circular microstructure of the cross microstructure mold in the upper and lower cavities (diameter of 180 m, microstructure height of 130 m feature dimension), the dynamic tensile was performed to adjust 1 mm, 3 mm, 5 mm, 7 mm, 9 mm, and 11 mm in various directions. The gasbag-assisted imprinting process reproduction forming operation window was established by multiple experiments. The experimental results showed that, under the axial operating conditions, the gas-assisted pressure within 0.08Mpa~0.15Mpa has good reproduction formability, as shown in Figure 23.

Conclusions
This study proposed a dynamic microstructure forming method for flexible polymer, and discussed the time-dependent creep behavior of dependent variables in the dynamic tensile process. Maxwell, Kelvin-Voight, Model A, and Model B were built for the cross microstructure mold. The findings showed that Model A has the steadiest performance. In terms of mold imprint uniformity, the flexible polymer mold imprint uniformity, mold pressure (stress, strain), and mechanical stress (von Mises stress and Tresca stress) effects in the tensile process were analyzed by Abaqus and ANSYS. The results showed that the system is uniform and stable. In terms of curved imprinting, for unequal curved imprinting, the geometric feature variation of the microstructure was calculated by MATLAB software. The results showed that the transverse deformation and longitudinal deformation are related to the location of the mold. In terms of asymmetric imprinting, an array cylinder with a diameter of 250um and a height of 75um was implemented. Different extensions were adjusted to find the similarities and differences in the forming of adjacent microstructures under curved imprinting. The single rectangular microstructure (1 × 1), array rectangular microstructure (4 × 4), array cylindrical microstructure, and array microlens hemispherical structure were calculated by asymmetric imprinting (MATLAB), as well as unequal rotation angle and tilt angle forming features, which are compared with the experiment to verify the feasibility of the prediction method. The result of imprint reproduction at rotation angle π/6° and tilt angle π/36° is close to the calculation result of the MATLAB software.

Conclusions
This study proposed a dynamic microstructure forming method for flexible polymer, and discussed the time-dependent creep behavior of dependent variables in the dynamic tensile process. Maxwell, Kelvin-Voight, Model A, and Model B were built for the cross microstructure mold. The findings showed that Model A has the steadiest performance. In terms of mold imprint uniformity, the flexible polymer mold imprint uniformity, mold pressure (stress, strain), and mechanical stress (von Mises stress and Tresca stress) effects in the tensile process were analyzed by Abaqus and ANSYS. The results showed that the system is uniform and stable. In terms of curved imprinting, for unequal curved imprinting, the geometric feature variation of the microstructure was calculated by MATLAB software. The results showed that the transverse deformation and longitudinal deformation are related to the location of the mold. In terms of asymmetric imprinting, an array cylinder with a diameter of 250um and a height of 75um was implemented. Different extensions were adjusted to find the similarities and differences in the forming of adjacent microstructures under curved imprinting. The single rectangular microstructure (1 × 1), array rectangular microstructure (4 × 4), array cylindrical microstructure, and array microlens hemispherical structure were calculated by asymmetric imprinting (MATLAB), as well as unequal rotation angle and tilt angle forming features, which are compared with the experiment to verify the feasibility of the prediction method. The result of imprint reproduction at rotation angle π/6 • and tilt angle π/36 • is close to the calculation result of the MATLAB software.