Thanks QShake for the hint. Currently I do not have an access to AASHTO reference. However my goal is to find out how the analysis is done, manually, and without using any charts or computer programs. Let us look at single moment case, eccentricty e=M/p, we compare this with L/6 value, where L is footing lenght in moment direction. if e<L/6, flexural equation is used (i.e. q=P/A+-M/S). on the other hand if e>L/6 flexural equation shall not be used and the static eqautions must be used (segma V=0, and segma M=0). Two equations for two unknowns, qmax and Leff. Now in the case of the biaxial bending, there are more unkowns than the number of equations you can generate. Any hint. thanks.