Tectonic and geological processes on Earth often result in structural anisotropy of the subsurface, which can be imaged by various geophysical methods. In order to achieve appropriate and realistic Earth models for interpretation, inversion algorithms have to allow for an anisotropic subsurface. Within the framework of this thesis, I analyzed a magnetotelluric (MT) data set taken from the Cape Fold Belt in South Africa. This data set exhibited strong indications for crustal anisotropy, e.g. MT phases out of the expected quadrant, which are beyond of fitting and interpreting with standard isotropic inversion algorithms. To overcome this obstacle, I have developed a two-dimensional inversion method for reconstructing anisotropic electrical conductivity distributions. The MT inverse problem represents in general a non-linear and ill-posed minimization problem with many degrees of freedom: In isotropic case, we have to assign an electrical conductivity value to each cell of a large grid to assimilate the Earth's subsurface, e.g. a grid with 100 x 50 cells results in 5000 unknown model parameters in an isotropic case; in contrast, we have the sixfold in an anisotropic scenario where the single value of electrical conductivity becomes a symmetric, real-valued tensor while the number of the data remains unchanged. In order to successfully invert for anisotropic conductivities and to overcome the non-uniqueness of the solution of the inverse problem it is necessary to use appropriate constraints on the class of allowed models. This becomes even more important as MT data is not equally sensitive to all anisotropic parameters. In this thesis, I have developed an algorithm through which the solution of the anisotropic inversion problem is calculated by minimization of a global penalty functional consisting of three entries: the data misfit, the model roughness constraint and the anisotropy constraint. For comparison, in an isotropic approach only the first two entries are minimized. The newly defined anisotropy term is measured by the sum of the square difference of the principal conductivity values of the model. The basic idea of this constraint is straightforward. If an isotropic model is already adequate to explain the data, there is no need to introduce electrical anisotropy at all. In order to ensure successful inversion, appropriate trade-off parameters, also known as regularization parameters, have to be chosen for the different model constraints. Synthetic tests show that using fixed trade-off parameters usually causes the inversion to end up by either a smooth model with large RMS error or a rough model with small RMS error. Using of a relaxation approach on the regularization parameters after each successful inversion iteration will result in smoother inversion model and a better convergence. This approach seems to be a sophisticated way for the selection of trade-off parameters. In general, the proposed inversion method is adequate for resolving the principal conductivities defined in horizontal plane. Once none of the principal directions of the anisotropic structure is coincided with the predefined strike direction, only the corresponding effective conductivities, which is the projection of the principal conductivities onto the model coordinate axes direction, can be resolved and the information about the rotation angles is lost. In the end the MT data from the Cape Fold Belt in South Africa has been analyzed. The MT data exhibits an area (> 10 km) where MT phases over 90 degrees occur. This part of data cannot be modeled by standard isotropic modeling procedures and hence can not be properly interpreted. The proposed inversion method, however, could not reproduce the anomalous large phases as desired because of losing the information about rotation angles. MT phases outside the first quadrant are usually obtained by different anisotropic anomalies with oblique anisotropy strike. In order to achieve this challenge, the algorithm needs further developments. However, forward modeling studies with the MT data have shown that surface highly conductive heterogeneity in combination with a mid-crustal electrically anisotropic zone are required to fit the data. According to known geological and tectonic information the mid-crustal zone is interpreted as a deep aquifer related to the fractured Table Mountain Group rocks in the Cape Fold Belt. / Tektonische und geologische Prozesse verursachen häufig eine strukturelle Anisotropie des Untergrundes, welche von verschiedenen geophysikalischen Methoden beobachtet werden kann. Zur Erstellung und Interpretation geeigneter, realistischer Modelle der Erde sind Inversionsalgorithmen notwendig, die einen anisotropen Untergrund einbeziehen können. Für die vorliegende Arbeit habe ich einen magnetotellurischen (MT) Datensatz vom Cape Fold Gürtel in Südafrika untersucht. Diese Daten weisen auf eine ausgeprägte Anisotropie der Kruste hin, da z.B. die MT Phasen außerhalb des erwarteten Quadranten liegen und nicht durch standardisierte isotrope Inversionsalgorithmen angepasst und ausgewertet werden können. Um dieses Problem zu beheben, habe ich eine zweidimensionale Inversionsmethode entwickelt, welche eine anisotrope elektrische Leitfähigkeitsverteilungen in den Modellen zulässt. Die MT Inversion ist im allgemeinen ein nichtlineares, schlecht gestelltes Minimierungsproblem mit einer hohen Anzahl an Freiheitsgraden. Im isotropen Fall wird jeder Gitterzelle eines Modells ein elektrischer Leitfähigkeitswert zugewiesen um den Erduntergrund nachzubilden. Ein Modell mit beispielsweise 100 x 50 Zellen besitzt 5000 unbekannte Modellparameter. Im Gegensatz dazu haben wir im anisotropen Fall die sechsfache Anzahl, da hier aus dem einfachen Zahlenwert der elektrischen Leitfähigkeit ein symmetrischer, reellwertiger Tensor wird, wobei die Anzahl der Daten gleich bleibt. Für die erfolgreiche Inversion von anisotropen Leitfähigkeiten und um die Nicht-Eindeutigkeit der Lösung des inversen Problems zu überwinden, ist eine geeignete Einschränkung der möglichen Modelle absolut notwendig. Dies wird umso wichtiger, da die Sensitivität von MT Daten nicht für alle Anisotropieparameter gleich ist. In der vorliegenden Arbeit habe ich einen Algorithmus entwickelt, welcher die Lösung des anisotropen Inversionsproblems unter Minimierung einer globalen Straffunktion berechnet. Diese besteht aus drei Teilen: der Datenanpassung, den Zusatzbedingungen an die Glätte des Modells und die Anisotropie. Im Gegensatz dazu werden beim isotropen Fall nur die ersten zwei Parameter minimiert. Der neu definierte Anisotropieterm wird mit Hilfe der Summe der quadratischen Abweichung der Hauptleitfähigkeitswerte des Modells gemessen. Die grundlegende Idee dieser Zusatzbedingung ist einfach. Falls ein isotropes Modell die Daten ausreichend gut anpassen kann, wird keine elektrische Anisotropie zusätzlich in das Modell eingefügt. Um eine erfolgreiche Inversion zu garantieren müssen geeignete Regularisierungsparameter für die verschiedenen Nebenbedingungen an das Modell gewählt werden. Tests mit synthetischen Modellen zeigen, dass bei festgesetzten Regularisierungsparametern die Inversion meistens entweder in einem glatten Modell mit hohem RMS Fehler oder einem groben Modell mit kleinem RMS Fehler endet. Die Anwendung einer Relaxationsbedingung auf die Regularisierung nach jedem Iterationsschritt resultiert in glatteren Inversionsmodellen und einer höheren Konvergenz und scheint ein ausgereifter Weg zur Wahl der Parameter zu sein. Die vorgestellte Inversionsmethode ist im allgemeinen in der Lage die Hauptleitfähigkeiten in der horizontalen Ebene zu finden. Wenn keine der Hauptrichtungen der Anisotropiestruktur mit der vorgegebenen Streichrichtung übereinstimmt, können nur die dazugehörigen effektiven Leitfähigkeiten, welche die Projektion der Hauptleitfähigkeiten auf die Koordinatenachsen des Modells darstellen, aufgelöst werden. Allerdings gehen die Informationen über die Rotationswinkel verloren. Am Ende meiner Arbeit werden die MT Daten des Cape Fold Gürtels in Südafrika analysiert. Die MT Daten zeigen in einem Abschnitt des Messprofils (> 10 km) Phasen über 90 Grad. Dieser Teil der Daten kann nicht mit herkömmlichen isotropen Modellierungsverfahren angepasst und daher mit diesen auch nicht vollständig ausgewertet werden. Die vorgestellte Inversionsmethode konnte die außergewöhnlich hohen Phasenwerte nicht wie gewünscht im Inversionsergebnis erreichen, was mit dem erwähnten Informationsverlust der Rotationswinkel begründet werden kann. MT Phasen außerhalb des ersten Quadranten können für gewöhnlich bei Anomalien mit geneigter Streichrichtung der Anisotropie gemessen werden. Um diese auch in den Inversionsergebnissen zu erreichen ist eine Weiterentwicklung des Algorithmus notwendig. Vorwärtsmodellierungen des MT Datensatzes haben allerdings gezeigt, dass eine hohe Leitfähigkeitsheterogenität an der Oberfläche in Kombination mit einer Zone elektrischer Anisotropie in der mittleren Kruste notwendig sind um die Daten anzupassen. Aufgrund geologischer und tektonischer Informationen kann diese Zone in der mittleren Kruste als tiefer Aquifer interpretiert werden, der im Zusammenhang mit den zerrütteten Gesteinen der Table Mountain Group des Cape Fold Gürtels steht.

The goal of this research is to develop a robust, continuum-based approach for a three-dimensional, Finite Element Method (FEM) simulation of layered geomaterials. There are two main approaches to the numerical modeling of layered geomaterials; discrete or discontinuous techniques and an equivalent continuum concept. In the discontinuous methodology, joints are explicitly simulated. Naturally, discrete techniques provide a more accurate description of discontinuous materials. However, they are complex and necessitate care in modeling of the interface. Also, in many applications, the definition of the input model becomes impractical as the number of joints becomes large. In order to overcome the difficulties associated with discrete techniques, a continuum-based approach has become popular in some application areas. When using a continuum model, a discrete material is replaced by a homogenized continuous material, also known as an 'equivalent continuum'. This leads to a discretization that is independent of both the orientation and spacing of layer boundaries. However, if the layer thickness (i.e., internal length scale of the problem) is large, the classical continuum approach which neglects the effect of internal characteristic length can introduce large errors into the solution. In this research, a full 3D FEM formulation for the elasto-plastic modeling of layered geomaterials is proposed within the framework of Cosserat theory. The effect of the bending stiffness of the layers is incorporated in the matrix of elastic properties. Also, a multi-surface plasticity model, which allows for plastic deformation of both the interfaces between the layers and intact material, is introduced. The model is verified against analytical solutions, discrete numerical models, and experimental data. It is shown that the FEM Cosserat formulation can achieve the same level of accuracy as discontinuous models in predicting the displacements of a layered material with a periodic microstructure. Furthermore, the method is capable of reproducing the strength behaviour of materials with one or more sets of joints. Finally, due to the incorporation of layer thickness into the constitutive model, the FEM Cosserat formulation is capable of capturing complicated failure mechanisms such as the buckling of individual layers of material which occur in stratified media.

No description available.

No description available.

In this Thesis the magnetic and structural properties of thin epitaxial Fe, Co, and Ni films are discussed. Some of the in-situ prepared samples were used to characterize the degree of circular polarization of the newly installed beamline D1011 at MAX-lab. By means of x-ray magnetic circular dichroism (XMCD) and utilizing the associated magneto optic sum rules, the orbital (ml) and spin (ms) moments are determined directly in mB/atom with elemental specificity. The extended x-ray absorption fine structure (EXAFS) measurements yield site specific information on the local crystallographic structure. These measurements were performed using the circular x-rays of several beamlines. The influence of the degree of spatial source coherence lspat of the x-rays was characterized by means of Fresnel diffractometry. A correlation between enhanced XAS white line intensities and higher values of lspat was established for 20 ML Fe, Co, and Ni films on Cu(100). The degree of circularly polarized x-rays (Pc) at beamline D1011 at MAX-lab was characterized by studying Fe films on Cu(100) by means of XMCD. The maximum value of Pc is experimentally determined to Pc =0.85. The Au/Co/Au trilayer system was studied as a function of Co thickness, temperature, and Au cap thickness. A 10 mono-layer (ML) Co film, with an Au cap of 20 Å, shows a spin reorientation transition (SRT) from an in-plane to an out-of-plane easy direction as the temperature is lowered from 300 K to 200 K. The magnetic properities of these Co films are very different to what is found for bulk samples due to, in particular, the broken symmetry at the interfaces. The thickness dependent spin reorientation transition in the Fe/Ag(100) system was characterized by means of XMCD and EXAFS measurements. 3 ML Fe films show an out-of-plane easy direction with an 125% enhanced orbital moment as compared to the 25 ML Fe in-plane film. Simulations of the Fe L-edge EXAFS indicate the bulk Fe bcc structure for film thicknesses of 6-25 ML Fe. For 3 ML Fe strong deviations from this bcc phase is observed. Ultrathin Co films deposited on flat and vicinal Cu(111) in the thickness region 1-25 ML were studied by means of XMCD and scanning tunneling microscopy (STM). The vicinal Cu(111) Co deposition leads to the formation of elongated islands preferentially oriented along the step edges. In connection to this particular Co growth mode we observe an increase of both the orbital and the spin moment on the vicinal Cu(111) of about 25% relative to what was observed for Co on flat Cu(111).

La motivació de la tesis deriva en el interès de la indústria aeronàutica a explotar, per mitjà d'un disseny adaptat, la utilització del gir d'esquerda per protegir els reforços situats davant una esquerda que s'està propagant en la xapa d'una estructura integral. L'objectiu principal és l'avaluació i predicció del gir d'esquerda en situacions de càrrega pròximes a Mode I, proporcionant una eina de modelització i un criteri confident. L'entorn industrial sota el qual s'ha realitzat aquest treball requereix una predicció ràpida del comportament estructural proporcionant informació útil als constructors. Per aquest motiu la predicció del gir d'esquerda s'ha investigat utilitzant la teoria linear elàstica de la mecànica de la fractura (LEFM) i l'anàlisi amb elements finits (FEA).Durant aquest treball s'ha demostrat la importància i necessitat de caracteritzar el camp de tensions a la punta de l'esquerda amb el factor d'intensitat de tensió (SIF) conjuntament amb un segon paràmetre. La tensió uniforme, no singular, normal a la línea de l'esquerda i dependent en la geometria i càrrega de la proveta, es a dir la tensió T, ha estat seleccionada com a segon paràmetre per dur a terme les prediccions del gir d'esquerda.El criteri més desenvolupat per predir el gir d'esquerda en situacions pròximes a Mode I és el proposat per Buczek, Herakovich, Boone et al., anomenat WEFO en la tesis. Aquest combina el criteri de tensió principal màxima amb la tensió T i considera efectes d'anisotropia. LEFM s'ha utilitzat també en la predicció del gir d'esquerda sota càrregues quasi estàtiques controlant en tot moment la plastificació del lligament.En la investigació d'eines de modelització/simulació s'ha tingut en compte les capacitats d'aquestes en el camp de la mecànica de la fractura, de disseny, d'implementació, així com la complexitat d'ús. Tot i que hi ha un gran ventall de Softwares que compleixen els requeriments assenyalats, només aquells que es trobaven a l'abast de l'autor s'han analitzat. StressCheck ha estat escollit com a resultat de la investigació. L'avaluació de la propagació de l'esquerda en provetes compactes en tensió (CT) i en provetes amb dos elements reforçants (2SP) sota els règims de Paris i Forman ha estat satisfactòria.Un pas important ha estat la implementació de la capacitat d'extracció de la tensió T. La demostració de la fiabilitat en el seu càlcul s'ha demostrat mitjançant resultats en la literatura i càlculs analítics en provetes de doble biga en volada (DCB). Un aspecte a tenir en compte és la importància en realitzar anàlisis no linears geomètrics pel càlcul del SIF i la tensió T.Prediccions en la trajectòria de l'esquerda s'han realitzat en base amb els resultats obtinguts en l'estudi de modelització. La millor trajectòria s'ha predit per mitjà del criteri WEFO. No obstant, les diferents trajectòries obtingudes per una esquerda propagant-se en la direcció T-L o L-T no són comparables amb els resultats experimentals.Aquestes deficiències estan relacionades en la definició del punt d'inestabilitat de l'esquerda. Algunes referències posen de manifest que hi ha experiències on l'esquerda es comporta de forma estable tot i mostrar T > 0. Per un altre banda, els criteris WEF i WEFO defineixen la inestabilitat dependent d'una distancia específica del material, rc. Però la seva definició no és única i no existeix cap acord sobre el seu càlcul.L'autor proposa un criteri derivat dels criteris existents i basant-se en els assajos, simulacions i resultats obtinguts. Aquest deriva del treball de Pettit i la tensió T normalitzada, TR, proposada per Pook. La fiabilitat d'aquest criteri es demostra amb la proveta DCB. Les prediccions de la trajectòria de l'esquerda en la proveta cruciforme no són tant satisfactòries. Tot i així, s'ha d'accentuar que el criteri desenvolupat proporciona la predicció més acurada. / The motivation of this thesis started from the interest of aeronautical industry to exploit the utilization of crack turning to protect stiffeners in front of an approaching skin-crack in integral structures by a tailored design. The main objective was to assess and predict crack turning under nearly Mode I situations on structures that reproduce aeronautical conditions by providing a modelling tool and a reliable criterion. The industrial environmental in which this work has been carried out requires a fast prediction of the structural behaviour to provide useful inputs to aircraft designers. It is for this reason that the crack turning prediction was investigated by means of LEFM and FEA. During this work it has been shown the importance and necessity of a second parameter for the characterisation of the stress field at the crack tip besides the SIF. Among the different proposed second parameters, the uniform non-singular stress, normal to the crack line and dependent on the type of loading and specimen geometry, i.e. the T-stress, was selected for crack turning predictions due to both calculation simplicity and its independence of the crack tip distance. The most developed criterion for crack turning predictions near Mode I loading is the criterion proposed by Buczek, Herakovich and Boone et al., called the WEFO-criterion. This is the Maximal Principal Stress criterion implemented with the T-stress and taking into account anisotropic effects. A challenge of this thesis was to overcome the lack of prediction on crack turning provided by this last criterion.Although the validity of LEFM is restricted, it was applied for the prediction of crack turning for quasi-static loading while paying attention to possible plastification. A screening of existent commercial and non commercial tools was carried out in respect to their fracture mechanics capabilities, their design abilities, implementation as well as their complexity. Although, there are many software possibilities, only those within the reach of the author were evaluated. This resulted in the selection of the commercial tool StressCheck®. The assessment of crack propagation on compact tension and two stringer specimens governed by the Paris and Forman regimes was satisfactory compared with experimental results using the material data from simple standard specimens.An important step was the implementation of the T-stress extraction facility in the tool and the evidence of its reliability. The latter was proved by literature and analytical calculations on DCB specimens. An important finding was the importance to perform geometric non-linear analyses for computing SIF and T-stress to find values comparable with literature data and analytical calculations. Taking into account the results obtained on the modelling study, crack path predictions were performed. The best prediction by means of existing criteria was reached by the WEFO-criterion. Different crack paths were predicted for a crack propagating in T-L or L-T directions. However, these predictions were not satisfactorily reliable: the point in the crack path where crack turning should take place was not predicted adequately. Additionally, the crack paths were similar for T L and L-T directions. These deficiencies are related with the definition of the crack path instability. Some literature results have shown that in some experiences the crack behaved in a stable manner even if T > 0. Moreover, WEF and WEFO criteria define crack instability to be related with a material specific distance, rc, but, there is no agreement about its definition. Based on tests, simulation results and observations noted during this work, a compilation criterion was proposed. This is based on the work of Pettit and the normalised T-stress, TR, proposed by Pook. Its reliability was successfully proved on the DCB. The crack path predictions on the CFS were not as satisfactory. But even at its worst the developed criterion was the most accurate.

There exists considerable motivation to reduce vehicle weight through the adoption of lightweight materials, such as aluminum alloys, while maintaining energy absorption and component integrity under crash conditions. To this end, it is of particular interest to study the crash behaviour of lightweight tubular hydroformed structures to determine how the forming behaviour affects the axial crush response. Thus, the current research has studied the dynamic crush response of both non-hydroformed and hydroformed EN-AW 5018 and AA5754 aluminum alloy tubes using both experimental and numerical methods. Experiments were performed in which hydroforming process parameters were varied in a parametric fashion after which the crash response was measured. Experimental parameters included the tube thickness and the hydroformed corner radii of the tubes. Explicit dynamic finite element simulations of the hydroforming and crash events were carried out with particular attention to the transfer of forming history from the hydroforming simulations to the crash models. The results showed that increases in the strength of the material due to work hardening during hydroforming were beneficial in increasing energy absorption during crash. However, it was shown that thinning in the corners of the tube during hydroforming decreased the energy absorption capabilities during axial crush. Residual stresses resulting from hydroforming had little effect on the energy absorption characteristics during axial crush. The current research has shown that, in addition to capturing the forming history in the crash models, it is also important to account for effects of material non-linearity such as kinematic hardening, anisotropy, and strain-rate effects in the finite element models. A model combining a non-linear kinematic hardening model, the Johnson-Cook rate sensitive model, and the Yld2000-2d anisotropic model was developed and implemented in the finite element simulations. This combined model did not account for the effect of rotational hardening (plastic spin) due to plastic deformation. It is recommended that a combined constitutive model, such as the one described in this research, be utilized for the finite element study of materials that show sensitivity to the Bauschinger effect, strain-rate effects, and anisotropy.

In this thesis the magnetic and structural properties of ultra-thin 3d transition metals films have been investigated, in particular Fe, Ni and Co films. X-ray Magnetic Circular Dichroism (XMCD) has provided element specific spin (ms ) and orbital (ml ) moments per atom by utilizing the magneto optic sum-rules. Element specific hysteresis curves have been measured by means of X-ray Resonant Magnetic Scattering (XRMS), and the local crystallographic structure has been investigated using Extended X-ray Absorption Fine Structure (EXAFS). By performing XMCD on Fe/Ag(100) we observe a spin reorientation from in-plane to out-of-plane as the Fe thickness is lowered. At temperatures below 300K it occurs around 5-7 mono layers (ML) of Fe. While reorienting the magnetization out-of-plane the orbital moment increases with 125% but only a minor increase (5%) of the spin moment is observed. Extended X-ray Absorption Fine Structure (EXAFS) measurements indicate that films 6 ML and thicker have a bulk-like bcc structure. For the thin out-of-plane films, the local crystallographic structure is more complicated. The spin reorientation of the Au/Co/Au tri-layer system has been studied as a function of temperature, Co layer and Au cap thickness. An unexpected behavior of the orbital moment upon spin reorientation is found in these systems. An ex-situ prepared sample shows a smooth spin reorientation from an in-plane to an out-of-plane easy magnetization direction as the temperature is lowered from 300K to 200K. In-situ prepared samples have also been investigated and a novel phase diagram has been identified. The Au/Co interface has been explored during the Au capping by means of photoemission measurements. In the bi- and tri-layer system of Fe and Ni we have been able to manipulate the spin reorientation by varying the Fe and Ni thickness. A novel non-collinear interlayer exchange interaction for 3d ferro magnets in direct contact has been discovered for a set of samples. This exchange interaction is found to be strongly dependant on the preparation conditions.

There exists considerable motivation to reduce vehicle weight through the adoption of lightweight materials, such as aluminum alloys, while maintaining energy absorption and component integrity under crash conditions. To this end, it is of particular interest to study the crash behaviour of lightweight tubular hydroformed structures to determine how the forming behaviour affects the axial crush response. Thus, the current research has studied the dynamic crush response of both non-hydroformed and hydroformed EN-AW 5018 and AA5754 aluminum alloy tubes using both experimental and numerical methods. Experiments were performed in which hydroforming process parameters were varied in a parametric fashion after which the crash response was measured. Experimental parameters included the tube thickness and the hydroformed corner radii of the tubes. Explicit dynamic finite element simulations of the hydroforming and crash events were carried out with particular attention to the transfer of forming history from the hydroforming simulations to the crash models. The results showed that increases in the strength of the material due to work hardening during hydroforming were beneficial in increasing energy absorption during crash. However, it was shown that thinning in the corners of the tube during hydroforming decreased the energy absorption capabilities during axial crush. Residual stresses resulting from hydroforming had little effect on the energy absorption characteristics during axial crush. The current research has shown that, in addition to capturing the forming history in the crash models, it is also important to account for effects of material non-linearity such as kinematic hardening, anisotropy, and strain-rate effects in the finite element models. A model combining a non-linear kinematic hardening model, the Johnson-Cook rate sensitive model, and the Yld2000-2d anisotropic model was developed and implemented in the finite element simulations. This combined model did not account for the effect of rotational hardening (plastic spin) due to plastic deformation. It is recommended that a combined constitutive model, such as the one described in this research, be utilized for the finite element study of materials that show sensitivity to the Bauschinger effect, strain-rate effects, and anisotropy.

An experimental and numerical study of the isothermal and non-isothermal warm formability of an AA3003 aluminum alloy brazing sheet is presented. Forming limit diagrams were determined using warm limiting dome height (LDH) experiments with in situ strain measurement based on digital image correlation (DIC) techniques. Forming limit curves (FLCs) were developed at several temperature levels (room temperature, 100ºC, 200ºC, 250ºC, and 300ºC) and strain-rates (0.003, 0.018, and 0.1s-1). The formability experiments demonstrated that temperature has a significant effect on formability, whereas forming speed has a mild effect within the studied range. Elevating the temperature to 250C improved the formability more than 200% compared to room temperature forming, while forming at lower speeds increased the limiting strains by 10% and 17% at room temperature and 250ºC, respectively. Non-isothermal deep draw experiments were developed considering an automotive heat exchanger plate. A parametric study of the effects of die temperature, punch speed, and blank holder force on the formability of the part was conducted. The introduction of non-isothermal conditions in which the punch is cooled and the flange region is heated to 250C resulted in a 61% increase in draw depth relative to room temperature forming. In order to develop effective numerical models of warm forming processes, a constitutive model is proposed for aluminum alloy sheet to account for temperature and strain rate dependency, as well as plastic anisotropy. The model combines the Barlat YLD2000 yield criterion (Barlat et al., 2003) to capture sheet anisotropy and the Bergstrom (1982) hardening rule to account for temperature and strain rate dependency. Stress-strain curves for AA3003 aluminum alloy brazing sheet tested at elevated temperatures and a range of strain rates were used to fit the Bergstrom parameters, while measured R-values were used to fit the yield function parameters. The combined constitutive model was implemented within a user defined material subroutine that was linked to the LS-DYNA finite element code. Finite element models were developed based on the proposed material model and the results were compared with experimental data. Isothermal uniaxial tensile tests were simulated and the predicted responses were compared with measured data. The tensile test simulations accurately predicted material behaviour. The user material subroutine and forming limit criteria were then applied to simulate the isothermal warm LDH tests, as well as isothermal and non-isothermal warm deep drawing experiments. Two deep draw geometries were considered, the heat exchanger plate experiments developed as part of this research and the 100 mm cylindrical cup draw experiments performed by McKinley et al. (2010). The strain distributions, punch forces and failure location predicted for all three forming operations were in good agreement with the experimental results. Using the warm forming limit curves, the models were able to accurately predict the punch depths to failure as well as the location of failure initiation for both the isothermal and non-isothermal deep draw operations.