Resin-impregnated paper (RIP) bushing has gained significant interest due to its extended application in Extra High Voltage (EHV) and Ultra High Voltage (UHV) electricity transmission systems. However, the design criterion of its overall structure, the geometry parameters of the condenser layers, and stress release devices, etc., are still not fully understood. This article proposes a unique electric field optimization technique to integrate both the analytical and the numerical methods. The charge simulation method (CSM) is employed to create the overall equipotential surface, within which the finite element analysis (FEA) is adapted to study the localized field enhancement effects, taking into consideration the multi-physics coupled fields. A case study is performed on an actual UHV bushing. The results are compared to the traditional methods, to demonstrate the benefit of the hybrid method.