《钢-混凝土组合梁离散模型的参数估算》由会员分享,可在线阅读,更多相关《钢-混凝土组合梁离散模型的参数估算(14页珍藏版)》请在金锄头文库上搜索。
1、 指导教师评定成绩(五级制):指导教师签字:附件C:译文钢混凝土组合梁离散模型的参数估算T. Wroblewski, S. Berczyn ski, M. Abramowicz(West Pomeranian University of Technology)摘要:本文介绍了对钢混凝土组合梁的自振数学模型的参数进行估算的方法。该模型使用刚性有限元法(RFE)建立,基于实际组合梁的试验研究结果,来对参数进行估算。估算的参数包括:连接件的剪力和轴向刚度,以及考虑纵筋时,钢筋混凝土板的等效纵向弹性模量。在钢混凝土组合梁模型参数的估算过程中,是否采用了适当的标准,对获得精确的结果有很大影响。选择估算标
2、准的目的是,在不同的振型下,都能使自振频率的试验值与数学模型的计算值一致。组合梁的模型和所有的计算,都是用MATLAB完成的。关键词:刚性有限元模型(RFE);钢混凝土组合梁;估算标准;模型参数;振动1.引言组合结构是把两个或两个以上材料性能不同的构件连接在一起而形成的。钢混凝土梁就是一种组合结构,它由一个工字钢和放在其上的混凝土板组成。在钢筋混凝土板浇筑之前,把钢连接件焊接在钢梁的上翼缘,从而使这两种组件能共同工作。这种组合梁通常用作民用和工业建筑工程中的楼盖,以及桥梁工程中的承重主梁。本文是Berczyski和Wrblewski的论文1和2的一个延伸,在那两篇论文中,他们提出了相关数学模型
3、,并给出了三根钢混凝土组合梁的试验研究成果。那三根组合梁的连接类型各不相同。梁B1和B2(连接件间距不同)通过栓钉连接,梁B3用的则是带孔钢板的刚性连接。他们的研究和分析,旨在确定组合梁的基本动力特性。本文旨在探索出一种组合梁模型参数的估算算法。这种要创建的估算方法基于试验的研究结果,并且,不仅仅使用梁的自振频率一个参数(正如2之前所做的那样),也使用在试验过程中所确定的振型。本文采用RFE方法惯用的平面离散模型进行分析。2.材料和方法2.1问题的分析Biscontin等人的3是最先研究钢-混凝土组合梁动力特性的作品之一。他们根据欧拉梁理论,创造了组合梁的一维连续计算模型,使人们有可能找到组合
4、梁特征值问题的解决方案。在这个模型中,考虑了沿梁轴向的连接柔度,即它的剪切刚度。这个分析模型可用于解释实际组合梁的试验数据。它的数值计算结果与数次试验的结果一致,因而可用于确定梁特定的物理参数和连接的特性。2.2梁的计算模型钢混凝土梁的离散数学模型,是用刚性有限元法RFE惯用的方法建立的。RFE方法是由Kruszewski等人在16和17中提出来的。它把实际结构划分为由弹簧阻尼元件SDS所连接的没有变形的刚性有限元件。建模开始时,把梁分割成相同(或接近)长度的区段。这是初次的划分。接着,在这次划分所形成的每个元素的重心,用一个SDE集中该区段所有的弹簧和阻尼特性。下一步,用RFEs连接所有的S
5、DEs。这是二次划分。在建模时,组合梁的钢和混凝土部分是分别对待的。因此,这才可能对连接参数进行分析。现在常用的组合梁模型示意图见图1。初次划分时,把长度为L的组合梁分成n个长度为L的区段,再在每个区段的中心,放置两个SDE:一个集中了钢区段的特性,而另一个集中了钢筋混凝土板区段的特性。接下来,在这些SDE之间,安置长度为二次划分长度的RFE。端部的划分长度为L/2,其余的为L。最后,添加模拟连接的SDE,连接钢筋混凝土板和工型钢的RFE边缘,就完成了模型。图1.RFE模型:(a)分割的梁段(初次划分); (b)RFE系统(第二次划分)。每一个编号为i的RFE都有它自己独立的坐标系。该坐标系与
6、给定RFE的惯性矩主轴相重叠。在这样假设条件下,质量和惯性矩就是描述RFE所需的唯一参数。这些参数量可以用一个对角质量矩阵的形式表示: (1) 矩阵的前三项等于RFE的质量,后三项是RFE绕轴线的质量惯性矩。对于平面模型,只考虑在二维平面的运动,质量矩阵降维成以下形式: (2) 在这个将结构的连续部分分割成段而建立的模型中,可以从下面的关系式计算出RFE的参数: (3) 其中,(i)RFE材料的质量密度,A(i)是横截面面积; I3i绕轴线X3i的截面惯性矩。钢筋混凝土板或I形钢的截面积A(i)和截面惯性矩I3i都是用来描述RFE的。每个编号为k的SDE都有沿自己主轴线的独立坐标系。SDE的主
7、轴线有这样的性质作用在SDE这些轴线方向的力只导致这个方向的平移变形。 描述SDE的主要参数是其弹簧和阻尼特性的确定系数。对结构自振问题,阻尼特性可以忽略不计。弹簧特性可以用两个矩阵来描述:平动刚度系数矩阵和转动刚度系数矩阵。两个矩阵都是33的对角矩阵。 (4) 这些矩阵,比如质量矩阵,在平面坐标系中分析时可以降维。在确定二维平面上的运动时,矩阵应采用以下形式: (5) 钢筋混凝土板和工型钢的SDE模型的平移和旋转刚度系数,可以用下面的关系式计算: (6) 其中,E(k)杨氏弹性模量,G(k)剪切模量(Kirchhoff模量);(k)考虑非均匀切向应力模式的横截面形状系数(Timoshenko
8、剪切系数)。系数的计算方法见文献1。SDE的连接模型矩阵中的元素等于剪切连接刚度Kh和轴向连接刚度Kv。在这种情况下,假设旋转刚度系数为零。刚度矩阵K是在平动刚度系数和转动刚度系数的基础上建立的。惯性矩阵M是在独立的质量矩阵的基础上建立的。矩阵建立的方法详见由Kruszewski等人的16和Wittbrodt等人的15。从一般的运动微分方程得到的自由振动的微分方程,在忽略了外界的对梁作用,并忽略了阻尼效应之后,可以写成以下形式: (7) 其中,q广义坐标向量。以上确定自由振动的频率和对应振型的方程解法,在文献18中有详细描述。2.3模型参数的估算需估算的参数包括:(i)连接的剪切刚度Kh,(i
9、i)轴向刚度Kv(垂直于梁的轴向方向)和(iii)钢筋混凝土楼板在考虑了纵向钢筋时的等效弹性模量Ec。模型中使用的其他参数均取自文献或梁的数据库。指数S是在试验得到的前n个弯曲振动频率的相对差的平方和。 (10) 其中,n所选的频率数量。指数Z是MAC的平均补值。因此,在进行优化时,考虑到了不同振型。 (11) MAC(模态保证准则)是使试验方法和数值模拟计算出的振型相符合的一个标准。当MAC值接近1时,振型结果匹配得更好。 (12) 参数估算时采用以下形式的优化标准:J=wfS+wZ (13) 这种方法允许同时优化指数S和Z。在估算的过程中,考虑到频率和振型,对频率(指数S)和振型(指数Z)
10、分别采取适当的权重函数wf和w。要使用这样的权重函数,需确定在估算的过程中振动频率和振型的重要性。类似的解决方法在Mordini的12中提出。已采用的附加标准,是纵向振动基本频率的试验值与数学计算值的兼容条件(14)。该算法假定的限制条件(14),对确定楼板弹性模量EC的影响最大。 (14) 为了确定能让振型和频率与试验结果有最佳映射关系的模型参数,指数J已经最小化了。为了解决这个问题,在MATLAB环境下实现优化程序。3.应用实例3.1组合梁研究所使用的组合梁是由S235JRG2钢制成的IPE160热轧工字钢,连接到一块由C25/30混凝土制成的横截面为60600mm的混凝土板。梁B1和B2
11、的韧性连接是,由S235J2G3钢制造的,直径为10mm,高为50mm的KOCO-SD型双头螺柱。在梁B1和B2中,分别间隔150mm和250mm放一对螺栓。梁B3通过穿孔的钢条板进行刚性连接。条板是由S235JRG2钢制成的厚度为10mm的钢板。条板通过4mm厚的角焊缝焊接于梁上。梁的总长是3200mm。三条组合梁中的连接件分布见图1、图2和图3。图2.组合梁中连接件分布: (a) B1; (b) B2; (c) B3.图3.组合梁剖面图3.2试验研究试验台由两个钢架组成。动态测试时,通过四束钢丝将梁悬挂在钢架上,形成一个可以自由运动的梁。试验台和悬挂梁见图4。图4.试验台和悬挂梁进行测试时
12、,使用脉冲激励型号为KISLER9726A20000的模型锤。通过KISLER 8692C5M1型三向轴压电传感器记录加速度。记录加速度的测量点,沿梁长(工型梁的底部凸缘有9个测量点和混凝土板的上表面有三行共27个测量点)每间隔375mm均匀分布。有关测试的详细描述见1。在分析钢筋混凝土或组合结构的动态特性时,应考虑动态的纵向弹性模量Ed,而不是静态的割线模量Ecm。模量Ed可以通过下式估算:Ed=0.8Ecm+C (15) 其中,C是等于19GPa的常数。这个表达式借用了英国标准:结构用混凝土CP110:1972。这个关系式中,研究使用的混凝土的动态弹性模量为33.1GPA。在此基础上,钢筋
13、混凝土板的等效弹性模量Ec考虑了板中纵筋的影响,它的计算值是36.0GPA。3.3模型描述现在还没有专业的RFE计算软件,于是我们在MATLAB环境中编写源程序来解决这个问题。该程序能解决组合梁的自振问题,并控制梁的所有其它参数。对梁B1和B2,采用的划分方法可以使得放置的SDE能体现RFE对称轴上螺栓的属性。 梁B1被分成22段:中间20段长150mm,两端的长100mm。模型最后由46个SDE(23个钢筋混凝土板的SDE和23个钢梁的SDE)和65个RFE组成。梁模型见图5。图5.梁B1的RFE模型:(a)梁的区段划分 (b)梁的RFE划分梁B2被分成26个段:中间24段长125mm,两端
14、的长100mm。在这种情况下,由于取区段数n=20,所以要增加每段螺栓的分布密度1。这种RFE密度划分方法,使得有代表性的连接螺栓放在每间隔一个SDE的位置。梁B3的连接可以视为均匀放置,而不是由连接板条的分布密度确定初次划分。因而梁被分为32段,每段100mm,模型最终由6 6个通过99个RFE连接的SDE组成。3.4估算结果采用标准(13),进行了以下的数值试验。在第一阶段的参数估算,标准(13)取最小值,假设wf = 1.0, w= 0.0。这意味着,在估算过程中只对自振频率进行比较。从估算结果可以看到,自振频率的试验值和数值计算结果的一致性较好(见表1)。梁B1和B3相差在1以下。对于梁B2,最大相差1.1。表1.自振频率的动态试验结果与数值计算结果比较第一次估算,取wf= 1.0,w=0.1通过引入不同的初始值,使用收敛迭代算法进行测试,来确定模型的参数。每次迭代得到的变量估算结果都很相似。图6给出了梁B1的第四振型。使用S表示钢部分和混凝土部分