/* test of a transformed orthotropic material */

options results

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

address ami2d

reset

call 'ami2d:rexx/cstrain'

tiso 1 181.0 10.3 0.23 0.70 7.17

lam 2 0 1 1/2 ' 45.0'
lam 2 0 1 1/2 '-45.0'

coord 'polar'

do r=10 to 10.25 by 0.125
    do ang=90 to 0 by -9
        node 0 r ang
    end
end

do j=0 to 4
    i = j*2
    'quad8('j + 1',2,'i + 1','i + 3','i + 25','i + 23','i + 2','i + 14','i + 24','i + 12')'
    'press('j + 1',1,-0.1)'
end
'dispx(1,0.0)'
'dispx(12,0.0)'
'dispx(23,0.0)'
'dispy(11,0.0)'
'dispy(22,0.0)'
'dispy(33,0.0)'
exit