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)