In this paper, a bivariate Lindley distribution using Morgenstern approach is proposed which can be used for modeling bivariate life time data. Some characteristics of the distribution like moment generating function, joint moments, Pearson correlation coefficient, survival function, hazard rate function, mean residual life function, vitality function and stress-strength parameter R = P r(Y < X), are derived. The conditions under which the proposed distribution is an increasing (decreasing) failure rate distribution and positive (negative) quadrant dependent is discussed. Also, the method of estimating model parameters and stress-strength parameter by maximum likelihood is elucidated. Numerical illustration using simulated data is carried out to access the estimates in terms of mean squared error and relative absolute bias.Keywords Farlie-Gumbel-Morgenstern family, maximum likelihood estimation, mean residual life, mean squared error, positive quadrant dependence, relative absolute bias, stress-strength parameter, vector hazard rate, vitality function AMS 2010 subject classifications 60E05