We introduce a new class of harmonic univalent functions by using a generalized differential operator and investigate some of its geometric properties, like, coefficient estimates, extreme points and inclusion relations. Finally, we show that this class is invariant under Bernandi--Libera--Livingston integral for harmonic functions.