Automated optimization of polymer extrusion dies using finite element analysis

An extrusion die is used to continuously produce parts with a constant cross section; such as sheets, pipes, tire components and more complex shapes such as window seals. The die is fed by a screw extruder when polymers are used. The extruder melts, mixes and pressures the material by the rotation of either a single or double screw. The polymer can then be continuously forced through the die producing a long part in the shape of the die outlet. The extruded section is then cut to the desired length. Generally, the primary target of a well designed die is to produce a uniform outlet velocity without excessively raising the pressure required to extrude the polymer through the die. Other properties such as temperature uniformity and residence time are also important but are not directly considered in this work. Designing dies for optimal outlet velocity variation using simple analytical equations are feasible for basic die geometries or simple channels. Due to the complexity of die geometry and of polymer material properties design of complex dies by analytical methods is difficult. For complex dies iterative methods must be used to optimize dies. An automated iterative method is desired for die optimization. To automate the design and optimization of an extrusion die two issues must be dealt with. The first is how to generate a new mesh for each iteration. In this work, this is approached by modifying a Parasolid file that describes a CAD part. This file is then used in a commercial meshing software. Skewing the initial mesh to produce a new geometry was also employed as a second option. The second issue is an optimization problem with the presence of noise stemming from variations in the mesh and cumulative truncation errors. In this work a simplex method and a modified trust region method were employed for automated optimization of die geometries. For the trust region a discreet derivative and a BFGS Hessian approximation were used. To deal with the noise in the function the trust region method was modified to automatically adjust the discreet derivative step size and the trust region based on changes in noise and function contour. Generally uniformity of velocity at exit of the extrusion die can be improved by increasing resistance across the die but this is limited by the pressure capabilities of the extruder. In optimization, a penalty factor that increases exponentially from the pressure limit is applied. This penalty can be applied in two different ways; the first only to the designs which exceed the pressure limit, the second to both designs above and below the pressure limit. Both of these methods were tested and compared in this work.

High cycle fatigue of AA6082 and AA6063 aluminum extrusions

The high cycle fatigue behavior of hollow extruded AA6082 and AA6063 aluminum extrusions has been studied. Hollow extruded aluminum profiles can be processed into intricate shapes, and may be suitable replacements for fatigue critical automotive applications requiring reduced weight. There are several features inherent in hollow aluminum extrusions, such as seam welds, charge welds, microstructural variations and die lines. The effects of such extrusion variables on high cycle fatigue properties were studied by taking specimens from an actual car bumper extrusion. It appears that extrusion die lines create large anisotropy differences in fatigue properties, while welds themselves have little effect on fatigue lives. Removal of die lines greatly increased fatigue properties of AA6082 specimens taken transverse to the extrusion direction. Without die lines, anisotropy in fatigue properties between AA6082 specimens taken longitudinal and transverse to the extrusion direction, was significantly reduced, and properties associated with the orientation of the microstructure appears to be isotropic. A fibrous microstructure for AA6082 specimens showed great improvements in fatigue behavior. The effects of elevated temperatures and exposure of specimens to NaCl solutions was also studied. Exposure to the salt solution greatly reduced the fatigue lives of specimens, while elevated temperatures showed more moderate reductions in fatigue lives.

Three dimensional finite element simulation of polymer melting and flow in a single-screw extruder : optimization of screw channel geometry

Single-screw extrusion is one of the widely used processing methods in plastics industry, which was the third largest manufacturing industry in the United States in 2007 [5]. In order to optimize the single-screw extrusion process, tremendous efforts have been devoted for development of accurate models in the last fifty years, especially for polymer melting in screw extruders. This has led to a good qualitative understanding of the melting process; however, quantitative predictions of melting from various models often have a large error in comparison to the experimental data. Thus, even nowadays, process parameters and the geometry of the extruder channel for the single-screw extrusion are determined by trial and error. Since new polymers are developed frequently, finding the optimum parameters to extrude these polymers by trial and error is costly and time consuming. In order to reduce the time and experimental work required for optimizing the process parameters and the geometry of the extruder channel for a given polymer, the main goal of this research was to perform a coordinated experimental and numerical investigation of melting in screw extrusion. In this work, a full three-dimensional finite element simulation of the two-phase flow in the melting and metering zones of a single-screw extruder was performed by solving the conservation equations for mass, momentum, and energy. The only attempt for such a three-dimensional simulation of melting in screw extruder was more than twenty years back. However, that work had only a limited success because of the capability of computers and mathematical algorithms available at that time. The dramatic improvement of computational power and mathematical knowledge now make it possible to run full 3-D simulations of two-phase flow in single-screw extruders on a desktop PC. In order to verify the numerical predictions from the full 3-D simulations of two-phase flow in single-screw extruders, a detailed experimental study was performed. This experimental study included Maddock screw-freezing experiments, Screw Simulator experiments and material characterization experiments. Maddock screw-freezing experiments were performed in order to visualize the melting profile along the single-screw extruder channel with different screw geometry configurations. These melting profiles were compared with the simulation results. Screw Simulator experiments were performed to collect the shear stress and melting flux data for various polymers. Cone and plate viscometer experiments were performed to obtain the shear viscosity data which is needed in the simulations. An optimization code was developed to optimize two screw geometry parameters, namely, screw lead (pitch) and depth in the metering section of a single-screw extruder, such that the output rate of the extruder was maximized without exceeding the maximum temperature value specified at the exit of the extruder. This optimization code used a mesh partitioning technique in order to obtain the flow domain. The simulations in this flow domain was performed using the code developed to simulate the two-phase flow in single-screw extruders.

Erupted magma volume estimates at Santiaguito and Pacaya Volcanoes, Guatemala using digital elevation models

High resolution digital elevation models (DEMs) of Santiaguito and Pacaya volcanoes, Guatemala, were used to estimate volume changes and eruption rates between 1954 and 2001. The DEMs were generated from contour maps and aerial photography, which were analyzed in ArcGIS 9.0®. Because both volcanoes were growing substantially over the five decade period, they provide a good data set for exploring effective methodology for estimating volume changes. The analysis shows that the Santiaguito dome complex grew by 0.78 ± 0.07 km3 (0.52 ± 0.05 m3 s-1) over the 1954-2001 period with nearly all the growth occurring on the El Brujo (1958-75) and Caliente domes (1971-2001). Adding information from field data prior to 1954, the total volume extruded from Santiaguito since 1922 is estimated at 1.48 ± 0.19 km3. Santiaguito’s growth rate is lower than most other volcanic domes, but it has been sustained over a much longer period and has undergone a change toward more exogenous and progressively slower extrusion with time. At Santiaguito some of the material being added at the dome is subsequently transported downstream by block and ash flows, mudflows and floods, creating channel shifting and areas of aggradation and erosion. At Pacaya volcano a total volume of 0.21 ± 0.05 km3 was erupted between 1961 and 2001 for an average extrusion rate of 0.17 ± 0.04 m3 s-1. Both the Santiaguito and Pacaya eruption rate estimates reported here are minima, because they do not include estimates of materials which are transported downslope after eruption and data on ashfall which may result in significant volumes of material spread over broad areas. Regular analysis of high resolution DEMs using the methods outlined here, would help quantify the effects of fluvial changes to downstream populated areas, as well as assist in tracking hazards related to dome collapse and eruption.

SO2 emissions at Volcan de Colima, 2003-2007

Volcán de Colima has been continuously erupting since the onset of dome growth in 1998. This period of unrest has had 4 prominent periods; 1998-1999, 2003, 2004-2005, and the current dome growth that began in February of 2007. Each of these episodes was marked by lava extrusion forming a dome and lava flows, followed by explosions that destroyed the dome. The Correlation Spectrometer (COSPEC) was used to determine SO2 emission rates on 164 days from May 2003 to February 2007, using both stationary ground based scans and some flight traverses. Scans were separated into the categories of explosive degassing and passive, or background degassing. These scans show variation in the SO2 flow rate from below detection limit (~3 t/d depending on environmental conditions) during background, passive emissions to a peak of 2949 t/d (34 kilograms/second) during an explosion on 9 October, 2004. Both passive and explosive degassing increased when there was lava extrusion in 2004 and with the increased explosive activity in 2005. These two different processes of degassing wax with each other when activity increases and wane together as well, indicating a parallel cyclicity in the volcanic eruption and degassing rates, where the conduit partially seals (pressurizes) between explosions. Colima’s gas and eruptive behavior is compared to similar systems such as Santiaguito and Soufrière Hills, Montserrat. About 2/3 of Colima’s SO2 degassing, amounting to 1.3 x 105 tonnes in 3.74 yrs has come in short lived small (VEI=0-1) vertical explosions that occurred at the rate of 100-3000explosions/ month, and the remaining third has occured in continuous passive degassing. Colima emits sulfur at a rate equivalent to about 0.04 to 0.08 wt % S, similar to other andesitic convergent plate boundary volcanoes. There has been an explosive destruction of the dome in every cycle for that past 5 years, and it is assumed that the current dome which began growth in February, 2007 (just at the end of this study) will be destroyed. Higher emission rates seen in the quiescence of 2006 may have eased the pressure at the time, resulting in the slow effusion of the current dome and lack of explosivity.

Intergranular Corrosion and Stress Corrosion Cracking of Extruded AA6005A

A research program focused on understanding the intergranular corrosion (IGC) and stress corrosion cracking (SCC) behavior of AA6005A aluminum extrusions is presented in this dissertation. The relationship between IGC and SCC susceptibility and the mechanisms of SCC in AA6005A extrusions were studied by examining two primary hypotheses. IGC susceptibility of the elongated grain structure in AA6005A exposed to low pH saltwater was found to depend primarily on the morphology of Cu-containing precipitates adjacent to the grain boundaries in the elongated grain structure. IGC susceptibility was observed when a continuous (or semi-continuous) film of Cu-containing phase was present along the grain boundaries. When this film coarsened to form discrete Cu-rich precipitates, no IGC was observed. The morphology of the Cu-rich phase depended on post-extrusion heat treatment. The rate of IGC penetration in the elongated grain structure of AA6005A-T4 and AA6005A-T6 extrusions was found to be anisotropic with IGC propagating most rapidly along the extrusion direction, and least rapidly along the through thickness direction. A simple 3-dimensional geometric model of the elongated grain structure was accurately described the observed IGC anisotropy, therefore it was concluded that the anisotropic IGC susceptibility in the elongated grain structure was primarily due to geometric elongation of the grains. The velocity of IGC penetration along all directions in AA6005A-T6 decreased with exposure time. Characterization of the local environment within simulated corrosion paths revealed that a pH gradient existed between the tip of the IGC path and the external environment. Knowledge of the local environment within an IGC path allowed development of a simple model based on Fick's first law that considered diffusion of Al3+ away from the tip of the IGC path. The predicted IGC velocity agreed well with the observed IGC velocity, therefore it was determined that diffusion of Al3+ was the primary factor in determining the velocity of IGC penetration. The velocity of crack growth in compact tensile (CT) specimens of AA6005A-T6 extrusion exposed to 3.5% NaCl at pH = 1.5 was nearly constant over a range of applied stress intensities, exposure times, and crack lengths. The crack growth behavior of CT specimens of AA6005A-T6 extrusion exposed to a solution of 3.5% NaCl at pH = 2.0 exhibited similar behavior, but the crack velocity was ~10.5X smaller than that those exposed to a solution at pH =1.5. Analysis of the local stress state and polarization behavior at the crack tip predicted that increasing the pH of the bulk solution from 1.5 to 2.0 would decrease the corrosion current density at the crack tip by approximately 11.8X. This predicted decrease in corrosion current density was in reasonable agreement with the observed decrease in SCC velocity associated with increasing the solution pH from 1.5 to 2.0. The agreement between the predicted and observed SCC velocities suggested that the electrochemical reactions controlling SCC in AA6005A-T6 extrusions are ultimately controlled by the pH gradient that exists between the crack tip and external environment.