【转载】ANSYS 动力分析 (9) - 瞬态动力分析 (1)

原文地址：http://htbbzzg.blog.163.com/blog/static/69725206201081663611208/

第四章   瞬态动力分析

    第一节：瞬态动力分析的定义和目的 
    第二节：瞬态分析状态的基本术语和概念 
    第三节：在ANSYS中如何进行瞬态分析 
    第四节：瞬态分析实例

 

第一节  定义和目的  
    什么是瞬态动力分析? 
        -   它是肯定随时间变化载荷（例如爆炸）做用下结构响应的技术； 
        -   输入数据： 
            *   做为时间函数的载荷 
        -   输出数据： 
            *   随时间变化的位移和其它的导出量，如：应力和应变。  
    瞬态动力分析能够应用在如下设计中： 
        -  承受各类冲击载荷的结构，如：汽车中的门和缓冲器、建筑框架以及悬挂系统等； 
        -  承受各类随时间变化载荷的结构，如：桥梁、地面移动装置以及其它机器部件； 
        -  承受撞击和颠簸的家庭和办公设备，如：移动电话、笔记本电脑和真空吸尘器等。

 

第二节  术语和概念  
    包括的主题以下： 
         -   运动方程 
         -   求解方法 
         -   积分时间步长

     (1)  运动方程 
         -    用于瞬态动力分析的运动方程和通用运动方程相同：

                   

                  

  
         -    这是瞬态分析的最通常形式，载荷可为时间的任意函数； 
         -    按照求解方法， ANSYS 容许在瞬态动力分析中包括各类类型的非线性- 大变形、接触、塑性等等。

      (2)  求解方法

                        

                      

  

   
      运动方程的两种求解法： 
            -    模态叠加法（在第六章中讨论） 
            -    直接积分法： 
                  *   运动方程能够直接对时间按步积分。在每一个时间点，需求解一组联立的静态平衡方程（F=ma）； 
                  *   ANSYS  采用 Newmark 法这种隐式时间积分法； 
                  *   ANSYS/LS-DYNA  则采用显式时间积分法； 
                  *   有关显式法和隐式法的讨论请参见第一章。  
            -   求解时便可用减缩结构矩阵，也可用完整结构矩阵； 
            -   减缩矩阵： 
                *   用于快速求解； 
                *   根据主自由度写出 [K]， [C]， [M] 等矩阵，主自由度是彻底自由度的子集； 
                *   减缩的 [K] 是精确的，但减缩的 [C] 和 [M] 是近似的。此外，还有其它的一些缺陷，但不在此讨论。 
            -    完整矩阵： 
                *   不进行缩减。 采用完整的 [K]， [C]， 和 [M] 矩阵； 
                *   在本手册中的所有讨论都是基于此种方法。

    (3)  积分时间步长  
            -   积分时间步长（亦称为 ITS 或 Dt）是时间积分法中的一个重要概念 
                *   ITS = 从一个时间点到另外一个时间点的时间增量 Dt ； 
                *   积分时间步长决定求解的精确度，于是其数值应仔细选取。 
            -   ITS 应足够小以获取下列数据： 
                *   响应频率 
                *   载荷突变 
                *   接触频率（若是存在的话） 
                *   波传播效应（若存在）
  
    (4)  响应频率 
            -  不一样类型载荷会在结构中激发不一样的频率（响应频率）； 
            -  ITS 应足够小以获取所关心的最高响应频率 （最低响应周期）； 
            -  每一个循环中有 20 个时间点应是足够的，即： 
                Dt = 1/20f 
            式中 ，f 是所关心的最高响应频率。

                      

  

    (5)  载荷突变 
        -   ITS  应足够小以获取载荷突变                   

                       

  
    (6)  接触频率 
        -   当两个物体发生接触，间隙或接触表面一般用刚度（间隙刚度）来描述； 
        -   ITS 应足够小以获取间隙“弹簧”频率； 
        -   建议每一个循环三十个点，这才足以获取在两物体间的动量传递，比此更小的 ITS  会形成能量损失，而且冲击可能不是彻底弹性的。

                        

   
    (7)  波传播 
        -   由冲击引发。在细长结构中更为显著（以下落时以一端着地的细棒） 
        -   须要很小的 ITS ，而且在沿波传播的方向须要精细的网格划分 
        -   显式积分法（在 ANSYS-LS/DYNA 采用）可能对此更为适用

 

                        

 

 

第三节  分析步骤  
    在此节中只讨论完整矩阵，五个主要步骤： 
        –  建模 
        –  选择分析类型和选项 
        –  规定边界条件和初始条件 
        –  施加时间历程载荷并求解 
        –  查看结果

     (1)   建模  
    模型 
        -  容许全部各类非线性 
        -  记住要输入密度！ 
        -  其他参见第一章建模所要考虑的问题
    典型命令流：  
        /PREP7 
        ET,... 
        MP,EX,... 
        MP,DENS,… 
        ! 创建几何模型 
        … 
        ! 划分网格 
        ... 
  
    (2)  选择分析类型和选项： 
        -  进入求解器并选择瞬态分析 
        -  求解方法和其它选项- 将在下面讨论 
        -  阻尼 – 将在下面讨论 
  
    典型命令: 
        /SOLU 
        ANTYPE，TRANS，NEW

                       

                        

 

  
    a.  求解方法 
     完整矩阵方法为缺省方法。容许下列非线性选项： 
        –  大变形 
        –  应力硬化 
        –  Newton-Raphson 解法 

    b.  集中质量矩阵 
         主要用于细长梁和薄壁壳或波的传播 

    c.  公式求解器 
         由程序自行选择

                      

 

    典型命令流：  
        TRNOPT，FULL 
        NLGEOM，… 
        SSTIF，… 
        NROPT，… 
        LUMPM，… 
        EQSLV，... 

 d.    阻尼 
        -   α 和 b 阻尼都可用； 
        -   在大多数状况下，忽略 α 阻尼（粘性阻尼），仅规定 b 阻尼（由滞后形成的阻尼）： 
                b = 2x/w 
        式中 x 为阻尼比，w 为主要响应频率 （rad/sec）。
  
    典型命令: 
        ALPHAD，… 
        BETAD，… 

                       

 

     (3)    设置边界条件和初始条件
     -   在这种状况下边界条件为载荷或在整个瞬态过程当中一直为常数的条件，例如： 
            *   固定点（约束） 
            *   对称条件 
            *   重力 
     -   初始条件将在下面讨论

                      

 

    典型命令流：  
        DK，… ! 或 D或 DSYM 
        DL，… 
        DA，… 
        ACEL，… 
        OMEGA，...

     初始条件 
        -   时间 t = 0 时的条件：u0， v0， a0 
        -   它们的缺省值为， u0 = v0 = a0 = 0 
        -   可能要求非零初始条件的实例： 
          *  飞机着陆 (v0 不等于 0) 
          *  高尔夫球棒击球 (v0 不等于 0) 
          *  物体跌落试验 (a0 不等于 0)

 

                     

   
    施加初始条件的两种方法： 
        -   以静载荷步开始 
            *  当只需在模型的一部分上施加初始条件时，例如，用强加的位移将悬臂梁的自由端从平衡位置“拨”开时，这种方法是有用的； 
            *  用于须要施加非零初始加速度时。 
        -  使用 IC  命令 
           *   Solution > Apply > Initial Condit’n > Define + 
           *   当需在整个物体上施加非零初始位移或速度时 IC  命令法是有用的。
  

 

第四节  实例
    实例 -  物体从静止状态下落 
        -  这种状况 a0=g （重力加速度）v0=0 
        -  采用静载荷步法 
        -  载荷步1： 
            *  关闭瞬态效应。用 TIMINT，OFF 命令或：
                Solution > Time/Frequenc > Time Integration... 
            *  采用小的时间间隔，例， 0.001； 
            *  采用 2  个子步， 分步加载（若是采用线性载荷或一个子步， v0  就将是非零的）； 
            *  保持物体静止，例如，固定物体的所有自由度； 
            *  施加等于 g  的加速度； 
            *  求解。

    典型命令流：  
           ! 载荷步 1 
           TIMINT，OFF          !  关闭瞬态效应 
           TIME，0.001           !  小的时间间隔 
           NSEL，…               !  选择全部小物体的全部节点 
           D，ALL，ALL，0    !  并在全部方向上定义固定约束 
           NSEL，ALL 
           ACEL，…              !  加速度值 
           NSUBST，2            ! 两个子步 
           KBC，1                  !  阶梯载荷 
           SOLVE 

  
    -   载荷步 2： 
        *  打开瞬态效应； 
        *  释放物体，例如， 删除物体上的 DOF 自由度约束； 
        *  规定终止时间，连续进行瞬态分析。

                     

 

    典型命令流：  
        ! 载荷步 2 
        TIMINT，ON                     !  打开瞬态效应开关 
        TIME，…                         !  指定载荷步实际的终点时刻 
        NSEL，…                        ! 选择全部小物体的全部节点 
        DDELE，ALL，ALL         !  并删除全部约束 
        NSEL，ALL 
        SOLVE 
        ... 
  
    实例 – 将悬臂梁的自由端从平衡位置“拨”开“ 
        -   这种状况时，在梁的自由端 u0 不等于 0 ，  v0=0； 
        -   用静载荷步法； 
        -   载荷步 1： 
            –  关闭瞬态效应。用 TIMINT，OFF 命令或 
                    Solution > Time/Frequenc > Time Integration... 
            –  采用小的时间间隔，例如， 0.001； 
            –  2 个子步， 分步加载（若是采用线性载荷或用一个子步，v0  就将是非零的）； 
            –  在梁的自由端施加所要求的非零位移； 
            –  求解。

     典型命令流：  
            ! 载荷步 1 
            TIMINT，OFF             !  关闭瞬态效应 
            TIME，0.001              !  小的时间间隔 
            D，…                         !  在指定节点定义强制位移 
            NSUBST，2               !  两个子步 
            KBC，1                     !  阶梯载荷步 
            SOLVE 
  
       - 载荷步2： 
          –  打开瞬态效应； 
          –  删除强加位移； 
          –  指定终止时间，连续进行瞬态分析。
  
      典型命令流：

          ! 载荷步 2 
          TIMINT，ON      ! 打开瞬态效应开关 
          TIME，…        ! 指定载荷步实际的终点时刻 
          DDELE，…       ! 删除全部强制位移 
          ... 
          SOLVE 
          ... 

  
      实例 -  高尔夫球棒端头的初速度 
            -  假定只对高尔夫球棒端头建模，而且整个端头运动，这时有初始条件 v0 不等于0。 同时又假定 u0 = a0 = 0； 
            -  在这种状况下使用 IC  命令法是方便的 
    1  选择球棒上的所有节点； 
    2  用 IC 命令施加初始速度或； 
            –  选择 Solution > Apply > Initial Condit’n > Define + 
             -   选用所有节点 
            –  选择方向并输入速度值 
    3  激活所有节点； 
    4  规定终止时间，施加其它载荷条件（若是存在的话），而后求解。

                       

 

    典型命令流：           NSEL，…          IC，…          NSEL，ALL          TIME，…          …          SOLVE         实例 – 承受冲击载荷的固定平板          -  此种状况下 u0 = v0 = a0 = 0；          -  这些初始条件都是 ANSYS 中的缺省初始条件值，因此这里没必要再规定它们！          -  只施加边界条件和冲击载荷，而后求解。