Optimización Topológica Isogeométrica Basada en Derivadas Topológicas

Información Básica

• ID del Artículo: 2509.09236
• Título: Optimización Topológica Isogeométrica Basada en Derivadas Topológicas
• Autores: Guilherme Henrique TEIXEIRA (Universidad Tecnológica de Graz), Nepomuk KRENN (Instituto Johann Radon), Peter GANGL (Instituto Johann Radon), Benjamin MARUSSIG (Universidad Tecnológica de Graz)
• Clasificación: math.NA cs.CE cs.NA math.OC
• Fecha de Publicación: 15 de octubre de 2025 (versión v3 de arXiv)
• Enlace del Artículo: https://arxiv.org/abs/2509.09236

Resumen

La optimización topológica es una herramienta importante en ingeniería para el diseño de estructuras optimizadas. Sin embargo, los cambios topológicos generalmente requieren remallado, lo que puede resultar desafiante. Este artículo propone un método de optimización topológica isogeométrica impulsado por derivadas topológicas. La combinación del método de conjuntos de nivel con un marco isogeométrico inmerso permite actualizaciones geométricas sin fisuras sin necesidad de remallado. Simultáneamente, las derivadas topológicas proporcionan modificaciones topológicas sin necesidad de definir agujeros iniciales. Se investigó el impacto de funciones base de orden superior en la representación de conjuntos de nivel y la aproximación de soluciones. Dos ejemplos numéricos validan el método propuesto, demostrando que la aproximación de soluciones con funciones base de orden superior mejora la precisión, mientras que las funciones base lineales siguen siendo suficientes para la representación de funciones de conjuntos de nivel.

Antecedentes de Investigación y Motivación

Definición del Problema

El problema central que aborda esta investigación es la dificultad de reconstrucción de mallas en la optimización topológica tradicional. Durante el proceso de optimización estructural, los cambios topológicos del dominio de diseño (como la generación, fusión o desaparición de agujeros) generalmente requieren la regeneración de la malla computacional, lo que no solo es costoso desde el punto de vista computacional, sino que también puede provocar inestabilidad numérica.

Importancia

La optimización topológica, como herramienta importante en el diseño estructural, se aplica ampliamente en aeronáutica, fabricación automotriz, ingeniería civil y otros campos. Resolver el problema de reconstrucción de mallas es significativo para mejorar la eficiencia de optimización y garantizar la estabilidad numérica.

Limitaciones de Métodos Existentes

1. Método de Conjuntos de Nivel Tradicional: Aunque puede evitar la reconstrucción explícita de mallas, es difícil generar automáticamente nuevos agujeros y depende de la configuración de agujeros de la geometría inicial
2. Método de Densidad: Aunque evita problemas de malla, puede producir elementos grises y patrones de tablero de ajedrez
3. Método de Derivadas de Forma: Solo puede manejar cambios de límites, no cambios topológicos

Motivación de la Investigación

Los autores proponen combinar análisis isogeométrico (IGA), método de conjuntos de nivel y derivadas topológicas con el objetivo de:

• Eliminar la necesidad de reconstrucción de mallas
• Generar automáticamente nuevos agujeros sin necesidad de predefinición
• Aprovechar la suavidad de orden superior de B-splines para mejorar la precisión
• Simplificar el proceso de actualización geométrica

Contribuciones Principales

1. Marco de Método Innovador: Primera combinación del método puramente impulsado por derivadas topológicas de Amstutz y Andrä con análisis isogeométrico
2. Optimización sin Reconstrucción de Mallas: Implementación de actualizaciones geométricas sin fisuras durante cambios topológicos mediante un marco isogeométrico inmerso
3. Investigación de Funciones Base de Orden Superior: Estudio sistemático del impacto de diferentes órdenes polinomiales en la representación de conjuntos de nivel y aproximación de soluciones
4. Validación Numérica: Verificación de la efectividad del método mediante ejemplos clásicos, demostrando que funciones base lineales son suficientes para representar conjuntos de nivel mientras que funciones base de orden superior mejoran la precisión de la solución

Explicación Detallada del Método

Definición de la Tarea

Entrada: Dominio de diseño D, condiciones de límite, cargas, propiedades de materiales Salida: Distribución óptima de material Ω ⊂ D Objetivo: Minimizar la flexibilidad satisfaciendo restricciones de volumen Restricciones: Ecuaciones de elasticidad lineal, condiciones de límite

Arquitectura del Modelo

1. Formulación del Problema de Elasticidad Lineal

Resolver en dominio fijo D:

-∇ · (αΩσ(u)) = 0 en D
u = 0 en ΓD
σ(u) · n = τ en ΓN
σ(u) · n = 0 en Γ0

donde el parámetro de material αΩ se define como:

αΩ = {
  αin   en Ω (región de material)
  αout  en D\Ω (región de vacío, αout ≪ 1)
}

2. Representación de Conjuntos de Nivel

El dominio Ω se representa implícitamente mediante la función de conjunto de nivel ϕ:

ϕ(x) < 0 ⟺ x ∈ Ω
ϕ(x) = 0 ⟺ x ∈ ∂Ω  
ϕ(x) > 0 ⟺ x ∈ D\Ω

La función de conjunto de nivel se discretiza usando B-splines:

ϕ(ξ,η) = Σᵢⱼ Bᵢ,d(ξ)Bⱼ,d(η)cᵢⱼ

3. Cálculo de Derivadas Topológicas

Para el problema de elasticidad lineal, la expresión analítica de la derivada topológica es:

dℐ(Ω)(x₀) = {
  -3αin((αout-αin)/(2αout+αin))σ(u):ε(u) - l  si x₀ ∈ Ω
  -3αout((αin-αout)/(2αin+αout))σ(u):ε(u) + l  si x₀ ∈ D\Ω
}

4. Algoritmo de Actualización de Conjuntos de Nivel

Se utiliza interpolación lineal esférica para actualizar el conjunto de nivel:

ϕⁱ⁺¹ = (1/sin θᵢ)[sin((1-κᵢ)θᵢ)ϕⁱ + sin(κᵢθᵢ)gΩᵢ]

donde θᵢ es el ángulo L² entre el conjunto de nivel actual y la derivada topológica generalizada gΩᵢ.

Puntos de Innovación Técnica

1. Método Isogeométrico Inmerso

• Uso de una malla de fondo B-spline fija, evitando reconstrucción de mallas
• Distinción entre regiones de material y vacío mediante el parámetro de material αΩ
• Mantenimiento de la suavidad geométrica de orden superior

2. Impulso Puramente por Derivadas Topológicas

• No requiere resolver la ecuación de Hamilton-Jacobi, depende únicamente de derivadas topológicas
• Genera automáticamente nuevos agujeros sin necesidad de predefinición de agujeros iniciales
• Simplifica el proceso de optimización y mejora la eficiencia computacional

3. Tratamiento de Elementos Cortados

• Uso de biblioteca de integración geométrica implícita para cálculo preciso de proporciones cortadas
• Adopción de integración gaussiana estándar para evitar inestabilidad en integración de regiones pequeñas
• Implementación de filtrado promediado para suavizar saltos de propiedades de materiales

Configuración Experimental

Ejemplos Numéricos

Se seleccionaron dos problemas de referencia clásicos de optimización topológica:

1. Problema de Viga Voladiza: Dominio rectangular 2×1, extremo izquierdo fijo, carga puntual en extremo derecho
2. Problema de Cuarto de Anillo: Cuarto de anillo con radio interior 1 y radio exterior 2, fondo fijo, carga puntual en parte superior

Configuración de Parámetros

• Malla: Malla de fondo 128×128
• Parámetros de Material: αin = 1, αout = 10⁻⁴
• Control de Volumen: l = 5
• Parámetro de Filtrado: γ = 10⁻⁴
• Propiedades de Material: E = 1, ν = 1/3
• Condición de Convergencia: εθ = 1, máximo 200 iteraciones

Configuración de Comparación

Se investigaron dos configuraciones:

1. Orden Idéntico: Conjunto de nivel y solución con el mismo orden polinomial (d = p)
2. Orden Mixto: Conjunto de nivel lineal (d = 1) con solución de orden superior (p > 1)

Indicadores de Evaluación

• Valor de función objetivo J (flexibilidad)
• Ángulo de convergencia θ
• Área relativa Ai/A₀
• Número de iteraciones de convergencia

Resultados Experimentales

Resultados Principales del Problema de Viga Voladiza

Rendimiento de Convergencia

• Ventaja de Orden Superior: p ≥ 2 converge más rápidamente que p = 1, con descenso más rápido de la función objetivo en etapas intermedias
• Convergencia Consistente: Todas las configuraciones convergen a soluciones óptimas similares (J ≈ 6.97-7.00)
• Eficiencia de Iteración: p = 2 y p = 4 requieren el mismo número de iteraciones, p = 3 requiere 29 pasos con conjunto de nivel de orden superior vs 24 pasos con conjunto de nivel lineal

Impacto del Orden del Conjunto de Nivel

• El conjunto de nivel lineal (d = 1) produce formas finales similares al conjunto de nivel de orden superior (d = p)
• Diferencias de paridad: Los órdenes impares comparten puntos de Greville causando efecto de promediado, produciendo diferencias leves

Resultados Principales del Problema de Cuarto de Anillo

Fenómeno de Óptimo Local

• Diferencias de Convergencia: Configuraciones de orden inferior p = 1,2 se detienen en ángulos más altos (θ = 4.72, 9.33)
• Solución Óptima: p = 3(d = 1) alcanza el valor de función objetivo más bajo de 4.12, p = 3(d = 3) es 4.13
• Balance de Eficiencia: p = 2(d = 2) alcanza 4.15 en 129 iteraciones, con eficiencia relativamente alta

Diversidad de Formas

Diferentes configuraciones producen diferentes estructuras topológicas, reflejando los desafíos inherentes de óptimos locales en optimización topológica.

Hallazgos Clave

1. Aproximación de Solución de Orden Superior Beneficiosa: Mejora la velocidad de convergencia y precisión de la solución
2. Conjunto de Nivel Lineal Suficiente: d = 1 es suficiente para representar con precisión los límites geométricos
3. Independencia de Malla: En malla fina 128×128, la ventaja de orden superior es limitada
4. Importancia del Filtrado: El filtrado promediado y el filtrado PDE son críticos para la calidad de resultados

Trabajo Relacionado

Clasificación de Métodos de Optimización Topológica

1. Método de Densidad: SIMP, RAMP, etc., pueden producir elementos grises
2. Método de Conjuntos de Nivel: Evita elementos grises pero requiere agujeros iniciales
3. Método de Derivadas Topológicas: Puede generar automáticamente agujeros, adoptado en este artículo

Optimización Topológica Isogeométrica

• Trabajos Tempranos: Dedè et al. (2012) modelo de campo de fase, Shojaee et al. (2012) funciones de base radial
• Método de B-splines: Jahangiry y Tavakkoli (2017), Khatibinia et al. (2020)
• Innovación del Artículo: Primera combinación de derivadas topológicas puras con método isogeométrico

Método Inmerso

Ampliamente aplicado en mecánica de fluidos, mecánica de sólidos y otros campos, este artículo lo introduce en optimización topológica para lograr optimización sin reconstrucción de mallas.

Conclusiones y Discusión

Conclusiones Principales

1. Se desarrolló exitosamente un marco de optimización topológica isogeométrica inmersa impulsada por derivadas topológicas
2. Las funciones base de orden superior mejoran la aproximación de soluciones pero las funciones base lineales son suficientes para representar conjuntos de nivel
3. El método evita reconstrucción de mallas y puede generar automáticamente agujeros
4. Se verificó la efectividad del método en problemas de referencia clásicos

Limitaciones

1. Problema de Óptimo Local: Diferentes configuraciones de parámetros pueden converger a diferentes soluciones óptimas locales
2. Tratamiento de Saltos de Material: La ventaja de funciones base de orden superior en la representación de saltos en la interfaz material-vacío no es evidente
3. Sensibilidad de Parámetros: Parámetros como el parámetro de filtrado γ tienen impacto importante en los resultados
4. Limitación a 2D: Actualmente solo se verifica en problemas 2D

Direcciones Futuras

1. Investigación profunda de representación de orden superior de saltos en interfaz de material
2. Extensión a problemas tridimensionales
3. Optimización topológica multiaterial
4. Condiciones de restricción más complejas como restricciones de esfuerzo

Evaluación Profunda

Ventajas

1. Fuerte Innovación de Método: Combinación orgánica de tres tecnologías avanzadas para resolver problemas de ingeniería práctica
2. Base Teórica Sólida: Derivación rigurosa de derivadas topológicas, fórmulas matemáticas completas
3. Diseño Experimental Razonable: Investigación sistemática de diferentes configuraciones de orden, comparación exhaustiva
4. Alta Practicidad de Ingeniería: Evita reconstrucción de mallas, adecuado para aplicaciones prácticas

Insuficiencias

1. Ausencia de Extensión a 3D: Solo verificado en 2D, la mayoría de aplicaciones de ingeniería son 3D
2. Orientación de Parámetros Insuficiente: Falta análisis profundo sobre selección de parámetros clave como parámetro de filtrado
3. Falta de Análisis de Eficiencia Computacional: Sin comparación de costos computacionales con métodos tradicionales
4. Verificación de Geometría Compleja Insuficiente: Los ejemplos son relativamente simples, falta verificación de estructuras de ingeniería complejas

Impacto

1. Contribución Académica: Proporciona nueva ruta técnica para el campo de optimización topológica
2. Valor de Ingeniería: Potencial aplicación en diseño estructural de aeronáutica, automoción y otros campos
3. Reproducibilidad: Basado en código GeoPDEs de código abierto, facilita reproducción y extensión

Escenarios Aplicables

1. Optimización de Estructura de Continuo: Especialmente adecuado para problemas de diseño que requieren cambios topológicos frecuentes
2. Etapa de Diseño Conceptual: Exploración rápida de diferentes configuraciones topológicas
3. Aplicaciones de Requisitos de Alta Precisión: Escenarios que aprovechan la suavidad de orden superior de B-splines
4. Diseño Automatizado: Optimización automática sin necesidad de definición manual de agujeros iniciales

Referencias

El artículo cita 34 referencias importantes que cubren trabajos clásicos y de vanguardia en los campos centrales de optimización topológica, análisis isogeométrico y método de conjuntos de nivel, proporcionando una base teórica sólida para la investigación. Las referencias clave incluyen el método de derivadas topológicas de Amstutz y Andrä (2006), la teoría fundamental de análisis isogeométrico de Hughes et al. (2005), entre otros.