In this paper semi-smooth Newton methods for optimal control problems governed by the dynamical Lamé system are considered and their convergence behavior with respect to superlinear convergence is analyzed. Techniques from Kröner, Kunisch, Vexler (2011), where semi-smooth Newton methods for optimal control of the classical wave equation are considered, are transferred to control of the dynamical Lamé system. Three different types of control actions are examined: distributed control, Neumann boundary control and Dirichlet boundary control. The problems are discretized by finite elements and numerical examples are presented. 2 U , subject to y = S(u), y ∈ Y, u ∈ U ad ⊂ U ω with control space U ω , state space Y and α > 0. The control-to-state operator S : U ω → Y is assumed to be affine-linear, the functional G : Y → R to be quadratic. The control and state space and the operators are defined in more detail in the next section. The choice of the control-to-state operator incorporates distributed as well as Neumann and Dirichlet boundary control problems of the dynamical Lamé system which we will consider later. The set of admissible controls is defined by U ad = { u ∈ U ω | u a ≤ u ≤ u b } for given u a , u b ∈ U ω. To specify the control-to-state operator we introduce the dynamical Lamé system. Let Ω ⊂ R d , d = 1, 2, 3, be a bounded domain with C 2-boundary (bounded