Reply To: How to apply displacement derivative dependent force on bunch of nodes in apdl?

[Mehdi Baghaee]

Since I am doing transient analysis, I am looking for the simplest way (commands) in apdl to apply the load formulation (P=dW(x,y,z,t)/dx+dW(x,y,z,t)/dt) on a node!

For example the node 11 in the following cantilevered beam.

/BATCH  
/input,menust,tmp,''
WPSTYLE,,,,,,,,0
/BATCH  
/input,menust,tmp,''
WPSTYLE,,,,,,,,0
/NOPR   
KEYW,PR_SET,1   
KEYW,PR_STRUC,1 
KEYW,PR_THERM,0 
KEYW,PR_FLUID,0 
KEYW,PR_ELMAG,0 
KEYW,MAGNOD,0   
KEYW,MAGEDG,0   
KEYW,MAGHFE,0   
KEYW,MAGELC,0   
KEYW,PR_MULTI,0 
/GO 
/PREP7  
ET,1,BEAM188
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,0
KEYOPT,1,4,0
KEYOPT,1,6,0
KEYOPT,1,7,0
KEYOPT,1,9,0
KEYOPT,1,11,0   
KEYOPT,1,12,0   
KEYOPT,1,15,0   
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,1e9
MPDATA,PRXY,1,,0.3  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2000 
SECTYPE,   1, BEAM, RECT, , 0   
SECOFFSET, CENT 
SECDATA,0.01,0.01,0,0,0,0,0,0,0,0,0,0   
N,1,0,0,0,,,,   
N,2,0.1,0,0,,,, 
N,3,0.2,0,0,,,, 
N,4,0.3,0,0,,,, 
N,5,0.4,0,0,,,, 
N,6,0.5,0,0,,,, 
N,7,0.6,0,0,,,, 
N,8,0.7,0,0,,,, 
N,9,0.8,0,0,,,, 
N,10,0.9,0,0,,,,
N,11,1,0,0,,,,  
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,   1 
TSHAP,LINE  
FLST,2,2,1  
FITEM,2,1   
FITEM,2,2   
E,P51X  
FLST,2,2,1  
FITEM,2,2   
FITEM,2,3   
E,P51X  
FLST,2,2,1  
FITEM,2,3   
FITEM,2,4   
E,P51X  
FLST,2,2,1  
FITEM,2,4   
FITEM,2,5   
E,P51X  
FLST,2,2,1  
FITEM,2,5   
FITEM,2,6   
E,P51X  
FLST,2,2,1  
FITEM,2,6   
FITEM,2,7   
E,P51X  
FLST,2,2,1  
FITEM,2,7   
FITEM,2,8   
E,P51X  
FLST,2,2,1  
FITEM,2,8   
FITEM,2,9   
E,P51X  
FLST,2,2,1  
FITEM,2,9   
FITEM,2,10  
E,P51X  
FLST,2,2,1  
FITEM,2,10  
FITEM,2,11  
E,P51X  
FLST,2,1,1,ORDE,1   
FITEM,2,1   
/GO 
D,P51X,ALL, , , , , , , , , ,  
FINISH  
/SOL
ANTYPE,4
TRNOPT,FULL 
LUMPM,0 
DELTIM,0.01,0,0 
OUTRES,ERASE
OUTRES,ALL,1
TIME,1  
FLST,2,1,1,ORDE,1   
FITEM,2,11  
FLST,2,1,1,ORDE,1   
FITEM,2,10  
IC,P51X,UZ,0.1,0.2, 
SOLVE   
FINISH