本論文は,様々なネットワーク最適化問題の数学的構造について様々な観点から調べたものである.主たる結果はネットワーク最適化問題の代表例である最大流問題に,関するいくつかの結果と,Min-Plus代数に値をもつ行列の固有値と固有ベクトルに関する特徴づけに関する結果からなっている. / 博士(理学) / Doctor of Philosophy in Science / 同志社大学 / Doshisha University

[pt] Neste trabalho apresentam-se três tipos de técnicas de análise do colapso estrutural através do método dos elementos finitos: análise linearizada da carga crítica, análise incremental da carga crítica e análise não linear completa. Na análise linearizada da carga crítica formulou-se um problema de autovalor empregando matrizes de rigidez baseadas na configuração indeformada da estrutura e materiais com comportamento linear elástico. No caso da análise incremental da carga crítica, o problema de autovalor foi formulado empregando matrizes de rigidez incrementais para levar em consideração os grandes deslocamentos e propriedades não lineares do material. Finalmente, na análise não linear completa a configuração deformada da estrutura e propriedades não lineares do material são atualizadas durante todo o processo incremental-iterativo até atingir a carga crítica. Desenvolveu-se uma implementação computacional para estudar as três técnicas de análise em estruturas planas como vigas, colunas, pórticos e arcos, empregando elementos isoparamétricos bidimensionais para estado plano de tensões. A configuração deformada da estrutura, devido aos grandes deslocamentos e rotações dos elementos, foi considerada através de uma formulação Lagrangeana Total, enquanto o comportamento inelástico do material foi modelado empregando um modelo elastoplástico de Von Mises (J2) com encruamento isotrópico. Nos exemplos apresentados mostrou-se a influência da não linearidade geométrica e física na estimativa de cargas críticas e no comportamento pós-crítico, podendo ocorrer bifurcações ao longo da trajetória de equilíbrio fundamental definida no espaço carga-deslocamentos. / [en] This work presents three kinds of techniques for collapse analysis using the finite element method: linear buckling analysis, nonlinear buckling analysis and full nonlinear analysis. The linear buckling analysis requires the definition of an eigenvalue problem using a stiffness matrix formulation based on the initial configuration of the structure and under the assumption of a linear elastic material behavior. In the case of nonlinear buckling analysis, the eigenvalue problem was formulated employing an incremental stiffness matrix in order to consider the effects of large displacements and nonlinear material properties in the critical load estimation. Finally, the full nonlinear analysis takes into account the deformed configuration and the nonlinear material properties of the structure, updating both of them through all the incremental-iterative process up to reaching the critical load. A Finite Element computational program, using plane stress isoperimetric bidimensional elements, was developed to study the three analysis techniques applied to plane structures such as beams, columns, frames and arches. The deformed configuration of the structure, due to large displacements and rotations, was considered through the Total Lagrangian formulation, whereas the inelastic material behavior was modeled using the Von Mises plasticity model with isotropic hardening. The examples presented in this article show the influence of geometric and material nonlinearity in the critical load estimation and the postcritical behavior, being this the reason for the potential occurrence of bifurcation points over the fundamental equilibrium path defined in the load-displacement space.

For noiseless sampled data, we describe the close connections between Prony--like methods, namely the classical Prony method, the matrix pencil method and the ESPRIT method. Further we present a new efficient algorithm of matrix pencil factorization based on QR decomposition of a rectangular Hankel matrix. The algorithms of parameter estimation are also applied to sparse Fourier approximation and nonlinear approximation.

El análisis y el control de las vibraciones cobra especial importancia en muchas ramas de la ingeniería, en especial la ingeniería mecánica, civil, aeronáutica y automovilística. Tal es así que prácticamente se identi¿ca como un área independiente dentro del análisis dinámico de estructuras. Desde los comienzos de esta teoría, las fuerzas disipativas o de amortiguamiento han sido uno de los fenómenos más difíciles de modelizar. El modelo viscoso, por su sencillez y versatilidad ha sido y sigue siendo el gran paradigma de los modelos de amortiguamiento. Sin embargo, como consecuencia de la aparición de materiales con memoria se introdujo el fenómeno de la viscoelasticidad; Esta, si bien está también 'íntimamente ligada ' a la velocidad de la respuesta, necesito de la introducción de las denominadas funciones hereditarias, que permiten poner a las fuerzas disipativas como función no solo de la velocidad instantánea sino de la historia de velocidades desde el comienzo del movimiento, de ahí el termino memoria. De forma natural, el avance teórico introducido en el modelo supone también una complicación computacional, pues donde antes teníamos un sistema lineal de ecuaciones diferenciales ahora tenemos un sistema de ecuaciones integro-diferenciales. El análisis de las vibraciones libres de los sistemas con amortiguamiento viscoelástico conduce a un problema nolineal de autovalores donde la característica principal es una matriz de amortiguamiento que depende de la frecuencia de excitación. El estudio de la solución de autovalores y autovectores de este problema es importante si se desean conocer los modos de vibración de la estructura o si se pretende obtener la respuesta en el dominio de la frecuencia del sistema. El objetivo fundamental de esta Tesis Doctoral es doble: Por un lado, profundizar en el conocimiento del problema de autovalores de sistemas viscoelásticos proponiendo para ello nuevos métodos numéricos de resolución. Por otro, desarrollar un nuevo modelo viscoso que, bajo ciertas condiciones, reproduzca la respuesta del modelo viscoelástico con su¿ciente aproximación. La Tesis se divide en ocho capítulos, de ellos el cuerpo principal se encuentra en los seis centrales (Capítulos 2 a 7. Todos ellos son artículos de investigación que, o bien han sido publicados, o bien están en proceso de revisión en revistas contenidas en el Journal Citation Reports (JCR). Por esta razón, todos los capítulos conservan la estructura intrínseca de un artículo, incluidas una introducción y una bibliografía en cada uno. Los cuatro primeros capítulos (Capítulos 2 a 5) se centran en el estudio del problema no lineal de autovalores. Se proponen dos metodologías de resolución: la primera es un procedimiento iterativo basado en el esquema del punto-¿jo y desarrollado para sistemas proporcionales o ligeramente no-proporcionales (aquellos en los que los modos se presentan desacoplados o casi desacoplados). La segunda metodología (presentada en dos capítulos diferentes), denominada paramétrica, permite obtener soluciones casi-analíticas de los autovalores, tanto para sistemas de un grado de libertad como para sistemas de múltiples grados de libertad y dentro de 'estos, para sistemas proporcionales y no proporcionales. El estudio del problema de autovalores se completa con un capítulo dedicado a los autovalores reales, también denominados autovalores no viscosos. En 'él se demuestra una nueva caracterización maten ática que deben cumplir dichos autovalores y que permite proponer un nuevo concepto: el conjunto no-viscoso. Los dos 'últimos capítulos (Capítulos 6 y 7) analizan el Modelo Viscoso Equivalente como propuesta para la modelización de la respuesta de sistemas viscoelásticos. El análisis se realiza desde el dominio de la frecuencia estudiando la función de transferencia. En una primera etapa (pen último capítulo), de naturaleza más maten ática, se demuestra que la función de transferencia exacta de un modelo viscoelástico se puede expresar como suma de una función de transferencia propia de un modelo viscoso más un término denominado residual, directamente dependiente del nivel de amortiguamiento inducido y del acoplamiento modal (noproporcionalidad de la matriz de amortiguamiento). En una segunda etapa ('ultimo capítulo), se desarrolla una aplicación para estructuras reales formadas por entramados planos de elementos 1D amortiguados con capas de material visco elástico. Este tipo de estructuras ha permitido usar una variante mejorada del método paramétrico para la obtención de los autovalores, de forma que en este 'ultimo capítulo ha servido como nexo de unión de las metodologías más importantes desarrolladas en la Tesis. / Lázaro Navarro, M. (2013). The Eigenvalue Problem in Linear Viscoelastic Structures: New Numerical Approaches and the Equivalent Viscous Model [Tesis doctoral no publicada]. Universitat Politècnica de València. https://doi.org/10.4995/Thesis/10251/30062 / TESIS

For noiseless sampled data, we describe the close connections between Prony--like methods, namely the classical Prony method, the matrix pencil method and the ESPRIT method. Further we present a new efficient algorithm of matrix pencil factorization based on QR decomposition of a rectangular Hankel matrix. The algorithms of parameter estimation are also applied to sparse Fourier approximation and nonlinear approximation.

Tato dizertační práce se zabývá dvěma případy interakce tělesa s tekutinou (FSI). První z nich se zabývá analýzou vzájemné interakce mezi rotorem čerpadla a kapalinou uvnitř těsnící spáry. Vliv těsnící spáry na dynamiku celého stoje je popsán pomocí dynamických parametrů, které jsou také označovaný jako přídavné účinky. V současnosti používané modely těsnících spár používají pro stanovení dynamických parametrů řadu zjednodušujících předpokladů. V této práci je prezentováno pět různých analýz dynamických parametrů těsnící spáry čerpadla na okysličovadlo. Každá z těchto pěti analýz používá jinou míru zjednodušení výpočetního modelu. V případě největšího zjednodušení je modelován pouze objem kapaliny uvnitř těsnící spáry. Nejkomplexnější analýza pro stanovení dynamických parametrů těsnící spáry používá pro výpočet model celého čerpadla s excentrickou polohou rotoru. Druhá část této dizertační práce definuje novou metodu pro řešení interakce kapaliny s pružným tělesem. Tato metoda využívá řešení inverzního problému kmitání. Přímý problém kmitání, který je také označován jako problém vlastních hodnot, používá jako vstupy pro řešení matice hmotnosti, tuhosti a tlumení, které jsou dohromady označovány jako koeficientové matice, na základě kterých je v nejobecnějším případě stanovena Jordanovská matice a také modální matice pravostranných a levostranných vlastních vektorů. Při řešení inverzního problému kmitání jsou stanoveny koeficientové matice na základě Jordanovské matice a modálních matic pravostranných a levostranných vlastních vektorů. Existují dva případy inverzního problému kmitání. V případě, že jsou známy všechny vstupní vlastní čísla a vlastní vektory, pak se jedná o tzv. plný problém. Naopak v případě, že alespoň 1 mód kmitání soustavy není znám, tak se jedná o tzv. částečný problém. V této práci je prezentováno 5 algoritmů pro řešení inverzního problému v kmitání. Nicméně pro každý typ inverzního problému kmitání je prezentován jeden univerzální algoritmus. Algoritmus pro řešení plných problémů byl poprvé prezentován v roce 1979 Otakarem Daňkem. Algoritmy pro řešení částečných problémů, které jsou prezentovány v této práci, jsou vůbec prvními algoritmy pro řešení tohoto typu inverzního problému kmitání. Univerzální algoritmus pro řešení částečných problémů je označován jako algoritmus pro řešení částečných problémů s volbou doplňkových vlastních hodnot. Aplikace těchto dvou univerzálních algoritmů pro řešení inverzního problému kmitání pro případ plných i částečných problémů je ukázána na řešení dvou případů interakce pružného tělesa s kapalinou.

Pensionsmyndigheten's (the Swedish Pensions Agency) current calculation of the internal rate of return for 7.7 million premium pension savers is both time and resource consuming. This rate of return mirrors the overall performance of the funded part of the pension system and is analyzed internally, but also reported to the public monthly and yearly based on differently sized data samples. This thesis aims to investigate the possibility of utilizing other approaches in order to improve the performance of these calculations. Further, the study aims to verify the results stemming from said calculations and investigate their robustness. In order to investigate competitive matrix methods, a sample of approaches are compared to the more classical numerical methods. The approaches are compared in different scenarios aimed to mirror real practice. The robustness of the results are then analyzed by a stochastic modeling approach, where a small error term is introduced aimed to mimic possible errors which could arise in data management. It is concluded that a combination of Halley's method and the Jacobi-Davidson algorithm is the most robust and high performing method. The proposed method combines the speed and robustness from numerical and matrix methods, respectively. The result show a performance improvement of 550% in time, while maintaining the accuracy of the current server computations. The analysis of error propagation suggests the output error to be less than 0.12 percentage points in 99 percent of the cases, considering an introduced error term of large proportions. In this extreme case, the modeled expected number of individuals with an error exceeding 1 percentage point is estimated to be 212 out of the whole population. / Pensionsmyndighetens nuvarande beräkning av internräntan för 7,7 miljoner pensionssparare är både tid- och resurskrävande. Denna avkastning ger en översikt av hur väl den fonderade delen av pensionssystemet fungerar. Detta analyseras internt men rapporteras även till allmänheten varje månad samt årligen baserat på olika urval av data. Denna uppsats avser att undersöka möjligheten att använda andra tillvägagångssätt för att förbättra prestanda för denna typ av beräkningar. Vidare syftar studien till att verifiera resultaten som härrör från dessa beräkningar och undersöka deras stabilitet. För att undersöka om det finns konkurrerande matrismetoder jämförs ett urval av tillvägagångssätt med de mer klassiska numeriska metoderna. Metoderna jämförs i flera olika scenarier som syftar till att spegla verklig praxis. Stabiliteten i resultaten analyseras med en stokastisk modellering där en felterm införs för att efterlikna möjliga fel som kan uppstå i datahantering. Man drar slutsatsen att en kombination av Halleys metod och Jacobi-Davidson-algoritmen är den mest robusta och högpresterande metoden. Den föreslagna metoden kombinerar hastigheten från numeriska metoder och tillförlitlighet från matrismetoder. Resultatet visar en prestandaförbättring på 550 % i tid, samtidigt som samma noggrannhet som ses i de befintliga serverberäkningarna bibehålls. Analysen av felutbredning föreslår att felet i 99 procent av fallen är mindre än 0,12 procentenheter i det fall där införd felterm har stora proportioner. I detta extrema fall uppskattas det förväntade antalet individer med ett fel som överstiger 1 procentenhet vara 212 av hela befolkningen.

Garantierte untere Eigenwertschranken (GLB) für elliptische Eigenwertprobleme partieller Differentialgleichungen sind in der Theorie sowie in praktischen Anwendungen relevant. Auf Grund des Rayleigh-Ritz- (oder) min-max-Prinzips berechnen alle konformen Finite-Elemente-Methoden (FEM) garantierte obere Schranken. Ein Postprocessing nichtkonformer Methoden von Carstensen und Gedicke (Math. Comp., 83.290, 2014) sowie Carstensen und Gallistl (Numer. Math., 126.1, 2014) berechnet GLB. In diesen Schranken ist die maximale Netzweite ein globaler Parameter, das kann bei adaptiver Netzverfeinerung zu deutlichen Unterschätzungen führen. In einigen numerischen Beispielen versagt dieses Postprocessing für lokal verfeinerte Netze komplett. Diese Dissertation präsentiert, inspiriert von einer neuen skeletal-Methode von Carstensen, Zhai und Zhang (SIAM J. Numer. Anal., 58.1, 2020), einerseits eine modifizierte hybrid-high-order Methode (m=1) und andererseits ein allgemeines Framework für extra-stabilisierte nichtkonforme Crouzeix-Raviart (m=1) bzw. Morley (m=2) FEM. Diese neuen Methoden berechnen direkte GLB für den m-Laplace-Operator, bei denen eine leicht überprüfbare Bedingung an die maximale Netzweite garantiert, dass der k-te diskrete Eigenwert eine untere Schranke für den k-ten Dirichlet-Eigenwert ist. Diese GLB-Eigenschaft und a priori Konvergenzraten werden für jede Raumdimension etabliert. Der neu entwickelte Ansatz erlaubt adaptive Netzverfeinerung, die für optimale Konvergenzraten auch bei nichtglatten Eigenfunktionen erforderlich ist. Die Überlegenheit der neuen adaptiven FEM wird durch eine Vielzahl repräsentativer numerischer Beispiele illustriert. Für die extra-stabilisierte GLB wird bewiesen, dass sie mit optimalen Raten gegen einen einfachen Eigenwert konvergiert, indem die Axiome der Adaptivität von Carstensen, Feischl, Page und Praetorius (Comput. Math. Appl., 67.6, 2014) sowie Carstensen und Rabus (SIAM J. Numer. Anal., 55.6, 2017) verallgemeinert werden. / Guaranteed lower eigenvalue bounds (GLB) for elliptic eigenvalue problems of partial differential equation are of high relevance in theory and praxis. Due to the Rayleigh-Ritz (or) min-max principle all conforming finite element methods (FEM) provide guaranteed upper eigenvalue bounds. A post-processing for nonconforming FEM of Carstensen and Gedicke (Math. Comp., 83.290, 2014) as well as Carstensen and Gallistl (Numer. Math., 126.1,2014) computes GLB. However, the maximal mesh-size enters as a global parameter in the eigenvalue bound and may cause significant underestimation for adaptive mesh-refinement. There are numerical examples, where this post-processing on locally refined meshes fails completely. Inspired by a recent skeletal method from Carstensen, Zhai, and Zhang (SIAM J. Numer. Anal., 58.1, 2020) this thesis presents on the one hand a modified hybrid high-order method (m=1) and on the other hand a general framework for an extra-stabilized nonconforming Crouzeix-Raviart (m=1) or Morley (m=2) FEM. These novel methods compute direct GLB for the m-Laplace operator in that a specific smallness assumption on the maximal mesh-size guarantees that the computed k-th discrete eigenvalue is a lower bound for the k-th Dirichlet eigenvalue. This GLB property as well as a priori convergence rates are established in any space dimension. The novel ansatz allows for adaptive mesh-refinement necessary to recover optimal convergence rates for non-smooth eigenfunctions. Striking numerical evidence indicates the superiority of the new adaptive eigensolvers. For the extra-stabilized nonconforming methods (a generalization of) known abstract arguments entitled as the axioms of adaptivity from Carstensen, Feischl, Page, and Praetorius (Comput. Math. Appl., 67.6, 2014) as well as Carstensen and Rabus (SIAM J. Numer. Anal., 55.6, 2017) allow to prove the convergence of the GLB towards a simple eigenvalue with optimal rates.