A large amount of low-rise unreinforced masonry structures (URM) can be foundin Groningen, the Netherlands. More and more induced earthquakes with shortduration have been detected in this region due to gas exploitation. Local unreinforcedmasonry (URM) buildings were initially not designed for withstanding seismicactions, so that unexpected damage may occur due to their vulnerability, raising insecurityamong residents. Existing low-rise masonry buildings in Groningen can bedivided into different categories based on their characteristics. Two types of residentialmasonry buildings that fulfil the prerequisites for performing non-linear seismicassessment are chosen to be studied in this thesis project, including the terracedhouse and the detached house.The seismic assessment of structures requires the use of both a discretization methodand a seismic assessment method. The discretization method is used to translate themechanical model into a finite element model used for the numerical analysis. Severalmethods have previously shown to be applicable for seismic assessment, but thiswork investigates the implications of using a continuum model (CM) and an equivalentframe model (EFM) approach to discretization in the general-purpose finiteelement package described in DIANA-FEA-BV (2017). The continuum model approachadopted was in a previous work by Schreppers et al. (2017) validated againstexperimental results and is as such deemed representative of the physical behaviourof the mechanical models investigated. An equivalent frame model approach to beused with DIANA is proposed in the work by Nobel (2017). The continuum modelapproach uses continuum elements with a constitutive model developed for the seismicassessment of masonry structures. This constitutive model captures both shearand flexural failure mechanisms. The equivalent frame model approach uses a combinationof numerically integrated beam elements and nodal interfaces, each witha distinct constitutive model, thus decoupling the description of the flexural andshear behaviour. This approach aims to capture the macro-behaviour at the structurallevel. The applicability of the proposed equivalent frame model approach isevaluated by how well it replicates the validated continuum model approach results.The two discretization methods described are evaluated using two types of seismicassessment methods. The first seismic assessment method used consists of first performinga quasi-static non-linear pushover analysis (NLPO) on the model. Thisresults in the pushover curve, which describes the global behaviour of the modelunder an equivalent lateral load based on the fundamental mode shape of the structure.The pushover curve is then used with the N2-method described in EN1998-1iii(2004) to assess at which peak ground acceleration (PGA) that the model reachesthe near-collapse (NC) limit state. The second seismic assessment method consistsof performing dynamic non-linear time-history analyses (NLTH). This method usesrecorded accelerograms to impose the inertial forces. The PGA for the accelerogramwhere the near-collapse limit state is reached is compared to the PGA fromthe use of the N2-method. The applicability of the pushover analysis in conjunctionwith the N2-method is evaluated by how well it replicates the PGA found from thetime-history analyses and by how well it replicates local failure mechanisms.Therefore, the main objectives of this project can be described by the following twoquestions:i. To what extent can the equivalent frame method be applicable as a properdiscretization method for pushover analyses and time-history analyses oflow-rise unreinforced masonry residential buildings in the Groningen region?ii. To what extent can the non-linear pushover method be adopted toassess the seismic behaviour of low-rise unreinforced masonry residentialbuildings in the Groningen region?The applicability of the equivalent frame model showed to vary. For describing localfailure mechanisms its applicability is poor. Further work on connecting the edgepiers to transverse walls is needed. For seismic assessment using the N2-method theapplicability of the equivalent frame model approach is sensible. The conservativedisplacement capacity counteracts the fact that it is worse at describing local unloading,which produced a larger initial equivalent stiffness of the bi-linear curvesin comparison to the continuum model. For seismic assessment using the timehistorysignals, its applicability is possible. While it could show different behaviourin terms of displacement and damping forces, it still showed a similar PGA at thenear-collapse limit state for the cases at hand.The seismic assessment of the terraced and detached houses by the N2-method issimilar to the seismic prediction by applying time-history analyses. However, thereare still some variations in the initial stiffness, force capacity and displacement capacitybetween these two assessment methods due to the assumptions and limitationsin this study. Overall, considering the pros and cons of the quasi-static pushovermethod, it is deemed applicable during the seismic assessment of the unreinforcedmasonry structures in the Groningen area.

This dissertation advances tools for the certification of unmanned aircraft system (UAS) flight controllers. We develop two thrusts to this goal: (1) the validation and improvement of an uncertain UAS framework based on integral quadratic constraint (IQC) theory and (2) the development of novel IQC theorems which allow the analysis of uncertain systems having time-varying characteristics. Pertaining to the first thrust, this work improves and implements an IQC-based robustness analysis framework for UAS. The approach models the UAS using a linear fractional transformation on uncertainties and conducts robustness analysis on the uncertain system via IQC theory. By expressing the set of desired UAS flight paths with an uncertainty, the framework enables analysis of the uncertain UAS flying about any level path whose radius of curvature is bounded. To demonstrate the versatility of this technique, we use IQC analysis to tune trajectory-tracking and path-following controllers designed via H2 or H-infinity synthesis methods. IQC analysis is also used to tune path-following PID controllers. By employing a non-deterministic simulation environment and conducting numerous flight tests, we demonstrate the capability of the framework in predicting loss of control, comparing the robustness of different controllers, and tuning controllers. Finally, this work demonstrates that signal IQCs have an important role in obtaining IQC analysis results which are less conservative and more consistent with observations from flight test data. With regards to the second thrust, we prove a novel theorem which enables robustness analysis of uncertain systems where the nominal plant and the IQC multiplier are linear time-varying systems and the nominal plant may have a non-zero initial condition. When the nominal plant and the IQC multiplier are eventually periodic, robustness analysis can be accomplished by solving a finite-dimensional semidefinite program. Time-varying IQC multipliers are beneficial in analysis because they provide the possibility of reducing conservatism and are capable of expressing uncertainties that have unique time-domain characteristics. A number of time-varying IQC multipliers are introduced to better describe such uncertainties. The utility of this theorem is demonstrated with various examples, including one which produces bounds on the UAS position after an aggressive Split-S maneuver. / Doctor of Philosophy / This work develops tools to aid in the certification of unmanned aircraft system (UAS) flight controllers. The forthcoming results are founded on robust control theory, which allows the incorporation of a variety of uncertainties in the UAS mathematical model and provides tools to determine how robust the system is to these uncertainties. Such a foundation provides a complementary perspective to that obtained with simulations. Whereas simulation environments provide a probabilistic-type analysis and are oftentimes costly, the following results provide worst-case guarantees—for the allowable disturbances and uncertainties—and require far less computational resources. Here we take two approaches in our development of certification tools for UAS. First we validate and improve on an uncertain UAS framework that relies on integral quadratic constraint (IQC) theory to analyze the robustness of the UAS in the presence of uncertainties and disturbances. Our second approach develops novel IQC theorems that can aid in providing bounds on the UAS state during its flight trajectory. Though the applications in this dissertation are focused on UAS, the theory can be applied to a wide variety of physical and nonphysical problems wherein uncertainties in the mathematical model cannot be avoided.

L'application d'Abel-Jacobi fait le lien entre la forme de Weierstrass d'une courbe elliptique définie sur C et le tore complexe qui lui est associé. Il est possible de la calculer en un nombre d'opérations quasi-linéaire en la précision voulue, c'est à dire en temps O(M(P) log P). Son inverse est donné par la fonction p de Weierstrass, qui s'exprime en fonction de thêta, une fonction importante en théorie des nombres. L'algorithme naturel d'évaluation de thêta nécessite O(M(P) sqrt(P)) opérations, mais certaines valeurs (les thêta-constantes) peuvent être calculées en O(M(P) log P) opérations en exploitant les liens avec la moyenne arithmético-géométrique (AGM). Dans ce manuscrit, nous généralisons cet algorithme afin de calculer thêta en O(M(P) log P). Nous exhibons une fonction F qui a des propriétés similaires à l'AGM. D'une façon similaire à l'algorithme pour les thêta-constantes, nous pouvons alors utiliser la méthode de Newton pour calculer la valeur de thêta. Nous avons implanté cet algorithme, qui est plus rapide que la méthode naïve pour des précisions supérieures à 300 000 chiffres décimaux. Nous montrons comment généraliser cet algorithme en genre supérieur, et en particulier comment généraliser la fonction F. En genre 2, nous sommes parvenus à prouver que la même méthode mène à un algorithme qui évalue thêta en O(M(P) log P) opérations ; la même complexité s'applique aussi à l'application d'Abel-Jacobi. Cet algorithme est plus rapide que la méthode naïve pour des précisions plus faibles qu'en genre 1, de l'ordre de 3 000 chiffres décimaux. Nous esquissons également des pistes pour obtenir la même complexité en genre quelconque. Enfin, nous exhibons un nouvel algorithme permettant de calculer une isogénie de courbes elliptiques de noyau donné. Cet algorithme utilise l'application d'Abel-Jacobi, car il est facile d'évaluer l'isogénie sur le tore ; il est sans doute possible de le généraliser au genre supérieur / The Abel-Jacobi map links the short Weierstrass form of a complex elliptic curve to the complex torus associated to it. One can compute it with a number of operations which is quasi-linear in the target precision, i.e. in time O(M(P) log P). Its inverse is given by Weierstrass's p-function, which can be written as a function of theta, an important function in number theory. The natural algorithm for evaluating theta requires O(M(P) sqrt(P)) operations, but some values (the theta-constants) can be computed in O(M(P) log P) operations by exploiting the links with the arithmetico-geometric mean (AGM). In this manuscript, we generalize this algorithm in order to compute theta in O(M(P) log P). We give a function F which has similar properties to the AGM. As with the algorithm for theta-constants, we can then use Newton's method to compute the value of theta. We implemented this algorithm, which is faster than the naive method for precisions larger than 300,000 decimal digits. We then study the generalization of this algorithm in higher genus, and in particular how to generalize the F function. In genus 2, we managed to prove that the same method leads to a O(M(P) log P) algorithm for theta; the same complexity applies to the Abel-Jacobi map. This algorithm is faster than the naive method for precisions smaller than in genus 1, of about 3,000 decimal digits. We also outline a way one could reach the same complexity in any genus. Finally, we study a new algorithm which computes an isogeny of elliptic curves with given kernel. This algorithm uses the Abel-Jacobi map because it is easy to evaluate the isogeny on the complex torus; this algorithm may be generalizable to higher genera

時間數列分析過程，常常發現其走勢，隨著時間過程而演變，應用傳統的線性模式來配適，往往很難獲得合適預測值。因此近幾年來，非線性時間數列結構性改變的研究越來越受到重視，也一直是時間數列及計量經濟學學者所熱衷的研究主題之一。本文利用模糊理論的觀念，以模糊炳找出配適ARCH模式數列之轉折區間，分別以轉折區問起始點及結束點作為截取點，去配適ARCH(1)模式，稱之為截取式自迴歸條件變異數分析法(Trimmed ARCH(1) model)。針對台幣對美元銀行間每日收盤匯率，分別以單變量ARIMA、ARCH(1)、Trtmmed ARCH(1)來建構模式，並做比較分析。比較結果發現，以轉折區間結束點作為截取點之Trimmed ARCH(1)模式，其預測值最為準確，大為改善了原來ARCH(1)模式之預測水準。 / In time series analysis, we often find the trend of which changing with time. Using the traditional model fitting can't get a good prediction. Hence the research of structure change of non-linear time series is attentive in recent years, and non-linear time series analysis is a research topic which the scholars of time series and econometrics are intent on. This article tries to use the theory of fuzzy ,to recognize the structure change period by the fuzzy classification, let the first point and the last point of the structure change period be the cute points, to fit ARCH(1) mod ie which we called the Trimmed ARCH(1) model. We use the data of the exchange rate between N.T dol liars and U.S dollars to compare the ARIMAwith ARCH(1) and Trimmed ARCH(1), the forcasting performance shows that Trimmed ARCH(1) model takes a better prediction result.

近幾年來，非線性時間數列分析一直是時間數列及計量經濟學者所熱衷的研究主題之一，而非線性時間數列結構改變的研究也越來越受到重視。其中的門檻自迴歸模式，雖具有線性模式所不能配適的特性，但模式建構的問題，一直是其在發展應用上的瓶頸。本研究擬以門檻自迴歸模式建構的流程並結合遺傳演算法的最佳化搜尋技術，架構出時間數列遺傳演算法，藉此演算法則及程序，全域性地搜尋最佳的門檻自迴歸模式。 / Non-linear time series analysis is a research topic which the schalors of time series and econometrics are intent on, and the research of structure change of non-linear time series is attentive. Threshold autoregressive model (TAR model) of non-linear time series has some characters which linear model fail to fit while the problem of how to find an appropriate threshold value is still attracted many researchers attention. In this paper, we present about searching the parameters for a TAR model by genetic algorithms.

當我們利用Shewhart管制圖(Shewhart control chart)或累積和管制圖(Cumulative-sum chart. CUSUM chart)來偵測製程時,通常假設製品係獨立取自一個服從均數μ和標準差為σ的獨立常態分配的管制下進行。但是若產品特性值呈現自相關時,這類管制圖就可能發生誤導的結果。本文利用時間序列模式來解決具相關變數的管制圖問題。並考慮利用非線性時間序列模式及特別原因管制圖(special-cause control chart)來檢視台灣經濟景氣指標是否處於控制中的狀態。並討論特別原因管制圖的連串長度分佈(run length distribution)。在最後的實例分析中,介紹自動控制的觀念。 / Traditionally, in the quality control process, such as: Shewhart control chart or CUSUM chart, it is assumed that the observation process follows an i.i.d normal distribution. If the assumption for independence fails, that is when the process exhibits type of autocorrelation, we need to find a more reliable decision method. In this paper, we will apply the time series analysis and structure changed concept to slove the serial correlation problem. The idea of automatic control can be applied in the explanation of this nonlinear process. Finally, a time series about the monitoring indicators of Taiwan is discussed in detail as an example.

While moment matching can effectively reduce the dimension of a linear, time-invariant system, it can simultaneously fail to improve the stable time-step for the forward Euler scheme. In the context of a semi-discrete heat equation with spatially smooth forcing, the high frequency modes are virtually insignificant. Eliminating such modes dramatically improves the stable time-step without sacrificing output accuracy. This is accomplished by modal filtration, whose computational cost is relatively palatable when applied following an initial reduction stage by moment matching. A bound on the norm of the difference between the transfer functions of the moment-matched system and its modally-filtered counterpart yields an intelligent choice for the mode of truncation. The dual-stage algorithm disappoints in the context of highly nonnormal semi-discrete convection-diffusion equations. There, moment matching can be ineffective in dimension reduction, precluding a cost-effective modal filtering step.

Due to the emergence of new materials and advanced engineering technology, slender footbridges are increasingly becoming popular to satisfy the modern transportation needs and the aesthetical requirements of society. These structures however are always &quotlively" with low stiffness, low mass, low damping and low natural frequencies. As a consequence, they are prone to vibration induced by human activities and can suffer severe vibration serviceability problems, particularly in the lateral direction. This phenomenon has been evidenced by the excessive lateral vibration of many footbridges worldwide such as the Millennium Bridge in London and the T-Bridge in Japan. Unfortunately, present bridge design codes worldwide do not provide sufficient guidelines and information to address such vibrations problems and to ensure safety and serviceability due to the lack of knowledge on the dynamic performance of such slender vibration sensitive bridge structures. A conceptual study has been carried out to comprehensively investigate the dynamic characteristics of slender suspension footbridges under human-induced dynamic loads and a footbridge model in full size with pre-tensioned reverse profiled cables in the vertical and horizontal planes has been proposed for this purpose. A similar physical suspension bridge model was designed and constructed in the laboratory, and experimental testings have been carried out to calibrate the computer simulations. The synchronous excitation induced by walking has been modelled as crowd walking dynamic loads which consist of dynamic vertical force, dynamic lateral force and static vertical force. The dynamic behaviour under synchronous excitation is simulated by resonant vibration at the pacing rate which coincides with a natural frequency of the footbridge structure. Two structural analysis software packages, Microstran and SAP2000 have been employed in the extensive numerical analysis. Research results show that the structural stiffness and vibration properties of suspension footbridges with pre-tensioned reverse profiled cables can be adjusted by choosing different structural parameters such as cable sag, cable section and pretensions in the reverse profiled cables. Slender suspension footbridges always have four main kinds of vibration modes: lateral, torsional, vertical and longitudinal modes. The lateral and torsional modes are often combined together and become two kinds of coupled modes: coupled lateral-torsional modes and coupled torsionallateral modes. Such kind of slender footbridges also have different dynamic performance in the lateral and vertical directions, and damping has only a small effect on the lateral vibration but significant effect on the vertical one. The fundamental coupled lateral-torsional mode and vertical mode are easily excited when crowd walking dynamic loads are distributed on full bridge deck. When the crowd walking dynamic loads are distributed eccentrically on half width of the deck, the fundamental coupled torsional-lateral mode can be excited and large lateral deflection can be induced. Higher order vertical modes and coupled lateral-torsional modes can also be excited by groups of walking pedestrians under certain conditions. It is found that the coupling coefficient introduced in this thesis to describe the coupling of a coupled mode, is an important factor which has significant effect on the lateral dynamic performance of slender suspension footbridges. The coupling coefficient, however, is influenced by many structural parameters such as cable configuration, cable section, cable sag, bridge span and pre-tensions, etc. In general, a large dynamic amplification factor is expected when the fundamental mode of a footbridge structure is the coupled lateral-torsional mode with a small coupling coefficient. The research findings of this thesis are useful in understanding the complex dynamic behaviour of slender and vibration sensitive suspension footbridges under humaninduced dynamic loads. They are also helpful in developing design guidance and techniques to improve the dynamic performance of such slender vibration sensitive footbridges and similar structures and hence to ensure their safety and serviceability.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-08-22T19:52:23Z No. of bitstreams: 1 ANA EMÍLIA VICTOR BARBOSA - DISSERTAÇÃO PPGCC 2007..pdf: 1669761 bytes, checksum: f47054507fe9200c8d1d56d2848ae276 (MD5) / Made available in DSpace on 2018-08-22T19:52:23Z (GMT). No. of bitstreams: 1 ANA EMÍLIA VICTOR BARBOSA - DISSERTAÇÃO PPGCC 2007..pdf: 1669761 bytes, checksum: f47054507fe9200c8d1d56d2848ae276 (MD5) Previous issue date: 2007-04-20 / Capes / Neste trabalho propomos uma técnica que visa detectar violações de propriedades comportamentais automaticamente durante a execução de sistema de software concorrentes. A técnica foi inspirada na metodologia de desenvolvimento Design by Contract (DbC). DbC permite que os desenvolvedores adicionem aos programas asserções para que sejam verificadas em tempo de execução. O uso de asserções para expressar propriedades de programas concorrentes (multithreaded)eparalelos, entretanto,não ésuficiente. Nesses sistemas,muitas das propriedades comportamentais de interesse, como vivacidade e segurança, não podem ser expressas apenas com asserções. Essas propriedades requerem o uso de operadores temporais. Neste trabalho, utilizamos Lógica Linear Temporal (Linear Time Logic - LTL) para expressar o comportamento desejado. Para dar suporte a checagem do comportamento dos programas em tempo de execução, propomos uma técnica baseada em Programação Orientada a Aspectos, que permite que o programa seja continuamente monitorado (o comportamento é checado através do uso de autômatos que permite a deteção de comportamentos inesperados). Associada a cada propriedade comportamental existe um conjunto de pontos de interesse do código-fonte que devem obedece-la. Esses pontos são então monitorados durante a execução do sistema através do uso de aspectos. Entre outros benefícios, a técnica permite que o sistema de software alvo seja instrumentado de maneira não intrusiva, sem alterar o código-fonte — particulamente, nenhum código do software alvo deve ser modificado para execução da monitoração. Para validar este trabalho, desenvolvemos como prova de conceitos um protótipo que implementa a técnica e permite a monitoração de programas Java multi-threaded, chamado DesignMonitor. Essa ferramenta é apresentada e discutida através de um estudo de caso para demonstrar a aplicação da técnica / In this work we propose and develop a technique that allows to detect the violation of behavior properties of concurrent systems. The technique was inspired by the Design by Contract (DbC) programming methodology, which proposes the use of assertions and their evaluation at runtime to check programs behavior. The use of simple assertions to express properties of concurrent and parallel programs, however, is not sufficient. Many of the relevant properties of those systems,s uch as liveness and security, can not be expressed with simple assertions. Thesepropertiesrequiretheuseof temporal operators. In our work, we used Linear Time Logic (LTL) to specify the expected behavior. To support the runtime checking of the program against the expected behavior, we propose a technique, based on Aspect-Oriented Programming, that allows the program to be continuously monitored (behavior is checked against automata that allows the detection of unexpected behaviors). Each property is mapped to a set of points of interest in the target program. Those points are then monitored during the system execution through aspects. Among other benefits, the technique allows the instrumentation of the target software to be performed automatically and in a non-intrusive way — in particular, no code must be changed toturn monitoring on or off. To validate the work, we developed a proof of concept prototype tool that implements the technique and allows the monitoring of multi-threaded Java programs, called DesignMonitor. The tool was used in case study that has allowed the evaluation and the discussion of practical issues related with the technique.

Orientador: Pedro Luis Dias Peres / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Eletrica e de Computação / Made available in DSpace on 2018-08-05T03:23:54Z (GMT). No. of bitstreams: 1 Leite_ValterJuniordeSouza_D.pdf: 4310851 bytes, checksum: bcd0414d19eb46e02496290857fdf9bc (MD5) Previous issue date: 2005 / Resumo: Este trabalho trata da aplica¸c¿ao de funcionais de Lyapunov e Lyapunov-Krasovskii dependentes de parâmetro a alguns problemas selecionados da área de controle robusto, a saber: D-estabilidade robusta de polítipo de matrizes, D-estabilidade robusta de politopos de polinômios matriciais, estabilidade robusta de sistemas neutrais com atrasos variantes no tempo e controle robusto H8 de sistemas discretos no tempo com atraso nos estados. É utilizada a representação politópica para as incertezas dos sistemas estudados. São obtidas formulações convexas, na forma de desigualdades matriciais lineares, suficientes para a solução dos problemas selecionados. Essas condições podem ser resolvidas numericamente de maneira eficiente por meio de algoritmos especializados baseados em pontos interiores. Os resultados obtidos são menos conservadores que os encontrados na literatura, baseados em geral na estabilidade quadrática, isto é, as matrizes dos funcionais são fixas e independentes da incerteza / Abstract: This work deals with the application of parameter dependent Lyapunov and Lyapunov-Krasovskii functionals to some selected problems of robust control: robust D-stability of polytopes of matrices, robust D-stability of polytopes of polynomial matrices, robust stability of uncertain neutral systems with timevarying delays and robust H8 control of uncertain discrete time delay systems. The polytopic representation is used to describe the uncertainties. Convex formulations are obtained, in terms of inear matrix inequalities, that are sufficient for the solution of the selected problems. Those conditions can be solved in a efficient way through specialized interior point algorithms. The obtained results are less conservative than those from the literature, in general based on quadratic stability, i.e., the matrices in the functionals are fixed and do not depend on the uncertainty / Doutorado / Automação / Doutor em Engenharia Elétrica