Page 144 - 《振动工程学报》2025年第9期
P. 144

2074                               振     动     工     程     学     报                     第 38 卷

                        (
                     w  c M ∂M   N ∂N     T  ∂T  )  dx          式中,A x 为等效单元的面积;l z 代表           z 方向的长度。
                                                                       y
                 δ =           +       +               (4)
                  F z  0 EI ∂F  EA ∂F    κGA ∂F cosθ
                                                                                                           x
                              z       z        z                其他等效剪切和拉伸弹性模量的推导过程与                       G z 的
                       (                         )
                     w
                      c M ∂M     N ∂N     T   ∂T   dx
                δ  =           +       +               (5)      推导过程类似,这里不再赘述。此外,等效模型的泊
                 M
                  0   0 EI ∂M   EA ∂M    κGA ∂M cosθ
                             0        0         0               松比   μ y 和材料密度     ρ e 的分析推导过程,可参考文
                                                                                    q
                                                                      x
              式中,EI、EA    和  κGA  分别为抗弯、抗拉和抗剪刚度;
                                                                献  [12, 18]。
              κ 为剪切修正系数,本文取           κ=5/6。
                  将式  (2) 代入式   (3)~(5),由于结构应处于纯剪状
                                                                2    波  纹  夹  芯  板  有  限  元  等  效  模  型
              态,最高点的转角        δ M 0 和竖直位移   δ F z 需为  0。经简单
              推导,可得:

                                   a   a   a                  2.1    三阶剪切变形理论

                                    11  12  13
                                   a 21  a 22  a 23

                                                                  为了构造波纹夹芯板等效有限元模型,首先需
                                1    a 31  a 32  a 33
                          δ =                        (6)
                               EA    a   a                    确定夹芯板结构的截面横向剪切变形。为此,国内
                           F x
                                      22  23
                                     a   a                    外学者已经提出了多种变形理论。对于工程结构中

                                       32  33
              其中:
                  w  [  (         ) 2                ]          大量使用的中厚板或厚板,MINDLIN               提出了考虑剪
                    c 12        π x        2(1+µ)      dx
                                                    2
                                        2
               a =       h+hcos     + cos θ+      sin θ    ,                            [19]                 [20]
                11
                    0  t  2      c            κ        cosθ     切变形的一阶剪切变形理论             。在此基础上,REDDY
                          [
                       w  c 12  (    π x  )                     进一步提出了三阶剪切变形理论,即假设变形前垂
               a = a =        h+hcos     (c− x)−
                12
                    21
                         0  t 2       c                         直于平板中面的直线,变形后为三阶曲线。与实际
                                           ]
                             2(1+µ)          dx
                    sinθcosθ +      sinθcosθ    ,               变形情况相比,三阶剪切变形理论能更好地反映中
                                κ           cosθ
                                                                厚板、复合材料层合板           [21]  等截面的剪切变形和剪应
                       w      (         )
                         c 12        π x  dx
               a = a =    −   h+hcos         ,                  力沿厚度方向的分布情况。本文将采用三阶剪切变
               13
                    31
                        0  t 2        c  cosθ
                     [
                   w  c 12            2(1+µ)     ]  dx          形理论描述波纹夹芯板的变形。
                             2
                                  2
                                               2
               a =      (x−c) +sin θ +      cos θ     ,
               22      2
                    0  t                 κ        cosθ              根据三阶横向剪切变形理论              [20] ,板内任意点的
                       w
                         c 12     dx
               a = a =      (x−c)    ,                          位移为:
               23   32  0 t 2    cosθ                                                            (       )
                   w  c 12 dx                                                                4Z  3   ∂w 0  ,
               a =                                     (7)        u(X,Y,Z) = u (X,Y)+Zθ (X,Y)−  2  θ +
                                                                            0
                                                                                                  Y
                                                                                      Y
               33     2                                                                      3H       ∂X
                    0 t cosθ
                                                                                                 (        )
                  在一个胞元内,单位载荷            F x 实际是由左右两部                                        4Z  3    ∂w 0
                                                                  v(X,Y,Z) = v (X,Y)−Zθ (X,Y)−   −θ +      ,
                                                                                                    X
                                                                            0
                                                                                     X
              分的半波纹芯共同承担的,如图                3  所示。由于这两                                       3H  2     ∂Y
              部分是对称的,因此在自由边的真实水平位移                      δ x 仅     w(X,Y,Z) = w (X,Y)                     (10)
                                                                             0
              为半波纹位移       δ F x 的一半,即:                         式中,θ X 为板中面法线绕          X  轴的转角;θ Y 为板中面
                                     1
                                                       (8)      法线绕    Y  轴的转角;H    为整个板的厚度。若不考虑
                                 δ x = δ F x
                                     2
                                                                u  和  v 表达式中的第三项,则三阶横向剪切变形退化
                                    2c
                                                                为一阶剪切变形模型。
                                      δ x         δ x
                                  F x
                                                                    对于正交各向异性板,根据材料的物理方程,板
                                                                内各点的应力-应变关系为:
                       2h
                                                                        (i)               (i)   
                                                                      σ      Q   Q     0   ε 
                                                                       X      11  12         X 
                                                                                              
                                                                     
                                                                                              
                                                                         
                                                                                                   
                                                                              
                       z                                              σ  =  Q    Q     0   ε       (11)
                                                                                             
                                                                              
                                                                       Y       21  22       Y 
                                                                                              
                                                                       τ XY    0   0   Q 66   γ  XY  
                                                                                              
                      o      x
                                                                         τ YZ     Q 44  0     γ YZ
                                                                       {   } (i)  [       ] (i) {  }
                  图 3 一个完整波纹芯胞元受         F x 作用的变形状况                     τ    =    0  Q       γ          (12)
                                                                          XZ            55     XZ
              Fig. 3 The  deformation  of  a  complete  corrugated  core  under
                                                                式中,i = 1, 2, 3  分别代表波纹夹芯板的上面板、波纹
                    the force F x
                                                                芯层等效板和下面板。各弹性系数分别为:

                  根据剪切模量的定义,波纹芯的等效剪切模量                                     E                  E µ
                                                                     (i)
                                                                                         (i)
                                                                                    (i)
                                                                   Q =       X   ,Q = Q =       Y  XY  ,
              G z 为:                                                 11  1−µ µ YX   12   21  1−µ µ YX
                x
                                                                                                 XY
                                                                            XY
                                   F                                 (i)   E Y      (i)       (i)
                                    x    F  x                      Q =           ,Q = G ,Q = G ,
                                                                                              55
                                                                                                   XZ
                                                                     22
                                                                                    44
                                                                                         YZ
                              τ    A xy  2cb   2h                       1−µ µ  YX
                                                                            XY
                         G =   zx  =  =     =          (9)           (i)
                          zx
                              γ zx  δ x  δ x  cδ F x               Q = G  XY                             (13)
                                                                     66
                                    l    2h                         以上本构方程的矩阵形式为:
                                    z
   139   140   141   142   143   144   145   146   147   148   149