Research Papers

A Theory in Finite Differences for Computation of Interchanging Principal Stresses and Principal Axes Rotation

[+] Author and Article Information
S. H. Stefanov

Department of Mechanical Engineering and Automation,
University of Forestry (LTU),
Faculty of Forest Industry,
10 Kliment Ohridski Blvd.,
Sofia 1756,Bulgaria
e-mail: stefanov_sh@yahoo.com;

Manuscript received February 27, 2012; final manuscript received September 19, 2012; accepted manuscript posted October 10, 2012; published online May 16, 2013. Assoc. Editor: Daining Fang.

J. Appl. Mech 80(4), 041009 (May 16, 2013) (8 pages) Paper No: JAM-12-1083; doi: 10.1115/1.4007796 History: Received February 27, 2012; Revised September 19, 2012; Accepted October 10, 2012

Suppose consecutive ordinates of three arbitrary and nonproportional stress-time functions of plane state of stress are entered into a computer by a little finite time step. The theory proposed solves the following problem: correct computation of the ordinates of the principal stress-time functions and the angle of principal axes rotation. This problem is not as simple as researchers approached it prior to the computer era. First of all, the correct solution for the principal stresses and the principal axes rotation require correct interchange of the principal stresses while computing them, i.e., correct interchange of the plus/minus signs in the well-known equations for them. For the interchange analysis, an ellipse of stress transformation in the three-dimensional stress-coordinate space is revealed. By changing a coordinate scale, the ellipse turns into a circumference that is an analog to, but different from, a Mohr circle. The correct solution also requires treating the principal axes rotation in little finite differences per little time differences during which little finite elements appear as building the stressing path in the three-dimensional stress-coordinate space. Based on the ellipse/circumference mentioned, three interchange conditions are revealed. The third one is the most important. And, a necessity is also revealed for dividing some stressing path's elements into two subelements. Based on all the findings, the main commands of an algorithm for computing the ordinates of the principal stress-time functions and the angle of principal axes rotation are presented. The correct solution of the problem has been achieved thanks to new notions taken from the so-called integration of damage differentials (IDD) theory. In fact, the paper presents a new contribution to the variable plane stress state analysis.

Copyright © 2013 by ASME
Your Session has timed out. Please sign back in to continue.


Beer, F., Johnston, E., and DeWolf, J., 2002, Mechanics of Materials, 3rd. ed., McGraw-Hill, New York.
Stefanov, S., 2007, “Strength of Materials,” Leso-Tehnicheski Universitet, Sofia, Bulgaria, (in Bulgarian).
Findley, W., Mathur, P., Szczepanski, E., and Temel, A., 1961, “Energy Versus Stress Theories for Combined Stress—A Fatigue Experiment Using a Rotating Disk,” ASME J. Basic Eng., 83, pp. 10–14. [CrossRef]
Tong, Z., Zhang, J., Yu, Y., and Zhang, G., 2010, “Drained Deformation Behavior of Anisotropic Sands During Cyclic Rotation of Principal Stress Axes,” J. Geotech. Geoenv. Eng., 136, pp. 1509–1518. [CrossRef]
Carpinteri, A., Spagnoli, A., Vantadori, S., and Viappiani, D., 2008, “A Multiaxial Criterion for Notch High-Cycle Fatigue Using a Critical-Point Method,” J. Eng. Fract. Mech., 75, pp. 1864–1874. [CrossRef]
Carpinteri, A., Spagnoli, A., and Vantadori, S., 2009, “Multiaxial Fatigue Life Estimation in Welded Joints Using the Critical Plane Approach,” Int. J. Fatigue, 31, pp. 188–196. [CrossRef]
Miller, W., Ohji, K., and Marin, J., 1966, “Rotating Principal Stress Axes in High Cycle Fatigue,” ASME Paper No. 66-WA/Met 9.
Stefanov, S., 2011, “Integration of Damage Differentials (IDD): General Presentation and Verification Results,” Int. J. Fatigue, 33, pp. 977–984. [CrossRef]
Stefanov, S., 2012, “IDD (Integration of Damage Differentials),” http://www.freewebs.com/fatigue-life-integral
Brighenti, R., Carpinteri, A., and Vantadori, S., 2012, “Fatigue Life Assessment Under a Complex Multiaxial Load History: An Approach Based on Damage Mechanics,” Fatigue Fract. Eng. M., 35, pp. 141–153. [CrossRef]
Brighenti, R., and Carpinteri, A., 2012, “A Notch Multiaxial-Fatigue Approach Based on Damage Mechanics” Int. J. Fatigue.39, pp. 122–133. [CrossRef]
Stefanov, S., 1997, “The Curvilinear Integral Method: A Question to the Critical Plane Concept,” Int. J. Fatigue, 19, pp. 101–107. [CrossRef]


Grahic Jump Location
Fig. 1

(a) An initial immovable cuboid of plane stress state with variable stresses σx = σx(t), σy = σy(t), and τxy = τxy(t); (b) another immovable cuboid at an angle α; (c) rotating principal cuboid at the angle α′ = α′(t)

Grahic Jump Location
Fig. 2

(a) Three arbitrary and nonproportional oscillograms σx(t), σy(t) and τxy(t); (b) σxyxy loading path called variant trajectory (Sv) and its finite element Δsv per Δt

Grahic Jump Location
Fig. 3

Transforming ellipse

Grahic Jump Location
Fig. 4

A view of the transforming ellipse from the arrowhead of the τxy axis

Grahic Jump Location
Fig. 5

Transforming circumference and principal axes

Grahic Jump Location
Fig. 6

Transformation of variant Δsv ≡ MpM finite element from σxyxy coordinate space into invariant Δs ≡ Mp,0M′ element in σ′-σ″-Δτ coordinate space, and other additional illustrations

Grahic Jump Location
Fig. 7

(a) Preceding (sp,0) stress state at tp = t − Δt on immovable cuboid that will be principal at t; (b) invariant (Δs) little finite stress state difference during Δt on the same cuboid; (c) invariant (s) stress state at t when the same cuboid is principal

Grahic Jump Location
Fig. 8

Circular transformation of (projection of) MpM variant element into invariant Mp,0M′, and other additional illustrations

Grahic Jump Location
Fig. 10

Variant (Sv) trajectory intersecting ηv axis (O point)

Grahic Jump Location
Fig. 9

Transformation of MpM variant element projection into invariant Mp,0M′, and other additional illustrations in Oξv*τxy coordinate plane



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In