BjornF
Mechanical
- Mar 18, 2010
- 9
Hi, everyone
I am trying to model a 2D cell (square) of fibres in shear. To do so, I am applying a prescribed displacement at the topside nodes, clamp the nodes at the bottomside nodes, and prescribe corresponding displacements at the sides of the cell.
The problem is, in order to get the right shape, I need to couple the topside nodes in the vertical direction to ensure the correct shape.
Here is the short code-snippet I have written about the BCs:
--------------------------------
! X = Xsheet, 0 - Apply BCs on side of the cell
nsel,s,loc,x,0
nsel,a,loc,x,Xsheet
*get,nnodes,node,0,count ! Count number of nodes in set
*get,nodenumber,node,0,num,min ! Get initial node number
*do,j,1,nnodes
*get,nodepos,node,nodenumber,loc,y ! Get nodal position in y
d,nodenumber,ux,(nodepos/(fiberL*10))*shear_strain ! Set displacement
*get,nodenumber,node,nodenumber,nxth ! Get next node number
*enddo
! Y = 0
nsel,s,loc,y,0
d,all,all
! Y = Ysheet
nsel,s,loc,y,Ysheet
d,all,ux,shear_strain
cp,1,uy,all
------------------------
However, the vertical displacements are not coupled for the selected nodes. Is there something in this code that overrides the CP command?
Thanks for any input
Björn
I am trying to model a 2D cell (square) of fibres in shear. To do so, I am applying a prescribed displacement at the topside nodes, clamp the nodes at the bottomside nodes, and prescribe corresponding displacements at the sides of the cell.
The problem is, in order to get the right shape, I need to couple the topside nodes in the vertical direction to ensure the correct shape.
Here is the short code-snippet I have written about the BCs:
--------------------------------
! X = Xsheet, 0 - Apply BCs on side of the cell
nsel,s,loc,x,0
nsel,a,loc,x,Xsheet
*get,nnodes,node,0,count ! Count number of nodes in set
*get,nodenumber,node,0,num,min ! Get initial node number
*do,j,1,nnodes
*get,nodepos,node,nodenumber,loc,y ! Get nodal position in y
d,nodenumber,ux,(nodepos/(fiberL*10))*shear_strain ! Set displacement
*get,nodenumber,node,nodenumber,nxth ! Get next node number
*enddo
! Y = 0
nsel,s,loc,y,0
d,all,all
! Y = Ysheet
nsel,s,loc,y,Ysheet
d,all,ux,shear_strain
cp,1,uy,all
------------------------
However, the vertical displacements are not coupled for the selected nodes. Is there something in this code that overrides the CP command?
Thanks for any input
Björn