Aerodynamic solver using the transonic small-disturbance (TSD) equation has frequently been used to perform practical aeroelastic analysis for many aircraft models. In the present study, the more accurate aeroelastic analysis solver using the TSD equation was developed by considering the viscous effects of the boundary-layer. The viscous effects were considered using Green's lag-entrainment equations and an inverse boundary-layer method. Through aerodynamic analyses for several aircraft wings, the viscous-inviscid interaction approach could improve the accuracy of the aerodynamic computation using the TSD equation. Finally, the aeroelastic characteristics were investigated using comparisons of the time responses between the inviscid and viscous flows.