The development of self-powered wireless sensor networks for structural and machinery health monitoring has attracted considerable attention in the research field during the last decade. Since the low-duty-cycle wireless sensor networks have significantly reduced the power requirements to the range of tens to hundreds of microwatts, it is possible to harvest environmental energy as the power supply instead of using batteries. Vibration energy harvesting using piezoelectric materials has become the most popular technique, which has a good potential to generate adequate power. However, there is a limitation for the conventional beam-shaped harvester designs in real applications due to their limited bandwidth. In order to overcome this limitation, the essential objective of this thesis is to develop harvesters with multi-resonance structures. The multi-resonance harvester with good broadband performance can achieve close resonance frequencies and relatively large power output in each vibration mode. The main tasks and contributions of this thesis are summarised as follows: • A parametric analysis is presented to determine how the modal structural and electromechanical performances of cantilevered beam harvesters are affected by two modal factors designated as mass ratio and electromechanical coupling coefficient (EMCC). The modal performance of using rectangular, convergent and divergent tapered configurations with and without extra masses are systematically analysed by geometric variation using the finite element analysis (FEA) software ABAQUS. • A modal approach using the two modal factors to evaluate the modal performance of harvesters is introduced and a configurational optimization strategy based on the modal approach is developed to pre-select the configurations of multi-resonance harvesters with better modal structural performance and close resonance frequencies in multiple modes. Using this optimization strategy obviates the need to run the full analysis at the first stage. • A novel two-layer stacked harvester, which consists of a base cantilevered beam that is connected to an upper beam by a rigid mass, is developed. By altering the dimensions and the locations of the masses, the two-layer harvester can generate two close resonance frequencies with relatively large power output. The effects of using rectangular, convergent and divergent tapered beam configurations are systematically analysed. • Multi-layer stacked harvesters with up to five layers are developed. The three-layer harvesters with different mass positions, which can generate three close resonance frequencies, are optimized using the configurational optimization strategy. • A novel doubly-clamped multi-layer harvester, which is able to generate five close resonance frequencies with relatively large power output, is developed and thoroughly analysed. • An experimental study of the multi-layer stacked harvester is presented to validate the simulated results and the configurational optimization strategy. • An experimental study of the two-layer stacked harvester using high performance single crystal piezoelectric material PIMNT is presented. The harvester using PIMNT can generate nearly 10 times larger power output and 3.5 times wider bandwidth than using PZT. Besides, by modifying the location of the piezoelectric layer, anti-resonances between two adjacent modes can be eliminated.

A supercavitating vehicle, a next-generation underwater vehicle capable of changing the paradigm of modern marine warfare, exploits supercavitation as a means to reduce drag and achieve extremely high submerged speeds. In supercavitating flows, a low-density gaseous cavity entirely envelops the vehicle and as a result the vehicle is in contact with liquid water only at its nose and partially over the afterbody. Hence, the vehicle experiences a substantially reduced skin drag and can achieve much higher speed than conventional vehicles. The development of a controllable and maneuvering supercavitating vehicle has been confronted with various challenging problems such as the potential instability of the vehicle, the unsteady nature of cavity dynamics, the complex and non-linear nature of the interaction between vehicle and cavity. Furthermore, major questions still need to be resolved regarding the basic configuration of the vehicle itself, including its control surfaces, the control system, and the cavity dynamics. In order to answer these fundamental questions, together with many similar ones, this dissertation develops an integrated simulation-based design tool to optimize the vehicle configuration subjected to operational design requirements, while predicting the complex coupled behavior of the vehicle for each design configuration. Particularly, this research attempts to include maneuvering flight as well as various operating trim conditions directly in the vehicle configurational optimization. This integrated approach provides significant improvement in performance in the preliminary design phase and indicates that trade-offs between various performance indexes are required due to their conflicting requirements. This dissertation also investigates trim conditions and dynamic characteristics of supercavitating vehicles through a full 6 DOF model. The influence of operating conditions, and cavity models and their memory effects on trim is analyzed and discussed. Unique characteristics are identified, e.g. the cavity memory effects introduce a favorable stabilizing effect by providing restoring fins and planing forces. Furthermore, this research investigates the flight envelope of a supercavitating vehicle, which is significantly different from that of a conventional vehicle due to different hydrodynamic coefficients as well as unique operational conditions.

<div>As physical testing does not always yield insight into the mechanistic cause of failures, computational modeling is often used to develop an understanding of the goodness of a design and to shorten the product development time. One common, and widely used analysis technique is the Finite Element Method. A significant difficulty with the finite element method is the effort required to generate an analysis-suitable mesh due to the difference in the mathematical representation of geometry CAD and CAE systems. CAD systems commonly use Non-Uniform Rational B-Splines (NURBS) while the CAE tools rely on the finite element mesh. Efforts to unify CAD and CAE by carrying out analysis directly using NURBS models termed Isogeometric Analysis reduces the gap between CAD and CAE phases of product development. However, several challenges still remain in the field of isogeometric analysis. A critical challenge relates to the output of commercial CAD systems. B-rep CAD models generated by commercial CAD systems contain uncoupled NURBS patches and are therefore not suitable for analysis directly. Existing literature is largely missing methods to smoothly couple NURBS patches. This is the first topic of research in this thesis. Fracture-caused failures are a critical concern for the reliability of engineered structures in general and semiconductor chips in particular. The back-end of the line structures in modern semiconductor chips contain multi-material junctions that are sites of singular stress, and locations where cracks originate during fabrication or testing. Techniques to accurately model the singular stress fields at interfacial corners are relatively limited. This is the second topic addressed in this thesis. Thus, the overall objective of this dissertation is to develop an isogeometric framework for parametric domain decomposition and analysis of singular stresses using enriched isogeometric analysis.</div><div><br></div><div>Geometrically speaking, multi-material junctions, sub-domain interfaces and crack surfaces are lower-dimensional features relative to the two- or three-dimensional domain. The enriched isogeometric analysis described in this research builds enriching approximations directly on the lower-dimensional geometric features that then couple sub-domains or describe cracks. Since the interface or crack geometry is explicitly represented, it is easy to apply boundary conditions in a strong sense and to directly calculate geometric quantities such as normals or curvatures at any point on the geometry. These advantages contrast against those of implicit geometry methods including level set or phase-field methods. In the enriched isogeometric analysis, the base approximations in the domain/subdomains are enriched by the interfacial fields constructed as a function of distance from the interfaces. To circumvent the challenges of measuring distance and point of influence from the interface using iterative operations, algebraic level sets and algebraic point projection are utilized. The developed techniques are implemented as a program in the MATLAB environment named as <i>Hierarchical Design and Analysis Code</i>. The code is carefully designed to ensure simplicity and maintainability, to facilitate geometry creation, pre-processing, analysis and post-processing with optimal efficiency. </div><div><br></div><div>To couple NURBS patches, a parametric stitching strategy that assures arbitrary smoothness across subdomains with non-matching discretization is developed. The key concept used to accomplish the coupling is the insertion of a “parametric stitching” or p-stitching interface between the incompatible patches. In the present work, NURBS is chosen for discretizing the parametric subdomains. The developed procedure though is valid for other representations of subdomains whose basis functions obey partition of unity. The proposed method is validated through patch tests from which near-optimal rate of convergence is demonstrated. Several two- and three-dimensional elastostatic as well as heat conduction numerical examples are presented.</div><div><br></div><div>An enriched field approximation is then developed for characterizing stress singularities at junctions of general multi-material corners including crack tips. Using enriched isogeometric analysis, the developed method explicitly tracks the singular points and interfaces embedded in a non-conforming mesh. Solution convergence to those of linear elastic fracture mechanics is verified through several examples. More importantly, the proposed method enables direct extraction of generalized stress intensity factors upon solution of the problems without the need to use <i>a posteriori</i> path-independent integral such as the J-integral. Next, the analysis of crack initiation and propagation is carried out using the alternative concept of configurational force. The configurational force is first shown to result from a configurational optimization problem, which yields a configurational derivative as a necessary condition. For specific velocities imposed on the heterogeneities corresponding to translation, rotation or scaling, the configurational derivative is shown to yield the configurational force. The use of configurational force to analyze crack propagation is demonstrated through examples.</div><div><br></div><div>The developed methods are lastly applied to investigate the risk of ratcheting-induced fracture in the back end of line structure during thermal cycle test of a epoxy molded microelectronic package. The first principal stress and the opening mode stress intensity factor are proposed as the failure descriptors. A finite element analysis sub-modeling and load decomposition procedure is proposed to study the accumulation of plastic deformation in the metal line and to identify the critical loading mode. Enriched isogeometric analysis with singular stress enrichment is carried out to identify the interfacial corners most vulnerable to stress concentration and crack initiation. Correlation is made between the failure descriptors and the design parameters of the structure. Crack path from the identified critical corner is predicted using both linear elastic fracture mechanics criterion and configurational force criterion. </div>