Incorporating moment-rotation springs to 1D FEA
Incorporating moment-rotation springs to 1D FEA
(OP)
Does anyone know how to include moment-rotation springs into a 1D FEA model, or know of a reference which explains this? I've got my own program to model a laterally loaded pile with p-y springs, but now I want to incorporate M-θ springs and was wondering how to implement them. Specifically, the effect of the M-θ springs on the load vectors: can I include the M-θ spring effect as an additional set of loads to the element load vector? Do all terms in the element load vector (F1, M1, F2, M2) get affected, or just the moment terms?
Thanks in advance for any help.
Thanks in advance for any help.





RE: Incorporating moment-rotation springs to 1D FEA
RE: Incorporating moment-rotation springs to 1D FEA
So:
(kp-y + kM-θ).[mM]
Rather than the current:
kp-y.[mM]
???
Would the units be correct? Moreover, the units of the M-θ gradient are different to the p-y stiffness.
RE: Incorporating moment-rotation springs to 1D FEA
The element stiffness matrix without spring is:
k11 k12 k13 k14 k15 k16
k21 k22 k23 k24 k25 k26
k31 k32 k33 k34 k35 k36
k41 k42 k43 k44 k45 k45
k51 k52 k53 k54 k55 k56
k61 k62 k63 k64 k65 k66
Add spring stiffness the matrix will be like
k11* k12 k13 k14 k15 k16
k21 k22* k23 k24 k25 k26
k31 k32 k33* k34 k35 k36
k41 k42 k43 k44* k45 k45
k51 k52 k53 k54 k55* k56
k61 k62 k63 k64 k65 k66*
k11*=k11+spring k22*=k22+spring
RE: Incorporating moment-rotation springs to 1D FEA
Nonetheless, from memory, the governing equation is:
EI.d4w/dx4 + P.d2w/dx2 + kp-y.dw/dx = q
Where:
EI = Flexural rigidity of pile
kp-y = Stiffness of p-y curve
P = Axial load
When solving with FEA and using matrix notation:
([kM] + [mM] + [gM]).{y} = {Q}
Where:
[kM] = Element stiffness matrix (for 'internal' pile forces caused by EI)
[mM] = Element mass matrix (for 'external' pile forces caused by soil reaction - currently just p-y springs incorporated; p values are distributed loads across the element length)
[gM] = Element geometric matrix (for inclusion of axial force across elements)
{y} = Displacement vector
{Q} = Load vector
Where do moment-rotation springs fit into the above governing equation? Moreover, what is the governing equation when including moment-rotation springs?
The units of the element stiffness matrix appear incompatible with the stiffness of a moment-rotation curve (units of moment-rotation stiffness would be kNm per radian for a lumped spring or kNm/m per radian for a distributed moment reaction)
RE: Incorporating moment-rotation springs to 1D FEA
Your moment rotation springs become an external force dependent on displacement, follow me for a single line of the matrix system
k_n*d_n = f_n
k = stiffness
d = displacement
f = force
n is the line in the matrix being represented.
f_n however is a function of d_n
f_n = m_n*d_n
m_n = rotational stiffness
Now substitute
k_n*d_n = m_n*d_n
There is now a problem. You can't solve this system since you need all d_n terms on the left side. So we rearrange
k_n*d_n - m_n*d_n = 0
combine terms to make your stiffness term more apparent
(k_n - m_n)*d_n = 0
So, to include your rotational stiffness you need to subtract off the stiffness of your M-θ springs.
RE: Incorporating moment-rotation springs to 1D FEA
Do you have a reference explaining the basis for what you discuss? I think you are both coming at this from a structural engineering perspective where there is only one element matrix (the element stiffness matrix). To model distributed loads along a pile you need other matrices, as I mentioned in my previous post.
RE: Incorporating moment-rotation springs to 1D FEA
Also, I'm not coming at this from a single element perspective,the concept is the same for all fea, whether performing a single element or many elements. The global stiffness matrix for a many element system is built entirely from local stiffness matrices of single elements.
RE: Incorporating moment-rotation springs to 1D FEA
RE: Incorporating moment-rotation springs to 1D FEA
I didn't say you were coming at this from a single element perspective. I presumed that you were more familiar with structural analysis where there is typically just the element stiffness matrix rather than the other two matrices I mentioned in my 3rd post in this thread.
The element mass matrix is not just used for dynamic analyses. A beam on an elastic foundation also has an element mass matrix too. For a 4DOF laterally loaded beam-column (pile):
[mM] =
CODE -->
multiplied by kp-y.L/420
The above is pretty standard for what I'm looking at and is referenced in several places, including 'Programming the Finite Element Method' by Smith, Griffiths and Margetts.
I suspect I could include the M-θ reaction in the element load vector, though I need to know how to resolve a distributed moment. I do wonder whether I need to add another matrix similar to the above, however.
RE: Incorporating moment-rotation springs to 1D FEA
Reference:
DAVIDSON, H. L., CASS, P. G., KHILJI, K. H. and MCQUADE, P. V., 1982. Laterally Loaded Drilled Pier Research, Volume 1: Design Methodology. Monroeville: GAI Consultants Inc.
Contacting the Electric Power Research Institute (whose address was 3412 Hillview Avenue, Palo Alto, California 94304 in 1982) will result in an auto-reply with a link to their portfolio. The above reference can be accessed from there.
RE: Incorporating moment-rotation springs to 1D FEA
1. If you are combining it with your stiffness matrix, you have inconsistent units
2. If you aren't combining with your stiffness matrix, you now must increase the size of all matrices involved to make the addition of the compatible with each other
Maybe we're perhaps speaking two different languages, but trying to say the same thing. I've always seen what you are calling the mass matrix called the foundational stiffness matrix, not mass. For the element equation you have shown above Mass*displacement =/= force, so you have inconsistent units there if this isn't a stiffness matrix of some sort.
I'm also not sure about your governing equation. It seems to me, if I'm thinking correctly on you p-y springs and they are foundational springs, then their stiffness should be multiplied by w instead of d2w/dx2 since it is related to the amount the pile displaces. Additionally, I would have to spend more time with it, but the axial term looks correct for a bar that is loaded only axially, but looks incorrect for a beam in bending.
RE: Incorporating moment-rotation springs to 1D FEA
I think this is almost certainly a case of different engineering fields having the same names for slightly different things. Indeed that was what I was alluding to before. Indeed mass doesn't actually come into the [mM] matrix at all for this sort of analysis, so I think the name was borrowed.
You are also quite right that the equation I gave in my third post in this thread should've actually been:
EI.d4w/dx4 + P.d2w/dx2 + kp-y.dw/dx = q
I have corrected the post above to reflect this. Thanks for pointing it out - I had implemented it as above, though I didn't write it down here correctly.
The axial term is for axial loading of an element. Is there a different term for axial + bending? Is this where I need to upgrade it to a Timoshenko beam (i.e. to account for shear deformations)? That is on my to-do list. Or is there some other matrix that is used instead?
RE: Incorporating moment-rotation springs to 1D FEA
EI has units of kN-m2 not kN/m2.
E = kN/m2
I = m4
EI = kN-m2
I also don't understand your division by y on the mass matrix though it may be my lack of understanding of p-y springs. It may also be not understanding the terms you have contained in your y vector. Going back to your governing equation w(x) has units of m dw/dx is unitless since it's an angle. Units look something like this for derivatives of w(x):
w(x) = m
dw/dx = unitless
d2w/dx2 = 1/m
d3w/dx3 = 1/m2
d4w/dx4 = 1/m3
From what I see in your governing equation you have the first term having units of kN/m, the second term, if P = AE/L, has units of kN/m2 and the third term has units of kN/m2 based on what you said about the k-y springs.
As for the timoshenko beam, it's been a while since I looked at it, and it's not as easy to remember as classical beam theory, so I would have to revisit that, which I can't do at the moment.
RE: Incorporating moment-rotation springs to 1D FEA
The governing equation and matrices are published in a number of places, so please don't rely on what I've put in here - there's every chance I've made other typos.