基于ANSYS的水电站厂房刚度中心数值分析.pdf
《基于ANSYS的水电站厂房刚度中心数值分析.pdf》由会员分享,可在线阅读,更多相关《基于ANSYS的水电站厂房刚度中心数值分析.pdf(7页珍藏版)》请在咨信网上搜索。
1、中国农村水利水电China Rural Water and H基于ANSYS的水电站厂房刚度中心数值分析赵廷红,陈健(兰州理工大学能源与动力工程学院,甘肃 兰州 730050)摘要:为了分析水电站厂房的动力特性、求解其刚度中心的位置,在传统刚度中心计算方式的基础上,对计算公式中柔度系数的获取方式做出了改变,修订了水电站厂房刚度中心的计算方法,使刚度中心求解方法的适用性更为广泛。以某水电站厂房为计算模型,应用有限元软件ANSYS进行了刚柔耦合仿真分析,求解得出了该结构的刚度中心,此外还对使用ANSYS求解刚度中心时影响计算结果精确性的因素进行了分析。结果显示:本文修订的刚度中心计算方法结合ANS
2、YS软件不仅可求解得出复杂结构的刚度中心,而且求解速度与精度皆高于传统计算方法,但在求解过程中,网格的不均匀加密将引起计算模型的抗侧刚度发生不均匀变化,对刚度中心数值计算结果的精确性存在一定影响。关键词:水电站厂房;刚度中心;ANSYS;刚柔耦合;数值分析中图分类号:TU311.41 文献标识码:A DOI:10.12396/znsd.222095赵廷红,陈 健.基于ANSYS的水电站厂房刚度中心数值分析 J.中国农村水利水电,2023(8):241-247.DOI:10.12396/znsd.222095.ZHAO T H,CHEN J.Numerical analysis of stiff
3、ness center of hydropower house based on ANSYS J.China Rural Water and Hydropower,2023(8):241-247.DOI:10.12396/znsd.222095.Numerical Analysis of Stiffness Center of Hydropower House Based on ANSYSZHAO Ting-hong,CHEN Jian(College of Energy and Power Engineering,Lanzhou University of Technology,Lanzho
4、u 730000,Gansu Province,China)Abstract:To analyze the dynamic characteristics of the hydropower plant and solve the position of the stiffness center,based on the traditional calculation method of the stiffness center,the method of obtaining the flexibility coefficient in the calculation formula is c
5、hanged,and the calculation method of the stiffness center of the hydropower plant is revised.By taking the powerhouse of a hydropower station as the calculation model,this paper carries out the rigid-flexible coupling simulation analysis by using the finite element software ANSYS,and the stiffness c
6、enter of the structure is obtained.In addition,the factors affecting the accuracy of the calculation results are analyzed when ANSYS is used to solve the stiffness center.The results show that the revised stiffness center calculation method combined with ANSYS software can solve the stiffness center
7、 of complex structure,the solution speed and accuracy are also higher than the traditional calculation method.However,in the solution process,the uneven refinement of the grid will cause the lateral stiffness of the calculation model to change unevenly,which has a certain influence on the accuracy o
8、f the numerical calculation results of the stiffness center.Key words:hydropower house;stiffness center;ANSYS;rigid-flexible coupling;numerical calculation0引 言在结构动力分析的过程中,刚度中心是一项极为重要的参数,大量研究表明1,2它与质量中心的相对位置决定了结构在地震等动力作用下发生平扭耦联等不良现象的剧烈程度,随着二者偏心距离的增大,结构的动位移与动应力也随之增大,控制偏心距可有效避免结构在地震作用下因扭转响应过大而发生扭转破坏。对此
9、,我国相关规范3也对结构的刚度分布进行了一定的约束,从而保证结构在设计层面的合理性,同时防止因建筑物的刚度中心与质量中心之间的距离过大而出现的动文章编号:1007-2284(2023)08-0241-07收稿日期:2022-10-17基金项目:国家自然科学基金项目(51969012)。作者简介:赵廷红(1974-),女,副教授,博士,主要研究方向为水工结构动力分析。E-mail:。水力发电241基于ANSYS的水电站厂房刚度中心数值分析 赵廷红 陈健力响应异常增大、乃至危及结构体系安全的问题。对于单一单层建筑结构,通常直接采用静力学法来求解其刚度中心,即先分析整体结构的刚度分布情况,而后直接应
10、用静力学公式求解出该结构抗侧反力的合力点4,即得出其刚度中心。该方法也可用于求解大型多层复杂结构刚度中心的计算,但其求解过程将非常复杂,而且也很难得到理想的结果。为了克服这种现象,研究者5,6提出了一些求解多层复杂结构刚度中心的方法,这些方法虽然各有特色,但它们的主要思路都是:当整体结构只在某层刚度中心作用水平力时,该层结构只有平动而不发生扭转,而其他层可以有扭转和平动,而后通过观测结构各部分在荷载作用下的响应去求解刚度中心。按照该思路,该层的刚度中心既反映了本层构件的刚度分布又反映了整体结构中其他各层构件的刚度分布情况。目前,水电站厂房抗震研究已有大量成果并逐渐聚焦于流固耦合模拟以及非线性分
11、析等方面,但从整体来看,其地震响应研究是滞后的7,关于水电站厂房平扭耦联效应的分析以及相应的刚度中心计算求解尚未获得关注。随着计算机技术的飞速发展,结构仿真模拟软件也应运而生,结合有限元技术,研究人员可以非常方便地通过数值计算近似地得到结构的刚度中心,但现阶段具备刚度中心自动求解功能的工程软件局限性太大,不适用于除土木工程外其他行业的应用,无法很好利用有限元方法的优势。水电站厂房与土木工程中的框架结构类似,但因其内部设有廊道、流道、孔洞等功能区域,而这些因素使得它的几何外形较土木工程中的框架结构体系要复杂得多,若将其内部的各部分构件按梁、板、柱等形状规则的构件进行简化再用土木行业的数值模拟软件
12、去求解刚度中心,必然使得数值计算的准确性大打折扣,基于此,有必要将刚度中心的求解方法进行扩展。ANSYS是目前应用最为广泛的一款有限元分析软件,可以发挥出有限元方法最大的优势,任意形状的结构均可在精确三维建模后,通过增大结构的离散程度获取可接受的数值结果。本文拟应用ANSYS软件结合现有求解土木工程结构刚度中心的方法来求解水电站厂房的刚度中心;此外,因考虑到有限元网格的划分是影响数值计算结果精度的主要因素之一,本文拟将研究网格的划分对刚度中心计算精度的影响,以期能为水工结构求解刚度中心时提供参考。1刚度中心的计算原理1.1基本假设在求解建筑结构的刚度中心时,必须遵循两条基本假设8:假定楼板在平
13、面内为无限刚性,即受到平面内作用力时仅有刚体位移产生而无弹性变形;刚性楼板上存在这样一个点,当平面内作用力的合力通过此点时,楼板只有平动位移而不发生扭转位移,否则反之。1.2理论公式根据相关文献 6,基于结构整体各层相互影响下刚度中心求解方法的理论推导过程如下:如图1,设某层刚度中心的坐标值为R(xr,yr),在质心或任意非刚度中心点O(x,y)处施加平行于x轴方向的单位力(无量纲数1),记为fx,则刚度中心R(xr,yr)处将随之产生抗侧合力,记为p,则fx与p二者同时形成一组力偶mx,楼板会由此产生扭转位移x,若在P(x,y)处作用反向力偶mz,则楼板由此所产生扭转位移与x方向相反,二者可
14、相互抵消,当mz=mx时,该层楼板的扭转位移归零。基于上述思路,以O(x,y)为原点,列出结构由x、y及三个方向为分量且以柔度矩阵(即刚度矩阵的逆矩阵)表示的受力平衡方程,如下式(扭转位移与力矩以逆时针为正):()fxfymz=()xy(1)也即:()111213212223313233()fxfymz=()xy(2)其中,为结构的柔度矩阵,11至33为对应单位荷载作用下相应的位移,即柔度系数。由前述有:()111213212223313233()10-()yr-y=()xy(3)求解=0,可得下式(力偶与扭转位移以逆时针为正):31-33(yr-y)=0(4)此时注意到:31=x,即上述推导
15、中施加单位力后楼板的扭转位移;而33为单位扭矩作用下楼板的扭转位移,记为z。所以有:x-z(yr-y)=0(5)上式中仅有刚度中心的坐标值yr未知,故可求得。同理,可求得刚度中心的另一个坐标值xr。经整理,考虑结构整体各层相互影响时刚度中心的求解公式如下:yr=y+xz(6)xr=x-yz(7)由于在ANSYS中无法添加单位荷载,所以柔度系数并不能够直接从柔度矩阵中获取,而且对于节点数量庞大的有限元模型而言,直接分析其刚度矩阵的组成也是不现实的。根据柔度系数的定义,x、y与z的值可由点O处作用实际的荷载Fx、Fy及Mz后产生实际扭转位移x、y及z求解其比值获得,所以可将计算公式可变换为:图1计
16、算原理Fig.1 Theory of calculation242基于ANSYS的水电站厂房刚度中心数值分析 赵廷红 陈健yir=yi+i xMizizFi x(8)xir=xi-i yMizizFi y(9)式中:xir、yir为第i层楼板刚度中心的坐标值;xi、yi为第i层楼板荷载作用点的坐标值;ix为第i层楼板荷载作用点沿x轴正方向施加作用力时楼板的扭转位移;iy为在第i层楼板荷载作用点沿y轴正方向施加作用力时楼板的扭转位移;iz为第i层楼板荷载作用点沿逆时针方向施加力矩时楼板的扭转位移;Miz为在第i层楼板荷载作用点沿逆时针方向作用的力矩;Fix为第i层楼板荷载作用点沿x正方向作用的水
17、平力;Fiy为第i层楼板荷载作用点沿y正方向作用的水平力。1.3刚性楼板的处理在数值计算中,不同构件之间刚度差异太大将使得计算结果误差异常增大,甚至引起求解无法收敛的情况9。利用公式(8)、(9)求解结构刚度中心需要将楼板视为刚体,而建筑结构的其余部分皆为柔性体,为克服上述困难,在计算时采用 ANSYS中的刚柔耦合仿真分析方法来进行处理,即通过接触算法将刚体与柔性体接触面上的节点进行耦合,从而将二者连接在一起,该方法广泛应用于机械制造、轨道交通等研究领域,是多体系统动力学与有限元分析方法相结合的表现10。结构刚度中心其实就是约束刚性楼板位移的柔性体所产生抗侧反力的合力点,所以在刚柔耦合分析过程
18、中并不需要进行瞬态动力学分析,这类问题只进行静力学分析便能满足要求11。在几何建模阶段将楼板单独剖分出来进行建模,计算时将其刚度行为设置为刚性,则内部程序将会在楼板的质心处创建一个包含其整体质量、质心坐标值以及转动惯量等动力特性的质量单元,并作为导向节点来表达楼板的刚体运动,故而不必再划分应力应变的网格。但为了让楼板在后处理中也能显示出来,通过设置让程序划分出表示楼板几何形状的网格。刚体与柔性体(楼板与结构主体)之间必须设置接触来进行连接,二者不能合并在一个整体(part)内,而剩余的柔性体(结构主体与地基)将合并为一个整体,使所有柔性体连续、共用节点。由于本研究并不需要进行碰撞、摩擦等高度非
19、线性的接触分析,楼板与结构主体在分析过程中始终保持紧密连接,所以将刚体与柔性体的接触类型设置为绑定连接。刚体与柔性体的接触算法采用多点约束法,该方法是ANSYS众多连接方法中适用性较为广泛、计算精度较高的一种方法,多点约束算法不仅适用于刚柔耦合分析,也可以将不同自由度的单元通过自动创建约束方程联系起来,所以用多点约束法对水电站厂房这类复杂建筑结构进行模拟前处理,将会大大减小数值模拟前处理的工作量12。多点约束法的本质是通过一个或多个节点的自由度为标准值,令其余节点的自由度与之通过数学方程建立耦合关系12,13,如式(10)所示,使用该方法将楼板处理为刚体后如图2所示。Uj+i=1nCi Ui=
20、C0(10)式中:j为主节点的节点编号;i为从节点的节点编号;U为自由度,如位移、温度等,在本文中为位移值;Ci为加权系数;C0为常数项。2工程实例2.1工程概况某水电站是河段上游开发规划中的第六个梯级电站,该水电站以发电为主,同时兼顾防洪、航运等功能。工程规模等级为等大(2)级,按百年一遇的洪水进行设计,千年一遇的洪水进行校核,对应的洪峰流量为30 000和38 100 m3/s,坝址处的多年平均流量732 m3/s,正常蓄水位对应的水库库容为1.72 亿m3。厂内溢流式水电站厂房与混凝土重力坝在结构形式、应力应变行为等方面十分相近14,后文以坝体、坝段进行表述。从右岸至左岸分别为右副坝、航
21、运坝段、泄洪闸坝段、电站厂房坝段、安装间坝段以及左副坝。其中,电站厂房坝段每两个发电机组由两侧的横缝分割为一个机组坝段,每个机组坝段的坝体下部是发电机组的过流通道,坝体上部是表孔溢洪通道,发电机组与表孔泄洪的流道交叠布置,水电站厂房在正常工作阶段有蓄水功能,当洪水来临时启用表孔进行泄洪,过流形式为厂内溢流式,结构体系的质量与刚度分布错综复杂。2.2坐标系设定ANSYS中的坐标系遵循右手定则,全局坐标系原点位于坝体上游迎水面与地基以及左岸相交处,左岸至右岸的垂直方向为x轴正方向、上游至下游的垂直方向为z轴正方向。在全局坐标系中,整体模型 x方向宽 35 m、y方向高 225 m、z方向长 369
22、 m;坝体x方向宽35 m、y方向高75 m、z方向总长71.3 m。为便于求解刚度中心,特建立如图3所示的局部坐标系,其原点位于全局坐标系的(35,70,-4.5)处,右岸至左岸的垂直方向平行于x轴正方向、上游至下游的垂直方向平行于y轴正方向。本文以下所有计算及分析过程均基于该局部坐标系。2.3计算模型在以往与本文研究对象相似的动静力研究中15,16,取典型坝段进行分析是普遍做法,本文基于这种常规的处理方法去求解结构的刚度中心。以一段坝段为研究对象,将坝体顶部的楼板按第一节的方法设置为刚性,采用无质量地基法来考虑结构和坝基的相互作用,坝基沿上、下游及深度方向延伸1.5至2倍坝高,坝体结构密度
23、为 2.5103 kg/m3、弹性模量为 2.81010 Pa、泊松比为0.167;坝基的弹性模量为3109 Pa、泊松比为0.3。坝图2基于多点约束法的刚性楼板Fig.2 Rigid floor based on multi-point constraint method243基于ANSYS的水电站厂房刚度中心数值分析 赵廷红 陈健基底部为全约束,坝基四周为法相约束,由于坝段两侧由横缝分割,本文视坝段间无相互作用,各个坝段独立承受荷载作用15,坝体两侧视为自由面。需要注意坝基的属性尺寸、边界条件的设定等各方面因素均对结构刚度中心的求解结果产生影响。有限元网格的类型对计算结果的精度有所影响,为
24、得到质量较好的有限元网格,将几何模型切分为各个可以扫掠的体,划分为八节点六面体以及六节点五面体的网格单元。远离坝体一定范围外的坝基应力应变程度很小,将此部分网格尺寸增大以减轻计算成本。在实际工程中,对于同一结构体系,数值计算往往需要划分不同的网格方案以获得具有相当精度的计算结果,而在网格整体加密的过程中,结构的应力应变将逐步收敛于真实值,即结构的刚度(柔度)发生了变化,所以刚度中心的坐标值也将随之改变,为分析刚度中心的求解在网格加密过程中的变化,对有限元模型进行整体加密,改变网格尺寸大小以及扫掠路径的单次长短,使坝体与坝基同步加密。本文共设置五组不同的网格方案进行对比分析,各方案的网格参数如表
25、1所示,因篇幅所限,在此只列出方案3的网格划分图(图4)。2.4刚度中心的求解及验证因篇幅所限,在此仅以方案3为例叙述刚度中心的求解以及验证过程。将荷载作用点设为(35,0,0),作用平行于x轴的作用力Fx=1 N、平行于y轴的作用力Fy=1 N、绕z轴的力矩Mz=1 Nm,边界条件的设定与荷载的施加如图5所示。结构变形如图6,在Fx作用下结构发生明显的不均匀位移变形,结构整体向x轴正方向弯曲,受荷载直接作用的上游侧总体位移量比下游侧大,约为1.5至1.8倍,扭转位移强烈;在Fy作用下结构的变形位移规律与Fx作用时基本一致,结构整体向y轴正方向弯曲,受荷载直接作用的左岸侧的总体位移量比右岸侧大
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 ANSYS 水电站 厂房 刚度 中心 数值 分析
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。