【求助】如何定义角速度

zll2001

我用ls－dyna模拟车轮旋转并前进的过程，如何给车轮上的节点加上角速度？
请指教具体操作过程。

hekie

命令omega,x,y,z(x,y,z分别是相对整体坐标的三个方向的角速度值)

zll2001

是的，我用如下定义了车轮节点的角速度：
EDLOAD,ADD,omgx,,ROLL,TIME1,ROTZ1
但模拟的结果车轮不但没有旋转，而且单元发生了畸变，
好像车轮被压扁了一样变大、变薄了。

willow

Roll這個component 是不是包含旋轉物件的所有node
也許應該試試只加在中心點

zll2001

是的，ROLL是包括旋转体上的所有节点。
为什么只加在中心点上，如果找不到中心点位置呢？
而且中心点指什么位置？应该是中心线吗？

liu02

是不是你的约束条件不合理？

zll2001

不是，除了简单的对称面约束，只有加速度一种载荷

liu02

可能小鸟撞击螺旋桨的例子会对你有所启发，在ansys/ls-dyan的2张cd里有这个例子，找到后麻烦给我寄一份:liuchaoxun@sohu.com

willow

似乎加載於非rigid body的角加速度皆會換算成離心力施加於其上
所以分析的過程看不到物體自轉的情況
只看得到力的作用不知道是不是這樣

zll2001

willow,如果象您所说的这样，对于非rigid body如何实现使其旋转前进的状态。

willow

在part上貼上rigid element
或是在剛體與彈性體間切換
  
可以做得出來自轉與前進　但不知是否是正確的做法
我只做自轉　圖檔太大被我切了　原始可以轉完一圈

zll2001

我看到了您做的图象，您是否是在内圈贴了一层单元？能否教给我具体的操作过程。
另外，我也想到了一个与您想法类似的方法，但没有试验成功，请您指教问题大约是什么？我在内圈定义了一个刚性轴，使用shell163单元，与外圈弹性体glue到一起，然后定义刚性轴转动，但没能带动外圈弹性体一起转。

willow

用glue 就不對了
  
直接建立一個Cylinder 然後先做solid mesh
在使用shell element mesh前後任一面   
實作應該要避開會contact 的面
注意shell element 所使用的材料性質必須是 rigid  
然後把角速度家在shell element 的part上
就可以了

zll2001

很遗憾，我没有理解“在使用shell element mesh前後任一面實作應該要避開會contact 的面 ”，能否详细解释一下，或者将您上面做的题目的命令流发给我。zll_2001@fescomail.net，谢谢！
还有一个问题请教，用solid164划分阶梯轴不能实现，长以三角形单元代替，但这样会出一些警告，如何解决？

liu02

willow是高手。

zll2001

willow,难道除了加贴刚性单元的方法，弹性体就无法转动吗？！
真是郁闷。

willow

是不是只有這方法我不知道
命令流得整理一下才有  因為是用GUI做的
跟其他job混在一起

zll2001

那您能否重新解释一下“在使用shell element mesh前後任一面實作應該要避開會contact 的面”

willow

我的意思是
要選取圓柱的正面或反面貼上rigid element 都沒關係
但是如果此面會與其他part contact的話就不適合了
  
事實上  我也還在試  因為這個模型的hourglass energy粉大
還在想是什麼原因
我也覺得只能施加在剛體這個情況不太對

zll2001

以下是我做成功的一个试验题目，你可以看看
ET,1,SOLID164  
ET,2,SOLID164   
  
MPMOD,1,19   
MP,EX,1,210E3   
MP,DENS,1,7.82E-5   
MP,NUXY,1,0.3
TB,PLAW,1, , ,1  
TBDAT,1,0.286,   
TBDAT,2,0.00689,
TBDAT,3,1.0,
TBDAT,4, ,   
TBDAT,5, ,   
TBDAT,6, ,   
edmp,hgls,1,4,,,,,
edhgls,0.12
  
MP,DENS,2,0.0000782,   
MP,EX,2,210000,   
MP,NUXY,2,0.3,  
EDMP,RIGID,2,7,7  
  
BLOCK,0,1,0,1,0,10,
  
wpof,,2.000000
wpro,,,90.000000
wpof,-1.000000  
  
CYLIND,0.6,1,-0.2,0.3,0,360,   
CYLIND,0.6,1.2,-0.7,-0.2,0,360,   
   
FLST,2,2,6,ORDE,2   
FITEM,2,2   
FITEM,2,-3   
VADD,P51X   
  
TYPE,1
MAT,1
CM,_Y,VOLU   
VSEL, , , ,       4  
CM,_Y1,VOLU  
CMSEL,S,_Y   
!*   
CMSEL,S,_Y1  
VATT,       1, ,   1,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
!*   
ESIZE,0.5,0,
MSHAPE,1,3D  
MSHKEY,0
!*   
CM,_Y,VOLU   
VSEL, , , ,       4  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y   
!*   
VMESH,_Y1   
!*   
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2     
  
TYPE,2
MAT,2  
CM,_Y,VOLU   
VSEL, , , ,       1  
CM,_Y1,VOLU  
CMSEL,S,_Y     
CMSEL,S,_Y1  
VATT,       2, ,   2,       0   
CMSEL,S,_Y   
CMDELE,_Y   
CMDELE,_Y1   
FLST,5,1,4,ORDE,1   
FITEM,5,10   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
LESIZE,_Y1, , ,10, , , , ,1   
FLST,5,1,4,ORDE,1   
FITEM,5,6   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
LESIZE,_Y1, , ,2, , , , ,1     
FLST,5,1,4,ORDE,1   
FITEM,5,7   
CM,_Y,LINE   
LSEL, , , ,P51X  
CM,_Y1,LINE  
CMSEL,,_Y   
LESIZE,_Y1, , ,2, , , , ,1     
MSHAPE,0,3D  
MSHKEY,1   
CM,_Y,VOLU   
VSEL, , , ,       1  
CM,_Y1,VOLU  
CHKMSH,'VOLU'   
CMSEL,S,_Y     
VMESH,_Y1      
CMDELE,_Y   
CMDELE,_Y1   
CMDELE,_Y2  
/SOLU  
EDPART,CREATE
  
ASEL,S,,,9                  
ASEL,A,,,10
NSLA,R,1
CM,WORKSURF,NODE
ASEL,ALL
NSEL,ALL
  
FLST,2,1,5,ORDE,1   
FITEM,2,13   
!*   
/GO  
DA,P51X,UX,0
  
FLST,2,1,5,ORDE,1   
FITEM,2,3   
/GO  
DA,P51X,UX,0
  
FLST,2,1,5,ORDE,1   
FITEM,2,3   
/GO  
DA,P51X,Uy,0
  
FLST,2,1,5,ORDE,1   
FITEM,2,3   
/GO  
DA,P51X,Uz,0
  
FLST,2,1,5,ORDE,1   
FITEM,2,2   
/GO  
DA,P51X,UZ,0
FLST,2,1,5,ORDE,1   
FITEM,2,1   
/GO  
DA,P51X,UZ,0
  
VSEL,S,,,4
NSLV,R,1
CM,roll,NODE
VSEL,ALL
NSEL,ALL
*DIM,TIME1,,2                              
*DIM,ROTz1,,2
TIME1(1)=0,4.00001
ROTz1(1)=6.28,6.28
cgloc,0.3,2,1
edload,add,omgx,,roll,time1,rotz1
!!定义摩擦                           
EDCGEN,sTS,WORKSURF,2,0.45,0.45
EDCONTACT,0.1,,2,2,4,1,2,
!定义车轮重量
NSEL,r, , ,     117  
CM,zh,NODE  
NSEL,ALL
vsel,all
*DIM,TIME2,,2                              
*DIM,ROTz2,,2
TIME2(1)=0,4.0000
ROTz2(1)=-0.14,-0.14
!EDLOAD,ADD,ay,,zh,time2,rotz2
!///////////
TIME,1,
EDRST,40,   
EDHTIME,40,