naturepointing
Marine/Ocean
- Mar 5, 2014
- 4
Hello, everyone.
Here I am trying to loop to add some geometrical imperfection to my components by name index, but it does not work and I can not see any deflection from the element. and I do not know why .
modmsh,detach ! detach FE-model from geometry needed to move nodes
pi=acos(-1.0)
*dim,p1node,,100000
*do,j,1,ns
allsel
*do,i,1,3
nump=99+i+14*j
nsel,s,node,,plate_%99+i+14*j%
*get,nmax,node,plate_%99+i+14*j%,count
*get,minno,node,plate_%99+i+14*j%,num,min
n0=minno
p1node(1)=n0
*do,m,2,nmax !to get every node number in the selected plate
n0=ndnext(n0)
p1node(m)=n0
*enddo
*do,in,1,nmax ! start loop over all nodes in one plate
*get,xno,node,p1node(in),loc,x
*get,yno,node,p1node(in),loc,y
*get,zno,node,p1node(in),loc,z
x0=100+(j-1)*a
dy= wopl*(b1*sin(pi*(xno-x0)/a) !all these parameter has been given in the previous code
n,p1node(in),xno,yno+dy,zno
*enddo
*enddo
*enddo
Here I am trying to loop to add some geometrical imperfection to my components by name index, but it does not work and I can not see any deflection from the element. and I do not know why .
modmsh,detach ! detach FE-model from geometry needed to move nodes
pi=acos(-1.0)
*dim,p1node,,100000
*do,j,1,ns
allsel
*do,i,1,3
nump=99+i+14*j
nsel,s,node,,plate_%99+i+14*j%
*get,nmax,node,plate_%99+i+14*j%,count
*get,minno,node,plate_%99+i+14*j%,num,min
n0=minno
p1node(1)=n0
*do,m,2,nmax !to get every node number in the selected plate
n0=ndnext(n0)
p1node(m)=n0
*enddo
*do,in,1,nmax ! start loop over all nodes in one plate
*get,xno,node,p1node(in),loc,x
*get,yno,node,p1node(in),loc,y
*get,zno,node,p1node(in),loc,z
x0=100+(j-1)*a
dy= wopl*(b1*sin(pi*(xno-x0)/a) !all these parameter has been given in the previous code
n,p1node(in),xno,yno+dy,zno
*enddo
*enddo
*enddo