冲击载荷下柔性储液罐动态响应数值模拟及规律分析.pdf

《冲击载荷下柔性储液罐动态响应数值模拟及规律分析.pdf》由会员分享，可在线阅读，更多相关《冲击载荷下柔性储液罐动态响应数值模拟及规律分析.pdf（7页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第3 1 卷第5 期2 0 1 1 年9 月爆炸与冲击E X P L O S I O NA N DS H O C KW A V E SV 0 1 3 1，N O 5S e p t，2 0 1 1文章编号：1 0 0 1 1 4 5 5(2 0 1 1)0 5 0 4 6 9 0 6冲击载荷下柔性储液罐动态响应数值模拟及规律分析曹源，金先龙，李政(上海交通大学机械系统与振动国家重点实验室，上海2 0 0 2 4 0)摘要：针对柔性材料大变形、非线性等特点，对柔性储液罐的冲击响应规律进行研究，依据实验结果，采用A L E 有限元方法对冲击情况下柔性罐状储液容器及其中流体的动态响应进行了数值模拟。所

2、得结果与实验结果具有较好的一致性。在此基础上，深入分析了空投高度、装水量等对柔性储液罐动态响应的影响规律。重点关注了表面一些关键部位的变形和应力变化情况。结果表明，柔性储液罐的动态应力和变形均与抛落高度成正比。但表面最大应力与装水量呈“V”型关系。关键词：流体力学；动态响应；A L E 方法；柔性储液罐；流固耦合；冲击载荷中图分类号：0 3 5 3 1国标学科代码：1 3 0 2 0 7 0文献标志码：A储液容器广泛应用于储运、航空等领域，随着材料性能的提高，在许多抗冲击场合，如空投、爆炸中，柔性储液容器开始被越来越多的应用。许多学者已对容器内液面晃动、容器冲击响应以及流固耦合等问题进行了长期

3、深入的研究 1。3 。与早期研究相比，柔性容器冲击响应问题的难点是大变形和非线性，使得其在冲击过程中的动态响应规律更加复杂。对于这类问题的研究，目前还是以实验为主。随着计算方法的改进和计算机性能的提高，数值模拟手段已有可能发挥更大的作用。对于大变形的流固耦合问题，很难用单一的L a g r a n g e 方法或E u l e r 方法进行研究。所以A L E 方法开始应用于流固耦合分析H 5 。Z H A N GA i n i a n 等口3 分别使用附加质量法、L a g r a n g e 方法和A L E方法分析了船舶碰撞时流固耦合作用，并认为与A L E 方法相比，附加质量方法和L

4、a g r a n g e 方法低估了液体与结构之间相互作用的影响。M A n g h i l e r i 等口 研究了储液箱体与地面碰撞时，箱内液面大幅晃动及箱体变形的问题，并将A L E 方法与其他数值模拟方法进行了比较，并且与实验结果吻合很好。本文中针对跌落或空投中的柔性储液罐动态响应进行实验研究和数值模拟。采用基于L a g r a n g e方法描述下的B e l y t s c h k o T s a y(B T)膜单元建立柔性储液容器模型，流体部分则是基于A L E 方法建模，并利用罚函数法实现流固耦合。在实验基础上，使用数值模拟的方法深入分析了空投高度、装水量等因素对容器动态响

5、应的影响规律。1 流固耦合方法与原理利用A L E 方法建模与罚函数法相结合对冲击过程中的流固耦合现象进行数值模拟。罚函数法的优点在于保证了耦合接触过程中的能量守恒。首先，对于每个结构节点，搜索包含该节点的A L E 单元。然后，通过计算穿透速度与时间来计算流体节点的穿透深度。最后，通过正比穿透深度计算基于罚函数的耦合力。耦合力可以看作是节点总力中的一个外部力。在此基础上，通过计算每一个时间步长的总节点力，可以得到流体和结构体耦合界面上每一点的速度、位移和应力等物理量，从而实现流固耦合界面上的相互作用，即耦合计算。A L E 方法则综合了L a g r a n g e 方法与E u l e r

6、 方法的优点，在材料域*收稿日期：2 0 1 0 0 7 1 2；修回日期：2 0 1 1-0 1 1 7基金项目：国家自然科学基金项目(1 1 0 7 2 1 5 0)；教育部博士点基金项目(0 0 7 0 2 4 8 1 1 0)；机械系统与振动国家重点实验室基金项目(M S V-M S2 0 1 1 0 7)作者简介：曹源(1 9 7 6 一)，男，博士，讲师。万方数据4 7 0爆炸与冲击第3 1 卷与空间域外引入了参考域，并在参考域网格上求解，既解决了L a g r a n g e 方法中材料严重扭曲的问题，又解决了E u l e r 方法中移动边界引起的复杂性问题。基于A L E 方

7、法的连续性方程和动量方程为弓；I，+c i。；A f-p。3。7 )i j=o(1)I D 等l，+阳瓦O v f _ 惦3 a o+6 i(2)式中：?(为A L E 坐标；u i 为水的流动速度；z i 和z J 为空间坐标；C r 和c，均为对流速度；P 为流体密度；b：为流体体力，为应力张量。对储液容器变形的描述使用弹性体的连续方程i(3)8 万Ix 一茜十i式中：X 是L a g r a n g e 坐标，p。为容器密度，为体力，“为固体结构位移。对于水这种弱可压缩流体，引入线性牛顿流体本构方程一一p 8#+2 口(4)式中：P 为水的静态压力，t l=-(s 4)为动力粘性系数，s

8、#为应变率张量，岛为克罗内克常数。流固耦合需在液体和结构界面上满足几何相容条件和力的平衡条件，即一口。一譬f，F+F。一0Vit(5)5 口s il，r 十s2L o J“6X式中m 和u。分别为流固耦合界面上水的流动速度和固体材料速度，F。和F。分别是流体和固体结构作用在流固耦合界面上的力，这些力可由罚函数方法计算获得F。一万d Z Z+面d Z+2 z(6)式中：Z 是穿透量；车是阻尼系数，当2 时，振荡趋于收敛，因此可用于消除计算中的振荡。定义叫一k 些，m s 和m f 分别为固体和流体质量；k 一掣，K i 为单元容积模量，A i 代表单元面积，V。、优。一m fV i是单元体积，口

9、为缩放系数。根据上述计算式可以得到阻尼系数e，取一o 5。2实验简介为验证本文中方法和模型的正确性，首先对装水量刀一8 0 的柔性储液罐进行1 0 0m 高度的空投冲击实验。所使用的储液罐为矽o 5I T I 0 8m 的圆柱体，如图1 所示。储液罐壁厚为1 0m m。罐体采用帘布材料，该帘布采用尼龙织物，基体采用天然橡胶。其顶端由钢盖密封。在1 0 0m 高度，以零初始速度垂直抛落，整个抛落过程由高速照相机记录。实验中，罐体肩部(位置1)破裂，储液罐内的水溅出。所以，罐体肩部被确定为脆弱位置之一。此外，罐体中部(位置2)和底部(位置3)也是易受冲击和易损部位，这些部位也是数值模拟时关注的重点

10、。3 数值模拟与分析厂、A i rW a t e r叉C o n t a i n e rb o d y(P o s i t i o n2、图1 储液罐结构F i g 1S t u c t u r eo fc o n t a i n e r3 1 数值模型以实验为基础，建立三维有限元数值模型。该模型包括：地面、罐状柔性储液容器、水和容器内空万方数据第5 期曹源等：冲击载荷下柔性储液罐动态响应数值模拟及规律分析4 7 1气。相对于柔性容器，冲击过程中混凝土地面变形较小，模型中将地面材料作为弹塑性体处理。容器的钢盖处理为弹性体。各项参数如表1 所示，其中，l D 为密度，E 为弹性模量，v 为泊松比

11、，盯，为屈服应力，G为剪切模量，p 为硬化系数。容器主体复合材料结构的力学性能是数值模拟的关键，本文中采用X I AZ h i h u i 等嘲提出的周期性边界加载方法计算获得帘布材料的力学性能参数，如表2表1 材料参数T a b l e1M a t e r i a lp o r p e r t i e sa n dp a r a m e t e r s所示。其中E，、E z 和E 3 分别为材料3 个方向的弹性模量；G，。、G。和G。分别为材料3 个方向的剪切模量m z、地和蛐分别为材料3 个方向的泊松比。由表2 可知材料呈横观各向异性。由于容器壁较薄，所以采用基于L a g r a n g

12、 e 方法描述的B T 膜单元，沿膜厚度方向取3 个积分点。表2 容器材料力学参数T a b l e2M e c h a n i c a lp r o p e r t i e so ft h ec o m p o s i t e为了节省计算时间，忽略罐体坠落过程中的空气阻力，将抛落高度换算为罐体冲击地面时的速度，作为初始条件。对于内部流体采用基于A L E 方法描述的六面体单元。为了确保结果的准确性，所有单元尺寸和形状均仔细控制。整个模型单元数为2 2 80 3 3，节点数为2 3 84 8 4。数值模拟采用基于显示计算方法的A N S Y S L S-D Y N A 有限元软件，为确保计算的

13、稳定性，时间步长由模型的最小单元尺寸控制，为约2 0 弘s。在流固耦合的模拟中，过多的耦合点将导致数值计算不稳定，而耦合点较少则易产生泄漏现象。由于柔性容器为曲面结构，存在较小的夹角，故耦合点数取为4。此外，耦合最小体积参数决定了流体物质达到某个单元体积时，流体与结构发生耦合作用。由于空投物体在触地时属于高速碰撞冲击问题，极易产生泄漏现象，所以计算中将耦合最小体积参数值设为0 1，可以有效地减小流体的泄漏现象。3 2 数值模拟与实验结果比较对相同条件下的数值模拟和实验结果进行比较，以验证方法和模型的有效性。定义罐体直径变化量A d 来评估柔性罐体的变形情况。所谓直径变化量，是指冲击过程中不同时

14、刻罐体中间直径d 与盛有相应水量的罐体静置于地面时罐体中间直径d。之差f OA d=【d 一巩对罐体径向变形量的实验结果和柔性罐状容器的数值模拟结果相比较发现，数值模拟结果相对于物理实验结果要小一些，但仍处于可接受的范围，如图2 所示。不同时刻下柔性容器应力分布和内部流体的变化，如图3 所示。从图中可以看出，流体形态与柔性容器变形一致，其内部的空气先是被压缩，随着柔性容器的反弹，空气完全被水包围。当冲击开始时，柔性容器表面最大应力出现在最先与地面接触的底部；反弹时，由于钢盖惯性较大，在顶部形成“弹坑”形状，此时，容器肩部变形最大，应力也最大。t=0t 0图2 变形量比较F i g 2T h e

15、c o m p a r i s o no ft h ed e f o r m a t i o n(8)万方数据口争。蔫第3 1 巷图3 柔性容器表面应力分布及漉体形态F i z 3S tr e s sd i s t r i b u t i o no ft h ec o n t a i n e ra n dc o n f i g u r a t i o no ft h ef l u i d柔性容器上3 个位置主应力如变化曲线如图4 所示。从曲线中可以看出3 个位置的应力依次达到最大值，与这些位置受冲击的先后顺序一致。其中，肩部的最大应力几乎是其他2 个位置的2 倍。这也与实验中罐体反弹过程中，肩

16、部破裂的结果十分吻合。3 3 响应规律分析及讨论从抛落高度和装水量等2 个方面对容器的冲击响应进行研究。对装水量日；8 0 时，1 0 0、2 0 0 和4 0 0m 等3 个抛落高度进行数值模拟。图5 显示了这3 种高度下，3 个位置上最大主应力J。的比较。图4 关键位置处主应力曲线F i glI r i n C i l l,q】s l r c s s I【k t-vp o s i l i o n s图5 不同高度下最大应力F i g5P r i n c i p a ls r c s s e sa Id i f I e r e n th e i g h t s万方数据第5 期曹源等：冲击载荷

17、下柔性储液罐动态响应数值模拟及规律分析4 7 3柔性容器肩部的应力一直大于另外2 个位置的，而随着高度的增加，容器中部的应力增加最快。当高度达到4 0 0m 时，该位置的最大应力几乎与肩部相同。所以，随着高度的增加，容器中部也成为危险部位。另一方面，容器底部的应力随高度的增加，升高缓慢。特别是当高度低于2 0 0m 时，容器底部应力变化不大。图6 给出了3 个抛落高度下柔性容器直径的变化曲线。容器直径变化量和最大幅度都随高度的增加而变大。其中容器直径变化量的最大值从3 0 6m m 增加到6 4m m。图中还可以看出，尽管抛落高度不同，但振荡收敛的时间没有大的改变。图7 显示的是空抛落高度1

18、0 0m 时，装水量刁一3 0，5 0，8 0 和1 0 0 等4 种情况下容器上3 个位置的最大应力。在装水量不同的情况下仍然是容器肩部应力最大，且与装水量呈“V”型关系。当叩=8 0 时，肩部应力是4 种情况下最小的。原因可能是：当装水量少时，钢盖可以上下大幅运动，导致容器肩部产生大变形。随着装水量的增加，肩部变形减小，应力也降低。但是当容器1 0 0 充满水时，冲击产生的力和能量可以通过水直接作用在肩部。由于水的可压缩性很小，刚度又大于柔性材料，减震效果小于容器本身，所以导致肩部的应力急剧增加。与肩部应力的变化不同，容器底部和容器中部的应力与装水量的变化的关系不大。还发现容器的变形量随着

19、装水量增加而呈线性增加。图6 不同高度下柔性容器直径变化F i g 6C o n t a i n e rd e f o r m a t i o n sa td i f f e r e n th e i g h t s圆S h o u1 d e,团C o n tU n er b o d y囚B o t t(I m蘑雾耄鋈薹鋈亵薹莲4 06 08 01 0 0n，图7 不同装水量下最大应力F i g 7P r i n c i p a ls t r e s s e su n d e rd i f f e r e n ta m o u n to fw a t e r4结论在对柔性储液罐进行空投实验的基

20、础上，采用基于A L E 方法描述的流固耦合状态和三维有限元模型对柔性储液容器的冲击响应进行研究。实验和数值模拟结果都显示，柔性储液容器的肩部是其薄弱部位。当达到一定的抛落高度时(如4 0 0m)，容器中部也会变成另外一个薄弱部位。容器的装水量对其抗冲击性有很大影响。其中肩部的应力与装水量呈“V”型关系。而柔性容器的变形则基本与装水量呈线性关系。在其他条件不变的情况下，装水量为8 0 时，柔性储液容器的抗冲击性最好。参考文献：1 3M o i s e e vNN I n t r o d u c t i o nt ot h et h e o r yo fo s c i l l a t i o n

21、 so fl i q u i dc o n t a i n i n gb o d i e s J A d v a n c e si nA p p l i e dM e c h a n i c s，1 9 6 4，8(5)：2 3 3 2 8 9 2 M aDC，S uTC F l u i d-s t r u c t u r eV i b r a t i o na n dL i q u i dS l o s h i n g C f fP r o c e e d i n g so fP r e s s u r eV e s s e l sa n dP i p i n gC o n f e r e

22、n c e N e wY o r k：A m e r i c a nS o c i e t yo fM e c h a n i c a lE n g i n e e r s，19 8 7 3 3S h i m i z uN A d v a n c e sa n dt r e n d si ns e i s m i cd e s i g no fc y l i n d r i c a ll i q u i ds t o r a g et a n k s J J a p a nS o c i e t yo fM e c h a n i c a lE n g i n e e r sI n t e r

23、 n a t i o n a lJ o u r n a l，1 9 9 0，3 3(2)：1 1 1-1 2 4 4 3N i t i k i t p a i b o o nC，B a t h eKJ A na r b i t r a r yL a g r a n g i a n-E u l e r i a nv e l o c i t yp o t e n t i a lf o r m u l a t i o nf o rf l u i d s t r u c t u r ei n t e r a e t i o n J C o m p u t e r s&S t r u c t u r e

24、 s，1 9 9 3，4 7(4 5)：8 7 1-8 9 1 万方数据4 7 4爆炸与冲击第3 1 卷 5 3 6 7 8 N o m u r aT A L Ef i n i t ee l e m e n tc o m p u t a t i o n so ff l u i d-s t r u c t u r ei n t e r a c t i o np r o b l e m s -J C o m p u t e rM e t h o d si nA p p l i e dM e c h a n i c sa n dE n g i n n e e r i n g，1 9 9 4，1 1

25、2(1 2 3 4)：2 9 1 3 0 8 Z H A N GA i-n i a n，S u z u k iK Ac o m p a r a t i v es t u d yo fn u m e r i c a ls i m u l a t i o n sf o rf l u i d-s t r u c t u r ei n t e r a c t i o no fl i q u i d f i l l e dt a n kd u r i n gs h i pc o l l i s i o nI-J O c e a nE n g i n e e r i n g，2 0 0 7，3 4(5 6

26、)：6 4 5 6 5 2 A n g h i l e r iM，L u i g iML，C a s t e l l e t t i，e ta 1 F l u i d s t r u c t u r ei n t e r a c t i o no fw a t e rf i l l e dt a n k sd u r i n gt h ei m p a c tw i t ht h eg r o u n dE J 3 I n t e r n a t i o n a lJ o u r n a lo fI m p a c tE n g i n e e r i n g，2 0 0 5，3 1(3)：2

27、 3 5 2 5 4 X I AZ h i h u i，Z H A N GY u n-f a n g，E l l y i nF Au n i f i e dp e r i o d i c a lb o u n d a r yc o n d i t i o n sf o rr e p r e s e n t a t i v ev o l u m ee l e r n e n t so fc o m p o s i t e sa n da p p l i c a t i o n s l,J I n t e r n a t i o n a lJ o u r n a lo fS o l i d sa

28、n dS t r u c t u r e，2 0 0 3，4 0(8)：1 9 0 7 1 9 2 1 D y n a m i ca n a l y s i so ff l e x i b l ec o n t a i n e r su n d e ri m p a c t。C A OY u a n，JI NX i a n-l o n g，L IZ h e n g(S t a t eK e yL a b o r a t o r yM e c h a n i c a lS y s t e ma n dV i b r a t i o n，S h a n g h a iJi a o t o n gU

29、 n i v e r s i t y，S h a n g h a i2 0 0 0 4 0，C h i n a)A b s t r a c t：T ou n d e r s t a n da n di m p r o v et h ec r a s h w o r t h i n e s so ff l e x i b l ec o n t a i n e r s，t h ed y n a m i cb e h a v i o r so ft h ec o n t a i n e ru n d e ri m p a c tl o a d i n gw e r en u m e r i c a

30、l l ys i m u l a t e d T h ew a t e rm o v e m e n ta n ds l o s h i n gi nt h ec o n t a i n e rw e r em o d e l e db yu s i n gt h eA L Em e t h o d T h ei n t e r a c t i o nb e t w e e nt h ec o n t a i n e ra n df l u i dw a ss t u d i e db ya p p l y i n gt h ep e n a l t ym e t h o d T h em o

31、 d e la n da p p r o a c h e sw e r ev a l i d a t e db yc o m p a r i n gn u m e r i c a lr e s u l t sw i t he x p e r i m e n t a ld a t a B a s e do nt h ea b o v ei n v e s t i g a t i o n，n u m e r i c a ls i m u l a t i o nw a sp e r f o r m e dt oe x p l o r et h ei n f l u e n c e so fd r o

32、 ph e i g h ta n dw a t e rq u a n t i t yo nt h ed y n a m i cr e s p o n s eo ft h ec o n t a i n e r A n dt h ed e f o r m a t i o na n ds t r e s sd i s t r i b u t i o na tt h ek e yp o s i t i o n sw e r ea n a l y z e d R e s u l t sd i s p l a yt h a tt h es h o u l d e ro ft h ec o n t a i

33、n e ri st h em o s tv u l n e r a b l ep a r td u r i n gi m p a c t；t h a tb o t hd r o ph e i g h ta n dw a t e rq u a n t i t yh a v eg r e a te f f e c t so nt h ec r a s h w o r t h i n e s so ft h ec o n t a i n e r；t h a tt h em a x i m u md y n a m i cs t r e s so ft h ec o n t a i n e ri sp

34、 r o p o r t i o n a lt od r o ph e i g h ta n dt h e r el i e saV s h a p e dr e l a t i o nb e t w e e nt h em a x i m u md y n a m i cs t r e s sa n dw a t e rq u a n t i t y K e yw o r d s：f l u i dm e c h a n i c s；d y n a m i cr e s p o n s e；A L E；f l e x i b l ec o n t a i n e r；f l u i d s

35、t r u c t u r ei n t e r a c t i o n；i m p a c t1 0 a dR e c e i v e d1 2J u l y2 0 1 0；R e v i s e d1 7J a n u a r y2 0 1 1S u p p o r t e db yt h eN a t i o n a lN a t u r a lS c i e n c eF o u n d a t i o no fC h i n a(1 1 0 7 2 1 5 0)C o r r e s p o n d i n ga u t h o r：C A OY u a n，c a o y u a

36、n s j t u e d u c n(责任编辑王易难)万方数据冲击载荷下柔性储液罐动态响应数值模拟及规律分析冲击载荷下柔性储液罐动态响应数值模拟及规律分析作者：曹源，金先龙，李政，CAO Yuan，JIN Xian-long，LI Zheng作者单位：上海交通大学机械系统与振动国家重点实验室,上海,200240刊名：爆炸与冲击英文刊名：Explosion and Shock Waves年，卷(期)：2011,31(5)参考文献(8条)参考文献(8条)1.Moiseev N N Introduction to the theory of oscillations of liquid conta

37、ining bodies 1964(05)2.Ma D C;Su T C Fluid-structure Vibration and Liquid Sloshing 19873.Shimizu N Advances and trends in seismic design of cylindrical liquid storage tanks 1990(02)4.Nitikitpaiboon C;Bathe K J An arbitrary Lagrangian-Eulerian velocity potential formulation for fluid-structure intera

38、ction 1993(4/5)5.Nomura T ALE finite element computations of fluid-structure interaction problems 1994(1/2/3/4)6.ZHANG Ai-nian;Suzuki K A comparative study of numerical simulations for fluid-structure interaction ofliquidfilled tank during ship collision 2007(5/6)7.Anghileri M;Luigi M L;Castelletti Fluid-structure interaction of water filled tanks during the impact withthe ground外文期刊 2005(03)8.XIA Zhi-hui;ZHANG Yun-fang;Ellyin F A unified periodical boundary conditions for representative volumeelements of composites and applications外文期刊 2003(08)本文链接：http:/