The paper describes and compares in detail, directly accumulation method and global stiffness matrix method. Taking efficiency, stability and convergence into consideration, the paper puts forward an explicit-implicit combined method which is suitable for force models.