Determination of collapse load-carrying capacity of elasto-plastic material is very important in designing structure. The problem is commonly solved by elasto-plastic finite element method （FEM）. In order to deal with material nonlinear problem involving strain softening problem effectively, a new numerical method-damped Newton method was proposed. The iterative schemes are discussed in detail for pure equilibrium models. In the equilibrium model, the plasticity criterion and the compatibility of the strains are verified, and the strain increment and plastic factor are treated as independent unknowns. To avoid the stiffness matrix being singularity or condition of matrix being ill, a damping factor a was introduced to adjust the value of plastic consistent parameter automatically during the iterations. According to the algorithm, the nonlinear finite element program was complied and its numerical example was calculated. The numerical results indicate that this method converges very fast for both small load steps and large load steps. Compared with those results obtained by analysis and experiment, the predicted ultimate bearing capacity from the proposed method is identical.
A new node-pairs contact algorithm is proposed to deal with a composite material or bi-material interface crack face contact and friction problem （e.g., resistant coating and thermal barrier coatings） subjected to complicated load conditions. To decrease the calculation scale and calculation errors, the local Lagrange multipliers are solved only on a pair of contact nodes using the Jacobi iteration method, and the constraint modification of the tangential multipliers are required. After the calculation of the present node-pairs Lagrange multiplier, it is turned to next contact node-pairs until all node-pairs have finished. Compared with an ordinary contact algorithm, the new local node-pairs contact algorithm is allowed a more precise element on the contact face without the stiffness matrix singularity. The stress intensity factors （SIFs） and the contact region of an infinite plate central crack are calculated and show good agreement with those in the literature. The contact zone near the crack tip as well as its influence on singularity of stress fields are studied. Furthermore, the frictional contacts are also considered and found to have a significant influence on the SIFs. The normalized mode-II stress intensity factors KII for the friction coefficient decrease by 16% when f changes from 1 to 0.
