/* rotating beam */
options results

if ~show('ports', "AMI2D") then do
    address command 'run ami2d:bin/ami2d'
    do while ~show('ports', "AMI2D")
    end
end

address ami2d

call 'ami2d:rexx/plstress'

if ~show('l', "rexxmathlib.library") then do
    check = addlib('rexxmathlib.library',0,-30,0)
end

h = 0.5
L = 10.0
E1 = 1000
nu1 = 0.3000

reset
iso 1 E1 nu1
do i=0 to 20
    xx = L*i/20
    node i+1 xx 0.0
    node i+31 xx h
end
do i=0 to 19
    quad4 i+1 1 i+1 i+2 i+32 i+31
end

dispx 1 0
dispy 1 0
alias 'rotx=pbc(a1,1,0,1,10*cos(a2)-0.5*sin(a2)-10)'
alias 'roty=pbc(a1,1,0,2,10*sin(a2)+0.5*cos(a2)-0.5)'
rotx 51 45
roty 51 45
exit