In this dissertation, a novel fast modeling technique called peri-ultrasound that can model both linear and nonlinear ultrasonic behavior of materials is developed and implemented. Nonlinear ultrasonic response can detect even very small material non- linearity. Quantification of the material nonlinearity at the early stages of damage is important to avoid catastrophic failure and reduce repair costs. The developed model uses the nonlocal continuum-based peridynamic theory which was found to be a good simulation tool for handling crack propagation modeling, in particular when multiple cracks grow simultaneously. The developed peri-ultrasound modeling tool has been used to model the ultrasonic response at the interface of two materials in presence of an interface crack. Also, the stress wave propagation in a half-space (or half-plane for a 2-dimensional problem) with boundary loading is investigated using peri-ultrasound modeling. In another simulation, well-established two-dimensional Lamb's problem is investigated where the results are verified against available analytical solution. Also, the interaction between the surface wave and a surface breaking crack is studied.