The effect of eight different cost functions on trunk muscle forces, spinal loads and stability was investigated. Kinematics-based approach combined with nonlinear finite element modeling and optimization were used to model in vivo measurements on isometric forward flexions at similar to 40 degrees and -65 degrees in sagittal plane with or without a load of 180N in hands. Four nonlinear (Sigma stress', Sigma stress(2), Sigma force' and muscle fatigue) and four linear (Sigma stress, Sigma force, axial compression and double-linear) criteria were considered. Predicted muscle activities were compared with measured EMG data. All predictions, irrespective of the cost function used, satisfied required kinetic, kinematics and stability conditions all along the spine. Four criteria (Sigma stress', Sigma stress(2), fatigue and double-linear) predicted muscle activities that qualitatively matched measured EMG data. The fatigue and double-linear criteria were inadequate in predicting greater forces in larger muscles with no consideration for their moment arms. Nearly the same stability margin was computed under these four cost functions. At the lower lumbar levels, the compression forces differed by < 20% and the shear forces by < 14% as various cost functions were considered. Smaller axial compression and anterior shear forces (by Z6%) were computed when only the active components rather than the total muscle forces were taken as unknown in the Tstress' cost function. Overall, one single cost function of Estress(2) or Estress' rather than a multi-criteria one was found sufficient and adequate in yielding plausible results comparable with measured EMG activities and disc pressure. (c) 2005 IPEM. Published by Elsevier Ltd. All rights reserved.