Optimization of the heating system by electrical resistances in a rapid thermal response mold based on MSR-PSO-FEM

==Resumen==

<span id='OLE_LINK1'></span><span id='OLE_LINK2'></span>
En el presente trabajo se presenta la optimización del sistema de calentamiento por resistencias eléctricas para un molde de respuesta térmica rápida (MRTR), usando la metodología de superficie de respuesta (MSR) así como una técnica de optimización. Aplicando la técnica de diseño de experimentos Box-Behnken, se diseñó una matriz de experimentos con cuatro factores y tres niveles. Las variables de diseño que se emplean para describir el diseño y la forma del sistema de calentamiento, son el flujo de calor por unidad de área de la resistencia eléctrica, las distancias desde el centro de los calentadores a la superficie de la cavidad, la distancia entre las resistencias eléctricas adyacentes y el espesor de la placa de calentamiento. El tiempo de calentamiento, la variación de la temperatura en la cavidad y los esfuerzos de Von-Mises fueron considerados como las variables del modelo. Se realizan análisis térmicos y de resistencia térmico-estructural del modelo basado en método de elemento finito (FEM por sus siglas en inglés) para adquirir las variables objetivas. Además, se desarrollan modelos matemáticos de superficie de respuesta mediante el modelo de regresión mixta y el modelo de superficie de respuesta y se emplea el método de análisis de varianza (ANOVA por sus siglas en inglés), para comprobar la exactitud de dichos modelos matemáticos. Con los modelos obtenidos, se optimiza la posición de las resistencias eléctricas y se reduce la relación entre masa y el volumen del inserto de la cavidad para minimizar el tiempo de calentamiento dentro de una distribución de temperatura y resistencia estructural razonables, acoplando los modelos de superficie de respuesta desarrollados con el método de optimización de enjambre de partículas (PSO por sus siglas en inglés). Los resultados obtenidos muestran que el tiempo de calentamiento requerido en la superficie de la cavidad se puede reducir significativamente en el ciclo de moldeo, demostrando con estos la efectividad del sistema de calentamiento.

'''Palabras Clave''': Molde de Respuesta Térmica Rápida (MRTR), moldeo por inyección, Optimización de Enjambre de Partículas (PSO), análisis térmico

==Abstract==

The aim of this work to optimization of the heating system by electrical resistances for a rapid thermal response mold (MRTR), using the response surface methodology (MSR). Applying the design technique of the Box-Behnken experiments, a matrix of experiments with four factors and three levels was designed. The design variables that are used to describe the design and shape of the heating system are the heat flux per unit area of the electrical resistance, the distances from the center of the heaters to the surface of the cavity, the distance between the adjacent electrical resistors and the thickness of the heating plate. The heating time, the variation of the temperature in the cavity and the Von-Mises stress were considered as the variables of the model. Thermal and thermal-structural resistance analyzes of the model based on finite element method (FEM) are performed to acquire the objective variables. Mathematical response surface models are developed using the mixed regression model and the response surface model and the variance analysis method (ANOVA) is used to verify the accuracy of these mathematical models. With the obtained models, the position of the electric resistances is optimized and the ratio between mass and volume of the cavity insert is reduced to minimize the heating time within a reasonable temperature distribution and structural strength, coupling the surface models of response developed with the method of particle swarm optimization (PSO). The results obtained indicate that the required heating time on the surface of the cavity can be significantly reduced in the molding cycle, demonstrating with these the effectiveness of the heating system.

'''Keywords''': Rapid thermal response molding (RTRM), injection molding, Particle Swarm Optimization (PSO), thermal analysis

==1. Introducción==

Recientemente, la industria de los polímeros está creciendo de manera sorprendente, sobre todo en las últimas décadas, por lo que es necesario estudiar nuevas técnicas que mejoren los procesos de manufactura de partes plásticas, como inyección y extrusión. Las mejoras al proceso se han visto en tres principales áreas, las máquinas de inyección, los moldes y los polímeros. La temperatura del molde es un parámetro de procesamiento importante en el tiempo de ciclo, además, de afectar significativamente la calidad superficial de la pieza moldeada. Si la temperatura del molde es baja en el proceso de moldeado, se acorta el tiempo de ciclo, sin embargo, esto provocará que polímero fundido solidifique prematuramente, por lo tanto, da lugar a una calidad superficial pobre en la superficie de la pieza, rugosidad, líneas de soldadura, rechupes, entre otros defectos [1,2]. Para minimizar estos problemas desde el punto de vista de proceso, se requiere aumentar la presión y/o la velocidad de inyección [3]. Si la parte de plástico requiere de un alto aspecto de brillo, los defectos superficiales requieren ser eliminados mediante operaciones secundarias como el lijado y pulido. Estas operaciones aumentan el costo de producción, el tiempo total del proceso de producción y requieren energía consumida en procesos secundarios. Aumentar la temperatura en la superficie de la cavidad del molde hasta la temperatura de cristalización (<math display="inline">{T}_{c}</math>) mejora la calidad de la superficie de la pieza [4], pero esto extenderá el tiempo de ciclo en el proceso de moldeo regular.

Por esta razón, recientemente se ha desarrollado una nueva tecnología de moldeo llamada MRTR para satisfacer las exigencias cada vez más estrictas de las piezas de plástico. En este proceso de moldeo, la temperatura de la cavidad del molde se cambia alternativamente mediante un rápido calentamiento y enfriamiento en cada ciclo de moldeo. Debido a la elevada temperatura de la cavidad del molde, el polímero fundido llena fácilmente la superficie de la cavidad del molde y por lo tanto pueden moldearse piezas de alto brillo. Para variar la temperatura de la cavidad del molde sin un aumento significativo en el tiempo de ciclo, se requiere una técnica rápida de calentamiento del molde. En los últimos años, se han propuesto algunos métodos rápidos de calentamiento del molde. Estos métodos se pueden dividir principalmente en dos categorías: calefacción externa y calefacción interna. La primera categoría incluye principalmente calentamiento por inducción [5], calentamiento asistido por gas [6], calentamiento con alta frecuencia [7], calentamiento por infrarrojos [8], entre otros métodos de calentamiento. En cuanto a la segunda categoría, se tiene el calentamiento asistida por vapor [9], calentamiento por fluido [10-12] y calentamiento por resistencias eléctricas [13-15] son los métodos más utilizados.

Los métodos de calentamiento externos tienen una buena eficiencia en el calentamiento, así como un costo elevado, además la estructura compleja de los moldes limita sus aplicaciones en la producción real. En aplicaciones prácticas, los métodos de calentamiento interno son preferidos en términos de su mejor adaptabilidad y estabilidad, particularmente para el molde con un volumen relativamente grande. Algunos canales se perforan dentro del molde para transportar el medio caliente o instalar los elementos de calentamiento eléctrico.Los métodos de calentamiento por fluido [10-12] han mostrado una eficiencia mayor en el calentamiento; sin embargo, la infraestructura y costo para su funcionamiento es mayor que la que se requiere en un sistema de calentamiento por resistencias eléctricas, éstas presentan además de las ventajas mencionadas, el mantenimiento casi nulo y el reemplazo de cualquier componente en caso de avería es mínimo y económico. Para obtener el mejor rendimiento del sistema de calentamiento en términos de la distribución de la temperatura en la superficie de la cavidad del molde y la eficiencia de calentamiento, es necesario optimizar la posición de los elementos de calentamiento en el molde. Para solucionar este problema, se han propuesto algunos métodos de diseño de optimización. Li et al. [16], propusieron la estrategia de optimización combinando la metodología de la superficie de respuesta (MSR) con el algoritmo genético (AG) para optimizar el diseño de los canales de calentamiento para un molde de inyección de calentamiento rápido con calentamiento asistido por vapor. Una red neuronal artificial (RNA) [17], suele construir para las variables definidas de un proceso o sistema ajustando una ecuación polinomial, muy semejante a la obtenida por la MSR y estar acoplada a un AG para su optimización [18]. Sin embargo, se ha demostrado que el método de búsqueda estocástico para la optimización, llamado PSO resulta ser más rápido y mejor que las redes neuronales y los algoritmos genéticos [19-22].

En este artículo se presenta la optimización del sistema de calentamiento por resistencias eléctricas para un molde de respuesta térmica rápida (MRTR), usando la metodología de superficie de respuesta (MSR) así como una técnica de optimización. Aplicando la técnica de diseño de experimentos Box-Behnken, se diseñó una matriz de experimentos con cuatro factores y tres niveles. Las variables de diseño que se emplean para describir el diseño y la forma del sistema de calentamiento, son el flujo de calor por unidad de área de la resistencia eléctrica, las distancias desde el centro de los calentadores a la superficie de la cavidad, la distancia entre las resistencias eléctricas adyacentes y el espesor de la placa de calentamiento. El tiempo de calentamiento, la variación de la temperatura en la cavidad y los esfuerzos de Von-Mises fueron considerados como las variables del modelo. Se realizan análisis térmicos y de resistencia térmico-estructural del modelo basado en método de elemento finito (MEF) para adquirir las variables objetivas. Además, se desarrollan modelos matemáticos de superficie de respuesta mediante el modelo de regresión mixta y el modelo de superficie de respuesta y se emplea el método de análisis de varianza (ANOVA), para comprobar la exactitud de dichos modelos matemáticos. Con los modelos obtenidos, se optimiza la posición de las resistencias eléctricas y se reduce la relación entre masa y el volumen del inserto de la cavidad para minimizar el tiempo de calentamiento dentro de una distribución de temperatura y resistencia estructural razonables, acoplando los modelos de superficie de respuesta desarrollados con el método de optimización de enjambre de partículas (PSO, por sus siglas en inglés). Los resultados obtenidos muestran que el tiempo de calentamiento requerido en la superficie de la cavidad se puede reducir significativamente en el ciclo de moldeo, demostrando con estos la efectividad del sistema de calentamiento.

==2. Modelo físico y matemático del molde de respuesta térmica rápida (MRTR)==

En la presente sección se presenta la estructura física del molde a evaluar, se muestra la cavidad del molde, las resistencias eléctricas, los canales de enfriamiento, los materiales utilizados y las condiciones de frontera, además, los fundamentos matemáticos necesarios para la evaluación de la respuesta térmica en el molde de inyección.

La principal diferencia entre el moldeo por inyección convencional (MIC) y el moldeo con un MRTR es la estrategia de control de la temperatura del molde. En el MRTR, el molde debe calentarse y enfriarse rápidamente para lograr un moldeo por inyección a alta temperatura sin aumentar demasiado el tiempo de ciclo de moldeo. La Figura 1, muestra de manera esquemática el MRTR que se diseña y propone en este trabajo, en la Figura 1(a), se observan los componentes principales del molde, canales de enfriamiento, resistencias eléctricas, placa retenedora, inserto, entre otros, el molde propuesto está bajo dimensiones estandarizadas de la empresa HASCO®, de <math>196 mm \times 296 mm</math>. La Figura 1 (b), se especifica las dimensiones de la pieza final, así como la cavidad del MRTR. Para este trabajo se propuso una probeta de tensión en base a la norma ASTM D638-03 [23], tipo I. La Figura 1(c), muestra un plano de corte trasversal del MRTR, el corte se realizó en el centro de la probeta de tensión.

Un control dinámico de la temperatura del molde se implementa calentando y enfriando el molde en cada ciclo del proceso, por lo que la transferencia de calor es muy diferente a la del proceso MIC. En la etapa de calentamiento del molde, este debe ser calentado hasta la temperatura preestablecida utilizando resistencias eléctricas. Wang propone que, para una mejor apariencia de la pieza plastica, la temperatura en la superficie de la cavidad del molde debe ser en promedio la temperatura de cristalización del polímero en la etapa del llenado [4]. El calor generado se transfiere principalmente al acero del molde a través de la conducción del calor, solamente una pequeña cantidad de calor se acumula en las resistencias eléctricas así como el marco del molde y el resto se disipa al ambiente.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image1.png|center|382px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 1'''. Esquema del MRTR para una probeta de tensión. (a) Isométrico del MRTR. (b)  Probeta de tensión, unidades en mm, (c) Corte de sección transversal del MRTR
|}

La Figura 2(a), ilustra esquemáticamente la transferencia de calor en la etapa de calentamiento del MRTR con calentamiento eléctrico. Basándose en la ley de conservación de energía, la ecuación de equilibrio de transferencia de calor del molde durante la etapa de calentamiento puede expresarse:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{Q}_{e}={Q}_{m}+{Q}_{c}+{Q}_{a}^{c}+{Q}_{i}^{c}</math>
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(1)
|}

donde <math display="inline">{Q}_{e}</math> es el calor generado por las resistencias eléctricas; <math display="inline">{Q}_{m}</math> es el calor transferido al acero del molde; <math display="inline">{Q}_{c}</math> indica el calor acumulado en las resistencias; <math display="inline">{Q}_{a}^{c}</math> y <math display="inline">{Q}_{i}^{c}</math> representan el calor disipado al ambiente y al marco del molde en la etapa de calentamiento del molde, respectivamente.

Los <math display="inline">{Q}_{c}</math>, <math display="inline">{Q}_{a}^{c}</math> y <math display="inline">{Q}_{i}^{c}</math> son bastante pequeños en comparación con el <math display="inline">{Q}_{m}</math> y por lo tanto tienen efectos despreciables sobre el campo de temperatura del molde durante el calentamiento. Así, la ec. (1) se puede simplificar como [24,25]:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{Q}_{e}\approx {Q}_{m}</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(2)
|}

donde <math display="inline">{Q}_{e}</math> y <math display="inline">{Q}_{m}</math> puede ser representado como:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{Q}_{e}={P}_{r}{t}_{cr}</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(3)
|}


{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{Q}_{m}={C}_{m}{\rho }_{m}{V}_{m}(\overline{{T}_{m}}-\overline{\overline{{T}_{m}^{i}}})</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(4)
|}

donde <math display="inline">{P}_{r}</math> es la potencia total de las resistencias eléctricas; <math display="inline">{t}_{cr}</math> es el tiempo de calentamiento requerido; <math display="inline">\overline{{T}_{m}^{i}}</math> es la temperatura inicial del molde;  <math display="inline">\overline{{T}_{m}}</math> es la temperatura media del molde en el tiempo <math display="inline">{t}_{cr}</math>; <math display="inline">{C}_{m}</math>, <math display="inline">{\rho }_{m}</math> y <math display="inline">{V}_{m}</math> representan el calor especifico, densidad y volumen del molde, respectivamente.

Sustituyendo las ecs. (3) y (4) en (2) se obtiene la siguiente expresión:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{t}_{cr}={[C}_{m}{\rho }_{m}{V}_{m}(\overline{{T}_{m}}-\overline{{T}_{m}^{i}})]/{P}_{r}</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(5)
|}

Cuando se moldea una pieza de plástico específica usando un MRTR, las  <math display="inline">\overline{{T}_{m}}</math> y <math display="inline">\overline{{T}_{m}^{i}}</math> están generalmente predeterminadas y por lo tanto pueden considerarse constantes. De aquí que, <math display="inline">{t}_{cr}</math> depende principalmente de los parámetros <math display="inline">{C}_{m}</math>, <math display="inline">{\rho }_{m}</math>, <math display="inline">{V}_{m}</math> y <math display="inline">{P}_{r}</math>. La eficiencia de calentamiento del molde se puede mejorar, es decir, el <math display="inline">{t}_{cr}</math> en la ec. (5) se puede minimizar reduciendo los valores de <math display="inline">{C}_{m}</math>, <math display="inline">{\rho }_{m}</math>, <math display="inline">{V}_{m}</math> y/o aumentando el valor del <math display="inline">{P}_{r}</math>.

El producto de <math display="inline">{C}_{m}</math> y <math display="inline">{\rho }_{m}</math> para la mayoría de los materiales del molde tiene poca diferencia, por lo que la forma de seleccionar el material del molde con menor <math display="inline">{C}_{m}</math> y <math display="inline">{\rho }_{m}</math> para mejorar la eficiencia de calefacción del molde no es la adecuada [26]. Por lo tanto, bajo la condición de resistencia y rigidez estructural del molde adecuadas, el volumen del molde calentado debe reducirse tanto como sea posible para reducir la capacidad calorífica del molde. Además, aumentar el número y/o potencia de las resistencias eléctricas adecuadamente puede mejorar aún más la eficiencia de calentamiento del molde.

En la Figura 2(b), se muestra la etapa de enfriamiento del molde, en esta, se bombea refrigerante a los canales de enfriamiento del molde para enfriar rápidamente el molde y la pieza moldeada. El calor acumulado en el molde de acero y las resistencias eléctricas en la etapa de calentamiento, así como el calor liberado de la masa fundida de polímero, durante el enfriamiento son eliminados por el refrigerante en circulación; además, una pequeña cantidad de calor se disipa al ambiente y al marco del molde.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image2.png|center|358px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 2'''. Esquemas de transferencia de calor del MRTR con calentamiento eléctrico. (a) Calentamiento. (b) Enfriamiento
|}


La ecuación de equilibrio de transferencia de calor del molde durante el enfriamiento puede expresarse como:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{Q}_{p}+{Q}_{m}+{Q}_{c}={Q}_{r}+{Q}_{a}^{e}+{Q}_{i}^{e}</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(6)
|}

donde <math display="inline">{Q}_{p}</math> es el calor liberado del polímero fundido; <math display="inline">{Q}_{r}</math> representa el calor removido por el refrigerante; <math display="inline">{Q}_{a}^{e}</math> y <math display="inline">{Q}_{i}^{e}</math> representan el calor disipado al ambiente y al marco del molde en la etapa de enfriamiento del molde, respectivamente [27,28].

El <math display="inline">{Q}_{p}</math> y <math display="inline">{Q}_{r}</math> pueden ser representados por las siguientes expresiones:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math display="inline">{Q}_{p}={C}_{p}{\rho }_{p}{V}_{p}(\overline{{T}_{p}^{i}}-</math><math>\overline{{T}_{p}})</math>
|}
|  colspan='2'  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(7)
|-
|  colspan='2'  style="text-align: right;vertical-align: top;text-align: center;white-space: nowrap;"|<math>{Q}_{r}={h}_{c}{A}_{c}\int_{0}^{{t}_{er}}\left[ {\overline{T}}_{cs}\left( t\right) -{T}_{r}\right] dt</math>
|  colspan='2'  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(8)
|}

donde <math display="inline">{C}_{p}</math>, <math display="inline">{\rho }_{p}</math> y  <math display="inline">{V}_{p}</math> representan el calor específico, la densidad y el volumen de la masa fundida de polímero, respectivamente; <math display="inline">\overline{{T}_{p}^{i}}</math> es la temperatura inicial del polímero fundido; <math display="inline">\overline{{T}_{p}}</math> es la temperatura de eyección de la pieza; <math display="inline">{h}_{c}</math> es el coeficiente de transferencia de calor por convección en la interfase entre el refrigerante y los canales de enfriamiento; <math display="inline">{A}_{c}</math> es la superficie interna de los canales de refrigeración; <math display="inline">{t}_{er}</math> representa el tiempo de enfriamiento requerido; <math display="inline">{\overline{T}}_{cs}</math> indica la temperatura media de la pared de los canales de refrigeración en el instante t; y <math display="inline">{T}_{r}</math> es el valor promedio de la temperatura del refrigerante.

Al sustituir las ecs. (4), (7) y (8) en (6), se puede obtener la siguiente expresión:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{C}_{p}{\rho }_{p}{V}_{p}\left( \overline{{T}_{p}^{i}}-\overline{{T}_{p}}\right) +</math><math>{C}_{m}{\rho }_{m}{V}_{m}\left( \overline{{T}_{m}}-\overline{{T}_{m}^{i}}\right) +</math><math>{Q}_{c}={h}_{c}{A}_{c}\int_{0}^{{t}_{er}}\left[ {\overline{T}}_{cs}\left( t\right) -{T}_{r}\right] dt+</math><math>{Q}_{a}^{e}+{Q}_{i}^{e}</math>
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(9)
|}


Los efectos de <math display="inline">{Q}_{c}</math>, <math display="inline">{Q}_{a}^{e}</math> y <math display="inline">{Q}_{e}^{e}</math> en el campo de temperatura del molde durante el enfriamiento pueden ser despreciables, de modo que la ec. (10) se puede expresar como:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{C}_{p}{\rho }_{p}{V}_{p}\left( \overline{{T}_{p}^{i}}-\overline{{T}_{p}}\right) +</math><math>{C}_{m}{\rho }_{m}{V}_{m}\left( \overline{{T}_{m}}-\overline{{T}_{m}^{i}}\right) =</math><math>{h}_{c}{A}_{c}\int_{0}^{{t}_{er}}\left[ {\overline{T}}_{cs}\left( t\right) -{T}_{r}\right] dt</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(10)
|}


Mediante el uso del teorema del valor medio de las integrales, la ec. (10) se puede representar como la siguiente expresión:

{| class="formulaSCP" style="width: 49%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>{t}_{er}=\frac{{C}_{p}{\rho }_{p}{V}_{p}\left( \overline{{T}_{p}^{i}}-\overline{{T}_{p}}\right) +{C}_{m}{\rho }_{m}{V}_{m}\left( \overline{{T}_{m}}-\overline{{T}_{m}^{i}}\right) }{{h}_{c}{A}_{c}\left[ {\overline{T}}_{cs}\left( \epsilon \right) -{T}_{r}\right] }</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(11)
|}
para valores  <math display="inline">0\leq \epsilon \leq {t}_{er}</math>.

Para el moldeo de una parte de plástico específica, se predetermina el material polimérico correspondiente, el volumen de pieza y los parámetros de procesamiento, es decir, el <math display="inline">{C}_{p}</math>, <math display="inline">{\rho }_{p}</math>, <math display="inline">{V}_{p}</math>, <math display="inline">\overline{{T}_{p}^{i}}</math>, <math display="inline">{T}_{p}</math>, <math display="inline">\overline{{T}_{m}}</math>, <math display="inline">\overline{{T}_{m}^{i}}</math> en la ec. (11) pueden considerarse constantes. Por lo tanto, el <math display="inline">{t}_{er}</math> depende principalmente de los parámetros <math display="inline">{C}_{m}</math>, <math display="inline">{\rho }_{m}</math>, <math display="inline">{V}_{m}</math>, <math display="inline">{A}_{c}</math>  y <math display="inline">{T}_{r}</math>. Para mejorar la eficiencia de enfriamiento del molde, es decir, acortar el valor de <math display="inline">{t}_{er}</math> en la ec. (11), se necesita disminuir los valores de <math display="inline">{C}_{m}</math>, <math display="inline">{\rho }_{m}</math>, <math display="inline">{V}_{m}</math> y <math display="inline">{T}_{r}</math> y/o aumentar los valores de <math display="inline">{h}_{c}</math> y <math display="inline">{A}_{c}</math>.  La opción de minimizar el volumen del molde para mejorar la eficiencia de calentamiento del molde también puede mejorar la eficiencia de enfriamiento del molde. Además, la eficiencia de enfriamiento del molde puede mejorarse adicionalmente aumentando el caudal de refrigerante, bajando la temperatura del refrigerante y aumentando el área de enfriamiento del molde.

==3. Análisis térmico y estructural (Análisis térmico-estructural)==

El análisis térmico transitorio y el análisis estructural se implementaron con base en un software comercial de ANSYS® para analizar la respuesta térmica y la resistencia estructural del MRTR. Para este molde, el sistema de calentamiento por resistencias eléctricas se puede utilizar para un calentamiento uniforme y eficiente. La geometría propuesta para la cavidad del molde es una probeta de tensión, la sección de mayor importancia es donde se concentra la carga y se provoca la falla (cuello de la probeta) y por consecuencia es el área de análisis de este trabajo. La Figura 3(a), muestra la geometría total del inserto de la cavidad que será calentado, en la Figura 3 (b), se observa la simplificación de la geometría para la simulación térmica-estructural, esta simplificación es el resultado de parte central de la probeta de tensión, que es donde se concentra los esfuerzos y la falla mecánica.

En esta ilustración <math>a, h</math>, y <math>q</math> representan la distancia de centro a centro de las resistencias, la altura del centro de la resistencia a la superficie de la cavidad y el flujo de calor por unidad de área, respectivamente.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image3.png|center|382px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 3'''. Geometría del inserto de la cavidad; (a) Geometría simétrica para la simulación; (b) Geometría total
|}


El material del molde es un acero pre-endurecido (AISI P20) y sus propiedades físicas se describen en la Tabla 1.


<div class="center" style="font-size: 75%;">'''Tabla 1'''. Propiedades físicas del inserto del molde [29]
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:30%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|'''Propiedad física'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor'''
|-
|  style="border-top: 1pt solid black;vertical-align: top;"|Material
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|AISI P-20
|-
|  style="vertical-align: top;"|Densidad
|  style="text-align: center;vertical-align: top;"|7850 kg/ <math display="inline">{m}^{3}</math>
|-
|  style="vertical-align: top;"|Calor específico
|  style="text-align: center;vertical-align: top;"|460 kJ/°C
|-
|  style="vertical-align: top;"|Conductividad térmica 
|  style="text-align: center;vertical-align: top;"|34 W/m°C
|-
|  style="vertical-align: top;"|Módulo de elasticidad
|  style="text-align: center;vertical-align: top;"|205 GPa
|-
|  style="border-bottom: 1pt solid black;vertical-align: top;"|Coeficiente de Poisson
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.285
|}


===3.1. Análisis de respuesta térmica ===

El modelo de malla y las condiciones de frontera en el análisis transitorio de transferencia de calor se muestran en la Figura 4. Se considera la temperatura inicial de 60°C, equivalente a la temperatura de trabajo en un proceso industrial en un periodo transitorio estable.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image4.jpeg|center|394px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 4'''. Modelo de malla y condiciones de frontera para análisis transitorio de transferencia de calor
|}


El ciclo de moldeo en el MRTR inicia abierto y sin fluido refrigerante en el sistema de enfriamiento, el sistema de enfriamiento y calentamiento está dividido por una delgada capa de aislamiento (Figura 3c),  se eleva la temperatura del molde por medio de las resistencias eléctricas hasta alcanzar la <math display="inline">{T}_{c}</math> promedio en la superficie de la cavidad, la condición de convección de aire libre se cargó en la superficie externa del molde y el coeficiente de transferencia de calor convectivo correspondiente fue de 20 <math display="inline">w/{m}^{2}</math>°C [30]. Para la parte lateral izquierda se implementó una condición de simetría que permite reducir el modelo de malla, en cuanto a la parte lateral derecha, se aplicó un flujo de calor de cero por defecto debido a la capa de aislamiento. El medio para calentar el inserto del molde son resistencias eléctricas.

===3.2. Análisis térmico-estructural===

En la fase de llenado y empaquetamiento, la superficie de la cavidad del molde sufrira una presión de cavidad grande. Por lo tanto, el molde debe ser lo suficientemente fuerte para reducir la deformación de la superficie de la caviad y mejorar su durabilidad. La posicion de las resistencias electricas tiene un gran efecto sobre la eficiencia de la respuesta termica del molde, ademas de tener un efecto significativo en su resistencia estructural. Por ejemplo, situar las resistencias eléctricas mucho cerca de la supercie del molde puede aumentar la eficiencia de calentamiento, pero disminuir la resistencia estructural del molde. El volumen de la placa de acero que se desea calentar es otro factor importante que se requiere estudiar para reducir el tiempo de calentamiento en el ciclo de moldeo, al reducir el volumen de la placa del sistema de calentamiento, en consecuencia, el tiempo de calentamiento se reducirá. Como resultado, es necesario estudiar el espesor de la placa del sistema de calentamiento y el efecto de la posición de las resistencias eléctricas sobre la resistencia estructural del molde. La presión del polímero causada por la etapa de llenado se calculó mediante una simulación de moldeo por inyección usando Moldex3D. Los parámetros de moldeo por inyección se eligieron considerando el uso de un polietileno de alta densidad (PEMEX, PADMEX 60120-U) con condiciones de  trabajo mostradas en la Tabla 2.

<div class="center" style="font-size: 75%;">'''Tabla 2'''. Parámetros de inyección
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:45%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;vertical-align: top;"|'''Parámetro '''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor'''
|-
|  style="border-top: 1pt solid black;vertical-align: top;"|Temperatura de  inyección
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|200 <sup>0</sup>C
|-
|  style="vertical-align: top;"|Temperatura del molde
|  style="text-align: center;vertical-align: top;"|60 <sup>0</sup>C
|-
|  style="vertical-align: top;"|Tiempo de inyección de moldeo 
|  style="text-align: center;vertical-align: top;"|1 s
|-
|  style="vertical-align: top;"|Llenado de cavidad  
|  style="text-align: center;vertical-align: top;"|98% del volumen total de llenado
|-
|  style="vertical-align: top;"|Presión de empaquetamiento 
|  style="text-align: center;vertical-align: top;"|80MPa
|-
|  style="border-bottom: 1pt solid black;vertical-align: top;"|T. de presión de empaquetamiento
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|5.13 s
|}


En la Figura 5, se muestra el perfil de presión de la cavidad a lo largo del tiempo. La presión de inyección máxima, alcanzada en correspondencia con el corredor principal, fue de 65 MPa. A través de esta simulación por volumenes finitos, la distribución de carga se determinó y se usó como condición de frontera en la simulación estructural subsiguiente. Sin embargo, para simplificar la configuración de la simulación estructural y por razones de seguridad, se consideró una presión de inyección de 90 MPa aplicada en la superficie de la cavidad, aunque el valor maximo de la presion de inyección se alcanzó solo en la posición central, debajo del bebedero.

La Figura 6, muestra el modelo de malla y las condiciones de frontera para el análisis de resistencia estructural. En esta ilustración, P representa la presión de inyección.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image5.png|center|482px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 5'''. Perfil de presión de la cavidad del molde en la etapa de llenado
|}

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image6.jpeg|center|494px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 6'''. Modelo de malla y condiciones del análisis térmico-estructural
|}


En el límite izquierdo del modelo, la condición de simetría es implementada, en le límite derecho, el desplazamiento en la dirección normal es cero. Dado que los insertos de molde son soportados por la placa de retención para la fijación, el desplazamiento en el borde inferior del modelo es también cero en la dirección vertical.

==4. Trabajo experimental ==

Para estudiar la relación entre el espesor del inserto del sistema de calentamiento, la posición de las resistencias eléctricas, la eficiencia de respuesta térmica, uniformidad de temperatura y resistencia estructural del MRTR, se implementó la metodología de la superficie de respuesta combinada con la regresión mixta para obtener los correspondientes modelos matemáticos.

===4.1. Modelo de superficie de respuesta===

El modelo de superficie de respuesta es una función analítica útil para describir la relación numérica entre las variables experimentales controlables y las variables de respuesta. En MSR, la forma cuantitativa de relación entre las variables experimentales y la variable de respuesta deseada puede expresarse mediante la siguiente expresión:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>Y=f\left( {x}_{1},{x}_{2},\ldots ,{x}_{k}\right) +\epsilon</math> 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(12)
|}

donde ''Y'' representa la variable de respuesta, ''f'' es la función de respuesta desconocida, <math display="inline">{x}_{1},{x}_{2},\ldots ,{x}_{k}</math> son las variables experimentales independientes y <math display="inline">\epsilon</math>  es el error de ajuste que representa otras fuentes de variabilidad no contabilizadas para ''f''. El objetivo principal del MSR es establecer una función aproximada entre los factores y las respuestas a través de una serie especial de experimentos y análisis estadístico. Para MSR, hay muchos modelos que se pueden utilizar para establecer las funciones para aproximaciones, pero el modelo más común es el polinomio basado en la expansión de la serie de Taylor. En este estudio, la aproximación de ''f'' será propuesta utilizando el modelo de regresión polinomial de segundo orden, que puede escribirse de la siguiente manera:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>Y={\beta }_{0}+\sum _{i=1}^{k}{\beta }_{i}{x}_{i}+\sum _{i=1}^{k}{\beta }_{ii}{{x}_{i}}^{2}+</math><math>\sum _{i=1}^{k-1}\sum _{j=2}^{k}{\beta }_{ij}{x}_{i}{x}_{j}+\epsilon</math> 
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(13)
|}

donde <math display="inline">{\beta }_{0}</math> es la constante del modelo; <math display="inline">{\beta }_{i}</math> representa el término lineal; <math display="inline">{\beta }_{ii}</math> denota el término cuadrático; <math display="inline">{\beta }_{ij}</math> es el término de interacción; ''k'' es el número de factores o variables; y <math display="inline">\epsilon</math>  es el error estadístico. Se debe implementar una serie de experimentos reales o de simulación para obtener los valores de la variable de respuesta, que se pueden registrar en el vector ''Y''. La función aproximada entre las respuestas y los términos del modelo basados en el modelo polinomial cuadrático puede escribirse en la forma matricial de la siguiente manera:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>Y=X\beta +\epsilon</math> 
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(14)
|}

donde ''X'' es una matriz de términos de modelo incluyendo los términos lineales, términos cuadrados y términos de interacción; <math display="inline">\beta</math>  es un vector de columna que contiene los coeficientes de la pendiente; y <math display="inline">\, \epsilon</math>  es el vector de error residual. El MSR se basan en tres tipos de matrices de diseño de experimentos (DOE), incluyendo diseños de compuestos centrales (DCC), diseño de Box-Behnken (DBB) y estimación de error cuadrático medio integrado óptimo (EIMSE-óptimo). Para evaluar el vector de coeficiente de regresión <math display="inline">\beta</math>  del modelo predictivo, se puede utilizar el método de mínimos cuadrados. Se puede describir con la siguiente ecuación en la forma matriz:

{| class="formulaSCP" style="width: 47%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>\beta ={({X}^{T}X)}^{-1}{X}^{T}Y</math>
|}
|  style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(15)    
|}


donde <math display="inline">{X}^{T}</math> es la transpuesta de la matriz ''X''. Para evitar la singularidad de la matriz <math display="inline">{X}^{T}X</math>, debe observarse que el número de puntos de diseño experimental debe ser mayor que el número de los coeficientes de la pendiente.

===4.2. Diseño Box-Behnken===

Para el diseño experimental, se seleccionó el método DBB porque ofrece un equilibrio razonable entre precisión de predicción y economía de ejecución. Requiere tres niveles de cada factor y es siempre de forma esférica. El diseño Box-Behnken fue específicamente seleccionado ya que requiere menos ejecuciones que un DCC en casos de tres o cuatro variables. Este diseño cúbico se caracteriza por un conjunto de puntos que se encuentran en el punto medio de cada borde de un cubo multidimensional y el punto central se replica (<math>n=5</math>) mientras que las "esquinas perdidas" ayudan al experimentador a evitar los factores extremos combinados. Esta propiedad evita una posible pérdida de datos en esos casos. Como se muestra en Figura 3(a), el sistema de calentamiento está definido por cuatro variables: la distancia del centro de la resistencia a la superficie de la cavidad, ''h''; la distancia de centro a centro de las resistencia, ''a'';  el flujo de calor por unidad de área, ''q''; y el espesor de la placa del sistema de calentamiento, ''b''. Por lo tanto, se eligió un diseño de ensayo Box-Behnken con 29 ensayos para cuatro factores y tres niveles. Entre los ensayos experimentales, 24 ensayos estan en el borde del espacio experimental y los otros se replicaron en el centro del cubo espacial. Las rangos de los cuatro factores se pueden definir basándose en la experiencia previa y el análisis de ingeniería.

Los tres niveles de los factores se enumeran en la Tabla 3. -1, 0 y +1 representan los niveles bajo, medio y alto de los factores, respectivamente. El diseño del experimento basado en DBB se puede ver en la Tabla 4.

<div class="center" style="font-size: 75%;">'''Tabla 3'''. Niveles de los factores
</div>

{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:40%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: left;vertical-align: top;"|'''Factor'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Nivel -1'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Nivel 0'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Nivel +1'''
|-
|  style="border-top: 1pt solid black;vertical-align: top;"|h (''mm'')
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|8.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|10.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|12.00
|-
|  style="vertical-align: top;"|a (''mm'')
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|22.00
|-
|  style="vertical-align: top;"|''q'' (''w/'' <math display="inline">{m}^{2}</math>)
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|3.00
|-
|  style="border-bottom: 1pt solid black;vertical-align: top;"|b (mm)
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|20.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|22.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|24.00
|}


Para evaluar la influencia del sistema en la eficiencia de calentamiento, la uniformidad de la temperatura y la resistencia estructural, se definen tres variables de respuesta. Se utilizó el tiempo de calentamiento requerido, <math display="inline">{t}_{cr}</math>, para que toda la superficie del molde se calentara a la temperatura designada.

<div class="center" style="font-size: 75%;">'''Tabla 4'''. Resultados del análisis de los puntos experimentales seleccionados
</div>

{| style="width: 80%;margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''No. de exp.'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''h (mm)'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''a (mm)'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''q (w/''' <math display="inline">{\mathit{\boldsymbol{m}}}^{\mathit{\boldsymbol{2}}}</math>''')'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''b (mm)'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{cr}}}</math>'''  (s)'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{d}}}</math>'''   (s)'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math display="inline">{\mathit{\boldsymbol{\sigma }}}_{\mathit{\boldsymbol{v}}}</math>''' (MPa)'''
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|1
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|10.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|14.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|2.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|20.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|17.8
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|10.68
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|448
|-
|  style="text-align: center;vertical-align: top;"|2
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|28.9
|  style="text-align: center;vertical-align: top;"|3.564
|  style="text-align: center;vertical-align: top;"|545
|-
|  style="text-align: center;vertical-align: top;"|3
|  style="text-align: center;vertical-align: top;"|8.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|20.00
|  style="text-align: center;vertical-align: top;"|18.85
|  style="text-align: center;vertical-align: top;"|4.564
|  style="text-align: center;vertical-align: top;"|499
|-
|  style="text-align: center;vertical-align: top;"|4
|  style="text-align: center;vertical-align: top;"|8.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|49.7
|  style="text-align: center;vertical-align: top;"|2.665
|  style="text-align: center;vertical-align: top;"|525
|-
|  style="text-align: center;vertical-align: top;"|5
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|20.00
|  style="text-align: center;vertical-align: top;"|47.25
|  style="text-align: center;vertical-align: top;"|3.921
|  style="text-align: center;vertical-align: top;"|511
|-
|  style="text-align: center;vertical-align: top;"|6
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|28.45
|  style="text-align: center;vertical-align: top;"|5.756
|  style="text-align: center;vertical-align: top;"|505
|-
|  style="text-align: center;vertical-align: top;"|7
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|29.4
|  style="text-align: center;vertical-align: top;"|1.891
|  style="text-align: center;vertical-align: top;"|543
|-
|  style="text-align: center;vertical-align: top;"|8
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|20.00
|  style="text-align: center;vertical-align: top;"|13.7
|  style="text-align: center;vertical-align: top;"|7.22
|  style="text-align: center;vertical-align: top;"|505
|-
|  style="text-align: center;vertical-align: top;"|9
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|15.8
|  style="text-align: center;vertical-align: top;"|6.083
|  style="text-align: center;vertical-align: top;"|502
|-
|  style="text-align: center;vertical-align: top;"|10
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|16.95
|  style="text-align: center;vertical-align: top;"|7.814
|  style="text-align: center;vertical-align: top;"|530
|-
|  style="text-align: center;vertical-align: top;"|11
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|20.00
|  style="text-align: center;vertical-align: top;"|24.45
|  style="text-align: center;vertical-align: top;"|2.892
|  style="text-align: center;vertical-align: top;"|538
|-
|  style="text-align: center;vertical-align: top;"|12
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|23.4
|  style="text-align: center;vertical-align: top;"|5.716
|  style="text-align: center;vertical-align: top;"|500
|-
|  style="text-align: center;vertical-align: top;"|13
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|23.4
|  style="text-align: center;vertical-align: top;"|5.713
|  style="text-align: center;vertical-align: top;"|499
|-
|  style="text-align: center;vertical-align: top;"|14
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|48.65
|  style="text-align: center;vertical-align: top;"|5.795
|  style="text-align: center;vertical-align: top;"|473
|-
|  style="text-align: center;vertical-align: top;"|15
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|23.4
|  style="text-align: center;vertical-align: top;"|5.714
|  style="text-align: center;vertical-align: top;"|499
|-
|  style="text-align: center;vertical-align: top;"|16
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|23.4
|  style="text-align: center;vertical-align: top;"|5.711
|  style="text-align: center;vertical-align: top;"|499
|-
|  style="text-align: center;vertical-align: top;"|17
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|20.00
|  style="text-align: center;vertical-align: top;"|22.55
|  style="text-align: center;vertical-align: top;"|7.166
|  style="text-align: center;vertical-align: top;"|578
|-
|  style="text-align: center;vertical-align: top;"|18
|  style="text-align: center;vertical-align: top;"|8.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|12.5
|  style="text-align: center;vertical-align: top;"|3.882
|  style="text-align: center;vertical-align: top;"|502
|-
|  style="text-align: center;vertical-align: top;"|19
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|57.6
|  style="text-align: center;vertical-align: top;"|1.431
|  style="text-align: center;vertical-align: top;"|542
|-
|  style="text-align: center;vertical-align: top;"|20
|  style="text-align: center;vertical-align: top;"|8.00
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|15.95
|  style="text-align: center;vertical-align: top;"|9.935
|  style="text-align: center;vertical-align: top;"|449
|-
|  style="text-align: center;vertical-align: top;"|21
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|22.85
|  style="text-align: center;vertical-align: top;"|9.148
|  style="text-align: center;vertical-align: top;"|482
|-
|  style="text-align: center;vertical-align: top;"|22
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|17.85
|  style="text-align: center;vertical-align: top;"|2.416
|  style="text-align: center;vertical-align: top;"|548
|-
|  style="text-align: center;vertical-align: top;"|23
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|3.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|12.15
|  style="text-align: center;vertical-align: top;"|11.86
|  style="text-align: center;vertical-align: top;"|447
|-
|  style="text-align: center;vertical-align: top;"|24
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|58.45
|  style="text-align: center;vertical-align: top;"|2.983
|  style="text-align: center;vertical-align: top;"|515
|-
|  style="text-align: center;vertical-align: top;"|25
|  style="text-align: center;vertical-align: top;"|8.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|21.8
|  style="text-align: center;vertical-align: top;"|3.564
|  style="text-align: center;vertical-align: top;"|507
|-
|  style="text-align: center;vertical-align: top;"|26
|  style="text-align: center;vertical-align: top;"|12.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|1.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|55.85
|  style="text-align: center;vertical-align: top;"|3.571
|  style="text-align: center;vertical-align: top;"|509
|-
|  style="text-align: center;vertical-align: top;"|27
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|18.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|22.00
|  style="text-align: center;vertical-align: top;"|23.4
|  style="text-align: center;vertical-align: top;"|5.716
|  style="text-align: center;vertical-align: top;"|499
|-
|  style="text-align: center;vertical-align: top;"|28
|  style="text-align: center;vertical-align: top;"|10.00
|  style="text-align: center;vertical-align: top;"|14.00
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|24.00
|  style="text-align: center;vertical-align: top;"|21.3
|  style="text-align: center;vertical-align: top;"|9.462
|  style="text-align: center;vertical-align: top;"|450
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|29
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|8.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|22.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|22.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|24.65
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.814
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|540
|}


El polímero inyectado es un PADMEX 60120-U suministrado por la empresa PEMEX, el cual presenta una <math display="inline">{T}_{c}</math> de 116 [31,32], temperatura que se utiliza para evaluar la eficiencia de calentamiento en el presente trabajo. El esfuerzo máximo de Von Mises, <math display="inline">{\sigma }_{v}</math>, ejercido sobre el molde se utilizó para evaluar la resistencia estructural. La diferencia de las temperaturas máxima y mínima superficial del molde, <math display="inline">{t}_{d}</math>. Los resultados de la simulación con los parámetros iniciales, es decir, todos los factores en nivel cero, muestran que la mayoría de la superficie de la cavidad se ha calentado a la temperatura designada después de un calentamiento de 23.4 s.

La Figura 7 muestra la gráfica de contorno de la distribución de la temperatura en la superficie de la cavidad cuando se termina el proceso de calentamiento para el sistema de calentamiento del molde con los parámetros iniciales. El largo tiempo de calentamiento y diferencia de temperatura de la superficie de la cavidad fuera de norma conducirá a un largo tiempo de ciclo de moldeo y baja calidad del producto.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image7.png|center|482px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 7'''. Gráfica de contorno de distribución de temperaturas con parámetros iniciales
|}


La Figura 8 muestra la distribución de los esfuerzos de Von Mises en la placa de calentamiento en los parámetros iniciales del diseño de experimentos, se observa en la posición de las resistencias eléctricas es donde se concentra los esfuerzos máximos.

<span id='_Ref513457402'></span>
{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image8.png|center|482px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 8'''. Esfuerzos de Von Mises con parámetros iniciales
|}

==5. Resultados y discusión==

De acuerdo con el resultado de DBB, se realizaron análisis de respuesta térmica y análisis de resistencia estructural basado en análisis de elementos finitos (FEA, por sus siglas en inglés) para obtener los valores de las diferentes variables de respuesta. Los resultados experimentales se presentan en la Tabla 4. Las variables de respuesta obtenidas de los experimentos se analizaron por regresión múltiple lineal para determinar los modelos matemáticos con mejores resultados. La adecuación y fiabilidad de los modelos de regresión también se probaron mediante el análisis de varianza (ANOVA).

===5.1 Análisis de varianza===

Se aplicó un ANOVA para examinar la importancia de los términos del modelo cuadrático. Los resultados de la ANOVA para el modelo cuadrático del tiempo de calentamiento requerido, <math display="inline">{t}_{cr}</math>, se presentan en la Tabla 5.


<div class="center" style="font-size: 75%;">'''Tabla 5'''. ANOVA para el tiempo de calentamiento requerido
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:70%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Fuente'''
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Sum. de cuadrados'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''DF'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''M. de Cuadrado'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor F'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor P (Prob>F)'''
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|Modelo
|  colspan='2'  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|5511.50
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|9
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|612.39
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|2019.83
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|A-h
|  colspan='2'  style="text-align: center;vertical-align: top;"|85.87
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|85.87
|  style="text-align: center;vertical-align: top;"|283.21
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|B-a
|  colspan='2'  style="text-align: center;vertical-align: top;"|162.44
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|162.44
|  style="text-align: center;vertical-align: top;"|535.76
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|C-w
|  colspan='2'  style="text-align: center;vertical-align: top;"|4352.93
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|4352.93
|  style="text-align: center;vertical-align: top;"|14357.1
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|D-b
|  colspan='2'  style="text-align: center;vertical-align: top;"|78.03
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|78.03
|  style="text-align: center;vertical-align: top;"|257.36
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|AB
|  colspan='2'  style="text-align: center;vertical-align: top;"|1.76
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|1.76
|  style="text-align: center;vertical-align: top;"|5.79
|  style="text-align: center;vertical-align: top;"|0.0265
|-
|  style="text-align: center;vertical-align: top;"|AD
|  colspan='2'  style="text-align: center;vertical-align: top;"|2.18
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|2.18
|  style="text-align: center;vertical-align: top;"|7.18
|  style="text-align: center;vertical-align: top;"|0.0149
|-
|  style="text-align: center;vertical-align: top;"|BC
|  colspan='2'  style="text-align: center;vertical-align: top;"|2.64
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|2.64
|  style="text-align: center;vertical-align: top;"|8.71
|  style="text-align: center;vertical-align: top;"|0.0082
|-
|  style="text-align: center;vertical-align: top;"|CD
|  colspan='2'  style="text-align: center;vertical-align: top;"|20.70
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|20.70
|  style="text-align: center;vertical-align: top;"|68.28
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|<math>{C}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|804.97
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|804.97
|  style="text-align: center;vertical-align: top;"|2655.02
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Residual
|  style="text-align: center;vertical-align: top;"|5.76
|  style="text-align: center;vertical-align: top;"|19
|  style="text-align: center;vertical-align: top;"|0.30
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|F. de ajuste
|  style="text-align: center;vertical-align: top;"|5.76
|  style="text-align: center;vertical-align: top;"|15
|  style="text-align: center;vertical-align: top;"|0.38
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  style="text-align: center;vertical-align: top;"|Error puro
|  colspan='2'  style="text-align: center;vertical-align: top;"|0.000
|  style="text-align: center;vertical-align: top;"|4
|  style="text-align: center;vertical-align: top;"|0.000
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  style="text-align: center;vertical-align: top;"|Total
|  colspan='2'  style="text-align: center;vertical-align: top;"|5517.27
|  style="text-align: center;vertical-align: top;"|28
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Desv. Est.
|  style="text-align: center;vertical-align: top;"|0.55
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada''
|  style="text-align: center;vertical-align: top;"|0.9990
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Media
|  style="text-align: center;vertical-align: top;"|27.60
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada ajustada''
|  style="text-align: center;vertical-align: top;"|0.9985
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|C.V.%
|  style="text-align: center;vertical-align: top;"|2.00
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada de la predicción''
|  style="text-align: center;vertical-align: top;"|0.9969
|-
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|PRESS
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|17.20
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''Precisión adecuada''
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|141.596
|}


La suma de cuadrados en la tabla se refiere a la suma de cuadrados para cada término, que mide la variabilidad en los datos aportados por el término. DF representa el número de términos que contribuirán a la predicción del error. El valor del modelo F de 2019.83 indica que el modelo se considera estadísticamente significativo. Los valores de Valor P (Prob>F) inferiores a 0.0500 indican que los términos del modelo son significativos. En este caso A, B, C, D, AB, AD, BC, CD, <math display="inline">{C}^{2}</math> son términos de modelo significativos. Los valores superiores a 0.1000 indican que los términos del modelo no son significativos. La ''R-Cuadrada de la predicción'' es de 0.9969 está en un acuerdo razonable con la ''R-Cuadrada ajustada'' de 0.9985. La ''Precisión adecuada'' mide la relación señal/ruido. Es deseable una relación mayor que 4. La relación es de 141.596 lo que indica una señal adecuada. Este modelo se puede utilizar para navegar por el espacio de diseño.

Los resultados del análisis de ANOVA para la diferencia de las temperaturas máxima y mínima superficial del molde, <math display="inline">{t}_{d}</math>, se presentan en la Tabla 6. El valor del modelo F es de 188.18 implica que el modelo es significativo. Los valores de Valor P (Prob>F) inferiores a 0.0500 indican que los términos del modelo son significativos.  En este caso A, B, C, D, AB, AC, BC, <math display="inline">{A}^{2}</math>, <math display="inline">{B}^{2}</math> y <math display="inline">\, {C}^{2}</math> son términos de modelo significativos. La ''R-Cuadrada de la predicción'' de 0.9557 está en un acuerdo razonable con la ''R-Cuadrada ajustada'' de 0.9853. La ''Precisión adecuada'' mide la relación señal/ruido. Es deseable una relación mayor que 4. La relación de 54.749 indica una señal adecuada. Este modelo se puede utilizar para navegar por el espacio de diseño.


<div class="center" style="font-size: 75%;">'''Tabla 6'''. ANOVA para la diferencia de las temperaturas máxima y mínima superficial del molde
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:70%;"  
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Fuente'''
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Sum. de cuadrados'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''DF'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''M. de Cuadrado'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor F'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor P (Prob>F)'''
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|Modelo
|  colspan='2'  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|224.98
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|10
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|22.50
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|188.18
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|A-h
|  colspan='2'  style="text-align: center;vertical-align: top;"|11.20
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|11.20
|  style="text-align: center;vertical-align: top;"|93.71
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|B-a
|  colspan='2'  style="text-align: center;vertical-align: top;"|160.43
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|160.43
|  style="text-align: center;vertical-align: top;"|1341.9
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|C-w
|  colspan='2'  style="text-align: center;vertical-align: top;"|29.80
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|29.80
|  style="text-align: center;vertical-align: top;"|249.22
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|D-b
|  colspan='2'  style="text-align: center;vertical-align: top;"|3.75
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|3.75
|  style="text-align: center;vertical-align: top;"|31.37
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|AB
|  colspan='2'  style="text-align: center;vertical-align: top;"|3.13
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|3.13
|  style="text-align: center;vertical-align: top;"|26.16
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|AC
|  colspan='2'  style="text-align: center;vertical-align: top;"|2.29
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|2.29
|  style="text-align: center;vertical-align: top;"|19.15
|  style="text-align: center;vertical-align: top;"|0.0004
|-
|  style="text-align: center;vertical-align: top;"|BC
|  colspan='2'  style="text-align: center;vertical-align: top;"|6.45
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|6.45
|  style="text-align: center;vertical-align: top;"|53.96
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|<math>{A}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|1.29
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|1.29
|  style="text-align: center;vertical-align: top;"|10.76
|  style="text-align: center;vertical-align: top;"|0.0042
|-
|  style="text-align: center;vertical-align: top;"|<math>{B}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|1.64
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|1.64
|  style="text-align: center;vertical-align: top;"|13.73
|  style="text-align: center;vertical-align: top;"|0.0016
|-
|  style="text-align: center;vertical-align: top;"|<math>{C}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|4.17
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|4.17
|  style="text-align: center;vertical-align: top;"|34.89
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|Residual
|  colspan='2'  style="text-align: center;vertical-align: top;"|2.15
|  style="text-align: center;vertical-align: top;"|18
|  style="text-align: center;vertical-align: top;"|0.12
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|F.  de ajuste
|  style="text-align: center;vertical-align: top;"|2.15
|  style="text-align: center;vertical-align: top;"|14
|  style="text-align: center;vertical-align: top;"|0.15
|  style="text-align: center;vertical-align: top;"|34158
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Error puro
|  style="text-align: center;vertical-align: top;"|1.800E-005
|  style="text-align: center;vertical-align: top;"|4
|  style="text-align: center;vertical-align: top;"|4.5 E-6
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Total
|  style="text-align: center;vertical-align: top;"|227.13
|  style="text-align: center;vertical-align: top;"|28
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Desv. Est.
|  style="text-align: center;vertical-align: top;"|0.35
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada''
|  style="text-align: center;vertical-align: top;"|0.9905
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Media
|  style="text-align: center;vertical-align: top;"|5.44
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada ajustada''
|  style="text-align: center;vertical-align: top;"|0.9853
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|C.V.%
|  style="text-align: center;vertical-align: top;"|6.36
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada de la predicción''
|  style="text-align: center;vertical-align: top;"|0.9557
|-
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|PRESS
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|10.07
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''Precisión adecuada''
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|54.749
|}


Los resultados del análisis de ANOVA para la los Esfuerzos de Von Mises se presentan en la Tabla 7. El valor del modelo F es de 26.89 implica que el modelo es significativo. Los valores de Valor P (Prob>F) inferiores a 0.0500 indican que los términos del modelo son significativos.


<div class="center" style="font-size: 75%;">'''Tabla 7'''. ANOVA para el modelo cuadrático de los Esfuerzos de Von Mises
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:75%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Fuente'''
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Sum. de cuadrados'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''DF'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''M. de Cuadrado'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor F'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Valor P (Prob>F)'''
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|Modelo
|  colspan='2'  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|27053.77
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|8
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|3381.72
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|26.89
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|A-h
|  colspan='2'  style="text-align: center;vertical-align: top;"|1344.08
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|1344.08
|  style="text-align: center;vertical-align: top;"|10.69
|  style="text-align: center;vertical-align: top;"|0.0038
|-
|  style="text-align: center;vertical-align: top;"|B-a
|  colspan='2'  style="text-align: center;vertical-align: top;"|21420.75
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|21420.75
|  style="text-align: center;vertical-align: top;"|170.31
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  style="text-align: center;vertical-align: top;"|C-w
|  colspan='2'  style="text-align: center;vertical-align: top;"|140.08
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|140.08
|  style="text-align: center;vertical-align: top;"|1.11
|  style="text-align: center;vertical-align: top;"|0.3038
|-
|  style="text-align: center;vertical-align: top;"|D-b
|  colspan='2'  style="text-align: center;vertical-align: top;"|270.75
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|270.75
|  style="text-align: center;vertical-align: top;"|2.15
|  style="text-align: center;vertical-align: top;"|0.1579
|-
|  style="text-align: center;vertical-align: top;"|AC
|  colspan='2'  style="text-align: center;vertical-align: top;"|484.00
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|484.00
|  style="text-align: center;vertical-align: top;"|3.85
|  style="text-align: center;vertical-align: top;"|0.0639
|-
|  style="text-align: center;vertical-align: top;"|AD
|  colspan='2'  style="text-align: center;vertical-align: top;"|1640.25
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|1640.25
|  style="text-align: center;vertical-align: top;"|13.04
|  style="text-align: center;vertical-align: top;"|0.0017
|-
|  style="text-align: center;vertical-align: top;"|<math>{A}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|992.00
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|992.00
|  style="text-align: center;vertical-align: top;"|7.89
|  style="text-align: center;vertical-align: top;"|0.0109
|-
|  style="text-align: center;vertical-align: top;"|<math>{B}^{2}</math>
|  colspan='2'  style="text-align: center;vertical-align: top;"|529.80
|  style="text-align: center;vertical-align: top;"|1
|  style="text-align: center;vertical-align: top;"|529.80
|  style="text-align: center;vertical-align: top;"|4.21
|  style="text-align: center;vertical-align: top;"|0.0534
|-
|  style="text-align: center;vertical-align: top;"|Residual
|  colspan='2'  style="text-align: center;vertical-align: top;"|2515.47
|  style="text-align: center;vertical-align: top;"|20
|  style="text-align: center;vertical-align: top;"|125.77
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|F. de ajuste
|  style="text-align: center;vertical-align: top;"|2514.67
|  style="text-align: center;vertical-align: top;"|16
|  style="text-align: center;vertical-align: top;"|157.17
|  style="text-align: center;vertical-align: top;"|785.83
|  style="text-align: center;vertical-align: top;"|< 0.0001
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Error puro
|  style="text-align: center;vertical-align: top;"|0.80
|  style="text-align: center;vertical-align: top;"|4
|  style="text-align: center;vertical-align: top;"|0.20
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  style="text-align: center;vertical-align: top;"|Total
|  colspan='2'  style="text-align: center;vertical-align: top;"|29569.24
|  style="text-align: center;vertical-align: top;"|28
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|  style="text-align: center;vertical-align: top;"|
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Desv. Est.
|  style="text-align: center;vertical-align: top;"|11.21
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada''
|  style="text-align: center;vertical-align: top;"|0.9149
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|Media
|  style="text-align: center;vertical-align: top;"|506.52
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada ajustada''
|  style="text-align: center;vertical-align: top;"|0.8809
|-
|  colspan='2'  style="text-align: center;vertical-align: top;"|C.V.%
|  style="text-align: center;vertical-align: top;"|2.21
|  style="text-align: center;vertical-align: top;"|
|  colspan='2'  style="text-align: center;vertical-align: top;"|''R-Cuadrada de la predicción''
|  style="text-align: center;vertical-align: top;"|0.7453
|-
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|PRESS
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|7531.18
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|
|  colspan='2'  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|''Precisión adecuada''
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|17.898
|}


En este caso A, B, C, D, AC, AD, <math display="inline">{B}^{2}</math> y <math display="inline">{A}^{2}</math> son términos  significativos del modelo. Los valores superiores a 0.1000 indican que los términos del modelo no son significativos. La ''R-Cuadrada de la predicción'' de 0.753 está en un acuerdo razonable con ''la R-Cuadrada ajustada'' de 0.8809. La ''Precisión adecuada'' mide la relación señal/ruido. Es deseable una relación mayor que 4. La relación de 17.898 indica una señal adecuada. Este modelo se puede utilizar para navegar por el espacio de diseño.

===5.2. Modelo de regresión de las respuestas===

Los modelos de regresión de la superficie de respuesta para evaluar el tiempo de calentamiento requerido, la variación de la temperatura superficial del molde y los esfuerzos máximos de Von Mises se expresan como:

{| class="formulaSCP" style="width: 48%;border-collapse: collapse;width: 100%;text-align: center;" 
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>Tcr=14.36691-1.22813\ast h+2.15417\ast a\,</math>

<math>\quad \quad \quad -33.15380\, \ast \, w\, +\, 1.70625\, \ast b</math>

<math>\quad \quad \quad -0.082813\ast h\ast a\, +\, 0.18438\, \ast \, h\, \ast b</math>

<math>\quad \quad \quad -0.20312\, \ast a\ast w-1.13750\ast w\ast b\,</math>

<math display="inline">\quad \quad \quad +10.69730\ast {w}^{2}</math> 
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(16)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>Td=32.31387-0.076823\ast h-2.49598\ast a\,</math>

<math>\quad \quad \, \, +6.65814\ast w-0.27954\ast b\,</math>

<math>\quad \quad \, \, +0.11053\ast h\ast a+0.37825\, \ast h\ast w\,</math>

<math>\quad \quad \, \, -0.31750\ast a\ast w\, -0.10931\ast \, {h}^{2}</math>

<math>\quad \quad \, \, +0.030877\ast {a}^{2}-0.78747\ast {w}^{2}</math>
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(17)
|-
| 
{| style="text-align: center; margin:auto;width: 100%;"
|-
| <math>EVM=-560.22822+45.72348\ast h\,</math>

<math>\quad \quad \quad \, \, +30.2755\ast a-58.41667\, \ast w</math>

<math>\quad \quad \quad \, \, +48.250\ast b+5.50\ast h\ast w</math>

<math>\quad \quad \quad \, \, -5.06250\ast h\ast b+2.99716\ast {h}^{2}</math>

<math>\quad \quad \quad \, \, -0.54759\ast {a}^{2}</math>
|}
|  style="vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(18)
|}


Los modelos matemáticos de las ecs. (16)-(18) pueden usarse para predecir el tiempo de calentamiento requerido, <math display="inline">{t}_{cr}</math>, la variación de la temperatura superficial del molde, <math display="inline">\, {t}_{d}</math>, los esfuerzos máximos de Von Mises, <math display="inline">{\sigma }_{v}</math>. Los modelos matemáticos desarrollados anteriormente producen una forma simple y eficiente de predecir los valores de las variables de respuesta dentro de los límites de los factores estudiados.

Esto será útil para la optimización de parámetros de diseño en el siguiente paso. Antes de la aplicación de los modelos polinomiales cuadráticos, es necesario verificar su precisión confirmando que los modelos han extraído toda la información relevante de los resultados experimentales. Por esta razón, se seleccionó el análisis residual como la herramienta de diagnóstico [33]. Los diagramas de probabilidad normal de los residuales para el tiempo de calentamiento requerido (<math display="inline">{t}_{cr}</math>), diferencia de temperatura de superficie de cavidad máxima (<math display="inline">{t}_{d}</math>) y los esfuerzos de Von Mises (<math display="inline">{\sigma }_{v})</math> se muestran en la Figura 9(a)-(c), respectivamente.  Se puede encontrar que los residuos para cada variable de respuesta generalmente caen en una línea recta que implica que los errores se distribuyen normalmente y las adecuaciones de estos modelos matemáticos se ajustan mediante técnicas de regresión de mínimos cuadrados.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image9.png|center|382px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| <span id='_Ref511728576'></span><span id='_Ref511728571'></span>'''Figura 9'''. Gráficas de probabilidad normal. (a) Tiempo de calentamiento requerido. (b) Diferencia máxima de temperaturas de la superficie de la cavidad. (c) Esfuerzos de Von Mises
|}


Se realizó una prueba de confirmación con un conjunto de parámetros de diseño aleatorios para verificar los modelos matemáticos desarrollados. La Tabla 8 muestra los conjuntos de datos aleatorios para la verificación y la comparación de los resultados de cálculo por los modelos desarrollados y la simulación numérica. Se puede observar que los errores porcentuales calculados son muy pequeños. Las desviaciones porcentuales para las tres variables de respuesta entre los resultados estimados y los experimentales son de 1.85%,0.095% y 0.308%, respectivamente.

Estos pequeños errores indican que los modelos polinomiales desarrollados son adecuados y pueden usarse para predecir las variables objetivas con buena precisión. Se puede deducir, que es razonable usar estos modelos matemáticos para la optimización.

<div class="center" style="font-size: 75%;">'''Tabla 8'''. Comparación de los modelos polinomiales cuadráticos desarrollados y simulación con variables de diseño aleatorias
</div>

{| style="margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;width:75%;"
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Parámetros'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''h''' 

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''a''' 

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''q'''w/ <math display="inline">{m}^{2}</math>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''b''' 

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{cr}}}</math>

s
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{d}}}</math>

°C
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{\sigma }}}_{\mathit{\boldsymbol{v}}}</math>

MPa
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|V. modelo
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|11.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|19.50
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|1.50
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|23.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|40.14
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|3.760
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|519.6
|-
|  style="text-align: center;vertical-align: top;"|V. simulación
|  style="text-align: center;vertical-align: top;"|11.00
|  style="text-align: center;vertical-align: top;"|19.50
|  style="text-align: center;vertical-align: top;"|1.50
|  style="text-align: center;vertical-align: top;"|23.00
|  style="text-align: center;vertical-align: top;"|39.40
|  style="text-align: center;vertical-align: top;"|3.757
|  style="text-align: center;vertical-align: top;"|518
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|% Error
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|1.85
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.095
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|0.308
|}


Para visualizar el efecto de las variables en las respuestas, los diagramas de superficie de respuesta tridimensional se pueden usar para describir los modelos matemáticos extraídos. Para la respuesta establecida del tiempo de calentamiento requerido las superficies de respuesta correspondientes se muestran en la Figura 10. La Figura 10(a) muestra la gráfica de superficie de respuesta del efecto del término interactivo, la distancia de centro a centro de las resistencias y  la altura del centro de la resistencia a la superficie de la cavidad (ah) en el tiempo de calentamiento requerido, se puede observar  el tiempo de calentamiento requerido se puede reducir a medida que se reduce la distancia de centro a centro de las resistencias y  la altura del centro de la resistencia a la superficie y aumenta de manera contraria.

El tiempo de calentamiento mínimo requerido sería de 12.15 s, correspondiente a la distancia más pequeña de la distancia de centro a centro de las resistencias (<math>a = 14 mm</math>) y la menor la altura del centro de la resistencia a la superficie de la cavidad (<math>h = 8 mm</math>). En la Figura 10(b) se observa la gráfica de superficie de respuesta del efecto del término interactivo, la distancia de centro a centro de las resistencias y el flujo de calor por unidad de área (aq) en el tiempo de calentamiento requerido, se puede interpretar que el aumento del flujo de calor da como resultado la reducción del tiempo de calentamiento requerido, especialmente cuando la distancia de centro a centro de las resistencias (a) está en su nivel bajo.

La Figura 10(c) muestra la gráfica de superficie de respuesta del efecto del término interactivo, el espesor de la placa de calentamiento y el flujo de calor por unidad de área en el tiempo de calentamiento requerido, el tiempo mínimo de calentamiento requerido se puede obtener con el nivel bajo del espesor de la placa de calentamiento (<math>b= 20 mm</math>) y el nivel alto del flujo de calor por unidad de área (<math display="inline">q= 3 w/{m}^{2}</math>).

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image10.png|center|382px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 10'''. Gráficas de superficie de los efectos combinados de (a) la distancia de centro a centro de las resistencias y la altura del centro de la resistencia a la superficie de la cavidad (ah), (b) la distancia de centro a centro de las resistencias y el flujo de calor por unidad de área (aq) y (c) el espesor de la placa de calentamiento y el flujo de calor por unidad de área (bq) para el tiempo de calentamiento
|}

==6. Optimización del sistema de calentamiento==

Según la discusión anterior, se podría concluir que el espesor de la placa y la posición las resistencias eléctricas tienen efectos significativos y complejos sobre la eficiencia de calentamiento, la uniformidad de temperatura en cavidad del molde y la resistencia del molde.

Dado que la productividad es un indicador importante en la industria de moldeo por inyección, se desea que el tiempo de calentamiento sea lo más corto posible. Por lo tanto, el propósito de esta optimización es minimizar el tiempo de calentamiento requerido para mejorar la eficiencia de calentamiento. Al mismo tiempo, la variación  de la temperatura de la superficie del molde de acuerdo a la norma ASTM-D 3641–02 y D 4703–03 [34,35], debe ser inferior a 5°C para garantizar la calidad de la superficie de la pieza. Los esfuerzos máximos de Von Mises no deben ser superiores a 500 MPa. Además, los valores de las tres variables de diseño están todos en el espacio de diseño designado, como se muestra en la Tabla 3. El problema de optimización se puede representar mediante las siguientes ecuaciones:

* Encontrar: ''h'', ''a'', ''q'', ''b''

* Minimizar:

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;" | <math>{t}_{cr}(h,a,q,b)</math>
|}
| style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(19)
|}

* Sujeto a restricciones: 

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;"|<math display="inline">{t}_{d}\left( h,a,q,\, b\right) \leq 5\hbox{°C}d\qquad \hbox{y}\qquad  {\sigma }_{v}\left( h,a,q,b\right) \leq 500\, MPa</math>
|}
| style="text-align: right;vertical-align: top;width: 5px;text-align: right;white-space: nowrap;"|(20)
|}

* Dentro de los rangos: 

{| class="formulaSCP" style="width: 100%; text-align: left;" 
|-
| 
{| style="text-align: left; margin:auto;width: 100%;" 
|-
| style="text-align: center;"| <math display="inline">8mm\leq h\leq 12mm</math>; <math display="inline">14mm\leq a\leq 22mm</math>
|-
| style="text-align: center;"|<math display="inline">1w/{m}^{2}\leq q\leq 3w/{m}^{2}</math>; <math display="inline">20m\leq b\leq 24mm</math> 
|}
| style="text-align: right;vertical-align: center;width: 5px;text-align: right;white-space: nowrap;"|(21)
|}

El problema de optimización expresado con las ecuaciones anteriores se resolvió acoplando los modelos matemáticos desarrollados con el algoritmo PSO. El diagrama de flujo del diseño de optimización se muestra en la Figura 11.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image11.png|center|650px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 11'''. Metodología de diseño óptimo basado en DOE, MSR y PSO
|}

Los valores de las variables dependientes y las variables de respuesta correspondientes antes y después de la optimización se presentan en la Tabla 9.

<div class="center" style="font-size: 75%;">'''Tabla 9'''. Comparación de los valores iniciales y óptimos de diseño
</div>

{| style="width: 50%;margin: 1em auto 0.1em auto;border-collapse: collapse;font-size:85%;" 
|-
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''Parámetros'''
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''h'''

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''a'''

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''q'''

w/ <math display="inline">{m}^{2}</math>
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|'''b'''

mm
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{cr}}}</math>

s
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{t}}}_{\mathit{\boldsymbol{d}}}</math>

°C
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|<math>{\mathit{\boldsymbol{\sigma }}}_{\mathit{\boldsymbol{v}}}</math>

MPa
|-
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|V. inicial
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|10.0
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|18.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|2.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|22.00
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|23.4
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|5.714
|  style="border-top: 1pt solid black;text-align: center;vertical-align: top;"|499
|-
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|V. optimizado
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|8.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|18.03
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|2.82
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|20.00
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|12.15
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|4.99
|  style="border-bottom: 1pt solid black;text-align: center;vertical-align: top;"|479.54
|}


En la Figura 12 se puede observar que el tiempo de calentamiento requerido optimizado en el centro de la superficie de la cavidad de 12.15 s representa una reducción de 92.59 % en comparación con el tiempo de calentamiento requerido inicial de 23.4 s.

La Figura 13 muestra la variación de la temperatura a lo largo de la superficie del molde se mejora aproximadamente 14.50 % del valor de referencia de 5.714 °C al valor optimizado de 4.999 °C, y el esfuerzo máximo de Von Mises disminuye 4.05 % con respecto al valor inicial de 499 MPa, quedando en 479.54 MPa, por debajo del valor predefinido en la optimización mediante PSO.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image12.png|center|482px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 12'''. Respuestas de temperatura en el centro de la superficie de la cavidad en etapa de calentamiento
|}


{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image13.png|center|482px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 13'''. Distribución de la temperatura en la superficie de la cavidad entre el caso de diseño óptimo y el inicial
|}


Se realizó la simulación respuesta térmica y el análisis térmico-estructural con el parámetro optimizado para validar los resultados de optimización. Los resultados de la simulación muestran que el tiempo de calentamiento requerido, la variación de la temperatura superficial del molde y los esfuerzos máximos de Von Mises son 12.70 s, 4.762 °C y 500 MPa, respectivamente.

Los valores en porcentaje de errores entre los valores optimizados y experimentales de <math display="inline">{t}_{cr}</math>, <math display="inline">{t}_{d}</math> y <math display="inline">{\sigma }_{v}</math> son 4.25%, 4.78 %, 4.09 %, respectivamente.

Las gráficas de contorno de la distribución de la temperatura y la distribución de esfuerzos de Von Mises basada en los parámetros óptimos se muestran en las Figura 14 y 15, respectivamente.

{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image14.png|center|494px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| '''Figura 14'''. Gráfica de contorno de distribución de temperaturas con valores optimizados
|}


{| style="text-align: center; border: 1px solid #BBB; margin: 1em auto; width: auto;max-width: auto;"
|-
|style="padding:10px;"| [[Image:Review_590513318777-image15.png|center|506px]]
|- style="text-align: center; font-size: 75%;"
| colspan="1" style="padding:10px;"| <span id='_Ref510439590'></span>'''Figura 15'''. Esfuerzos de Von Mises con valores optimizados
|}

==7. Conclusiones==

El moldeo por respuesta térmica rápida es una nueva tecnología de moldeo por inyección desarrollada en los últimos años. El ciclo de moldeo es relativamente largo debido a la baja eficiencia de calentamiento, el bajo rendimiento de producción es debido a la distribución de temperatura no uniforme en la superficie de la cavidad y vida de servicio relativamente corta debido a la fatiga termo-mecánica son tres preocupaciones principales para el MRTR. Teniendo en cuenta tantos factores y limitaciones en el diseño de MRTR, implica la aplicación de análisis y optimización computarizados avanzados para la optimización del diseño de MRTR.

En este estudio se investigó la respuesta térmica y la resistencia estructural del MRTR por análisis térmico y estructural. Los resultados del análisis térmico muestran que el sistema de calentamiento del molde de inyección tiene la capacidad potencial para ser calentado rápida y uniformemente usando resistencias eléctricas. Los resultados del análisis de campos acoplados térmico-estructural indican que la distribución de los esfuerzos en el bloque de la cavidad es muy desigual. Existen concentraciones de esfuerzos en los bordes de las resistencias eléctricas cerca de la superficie de la cavidad del molde. Tal concentración de esfuerzo y relajación durante el proceso de calentamiento y enfriamiento frecuente en el MRTR puede conducir a la iniciación de la microfisura en el orificio donde se posiciona la resistencia eléctrica y se podría extender hasta la fatiga térmico-mecánica después de ciertos ciclos de moldeo.

Por lo tanto, para el sistema de calentamiento del molde, el espesor de la placa y la posición de las resistencias eléctricas son muy cruciales para el MRTR y debe ser cuidadosamente diseñado considerando la eficiencia de la respuesta térmica, la uniformidad de temperatura, la resistencia estructural del MRTR. La metodología de superficie de respuesta (MSR) basada en el diseño experimental de Box-Behnken y el experimento numérico se utilizó para explorar la relación numérica entre los parámetros de diseño y las variables de respuesta. Se implementó el método PSO para acoplar los modelos matemáticos desarrollados y así encontrar las variables de diseño óptimas para minimizar el tiempo de calentamiento requerido y, al mismo tiempo, asegurar uniformidad de temperatura aceptable y resistencia estructural al molde. Se obtuvo la posición óptima de las resistencias eléctricas  del molde para probetas de tensión, se predijeron y verificaron los valores mínimos de <math display="inline">{t}_{cr}</math>, <math display="inline">{t}_{d}</math> y <math display="inline">{\sigma }_{v}\,</math> realizando experimentos de confirmación. Después de la optimización, el tiempo de calentamiento requerido se redujo en 92.59 %, la uniformidad de la temperatura superficial de la cavidad se mejoró en 14.50 %, y la resistencia estructural del molde se mejoró ligeramente 4.05 %. Los modelos de superficie de respuesta cuadrática desarrollados pueden obtener resultados razonables con errores de menos de 5 % entre los datos predichos y los experimentales en experimentos de confirmación. Estos modelos pueden utilizarse para describir la relación entre factores experimentales (<math>h, a, q</math> y <math>b</math>) y variables de respuesta (<math display="inline">{t}_{cr}</math>, <math display="inline">{t}_{d}</math> y <math display="inline">{\sigma }_{v}</math>).

El método de optimización desarrollado por el acoplamiento MSR y PSO ha demostrado ser eficaz y será una poderosa herramienta para guiar el diseño de optimización de la posición de las resistencias eléctricas del MRTR, reduciendo así el tiempo de ciclo de moldeo, mejorando el rendimiento de producción y la vida útil del MRTR.

==Agradecimientos==

Los autores agradecen al Consejo Nacional de Ciencia y Tecnología (CONACyT) por la beca doctoral No. 428731.

Agradecemos a CoreTech System Co., Ltd. (Moldex3D) por su generoso suministro de licencias para esta investigación a través de su distribuido en el Estado de Guanajuato Lic. Salvador López de 3DCAD.

==Referencias==

{| style="width: 100%;" 
|-
|  style="vertical-align: top;"|[1] 
|  style="vertical-align: top;"|Kutz M., Applied plastics engineering handbook. Processing, materials and aplications, Second ed., Elsevier, Ed., William Andrew, Oxford, 2017.
|-
|  style="vertical-align: top;"|[2] 
|  style="vertical-align: top;"|Kennedy P. y Zheng R., Flow analysis of injection molds., Second ed., Henser,Munich, 2013.
|-
|  style="vertical-align: top;"|[3] 
|  style="vertical-align: top;"|Shoemaker J., Moldflow design guide: a resource for plastics engineers, First Edition ed., Hanser Verlag, Massachusetts, 2006. 
|-
|  style="vertical-align: top;"|[4] 
|  style="vertical-align: top;"|Wang W., Zhao G., Guan Y., Wu X. and Hui Y., Effect of rapid heating cycle injection mold temperature on crystal structures, morphology of polypropylene and surface quality of plastic parts, Journal of Polymer Research, 22.5:84, 2015. 
|-
|  style="vertical-align: top;"|[5] 
|  style="vertical-align: top;"|Chen S. C., Jong W. R. and Chang J. A., Dynamic mold surface temperature control using induction heating and its effects on the surface appearance of weld line, Journal of Applied Polymer Science, 101.2:1174-1180, 2006. 
|-
|  style="vertical-align: top;"|[6] 
|  style="vertical-align: top;"|Chen S. C., Minh P. S. and Chang J. A., Gas-assisted mold temperature control for improving the quality of injection molded parts with fiber additives, International Communications in Heat and Mass Transfer, 38.3:304-312, 2011. 
|-
|  style="vertical-align: top;"|[7] 
|  style="vertical-align: top;"|Yao D., Kimerling T. E.  and Kim B., High‐frequency proximity heating for injection molding applications, Polymer Engineering & Science, 46.7:938-945, 2006. 
|-
|  style="vertical-align: top;"|[8] 
|  style="vertical-align: top;"|Chang P. C. and Hwang S. J., Experimental investigation of infrared rapid surface heating for injection molding, Journal of applied polymer science, 102.4:3704-3713, 2006. 
|-
|  style="vertical-align: top;"|[9] 
|  style="vertical-align: top;"|Wang G., Zhao G., Li H. and Guan Y., Research on a new variotherm injection molding technology and its application on the molding of a large LCD panel, Polymer-Plastics Technology and Engineering, 48.7:671-681, 2009. 
|-
|  style="vertical-align: top;"|[10] 
|  style="vertical-align: top;"|Xu R. X. and Sachs E., Rapid thermal cycling with low thermal inertia tools, Polymer Engineering & Science, 49.2:305-316, 2009. 
|-
|  style="vertical-align: top;"|[11] 
|  style="vertical-align: top;"|Collomb, J., Balland, P., Francescato, P., Gardet, Y., Leh, D., & Saffré, P. Thermo-mechanical simulations for the comparison of heating channels geometries for composite material injection tools. Applied Thermal Engineering, 126:204-211, 2017.
|-
|  style="vertical-align: top;"|[12] 
|  style="vertical-align: top;"|Shayfull, Z., Sharif, S., Zain, A. M., Saad, R. M., & Fairuz, M. A. Milled groove square shape conformal cooling channels in injection molding process. Materials and Manufacturing Processes, 28.8:884-891, 2013.
|-
|  style="vertical-align: top;"|[13] 
|  style="vertical-align: top;"|Zhao G., Wang G., Guan Y.  and Li H., Research and application of a new rapid heat cycle molding with electric heating and coolant cooling to improve the surface quality of large LCD TV panels, Polymers for Advanced Technologies, 22.5:476-487, 2011. 
|-
|  style="vertical-align: top;"|[14] 
|  style="vertical-align: top;"|Xiao C. L. and Huang H. X., Multiobjective optimization design of heating system in electric heating rapid thermal cycling mold for yielding high gloss parts, Journal of Applied Polymer Science, 131.6, 2014. 
|-
|  style="vertical-align: top;"|[15] 
|  style="vertical-align: top;"|Xiao C. L.  and Huang H. X., Development of a rapid thermal cycling molding with electric heating and water impingement cooling for injection molding applications, Applied Thermal Engineering, 73.1:712-722, 2014. 
|-
|  style="vertical-align: top;"|[16] 
|  style="vertical-align: top;"|Li X. P., Zhao G. Q., Guan Y. J.  and Ma M. X., Optimal design of heating channels for rapid heating cycle injection mold based on response surface and genetic algorithm, Materials & Design, vol. 30, nº 10, pp. 4317-4323, 2009. 
|-
|  style="vertical-align: top;"|[17] 
|  style="vertical-align: top;"|Huang H. X., Li J. C., Li D. and Huang G. Q., New strategies for predicting parison dimensions in extrusion blow molding, Polymer-Plastics Technology and Engineering, vol. 50, nº 13, pp. 1329-1337, 2011. 
|-
|  style="vertical-align: top;"|[18] 
|  style="vertical-align: top;"|Shen C., Wang L.  and Li Q., Optimization of injection molding process parameters using combination of artificial neural network and genetic algorithm method, Journal of Materials Processing Technology, vol. 183, nº 2, pp. 412-418, 2007. 
|-
|  style="vertical-align: top;"|[19] 
|  style="vertical-align: top;"|Hassan, R., Cohanim, B., De Weck, O., & Venter, G. A comparison of particle swarm optimization and the genetic algorithm. In 46th AIAA/ASME/ASCE/AHS/ASC structures, structural dynamics and materials conference, p. 1897, 2005, April.
|-
|  style="vertical-align: top;"|[20] 
|  style="vertical-align: top;"|Park, J. B., Lee, K. S., Shin, J. R., & Lee, K. Y. A particle swarm optimization for economic dispatch with nonsmooth cost functions. IEEE Transactions on Power systems, vol. 20, nº 1, pp. 34-42, 2005. 
|-
|  style="vertical-align: top;"|[21] 
|  style="vertical-align: top;"|Settles, M., Rodebaugh, B., & Soule, T. Comparison of genetic algorithm and particle swarm optimizer when evolving a recurrent neural network. In Genetic and Evolutionary Computation Conference. Springer, Berlin, Heidelberg, pp. 148-149, 2003, July. 
|-
|  style="vertical-align: top;"|[22] 
|  style="vertical-align: top;"|Eberhart, R. C., & Shi, Y. Comparison between genetic algorithms and particle swarm optimization. In International conference on evolutionary programming. Springer, Berlin, Heidelberg, pp. 611-616, 1998, March.
|-
|  style="vertical-align: top;"|[23] 
|  style="vertical-align: top;"|ASTM, ASTM D638, ASTM Standard test method for tensile properties of plastics, in D638-03., West Conshohocken: ASTM International., 2003. 
|-
|  style="vertical-align: top;"|[24] 
|  style="vertical-align: top;"|Menges G., Michaeli W., and Mohren P., How to Make Injection Molds. 3rd edn. Hanser Gardner Pubns. pp 328, 2001.
|-
|  style="vertical-align: top;"|[25] 
|  style="vertical-align: top;"|Rao, N. S., & Schumacher, G. Design formulas for plastics engineers. Carl Hanser Verlag GmbH Co KG. pp 161, 2014.
|-
|  style="vertical-align: top;"|[26] 
|  style="vertical-align: top;"|Yao D., Chen S. C. and Kim B. H., Rapid thermal cycling of injection molds: an overview on technical approaches and applications, Advances in Polymer Technology, vol. 27, nº 4, pp. 233-255, 2008. 
|-
|  style="vertical-align: top;"|[27] 
|  style="vertical-align: top;"|Park, S. J., & Kwon, T. H. Optimal cooling system design for the injection molding process. Polymer Engineering & Science, vol. 38, nº 9, pp. 1450-1462, 1998.
|-
|  style="vertical-align: top;"|[28] 
|  style="vertical-align: top;"|Rao, N. S., Schumacher, G., Schott, N. R., & O’brien, K. T. Optimization of cooling systems in injection molds by an easily applicable analytical model. Journal of reinforced plastics and composites, vol. 21, nº 5, pp. 451-459, 2002.
|-
|  style="vertical-align: top;"|[29] 
|  style="vertical-align: top;"|Xiao C. L. and Huang H. X., Development of a rapid thermal cycling molding with electric heating and water impingement cooling for injection molding applications, Applied Thermal Engineering, vol. 73, nº 1, pp. 712-722, 2014. 
|-
|  style="vertical-align: top;"|[30] 
|  style="vertical-align: top;"|Cengel Y.  and Boles M., Thermodynamics: An Engineering Approach, 4th ed., New York: McGraw-Hill, 2002. 
|-
|  style="vertical-align: top;"|[31] 
|  style="vertical-align: top;"|Dorigato A., D'Amato M.  and Pegoretti A., Thermo-mechanical properties of high density polyethylene-fumed silica nanocomposites: effect of filler surface area treatment, Journal of Polymer Research, vol. 19, nº 6, p. 9889, 2012. 
|-
|  style="vertical-align: top;"|[32] 
|  style="vertical-align: top;"|Al Ashraf A., Thermal Analysis of Polymer by DSC, Center for Advanced Materials, vol. 2, p. 5, 2015. 
|-
|  style="vertical-align: top;"|[33] 
|  style="vertical-align: top;"|Montgomery D. C., Introduction to statistical quality control, Seventh ed., John Wiley & Sons, 2009. 
|-
|  style="vertical-align: top;"|[34] 
|  style="vertical-align: top;"|<span id='_GoBack'></span>ASTM, ASTM D3641, ASTM Standard Practice for Injection Molding Test Specimens of Thermoplastic Molding and Extrusion Materials, in ASTM D3641-15, West Conshohocken: ASTM INTERNATIONAL, 2003. 
|-
|  style="vertical-align: top;"|[35] 
|  style="vertical-align: top;"|ASTM, ASTM D4703, ASTM Standard practice for compression molding thermoplastic materials into test specimens, plaques, or sheets, in D4703-10a., West Conshohocken: ASTM International., 2010. 
|}