position p (p ∈ {1,2} 解释

yiqing

nd_pos( ND_PTR np, INT p, INT dof ),position p (p ∈ {1,2} denotes current or reference position, respectively;
请问(p ∈ {1,2},这个是什么意思哦

yiqing

sp = s_head
loop while sp # null
        if s_type(sp) = 1 then
                ;第1节点相关数据
                ii=ii+1
                np = s_node(sp,1)
                 xx = nd_pos(np,2,1)
                 yy = nd_pos(np,2,2)
                zz = nd_pos(np,2,3)； 里面2是什么意思？
                FFx= sb_force(sp,1,1)                       
                    FFy= sb_force(sp,1,2)  ；里面的1又是什么意思？                     
                    FFz= sb_force(sp,1,3)                       
                    mmx= sb_mom(sp,1,1)                        
                    mmy= sb_mom(sp,1,2)                        
                    mmz= sb_mom(sp,1,3)
                xtable(9,ii)=xx
                ytable(9,ii)=ffy
               
        endif
        sp = s_next(sp)
endloop
end

ycding

当前座标和参考坐标，参考座标是不会变的，大变形模式下当前座标是会变的。

yiqing

sb_force( bp, e, d )force at end e of beamSEL (e ∈ {1,2}; d-component, d ∈ {1,2,3}). These forces are expressed in terms of the beamSEL coordinate system — see the figure below for sign convention.
sb_force，这里面 e好像不是当前座标和参考坐标，与ND_PTR np, INT p, INT dof ),position p (p ∈ {1,2这个意义不同

yiqing

;---------锚索类单元结构内力输出子程序------------------------
def force_outp
xtable(9,1) = 0.0
ii=0
sp = s_head
loop while sp # null
        if s_type(sp) = 2 then
                ;第1节点相关数据
                ii=ii+1
                np = s_node(sp,1)
                 xx = nd_pos(np,2,1)
                 yy = nd_pos(np,2,2)
                zz = nd_pos(np,2,3)
                FFx= sc_force(sp)                       
                     
                xtable(9,ii)=xx
                ytable(9,ii)=ffx
                ;第2节点相关数据
               
        endif
        sp = s_next(sp)
endloop
end

force_outp

set logfile fy.txt
set log on
set pagelength 10000
pri table 9
quit