您的当前位置:首页正文

薄壁杆件的有限元分析法

2024-10-18 来源:威能网
维普资讯 http://www.cqvip.com 科技慵报开发与经济 文章编号:1005—6033(2006)13—0140—04 SC1-TECH INFORMATION DEVELOPMENT&ECONOMY 2006年第16卷第13期 收稿日期:2006—02—28 国固 围曰圆 国 髓 牛庆芳 ,尹永青 (1.太原理工大学,山西太原,030024;2.山西省建筑设计研究院,山西太原,030013) 摘要:阐述了薄壁杆件理论的发展及经典的薄壁约束扭转理论的不足,对应用十分 广泛且基于势能驻值原理的一种数值解法——有限元理论进行了进一步的研究,提出 薄壁空间杆件单元中每个节点应增加一项表征截面翘曲的位移分量,即截面的翘曲角 ,以计及截面翘曲的影响。指出有限元法完善了薄壁杆件理论,对实际工程的结构设 计有一定的参考价值。 关键词:薄壁杆件;有限元分析;截面翘曲;位移分量 中图分类号:TU32 文献标识码:A 薄壁结构在当今桥梁与建筑工程中是一种比较普遍的结构形式,随 个待定常数m和oa: 着我国国民经济持续、快速发展,交通运输和城市基础设施的建设日益 u(x):口l (1) 受到重视。公路、铁路、城市道路和立交枢纽等各种工程建设正蓬勃发 展,方兴未艾。与此相应,在桥梁结构和各种建筑工程中薄壁结构的应用 ———_厂——一二二二 r——一也日益广泛。 X ’  通常把薄板、壳体和薄壁截面组成的结构称为薄壁结构。与其他结 图1拉压单元 构形式相比,在基于满足强度要求下,薄壁结构具有重量轻、强度大、能 待定常数与位移参数的关系,可根据端点位移条件来确定。端点位 充分利用材料的特点。 移条件为: 本文将研究一种特殊的薄壁结构——薄壁杆件结构,薄壁结构是一 当x=0时,U--Ul,al 种特殊的壳体。它具有如下的几何特征:横截面最大几何尺寸与梁的长 当x=l时,u=u2, +啦Z=I上2, ---Ul 度相比为小量级,杆件的壁厚与横截面最大几何尺寸相比为小量级。由 于薄壁结构同时具有壳体和粱的结构优点,从而使得它成为力学研究的 由此求得n 一个重要分支。早在2O世纪三四十年代,薄壁杆件的静力学和稳定性理 代入式(1),得 论就逐渐产生了,主要有符拉索夫的开口薄壁梁约束扭转理论、乌曼斯 基的闭口薄壁梁约束扭转理论和符拉索夫的广义坐标法。符拉索夫采用 u x)=(1号)UI+争 两个基本的假定,提出了满足工程需要的实用开口薄壁梁约束扭转理论 或写成: 体系。这一理论自建立以来,尽管本身并不是十分的完美,但有较好的计 ( )=∑ ( ) (2) 算精度,仍然被世界各国所普遍接受,并为更深入的研究奠定了基础,提 出了一套关于闭口薄壁梁约束扭转的实用理论体系。 式中, 经典的薄壁约束扭转理论中已经导出了薄壁杆件及结构的基本微 {I Ⅳ ( )=1一手 ‘ 分方程,但是直接解基本方程求解析解,有时是很困难的,因此,需采用 (2a) 数值方程。数值方法大体上可以分为两大类:一类是直接对微分方程求 l ( )=} 数值解;另一类是通过能量原理求数值解。当直接对微分方程求解时,一 由式(2)可知, ( )就是单元端点位移 =1引起的位移函数.即形 切解微分方程的数值方法都可以使用。最常用的有:数值积分法、差分 函数。 法、加权残数法等。有限元法是基于势能驻值原理的一种数值解法,此方 式(2)又可以写成: 法应用十分广泛。 有限元法的要点是:先把整体拆开,分解成若干个单元,这个过程称 ( ):[Ⅳ ( ) ( )]JL u2 J l (2b) 作离散化。然后再将这些单元按一定的条件集合成整体。在一分一合,先 或简写成: 拆后搭的过程中,把复杂结构的计算问题转化为简单单元的分析和集合 tt(x)=Ⅳ( ) (2c) 问题。有限元法取杆段作为一个单元,单元的刚度方程是按某种假定的 式中,Ⅳ( )称为单元的形函数矩阵: 位移状态推导出的,是近似的表达式,因而求出的是近似解,当单元划分 Ⅳ( )=[Ⅳ ( ) ( )] (3) 细时,近似解趋于精确解。 式(3)就是所求的由节点位移巩确定位移函数u( )的转换关系。凡 1形函数 是单元的位移状态只由两个位移参数决定的位移状态,都可以采用式 (2c)所示的直线式形函数表示。 在有限元位移法中采用节点位移作为基本未知量,为了得到节点之 1.2三次式形函数 间各截面的位移必须假定位移的变化规律,即位移函数。下面分几种情 以横向弯曲单元为例(见图2)推导出单元上任意一点 的横向位移 况来讨论。 ( ).表示为: 1.1直线式形函数 ( )=Ⅳ( ) (4) 以拉压单元为例(见图1)加以说明。由于单元的位移状态由两个端 式中,Ⅳ( )=[ ( ) ( ) ( ) ( )] (5) 点位移参数 和u:来决定,因此位移函数可设为多项式,其中只包括两 c.fF{ l, l , 2, 2 } 140 维普资讯 http://www.cqvip.com 牛庆芳,尹永青薄壁杆件的有限元分析法 y 图2横向弯曲单元 Ⅳl(小1—3等+2手 其中 c 一: ̄ x2 x3 M( )=3等一2等 Ⅳ.( )一 x2+手 , , ,Ⅳ.是4个形函数。 凡是单元的位移状态是由4个位移参数决定的位移状态,都可以采 用式(6)所示的三次式形函数表示。 1_3薄壁单元扭转的形函数 对于薄壁单元的扭转问题,单元上任一截面的扭转位移,可用端点1 和2处的4个端点位移(见图3)来表示。因为它也是4个位移参数决定 的位移状态,故也可用三次式的形函数表示如下: 口( )=Ⅳ( ) (7) 式中,/v(x)同式(6)。 0 2 Z.I———一 0 II V y 图3薄壁单元扭转位移示意图 1.4薄壁空间弯扭单元的形函数 一般的薄壁杆单元,其轴向位移见图1,可用式(2b)表示其位移,即 采用直线式形函数表示轴向位移的变化规律。 一般薄壁空间弯扭单元的变形.包括在:my平面内的横向弯曲位移 ( ),在 平面内的横向弯曲位移w(x),以及绕轴的扭转位移o(x),见 图4。 横向弯曲位移w(x),”( )是由4个位移参数决定的位移状态,因此 采用三次式的形函数。绕 轴的扭转位移o(x),在上节中已表示。将以上 两个方向的横向弯曲位移和绕 轴的扭转位移合写在一起,得一般薄壁 空间弯扭单元节点间位移用端点位移表示的关系式如下: 纠I『Ⅳ( )0 0 1 l I lw(x)I=l 0 N(x)0 l=ld。I (8) )J【0 0Ⅳ( )J【 J 1.5考虑截面翘曲时薄壁空间单元的形函数 在薄壁构件的有限元分析模型中,考虑其约束扭转效应,这样,薄壁 空间杆件单元中每个节点应增加一项表征截面翘曲的位移分量,即截面 的翘曲角 ,以及截面翘曲的影响。 薄壁空间梁单元中每个节点的位移向量共7个,即u W, , , , 比,节点力向量也为7个,即P P P , , ,B,每个位移分量的插值 函数选取为:单元轴向位移函数u采用直线式的位移模式,同时考虑了 双向弯曲和截面翘曲对杆件轴向变形的影响;截面的横向弯曲位移采 , W采用三次式的位移模式;扭角位移函数 , , 也采用三次式的位移 模式,而翘曲位移函数则引入一个表征翘曲沿杆长分布的翘曲函数以 本刊E—mail:bjb@mail.sxinfo.net科技研讨 x 2 0,2 W l 口 l y 图4薄壁空间弯扭单元变形 来代替扭转率 ,合写在一起得单元节点间位移用端点位移表示的关 系式如下: U du lv(x)0 0 0 0 0 0 由 l 0 lv(x)0 0 0 0 0 0 0 lv(x)0 0 0 0 dw 0 0 0 lv(x)0 0 0 d 0 0 0 0 Ⅳ( )0 0 0 0 0 0 0 lv(x)0 dO, 0 0 0 0 0 0 lv(x) dO, dO. 2单元分析 应用势能驻值原理,推导薄壁扭转单元的刚度矩阵与等效节点荷载 向量。 对薄壁弯扭单元(见图5),任意截面 的扭转位移口( )可用节点位移 ds=[O。,0。 , ,02 ] 表示,其关系如下: 口( )=Ⅳ( )ds (9) 单元应变能力 l-争』(EL 伽 )dx (10) 将位移表达式(9)代人上式后,得 (f 1 f{EL( ( )ds)2+GJ ̄(N ( ) ) } 也可改写成以下形式: ,』 I  、㈦=} (f (x)rEJdV'(x)dx+fⅣ (x) I(x)dx ds 可得出单元的刚度矩阵: (。)=I (x) (x) +I N (x)rEJ ̄N (x)dx (11) 5 6 0 l l 由式(6)知 Ⅳ (x)=[一6争+6 ・ 手+3等6V -6fi- 一2 +3等】 (x)=[一 +・ 争-4}+ 一, 争一 }+ 争】 代入式(11),得 (。)=E,_z丑+G zl (12) 141 维普资讯 http://www.cqvip.com 牛庆芳,尹永青薄壁杆件的有限元分析法 本刊E-maihbjb@mail.sxinfo.net科技研讨 式中, 一P 一 £! 一一 一 4一 一 一 “)-}dTz cTd=}drK"d 二 ‘ 式中, 为整体坐标系中的单元刚度矩阵,K z ;而 为单元坐 £!丁 下£!一P 一 一 一 下4一 6 1 1O 15 1 1O f 6 1 1O f 3O 1 1O 2f 15 标转换矩阵。式(2O)可以写成: (13) ∑} ‘ ’d 一 (23) 为使势能表达式中的位移参数得到统一,可将单元杆端位移向量d 转换为结构节点位移向量D,其转换关系式为: d=C‘ D (14) (24) 5f 1 5f 1 1O 6 Zll= 1O 6 式中,矩阵c“ 为节点位移与单元杆端位移向量的转换矩阵。以图6 所示的薄壁杆扭转问题为例,结构分为4个单元,结构节点位移向量D 和单元(3)的杆端位移向量d㈦分别为: D={0l, , ,02 ,03,03 } d㈣={02,02 ,03,03 } (1) (2) (3) (4) ~ 5f 1 1O 30 5f 1 1O 现考虑薄壁杆件弯扭势能。设薄壁杆件上作用着以下外荷载,作用 于截面剪力中心的线分布荷载gI,q ,m,杆端作用力与力N,Q|, , , , 和B,见图5。故总外力势能为: 1 2 3 E,p=一I(q ̄/'+qrr/+mO)dx一[^善+ , Q, 一坳 一^f, +^ +朋 ] (15) ; 图6薄壁杆扭转问题 它们之间的转换关系为: f 0 0 1 0 0 0 1 :式中, , , 分别为截面剪力中心沿x,y, 方向的位移,0为截面的 扭转角,上标“”’表示对:求导数。 因此,单元荷载势能为 f I 0o D 【0 0 0 0 0 1 J “ 一I mOdx一[M,O+BO ] 6 因此,单元(3)的转换矩阵C㈤=[o o】 式中,子块0和,都是2x2阶子矩阵。 将式(24)代入应变能表达式(23)中,即得: 将位移表达式(9)代入后,得 I ‘‘’一I m/V(x)cl ̄lx_[ 鼬 + 0t柑t0 ] 6 (16) v--y, DT(c㈦ ㈦c㈨D=loTy(c㈨) ㈦c㈦D 结构总体刚度矩阵为: 也可改写为以下形式: f “ =一 I mN(x) ̄dx-d 6 6 (17) =∑(c‘ ) ‘ C (25) 将式(17)对位移求导,得等效节点荷载 j 结构的外荷载势能求和得到: pc・)=I mN(x) ̄dx 5 其中,P6 :[肘 B M B2] 为单元的两端荷载向量。 (18) ∑ “ 各单元的外荷载势能 “ 可表示为: f “ 一 I mN(x) ̄dx-d 5 3整体分析 .将式(24)代入后即得: 用有限元法分析任何结构时,其计算都分为两步:第一步,把整体结 构离散为独立的单元进行单元分析;第二步,把这些单元重新集合为整 体,使能满足整体的平衡条件和变形协调条件,进行整体分析。 整体分析是形成结构刚度矩阵和结构节点荷载向量。常用的方法是 刚度集成法。下面从能量原理的角度进行整体分析。 结构的总势能可写成: Ⅱ= (19) ㈨,Up-一∑Dr(C“ ) I ( ) 一∑Dr(C“ ) ・ 等效节点荷载如下: 』 ∑(c‘‘ ) I mN(x) ̄dx 最后指出,式(25)中的(c‘‘ ) “ c‘ 实际上就是扩大的单元刚度矩阵 式中, 为结构的应变能,可由各单元的应变能求和,得到: 引入矩阵c 所起的作用,就是使单元刚度矩阵进行换码和扩大。 有了结构总体刚度矩阵 后,利用各单元的等效节点荷载P“ 计算 U--∑U‘‘ 各单元的应变能U“ 已由上节求得,可表示为: ‘c)= ‘ (2o) 各节点的荷载向量凡( 为节点编号),从而再组成结构的荷载向量R。 在此基础上,结构的平衡方程写成: KD=R (c (21) 式中, 为局部坐标系中单元的杆端位移向量; “ 为局部坐标系中 的单元刚度矩阵。 在整体分析中,各单元有不同的局部坐标系,所以应从符号上加以 说明。 为了求和计算的方便,取统一的整体坐标,将杆端位移经坐标变换 关系式: =z (22) 引人边界条件,求解即可得到节点位移D,再利用式(24)得出单元杆 端位移向量d,通过式 ,得出局部坐标系中的杆端位移向量,最后根据 .p 』_ ( 为局部坐标系中的单元刚度矩阵)计算出单元各节点的内力。 4结语 本文对薄壁结构理论进行了概述,提出了经典理论在实际工程应用 中的不足,有限元法是基于势能驻值原理的一种数值解法,考虑了截面 的翘曲函数,而在薄壁构件中截面的翘曲变形是不容忽略的。 可得: (责任编辑:刘翠玲) 142 维普资讯 http://www.cqvip.com 科技情报开发与经济 文章编号:1005—6033(2006)13-0143—03 SCI-TECH INFORMATION DEVELOPMENT&ECONOMY 2006年第16卷第13期 收稿日期:2006—03-06 侧方交会的精度分析 张文康 (中国铝业山西分公司石灰石矿,山西河津,043300) 摘要:以交会角每l5。的间隔变化及与测角ot, 的相对大小变化组成的45种图形结 构为基础,对侧方交会精度进行了具体分析,提出了一般规范应当修改及修正的正确 数值的建议,总结出了交会角确定时,点位中误差随测角大小变化的规律及侧方交会 时提高精度的具体操作法。 关键词:侧方交会;交会角;点位中误差;精度分析 中图分类号:TD872 ̄.5 文献标识码:A 在中国铝业股份有限公司龙门山石灰石矿日常测量工作中,常采用 经纬仪交会法来加密控制点(简称图根点)。根据已知控制点的条件和地 形情况,经纬仪交会法分为前方交会、侧方交会和后方交会3种。石灰石 矿的测量人员多用侧方交会,因为侧方交会可在待定点直接摆站,进行 验收等测量,从而节约时间,减少工作量。所以本文拟对侧方交会的精度 进行分析,并寻找规律,以期在日常工作中,运用侧方交会中的有效规 律,既能高精度地保证工作质量,又能有效地利用时间,提高工作效率。 一P 般规范规定:第一,选(图根控制)点时,所求点必须和3个以上的 图1侧方交会示意图 级网影响较小可略而不计)。 已知控制点通视。交会角 不小于3O。,不大于150。,即要求最不利的图 形的点位中误差,不大于交会角y=90。时点位中误差的两倍(容许误差)。 第二,选点时,要求点的位置要选在土质坚实、埋设标志后不易损坏和不 会变动的地方。 本文就与规范的第一条规定有关系的侧方交会精度进行分析。 、/ 或,,lp= 盟 (1) (2) 1侧方交会法及其精度 1.1侧方交会 1.2侧方交会的精度 由(1)、(2)式可以看出,除了测角精度对%有影响外,布网图形也 对%有影响。当为定值,y=90。时,%为极小值,即布设交会角y=90。的 图形最为理想,但是在作业中,布设 =9O。的侧方交会图形通常是困难 的,所以和前方交会一样,允许 角在一定范围内变动,即要求 角不小 于3O。或不大于150。,使得最大误差(容许误差)不大于交会角y=90o时 点位中误差的两倍。那么,要求最不利图形的点位中误差不大于交会角 9oo时点位中误差的两倍,是否 角就是不小于3Oo或不大于150o呢? 也就是说, 角究竟在什么样的范围内?我们将对侧方交会的精度进行 大学,助教,太原理工大学,山西省太原市迎泽大街79号.030024 图l为侧方交会示意图。图l中,A, 为两已知点,在已知点A和未 知点P上设站,观测角a和交会角 ,以决定P点位置的方法,叫侧方交 会。侧方交会求未知点的坐标和前方交会求未知点的坐标应用的是相同 的公式。因此,推导侧方交会点位精度公式思路和推导前方交会点位精 度公式的思路是一致的。即先求AP边边长b及方位角a 的真误差公 式,而后求未知点P的坐标真误差公式,再转化成中误差公式.最后导出 P点点位中误差的公式(这里暂时假定已知点的精度较高,它们对低一 第・作者简介:牛庆芳,女,1979年7月生,2002年毕业于太原理工 The Fillite Element Method for the Thin-walled Bar NIU Qing-fang,YIN Yong-qing ABSTRACT:This paper expounds the theory of the thin——walled bar and the deficiency of the theory of the thin——walled confined distortion,makes fu【rther research on the finite element the0r、 —一a kind of numeric solution with very widespread availability and based on potential energy standing—value theory,and points out that in the thin—walled space member bar unit each nodal point should add a displacement component characterizing the cross—sectional warp,the warped angle ,for considering the influence of the cross—sectional warp,and describes that the finite element theory. which perfects the theory of the thin—walled bar,possesses a certain reference signiicance ffor the structural design of the actual engineering. KEY WORDS:thin—walled bar;finite element analysis;cross—sectional warp;displacement component 143 

因篇幅问题不能全部显示,请点此查看更多更全内容