小球撞击平板.doc

学习CAE已有很长一段时间，在simwe论坛也学到了很多，但在论坛里一直都是活动在dyna版，鉴于我本身的学习的前处理为vpg，故将制作一些算例，以此来回报simwe，并提高eta版块的关注度，vpg还是挺强大的。
5 s& w) G3 ]& F0 |0 e- E
4 G) Q& ?. B* ?  X7 `- j7 v4 G提前声明：本人做的算例是些简单应用，目的是帮助新手了解vpg，慢慢熟悉vpg的过程，请高手不要鄙视！谢谢！如果有时间我会做些复杂供新手参考~$ T( [6 o5 s. w7 s: u9 Y
( n) _8 z# @4 c9 ~+ k6 z
, E4 j, q# J' c9 }5 @( d
小球撞击平板模型介绍：这个算例相信在很多CAE书里面都有，也有很多是别人的做过。我平时一般的建模都是从UG里建，然后导出IGES文件，并将其导入vpg做前处理，而我的这个模型的建模完全是在VPG里完成的，目的就是让新手熟悉vpg而已。: k% k3 n0 g* a; ~+ L1 K2 B. h( e
尺寸：小球直径20mm；平板：50*50mm9 K- r1 `) F- a# b6 @
完成后效果图：
下载 (183.49 KB)
2010-7-31 09:47

! _1 O: {3 \3 U7 b& d4 f4 W0 ~% v: o' D
建模过程如下：
/ \& z+ }4 R+ p' F2 R1 t1、首先建立两个part，分别命名为plate和ball。
, Z1 ?7 K$ m: @/ l2、将plate置为当前part，如图：
2 ^+ Q0 z6 C$ o% f% `0 K
下载 (12.46 KB)
2010-7-31 09:47
7 T. @6 v' @+ ^5 n/ {- R8 W' z2 e: i
3、在pre——line中选择generate lines 即新建线，如图：  A0 g4 q% n, t$ R$ e. P4 ~

下载 (18.29 KB)
2010-7-31 09:47

s- v4 |4 L2 S, b4、通过坐标建立四条直线：& {7 L. [  I2 R- B
（-25,-25,0）（25,-25,0）
; ~, [* W4 |) ^; _3 h  @: s! q( O/ t（25,-25,0）（25,25,0）6 E& j* b: b  ?9 Y/ L
（25,25,0）（-25,25,0）2 m7 E3 a2 B, z) I5 W3 P
（-25,25,0）（-25,-25,0)
* ^5 u; W; h8 c+ M2 Z8 [5、进入element划分网格：选择4 lines mesh，以此选择上述四条直线，输入划分数量：50,50.如图：) S% o6 L0 ?7 m: m4 z

下载 (12.01 KB)
2010-7-31 09:47
完成plate的网格划分；' Y! F6 ~3 C4 C7 R* A# [
6、将ball置为当前part，选择element中的ball mesh，进行网格划分：直径20mm，离地间距：20mm，如图：
: w+ `; T7 a2 y5 W3 t" E- o
下载 (182.64 KB)
2010-7-31 09:47
，输入数值后，网格划分完毕。
0 E$ \7 L, C# u! Q& E0 P  u* y
1 @5 S! ]- ^, {下面为材料定义、属性定义，以及边界条件定义、接触定义。" l$ O! _# l  R$ r: b. V
7、定义材料material：将小球定义为刚体，只释放Z方向的自由度：com=1，con1=4，con2=7；平板定义为弹性材料，取默认值；
2 b+ W  p) A  B! T4 a2 ^4 U8、单元属性定义element property：新建shll属性，将积分厚度ts1~ts4均设为1mm，并将其赋给上述两个part；
$ l5 W5 W0 N: Y9 H: ~9，边界条件：将平板的四个角固定：boundary condition中定义SPC，选择CREATE，node，将六个自由度都选中：
y. K" D2 t1 x/ l& u+ Z
下载 (13.86 KB)
2010-7-31 10:05
点击OK后选择plate的四个角，完成边界定义. n2 V4 a) T1 b, `1 F7 Y  _

下载 (73.96 KB)
2010-7-31 10:05
( |# B% R6 \" l! z8 x) S( g: U
10、定义node set ：选择set menu，中的node，将小球定义到node set，如set 1中，在下面的初速度定义，以及接触定义要用的；$ V7 h& L" z7 Z# J( w; U
11、定义初始速度：选择initial condition中的velocity定义小球的初速度：create选择node set，初速度：Z=1000，即：1m/s；
! A& Q1 n. r& o& u/ \: h  H12、定义接触：选择contact，create定义接触为：automatic nodes to surface，将plate定义为master，小球定义为slave；; J7 K- X- g" Q! Q5 N$ c' L& d
13、控制菜单定义：选择analysis中的control card定义控制关键字：& q( P1 X! k5 H$ D8 O( C+ w2 d. ^6 g& J
*CONTROL_CONTACT0 ^: h/ P( f4 p7 M, A4 E/ e! ^1 @
$   slsfac    rwpnal    islchk    shlthk    penopt    thkchg     orien    enmass
" c" F  B6 `9 w3 F/ ~8 ^% `   0.10000       0.0         1         0         1         0         1         0
  y1 u) t% {+ w- E6 u; O* K. o" j$   usrstr    usrfrc     nsbcs    interm     xpene     ssthk      ecdt   tiedprj1 p5 H' q3 J2 J" M! O  A
         0         0        10         0     4.000         0         0         0
9 [; L0 Y% ~4 k. N1 r$    sfric     dfric       edc       vfc        th     th_sf    pen_sf         
, R1 p& G2 I7 A6 ]6 n+ t$ r       0.0       0.0       0.0       0.0       0.0       0.0       0.0          $ M5 y; N. S1 K0 H/ t& f
$   ignore    frceng   skiprwg    outseg   spotstp   spotdel                    5 l* o$ D: I7 D9 @! a
         0         0         0         0         0         0                    
2 ?* j6 ^7 I5 O3 _, o+ d*CONTROL_TERMINATION" C+ i& g1 x* t  Y; {5 C% ^' l  S
$   endtim    endcyc     dtmin    endeng    endmas                              , P, u. w8 ~, ^3 c
   0.02000         0       0.0       0.0       0.0                              ' L9 _. g  A& z
14、定义输出关键字：database ASCII0 @" x2 [, r  r+ ^" ^/ ?

*DATABASE_GLSTAT, b9 x" g1 h( L
$       dt    binary                                                            
. A1 M7 c/ N8 r1 g9 \6 J6 D1 T  .1000E-2         1                                                            0 ?9 I( m8 G; T0 C' b" p4 e! W
*DATABASE_MATSUM. G; J0 \; P6 c: s7 N+ \# g, p- S
$       dt    binary                                                            0 D! ?, ?' [8 y% ^2 g( Z7 W( P
  .1000E-2         0                                                            
3 ]; Z, m& G* g- c1 P6 O/ a! G15、定义关键字：database binary  Z3 ]* w5 |9 V' m+ u
*DATABASE_BINARY_D3PLOT7 G$ R0 B, Q' `- \9 R
$  dt/cycl      lcdt      beam     npltc                                       
$ L1 l! G# N! K! H$ N  .1000E-2         0         0         0                                        & G: t: I* j, d" ~, ]# ]
$    ioopt                                                                      0 ]- \7 ~/ e. y/ Y
         0                                                                      ) I  V. p) _4 Q4 p
16、最后保存，并将其export，导入dyna中计算即可。1 J0 K4 G" m! F, }" T
说明：对于控制关键字可根据需要定义。