Criação de programa informático em ambiente Matlab® para estudo de estruturas planas pelo Método de Elementos Finitos.

Os modelos a ser analisados pelo Método de Elementos Finitos são cada vez mais complexos e, nos tempos que correm, seria impensável realizar tais análises sem um apoio computorizado. Existe para esta finalidade uma vasta gama de programas que permitem realizar tarefas que passam pelo desenho de estruturas, análise estática de cargas, análise dinâmica e vibrações, visualização do comportamento físico (deformações) em tempo real, que permitem a otimização da estrutura. Sob o pretexto de permitir a qualquer utilizador uma análise de estruturas simples com o Método dos Elementos Finitos, surge esta tese, onde se irá criar de raiz um programa com interface gráfica no ambiente MATLAB® para análise de estruturas simples com dois tipos de elemento finito, triangular de deformação constante e quadrangular de deformação linear. O software desenvolvido, verificado por comparação com um software comercial dedicado para o efeito, efetua malhagem com elementos bidimensionais triangulares e quadriláteros e resolve modelos arbitrados pelo Método de Elementos Finitos, representando estes resultados visualmente e em formato tabular.

Prediction of stress-strain behaviors in steels using an integrated constitutive, FEM and ANN model

Austenitic steels with a carbon content of 0.0037 to 0.79 wt% C are torsion tested and modeled using a physically based constitutive model and an Integrated Phenomenological and Artificial neural Network (IPANN) model. The prediction of both the constitutive and IPANN models on steel 0.017 wt% C is then evaluated using a finite element (FEM) code ABAQUS with different reduction in the thickness after rolling through one roll stand. It is found that during the rolling process, the prediction accuracy of the reaction force from FEM simulation for both constitutive and IPANN models depends on the strain achieved (average reduction in thickness). By integrating FEM into IPANN model and introducing the product of strain and stress as an input of the ANN model, the accuracy of this integrated FEM and IPANN model is higher than either the constitutive or IPANN model.

Three-dimensional Finite Element modelling of laser shock peening process

Laser shock peening (LSP) is an innovative surface treatment technique for metal alloys, with the great improvement of their fatigue, corrosion and wear resistance performance. Finite element method has been widely applied to simulate the LSP to provide the theoretically predictive assessment and optimally parametric design. In the current work, 3-D numerical modelling approaches, combining the explicit dynamic analysis, static equilibrium analysis algorithms and different plasticity models for the high strain rate exceeding 106s-1, are further developed. To verify the proposed methods, 3-D static and dynamic FEA of AA7075-T7351 rods subject to two-sided laser shock peening are performed using the FEA package–ABAQUS. The dynamic and residual stress fields, shock wave propagation and surface deformation of the treated metal from different material modelling approaches have a good agreement.

The effect of clamping sequence on dimensional variability in sheet metal assembly

Finite element (FE) modelling techniques have become a popular tool for exploring welding and clamping sequence dependence in sheet metal assemblies. In the current paper, the dimensional variability associated with different assembly clamping sequences is investigated with a FE contact modelling approach implemented in the commercial code Abaqus. A simplified channel section assembly consisting of a top hat and bottom plate is the case study investigated. Expected variation modes of bow and twist were used to simulate key variability sources in the main structural component under investigation; the top hat of the channel section. It was found that final assembly variability can change considerably depending on clamp sequence selection. It was also found that different clamp sequences can control particular modes of variation better than others, and that there is not one particular clamping sequence that is the best for containing all variation modes. An adaptable assembly process is therefore suggested, where given the shape of input components the best available clamping sequence is selected. Comparison of the performance of the proposed adaptable clamping sequence to traditional fixed clamping sequences shows improvements for the dimensional control of variability in non-rigid components. While introducing such a method in production would require inspection of each component being assembled and investigation of the alternative clamping sequences, given access to fast and detailed dimensional inspection technology such as optical coordinate measuring machines (OCMM's), the approach shows promise for future application.

Computational simulation of chassis design and manufacturing

Australian Universities including Deakin University are entering the Formula One SAE car competition this year which was initialised in USA many years ago and has attracted many universities to participate all over the world. There are many categories to be competed for novelty of the design. To achieve optimised design of the car, computational techniques have widely been used in all aspects of the design and manufacturing of the car for the team at Deakin University. In this work, the design and manufacturing of the chassis was computationally simulated to optimise the structure for stress and natural frequencies. The chassis is designed using solid modelling code IDEAS with the geometry exported to ABAQUS, a FEM software, for optimisation.

Simulation of three-point bending test of titanium foam for biomedical applications

Using Titanium (Ti) foam as an implant material is a new approach for biomedical applications and it is important to understand the mechanical behaviors of this new foam material. In the present study, the bending of the Ti foam has been simulated and compared against recently published data [1]. FE Analysis has been performed by Abaqus software. Stiffness and Yield strength of foams between 50% (cortical bone) to 80% (cancellous bone) porosity range were considered. This study showed that crushable foam material model in Abaqus, which has developed primarily for Aluminum (Al) foam alloys, is also valid for Ti Foam before any crack or damage occurs in the sample.

Development of material models to predict the crashworthiness of tubes

Metallic tubes have been extensively studied for their crashworthiness as they closely resemble automotive crash rails. Recently, the demand to produce light weight yet safer vehicles has led to the need to understand the behaviour of novel materials such as composites, metallic foams and sandwich structures durign a crash. This paper presents a method to predict the crashworthiness of structural components using material modes. The material factors that most affect the crushing response are determined and quantified by developing and validating the crushing of a square tube model in Abaqus. The inputs from the model are used to construct a simple, physically realistic constitutive model and new test methods for predicting the material behaviour at high strain rates using low test speeds. These material models enable a designer to predict the crash behaviour of a structure without the need to perform extensive physical tests, thus reducing the time and cost of development.

The effect of skin passing on the material behavior of metal strip in pure bending and tension

The metal strip used in roll forming has often been preprocessed by (tension or roller) leveling or by skin-pass rolling, and as a consequence, may contain residual stresses. These stresses are not well observed by the tensile test, but could have a significant effect on the bending and springback behavior. With the advent of improved process design techniques for roll forming, including advanced finite element techniques, the need for precise material property data has become important. The major deformation mode of roll forming is that of bending combined with unloading and reverse bending, and hence property data derived from bend tests could be more relevant than that from tensile testing.

This work presents a numerical study on the effect of skin passing on the material behavior of stainless steel strip in pure bending and tension. A two dimensional (2-D) numerical model was developed using Abaqus Explicit to analyze the affect of skin passing on the residual stress profile across a section for various working conditions. The deformed meshes and their final stress fields were then imported as pre-defined fields into Abaqus Standard, and the post-skin passing material behavior in pure bending was determined. The results show that a residual stress profile is introduced into the steel strip during skin passing, and that its shape and stress level depend on the overall thickness reduction as well as the number of rolling passes used in the skin passing process. The material behavior in bending and the amount of springback changed significantly depending on the skin pass condition.

Modelling the side impact of carbon fiber tubes

Metallic tubes have been extensively studied for their crashworthiness as they closely resemble automotive crash rails. Recently, the demand to improve fuel economy and reduce vehicle emissions has led automobile manufacturers to explore the crash properties of light weight materials such as fibre reinforced polymer composites, metallic foams and sandwich structures in order to use them as crash barriers. This paper discusses the response of carbon fibre reinforced polymer (CFRP) tubes and their failure mechanisms during side impact. The energy absorption of CFRP tubes is compared to similar Aluminium tubes. The response of the CFRP tubes during impact was modelled using Abaqus finite element software with a composite fabric material model. The material inputs were given based on standard tension and compression test results and the in-plane damage was defined based on cyclic shear tests. The failure modes and energy absorption observed during the tests were well represented by the finite element model.

The indentation behaviour of carbon fibre composite tubes-experiments & modelling

Metallic tubes have been extensively studied for their crashworthiness as they closely resemble automotive crash rails. Recently, the demand to produce lighter weight, yet safer vehicles has led to the need to understand the crash behaviour of novel materials, such as fibre reinforced polymer composites, metallic foams and sandwich structures. This paper discusses the static indentation response of Carbon Fibre Reinforced Polymer (CFRP) tubes. The side impact on a CFRP tube involves various failure mechanisms. This paper highlights these mechanisms and compares the energy absorption of CFRP tubes with similar Aluminium tubes. The response of the CFRP tubes during bending was modelled using ABAQUS finite element software with a composite fabric material model. The material inputs were given based on standard tension and compression test results and the in-plane damage was defined based on cyclic shear tests. The failure modes and energy absorption observed during the tests were well represented by the finite element model.

Spherulitic networks : from structure to rheological property

A finite element method based on ABAQUS is employed to examine the correlation between the microstructure and the elastic response of planar Cayley treelike fiber networks. It is found that the elastic modulus of the fiber network decreases drastically with the fiber length, following the power law. The power law of elastic modulus G′ vs the correlation length ξ obtained from this simulation has an exponent of −1.71, which is close to the exponent of −1.5 for a single-domain network of agar gels. On the other hand, the experimental results from multidomain networks give rise to a power law index of −0.49. The difference between −1.5 and −0.49 can be attributed to the multidomain structure, which weakens the structure of the overall system and therefore suppresses the increase in G′. In addition, when the aspect ratio of the fiber is smaller than 20, the radius of the fiber cross-section has a great impact on the network elasticity, while, when the aspect ratio is larger than 20, it has almost no effect on the elastic property of the network. The stress distribution in the network is uniform due to the symmetrical network structure. This study provides a general understanding of the correlation between microscopic structure and the macroscopic properties of soft functional materials.

Analytic study on pure bending of metal sheets

In this work, analytical models of pure bending are developed to simulate a particular type of bend test and to determine possible errors arising from approximations used in analyzing experimental data. Analytical models proposed for steels include a theoretical solution of pure bending and a series of finite element models, based on the von Mises yield function, are subjected to different stress and strain conditions. The results show that for steel sheets the difference between measured and calculated results of the moment-curvature behaviour is small and the numerical results from the finite element models indicate that experimental results obtained from the test are acceptable in the range of the pure bending operation. Further for magnesium alloys, which exhibit unsymmetrical yielding, the algorithm of the yield function with a linear isotropic hardening model is implemented by programming a user subroutine in Abaqus for bending simulations of magnesium. The simulations using the proposed user subroutine extract better results than those using the von Mises yield function.

Finite element analysis of an axi-symmetric forward spiral extrusion of Mg-1.75Mn alloy

A potential severe plastic deformation process known as axi-symmetrical forward spiral extrusion (AFSE) has been studied numerically and experimentally. The process is based on the extrusion of cylindrical samples through a die with engraved spiral grooves in a near zero shape change manner. The process was simulated using a three dimensional finite element (FE) model that has been developed using commercial software, ABAQUS. In order to verify the finite element results, hot rolled and annealed samples of the alloy were experimentally processed by AFSE. The required extrusion forces during the process were estimated using the FE model and compared with the experimental values. The reasonable agreement between the FE results and experimental data verified the accuracy of the FE model. The numerical results indicate the linear strain distribution in the AFSE sample is only valid for a core concentric while the strain distribution in the vicinity of the grooves is non axi-symmetric. The FE simulation results from this research allows a better understanding of AFSE kinematics especially near the grooves, the required extrusion force and the resultant induced strain distribution in the sample. To compare the mechanical properties of the Mg-1.75Mn alloy before and after the process, a micro shear punch test was used. The tests were performed on samples undergoing one and four passes of AFSE. After four passes of AFSE, it was observed that the average shear strength of the alloy has improved by about 21%. The developedfinite element model enables tool design and material flow simulation during the process.

A study of plastic deformation during axisymmetric forward spiral extrusion and its subsequent mechanical property changes

Axisymmetric forward spiral extrusion (AFSE) accumulates large strains in its sample while extruding it through a die with engraved spiral grooves. A three-dimensional finite element model of AFSE has been developed using ABAQUS to investigate the deformation mode in detail, including the effect of groove geometry and the heterogeneity of plastic deformation. The numerical results demonstrated that the strain distribution in the AFSE sample cross section is linear in the radial direction within a concentric core while the distribution, outside the core, in the vicinity of the grooves is non-linear and non-axisymmetric. Mechanical properties and grain structure changes of the deformed sample were investigated. Improvements of mechanical properties in the processed samples can be attributed to the domination of the shear deformation mode in a plane normal to the extrusion axis and consequently the elongation of grains in the tangential direction