kill(all);

assume(G>0,K>0,Cw>0,E>0);

here :
-- lowercase l is the overall, twice clamped beam HALF-LENGTH, 
   i.e. the length of an equivalent clamped-free beam
-- uppercase L is the overall, twice clamped beam length, 

L = 2l = eta*d

L:eta*d;

l:L/2;

here:
-- lowercase "theta" is the relative rotation between the mispan section, and the two terminals
-- uppercase "Theta" is the relative rotation between the two terminals

Theta = 2 theta

cross-section characteristic length with respect to Vlasov theory of restrained warping

Restrained warping effect is negligible if L/d >> 1 .

tentative solution

psi(z):=A*sinh(z/d)+B*z/d;

tentative solution check

fullratsimp(
    -d^2*diff(psi(z),z,4)+diff(psi(z),z,2)=0
);

imposed counterotation at the two terminals

eq1 : psi(l)-psi(-l)=Theta;

imposed zero warping at terminals

diff(psi(z),z,1);
eq2 : ev(%=0,z=l);

solution

linsolve([eq1,eq2],[A,B]),globalsolve=true;

torsional moments, De Saint Venant and Vlasov portions

T_dSV :         G*K*diff(psi(z),z,1);
T_Vla : - (d^2)*G*K*diff(psi(z),z,3);

plotme:fullratsimp(
[
    ev(T_Vla / (G*K*Theta/L),z=  l,infeval),
    ev(T_Vla / (G*K*Theta/L),z=  0,infeval),
    ev(T_dSV / (G*K*Theta/L),z=  0,infeval),
    1
]
);
plot2d(
    plotme,
    [eta,0.5,100],
    logx,
    [legend,"T_{tot} = T_{Vla}, ext.","T_{Vla}, midspan","T_{dSV}, midspan","1"], 
    [xlabel, "2l / d_{Vla}"],
    [ylabel, "(T * 2l)/(G * Kt * [2 theta])"],
    [title,"Torsional stiffness, "]  
);

B:-G*K*d^2*psi(z);

diff(B/Theta/G/K,z,2)$
ev(%,sols)$
ev(%,z=L/2);
wxplot2d(
    %,
    [eta,1,100],
    logx,logy,
    [legend,"B/(G * Kt * [2 theta])"], 
    [xlabel, "2l / d_{Vla}"],
    [ylabel,"B/(G * Kt * [2 theta])"],
    [title,"Scaling factor for the Vlasov induced stress comps at the end"] 
);

Definition of a stiffening factor due to the Vlasov effect

T : ev(T_Vla,z=l);

( T/(G*K*Theta/L));
Amp_Vla : ratsubst(tanh(eta/2)*cosh(eta/2),sinh(eta/2),%);

Newton-Raphson procedure for estracting the 

d/L = eta 

ratio from the Vlasov stiffening factor 

residual definition; 

1/x = 1/Amp_Vla

such 1/Amp_Vla function is well defined on the 0.001 <= eta <=1000 range

wxplot2d(1/Amp_Vla,[eta,0.001,1000],[logx]);

wxplot2d(1/Amp_Vla,[eta,0.001,1000],[logx],[logy]);

definition of a Newton procedure:
residual function and its derivatives

res : 1/x-1/Amp_Vla;

define(f(x,eta) ,
    expand(res)
);

define(
        df(x,eta),
        expand(diff(f(x,eta),eta))
    );

define(
    NR(x,eta),
    fullratsimp(eta - f(x,eta)/df(x,eta))
);

stiffening effect due to Vlasov

l=960mm, L=1920mm

amp: [restrained warping stiffness] / [free warping stiffness]

OTW
amp: (617330*2) / (27391.3*2);

CTW
amp: (1.76728e7*2) / (1.76127e7*2);

amp: (1.76728e7*2) / (1.76127e7*2);

1.0;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;
NR(amp,%),numer;

d:1920/%;