Utilização de planejamento em rede Simplex na previsão de propriedades mecânicas de peças cerâmicas vermelhas

The use of Simplex design to predict the mechanical properties of brick ceramic

(The use of Simplex design to predict the mechanical properties of brick ceramic)

J. Alexandre, R. Toledo, F. Saboya, M. Pereira

UENF Centro de Ciência e Tecnologia ¾ CCT

Universidade Estadual do Norte Fluminense

Av. Alberto Lamego, 2000, Horto, Campos dos Goytacazes, RJ, Brasil, 28015-620

[e-mail]

Resumo

As indústrias cerâmicas de um modo geral, no seu processo produtivo, utilizam matérias-primas provenientes de sedimentos argilosos e/ou alterações de minerais formadores de rochas. Trata-se portanto de um material com alto grau de heterogeneidade tanto do ponto de vista mineralógico quanto em relação ao grau de alteração dos seus minerais, dificultando sobremaneira a previsão do comportamento mecânico do produto final (tijolos, telhas, etc.). Uma prática que é comumente adotada durante a confecção de peças cerâmicas é a mistura de diferentes tipos de solo, realizada muitas vezes, de maneira intuitiva e empírica, onde as características do produto são fortemente dependentes das proporções usadas nessas misturas. O presente trabalho apresenta uma metodologia numérica, utilizando o planejamento em rede Simplex, na previsão das propriedades resultantes destas misturas. Dependendo do modelo escolhido (linear, quadrático, cúbico, etc.), o resultado de uma mistura genérica de até três componentes pode ser determinado numericamente a partir do conhecimento prévio das propriedades de cada componente original, bem como das misturas em determinadas proporções fundamentais. Estes resultados foram testados empiricamente em escala laboratorial e podem ser adaptados posteriormente ao processo produtivo, possibilitando a economia de matéria-prima e reduzindo a exploração predatória das jazidas.

Palavras-chave: misturas, rede Simplex, cerâmica vermelha.

Abstract

In general, the brick ceramic industry uses clayey soils from sedimentation basin and/or from weathered rock as raw material. Due to their nature, these soils, most of the time, are extremely heterogeneous at their natural state. This heterogeneity is not only related to their mineralogical composition but also related to their grading and the degree of mineral weathering. Therefore, an important question arises: how to predict a desirable specific behavior of these materials with such a variable characteristic ? It is very well known that ceramic workers apply an empirical procedure, combining different materials with known properties to manufacture these products. Thus, this paper presents a method to help in overcome this problem by applying a numerical routine which is used to predict specific final properties of different soil admixtures. It is also a helpful tool to determine the ideal admixture proportion that will show some required property, starting from the previously known behavior of the original components and some basic admixtures. The obtained results by carrying out laboratory tests in small ceramic bricks are compared to those predicted ones in order to verify the model suitability. This study is performed in a laboratory scale but it can be easily adapted to the industrial production process, enhancing the economy and preserving the environment from the predatory exploration of borrow areas.

Keywords: mixtures, Simplex lattice, clay.

INTRODUÇÃO

As cerâmicas estruturais em geral utilizam como matéria-prima para confecção de seus artefatos, materiais argilosos provenientes principalmente do intemperismo das rochas e de depósitos sedimentares. Após os processos de intemperismo e/ou sedimentação os minerais decompostos ou simplesmente desagregados passam a formar um aglomerado constituindo a jazida de material argiloso, que possui características próprias dependentes de sua constituição e das condições de sua formação. Devido à complexidade de formação dessa matéria-prima, as previsões que podem ser feitas antes do seu processamento tornam-se um grande desafio do ponto de vista mecânico. Prever o comportamento de artefatos cerâmicos baseando-se nas propriedades mineralógicas de seus constituintes pode vir a se constituir em uma importante ferramenta na busca da qualidade de peças cerâmicas vermelhas. Um método iterativo é apresentado neste trabalho, onde o solos puros e suas misturas correspondentes são analisadas a partir de um número mínimo de experimentos através de um modelo numérico denominado planejamento em rede Simplex [1].

METODOLOGIA

A metodologia proposta neste trabalho baseia-se em um modelo numérico que permita programar as misturas a serem utilizadas na confecção de produtos cerâmicos. O método adotado é o de planejamento em rede Simplex, onde modelos são ajustados para representar as respostas sobre uma superfície. Objetiva-se então encontrar um modelo que descreva o comportamento das misturas com a menor margem de erro possível, quando comparadas com os resultados reais dos experimentos, levando-se em conta o número mínimo necessário de experimentos.

O primeiro passo seria postular um ou mais modelos para descrever as propriedades de interesse (tensão de ruptura, porosidade, etc.) em função da composição das misturas. Em seguida, fazer um planejamento experimental, especificando as composições a serem estudadas. Uma avaliação da eficiência dos modelos pode ser feita através de métodos estatísticos e comparação entre os resultados experimentais e os obtidos com o modelo numérico.

EXPERIMENTOS

Para exemplificar a utilização do método, foram usadas três amostras de solo argiloso obtido na planície de inundação do rio Paraíba do Sul na região de Campos dos Goytacazes. As jazidas desta planície são basicamente constituídas de materiais argilosos cujo principal argilomineral é a caulinita, porém traços de gibsita e ilita são por vezes observados. Essas três amostras foram nomeadas de X₁, X₂ e X₃ e suas distribuições são mostradas no espaço Simplex que corresponde a um triângulo equilátero (Fig. 1). A soma das porcentagens dessas amostras (X= 0, 1/m, 2/m,..., ) deve atender a restrição da equação (A)

A distribuição simétrica das misturas dentro desse espaço é conhecida como rede (lattice), onde uma mistura pode ser um elemento puro ou interações binárias ou ternárias entre os mesmos, como mostrado na Fig. 2. Uma equação polinomial [2] de grau m , equação (B), pode ser usada para resolver o problema; isto é, para fazer uma análise de regressão que atenda à restrição da Equação (A) e passe a ter o número de termos igual ao número de pontos da rede Simplex (equação (C)). A partir da distribuição dos pontos na rede Simplex, atendendo à restrição da equação (A), pode-se através da equação polinomial canônica (equação (B)), encontrar as respostas desejadas (propriedades) em função das misturas.

Na Equação B, polinomial de grau m, h representa as respostas obtidas (observadas), b, os parâmetros do modelo e x as proporções dos componentes na mistura.

O número de pontos em uma rede Simplex {q,m} onde q é o número de componentes e m é o grau do polinômio é dado por O número de pontos em uma rede Simplex {q, m} para 3 £ q £ 10, 1 £ m £ 4 onde o número de níveis para cada componente é m+1 é mostrado na Tabela I.

Apenas a título de exemplificar aplicações do método, foi feito um estudo da variável tensão de ruptura a partir de três amostras (X₁, X₂, e X₃) queimadas à temperatura de 1200 ^ºC. Essas amostras foram moldadas e prensadas nas proporções estratégicas para que o grau dos modelos pudesse ser ajustado em função não apenas da aderência das respostas, mas também em função dos resultados de análise de erro. As interações entre as proporções amostradas são apresentadas na Fig. 3, onde o número de amostras para cada componente foi ensaiado em triplicata para que inferências estatísticas pudessem ser feitas ao modelo utilizando-se para isso dos critérios t de Student e Fisher. As análises estatísticas não foram apresentadas neste trabalho, por necessitarem de uma abordagem mais específica; mas podem ser vistas em detalhes, assim como outros estudos para outras propriedades e temperaturas em [2].

PREPARAÇÃO DAS AMOSTRAS

As amostras foram moldadas em forma de barras (Fig. 4) por prensagem, sob uma pressão de 25.0 MPa, com umidade de moldagem de 6,0%. Após moldagem, as amostras foram secas em estufa e posteriormente queimadas em forno elétrico-resistivo a uma temperatura de 1200 ^ºC sob taxa de aquecimento da ordem de 6 ^ºC/min e patamar de 3 horas. Após esta etapa, as amostras foram submetidas a ensaios de absorção de água, massa específica aparente, porosidade, variação linear e tensão de ruptura à flexão, onde neste trabalho só serão apresentados os resultados de tensão de ruptura à flexão, pois o principal objetivo é demonstrar a aplicabilidade do método.

TR - tensão de ruptura à flexão (MPa)

P - carga atingida no momento da ruptura (kN)

L ¾ distância entre os apoios do corpo de prova (m)

b ¾ largura do corpo de prova (m)

h - espessura do corpo de prova (m).

DETERMINAÇÃO DOS MODELOS EM FUNÇÃO DAS RESPOSTAS OBTIDAS

A partir dos resultados obtidos no experimento, uma avaliação estatística pode ser feita através da análise de variância (ANOVA) onde se determina um valor de R² que expressa o ajuste do modelo da fração descrita pela regressão, dada por:

Onde SQ_R e SQ_T são a soma quadrática dos resíduos e dos erros totais, respectivamente.

O primeiro modelo postulado foi um modelo linear (Fig. 5a) onde a distribuição da rede Simplex é formada por três pontos localizados nos vértices do triângulo.

Após a avaliação, foi observado que os resultados não foram satisfatórios, em virtude do valor de R² ser da ordem de 0,62, havendo a necessidade de aumentar o grau do modelo para um modelo quadrático (Fig. 5b). Este modelo é formado pelas três misturas puras que são os vértices do triângulo, mais três pontos sobre as arestas que são misturas binárias. Nestas, cada componente participa com 50% das proporções (1/2 de X₁, 1/2 de X₂ e 1/2 de X₃). Observando-se que o modelo quadrático, apresenta uma melhoria considerável no valor de R² que foi da ordem de 0,93, tentou-se um modelo mais complexo objetivando uma melhor aderência entre valores medidos e calculados. Considerou-se então um modelo cúbico especial (Fig. 5c), sendo composto por sete pontos, os seis do modelo de grau mais baixo e o sétimo formado por uma mistura ternária no centro do triângulo composta por 1/3 de X₁, 1/3 de X₂e 1/3 de X_3. Obteve-se então um R² de 0,94, que não se traduziu em uma melhoria significativa no resultado final quando comparado com o modelo quadrático.

O resultado obtido com os três modelos e a comparação feita entre estes e os resultados dos experimentos são mostrados na Tabela II. A Fig. 6 mostra a superfície de resposta para o modelo cúbico especial.

Os cálculos foram realizados com as seguintes equações

ANÁLISE DOS RESULTADOS E CONCLUSÃO

Os resultados obtidos com os modelos, quando comparados com os resultados dos ensaios em amostras em função da tensão de ruptura à flexão, mostram que os modelos são significativamente eficazes para previsão dos resultados dessa propriedade, embora o modelo linear, mais simples, mostrou pouca aderência entre resultados previstos e observados. Os resíduos, diferença entre os valores e , de uma maneira geral são aceitáveis, visto que a matéria-prima em estudo é de origem natural e heterogênea, geralmente com alto grau de variabilidade espacial . O grau do modelo necessário para previsão dos resultados, variou do modelo linear ao modelo cúbico especial para o caso das amostras queimadas a temperatura de1200 ^ºC, entretanto graus mais elevados ou mais baixos do modelo são geralmente funções das propriedades analisadas, da temperatura de queima e do grau de heterogeneidade das amostras.

Nota-se que um aumento no grau do modelo, consequentemente da complexidade e do número de ensaios requeridos, não se traduziu em uma melhor resposta para a variável ensaiada, onde o modelo quadrático foi o mais apropriado. Isto se deve ao fato de que erros sistemáticos referentes aos ensaios e a variabilidade das amostras, são do mesmo nível de significância dos coeficientes de grau mais elevados, tornando-os desprezíveis nos resultados previstos pelo modelo.

Outra importante observação a respeito da superfície de respostas, é a possibilidade de se utilizar algoritmos de otimização gráfica para a determinação da mistura adequada englobando todas as variáveis de interesse envolvidas no processo.

AGRADECIMENTOS

Os autores agradecem à Universidade Estadual do Norte Fluminense, ao técnico de laboratório André Luiz Flor Manhães e aos colegas Denise Ribeiro dos Santos e Ronaldo Lerner pelos ensaios executados.

REFERÊNCIAS

(Rec. 08/11/2000, Rev. 20/07/2001, Ac. 17/08/2001)

• [1] J. A. Cornell, "Experiments with mixtures: designs, models and the analysis of mixture data". 2. Ed., New York: John Wiley and Sons (1990).
• [2] J. Alexandre, "Análise de matéria-prima e composições de massa utilizada em cerâmicas vermelhas". Tese de Doutorado ž Campos dos Goytacazes ž RJ, Universidade Estadual do Norte Fluminense - UENF (2000).

Datas de Publicação

• Publicação nesta coleção
  04 Dez 2001
• Data do Fascículo
  Set 2001

Histórico

• Aceito
  17 Ago 2000
• Recebido
  08 Nov 2000
• Revisado
  20 Jul 2001