(*             CONICAL BAR IN TENSION             *)

ϕ=A*Log[R*(1+Cos[β])];
ω=B/R;
ψ=0;

SetDirectory[NotebookDirectory[]];
<<ABErtb.m;

c1=σβR/.β->β0
c2=σββ/.β->β0
t=σRR*Cos[β]-σβR*Sin[β]/.R->Rb

T=Simplify[Integrate[2*π*t*Rb^2*Sin[β],{β,0,β0}]]

Solve[{c1==0,c2==0,T==F},{A,B}]

R=Sqrt[r^2+z^2]
ϕ=A*Log[R+z]/.%9
ω=B/R/.%9

<<ABErtz.m;

f1=σzz*π*Tan[β0]^2/F/.{r->x*Tan[β0],z->1,ν->0.3};
f2=f1/.β0->π/3;
f3=f1/.β0->π/6;
f4=f1/.β0->π/12;
f5=f1/.β0->π/180;

Plot[{f2,f3,f4,f5},{x,0,1},PlotRange->{0,6}]