A second order accurate method in the infinity norm is proposed for general three dimensional anisotropic elliptic interface problems in which the solution and its derivatives, the coefficients, and source terms all can have finite jumps across one or several arbitrary smooth interfaces. The method is based on the 2D finite element-finite difference (FEFD) method but with substantial differences in method derivation, implementation, and convergence analysis. One of challenges is to derive 3D interface relations since there is no invariance anymore under coordinate system transforms for the partial differential equations and the jump conditions. A finite element discretization whose coefficient matrix is a symmetric semi-positive definite is used away from the interface; and the maximum preserving finite difference discretization whose coefficient matrix part is an M-matrix is constructed at irregular elements where the interface cuts through. We aim to get a sharp interface method that can have second order accuracy in the point-wise norm. We show the convergence analysis by splitting errors into several parts. Nontrivial numerical examples are presented to confirm the convergence analysis.