基于蒙特卡罗方法盾构刀盘结构可靠性分析.docVIP

基于蒙特卡罗方法盾构刀盘结构可靠性分析 　　摘要：盾构刀盘承担开挖土体作用，在整个掘进过程中刀盘工作环境恶劣，受力情况比较复杂，在考虑输入参数（如刀盘载荷、材料属性等）不确定性因素的基础上，对刀盘可靠性进行了研究，研究刀盘在不确定性因素的作用下刀盘结构的可靠性。该文首先通过三维建模软件Pro/E对盾构机刀盘进行建模，然后将模型导入到有限元软件ANSYS中，采用蒙特卡罗有限元方法对刀盘的可靠性进行分析，考察刀盘结构可靠性，确定刀盘的失效概率及影响其失效的显著因素。 　　关键词：盾构刀盘 可靠性 蒙特卡罗方法 ANSYS 　　中图分类号：U455.39 文献标识码：A 文章编号：1674-098X（2013）01（b）-00-03 　　随着人们对地下空间开发和利用，地下施工设备也在迅猛发展。作为地下施工的主要设备的盾构机正广泛的应用于各类城市的地铁、公路隧道、过江隧道及城市给排水系统的施工中[1]。 　　盾构刀盘作为盾构机的关键功能部件承担开挖土体作用，在整个掘进过程中刀盘工作环境恶劣，受力情况比较复杂，因此具有较高使用寿命刀盘对整个施工具有重要的影响[2]。 　　目前国内周奇才、刘义[3-4]等人对盾构刀盘的结构强度和振动方面已有了较为广泛的研究，得到刀盘了的应力、应变云图及各阶振型。但所其对刀盘定义材料属性和施加的载荷为常量，没有考虑输入参数的不确定性，无法获知在不确定因素影响下，刀盘的失效概率及是否具备良好的工作可靠性。 　　该文采用蒙特卡罗有限元方法对刀盘的可靠性进行分析，考察刀盘结构可靠性，确定刀盘的失效概率及影响其失效的显著因素。 　　1 结构可靠性及蒙特卡罗法 　　1.1 结构可靠性 　　可靠性分析法就是分析所建立模型中的一些输入参数和假设的不确定性对分析结构的影响，并用可靠性指标对结构进行判断，从而在不能完全消除输入参数的不确定性的情况下，提高分析的准确性和设计产品的质量与可靠性。同时，通过可靠性分析可以保证构件安全的基础上，避免不经济的设计[5]。 　　在结构可靠性分析中，极限状态表达式为：。其中随机矢量表征了工程中存在着的不确定信息，如载荷、材料参数的随机性等。 　　该文以刀盘材料的屈服极限、材料属性、载荷为输入变量来计算刀盘结构的可靠性。由强度理论可知，假设刀盘在工作过程中，不允许应力超过屈服强度的事故发生，如果应力超过了材料的屈服强度则认为失效，失效准则为： 　　则极限状态函数为：，为失效状态。其中Z为上式中所有不确定量的向量。 　　在该文中的可靠性就是要求的概率。 　　1.2 蒙特卡罗法 　　目前解决复杂结构可靠性的常用数值解法有蒙特卡罗法。蒙特卡罗法是一种采用数值模拟来解决与随机变量有关的实际工程问题的方法。对随机变量的数值模拟相当于一种实验，所以蒙特卡罗法又称为统计实验法。 　　在ANSYS中蒙特卡罗法分为直接法、拉丁超立方法、自定义方法3种，其中超拉丁方法的效率较直接法高（与直接法相比，避免了重复抽样）[6]。 　　2 可靠性分析文件建立 　　2.1 模型导入与定义参数变量 　　见图1。 　　将在Pro/E中建好的刀盘模型导入ANSYS中，导入后模型如图1。确定用于分析的输入变量和输出变量。随机输入变量有刀盘的载荷（刀盘推力、扭矩）、材料属性（弹性模量、泊松比）、刀盘材料的屈服极限。 　　输出变量有刀盘的最大等效应力、极限状态Z。具体输入变量参数如表1所示。将表1中参数定义在ANSYS中定义为参数变量，如图2所示。 　　2.2 有限元模型建立与加载 　　按表1定义材料单元类型、材料属性、划分网格、定义载荷和约束。采用soild185单元，采用自由划分网格。根据刀盘受力情况，刀盘法兰的连接处施加全 　　约束。 　　刀盘受到阻力视为均匀的分布在刀盘面板的6块辐条上，载荷大小为。在刀盘外径施加切向力。添加约束和载荷的模型如图3所示。 　　2.3 可靠性分析设置 　　进行求解，并提取刀盘的最大等效应力，命令流窗口中输入*get，maxstr，sort，max得到刀盘的最大节点应力，定义极限状态方程，在窗口中输入“Z=yeis-str”，得到分析文件。调入生成的分析文件，将输入变量设置成服从正态分布，并设置相应均值和标准差。确定输出变量为刀盘的最大等效应力和极限状态Z。最后选择拉丁超立方法抽样方法，做100次抽样模拟循环。 　　3 结果分析 　　查看，即刀盘失效的置信度为95%的概率，查取结果如图4所示。从图可以看出在置信度为95%的情况下，Z0的平均的概率为1.1927%，即说明刀盘强度的可靠性为98.8073%。 　　绘制Z的累积分布函数如图5所示。Z的分布区间为(～)，，与图5的结果相对应。直观的看Z的累积分布函数也可以看出Z0的占据了大部