PDF《管板组合结构等效力学模型》

第3 9卷第8期西南大学学报(自然科学版)

2 0

1 7年8月Vol.

39No.8JournalofS o u t h w e s tU n i v e r s i t y( N a t u r a lS c i e n c eE d i t i o n ) A u g .

2 0

1 7 D O I :

1 0 .

1 3

7 1

8 / j . c n k i . x d z k .

2 0

1 7 .

0 8 .

0 2

6 管板组合结构等效力学模型 ① 王玲, 董有恒, 王成军, 刘岩西南交通大学 力学与工程学院 应用力学与结构安全四川省重点实验室,成都

6 1

0 0

3 1 摘要:通过将管板组合结构等效为正交各向异性板的方法,探讨该类结构的力学响应的计算方法. 首先从平面问题 本构方程出发,应用变形协调条件得到其等效材料参数,然后应用刚度分析法得到其等效厚度与等效密度,最后通 过与3 D 有限元结果比较说明等效模型的正确性. 关键词:管板结构;

等效理论;

正交各向异性板 中图分类号:O

3 1

6 文献标志码:A 文章编号:1

6 7

3 9

8 6

8 (

2 0

1 7 )

0 8

0 1

8 0

0 5 管板组合是换热设备中常见的一种结构形式,广泛应用于冷凝器和锅炉中. 作为锅炉炉膛的焊制水冷 壁,由管板交替排列组成,在工作中承受热流、炉膛烟气压力、自重等多种载荷作用. 对管板组合结构的静动力学分析,前人做了一些研究工作. 文献[ 1-4] 针对焊制水冷壁结构提出了一 些具体的简化方法,并用于热应力分析. 文献[ 5-6 ] 利用二维有限差分和有限元法模拟了水冷壁中温度场 的分布. 文献[

7 ] 建立了锅炉水冷壁系统的有限元仿真模型. 但从现有文献看,有限元法是分析管板组合结 构的重要方法,但其复杂的建模,大的计算量,很难用于结构参数化设计. 由于焊制水冷壁由管板周期排列组合而成,使其整体表现出正交各向异性,文献[ 8-1 1] 提出用各向 异性板对这类结构进行简化,不仅可使计算量减小,还便于参数化设计. 本文将在文献[

1 ] 的基础上,将管板组合结构等效为正交各向异性板,通过变形协调方程和刚度分析 得到其相关的等效参数,通过与3 D 有限元结果比较说明方法的正确性.

1 理论推导

1 .

1 等效材料性质 考虑图1( a ) 所示周期性排列的管板 组合结构,取图1( b ) 所示截距为p的基本单元. 其中,鳍片宽度l,管子外径r 0,管子内径r i,鳍片厚度t f ,管子厚度t,管子中面距管心距离rm . 其等效正交各 向异性板如图2. 平面正交各向异性热力学问题,应力应变关系为 σ x σ y τ x y ì ? í ? ? ? ? ü ? ? ? ? ? ? = Ex x Ex y

0 Ex y Ey y

0 0

0 Gx y é ? ê ê ê ê ù ? ú ú ú ú * ε x ε y γx y é ? ê ê ê ê ù ? ú ú ú ú - α x α y

0 é ? ê ê ê ê ù ? ú ú ú ú T ? è ? ? ? ? ? ÷ ÷ ÷ (

1 ) 其中, σ x , σ y , σ x y 为应力分量, ε x , ε y , γx y 为应变分量, Ex x , Ey y , Ex y , Gx y 为弹性模量和剪切模量, α x = α y = α 为热膨胀系数. ① 收稿日期:2

0 1

5 0

1 2

9 基金项目:国家自然科学基金项目(

1 1

0 7

2 2

0 4 ) . 作者简介:王玲(

1 9

9 1 ) ,女,重庆璧山人,主要从事结构振动与控制的研究. 图1 管板结构模型与横截面模型尺寸图 图2 等效薄板示意图 平面问题的应变协调方程为 ?

2 ε x ? y

2 + ?

2 ε y ? x2 = ?

2 γx y ? x ? y (

2 ) 将应力应变关系( 1) 代入协调方程( 2) ,并利用σ x = ?

2 φ ? y

2 , σ y = ?

2 φ ? x2 , τ x y =- ?

2 φ ? x ? y ,得到用应力函数表示 的变形协调方程 Ey y Ex x Ey y - Ex y ( )

2 ・ ?

4 φ ? y

4 +

1 Gx y -

2 Ex y Ex x Ey y - Ex y ( )

2 ? è ? ? ? ÷・ ?

4 φ ? x2 ? y

2 + Ex x Ex x Ey y - Ex y ( )

2 ・ ?

4 φ ? x4 +α ?

2 T ? x2 + ?

2 T ? y

2 ? è ? ? ? ÷ =

0 (

3 ) 文献[

1 ] 给出了水冷壁结构在面内温度场作用下的变形协调方程为

1 E Fx * ?

4 F ? y

4 +

1 G Fx y -

2 μ E Fy ? è ? ? ? ÷ ?

4 F ? x2 ? y

2 +

1 E Fy * ?

4 F ? x4 +α ?

2 T ? x2 + ?

2 T ? y

2 ? è ? ? ? ÷ =0 (

4 ) 其中, F 为应力函数, T 为温度场, α为材料的热膨胀系数, E 为材料弹性模量, G 为剪切模量, μ为泊松比. Fx , Fy , Fz 为仅与结构几何有关的系数,表示为 Fx = p

1 2 rm t ? è ? ? ? ÷

3 ・ π

4 -

2 π ? è ? ? ? ÷ + π rm

4 t +

2 l t f Fy =

2 l t f +2 π rm t p , Fx y = p t f

2 l+ π rm t f

2 t 比较式(

3 ) 与式(

4 ) ,并利用Ex y =φ Ex x 以及F =φ・Fy 可得 Ex x = Fx Fy -μ

2 Fx ・E Ey y = Fy Fy -μ

2 Fx ・E Ex y =μ Ex x Gx y = Fx y Fy ・G ì ? í ? ? ? ? ? ? ? ? ? ? (

5 )

1 .

2 等效抗弯刚度

1 .

2 .

1 y 方向的抗弯刚度 从管板结构几何构型可知, y 方向的等效抗弯刚度为

2 西南大学学报( 自然科学版) h t t p : / / x b b j b . s w u . e d u . c n 第3 9卷Dy = Dy

1 +Dy

2 p (

6 ) 其中, Dy 1为管部分的抗弯刚度, Dy 2为鳍片部分的抗弯刚度. 由材料力学可得 Dy

1 = π( r

4 0 -r

4 I )

4 ・E Dy y = E h

3 1

2 ( 1-μ

2 ) = E t

3 f l

6 ( 1-μ

2 ) ì ? í ? ? ? ? ? ?

1 .

2 .

2 x 方向的抗弯刚度及扭转刚度 管板结构沿x 方向为变截面,整体抗弯能力主要由截面抗弯刚度最弱处决定. 由薄板弯曲理论得其等 效刚度为 Dx = Ex x h3

1 2 ( 1-μx y μy x ) (

7 ) 式中, μx y 和μy x 为等效后正交各向异性板的泊松比,且μy x Ex x =Ex y , μx y Ey y =Ex y , h 为其等效厚度, 将由 式(

8 ) 给出. 由文献[

8 ] 可得扭转刚度 Dx y 为Dx y = Dx Dy

1 .

3 等效厚度与等效密度 由薄板弯曲理论得y 方向的抗弯刚度为 Dy = Ey y h3

1 2 ( 1-μx y μy x ) 结合公式(

6 ) 及上述 Dy 的表达式得等效厚度为 h=

3 1

2 ( 1-μx y μy x ) Dy Ey y (

8 ) 由等效前后质量相等,得等效后密度ρ

1 为ρ1=ρ2ltf+πr

2 0 -r

2 i ( ) ( ) p h (

9 ) 其中, ρ 为管板材料密度.

2 算例为验证等效结果,分别用三统合有限元方法(

3 D - F EM) 和本文方法进行计算. 计算中取图1( b ) 所示长 为5 0倍节距长宽相等的管板组合结构为计算模型,基本单元参数如表1. 表1 基本单元尺寸及材料参数 几何参数/ mm 节距( p) 内径( r i) 外径( r0) 鳍片宽度( l) 鳍片厚度( t f )

3 4

9 .

4 1

2 .

3 4 .

8 6

37 4 材料参数 杨氏模量( E) / G P a 泊松比( υ) 剪切模型( G) / G P a 热膨胀系数( α) /

1 /℃ 密度( ρ) / ( k g・m-3 )

2 1

0 0 .

2 8

8 0 5*1 0-6

80 0

0 2 .

1 静态挠度 考虑四边简支,受单位均布载荷的正交各向异性板. 其挠度的经典 N a v i e r解[ 8] 为W=∑ ∞ m=1 ∑ ∞ n=1

4 ∫ a

0 ∫ b

0 s i n n π y b s i n n π y b d x d y π2 a b Dx m4 a4 +2 ( μ Dx +2 Dx y ) m2 n2 a2 b

2 +Dy n4 b

4 ? è ? ? ? ÷ s i n m π x a s i n n π y b

3 第8期王玲,等:管板组合结构等效力学模型 其中a, b 为板的长和宽. 计算中取上述级数前2 0项,结果如表2. 表2 板中心处挠度( *1 0-7 m)

3 D - F EM 当前解 相对误差/%

9 .

6 7

0 9 .

1 3

3 5 .

5 5

2 .

2 动态频率 对四边简支正交各向异性板,其频率公式为 fm n = π

2 Dx m a ? è ? ? ? ÷

4 +Dy n b ? è ? ? ? ÷

4 +2μ Dx +2 Dx y ( ) m n a b ? è ? ? ? ÷

2 m 其中, m 表单位面积质量,等于其等效密度与等效厚度之积. 计算管板结构前9阶固有频率,如表3. 表3 管板结构前9阶固有频率 f i j

3 D - F EM / H z 当前解/ H z 相对误差/ % f i j

3 D - F EM / H z 当前解/ H z 相对误差/ % f1

1 2

6 .

8 9

5 2

7 .

6 9

89 2 .

9 8

9 f1

3 1

8 5 .

2 3

6 1

8 6 .

5 4

01 0 .

7 0

4 f2

1 4

3 .

2 4

9 4

6 .

1 6

61 6 .

7 4

5 f2

3 2

0 7 .

3 3

8 2

1 1 .

5 7

54 2 .

0 4

4 f1

2 8

6 .

7 8

9 8

7 .

6 6

87 1

1 .

0 1

2 f3

3 2

4 0 .

7 5

5 2

4 9 .

2 4

03 3 .

5 4

5 f2

2 1

0 7 .

3 9

0 1

1 0 .

7 9

57 3 .

1 7

1 f1

4 3

2 2 .

0 6

6 3

2 4 .

7 6

76 0 .

8 3

9 f3

2 1

3 6 .

7 7

9 1

4 3 .

8 6

07 5 .

1 7

7 由表2及表3可见,使用本文方法得到的结构静力学解误差在6%以内,动态频率误差在7%以内. 说 明本文方法有效.

2 .

3 管板结构尺寸规模对等效结果精度的影响 以上算例是针对规模为5

0 p的结构进行的,实际工程中管板数量达到数百个,现讨论管板规模对等效 模型精度的影响. 图3给出了板中心受单位集中力该处挠度的相对误差随管板规模的变化情况. 图4给出 了四边简支管板前四阶固有频率的相对误差随管板规模的变化情况. 图3 管板中心挠度相对误差 图4 简支管板前四阶固有频率相对误差 由图

3、图4可见,随着管板规模的增加,等效板理论模型与3 D 有限元结果的相对误差逐渐减小,该 简化模型可用于大规模管板结构计算.

3 结论应用经典弹性理论将管板组合结构等效为正交各向异性板,得到其等效材料参量;

使用3 D 有限元法 分别从静力学和动态特性两方面对等效模型进行了验证;

讨论了等效模型对管板规模的适用性. 研究表明 该等效模型可应用于大规模管板组合结构的静动力学计算.

4 西南大学学报( 自然科学版) h t t p : / / x b b j b . s w u . e d u . c n 第3 9卷 参考文献: [

1 ] S E T OGU CH IK,WA D A H,M I YA Z O E M. S t r u c t u r a lA n a l y s i so fW e l d e d W a l l sf o rB o i l e rF u r n a c e s [ J ] . M i t s u b i s h i T e c hR e v ,

1 9

7 4,

1 1 (

2 ) :

1 3 7-1

4 5. [

2 ] S TUMP P W. T h e r m a lS t r e s s e s i nS t e a mB o i l e r sw i t h W e l d e dT u b eW a l l s[ R] . S u l z e rT e c h n i c a lR e v i e w,

1 9

6 6. [

3 ] S A I HO Y, T E R A D A K. T h eS t r e n g t ho fW e l d e d W a l lP a n e l s f o rB o i l e rF u r n a c e s[ J ] . T e c h n i c a lR e v i e w,

1 9

6 6, 3( 2) :

1 6-2 5. [

4 ] K I M T W,K I MJC,HWANGS H. S t r u c t u r a lA n a l y s i so fU t i l i t yB o i l e rW a t e r w a l l[ C] / / A S ME2

0 0 5I n t e r n a t i o n a l M e c h a n i c a lE n g i n e e r i n gC o n g r e s sa n dE x p o s i t i o n . Am 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 ,

2 0

0 5:

6 9 7-7

0 2. [

5 ] 叶科,吕俊复,张建胜,等. 膜式水冷壁换热系数的有限元分析 [ J ] . 清华大学学报( 自然科学版) ,

2 0

0 6,

4 6(

1 1) :

1 8

8 8-1

8 9 1. [

6 ] 李志宏,刘文铁,刘石. 膜式水冷壁壁温影响因素的数值分析 [ J ] . 热能动力工程,

2 0

0 3,

1 8 (

2 ) :

1 7 3-1

7 7. [

7 ] 黄永和,杨晨,何航行. 超临界直流锅炉水冷壁系统建模与仿真 [ J ] . 计算机仿真,

2 0

1 3,

3 0 (

1 ) :

3 1 4-3

1 9. [

8 ] 徐芝纶. 弹性力学 [ M] . 北京:高等教育出版社,

2 0

0 6. [

9 ] 周学军,张之峰. 压型钢板等效为正交各向异性板的有限元分析 [ J ] . 钢结构,

2 0

1 0,

2 5 (

8 ) :

3 5-3 8. [

1 0 ]苏文政,刘书田,张永存. 桁架板等效刚度分析 [ J ] . 计算力学学报,

2 0

0 7,

2 4 (

6 ) :

7 6 3-7

6 7. [

1 1 ]胡肇滋,钱寅泉. 正交构造异性板刚度计算的探讨 [ J ] . 土木工程学报,

1 9

8 7,

2 0 (

4 ) :

4 9-6 1. A nE q u i v a l e n tM e c h a n i c a lM o d e l f o ra T u b e - P l a t eC o m p o s i t eS t r u c t u r e WANG L i n g , D ONGY o u - h e n g , WANGC h e n g - j u n , L I U Y a n A p p l i e dM e c h a n i c sa n dS t r u c t u r eS a f e t yK e yL a b o r a t o r yo fS i c h u a nP r o v i n c e, S c h o o l o fM e c h a n i c sa n dE n g i n e e r i n g, S o u t h w e s t J i a o t o n gU n i v e r s i t y, C h e n g d u6

1 0

0 3 1, C h i n a A b s t r a c t : T h ec o m p u t a t i o n a lm e t h o do f t h em e c h a n i c a l r e s p o n s e so f t h et u b e - p l a t ec o m p o s i t es t r u c t u r e i s g i v e nb yp u t t i n gt h ec o m p o s i t es t r u c t u r et oa ne q u i v a l e n to r t h o t r o p i ca n i s o t r o p i cp l a t e . F i r s t , t h ee q u i v a - l e n tm a t e r i a l p a r a m e t e r sa r eo b t a i n e db yc o m b i n i n gt h ec o n s t i t u t i v ee q u a t i o n so ft h ep l a n ep r o b l e m sa n d t h ed e f o r m a t i o nc o m p a t i b i l i t yc o n d i t i o n s . T h e nt h ee q u i v a l e........