关于铁道车辆冲击out_of_range velocities

lqs10620

由于第一次实作比较大的模型，在计算的过程中（大概14ms）一直会出现30个NODE out_of_range velocities。我用ls-prepost观察了动画，显示在计算终止前有部分单元穿透刚性墙，然后接着丢失大片单元。我分析了下原因如下，期待得到您的指导。谢谢。
目标是建立半车模型施加初速度碰撞固定刚性墙。
1）单位制应该没有问题，我使用的是KG/MM/MS 应力单位输入Gpa
2）我在计算中使用了质量缩放，时间步改为0.67 控制了沙漏（最后会附上部分K文件）请指导。
3）整个模型有5000多个面组成，单元为shell163，壳单元设置均为默认。但处理连接关系的时候采用静力分析中的GLUE和overlap形式，没有预留壳单元偏置厚度。
由于这个原因造成了初始穿透有200多处，然后我在接触控制部分设置了ignore参数为1（忽略初始穿透）由于重新建模工作量很大，我想请问我是不是需要重新建模，然后在壳单元连接的地方预留1/2壳厚度的距离，然后用焊点约束定义连接关系？两个很长的面（8m）是不是也只需要定义两个焊点？
4）在模型中没有定义重力加速度。刚性墙处理为在ANSYS中先定义了四个节点定义单元属性生成单元，在K文件中刚性强定义中采用了这个节点。
5）网格划分控制的很均匀，已尽量采用映射网格划分技术，但有些单元为了处理其连接关系划分的比其他单元尺寸小很多。

材料模型关键字
*MAT_PLASTIC_KINEMATIC
         1 0.780E-05  207.      0.300000 0.345      0.00      0.00   
  0.00      0.00      0.00   
*MAT_RIGID
         2 0.700E-05  207.      0.300000       0.0       0.0       0.0
  1.00      7.00      7.00
参数控制关键字
*SET_SEGMENT
         2     0.000     0.000     0.000     0.000
    105000    105001    105002    105003     0.000     0.000     0.000     0.000
*CONTACT_AUTOMATIC_SURFACE_TO_SURFACE
         1         2         0         0         0         0         0         0
0.1000    0.1000     0.000     0.000     0.000             0 0.000    0.1000E+08
1.000     1.000     0.000     0.000     1.000     1.000     1.000     1.000   
*CONTACT_AUTOMATIC_SINGLE_SURFACE
$     SSID      MSID     SSTYP     MSTYP    SBOXID    MBOXID       SPR
         0         0         5         0                             0     
$       FS        FD        DC        VC       VDC    PENCHK        BT      
      0.20      0.20       1.0       0.0       0.0         0       0.0      
$      SFS       SFM       SST       MST      SFST      SFMT       FSF      
       0.0       0.0       1.0       0.0       0.0       0.0       1.0      
$     SOFT    SOFSCL    LCIDAB    MAXPAR     SBOPT     DEPTH     BSORT   
         1       0.0                 0.0      0.50       0.0         0   
$   PENMAX    THKOPT    SHLTHK     SNLOG      ISYM     I2D3D    SLDTHK    SLDSTF
       0.0         0         0         0         0         0                    
$     IGAP    IGNORE
         2         0   
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              CONTROL OPTIONS                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
*CONTROL_HOURGLASS
$      ihq        qh                                                            
         4   0.05000
*CONTROL_BULK_VISCOSITY
$       q2        q1                                                            
     1.500   0.06000
*CONTROL_CONTACT
$   SLSFAC    RWPNAL    ISLCHK    SHLTHK    PENOPT    THKCHG     ORIEN    ENMASS
      0.10       2.0         0         0         0         0         0         
$   USRSTR    USRFRC     NSBCS    INTERM     XPENE     SSTHK      ECDT   TIEDPRJ
         0         0        10         0       4.0         1         0         0
$    SFRIC     DFRIC       EDC       VFC        TH     TH_SF    PEN_SF
       0.0       0.0       0.0       0.0       0.0       0.0
$   IGNORE    FRCENG   SKIPRWG    OUTSEG   SPOTSTP   SPOTDEL
         1         0                                                              
*CONTROL_ENERGY
         2         2         2         2
*CONTROL_SHELL
$   wrpang    itrist     irnxx    istupd    theory       bwc     miter
  30.0             1        -1         1         2         2         1
*CONTROL_TIMESTEP
    0.0000    0.6700         0  0.00    -0.2100E-3   
*CONTROL_SUBCYCLE
*CONTROL_TERMINATION
  80.0             0   0.00000   0.00000   0.00000
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                                 TIME HISTORY                                 $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_BINARY_D3PLOT
0.8000   
*DATABASE_BINARY_D3THDT
0.4000   
*DATABASE_BINARY_D3DUMP
   3851473
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                              DATABASE OPTIONS                                $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*DATABASE_EXTENT_BINARY
         0         0         3         1         0         0         0         0
         0         0         4         0         0         0
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$                         INITIAL VELOCITY DEFINITIONS                         $
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
*INITIAL_VELOCITY
$     nsid    nsidex     boxid
                             1
$       vx        vy        vz      
      10.0       0.0       0.0
$
*DEFINE_BOX
$    boxid       xmm       xmx       ymn       ymx       zmn       zmx
         1      -100     11025        -5      1545      -850      2600
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$$$$  Rigidwalls
$
$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$
$
$...>....1....>....2....>....3....>....4....>....5....>....6....>....7....>....8
$  固定的、有限大小的刚性墙
*RIGIDWALL_PLANAR_FORCES
$     nsid    nsidex     boxid
         0         0         0
$
$       xt        yt        zt        xh        yh        zh      fric
   11030.0       0.0       0.0       0.0       0.0       0.0     0.000
$
$     soft      ssid      nid1      nid2      nid3      nid4
        10         0    105000    105001    105002    105003
$
$
*END

刚性墙nsid参数由于可能与刚性墙接触的节点比较多，不知道这样处理是否合适。

sillyoranger

可能与偏置1/2壳厚无关 我预留了1/2壳厚 仍出现你的问题
仅为参考

lqs10620

请问楼上你预留了之后怎么连起来的？

lqs10620

斑竹达人  过来解决这个问题

lqs10620

没人讨论不热烈啊

lqs10620

我期待的斑竹达人一直没有出现