Implementation of a partitioned algorithm for simulation of large CSI problems

The implementation of a partitioned numerical algorithm for determining the dynamic response of coupled structure/controller/estimator finite-dimensional systems is reviewed. The partitioned approach leads to a set of coupled first and second-order linear differential equations which are numerically integrated with extrapolation and implicit step methods. The present software implementation, ACSIS, utilizes parallel processing techniques at various levels to optimize performance on a shared-memory concurrent/vector processing system. A general procedure for the design of controller and filter gains is also implemented, which utilizes the vibration characteristics of the structure to be solved. Also presented are: example problems; a user's guide to the software; the procedures and algorithm scripts; a stability analysis for the algorithm; and the source code for the parallel implementation

Analysis, preliminary design and simulation systems for control-structure interaction problems

Software aspects of control-structure interaction (CSI) analysis are discussed. The following subject areas are covered: (1) implementation of a partitioned algorithm for simulation of large CSI problems; (2) second-order discrete Kalman filtering equations for CSI simulations; and (3) parallel computations and control of adaptive structures

Second-order discrete Kalman filtering equations for control-structure interaction simulations

A general form for the first-order representation of the continuous, second-order linear structural dynamics equations is introduced in order to derive a corresponding form of first-order Kalman filtering equations (KFE). Time integration of the resulting first-order KFE is carried out via a set of linear multistep integration formulas. It is shown that a judicious combined selection of computational paths and the undetermined matrices introduced in the general form of the first-order linear structural systems leads to a class of second-order discrete KFE involving only symmetric, N x N solution matrix

Parallel computations and control of adaptive structures

The equations of motion for structures with adaptive elements for vibration control are presented for parallel computations to be used as a software package for real-time control of flexible space structures. A brief introduction of the state-of-the-art parallel computational capability is also presented. Time marching strategies are developed for an effective use of massive parallel mapping, partitioning, and the necessary arithmetic operations. An example is offered for the simulation of control-structure interaction on a parallel computer and the impact of the approach presented for applications in other disciplines than aerospace industry is assessed

Discrete Kalman filtering equations of second-order form for control-structure interaction simulations

A second-order form of discrete Kalman filtering equations is proposed as a candidate state estimator for efficient simulations of control-structure interactions in coupled physical coordinate configurations as opposed to decoupled modal coordinates. The resulting matrix equation of the present state estimator consists of the same symmetric, sparse N x N coupled matrices of the governing structural dynamics equations as opposed to unsymmetric 2N x 2N state space-based estimators. Thus, in addition to substantial computational efficiency improvement, the present estimator can be applied to control-structure design optimization for which the physical coordinates associated with the mass, damping and stiffness matrices of the structure are needed instead of modal coordinates

Finding diamonds in the rough: Targeted Sub-threshold Search for Strongly-lensed Gravitational-wave Events

Strong gravitational lensing of gravitational waves can produce duplicate signals separated in time with different amplitudes. We consider the case in which strong lensing produces identifiable gravitational-wave events and weaker sub-threshold signals hidden in the noise background. We present a search method for the sub-threshold signals using reduced template banks targeting specific confirmed gravitational-wave events. We apply the method to all events from Advanced LIGO's first and second observing run O1/O2. Using GW150914 as an example, we show that the method effectively reduces the noise background and raises the significance of (near-) sub-threshold triggers. In the case of GW150914, we can improve the sensitive distance by $2.0\% - 14.8\%$. Finally, we present the top $5$ possible lensed candidates for O1/O2 gravitational-wave events that passed our nominal significance threshold of False-Alarm-Rate $\leq 1/30$ days

Extraction of Substructural Flexibility from Global Frequencies and Mode Shapes

A computational procedure for extracting substructure-by-substructure flexibility properties from global modal parameters is presented. The present procedure consists of two key features: an element-based direct flexibility method which uniquely determines the global flexibility without resorting to case-dependent redundancy selections; and, the projection of cinematically inadmissible modes that are contained in the iterated substructural matrices. The direct flexibility method is used as the basis of an inverse problem, whose goal is to determine substructural flexibilities given the global flexibility, geometrically-determined substructural rigid-body modes, and the local-to-global assembly operators. The resulting procedure, given accurate global flexibility, extracts the exact element-by-element substructural flexibilities for determinate structures. For indeterminate structures, the accuracy depends on the iteration tolerance limits. The procedure is illustrated using both simple and complex numerical examples, and appears to be effective for structural applications such as damage localization and finite element model reconciliation

Targeted Sub-threshold Search for Strongly-lensed Gravitational-wave Events

Strong gravitational lensing of gravitational waves can produce duplicated signals that are separated in time and with different amplitudes. We consider the case in which strong lensing produces identifiable gravitational-wave events together with weaker sub-threshold signals that are hidden in the noise background. We present a search method for the sub-threshold signals using reduced template banks targeting specific confirmed gravitational-wave events. We apply the method to an event from Advanced LIGO's first observing run O1, GW151012. We show that the method is effective in reducing the noise background and hence raising the significance of (near-) sub-threshold triggers. In the case of GW151012, we are able to improve the sensitive distance by 10%−25%. Finally, we present the 10 most significant events for GW151012-like signals in O1. Besides the already confirmed gravitational-wave detections, none of the candidates pass our nominal significance threshold of False-Alarm-Rate ≤ 1/30 days