毕业设计译文(土木工程专业).doc

《毕业设计译文(土木工程专业).doc》由会员分享，可在线阅读，更多相关《毕业设计译文(土木工程专业).doc（9页珍藏版）》请在得力文库 - 分享文档赚钱的网站上搜索。

1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流毕业设计译文(土木工程专业).精品文档.指导教师评定成绩(五级制)：指导教师签字：附件C：译文钢混凝土组合梁离散模型的参数估算T. Wroblewski, S. Berczyn ski, M. Abramowicz（West Pomeranian University of Technology）摘要:本文介绍了对钢混凝土组合梁的自振数学模型的参数进行估算的方法。该模型使用刚性有限元法（RFE）建立，基于实际组合梁的试验研究结果，来对参数进行估算。估算的参数包括：连接件的剪力和轴向刚度，以及考虑纵筋时，钢筋混凝土板的等效纵向弹性模量。在钢混凝

2、土组合梁模型参数的估算过程中，是否采用了适当的标准，对获得精确的结果有很大影响。选择估算标准的目的是，在不同的振型下，都能使自振频率的试验值与数学模型的计算值一致。组合梁的模型和所有的计算，都是用MATLAB完成的。关键词：刚性有限元模型（RFE）；钢混凝土组合梁；估算标准；模型参数；振动1.引言组合结构是把两个或两个以上材料性能不同的构件连接在一起而形成的。钢混凝土梁就是一种组合结构，它由一个工字钢和放在其上的混凝土板组成。在钢筋混凝土板浇筑之前，把钢连接件焊接在钢梁的上翼缘，从而使这两种组件能共同工作。这种组合梁通常用作民用和工业建筑工程中的楼盖，以及桥梁工程中的承重主梁。本文是Bercz

3、yski和Wrblewski的论文1和2的一个延伸，在那两篇论文中，他们提出了相关数学模型，并给出了三根钢混凝土组合梁的试验研究成果。那三根组合梁的连接类型各不相同。梁B1和B2（连接件间距不同）通过栓钉连接，梁B3用的则是带孔钢板的刚性连接。他们的研究和分析，旨在确定组合梁的基本动力特性。本文旨在探索出一种组合梁模型参数的估算算法。这种要创建的估算方法基于试验的研究结果，并且，不仅仅使用梁的自振频率一个参数（正如2之前所做的那样），也使用在试验过程中所确定的振型。本文采用RFE方法惯用的平面离散模型进行分析。2.材料和方法2.1问题的分析Biscontin等人的3是最先研究钢-混凝土组合梁动

4、力特性的作品之一。他们根据欧拉梁理论，创造了组合梁的一维连续计算模型，使人们有可能找到组合梁特征值问题的解决方案。在这个模型中，考虑了沿梁轴向的连接柔度，即它的剪切刚度。这个分析模型可用于解释实际组合梁的试验数据。它的数值计算结果与数次试验的结果一致，因而可用于确定梁特定的物理参数和连接的特性。2.2梁的计算模型钢混凝土梁的离散数学模型，是用刚性有限元法RFE惯用的方法建立的。RFE方法是由Kruszewski等人在16和17中提出来的。它把实际结构划分为由弹簧阻尼元件SDS所连接的没有变形的刚性有限元件。建模开始时，把梁分割成相同（或接近）长度的区段。这是初次的划分。接着，在这次划分所形成的

5、每个元素的重心，用一个SDE集中该区段所有的弹簧和阻尼特性。下一步，用RFEs连接所有的SDEs。这是二次划分。在建模时，组合梁的钢和混凝土部分是分别对待的。因此，这才可能对连接参数进行分析。现在常用的组合梁模型示意图见图1。初次划分时，把长度为L的组合梁分成n个长度为L的区段，再在每个区段的中心，放置两个SDE：一个集中了钢区段的特性，而另一个集中了钢筋混凝土板区段的特性。接下来，在这些SDE之间，安置长度为二次划分长度的RFE。端部的划分长度为L/2，其余的为L。最后，添加模拟连接的SDE，连接钢筋混凝土板和工型钢的RFE边缘，就完成了模型。图1.RFE模型：(a)分割的梁段（初次划分）；

6、 （b）RFE系统（第二次划分）。每一个编号为i的RFE都有它自己独立的坐标系。该坐标系与给定RFE的惯性矩主轴相重叠。在这样假设条件下，质量和惯性矩就是描述RFE所需的唯一参数。这些参数量可以用一个对角质量矩阵的形式表示： (1) 矩阵的前三项等于RFE的质量，后三项是RFE绕轴线的质量惯性矩。对于平面模型，只考虑在二维平面的运动，质量矩阵降维成以下形式： (2) 在这个将结构的连续部分分割成段而建立的模型中，可以从下面的关系式计算出RFE的参数： （3) 其中，（i）RFE材料的质量密度，A（i）是横截面面积; 绕轴线的截面惯性矩。钢筋混凝土板或I形钢的截面积A（i）和截面惯性矩都是用来描

7、述RFE的。每个编号为k的SDE都有沿自己主轴线的独立坐标系。SDE的主轴线有这样的性质作用在SDE这些轴线方向的力只导致这个方向的平移变形。 描述SDE的主要参数是其弹簧和阻尼特性的确定系数。对结构自振问题，阻尼特性可以忽略不计。弹簧特性可以用两个矩阵来描述：平动刚度系数矩阵和转动刚度系数矩阵。两个矩阵都是33的对角矩阵。 （4) 这些矩阵，比如质量矩阵，在平面坐标系中分析时可以降维。在确定二维平面上的运动时，矩阵应采用以下形式： （5) 钢筋混凝土板和工型钢的SDE模型的平移和旋转刚度系数，可以用下面的关系式计算： （6) 其中，E（k）杨氏弹性模量，G（k）剪切模量（Kirchhoff模

8、量）;（k）考虑非均匀切向应力模式的横截面形状系数（Timoshenko剪切系数）。系数的计算方法见文献1。SDE的连接模型矩阵中的元素等于剪切连接刚度Kh和轴向连接刚度Kv。在这种情况下，假设旋转刚度系数为零。刚度矩阵K是在平动刚度系数和转动刚度系数的基础上建立的。惯性矩阵M是在独立的质量矩阵的基础上建立的。矩阵建立的方法详见由Kruszewski等人的16和Wittbrodt等人的15。从一般的运动微分方程得到的自由振动的微分方程，在忽略了外界的对梁作用，并忽略了阻尼效应之后，可以写成以下形式： (7) 其中，q广义坐标向量。以上确定自由振动的频率和对应振型的方程解法，在文献18中有详细描

9、述。2.3模型参数的估算需估算的参数包括：（i）连接的剪切刚度Kh，（ii）轴向刚度Kv（垂直于梁的轴向方向）和（iii）钢筋混凝土楼板在考虑了纵向钢筋时的等效弹性模量Ec。模型中使用的其他参数均取自文献或梁的数据库。指数S是在试验得到的前n个弯曲振动频率的相对差的平方和。 （10) 其中，n所选的频率数量。指数Z是MAC的平均补值。因此，在进行优化时，考虑到了不同振型。 （11) MAC（模态保证准则）是使试验方法和数值模拟计算出的振型相符合的一个标准。当MAC值接近1时，振型结果匹配得更好。 （12) 参数估算时采用以下形式的优化标准：J=wfS+wZ （13) 这种方法允许同时优化指数S

10、和Z。在估算的过程中，考虑到频率和振型，对频率（指数S）和振型（指数Z）分别采取适当的权重函数wf和w。要使用这样的权重函数，需确定在估算的过程中振动频率和振型的重要性。类似的解决方法在Mordini的12中提出。已采用的附加标准，是纵向振动基本频率的试验值与数学计算值的兼容条件（14）。该算法假定的限制条件（14），对确定楼板弹性模量EC的影响最大。 （14) 为了确定能让振型和频率与试验结果有最佳映射关系的模型参数，指数J已经最小化了。为了解决这个问题，在MATLAB环境下实现优化程序。3.应用实例3.1组合梁研究所使用的组合梁是由S235JRG2钢制成的IPE160热轧工字钢，连接到一块

11、由C25/30混凝土制成的横截面为60600mm的混凝土板。梁B1和B2的韧性连接是，由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试验研究试验台由两个钢架组成。动态测试时,通过四束钢丝将梁悬挂在钢架上，形成一

12、个可以自由运动的梁。试验台和悬挂梁见图4。图4.试验台和悬挂梁进行测试时,使用脉冲激励型号为KISLER9726A20000的模型锤。通过KISLER 8692C5M1型三向轴压电传感器记录加速度。记录加速度的测量点，沿梁长（工型梁的底部凸缘有9个测量点和混凝土板的上表面有三行共27个测量点）每间隔375mm均匀分布。有关测试的详细描述见1。在分析钢筋混凝土或组合结构的动态特性时，应考虑动态的纵向弹性模量Ed,而不是静态的割线模量Ecm。模量Ed可以通过下式估算：Ed=0.8Ecm+C (15) 其中,C是等于19GPa的常数。这个表达式借用了英国标准：结构用混凝土CP110：1972。这个关

13、系式中,研究使用的混凝土的动态弹性模量为33.1GPA。在此基础上，钢筋混凝土板的等效弹性模量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）梁的区段划分

14、 （b）梁的RFE划分梁B2被分成26个段：中间24段长125mm，两端的长100mm。在这种情况下，由于取区段数n=20,所以要增加每段螺栓的分布密度1。这种RFE密度划分方法,使得有代表性的连接螺栓放在每间隔一个SDE的位置。梁B3的连接可以视为均匀放置，而不是由连接板条的分布密度确定初次划分。因而梁被分为32段，每段100mm，模型最终由6 6个通过99个RFE连接的SDE组成。3.4估算结果采用标准（13），进行了以下的数值试验。在第一阶段的参数估算，标准（13）取最小值，假设wf = 1.0, w= 0.0。这意味着，在估算过程中只对自振频率进行比较。从估算结果可以看到,自振频率的试

15、验值和数值计算结果的一致性较好（见表1）。梁B1和B3相差在1以下。对于梁B2,最大相差1.1。表1.自振频率的动态试验结果与数值计算结果比较第一次估算，取wf= 1.0，w=0.1通过引入不同的初始值,使用收敛迭代算法进行测试，来确定模型的参数。每次迭代得到的变量估算结果都很相似。图6给出了梁B1的第四振型。使用S表示钢部分和混凝土部分的垂直位移差。组合梁这两个部分沿垂直方向的位移差，是由连接在该方向的变形造成的，这将影响参数Kv的确定。这些差异在数学模型和进行的试验测试都有所表现。分析图6可以清楚地看到，在试验测试结果中，钢和混凝土在垂直方向上的协同工作情况比数值模型要好。因此，可以假定，

16、刚度Kv确定得太小，不能代表钢和混凝土在测试过程中观察到的协同水平（图6a）。图6.梁B1在第4振型时,试验和数值结果的S比较：（a）试验数据（b）振型在试验和数值的S比较（c）估算的数值数据第二阶段,取权重wf0.6, w= 1.0时,估算结果如下。当检查表2时,有人指出，对Eq引入MAC（13）导致频率的结果相差较大,并且,确定出的参数KV偏大。这意味着，在估算过程中,考虑了振型之后，更能精确地描述连接沿垂直方向的行为。表2.自振频率的动态试验结果与数值计算结果比较第二次估算，取wf0.6, w= 1.0在估算过程中，考虑了振型对MAC值的影响后，提高了问题中振型的一致性，如图7所示。从图

17、7（b）可以显示了，第一次和第二次数值估算的值与试验值的S。分析图7（b）可以清楚地看出，梁的第二次估算值的最大值与试验数据的最大值都差不多，更好地呈现了试验测试和RFE模型的一致性。在模型参数的估算中，引入MAC系数大大提高了振型的一致性，而略减小了（改变了2-3左右）自振频率的一致性。这种程度的一致性，似乎是频率和振型在解析模型与观察试验之间一致性的一个合理折衷。图7.第二次估算中，梁B1第4振型的试验和数值数据的S比较：（a）试验数据（b）振型在试验和数值的s比较（c）第二次估算的数值数据4.讨论我们经过分析，发现了当已经估算出的参数Kh、Kvi和Ec变化时，对所获得的结果的影响。我们对

18、每一个待定的参数，都进行了单独的分析。本文包含了梁B1的分析结果。参数是在估算的第二阶段确定的，即，将连接的刚度和钢筋混凝土板的动态弹性模量假定为初始值，表示为Kh0，Kv0和Ec0。在计算时，首先，把剪切刚度Khi作为唯一变量，让它从初始值Kv0的30到170范围内变化。然后，是轴向刚度连接Kv，也在相同范围内进行变化。最后，等效弹性模量Ec也从其初始值的60到140范围内变化。所得结果列在图8、图9和图10。图8.梁B1的自振频率随刚度Kh的变化情况图9.梁B1的自振频率随刚度Kv的变化情况图10.梁B1的自振频率随刚度弹性模量Ec的变化情况图8中曲线的分析结果表明，当连接的剪切刚度Kh增

19、加时,会增加所有参与分析的弯曲振动的频率。另一方面，轴向振动的频率不会受到Kh变化的影响。在轴向刚度KV变动时,可以观察到,自然弯曲振动的敏感度和频率之间有明确的关联。当刚度变化时,自然弯曲振动的频率越高，它就越敏感。轴向振动的基本频率，正如早期所观察到的一样,对刚度Kv的变化不敏感。图8和图9的曲线图的比较揭示，低频率的弯曲振动对刚度Kh的变化更敏感，而高频率的对轴向刚度KV的变化更为敏感。这是由于，弯曲振动频率的增加,伴随着使梁的钢和混凝土部分在垂直方向上相分离的惯性力的增加。这种情况下，当钢筋混凝土板的等效纵向弹性模量Ec改变（见图10）时，轴向振动的频率f1long对其最为敏感。这证明

20、,在参数估算过程中考虑标准（14）是一个正确的决定。参数Ec的确定起到了决定性的作用。应该指出的是,等效弹性模量的平均估算值，达到了Ec=33.8 GPA，略低于混凝土材料的试验值，即Ec =36.0 GPA（见第3.2章）。这百分之几的差异,是由测定模量Ed时使用Eq的简化方法（15）引起的。因为,对于水泥浆含量较高的混凝土,由于使用的混凝土总厚度不大于8mm,该方程中包含的常数C取19 GPa并不恰当。5.结论数次测试结果表明，该钢-混凝土组合梁的数学模型参数估算算法具有很好的收敛性。可以用模型代表实际的组合梁，并确定与试验过程中相似的动力特性。对组合梁模型的参数进行估算时，选择适当的权重

21、函数，对结果的准确性有很大影响。我们旨在使测试中观察得到的频率和振型与分析模型计算的结果一致。对于组合梁模型的参数，认定标准的选择是否恰当，对结果的精度会显著的影响。改变标准并确定标准成分的灵敏度之后，振型将有更好的一致性。然而，这种改变会使得自振频率的一致性有所降低。参考文献1 S. Berczyn ski, T. Wroblewski, 钢混凝土组合梁在Timoshenko梁模型中的振动研究，振动与控制学报11 (2005) 829848.2 S. Berczyn ski, T. Wroblewski, 钢混凝土组合梁的自振振型的试验验证，振动与控制学报16 (2010) 20572081

22、.3 G. Biscontin, A. Morassi, P. Wendel, 钢混凝土组合梁的振动, 振动与控制学报6 (2000) 691714.3 G. Biscontin, A. Morassi, P. Wendel, 钢混凝土组合梁的振动, 振动与控制学报6 (2000) 691714. 4 M. Dilena, A. Morassi, 钢混凝土组合梁连接处部分损坏的试验模型分析, 振动与控制学报 10 (2004) 897913.5 M. Dilena, A. Morassi, 振动在钢混凝土组合梁连接部分的退化和损伤检测中的应用，声音与振动学报320 (2009) 101124.

23、6 E.J. Sapountzakis, 对钢混凝土组合结构的可变形连接的动态分析, 计算机与结构 82 (2004) 717729.7 E.J. Sapountzakis, V.G. Mokos, 平行梁加劲板动态分析的改进模型, 工程结构30 (2008) 17201733.8 E.J. Sapountzakis, I.C. Dikaros, 平行梁加劲板的大挠度分析, 工程结构35 (2012) 254271.9 J. Maeck, M. Abdel Wahab, B. Peeters, G. De Roeck, J. De Visscher, W.P. De Wilde, J.-M.

24、Ndambi, J. Vantomme,识别钢筋混凝土结构损伤的动态刚度测定法, 工程结构22 (2000) 13391349.10 A. Teughels, J. Maeck, G. De Roeck, 有限元模型更新损伤函数进行损伤评估, 计算机与结构 80 (2002) 18691879.11 S. Maas, A. Zu rbes, D. Waldmann, M. Waltering, V. Bungard, G. De Roeck, 试验室动态测试法对混凝土结构的损伤评估, 工程结构 34 (2012) 351362.12 A. Mordini, K. Savov, H. Wenze

25、l, 有限元模型的修正：一个功能强大的结构健康监测工具, 国际结构工程4 (2007) 352358.13 A. Mordini, H. Wenzel, 基于VC + +的梁结构损伤检测, 结构工程电子期刊10 (2010) 1121.14 R.J. Allemang, 模态保证准则（MAC）：二十多年的使用和滥用, 声音与振动 37 (2003) 1421.15 E. Wittbrodt, I. Adamiec-Wojcik, S. Wojciech, 柔性多体系统动力学.刚性有限元法, 施普林格出版社，柏林，海德堡，纽约, 2006. 16 J. Kruszewski等, 刚性有限元法（波

26、兰）, 华沙，阿尔卡季, 1975.17 J. Kruszewski等,结构动力学的刚性有限元法（波兰），WNT，华沙, 1999.18 M. Rucka, K. Wilde,在MATLAB环境下的结构动力学实例（波兰）, WPG, 格但斯克, 2008.译文原文出处：T. Wrblewski, S. Berczyski, M .Abramowicz . Estimation of the parameters of the discrete model of a steelconcrete composite beamJ. elsevier, June 2013, 13(2):209219.