This paper proposes a voltage stability constrained optimal power flow (VSC-OPF) for an unbalanced distribution system with distributed generators (DGs) based on semidefinite programming (SDP). The AC optimal power flow (ACOPF) for unbalanced distribution systems is formulated as a chordal relaxation-based SDP model. The minimal singular value (MSV) of the power flow Jacobian matrix is adopted to indicate the voltage stability margin. The Jacobian matrix can be explicitly expressed by ACOPF state variables. The nonlinear constraint on the Jacobian MSV is then replaced with its maximal convex subset using linear matrix inequality (LMI), which can be incorporated in the SDP-based ACOPF formulation. A penalty technique is leveraged to improve the exactness of the SDP relaxation. Case studies performed on several IEEE test systems validate the effectiveness of the proposed method.