Why doesn't coriolis moment with Solid elems agree with hand calculations?
|
The attached input for the solid185 version of this test gives the expected value of MY 116. It uses a transient analysis and the stationary reference frame. Development suggested that the interpretation of the results from a static rotating reference frame was not clear and that the stationary frame should be used whenever possible. fini /clear md = 2 /PREP7 wpro,,,90 CYL4,0,0,10,,,,100 KWPLAN,-1,3,7,5 VSBW,1 LESIZE,9,,,10*2*md,,,,,1 ET,1,185,,2 ESIZE,8/md VSWEEP,ALL mp,ex,1,210000 mp,prxy,1,.3 mp,dens,1,7.8e-9 nsel,s,loc,x,-.001,.001 ! nsel,s,loc,x,50 d,all,all,0 alls CM,ROTSHAFT,ELEM *get,max_node,node,,num,max /solu ! antype,static ! coriolis,on antyp,trans coriolis,on,,,on time,1 delt,.5 cmomega,rotshaft,1000,0,0 *do,i,1,max_node ic,i,ux,,-10*ny(i) ! need this or moment is half expected value ic,i,uy,,10*nx(i) *enddo alls fcentrif,off solve !!!!!!!!! MY (THEORY) = wx*wz*Ixx = 1000*10*0.0116 = 116 /post1 set,last nsel,s,d,u,0 nforce |
||
|