Summary

Conhecimento Baseado Simulação Nuvem FE dos Processos de formação de folha de metal

Published: December 13, 2016
doi:

Summary

The following paper presents a novel FE simulation technique (KBC-FE), which reduces computational cost by performing simulations on a cloud computing environment, through the application of individual modules. Moreover, it establishes a seamless collaborative network between world leading scientists, enabling the integration of cutting edge knowledge modules into FE simulations.

Abstract

The use of Finite Element (FE) simulation software to adequately predict the outcome of sheet metal forming processes is crucial to enhancing the efficiency and lowering the development time of such processes, whilst reducing costs involved in trial-and-error prototyping. Recent focus on the substitution of steel components with aluminum alloy alternatives in the automotive and aerospace sectors has increased the need to simulate the forming behavior of such alloys for ever more complex component geometries. However these alloys, and in particular their high strength variants, exhibit limited formability at room temperature, and high temperature manufacturing technologies have been developed to form them. Consequently, advanced constitutive models are required to reflect the associated temperature and strain rate effects. Simulating such behavior is computationally very expensive using conventional FE simulation techniques.

This paper presents a novel Knowledge Based Cloud FE (KBC-FE) simulation technique that combines advanced material and friction models with conventional FE simulations in an efficient manner thus enhancing the capability of commercial simulation software packages. The application of these methods is demonstrated through two example case studies, namely: the prediction of a material’s forming limit under hot stamping conditions, and the tool life prediction under multi-cycle loading conditions.

Introduction

Finite Element (FE) simulations have become a powerful tool for optimizing process parameters in the metal forming industry. The reliability of FE simulation results is dependent on the accuracy of the material definition, input in the form of flow stress data or constitutive equations, and the assignment of the boundary conditions, such as the friction coefficient and the heat transfer coefficient. In the past few years, advanced FE simulations have been developed via the implementation of user-defined subroutines, which have significantly broadened the capability of FE software.

The use of such advanced FE simulations in the design of forming processes for structural components has been investigated by both the aviation and automotive industries, with the intention of producing lightweight structures that reduces operating costs and CO2 emissions. Particular focus has been placed on the replacement of steel components with lower density materials, such as aluminum alloys and magnesium alloys. However, these alloys, especially the stronger variants, offer limited formability at room temperature and thus complex-shaped components cannot be manufactured using the conventional cold stamping process. Therefore, advanced high temperature forming technologies, such as warm aluminum forming 1-4, hot stamping of aluminum alloys 5-9 and hot stamping of high strength steels 10, have been developed over the past decades to enable complex-shaped components to be formed. In general, high temperature forming processes involve significant temperature variations, strain rate and loading path changes 11, which would, for instance, cause inevitable viscoplastic and loading history dependent responses from the work piece materials. These are intrinsic features of high temperature forming processes and may be difficult to represent using conventional FE simulation techniques. Another desirable feature would be the ability to predict the tool life over multiple forming cycles in such processes, since they require low friction characteristics achieved through coatings that degrade with each forming operation. To represent all these features via the implementation of user-defined subroutines would be computationally very expensive. Moreover, the development and implementation of multiple subroutines would require excessive multi-disciplinary knowledge from an engineer conducting the simulations.

In the present work, a novel Knowledge Based Cloud FE (KBC-FE) simulation technique is proposed, based on the application of modules on a cloud computing environment, that enables an efficient and effective method of modeling advanced forming features in conjunction with conventional FE simulations. In this technique, data from the FE software is processed at each cloud module, and then imported back into the FE software in the relevant consistent format, for further processing and analysis. The development of these modules and their implementation in the KBC-FE is detailed.

Protocol

1. Desenvolvimento de uma alta temperatura de formação limite modelo de previsão Laser cortar as amostras para testes de formabilidade das chapas de liga de alumínio AA6082 (espessura 1,5 mm) para as geometrias selecionadas 12. Grave um padrão de grelha, composto por pontos circulares 0,75 milímetros de diâmetro com um espaçamento regular de 1 mm, sobre a superfície das amostras utilizando um método electrolítico 13. aplicar manualmente graxa grafite como um lubrificante sobre o lado não-gravadas. Monte o equipamento de teste de cúpula em uma prensa hidráulica elevada taxa de 12. Use uma máquina universal de ensaios hidráulico 250 kN. Aqueça o equipamento de teste de cúpula a uma temperatura de teste e definir o soco em uma velocidade constante movimento. Em seguida, iniciar o teste. Nota: As temperaturas são de teste 300, 400, e 450 ° C, respectivamente. As velocidades de teste incluem 75, 250, e 400 mm / s. Parar o teste na primeira ocorrência de estiramento. Nota: A imprensa Stroke (ou seja, a altura da amostra final) é definida de tal forma que estiramento é apenas observado na amostra formados. Medir a altura da amostra final utilizando um medidor de altura, e calcular as tensões e as taxas de tensão máxima (a taxa de mudança de tensão em relação ao tempo), utilizando um sistema de formação de 3D óptico análise. Analisar as alterações no espaçamento da grade para calcular as estirpes em cada ponto da amostra formados. Certifique-se de que o 3D óptico formando sistema de análise inclui uma câmera, o espécime formados e barras de escala de calibração 14. Nota: A amostra é colocada no centro de uma mesa giratória e fechado com as barras de escala, e as suas posições relativas são mantidas fixas para a duração da análise. Definir a câmara a uma altura fixa (por exemplo, 50 cm) e o ângulo (por exemplo, 30, 50, ou 70 °) para a amostra, e tirar fotografias ao longo de uma rotação completa (360 °) do prato rotativo, em incrementos de 15 ° . Nota: No presetrabalho nt, três conjuntos de imagens foram adquiridas de várias elevações e ângulos de câmera, a fim de mapear as cepas ao longo de toda a amostra 15. Carregar as imagens para o software de análise de formação de 3D óptico, e prossiga para calcular as tensões. Faça isso clicando na função "elipses de computação e pacote ', que detecta os pontos da grelha, seguido clicando a função' pontos 3D de computação e grade 'que edifica a grade. Nota: Calcule as tensões e visualizá-lo no modo de avaliação. Saída as distribuições de deformação para determinar as cepas limite para cada amostra com base na ISO 12004 16, e traçar os diagramas limite de conformação para diferentes velocidades de moldagem e temperaturas de formação. Calibrar um material modelo para AA6082 a diferentes temperaturas de 300 a 500 ° C e taxas de deformação de 0,1 a 10 s -1. Nota: O modelo de material e suas constantes para AA6082estão detalhados na referência 17. Implementar e unificar a função de rendimento anisotrópica Hosford 18, Marciniak-Kuczynski (MK) Teoria 19 e o modelo de material no passo 1,12 em um algoritmo de integração, a fim de formular o modelo de previsão de limite de conformação. Nota: O modelo é descrito na referência 11. Calibrar e verificar o modelo desenvolvido para o passo 1.13 a partir dos resultados experimentais obtidos no passo 1.11. Prever que formam os limites através do modelo verificada 11 a partir do passo 1.14. Nota: A Figura 1 mostra as previsões do modelo resultante a diferentes temperaturas, a uma velocidade de formação a 250 mm / s, ou equivalentemente, uma taxa de deformação de 6,26 s -1. 2. Desenvolvimento de um modelo interativo Friction / desgaste Realizar testes bola-on-disco para espécimes revestidos (disco) Prepare nitreto de titânio (TiN) revestimentos em aço para rolamentosGCr15 disco utilizando arco cátodo e cintilação de magnetrão de frequência média, com os parâmetros de deposição indicado na referência 20. Usando um microscópio eletrônico de varredura (SEM), obter / cross-section superfície topografia da amostra revestida. Medir a espessura do revestimento TiN através das imagens de SEM, comparando a topografia (brilho e contrato) da base e do revestimento materiais. Nota: Os procedimentos experimentais podem ser encontrados na referência 20. Usar uma luz de superfície inter-ferometric pró-filometer branco para obter a rugosidade da superfície da amostra. Colocar a amostra sob a lente e ajustar o microscópio para obter a estrutura de superfície clara. Iluminar a amostra e ajustar os ângulos de eixos X e Y para observar tiras clara interferência (que podem ser monitorados a partir da tela). Definir profundidade bruta no software e começar a medição. verificar automaticamente a superfície da amostra e calcular a rugosidade da superfície. Avaliar a força aderente of a amostra utilizando um aparelho de micro-zero. Aplicar uma carga crescente (máximo de 50 N) e uma distância de zero (no máximo 5 mm) sobre o revestimento TiN. Determinar a falha causando carga crítica do revestimento e obter as curvas de micro-zero 20. Avaliar a dureza da amostra, utilizando um penetrador dureza. Aplicar uma carga estática de 20 N na amostra por 15 s. Medir a diagonal da impressão causada pelo penetrador, em seguida, obter os valores de dureza do testador. Realizar testes bola-de-disco num tribômetro num ambiente ambiente (temperatura 25 ° C, humidade 30%). Usar um WC-6% bola (micro-dureza HV 1780, resistência à abrasão 1,380 N / cm, módulo de elasticidade 71 GPa) 6 mm de diâmetro como o homólogo contra o disco revestido. Ajustar a velocidade de deslizamento em relação a 5 mm / s. Aplicar uma carga normal de 200 N. Inicie os valores do motor e fricção registro usando o tribômetro. Interromper o teste a 180 s, 350 s, 400 s, e 450 s, respectivamente, para analisar a pista de desgaste utilizando um Omicroscópio ptical 20. Medir a topografia da superfície desgastada utilizando um perfilómetro de luz de interferometria superfície branca após o teste. Repetir os testes (Passo 2.1.6) com diferentes cargas normais (300 N, 400 N). Determinar a evolução do coeficiente de atrito até que a desagregação do revestimento duro, caracterizado por um aumento acentuado no coeficiente de atrito Traçar a evolução do coeficiente de atrito contra o tempo depois de gravar os valores de atrito no Passo 2.1.6. Nota: A evolução do coeficiente de atrito é apresentada na referência 20. Avaliar a evolução do coeficiente de atrito em termos de comportamento de desgaste e dos mecanismos associados. Nota: A evolução de atrito é caracterizado em três fases diferentes: (i) fase de baixa fricção, (ii) arar fase de atrito, e (iii) a fase de desagregação de revestimento 20,21. Avaliar o desgaste sTates a 180 s por interromper o teste manualmente, em seguida, analisar a pista de desgaste usando um microscópio óptico. Nota: Este passo é o de investigar os detritos de desgaste para o estágio de baixa fricção, tal como descrito no passo 2.2.2. Repita o passo 2.2.3 em 350 s, 400 s, e 450 s, respectivamente. Desenvolver o modelo de atrito interativo Caracterizar o coeficiente de atrito μ global combinando o atrito inicial μ α com o atrito das partículas de hardware que ara u PC (como mostrado na Eq. (1)) 20. (1) Combine o atrito arar entre a bola eo substrato (μ Ps) com a espessura do revestimento instantâneo (h) para modelar a repartição de revestimento induzida forte aumento do atrito arar μ Pc (Eq. (2)). Nota: Neste caso, é igual a μ μ PC PS quando a espessura do revestimento restante é zero (indicando o colapso completo do revestimento duro). (2) onde λ 1 e λ 2 são parâmetros do modelo introduzido para representar o significado físico do processo de desgaste. λ 1 descreve a influência de grandes partículas de desgaste aprisionadas, e λ 2 representa a intensidade do efeito de fricção arar, que se caracteriza pela inclinação do coeficiente de atrito. Usar um algoritmo de integração de tempo com base para obter a evolução da espessura do revestimento restante e modelar o desgaste acumulado sob condições de contacto variáveis. Atualizar a espessura do revestimento em cada loop de cálculo pela Eq. (3). (3) <img alt="equação 3" src="/files/ftp_upload/53957/53957eq3.jpg"/> em que h é 0 a espessura do revestimento inicial e a taxa de desgaste está dependente do tempo do revestimento. Modificar de Archard lei desgaste 22 (Eq. (4)) e aplicá-lo no presente modelo. (4) em que K é o coeficiente de desgaste, P é a pressão de contacto, v é a velocidade de deslizamento, e H representa a dureza c combinado do revestimento e o substrato. Use o modelo de Korsunsky para calcular a dureza combinada (Eq. (5)). (5) H, onde S é a dureza do substrato, α é a relação entre a dureza do revestimento e o substrato e é o coeficiente de β influência da espessura. Representar a carga de parâmetros dependentes X 1 e K pelo poder lequações aw. (6) (7) onde λ1 κ, κ λ1 K, Ν e Ν K são constantes materiais relacionadas com a evolução de atrito 20. Montar o modelo de atrito interativo para os resultados experimentais utilizando um algoritmo de integração desenvolvido no grupo dos autores para determinar os parâmetros do modelo. 3. Estudos KBC-FE Simulação de casos KBC-FE estudo de caso de simulação 1: previsão de limite de conformação sob condições hot stamping Criar e nomear um novo projeto de simulação no software de simulação FE. Selecione o processo como "formação de estampagem a quente" e do tipo solver como 'PAM-AutoStamp "quandosalvar o projeto. Importar o cunho interior porta clicando no 'Import ferramentas CAD' e depois em 'Import & transfe r' a porta interna 'IGS' arquivo de geometria na interface gráfica do software de simulação FE. Selecione a estratégia de 'Hot formando' para engrenar de ferramentas. Nomeie o objeto importado como 'Die'. Repita o passo 3.1.2 e «importação» os objetos de Punch and Blankholder, respectivamente. Clique em 'branco' no separador 'Set-up'. Clique em "Adicionar em branco" no "editor em branco", e definir o 'New objeto' como 'em branco'. Em seguida, selecione o tipo "superfície em branco '. Escolha 'esboço' para o tipo de definição e importar a forma em branco bclicando y em "Importar do arquivo CAD. Definir "Requinte" como "nível imposto 'e selecione o nível 1 em" opções de malha'. Desligue 'articulada Automatic' e definir «Malhagem» a 4 mm. Definir as propriedades dos materiais em 'editor em branco'. Clique em 'Carregar um material "sob o separador' Material '. Selecione o 'AA6082' (unidade: mm · kg · ms · C) material que as propriedades do material. Defina a "direção rolando 'para' x = 1 '. Defina a 'espessura em branco' para 2 mm, eo 'Temperatura inicial' em branco para 490 ° C. Nota: As propriedades do material e modelo de material estão descritas na referência 17. Clique em "Processo </strong> 'em' 'guia e selecione a opção' Set-up ícone + 'para carregar uma nova macro. Navegue até ' Stamp Hotforming "e selecione" HF_Validation_DoubleAction_GPa.ksa'. Na caixa de diálogo 'Personalizar', ativar o branco, Die, Punch, e Blankholder. Sob a aba "fases", ative gravidade, Segurar, Estampagem e Têmpera. Defina todos os parâmetros nos "Objetos atributos 'no separador' Set-up 'para corresponder com a configuração experimental real (em branco força de retenção = 50 kN, formando speed = 250 mm / s, coeficiente de atrito = 0,1, transferência de calor coeficientes de 23 como um função da abertura e contacto de pressão). Clique em "Verificar" ícone para verificar a simulação set-up e garantir que não haja erros nas definições acima. Clique no ícone 'Computation "para iniciar a siformula-. Nota: O software registra 11 estados durante a simulação em um computador host. Após a conclusão da simulação, observar os resultados da simulação na interface gráfica do software de simulação FE, e prossiga para gravar um 'script' para uma ação exportar os valores de contorno, ou seja, grande pressão (membrana), tensão menor (membrana) e temperatura de todos os elementos em branco, para um estado de simulação especificado. Clique em 'record' e valores de contorno de exportação manualmente. Clique em 'stop' para parar a gravação. Salve o script de modo a repetir a mesma ação para todos os 11 estados de simulação. Clique em 'play' ícone para carregar o script, clique em 'fazer tudo' para exportar os valores de contorno. Nota: Para cada contorno / estado individual, o software exporta automaticamente os valores nos arquivos 'ASCII' em 'major_strain_statenferrugem ',' minor_strain_statenumber ', e' temperature_statenumber ', respectivamente. Salve todos os arquivos exportados para um computador em nuvem. Execute o 'modelo de previsão de carícias "(ou seja, nuvem código do módulo), juntamente com todos os arquivos exportados no computador em nuvem. Prever o aparecimento de estiramento por meio do uso de formar modelo de previsão limite no computador nuvem. Nota: Este modelo 11 oferece aos usuários a opção de executar o modelo de previsão em um elemento individual ou de todos os elementos do branco. Inserir manualmente os detalhes de simulação / parâmetros no "modelo de previsão de carícias". Introduza o número de estados na simulação (estado 11), acidente vascular cerebral total do processo de estampagem (157 mm), estamparia de velocidade (250 mm / s), gama de tensão de interesse (o critério de seleção elemento, por exemplo, tensão> 0,2) e todos os elementos. Nota: O strana faixa limita os elementos para os quais carícias podem ocorrer, definindo um critério elemento, por exemplo, só os elementos com uma grande pressão final superior a 0,2 são selecionados para uma avaliação mais aprofundada no módulo. Depois de completar o cálculo do módulo no computador em nuvem, salvar automaticamente todos os dados (estreitamento resultados de previsão) em arquivos formatados 'ASCII'. Carregar o estado final dos resultados da simulação FE. No separador 'Contornos', clique em 'importado' e depois 'valores escalares'. Selecione o arquivo 'ASCII' obtido a partir do passo acima. Exibir as carícias resultados de previsão no software de simulação FE. KBC-FE estudo de caso de simulação 2: previsão de vida útil sob condições de carregamento multi-ciclo Criar e nomear um novo projeto de simulação no software de simulação FE. Selecione o process como 'Padrão carimbar' eo tipo solver como 'PAM-AutoStamp' ao salvar o projeto. Importar a geometria die clicando no 'Import ferramentas CAD' e depois em 'Import & transferência' do U die arquivo de geometria 'IGS' na interface gráfica do software de simulação FE. Selecione a estratégia de «validação» para engrenar de ferramentas. Nomeie o objeto importado como 'Die'. Repita o Passo 3.2.2 para importar os objetos de Punch and Blankholder, respectivamente. Clique em 'branco' no separador 'Set-up'. 'Adicionar em branco "no" editor em branco', defina o 'New objec t' como 'em branco' e selecione o tipo 'Superfície em branco'. Escolha 'Quatro points 'para o tipo de definição e definir o tamanho em branco para 120 × 80 mm2. Definir "Requinte" como "nível imposto ': Nível 1, sob" opções de malha'. Desligue 'articulada Automatic' e definir «Malhagem» a 1,5 mm. Definir as propriedades dos materiais em 'editor em branco'. Clique no 'Load um material "sob a aba' Material '. Selecione o 'AA5754-H111' (unidade: mm · kg · ms · C) material que as propriedades do material. Defina a "direção rolando 'para' x = 1 '. Defina a 'espessura em branco' para 1,5 mm. Clique em "processo" no separador 'Set-up' e selecione o ícone '+' para carregar uma nova macro. Navegue até' Stamp Viabilidade "e selecione" SingleActioin_GPa.ksa'. Na caixa de diálogo 'Personalizar', ativar o branco, Die, Punch, e Blankholder. Sob "fases", ativar a gravidade, Segurar, e estampagem. Defina todos os 'parâmetros' na simulação para corresponder com a configuração real experiência (forças em branco Prendendo = 5, 20, 50 kN, respectivamente, formando speed = 250 mm / s, coeficiente de atrito = 0,17). "Verificar" a simulação set-up e garantir que não haja erros nas definições acima. Clique no ícone "Computação" e começar a "Computação" para um 11-state U flexão simulação em um computador host. Após a conclusão da simulação, a exportação 'coordenar' de dados e dados de pressão de contacto de automaticamente para a peça de trabalho eferramentas (soco, morrer e titular em branco) como arquivos ASCII '' (como por Passos 3.1.11 e 3.1.12). Salve todos os arquivos exportados para um computador em nuvem. Execute o 'ferramenta módulo de previsão de vida ", juntamente com todos os arquivos exportados no computador em nuvem. Inserir manualmente formando parâmetros no "módulo de previsão de vida útil da ferramenta. Introduza os seguintes parâmetros: número de estados (estado 11), curso total (70 mm), estamparia de velocidade (250 mm / seg) e coeficiente de atrito inicial (0,17). Selecione a ferramenta (soco, morrer, ou titular em branco), e em seguida, iniciar o cálculo para um único elemento ou todos os elementos. Depois de se completar o cálculo do módulo no computador nuvem, gravar automaticamente todos os dados (incluindo a espessura de revestimento instantânea restante e o coeficiente de atrito) em ficheiros ASCII formatados ''. Carregar e exibir a espessura do revestimento restante e frictisobre o coeficiente para os elementos relevantes no software de simulação FE (como por Passo 3.1.17).

Representative Results

Simulação KBC-FE para Estiramento Prediction Em um processo de estampagem a quente, o uso de um espaço em branco otimizado-forma não só irá poupar custos de material, mas também ajudar a reduzir a presença de defeitos, tais como carícias, rachaduras e enrugamento. A forma em branco inicial afecta significativamente o fluxo de material durante a formação e, portanto, um desenho sensível da forma em branco é crítico para o sucesso do processo de estampagem a quente e a qualidade dos produtos finais. Para reduzir os esforços de experiências de tentativa e erro para determinar a geometria em branco ideal, simulação KBC-FE foi provado ser um método altamente eficiente e eficaz para minimizar as áreas com carícias. Usando esta técnica, cada simulação demora cerca de 2 horas, enquanto o paralelo computação módulo de nuvem para carícias previsão é concluída dentro de 4 horas. <p class="jove_content" fo:keep-together.withiN-page = "1"> A Figura 4 mostra a evolução da forma em branco usado na estampagem a quente, um exemplo de componente interno da porta automóvel. A forma em branco inicial, adotada a partir de um processo de estampagem a frio convencional, foi utilizado pela primeira vez na simulação KBC-FE. Os resultados experimentais na Figura 4 (a) mostram que a grande falha (rachaduras ou estiramento) áreas são visíveis após a estampagem a quente. Depois de uma iteração da optimização forma em branco, que pode ser visto na Figura 4 (b), que um painel de quase totalmente bem sucedida é formado com muito menos estiramento, comparado com a utilização da forma inicial em branco. Pode ver-se que existe ainda uma indicação de estiramento nos bolsos na parte superior direita e cantos esquerdo do painel. Após optimização adicional na Figura 4 (c), a forma optimizada em branco foi finalmente obtido sem estreitamento visível no painel. A forma em branco otimizado determinada pela simulação KBC-FE foi verificado experimentalmente através de estampagem a quenteEnsaios realizados sobre uma linha de produção totalmente automatizada oferecido por um fabricante do sistema de produção. Simulação KBC-FE de Previsão de Vida Ferramenta FE simulações convencionais de processos de formação de metal são realizados por um único ciclo. No entanto, num ambiente de produção, ciclos múltiplos de formação são executadas sobre uma dada ferramenta, onde se verificar que um aumento no número de ciclos que formam resulta num aumento da diferença entre os componentes formados. Esta variação durante o carregamento ferramenta multi-ciclo é o resultado de mudar a topografia da superfície. Por exemplo, o carregamento de múltiplos ciclos de formação de ferramentas com revestimentos funcionais irá levar a uma redução de espessura de revestimento devido ao desgaste. Além disso, a desagregação do revestimento, também irá ser influenciada por formação de parâmetros, tais como a carga / pressão, formando velocidades, etc A técnica KBC-FE permite osimulação de folha de metal formando processos sob condições de carga multi-ciclo, o que é essencial para a previsão da vida em serviço de formação de ferramentas com revestimentos funcionais avançados. Para investigar os efeitos da força de retenção em branco na vida útil da ferramenta, em branco valores força de retenção de 5, 20 e 50 kN foram examinados para uma velocidade de formação constante de 250 mm / s. A Figura 5 mostra a distribuição da espessura do revestimento da ferramenta restante com diferentes forças de retenção em branco depois de 300 ciclos de formação. Isto indica claramente que a espessura do revestimento restante diminui com um aumento da força de retenção em branco. A Figura 6 mostra a distribuição da espessura de revestimento e a pressão remanescente com as forças de retenção em branco de 5, 20 e 50 kN, respectivamente, ao longo da distância curvilínea do molde após 300 ciclos de formação. Como a região AB representa o ent dierance região durante o processo de dobragem em forma de U, a pressão e a distância em relação ao desgaste nesta região eram muito maiores do que outras regiões da matriz. Por conseguinte, o desgaste do revestimento ocorreu principalmente nesta área. Existem dois valores de pico de redução da espessura de revestimento a 20 kN e 50 kN, que correspondem aos dois picos, sob a pressão. Enquanto isso, a espessura do revestimento restante diminui com o aumento da força de retenção em branco. Os menores espessuras de revestimento restante com as forças de exploração em branco de 5, 20 e 50 kN, foram 0,905, 0,570 e 0,403 microns, respectivamente, em que a espessura do revestimento inicial foi de 2,1 microns. Figura 1: Comparação entre as estirpes experimentais e previstos, formando-limite em diferentes temperaturas. As estirpes limite de conformação aumentar à medida que a temperatura sobe, a uma velocidade constante de 250 mm/ s, ou equivalentemente, a taxa de deformação de 6,26 s-1. Por favor clique aqui para ver uma versão maior desta figura. Figura 2: Gráfico esquemático para o conhecimento baseado em nuvem FE simulação de um processo de chapa. software comercial de simulação FE, é usado para executar a simulação e exportar os resultados necessários para os módulos individuais. Os módulos, por exemplo, maleabilidade, transferência de calor, força pós-formação (microestrutura), previsão de vida útil da ferramenta, da ferramenta de projeto, etc., trabalham simultaneamente e independentemente na nuvem, portanto, permitindo a integração de conhecimentos de ponta de várias fontes em simulações FE . Por favor click aqui para ver uma versão maior desta figura. Figura 3: Geometria da peça de trabalho e ferramentas para a forma de U flexão simulação. As ferramentas, ou seja, soco, titular em branco e morrem, são modelados usando elementos rígidos. elementos de casca são usadas para elementos da peça de trabalho (em branco). Por favor clique aqui para ver uma versão maior desta figura. Figura 4: Evolução da forma em branco para estampagem a quente de uma porta de painel interno (exibido na simulação FE). Esquerda: Os números em molduras verdes representam formas em branco em cada fase de otimização, e os de vermelhoquadros de correspondência com a forma em branco antes da sua optimização. Direita: estiramento resultados de previsão em cada fase de otimização. (A) Os resultados iniciais com grande fracasso (craqueamento / estiramento mostrado na cor vermelha), (b) Redução do fracasso com algum estiramento após a primeira etapa de otimização, (c) forma em branco otimizado final sem estiramento visível. Por favor clique aqui para ver uma versão maior desta figura. Figura 5: A distribuição da espessura do revestimento restante (exibido na simulação FE) com forças que prende o de: (a) 5 kN, (b) de 20 kN, e (c) 50 kN, após 300 ciclos que formam a uma velocidade de estampagem constante de 250 mm / s. Por favorclique aqui para ver uma versão maior desta figura. Figura 6: Previsão da pressão de contacto e a espessura do revestimento restante com as forças de retenção em branco de: (a) 5 kN, (b) de 20 kN, e (c) 50 kN, ao longo da distância curvilínea da fieira a uma velocidade de estampagem constante de 250 mm / s. Por favor clique aqui para ver uma versão maior desta figura.

Discussion

A técnica de simulação KBC-FE permite simulações avançadas a ser realizado fora do local usando os módulos dedicados. Ele pode ser executado módulos funcionais em um ambiente de nuvem, que apontam-se os nós de diferentes especializações, para garantir que as simulações de processos são conduzidos com a maior precisão possível. Os aspectos críticos na simulação KBC-FE pode envolver a independência dos códigos FE, a eficiência da computação, ea precisão dos módulos funcionais. A realização de cada função avançada em um módulo seria contar com o desenvolvimento de um novo modelo e / ou uma nova técnica experimental. Por exemplo, o módulo de limite de conformação é desenvolvido com base no novo modelo de previsão de limite de conformação unificado 11, eo módulo de previsão de vida útil da ferramenta atrito atualmente tem sido desenvolvido pela implementação do modelo de atrito interativa 20. A técnica de simulação KBC-FE também oferece a função de cálculo selectiva, isto é, apenas os elementos cumprindo a selecçãocritérios são selecionados para uma avaliação mais aprofundada nos módulos individuais. Por exemplo, o módulo de previsão de vida da ferramenta selecciona automaticamente os elementos para os quais o revestimento duro tende a desagregação, pela classificação da taxa de desgaste de todos os elementos na primeira formando ciclo, assim, normalmente menos de 1% dos elementos serão seleccionados para posterior avaliações de vida útil da ferramenta sob condições de carga multi-ciclo. Na presente pesquisa, a previsão de vida útil da ferramenta depois de 300 ciclos de formação pode ser concluído dentro de 5 min.

Ao realizar os testes pertinentes e calibração em conformidade, o modelo de formação de limite poderia ser aplicada para a formação de simulações de processos para determinar consequentemente, os parâmetros ideais para a produção de um componente a partir de tais ligas com sucesso, e sem ocorrência de estiramento. O modelo de previsão limite de conformação foi desenvolvido como um módulo de nuvem que era independente do software FE ser utilizado, e pode ser aplicado a qualquer software FE para avaliar a formabilidade de um material duranteformando, sem sub-rotinas complicadas 17. Ao importar os dados relevantes para o modelo, os cálculos podem ser realizados para determinar se a falha poderia ocorrer, em regiões do componente que o usuário pode especificar, economizando recursos computacionais. No entanto, deve notar-se que, como as curvas de tensão-deformação são inseridos no software FE por meio de uma simples tabela de consulta, que podem ser difíceis de representar totalmente as propriedades do material a várias temperaturas e taxas de deformação durante a simulação.

No módulo de previsão de vida útil da ferramenta, o comportamento de atrito durante a formação pode ser previsto através da importação dos dados do histórico de deformação necessários para o módulo de atrito verificados 20, e em seguida, importando os pontos de dados discretos calculados pelo módulo de nuvem para cada elemento de volta para o software FE. Isso garante que o módulo de atrito avançado pode ser utilizado por todos os códigos de FE, independentemente da sua capacidade para incorporar utilizadores-sub-rotinas. Além disso, a modificaçãoule poderia ser executado em paralelo para reduzir ainda mais o tempo de computação. O modelo de atrito / desgaste interativo assumiu a ausência de partículas de desgaste durante a inicial de correr, e, como resultado, seria razoável esperar um valor inicial constante de coeficiente de atrito de 0,17 20. Embora este modelo revelou a evolução da distribuição de atrito, o comportamento de atrito durante um processo de formação é muito complicado e é difícil de integrar completamente o comportamento de atrito complexo a partir do módulo de nuvem para a simulação FE.

Enquanto tecnologia do futuro, a simulação KBC-FE contará com o desenvolvimento de pacotes de software de simulação FE baseados na Internet dedicado e robusto, o que exigiria um modelo de negócio altamente rentável, mas completamente diferente a ser estabelecida pelos desenvolvedores de software. Além disso, uma rede interna dedicada precisa ser construída dentro dos partidos de colaboração para garantir a segurança dos dados ea confiabilidade do sistema de controle industrial. </p>

Disclosures

The authors have nothing to disclose.

Acknowledgements

The financial support from Innovate UK, Ultra-light Car Bodies (UlCab, reference 101568) and Make it lighter, with less (LightBlank, reference 131818) are gratefully acknowledged. The research leading to these results has received funding from the European Union’s Seventh Framework Program (FP7/2007-2013) under grant agreement No. 604240, project title ‘An industrial system enabling the use of a patented, lab-proven materials processing technology for Low Cost forming of Lightweight structures for transportation industries (LoCoLite)’. Significant support was also received from the AVIC Centre for Structural Design and Manufacture at Imperial College London, which is funded by Aviation Industry Corporation of China (AVIC).

Materials

AA6082-T6 AMAG Material
AA5754-H111 AMAG Material
 1000 kN high-speed press ESH Forming press
ARGUS GOM Optical forming analysis
PAM-STAMP 2015 ESI FE simulation software
Matlab  MathWorks Numerical calculation software
Gleeble 3800 DSI Uniaxial tensile test
High Temperature Tribometer (THT) Anton Paar Friction property test
NewViewTM 7100 ZYGO Surface profilometer
 Magnetron sputtering equipment Coating deposition
Microhardness tester Wolpert Wilson Instruments
Nano-hardness indenter  MTS

References

  1. Miller, W. S., et al. Recent development in aluminium alloys for the automotive industry. Mater. Sci. Eng. A. 280 (1), 37-49 (2000).
  2. Bolt, P. J., Lamboo, N. A. P. M., Rozier, P. J. C. M. Feasibility of warm drawing of aluminium products. J. Mater. Process. Tech. 115 (1), 118-121 (2001).
  3. Li, D., Ghosh, A., et al. Effects of temperature and blank holding force on biaxial forming behavior of aluminum sheet alloys. J. Mater. Eng. Perform. 13 (3), 348-360 (2004).
  4. Toros, S., Ozturk, F., Kacar, I. Review of warm forming of aluminum-magnesium alloys. J. Mater. Process. Tech. 207 (1-3), 1-12 (2008).
  5. Wang, L., Strangwood, M., Balint, D., Lin, J., Dean, T. A. Formability and failure mechanisms of AA2024 under hot forming conditions. Mater. Sci. Eng. A. 528 (6), 2648-2656 (2011).
  6. Wang, L., et al. . TTP2013 Tools and Technologies for Processing Ultra High Strength Materials. , (2013).
  7. El Fakir, O., et al. Numerical study of the solution heat treatment, forming, and in-die quenching (HFQ) process on AA5754. Int. J. Mach. Tool. Manu. 87, 39-48 (2014).
  8. Raugei, M., El Fakir, O., Wang, L., Lin, J., Morrey, D. Life cycle assessment of the potential environmental benefits of a novel hot forming process in automotive manufacturing. J. Clean. Prod. 83, 80-86 (2014).
  9. Liu, J., Gao, H., Fakir, O. E., Wang, L., Lin, J. HFQ forming of AA6082 tailor welded blanks. MATEC Web of Conferences. 21 (05006), (2015).
  10. Karbasian, H., Tekkaya, A. E. A review on hot stamping. J. Mater. Process. Tech. 210 (15), 2103-2118 (2010).
  11. El Fakir, O., Wang, L., Balint, D., Dear, J. P., Lin, J. Predicting Effect of Temperature Strain Rate and Strain Path Changes on Forming Limit of Lightweight Sheet Metal Alloys. Procedia Eng. 81, 736-741 (2014).
  12. Shi, Z., et al. . the 3rd International Conference on New Forming Technology. , 100-104 (2012).
  13. Mohamed, M. S., Foster, A. D., Lin, J., Balint, D. S., Dean, T. A. Investigation of deformation and failure features in hot stamping of AA6082: Experimentation and modelling. Int. J. Mach. Tool. Manu. 53 (1), 27-38 (2012).
  14. Hosford, W. F. Comments on anisotropic yield criteria. Int. J. Mech. Sci. 27 (7), 423-427 (1985).
  15. Marciniak, Z., Kuczyński, K. Limit strains in the processes of stretch-forming sheet metal. Int. J. Mech. Sci. 9 (9), 609-620 (1967).
  16. Ma, G., Wang, L., Gao, H., Zhang, J., Reddyhoff, T. The friction coefficient evolution of a TiN coated contact during sliding wear. Appl. Surf. Sci. 345, 109-115 (2015).
  17. Põdra, P., Andersson, S. Simulating sliding wear with finite element method. Tribol. Int. 32 (2), 71-81 (1999).
  18. Archard, J. F. Contact and Rubbing of Flat Surfaces. J. Appl. Phys. 24 (8), 981-988 (1953).
  19. Liu, X., et al. Determination of the interfacial heat transfer coefficient in the hot stamping of AA7075. MATEC Web of Conferences. 21 (05003), (2015).

Play Video

Cite This Article
Zhou, D., Yuan, X., Gao, H., Wang, A., Liu, J., El Fakir, O., Politis, D. J., Wang, L., Lin, J. Knowledge Based Cloud FE Simulation of Sheet Metal Forming Processes. J. Vis. Exp. (118), e53957, doi:10.3791/53957 (2016).

View Video