In this paper a method for simultaneous swift non-linear analysis and optimal design/posture of mechanical/biomechanical systems is presented. The method is developed to get advantages of iterations in non-linear analysis and/or generations in genetic algorithm (GA) for the purpose of efficient analysis within the optimal design/posture. The method is applicable for both size and geometry optimizations wherein material and geometry non-linearity are present. In addition to established mechanical systems, the method can solve biomechanical models of human musculoskeletal system. Optimization-based procedures are popular methods for resolving the redundancy at joints wherein the number of unknown muscle forces is far more than the number of equilibrium equations. These procedures involve optimization of a cost function(s) which is assumed to be consistent with the central nervous system's strategy when activating muscles to assure equilibrium. However, because of the complexity of biomechanical problems (i.e., due to non-linear biomaterial, large deformation, redundancy of the problem and so on) efficient analysis are required within optimization procedures as suggested in this paper.