En este documento, proponemos un método para derivar esquemas numéricos que conservan propiedades energéticas para sistemas acoplados de dos sistemas naturales dados. Consideramos el caso en el que los dos sistemas están interconectados por la ley de acción-reacción. Aunque los esquemas derivados se basan en el método del gradiente discreto, en el caso bajo consideración, la ecuación de movimiento no tiene la forma usual representada mediante el uso de la matriz antisimétrica. Por lo tanto, los esquemas que conservan propiedades energéticas no pueden obtenerse utilizando directamente el método del gradiente discreto. Mostramos resultados numéricos para dos sistemas acoplados como ejemplos; el primer sistema es una combinación de la ecuación de onda y la ecuación elástica, y el segundo es del sistema masa-resorte y la ecuación elástica.