!   trw_rotor_2d.mac
!
rotor_cover_315a139_2d,1
rotor_sleeve_375a108_2d
!
init_entity_num,1
aglue,all
init_entity_num,1
!
!
*ask,pload_,Pressure Load,1
*ask,slv_,Solve this model now (True=1),0
!
/title,Rotor 551A182 @ %pload_% psi External
!
view,'f'
!
asel,s,area,,am-2
msh_area,0,3,1,0.003,0.003
init_entity_num,1
!
asel,s,area,,am
msh_area,0,3,1,0.010,0.010
init_entity_num,1
!
eplot
!
!
*if,slv_,eq,1,then
    lsel,s,line,,lm-9
    nsll,s
    d,all,uy
    init_entity_num,1
    !
    lsel,u,line,,lm-30
    lsel,u,line,,lm-43
    ln_group,km-2,km-3,km-10
    !
    sfl,all,pres,pload_
    !
    init_entity_num,1
    !
    finish
    /config,fsplit,250
    /solu
    eqslv,iter,5,3
    solve
    !
    finish
    !
    /post1
    /eface,1
    avprin,0,0,
    !
    plnsol,s,eqv,0,1
    /eface,1
    avprin,0,0,
    !
    *get,Smax,plnsol,,max
    !
    round_num,_Sy/Smax,3
    SFy=num_                ! calculated safety factor @ yield
    round_num,_Su/Smax,3
    SFu=num_                ! calculated safety factor @ ultimate
    !
    round_num,(_Sy/(_FSyR*Smax))-1,3
    MSy=num_                ! calculated margin of safety @ yield
    round_num,(_Su/(_FSuR*Smax))-1,3
    MSu=num_                ! calculated margin of safety @ ultimate
    !
    round_num,Smax,0
    Smax=num_
    !
    s=Smax
    p=pload_
    y=SFy
    u=SFu
    !
    /title,%p% psi External,  Smax=%s%,  SFy=%y%,  MSy=%MSy%,  SFu=%u%,  MSu=%MSu%
    /post1
    /eface,1
    avprin,0,0,
    !
    view,'f'
*endif