Перейдём к вычислению матрицы демпфирования для трёхмерного случая. Матрица демпфирования связана с интегралом квадрата пробной функции. Рассмотрим интеграл для одного тетраэдра с вершинами

Пробная функция на тетраэдре имеет вид $\upsilon_{(i)(i+3)}(x, y, z) = q_i \cdot \phi_i + q_{i+1} \cdot \phi_{i+1} + q_{i+2} \cdot \phi_{i+2} + q_{i+3} \cdot \phi_{i+3}$. Аналогично двумерному случаю, функции <<крышек>> для тетраэдра имеют линейный вид согласно (6.14$$\begin{cases} \phi_i(x, y, z) = a_i + b_i \cdot x + c_i \cdot y + d_i \cdot z\\ \phi_{i+1}(x, y, z) = a_{i+1} + b_{i+1} \cdot x + c_{i+1} \cdot y + d_{i+1} \cdot z\\ \phi_{i+2}(x, y, z) = a_{i+2} + b_{i+2} \cdot x + c_{i+2} \cdot y + d_{i+2} \cdot z\\ \phi_{i+3}(x, y, z) = a_{i+3} + b_{i+3} \cdot x + c_{i+3} \cdot y + d_{i+3} \cdot z \end{cases}$$)

Подставим пробную функцию в (6.31$$\int_{\text{тет}} \upsilon^2 \,dV.$$)

Для линейных функций <<крышек>> на тетраэдре справедливы следующие соотношения

где $V_{\text{тет}}$ — объём тетраэдра, который вычисляется по формуле (6.16$$V_{\text{тет}} = \frac{\displaystyle |\Delta|}{\displaystyle 6},$$).

Раскроем квадрат пробной функции

Применим формулы (6.33$$\begin{split} &\int_{\text{тет}} \phi_m \cdot \phi_n \,dV = \begin{cases} \frac{\displaystyle V_{\text{тет}}}{\displaystyle 10}, & m = n\\ \frac{\displaystyle V_{\text{тет}}}{\displaystyle 20}, & m \neq n \end{cases} \end{split}$$)

Упростим выражение

Введём обозначения для элементов локальной матрицы демпфирования тетраэдра

Таким образом, локальная матрица демпфирования для тетраэдрального элемента имеет вид

Локальная матрица демпфирования является симметричной. Глобальная матрица демпфирования $\mathbf{C}$ получается путём суммирования вкладов от всех тетраэдральных элементов сетки методом сборки: элементы локальных матриц добавляются к соответствующим элементам глобальной матрицы согласно глобальной нумерации узлов. Размерность глобальной матрицы демпфирования равна $N \times N$, где $N$ — общее количество узлов сетки.