可靠度大作业

工程结构可靠度分析学院：姓名：学号：指导教师：目录一、工程概况..............................................................................................................1二、计算依据和参考..................................................................................................3三、分析模型和参数..................................................................................................43.1计算模型..........................................................................................................43.2计算模型参数..................................................................................................5四、可靠指标计算....................................................................................................124.1MonteCarlo方法简介...................................................................................124.2响应面方法简介............................................................................................124.3可靠指标的计算............................................................................................13附录：ANSYSAPDL命令流....................................................................................251一、工程概况冷水江市资江二桥为一独塔双索面、完全对称、塔梁墩固结的斜拉桥，设置有22对、上下游2个索面斜拉索，共计88根。试考虑部分参数的随机性，应用ANSYSPDS模块计算在各施工工况（最多悬臂施工至第10对索已经张拉），由于突发断索导致的主梁、索塔、桥墩、斜拉索抗弯、抗拉、抗裂性等指标失效概率的变化（响应面法+MonteCarlo）。主梁标准断面如图1-1所示，索塔断面如图1-2所示。A-A342252,5图1-1主梁标准断面（单位：m）Ⅲ-ⅢⅣ-Ⅳ909010010022090909047090140140909040022090909047090400R90R90图1-2-a索塔空心断面（单位：cm）2Ⅴ-ⅤⅥ-Ⅵ2209090904709040032090909057090500R90R90图1-2-b索塔实心断面（单位：cm）题目编号2：工况：断索位置，仅一侧（9#拉索）；计算内容（可靠指标）：拉索抗拉能力（8#、10#），索塔抗弯、抗裂性。3二、计算依据和参考1、《公路工程技术标准》JTGB01-20032、《公路桥涵设计通用规范》JTGD60-20043、《公路钢筋混凝土及预应力混凝土桥涵设计规范》JTGD62-20044、《公路工程结构可靠度设计统一标准》GB/T50283-19994三、分析模型和参数3.1计算模型主梁、索塔、桥墩建模时可考虑尽量采用标准节段，主梁横隔板按节点荷载计入。拉索锚固点塔端以钢锚梁上锚固坐标为准，锚点与建模节点间通过刚臂连接。索塔内钢锚梁不必建模，按节点荷载计入。主梁预应力必须模拟，尽量按直线段处理。主梁按全预应力结构设计，墩塔按钢筋混凝土结构设计，塔梁墩抗裂性要求按规范D62进行。挂篮作用在主梁最前端，重600kN。为缩短计算时间，可以按10个节段悬臂施工完毕后的索力按结构一次形成的方式施加到主梁及索塔上，然后进行后续分析。主梁、索塔及桥墩采用beam188建模，钢横梁采用beam4建模，斜拉索采用link10建模。根据设计图纸采用通用计算软件ANSYS14.5建立结构有限元模型，模型如图3-1所示。图3-1ANSYS14.5结构模型图53.2计算模型参数考虑部分参数的随机性，如表3-1所示。参数统计特征表3-1编号参数均值变异系数分布1施工索力T12707kN0.05正态2施工索力T22352kN0.05正态3施工索力T32506kN0.05正态4施工索力T42715kN0.05正态5施工索力T52898kN0.05正态6施工索力T63011kN0.05正态7施工索力T73284kN0.05正态8施工索力T83434kN0.05正态9施工索力T93588kN0.05正态10施工索力T103694kN0.05正态11混凝土弹性模量E3.5E7kPa0.1正态12混凝土密度26.0kN/m30.1对数正态标准差与均值及变异系数之间的关系为：xxxV，正态分布又称为高斯分布，各参数的统计特征如图3-2所示。6图3-2-1F1统计特征图3-2-2F2统计特征7图3-2-3F3统计特征图3-2-4F4统计特征8图3-2-5F5统计特征图3-2-6F6统计特征9图3-2-7F7统计特征图3-2-8F8统计特征10图3-2-9F9统计特征图3-2-10F10统计特征11图3-2-11YOUNG统计特征图3-2-12DENSITY统计特征12四、可靠指标计算4.1MonteCarlo方法简介MonteCarlo法又名随机模拟法或统计试验法，将随机变量ix代入功能函数)(xg，得出一个函数值。若0)(xg，则在计算机程序中计入一次失效的实现；若0)(xg，则不计入。这样就完成了一次计算，再产生下一个随机数，重复上面的计算，直至完成预定的实验次数为止。此时，失效概率为nkxgPPnflim]0)([MonteCarlo法可选择直接抽样法、拉丁超立方抽样（LHS）和用户抽样处理，本例中采用拉丁超立方抽样（LHS）。其特点如下：①比直接抽样法更先进、更有效。②对抽样过程有“记忆”功能，可避免直接抽样法数据点集中而导致的仿真循环重复问题。③强制抽样过程中抽样点必须离散分布于整个抽样空间。④LHS抽样法比直接抽样法要少20%~40%的仿真循环资料。⑤需要指定仿真循环次数、重复次数、样本分布位置、循环终止准则（均值和标准方差精度等）和随机输入参数样本种子值。4.2响应面方法简介其基本思想就是选用一个适当的具有明确表达式的函数来近似代替一个不能明确表达的函数，对于可靠度分析来说，就是尽可能通过一系列确定性的试验即有限元数值计算结果来拟合一个响应面（明确表达式的函数）以代替未知的真实的极限状态曲面。给定一组结构性能参数，几何参数和荷载的取值，应用确定的有限元数值计算，就可以得到结构的一个响应值，取n组结构性能参数，几何参数和荷载的取值，就能得到相应的n个响应值，根据这n个响应值拟合的函数表示的曲面就叫响应面，以该响应面代替未知的真实的极限状态曲面进行可靠度分析的方法就响应面方法。13响应面的设计实质也就是响应面函数形式的确定。响应面函数形式的确定应满足的两个要求：（1）响应面函数数学表达式在基本能够描述真实函数的前提下应尽可能的简单，以方便可靠度分析；（2）响应面函数中的待定系数应尽可能的少，以便减少需要确定待定系数的结构有限元分析的工作量。本例取不含交叉项二次多项式为响应函数：211iniiiniiXcXbaZ式中：iX为基本随机变量；ia、ib、ic为响应面函数中的待定系数。4.3可靠指标的计算4.3.1拉索抗拉能力计算斜拉索抗力均值为1581MPa，变异系数为0.05，服从对数正态分布。9#斜拉索断裂后，8#和10#斜拉索应力水平的特征值如图4-1、4-2所示。MEAN0.50114E+09STDEV0.55247E+08SKEW0.81320E-01KURT0.16164E+01MIN0.35560E+09MAX0.67221E+09ConfidenceLimit95.00%.420E+09.440E+09.460E+09.480E+09.500E+09.520E+09.540E+09.560E+09.580E+09.600E+09.620E+09113263850SMAX8NumberofsamplesResultSetLHSLSJZJMeanValuesofSamples图4-1-18#斜拉索应力均值14MEAN0.50114E+09STDEV0.55247E+08SKEW0.81320E-01KURT0.16164E+01MIN0.35560E+09MAX0.67221E+09ConfidenceLimit95.00%.250E+08.375E+08.500E+08.625E+08.750E+08.875E+08.100E+09.113E+09.125E+09.138E+09.150E+09113263850SMAX8NumberofsamplesResultSetLHSLSJZJStandardDeviationofSample图4-1-28#斜拉索应力标准差MEAN0.50114E+09STDEV0.55247E+08SKEW0.81320E-01KURT0.16164E+01MIN0.35560E+09MAX0.67221E+09ConfidenceLimit95.00%0102030405060708090100.350E+09.375E+09.400E+09.425E+09.450E+09.475E+09.500E+09.525E+09.550E+09.575E+09.600E+09.625E+09.650E+09.675E+09Probabilityin%SMAX8ResultSetLHSLSJZJCumulativeDistributionFunction图4-1-38#斜拉索应力累积分布图15Significant:DENSITYF8Insignificant:F1F2F3F4F5F6F7F9F10YOUNGSignificancelevel:2.500%0.1.2.3.4.5.6.7.8.9ResultSetLHSLSJZJOutputParameterSMAX8Rank-OrderCorrelationSensitivities图4-1-48#斜拉索应力概率参数的灵敏度图MEAN0.50114E+09STDEV0.55247E+08SKEW0.81320E-01KURT0.16164E+01MIN0.35560E+09MAX0.67221E+090.05.1.15.2.25.3.35.4.45.5.356E+09.408E+09.461E+09.514E+09.567E+09.619E+09.672E+09RelativeFrequencySMAX8ResultSetLHSLSJZJHistogram图4-1-58#斜拉索应力概率参数柱状分布图16MEAN0.55876E+09STDEV0.57646E+08SKEW0.31592E+00KURT-0.39735E+00MIN0.44707E+09MAX0.69391E+09ConfidenceLimit95.00%.460E+09.480E+09.500E+09.520E+09.540E+09.560E+09.5