This paper presents a curvature method for analysis of beam-columns with different materials and arbitrary cross-section shapes and subjected to combined biaxial moments and axial load. Both material and geometric nonlinearities (the p-delta effect in this case) were incorporated. The proposed method considers biaxial curvatures and uniform normal strains of discrete cross-sections of beam-columns as basic unknowns, and seeks for a solution of the column deflection curve that satisfies force equilibrium conditions. A piecewise representation of the beam-column deflection curve is constructed based on the curvatures and angles of rotation of the segmented cross-sections. The resulting bending moments were evaluated based on the deformed column shape and the axial load. The moment curvature relationship and the beam-column deflection calculation are presented in matrix form and the Newton-Raphson method is employed to ensure fast and stable convergence. Comparison with results of analytic solutions and eccentric compression tests of wood beam-columns implies that this method is reliable and effective for beam-columns subjected to eccentric compression load, lateral bracings and complex boundary conditions.