Research Article Unified Trajectory Planning Algorithms for Autonomous Underwater Vehicle Navigation

Similar documents
Stable Trajectory Design for Highly Constrained Environments using Receding Horizon Control

Research Article MFT-MAC: A Duty-Cycle MAC Protocol Using Multiframe Transmission for Wireless Sensor Networks

Optimal Path Finding for Direction, Location and Time Dependent Costs, with Application to Vessel Routing

Path Planning. Marcello Restelli. Dipartimento di Elettronica e Informazione Politecnico di Milano tel:

AUV Cruise Path Planning Based on Energy Priority and Current Model

Visibility Graph. How does a Mobile Robot get from A to B?

Path Planning and Decision-making Control for AUV with Complex Environment

Autonomous Mobile Robots, Chapter 6 Planning and Navigation Where am I going? How do I get there? Localization. Cognition. Real World Environment

Waypoint Navigation with Position and Heading Control using Complex Vector Fields for an Ackermann Steering Autonomous Vehicle

Three-Dimensional Off-Line Path Planning for Unmanned Aerial Vehicle Using Modified Particle Swarm Optimization

Robust Controller Design for an Autonomous Underwater Vehicle

Evolutionary Path Planning and Navigation of Autonomous Underwater Vehicles

LAIR. UNDERWATER ROBOTICS Field Explorations in Marine Biology, Oceanography, and Archeology

Prediction-Based Path Planning with Obstacle Avoidance in Dynamic Target Environment

Research Article Regressive Structures for Computation of DST-II and Its Inverse

Time Constrained Randomized Path Planning Using Spatial Networks

Elastic Bands: Connecting Path Planning and Control

Introduction to Information Science and Technology (IST) Part IV: Intelligent Machines and Robotics Planning

Announcements. Exam #2 next Thursday (March 13) Covers material from Feb. 11 through March 6

Engineers and scientists use instrumentation and measurement. Genetic Algorithms for Autonomous Robot Navigation

Safe Prediction-Based Local Path Planning using Obstacle Probability Sections

EE631 Cooperating Autonomous Mobile Robots

Research Article Path Planning Using a Hybrid Evolutionary Algorithm Based on Tree Structure Encoding

Adapting the Wavefront Expansion in Presence of Strong Currents

Navigation and Metric Path Planning

Research Article Implementation of Personal Health Device Communication Protocol Applying ISO/IEEE

Research Article A Novel Metaheuristic for Travelling Salesman Problem

Neural Networks for Obstacle Avoidance

Variable-resolution Velocity Roadmap Generation Considering Safety Constraints for Mobile Robots

EE631 Cooperating Autonomous Mobile Robots

CONSTRUCTION OF THE VORONOI DIAGRAM BY A TEAM OF COOPERATIVE ROBOTS

Mobile Robot Path Planning in Static Environments using Particle Swarm Optimization

UNMANNED UNDERWATER VEHICLE SIMULATOR ENABLING THE SIMULATION OF MULTI- ROBOT UNDERWATER MISSIONS WITH GAZEBO

Interval based 3D-mapping of an island with a sector-scan sonar

Fast Local Planner for Autonomous Helicopter

Ian Mitchell. Department of Computer Science The University of British Columbia

A New Algorithm for Automated Aircraft Conflict Resolution

CS4733 Class Notes. 1 2-D Robot Motion Planning Algorithm Using Grown Obstacles

S-SHAPED ONE TRAIL PARALLEL PARKING OF A CAR-LIKE MOBILE ROBOT

A Reactive Bearing Angle Only Obstacle Avoidance Technique for Unmanned Ground Vehicles

Go with the Flow: Optimal AUV Path Planning in Coastal Environments

REINFORCEMENT LEARNING: MDP APPLIED TO AUTONOMOUS NAVIGATION

Agent Based Intersection Traffic Simulation

Motion Planning. Howie CHoset

Research Article Average Bandwidth Allocation Model of WFQ

Transactions on Information and Communications Technologies vol 16, 1996 WIT Press, ISSN

Efficiency. Narrowbanding / Local Level Set Projections

A unified motion planning method for a multifunctional underwater robot

Research Article Cross Beam STAP for Nonstationary Clutter Suppression in Airborne Radar

Real-time Obstacle Avoidance and Mapping for AUVs Operating in Complex Environments

inter.noise 2000 The 29th International Congress and Exhibition on Noise Control Engineering August 2000, Nice, FRANCE

Robotics Project. Final Report. Computer Science University of Minnesota. December 17, 2007

DETERMINATION OF THE OPTIMUM PATH ON THE EARTH S SURFACE. (extended abstract) Abstract

Geometry Definitions and Theorems. Chapter 9. Definitions and Important Terms & Facts

Optimal Trajectory Generation for Nonholonomic Robots in Dynamic Environments

Research Article Well-Suited Similarity Functions for Data Aggregation in Cluster-Based Underwater Wireless Sensor Networks

Functional Discretization of Space Using Gaussian Processes for Road Intersection Crossing

Research Article Data Visualization Using Rational Trigonometric Spline

arxiv: v1 [cs.ro] 2 Sep 2017

Vision-based Localization of an Underwater Robot in a Structured Environment

Obstacle Avoidance Project: Final Report

Graphs, Search, Pathfinding (behavior involving where to go) Steering, Flocking, Formations (behavior involving how to go)

Lab 6 - Ocean Acoustic Environment

Cluster Subgraphs Example, With Tile Graphs. Alternatives. Cluster Subgraphs. Cluster Subgraphs Example, With Tile Graphs

Open Access The Kinematics Analysis and Configuration Optimize of Quadruped Robot. Jinrong Zhang *, Chenxi Wang and Jianhua Zhang

Research Article Modeling and Simulation Based on the Hybrid System of Leasing Equipment Optimal Allocation

Research Article Secure Data Aggregation in Wireless Multimedia Sensor Networks Based on Similarity Matching

Safe prediction-based local path planning using obstacle probability sections

A new method for determination of a wave-ray trace based on tsunami isochrones

Planning With Uncertainty for Autonomous UAV

Research Article Cooperative Signaling with Soft Information Combining

A New Performance-Based Motion Planner for Nonholonomic Mobile Robots

Other Voronoi/Delaunay Structures

Robot Path Planning Method Based on Improved Genetic Algorithm

Probabilistic Methods for Kinodynamic Path Planning

CMU-Q Lecture 4: Path Planning. Teacher: Gianni A. Di Caro

DVFH - VFH*: Reliable Obstacle Avoidance for Mobile Robot Navigation Coupled with A*Algorithm Through Fuzzy Logic and Knowledge Based Systems

Mathematics in Orbit

1 Mission Level Design. of Autonomous Underwater Vehicles

UAVs Task and Motion Planning in the Presence of Obstacles and Prioritized Targets

Autonomous Robotics 6905

Multi-Objective Path Planning using Spline Representation

Planning for Simultaneous Localization and Mapping using Topological Information

Final Report: Dynamic Dubins-curve based RRT Motion Planning for Differential Constrain Robot

Robot Motion Planning

by Flow Simulation Recontouring of Jet Engine Compressor Blades

Visual Navigation for Flying Robots. Motion Planning

Robot Motion Control and Planning

Application of Hybrid A* to an Autonomous Mobile Robot for Path Planning in Unstructured Outdoor Environments

Lab #5 Ocean Acoustic Environment

INTEGRATING LOCAL AND GLOBAL NAVIGATION IN UNMANNED GROUND VEHICLES

Autonomous navigation in industrial cluttered environments using embedded stereo-vision

Vision Guided AGV Using Distance Transform

Detection and Motion Planning for Roadside Parked Vehicles at Long Distance

Humanoid Robotics. Path Planning and Walking. Maren Bennewitz

Mapping Contoured Terrain Using SLAM with a Radio- Controlled Helicopter Platform. Project Proposal. Cognitive Robotics, Spring 2005

Sung-Eui Yoon ( 윤성의 )

Trajectory Optimization

ON THE DUALITY OF ROBOT AND SENSOR PATH PLANNING. Ashleigh Swingler and Silvia Ferrari Mechanical Engineering and Materials Science Duke University

Cooperative Relative Localization for Moving UAVs with Single Link Range Measurements

Transcription:

ISRN Robotics Volume 213, Article ID 329591, 6 pages http://dx.doi.org/1.542/213/329591 Research Article Unified Trajectory Planning Algorithms for Autonomous Underwater Vehicle Navigation Oren Gal Technion, Israel Institute of Technology, 32 Haifa, Israel Correspondence should be addressed to Oren Gal; orengal@tx.technion.ac.il Received 9 April 213; Accepted 19 May 213 Academic Editors: G. C. Gini, J.-S. Liu, and M. Oussalah Copyright 213 Oren Gal. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. This paper presents two efficient methods for obstacle avoidance and path planning for Autonomous Underwater Vehicle (AUV). These methods take into account the dynamic constraints of the vehicle using advanced simulator of AUV considering low level control and stability effects. We present modified visibility graph local avoidance method and a spiral algorithm for obstacle avoidance. The algorithms were tested in challenged scenarios demonstrating safe trajectory planning. 1. Introduction Path planning and obstacle avoidance are an important issues for Autonomous Underwater Vehicles (AUVs), and as can be noticed lately, these fields are extensively studied. Currents disturbances can have a big influence in a different water depth and should be considered in underwater environments [1]. For the first time, Soulignac et al. [2] faced with strong current situations. Online replanning in 3D underwater environments with strong, dynamic, and uncertain currents was presented in [3]. AUV also suffers from a limited energy source, which can be minimized by planning optimal trajectories, extending the working time of the vehicle. Optimal path planning using A algorithm search was presented by Carroll et al. [4]. Garau et al. [5]. Pêtrès et al. [6, 7] used similar concepts using BF (breadth first) search named FM and FM. As the presented methods in this paper, other methods do not use grid-based search. The path is presented with a series of points, which are connected one by one. Path planning problem transformed to a constrained optimization problem in terms of the coordinates of these points, generating optimal paths considering AUV s dynamic constraints [1, 8 1]. One of the most known limitations in motion planning algorithms related to real-time computation ability. Planning methods based on local perceptions are computationally less expensive and thus time efficient. Bui and Kim [11] and Kanakakis et al. [12] apply fuzzy logic approaches to AUV pathplanning.antonellietal.[13, 14] integratevirtualforce field (VFF); all of these are not optimal planning paths methods. This paper presents several different AUV path planning algorithms avoiding obstacle based on local perception abilities based on forward looking sonar. The introduced algorithms inherently take into account AUVs dynamic and kinematic constraints. AUV trajectory is simulated as described later. Simulations in typical underwater environments are presented, demonstrating algorithm s capabilities. 2. Vehicle Simulator AUV platforms are known as underactuated vehicles models; these kinds of models generate dynamic constraints. In such cases, obstacle avoidance algorithms must consider the dynamic envelop of the vehicles not entering unstable states.wedevelopedandimplementedauvsimulator.the simulator takes into account the dynamic model of the AUV and the behavior of the sensors. The simulator uses only Forward Looking Sonars (FLS) for obstacle detection. In addition to FLS we have Side Scan Sonars (SSS) for detail

2 ISRN Robotics 1 2 3 8 1 4 5 6 7 9 Figure 1: Block diagram of the AUV simulator [Braginsky, BGU]. scanning of the environment. Block diagram of simulator is presented in Figure 1. 3. Path Planning Algorithms 3.1. Visibility Graph Algorithm 3.1.1. General Description. The main idea behind our algorithm is the well-known visibility graph method. However, visibility graph algorithm is an abstract one that cannot be used on-line in real time applications, so many changes and adjustments are needed. Basically, visibility graph describes the concept of connecting nodes from start point to the goal; nodes are connected to each other if and only if there are no obstacles between them, and the next node is visible from the current node. Another major issue of the algorithm logic is concerned with node s creation. One of the assumptions of the abstract visibility graph algorithm is that the obstacles are polygonal. On the real world we cannot assume that every obstacle is polygonal. Although, we can trap the obstacle in a polygonal shapebutperformancesarelikelytobeveryconservative.we propose a new concept constructing nodes around obstacles, enabling safety motion of the vehicle. After constructing a full graph of nodes and arcs, search overthegraphcanbeusedforoptimaltrajectory,suchasa algorithm. Algorithm completeness is based on the fact that the marine environment is quasi dynamic and a free visible node can be found at each time step, taking into account kinematic and dynamic models. 3.1.2. Algorithm Description. Our algorithm can be divided into three steps: (1) building visibility graph; (2) connecting nodes; and (3) searching over the graph for the best solution. (1) Building the graph s nodes: first, it should be noted that the visibility graph is a global algorithm and thus is not an efficient computation time one. Therefore, we try to minimize node s number to achieve fast calculation times. Unseen node Figure 2: Obstacle s nodes. Seen node that will be left out due to foreseen collision with the obstacle Legal node Figure 3: Risky nodes being reduced due to maneuver constraints. We divide the nodes into two different cases. (a) Vehicle nodes: we add nodes related to the vehicle. For example, a node straight ahead of the vehicle (that represents straight path without turnings), a node to the right of the vehicle (representing a right turn), and a node to the left. The easiest way is to choose a radius that will allow turning to the node and take nodes on a circle around the start position with this chosen radius. That way we have nodes that represent all turning directions as well as slowing down if needed. We assume that the obstacles are round by spheres generating smooth trajectories due to dynamic profile andauvsturningradius.incaseofanobstaclewhich is not round, we bound the obstacle with a circle, considering safety.

ISRN Robotics 3 12 96 1 94 8 92 6 9 4 89 2 86 84 2 2 2 4 6 8 1 12 6 4 2 2 4 6 (a) (b) 26 24 22 2 18 16 14 2 2 4 6 8 1 12 14 16 (c) Figure 4: OA with spiral algorithm with typical scenarios (a), (b), and (c). 7 6 5 4 3 2 1 1 Obstacle vessel path vessel to target path Vessel Waypoint 1 1 2 3 4 5 6 7 8 9 1 2 Obstacle 15 1 5 vessel path 5 α 1 5 1 15 2 25 3 35 4 45 5 Figure 5: Vehicle moving to a waypoint changing vehicle course. Figure 6: Calculating alpha angle in a collision course. The next step is to choose the location of the nodes. For each obstacle we calculate the tangent line to the line that connects the obstacle with the robot and add two nodes on this line (the tangent line) right outside of the circle of the obstacle, as can be shown in Figure 2.

4 ISRN Robotics 2 15 1 Obstacle Nodes are denoted as blue squares. A square, with an X inside it, is a node that was passed by the vehicle. In (b) and (c), we can see a close up of two types of nodes. In (c) an obstacle node and in (b) a vehicle related node are presented. 5 5 1 vessel path 15 New waypoints 5 1 15 2 25 3 35 4 45 5 Figure 7: Added waypoints avoiding an obstaces. (b) Goal nodes: we add nodes that will help achieve the planned mission. For example, if the vehicle mission is to follow a line, adding nodes on that line will help the vehicle stay close to the original line, defining these nodes s cost as favor one when there are no obstacles. (2) Connecting nodes with arcs: weconnectthenodes usingthesimpleideaofvision.ifonenodecanbeseenfrom the other node, these nodes are connected, as can be seen in Figure 3. We add a simple measure considering the turning radius of the vehicle considering the maneuver constraints based on the turning radius of the vehicle. Otherwise, we do not connect this node. (3) Search an optimal path: the last step is to search the optimal path over the graph from the current position of thevehicletothegoal.asmentionedabove,weusedana algorithm to search the shortest path in the visibility graph. We define a cost function, taking into account and the shortest distance from the line connecting targets (the line connectingthefirstnodeandthelastnodeateachiteration). For the two nodes n and m,wesetthecostofthearcthat connects them as F (n, m) = Dist (n, m). (1) Then we define the price for the node n as: P (n) = Dist From Target Line (n, s, t) + Cost Of Path (n), (2) where s is the start node, t is the target node, and the Cost Of Path function is the sum of the costs of the arc connecting this nodetothepreviousnodeandthecostoftheprevious.based on that, the shortest path between start to goal is calculated. 3.1.3. Algorithm Simulations. In Figure 4, we present dense environment demonstrating visibility graph algorithm. Start point is located in the bottom left and goal point is located in the bottom right. It can be seen in (a) the nodes that were added during the time the vehicle was advancing its path. 3.2. Spiral Algorithm 3.2.1. General Description. The spiral algorithm presented for the first time in this paper, as its name may suggest, prefers to stay away from the vehicle s predetermined path and locally avoids an obstacle when there is an imminent danger of a collision. In such a case, circumnavigate of an obstacle is done by adding and placing new waypoints in the vehicle s path, in a semicircle around the center of the original, and in such a radius the vehicle can pass safely by the obstacle. For simplicity, we assume all obstacles are spheres. In the next section, a detailed algorithm stages are described. 3.2.2. Algorithm Stages. In this section, we introduce a more detailed description of the algorithm. We distinguish between two basic situations and the relevant action in a case of an obstacle that may cause collision. (A) No obstacles were detected: in general, if there are no obstacles, the algorithm output would be diving and moving at nominal velocity at the same direction of the next waypoint in the vehicle s path. The nominal velocity and course can be easily converted to a waypoint. (B) An obstacle was detected:whenanobstacleisdetected, several calculations must be done in order to decide whether the obstacle is endangering the vehicle: (i) First, we check if there is a collision between the course of the vehicle and the obstacle. CheckIntersect is a subalgorithm that performs a set of vector calculations to decide whether or not a sphere and a line segment intersect. The intersection points of the sphere and infinite line containing the line segment, if any, are given by (5) based on(3) and(4): discrim =((obs loc vessel loc ) vessel path ) 2 obs loc vessel loc 2 +(obs rad vessel loc )2 (3) sol 1/2 = (obs loc vessel loc ) vessel path ± sqrt (discrim) vessel loc 2 (4) intersect point1/2 = vessel loc + sol 1/2 (vessel path ) (5) CheckIntersect variable returns a true value if an intersection occurs within the finite line segment and false otherwise, where obstacle is the obstacle s center coordinates and radius and vessel path is the line segment starting in the vehicle s location, pointing to vehicle s heading. We limit vessel path s length to be the distance given by min (5, vessel loc next targetloc ). (6)

ISRN Robotics 5 y 1 9 8 7 6 5 y 1 9 8 7 6 5 4 4 3 3 2 2 1 1 2 4 x 6 8 1 2 4 6 x 8 1 12 (a) (b) 1 9 8 7 6 y 5 4 3 2 1 2 4 x 6 8 1 (c) Figure 8: OA with spiral algorithm with typical scenarios (a), (b), and (c). Obstacles detection ability is related to the vehicle perception abilities, and only local avoidance can be done. (ii) In the next stage, we calculate Dist (obstacle loc, vessel loc ) obs radius < safe marge, (7) where Dist(obstacle loc,vessel loc ) is the Euclidean distance between the obstacle and the vehicle, obs radius is the obstacle s radius, and safe marge is the minimum distance at which we would like to consider the obstacle as a dangerous one. (iii) Finally, we calculate CheckIntersect (obstacle, vessel to target path) function value, where vessel to target path is the line segment connecting the vehicle to the next waypoint as can be seen in Figure 5. The rational is not to consider obstacles which are not cause collision, in a case of vehicle s turns. (C) Avoiding an obstacle: in a case of a collision course between the obstacle and the vehicle, we add number of waypoints around the obstacle in a spiral trajectory and locally avoid the obstacle, in the following way. (i) First, we determine whether the center of the obstacle is to the left or the right of the vehicle s heading. We choose to pass the obstacle on the side that minimizes the distance between the updated and the original path. (ii) Next, a trigonometric calculation is done, to determine how wide is the bite taken out of the vehicle s path by the obstacle. This is done in terms of the angle alpha as can be seen in Figure 6. Alpha is given by α=2 sin 1 ( intersect point 1 intersect point2 2 obs rad ). (8)

6 ISRN Robotics If 135 < alpha 18, five waypoints are added on the boundary of the sphere centered located at the obstacle s center, with the augmented radius mentioned above. The waypoints are added so that the difference in degree between themisequal,ascanbeseeninfigure7. If 9 < alpha 135, 4 waypoints are added in the same fashion, and if alpha 9, 3 are added. 3.2.3. Algorithm Simulations. The algorithm has been tested using a simulation of a real vehicle simulator as described above. The scenarios chosen were 1 1km squares, in which were scattered 5 to 7 obstacles of radii 15 meters. The scenarios are modeled after underwater minefields. The algorithmhasshowngoodresultssofar,bothinobstacle avoidance and in adhering to a course. Simulations can be seen in Figure 8. 4. Conclusions In this research we presented two different methods for AUV, basedonlocalandglobalplanningmethods.oneofthe major challenges in AUV obstacle avoidance is related to the underactuated model of the vehicle. We challenge these constraints by simulating vehicle trajectory using advanced simulator, modeling the perception abilities of the forward looking sonar, and the obstacles detection in realtime. We presented a general description for each one of the algorithm and a simulation results with a typical search fields of AUVs. Further research includes sea test with REMUS 1 AUV, testing the presented algorithms. [6] C. Pêtrès,Y.Pailhas,Y.Petillot,andD.Lane, Underwaterpath planing using fast marching algorithms, in Proceedings of the Oceans Europe,vol.2,pp.814 819,June25. [7] C. Pêtrès,Y.Pailhas,P.Patrón,Y.Petillot,J.Evans,andD. M. Lane, Path planning for autonomous underwater vehicles, IEEE Transactions on Robotics,vol.23,no.2,pp.331 341,27. [8] C. W. Warren, Technique for autonomous underwater vehicle route planning, IEEE Oceanic Engineering,vol.15,no. 3, pp. 199 24, 199. [9] A. Alvarez and A. Caiti, A genetic algorithm for autonomous underwater variability, in Proceedings of International Federation Automatic Control (IFAC) Conference Control Applications Marine Systems,Glasgow,UK,June21. [1] A. Alvarez, A. Caiti, and R. Onken, Evolutionary path planning for autonomous underwater vehicles in a variable ocean, IEEE Oceanic Engineering,vol.29,no.2,pp.418 429,24. [11] L.-D. Bui and Y.-G. Kim, An obstacle-avoidance technique for autonomous underwater vehicles based on BK-products of fuzzy relation, Fuzzy Sets and Systems,vol.157,no.4,pp.56 577, 26. [12] V. Kanakakis, K. P. Valavanis, and N. C. Tsourveloudis, Fuzzylogic based navigation of underwater vehicles, Intelligent and Robotic Systems, vol. 4, no. 1, pp. 45 88, 24. [13] G. Antonelli, S. Chiaverini, R. Finotello, and E. Morgavi, Realtime path planning and obstacle avoidance for an autonomous underwater vehicle, in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA 99), pp.78 83, Detroit, Mich, USA, May 1999. [14] G. Antonelli, S. Chiaverini, R. Finotello, and R. Schiavon, Real-time path planning and obstacle avoidance for RAIS: an autonomous underwater vehicle, IEEE Oceanic Engineering,vol.26,no.2,pp.216 227,21. References [1] D. Kruger, R. Stolkin, A. Blum, and J. Briganti, Optimal AUV path planning for extended missions in complex, fastflowing estuarine environments, in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA 7),pp.4265 427,April27. [2] M. Soulignac, P. Taillibert, and M. Rueher, Adapting the wavefront expansion in presence of strong currents, in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA 8), pp. 1352 1358, May 28. [3] D. R. Thompson, S. Chien, Y. Chao et al., Spatiotemporal path planning in strong, dynamic, uncertain currents, in Proceedings of the IEEE International Conference on Robotics and Automation (ICRA 1), pp. 4778 4783, May 21. [4]K.P.Carroll,S.R.McClaran,E.L.Nelson,D.M.Barnett, D. K. Friesen, and G. N. Williams, AUV path planning: an A approach to path planning with consideration of variable vehicle speeds and multiple, overlapping time-dependent exclusion zones, in Proceedings of the Symposium on Autonomous Underwater Vehicle Technology (AUV 92),pp.3 8,1992. [5] B. Garau, A. Alvarez, and G. Oliver, Path planning of autonomous underwater vehicles in current fields with complex spatial variability: an A approach, in Proceedings of the IEEE International Conference on Robotics and Automation,pp.194 198, Barcelona, Spain, April 25.

Rotating Machinery Engineering The Scientific World Journal Distributed Sensor Networks Sensors Control Science and Engineering Advances in Civil Engineering Submit your manuscripts at Electrical and Computer Engineering Robotics VLSI Design Advances in OptoElectronics Navigation and Observation Chemical Engineering Active and Passive Electronic Components Antennas and Propagation Aerospace Engineering Modelling & Simulation in Engineering Shock and Vibration Advances in Acoustics and Vibration