 Solution

Use the same real elastic stiffness tensor for both fields, including its minor and major symmetries, and let $t_i^a={\hat{\sigma }}_{ij}^a{n}_j$, with outward unit boundary normal $n_j$. No complex conjugation is needed for the initial bilinear identity. DefineThe divergence has two force-balance terms and two cross-gradient terms. The latter cancel becausethis follows by exchanging the two strain-gradient factors using major symmetry and the minor pair symmetries. The frequency-domain equations then giveIntegrating and applying the divergence theorem provesThis is the Betti identity for elastodynamic fields. Spatially varying stiffness is allowed because its derivatives remain inside the divergence. The original PDF's constitutive contraction is $c_{ijpq}{\partial }_p{\hat{u}}_q$; the TeX's repeated displacement field index is erroneous.