Matriz de Rigidez Elasto-Plastica

Boa Tarde,
Queria fazer um programa de Analise Estática não Linear Geométrica e Material. Para isso necessitava de uma Matriz de Rigidez Elasto-Plástica de uma Viga Bi-encastrada com as coordenadas (N1,Q1,M1,N2,Q2,M2).
Alguém me pode ajudar sobre este assunto. Uma folha PDF por exemplo, uma referência.

Bom Dia,
Consegui arranjar na Net uma Matriz de Rigidez Elasto-Plástica.
Agora só Falta Calcular as Propriedades Geométricas elasto-plasticas, D1m, D2m, D3m.

KMep extremidade (rigido-rigido)

D1m/Lr;0;D2m/Lr;-D1m/Lr;0;D2m/Lr

Simétrica;12*D3m/(Lr*Ld^2);6*D3m/(Lr*Ld);0;-12*D3m/)Lr*Ld^2);6*D3m/(Lr*Ld)

Simétrica;Sim;4*D3m/Lr;D2m/Lr;-6*D3m/(Lr*Ld);4*D3m/Lr*1/2

Simétrica;Sim;Sim;D1m/Lr;0;-D2m/Lr

Simétrica;Sim;Sim;Sim;12*D3m/(Lr*Ld^2);-6*D3m/(Lr*Ld)

Simétrica;Sim;sim;Sim;Sim;4*D3m/Lr

D1m,D2m,D3m – propriedades geométricas Elasto-plásticas
Dkm=0,5*(Dka+Dkb) com 1<=k<=3
Dkj=Σ[Dij*(yci)^(k-1)*dAn] indice j (a ou b)