Familia automatizada y convergencia

Paso 7 — Limpiar sin destruir la geometría

/PREP7
ALLSEL,ALL
DDELE,ALL,ALL
FDELE,ALL,ALL
CMDELE,fixed_nodes
CMDELE,tip_nodes
VCLEAR,ALL

VCLEAR,ALL elimina los elementos y nodos asociados al volumen, pero conserva el modelo sólido y sus atributos. Cambiar ESIZE sin limpiar y volver a mallar no refina una malla existente: solo cambia una preferencia para operaciones futuras.

Paso 8 — Automatizar la familia

*DIM,mesh_size,ARRAY,3
mesh_size(1)=0.0500
mesh_size(2)=0.0250
mesh_size(3)=0.0125

*DO,i,1,3
  ! divisiones, mallado, auditoria, medida y limpieza
*ENDDO

El bucle no crea tres estudios independientes escritos a mano. Obliga a que todos recorran exactamente la misma secuencia y reduce las oportunidades de cambiar una condición sin querer.

Paso 9 — Medir la convergencia

change_pct(1)=ABS(uy_tip(2)-uy_tip(1))/ABS(uy_tip(2))*100
change_pct(2)=ABS(uy_tip(3)-uy_tip(2))/ABS(uy_tip(3))*100
change_pct(3)=-1

Cada malla se compara con la siguiente más fina, utilizada como estimación mejorada. La última no puede compararse consigo misma; el valor -1 significa «sin comparación», no convergencia perfecta.

Paso 10 — Tomar una decisión honesta

El algoritmo recorre las comparaciones desde la malla más económica:

selected_case=0

! Si cambio 1→2 < 2 % y ambas mallas pasan:
!   selected_case=1
! En caso contrario, probar cambio 2→3.
! Si tampoco cumple:
!   selected_case=0

Nunca elegimos automáticamente la malla 3 por ser la más fina disponible. Si la última comparación supera el umbral, el estudio es inconcluso y necesita otra malla.

Por qué no usamos la tensión máxima

El empotramiento ideal introduce una transición abrupta entre nodos completamente restringidos y material libre. Una tensión puntual junto a esa frontera puede crecer o desplazarse al refinar, incluso cuando el desplazamiento global ya es estable. Elegirla como único criterio confundiría convergencia de una magnitud global con comportamiento local potencialmente singular.

En M07 compararemos tensiones en regiones justificables y separaremos tensión nominal, concentración física y singularidad numérica.

El CSV como acta del estudio

04_mesh_convergence.mac genera m04_mesh_study.csv:

case,mesh_h,nx,ny,nz,n_nodes,n_elements,n_type1,n_mat1,uy_tip,change_to_next_pct,mesh_pass,selected

Debe contener tres filas. Solo una puede tener selected=1; también es válido que todas sean cero si ninguna comparación cumple. El CSV registra intención, coste, auditoría, respuesta y decisión en una misma evidencia reproducible.