In this study, we introduce a two-dimensional (2D) human skeletal model coupled with knee, spine, and shoulder exoskeletons. The primary purpose of this model is to predict the optimal lifting motion and provide torque support from the exoskeleton through the utilization of inverse dynamics optimization. The kinematics and dynamics of the human model are expressed using the Denavit-Hartenberg (DH) representation. The lifting optimization formulation integrates the electromechanical dynamics of the DC motors in the exoskeletons of the knee, spine, and shoulder. The design variables for this study include human joint angle profiles and exoskeleton motor current profiles. The optimization objective is to minimize the squared normalized human joint torques, subject to physical and task-specific lifting constraints. We solve this optimization problem using the gradient-based optimizer SNOPT. Our results include a comparison of predicted human joint angle profiles, joint torque profiles, and ground reaction force (GRF) profiles between lifting tasks with and without exoskeleton assistance. We also explore various combinations of exoskeletons for the knee, spine, and shoulder. By resolving the lifting optimization problems, we designed the optimal torques for the exoskeletons located at the knee, spine, and shoulder. It was found that the support from the exoskeletons substantially lowers the torque levels in human joints. Additionally, we conducted experiments only on the knee exoskeleton. Experimental data indicated that using the knee exoskeleton decreases the muscle activation peaks by 35.00%, 10.03%, 22.12%, 30.14%, 16.77%, and 25.71% for muscles of the erector spinae, latissimus dorsi, vastus medialis, vastus lateralis, rectus femoris, and biceps femoris, respectively.