The cone penetration tests have been employed extensively in both onshore and offshore site investigations to obtain the strength properties of soils. Interpretation of effective internal friction angle gyp' becomes complicated for cones in silty clays or clayey silts, since the soil around the advancing cone may be under partially drained conditions. Although there exist several robust methods to estimate gyp ' , the pore pressure at the cone shoulder has to be measured to represent the drainage conditions. Many cone penetrometers in practice are not equipped with a pore pressure transducer. Even for a piezocone, the pore pressure recorded in-situ may be unreliable due to the poorly saturated or clogged filter. These limitations prohibit the application of existing methods. Large deformation finite element analyses were carried out within the formula of effective stress to reproduce the cone penetrations under various drainage conditions. The numerical approach was validated against the existing model tests in centrifuge and chamber, with wide ranges of penetration rates and soil types. A backbone curve is proposed to estimate the normalized cone resistance varying with the normalized penetration rate. Based on the backbone curve, a procedure is developed to predict gyp' of cohesive soils under undrained or partially drained conditions, replacing the pore pressure with the normalized penetration rate. The procedure can be used for soils with an overconsolidation ratio no larger than 5.