基于UG的“火轮”汽车轮毂三维建模及有限元应力分析

毕业设计(论文) 题 目： 基于UG的“火轮”汽车轮毂 三维建模及有限元应力分析 论文摘要： 车轮是左右整车性能最重要的安全部件，不仅要承受静态时车辆本身垂直方向的自重载荷，更需要承受车辆行驶中来自各个方向因起动、制动、转弯、路面凹凸不平等各种动态载荷所产生的不规则应力之考验。它的轴向跳动和径向跳动精度，又直接影响到整车行驶中的平稳性、抓地性、偏摆性、制动性等行驶性能。随着汽车工业的迅猛的发展，铝合金轮毂的应用越来越普遍。铝合金轮毂具有重量轻，降低油耗；散热性好，提高轮胎寿命；缓冲和吸震性好；造型美观，易加工，耐腐蚀等优点。但中国铝合金车轮行业普遍存在设计周期长，制造成本高等问题。在汽车设计制造中计算机辅助设计是必然趋势，因此采用先进的三维CAD软件和大型CAE软件对汽车铝合金轮毂进行结构设计，有利于缩短设计周期，提高产品质量。目前，在中国用有限单元法对铝合金车轮进行研究还处在起步阶段。只有少数的科研院所和高校对钢制车轮进行有限元分析研究。因此，有必要把有限元技术应用到铝合金车轮上，以解决生产实际问题。 本课题：基于UG的“火轮”汽车轮毂三维建模及有限元应力分析，是从汽车轮毂基本特征入手并结合工业设计美学，打造出具有中国元素的车轮，应用UG软件三维建模并完成网格划分，用有限元分析模块做出强度分析。 通过此次设计学习并熟练掌握UG软件的建模与有限元分析功能，对车轮做强度分析，为车轮结构优化设计提供依据。 关键词：铝合金车轮 UG 有限元 Abstract Wheel is the most important safety component of the vehicle performance, Not only to bear the load of the vehicle itself static weight on vertically. And withstand the test of vehicles from all directions, starting、braking、cornering、uneven surface and other dynamic loads generated by the irregular stress. Its axial and radial runout accuracy the smooth, grip, swing, braking and driving performance. By the rapid development of automobile industry, the increasing application of aluminum alloy wheels. Aluminum alloy wheels with light weight, lower fuel consumption; good , improve tire life; buffering and shock absorption; attractive appearance, ease of processing, corrosion resistance, etc.But in China aluminum alloy wheel industry cycles, is an inevitable trend in automotive design and manufacturing.Therefore the use of 3D CAD software and large CAE software for car aluminum wheel structure design, is benefit to shorten the design cycle and improve the quality of products.At present, in China adopt the finite element method on the aluminum alloy wheels research is still in its initial stage.Therefore, it is necessary to adopt the finite element technique on aluminum alloy wheel, to solve the practical problems of production. This project:Based on UG establish "steamer" alloy wheel 3D modeling and finite element stress analysis, Start from the basic characteristics of the car wheels and combine the aesthetics of industrial design,to create a wheel with Chinese elements.using UG software establish 3D modeling and complete the meshing finite element analysis model with a strength analysis. Through this design study and master UG software 3D modeling and finite element analysis function.Do strength analysis of the wheel.Provide the basis for the optimized design of the wheel structure. Key Words: alloy wheel Unigraphics NX Finite element 目录 第一章 绪论 1 1.1. 课题研究目的和意义 1 1.2. 国内外研究动态及现状分析 2 1.2.1. 车轮的发展趋势 2 1.2.2. 车轮疲劳分析研究 3 1.3. 论文研究目标 4 第二章 理论基础与模型建立 5 2.1 有限元技术及UG软件 5 2.1.1 有限元法基本原理 5 2.1.2 有限元法分析过程 5 2.1.3 UG软件介绍 7 2.2 车轮模型建立 8 2.2.1车轮的结构设计的基本步骤： 8 2.2.2 车轮参数确定及建立模型 9 2.3 本章小结 13 第三章 车轮径向疲劳的有限元分析 14 3.1 车轮径向疲劳试验原理 14 3.2 车轮径向载荷疲劳有限元分析 15 3.2.1 车轮有限元模型建立 15 3.2.2 径向载荷疲劳有限元分析 18 3.3本章小结 22 第四章 车轮弯曲疲劳的有限元分析 23 4.1 车轮弯曲疲劳试验原理 23 4.2 车轮弯曲疲劳有限元分析 24 4.2.1 车轮有限元模型建立 24 4.2.2 车轮弯曲疲劳有限元分析 27 4.3 本章小结 29 第五章 总结与展望 30 5.1 毕业设计总结 30 5.2 未来工作展望 31 参考文献 32 致谢 33 第1章　 绪论 1.1. 课题研究目的和意义 本课题研究从汽车轮毂的基本特征入手，应用UG三维软件对从工业设计美学【1】角度设计出的具有中国元素的汽车车轮建立模型，并用有限元法完成强度分析。通过本课题的研究，熟练掌握并应用UG三维软件，综合应用专业知识能够独立解决实际的工程问题，作为我们从高校迈向社会重要的一步。 随着汽车工业市场全球化，汽车产业已经成为国民经济的一项支柱产业。目前世界汽车市场上欧美、日韩车系独领风骚，很多国家工业起步早，汽车设计制造产业基础牢固，技术力量雄厚，汽车以其优异的性能与独特的设计风格领跑全球汽车市场，从而抑制了其他国家车辆行业的发展，就我国而言，中国的制造业日益成为世界性制造基地。国外各大汽车集团和著名摩托车生产企业，或在中国投资建厂，或在加大采购份额，国内一些技术先进、管理严谨、质量过硬的铝合金车轮生产企业，已经成为他们的全球供货商。显然，这是千载一遇的机遇，但又是新的挑战。我们工业起步晚、自主研发能力差、汽车制造业急于求成，过度依赖国外技术，直接导致我国汽车制造业发展一直很缓慢。就我们多数汽车生产厂的制造技术水平而言，与发达国家的同行比较存在着很大的差距。每个企业都面临着一个提升核心竞争力的问题，掌握先进制造技术是提升核心竞争力的重要方面。过度依赖国外先进、成熟的汽车技术而不是借鉴，更缺乏自主创新能力，这样便不能铸就我国的民族品牌，这就需要我们独立自主的研究本国特色的汽车，争取国产品牌跻身世界汽车产业前列，在车辆工业上我们任重而道远。 车轮是汽车行驶系统最重要的部件之一，车轮的优劣是衡量整车质量和档次的最主要象征之一。车轮设计是汽车外观造型与整体性能的一个重点，车轮不仅要承受静态时车辆本身垂直方向的自重载荷、通风（制动器散热）的需要，各种动态载荷所产生不规则应力的考验，而且车轮造型直接关系着汽车车身设计的品位和档次。因此，如何进行车轮的造型设计，如何设计出有风格特色和审美情趣的车轮已成为设计者们最关心的问题，另外全球汽车数量的迅速增长对能源和环境形成巨大挑战，解决环境和能源问题已经成为世界各个国家的基本共识，实现汽车轻量化即减轻汽车的自重成为首选措施，对铝合金车轮的节能、安全、环保提出了全新的要求，车轮结构轻量化成为设计重点。眼下越来越多的新技术正逐渐被应用到车轮的设计和开发上来。本课题通过对融汇中国传统文化要素的汽车轮毂建模，并用有限元分析技术使其不仅做到结构轻量化同时又满足人们对审美的需求，为汽车轮毂的结构优化设计提供好的 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 【2】。 1.2国内外研究动态及现状分析 1.2.1车轮的发展趋势 最早的铝合金车轮诞生是一些热衷于赛车的爱好者，为了能使车辆更轻以提高赛车的速度，便想方设法对车辆各零部件作 "轻量化"的设计，其中车轮是重点减轻的主要对象。20世纪30年代联邦德国汽车联合会、拜尔(BMW)发动机公司及戴姆勒一奔驰汽车公司，正式将钢制辐条式轮毅与铝质轧制轮辆相结合的车轮装上了汽车，为铝合金车轮的发展奠定了基础。二次世界大战和世界性的能源危机大大刺激了汽车商的"轻量化"需求。1945年汽车厂商纷纷开展批量生产铝合金车轮的研究，主要集中在铝合金车轮材质和成形工艺方面，但由于车轮的特殊安全要求，仍末能实施批量生产。1970年后，拜尔发动机公司首先将铸造铝合金车轮作为特殊部件装到了2002型轿车上，1972年又在双门小轿车上成批装上了铸造铝合金车轮，开创了铸造铝合金车轮批量用于轿车的新局面【3】。 我国的铝合金车轮工业起步较晚，最早使用铝合金车轮是在20世纪80年代初，国营洪都机械厂用砂型铸造生产铝合金车轮，起初用在三轮摩托车上，但数量很少，也末形成气候【8】。到20世纪80年代末，我国出现了第一个具有现代规模的戴卡轮毂制造有限公司，其规模和设备都步入了世界先进行列。紧接着又出现了南海中南铝合金轮载有限公司，这两个生产厂的生产装备部达到国际水平。当时国内汽车和摩托车对铝合金车轮的装车欲望还很低，钢圈还占据着绝对统治市场的地位。随着我国公路设施的飞速发展，铝合金车轮开始以极其迅猛之势在全国得到推广，生产铝合金车轮的工厂也像雨后春笋，在全国迅速蔓延【3】。 汽车的高速化迫使车轮朝 "三化"(扁平化、子午线化、无内胎化)迅猛发展。国外车轮规格趋向大直径、高强度、轻量化、更美观，直径大、轮辋宽使得轮胎与地面接触面积增大从而增大了与地面的附着力，使汽车操纵性更好，安全性又进一步提高。大直径铝合金车轮是现代车辆 "轻量化"、"高速化"、"现代化"的必然产物。 结构上除特殊场合装用二片式和三片式的复合车轮外，整体铸造的铝合金车轮已成为当今世界整个轻合金车轮的主流。锻造铝合金汽车轮毂自2003年在戴卡公司上线以来，目前在我国已崭露头角。以铸锻结合及液锻方式生产汽车轮毂在我国尚属空白。铸造铝轮毂由于起步较早，目前生产工艺已经相当成熟，轮毂所能承受的载荷也大幅提高，由于能源问题．汽车的轻量化对耗油量的影响又不容忽视．车轮作为汽车的重要部件，能否减重也是汽车制造业必须要考虑的。锻造轮毂由于具有很高的比强度，能够在较轻的自重下承载很高的载荷，其安全性更高。此外，锻造轮毂的高强度、高载荷的特性，更适合于重载卡车和豪华大巴车。而将锻造铝合金轮毂应用于重载卡车和豪华大巴在国外也仅仅处于起步推广阶段，在国内来看基本上还是一片空白。虽然铸造轮毂具有很多优势．但毕竟受工艺条件所限，很难在载重量和自身重量上继续突破。而锻造轮毂虽然具有很好的强度和安全性，但是由于受工艺条件所限．其材料利用率较低(与铸造相比)、模具制造费高(多套模具生产一个轮子)、原辅材料投入较多等等．使其生产成本较高。铸锻结合或液体模锻则刚好弥补了铸造与锻造的不足．采用铸锻结合或液体模锻工艺所生产出来的产品，不仅能够实现与铸造产品完全相同的复杂造型，同时对产品的机械性能又有大幅提高。对于工艺复杂的轮型可以采用铸锻结合或液体模锻的方式。铸锻结合和液体模锻在西方少数发达国家已经有所运用，我国尚属空白，因此具有良好地发展空间【4】。 作为象征整车档次之一的车轮外观，在点缀整车的"时装化"作用中越来越向着 "艺术化"方向发展，多变且时肇的铝合金车轮轮辐形态和迷人的色泽成了抢占市场的主要手段之一。传统单调的“辐条式”演变成了 "辐板式"；平面的辐板变成了带空间曲面和弧形面状态，甚至由中心对称演变至不对称艺术图案，以增加模仿制造难度，提高竞争力。对于车轮与整车的匹配和色泽的协调，通过表面处理获取与众不同的效果来竞争币场 (如全涂装亚光色、抛光轮、电镀轮、真空镀膜轮等)。对此，欧洲和亚洲因文化差别和不同的审美观点，形成了千姿百态的世界汽车铝合金车轮市场。 镁和铝是最适宜制造车轮的理想轻合金材料，它们有着传统钢车轮所无法比拟的许多优点，更能适应整车高速化、节能化、现代化、高档次发展的需要。国际上车轮材料用得最为普遍的材料是A356(相当我国的ZLlOlA)合金，也是我国铝合金车轮行业普遍使用的材料。材料方面有向镁合金发展的趋势，镁合金有着极高的比强度、疲劳强度和比弹性模量，是极为理想的车轮制造材料。使用镁合金制造的车轮具有极佳的减振性能。但目前除飞机、跑车等一些特殊场合外，镁合金车轮还末能像铝合金车轮那样广泛用于市场，主要原因是镁合金极易氧化的特点给熔炼、铸造所涉及的环保和生产成本等因素，带来了诸多暂时难于大量生产的难题，而且还有许多产业化的问题有待我们去研究解决。许多学者正研究使镁合金能量产的工艺和设备【5】。 1.2.2车轮疲劳分析研究 车轮的疲劳性能长期以来都是车轮性能研究的重点，20世纪70年代开始，欧美以及日本等国在铝合金车轮的性能试验方面就进行了研究，早在20世纪80年代美国和德国的汽车研究所就开发了FLIPS系统，利用有限元分析预测轮毂疲劳寿命。从此，有限元分析技术在车轮设计、研发中应用越来越广泛，在80到90年代，国外在车轮结构设计优化、疲劳失效等方面的研究全面展开【6】。英国MSC工程软件公司应Goodrich Aerospace公司所求，合作开发了MSC.FATIGUE下的“WHEELS”模块。该模块专用于飞机轮子和其它旋转体的疲劳寿命，可以计算轮上每一节点的应力，并在损伤最严重的表面域为每一节点给出其疲劳寿命和疲劳损伤云图。作为形状和受力情况极为相似的汽车车轮，这一模块同样适用。土耳其学者用有限元仿真方法对钢制冲压车轮的弯曲疲劳试验进行静态仿真分析，得出应力集中处的塑性应变，结合材料的应力-应变曲线采用应力应变法和Miner累积损伤理论测其寿命，并取得很好的结果。美国Akron大学教授用有限元法模拟铝合金车轮的径向疲劳试验，并考虑了轮胎及气压的影响，建立有效地力学模型，求出车轮的应力应变分布【7】【8】。 我国在80年代开始引进有限元分析技术和软件，车轮行业方面有限元研究不多，较为突出的是“汽车车轮结构分析方法与应用”课题，主要研究了车轮的静力学分析、车轮疲劳寿命估算和可靠性分析、结构抗疲劳设计等方面，建立了基于有限元结构分析的车轮优化设计方法及基于车轮台架试验 标准 excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载 的计算机模拟和仿真系统，此外还有很多人对车轮结构优化设计进行研究，通常采用有限元软件，以轮辋和轮辐的厚度为设计变量，以结构轻量化为目标，结果表明优化后结构应力接近于材料的强度极限，材料性能得到充分利用，结构重量有所降低，因此车轮强度分析与结构轻量化设计仍将是研究的重点【9~12】。 1.3论文研究目标 （1）学习和掌握UG软件的3D建模、有限元分析功能。 （2）针对实际工程结构，完成其3D几何建模、网格划分、强度分析。 （3）主要内容是针对对具有“中国风”特色的“火轮”汽车轮毂，用UG软件完成其3D几何建模、网格划分、强度分析。 （4）综合运用所学专业知识，独立完成工程结构件的强度分析,能够独立解决实际的工程结构件优化问题，作为我们从高校迈向工作岗位重要的一步。 第二章 理论基础与模型建立 2.1 有限元技术及UG软件 2.1.1 有限元法基本原理 计算机辅助工程CAE(Computer Aid2ed Engineering) 指工程设计中的分析计算与分析仿真, 而有限元法FEM( FiniteElement Method) 是计算机辅助工程CAE中的一种, 另外CAE还包含了边界元法BEM(Boundary Element Method) 和有限差分法FDM( Finite Difference Method) 等。这几种方法各有其优缺点, 各有其应用领域,但有限元法的应用最广。 有限元法是求解数理方程的一种数值计算方法，是将弹性理论、计算数学和计算机软件有机结合在一起的一种数值分析技术，是解决工程实际问题的一种有力的数值计算工具。有限元是一种离散化的数值方法。离散后的单元与单元间只通过节点相联系, 所有力和位移都通过节点进行计算。对每个单元选取适当的插值函数,使得该函数在子域内部、子域分界面上(内部边界) 以及子域与外界分界面(外部边界) 上都满足一定的条件。然后把所有单元的方程组合起来, 就得到了整个结构的方程。求解该方程,就可以得到结构的近似解。离散化是有限元方法的基础。必须依据结构的实际情况,决定单元的类型、数目、形状、大小以及排列方式。这样做的目的是将结构分割成足够小的单元,使得简单位移模型能足够近似地表示精确解【13】。 因次它可以对各种类型的工程和产品的物理力学性能进行分析、模拟、预测、 评价 LEC评价法下载LEC评价法下载评价量规免费下载学院评价表文档下载学院评价表文档下载 和优化,以实现产品技术创新, 故已广泛应用于各种力学、电学、磁学及很多结合学科领域; 同时, 由于它能够处理耦合问题, 使得其有更大的应用前景。你可以从专业的角度理解有限元:包括变分原理、等效积分和加权余量法等, 也可以从直观的意义上理解有限元: 把连续体划分为足够小的单元, 这些单元通过节点和边连接起来,通过选择简单函数(比如线形函数) 来近似表达位移或应力的分布或变化, 从而得到整个连续体物理量的分布和变化【14】。 2.1.2 有限元法分析过程 所谓有限元法（FEA）基本思想是把连续的几何机构离散成有限个单元，并在每一个单元中设定有限个节点，从而将连续体看作仅在节点处相连接的一组单元的集合体，同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律，再建立用于求解节点未知量的有限元方程组，从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元，选取适当的插值函数，使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时，列出一系列以节点、位移为未知量的线性方程组，利用计算机解出节点位移后，再用弹性力学的有关 公式 小学单位换算公式大全免费下载公式下载行测公式大全下载excel公式下载逻辑回归公式下载 ，计算出各单元的应力、应变，当各单元小到一定程度，那么它就代表连续体各处的真实情况【15】。本课题中主要应用有限元法完成铝合金车轮的强度分析，依据有限元法的基本思想对车轮网格划分，就是将车轮几何模型离散化，分成有限个细小的单元。 具体分析过程可以概括为六个步骤【16】： 1. 结构的离散化 结构的离散化是有限元分析的第一步，关系到计算角精度。离散化的过程简单地说是将分析的结构物划分成有限个单元体，是力学模型变成离散模型，以代替原来的连续体。为了有效地逼近实际的连续体和保证计算精度，就需要考虑选择单元的形状、确定单元的数目和确定划分 方案 气瓶 现场处置方案 .pdf气瓶 现场处置方案 .doc见习基地管理方案.doc关于群访事件的化解方案建筑工地扬尘治理专项方案下载 等问题。然后求解的问题就转变为求有限个自由度的节点位移。有限元法计算精度取决于划分单元的形状、大小、数量和分布情况，通常划分的单元愈多、愈密集、也就愈能反映实际结构状况，计算精度愈高，但计算工作量增大，时间增长，因此必须两方面兼顾，在满足精度的要求下尽可能减少单元数目。 2. 位移模式的选择 结构离散化完成后，就可以对典型单元进行特征分析。为了能用节点位移表示单元体的位移、应力、应变，就必须对单元中位移的分布做出一定的假设，也就是假定位移是坐标的某种简单的函数，这种函数称为唯一模式或位移函数。在有限元法应用中普遍选择多项式作为位移模式。其原因是多项式的数学运算简便，并且从所有光滑函数的局部来看都可以用多项式逼近，即不完全的泰勒级数。 根据所选定的位移模式，就可以导出用节点位移表示单元内任一点位移的关系式，其矩阵式是： (2-1) 式中 为单元内任一点的位移列阵； 为单元的节点位移列阵； 为形函数矩阵，它的元素是位置坐标的函数。 3. 单元力学特性的分析 位移模式选定以后进行单元力学特性的分析，包括三个部分内容。 1　 利用几何方程，由位移表达式(2-1)导出用节点位移表示单元应变的关系式； (2-2) 式中 是单元内任一点的应变列阵； 为单元应变矩阵。 2　 由应变表达式(2-2)导出用节点位移表示单元应力的关系式； (2-3) 式中 为单元内任一点的应力列阵； 是与单元材料有关的弹性矩阵。 3　 由虚功原理建立作用于单元上的节点力和节点位移之间的关系式，单元的刚度方程； (2-4) 式中 为单元刚度矩阵，可以推导出： (2-5) 上式的积分遍及整个单元的体积。 4. 等效节点力的计算 弹性体经过离散化后，假定力是通过节点从一个单元传递到另一个单元，但是作为实际的连续体，力是通过公共边界传递。作用在单元边界上的表面力和作用在单元上的体积力、集中力等都需要等效移置到节点上去，用等效的节点力来替代所有作用在单元上的力。移置的方法是按照作用在单元上的力与等效节点力，在任何虚伪位移上的虚功都相等的原则进行。 5. 单元迭加建立整个结构的平衡方程 集合所有单元的平衡方程，建立整个结构的平衡方程组集总体刚度矩阵 、载荷列阵 以及节点位移列阵 形成的整个结构的平衡方程为： (2-6) 6. 节点位移的求解和单元应力的计算 由结构的平衡方程组(2-6)求出节点位移，再利用公式(2-3)和求出的节点位移来计算各单元的应力，并加以整理得出所求的结果。 2.1.3 UG软件介绍 UG（Unigraphics NX）是Siemens PLM Software公司出品的一个产品工程解决方案，它为用户的产品设计及加工过程提供了数字化造型和验证手段。使企业能够通过新一代数字化产品开发系统实现向产品全生命周期管理转型的目标。 且包含了企业中应用最广泛的集成应用套件，用于产品设计、工程和制造全范围的开发过程。UG针对用户的虚拟产品设计和工艺设计的需求，提供了经过实践验证的解决方案。这是一个交互式CADCAM(计算机辅助设计与计算机辅助制造)系统，它功能强大，可以轻松实现各种复杂实体及造型的建构。 NX 为那些培养创造性和产品技术革新的工业设计和风格提供了强有力的解决方案。利用 NX 建模，工业设计师能够迅速地建立和改进复杂的产品形状， 并且使用先进的渲染和可视化工具来最大限度地满足设计概念的审美要求。 NX 包括了世界上最强大、最广泛的产品设计应用模块。 NX 具有高性能的机械设计和制图功能，为制造设计提供了高性能和灵活性，以满足客户设计任何复杂产品的需要。 NX 优于通用的设计工具，具有专业的管路和线路设计系统、钣金模块、专用塑料件设计模块和其他行业设计所需的专业应用程序。 UG的开发始于1990年7月，它是基于C语言开发实现的。UG NX是一个在二和三维空间无结构网格上使用自适应多重网格方法开发的一个灵活的数值求解偏微分方程的软件工具。其设计思想足够灵活地支持多种离散方案。因此软件可对许多不同的应用再利用。一个给定过程的有效模拟需要来自于应用领域(自然科学或工程)、数学(分析和数值数学)及计算机科学的知识。然而，所有这些技术在复杂应用中的使用并不是太容易。这是因为组合所有这些方法需要巨大的复杂性及交叉学科的知识。最终软件的实现变得越来越复杂，以致于超出了一个人能够管理的范围。一些非常成功的解偏微分方程的技术，特别是自适应网格加密（adaptivemeshrefinement）和多重网格方法在过去的十年中已被数学家研究，同时随着计算机技术的巨大进展，特别是大型并行计算机的开发带来了许多新的可能UG的目标是用最新的数学技术，即自适应局部网格加密、多重网格和并行计算，为复杂应用问题的求解提供一个灵活的可再使用的软件基础【17】。本课题是基于UG软件建模并有限元分析，在建模环境中完成铝合金车轮3D模型，进入高级仿真环境，对车轮划分网格，这是有限元法中结构离散化的应用，将力学模型转化为离散模型，指派材料，施加边界条件，通过求解完成应力分析。 2.2 车轮模型建立 2.2.1车轮的结构设计的基本步骤： 1. 根据车轮的装车情况对设计的车轮进行归类并初定出车轮的各种装配参数。 2. 车轮装配参数确定后，根据车轮的外观形状和偏距、装配情况等要素来确定车轮选用正向轮辋还是反向轮辋及其形状。 3. 轮辋确定后就可以根据装车的要求、车轮形状和载荷等确定出车轮安装和轮辐的初步设计。 4. 根据轮辐形状、轮辋形状和车轮的规格定义出机加余量和掏料结构的设计。 5. 对做出的车轮进行重量计算和有限元分析，确定是否需要进一优化车轮的产品结构。 6. 如果重量计算和有限元分析合理，即可完成车轮的结构设计并输出工程图纸，如果不合理则重复上面步骤进行更改。 2.2.2 车轮参数确定及建立模型 “火轮”的汽车轮毂轮毂材料为A356(ZAlSiMg)，国内牌号为ZLl01A。美国铝业协会标准对A356在不同时期的化学成分有不同的要求，A356合金锭称作A356.2,熔炼过程中的合金为A356.1,最终车轮成品的合金为A356.0，其主要化学成分如表2-1所示[18]： 表2-1 材料A356.0的主要化学成分要求 化学成分 Si Mg Fe Cu Mn Zn Ti 百分含量(%) 材料属性：密度: ；弹性模量: ；泊松比: ；许用应力 。本车轮模型采用整体式正向轮辋，铝合金压铸。 车轮规格：16×7J； 16—名义直径(in)； ×—表示整体式轮辋； 7—表示名义宽度(in)； J—轮辋类型，指深槽轮辋； 根据设计指导书车轮名义直径减去0.8mm，即：405.6x178(mm) 【19】； ET：45mm； PCD：4x100； 安装盘直径：150mm； 设计载荷：1600(LBS)(1LBS=0.454kg)。 根据GBT3487-2005 J型轮辋轮廓应符合如图2.1的规定[19]： 图2.1 5°J型轮辋轮廓 Fig.2.1 5-degree J-type rim contour 这些尺寸给出了轮辋槽的最小尺寸及其位置，以保证轮胎顺利装、拆。较大的槽底深度更方便无内胎气门嘴的安装，5°J型轮辋相关尺寸要求如下表2-1所示【19】： 表2-1 5°J型轮辋主要尺寸要求（单位:mm） 轮辋轮廓 A L(量规) 7J 178.0 19.5 19.5 22.0 45.0 轮辋直径代号 标定直径 平峰周长 圆峰周长 16 405.6 1274.2 1276.4 车轮胎圈座根据装胎方向不同也有不同的要求，如表2-2所示： 表2-2 胎圈座组合 组合状态 内侧 外侧 标准（适用正向轮辋） 圆峰 平峰 选用（适用反向轮辋） 平峰 圆峰 选用 圆峰 圆峰（凸峰圆角R4取3～6mm） 这里采用正向轮辋，即轮辋内侧圆峰、外侧平峰。 依国标及设计指导书轮辋内外侧凸峰的尺寸设计如图2.2所示： 内侧圆峰结构及尺寸 外侧平峰结构尺寸 图2.2 胎圈座设计 Fig.2.2 Bead seat design 根据设计指导书的要求气门孔结构如图2.3： 图2.3 气门孔结构 Fig2.3 Structure of valve 为了凸显艺术外观，螺栓孔与通风孔结合的表现形式，本车轮安装盘设计为裸露螺栓孔，以及根据 ；安装盘直径为 ；为达到更好的装配要求，安装面设计为有0.2°角的锥面，并且不进行涂装与电镀处理，安装盘如图2.5所示： 图2.5 安装盘结构 Fig.2.5 Structure of installation disk 轮辐设计的顺序是先根据车轮的外观设计出轮辐正面的轮廓线，再设计轮辐背面的轮廓线，根据装车空间的需要进行调整两条线的位置，确定轮辐剖截面的设计，考虑轮辐的受力分布、铝液顺序凝固和补缩通道的需要，靠近轮辋部位的轮辐厚度，在20mm到26mm之间。而靠近安装盘部位的轮辐厚度，在25mm到32mm之间。本车轮轮辐靠近轮辋部位较薄，厚度为21mm，靠近安装盘部位初步设计为30mm,强度分析前未进行轮辐掏料。 中心孔结构的设计主要是需要更多的考虑装配的关系，即中心孔与车轴之间、中心孔与卡口盖之间、中心孔与杯盖之间的配合。当车轴与中心孔之间进行配合时，除了保证车轮中心孔与车轴之间的配合为间隙配合外，还要求中心孔有足够的高度，以保证中心孔或装饰盖和车轴顶部不会干涉。本车轮中心孔结构如图2.5所示： 图2.6 中心孔结构 Fig.2.6 Structure of center 1　 施加载荷 由车轮径向疲劳试验原理可知，作用在车轮上的载荷主要有轮胎气压载荷和汽车负载。径向载荷疲劳试验的有限元分析模型中考虑径向载荷作用的同时，必颏考虑充气压力的影响。在一定范围(120°夹角)内轮辋胎圈座上施加呈半余弦弦函数分布的径向载荷，根据该车轮结构的结构特点，取在两个关键部位加载的情况分别进行分析，最大压力加载位置对窗口位置为Load1，也就是本模型通风孔靠近轮辋的小圆角部位，最大压力加载位置对辐条部分设为Load2，，并在整个外轮辋上施加充气压力,载荷情况如图3.4所示[21]： Load1径向载荷(对窗口) Load2径向载荷（对辐条） 气压载荷 图3.4 车轮径向载荷疲劳试验有限元模型 Fig.3.4 FEA model of wheel's radial load of fatigue test 其中径向载荷加载函数计算公式如下[21][22]： (3-2) (3-3) 式中， 为试验要求径向载荷， 分别为内、外轮辋胎圈座宽度， 分别为 内，外轮辋胎圈座半径， 为加载范围的一半，即120°的 ，等于60°,充气压力为实际试验时装胎车轮的轮胎充气压力，本模型气压值为2.5Bar，即250kPa．实验要求径向载荷 可根据公式(3-1): 得出，车轮设计载荷为1600LBS，1600LBS×0.454㎏=727㎏， ，取强化试验系数 =1.6，试验载荷 ，由于轮辋所受的力是二次曲线而不是均布的，车轮120°内轮辋网格分为11部分，所以根据有限元思想在网格单元加载荷可以近似为按二次曲线分布的径向载荷，则径向载荷 按加载角度变化值如下表3.1： 表3.1 轮辋的单元压力值 （单位： ） -50° -40° -30° -20° -10° 0° 10° 20° 30° 40° 50° 0.112 0.210 0.289 0.347 0.383 0.395 0.383 0.347 0.289 0.210 0.112 0.107 0.200 0.276 0.331 0.365 0.377 0.365 0.331 0.276 0.200 0.107 这样完成了车轮有限元模型的建立,下一步进行求解。 3.2.2 径向载荷疲劳有限元分析 进入UG求解模块(Solusion),求解器类型为NX-NASTRAN,两种加载载荷情况下所得变形情况如图3.5所示： Load1变形图 Load2 变形图 图3.5 径向载荷疲劳试验变形 Fig.3.5 Deformation of radial load fatigue test 在两种加载条件下(Load1、Load2)，在正对通风孔小圆角部位的轮辋胎圈座上施加径向载荷时，车轮结构的最大变形量为O．497mm。在正对辐条的轮辋胎圈座上施加径向载荷时，车轮结构的最大变形量为O.547mm；两种载荷条件下的最大变形均发生在车轮承受径向载荷的内轮辋区域，整个车轮结构的轮辋部分都有所变形，内轮辋的变形明显比外轮辋的变形量大，并且对于正对辐条的轮辋胎圈座加载时，轮辋变形量大于正对通风孔加载的变形．总体看来，车轮结构的变形量不大，车轮结构设计的刚度性能完全满足要求。 再来看两种载荷加载情况下应力分布如图3.6所示： Load1应力分布 Load2应力分布 图3.6 应力分布 Fig.3.6 The stress distribution Load1应力集中部分 Load2应力集中部分 图3.7 应力集中 Fig.3.7 The stress concentration 两种加载情况下的应力分析结果分别如上图．Load1给出了载荷正对通风孔小圆角部位时的应力分布情况，分析结果表明最大应力发生在通风孔小圆角处，最大应力 Load2给出载荷正对辐条位置加载时车轮结构内的应力分布情况，分析结果表明最大应力同样发生在通风孔小圆角处，靠近轮辐内侧与轮辋交界的部位，最大应力值为49.74MPa．从应力分析结果看来，在径向载荷作用下车轮结构内的应力水平并不高，最大应力值远低于铝合金车轮材料对称弯曲疲劳极限 ，完全满足强度要求。 3.3本章小结 本章的主要内容是车轮出厂测试的一项重要内容，车轮弯曲径向疲劳的有限元分析，主要阐述了径向疲劳测试的试验原理及方法，然后在UG的FEA(Finite Element Analysis）模块对车轮指派材料属性、网格划分、边界条件加载，来模拟车轮实际的工况，求解出车轮在模拟工况中各个部分的应力分布情况及最大变形量。 第四章 车轮弯曲疲劳的有限元分析 4.1 车轮弯曲疲劳试验原理 车轮弯曲疲劳试验也称动态横向疲劳试验，该试验是使车轮承受一个旋转的弯矩，模拟车轮在行车中承受弯矩负荷，标准要求车轮在试验弯矩下经历一定的疲劳循环后不得出现裂纹等破坏现象。 目前，国内外最常采用的车轮弯曲疲劳试验方法是采取车轮固定，由电机带动偏心块进行旋转，通过加载轴对车轮施加一个旋转的弯矩[24]。其结构原理如图4.1所示；根据车轮弯曲疲劳试验的工作原理，有限元分析的边界条件设置如下：因为车轮内轮辋边缘部分被试验台夹具压紧固定，不能旋转和移动，所以对内轮辋边缘施加全约束，即六个自由度全部被约束．车轮承受的弯矩是通过螺拴使车轮与安装盘、加载轴相连施加的，因此，在车轮结构的弯曲试验分析中需在车轮结构上添加安装盘、加载轴和连接螺栓，以模拟实际的加载方式，在加栽轴末端施加力，该力等于车轮试验弯矩除以建立的加载力臂的长度[25]。 图4.1 弯曲疲劳试验工作原理 Fig.4.1 Principle of Bending fatigue test 车轮试验弯矩 （ ）用 表示，由式（4-1）确定： （4-1） 式中： ——轮胎和道路间的设定摩擦系数； ——静负荷半径，车轮配用最大轮胎负荷下的静半径，单位m; ——车轮的内偏距或外偏距，单位m； ——车轮额定负载，单位N； ——强化系数。 4.2 车轮弯曲疲劳有限元分析 4.2.1 车轮有限元模型建立 1　 车轮几何模型前面已完成，现在需要建立有限元模型，在仿真环境中指派材料属性，划分网格、给定边界条件、施加载荷，完成有限元模型建立。首先进入装配环境，给车轮装配加载轴、螺栓、螺母，完成模型装配。然后进入高级仿真环境，对装配模型进行网格划分如图如图4.2所示： 图4.2车轮网格模型 Fig.4.2 Wheel mesh 2　 边界条件 车轮内轮辋边缘部分被试验台夹具压紧固定，车轮承受的弯矩是通过螺拴使车轮与安装盘、加载轴相连施加的，轮辋不能旋转和移动，所以对内轮辋边缘施加全约束，六个自由度全部被约束，约束如图4.3所示[26]： 图4.3 边界条件 Fig.4.3 Boundary condition 3　 施加载荷 车轮弯曲疲劳试验中螺栓预紧力对分析结果没有太大影响，所以这里不考虑施加螺栓预紧力。取轮胎与道路间摩擦系数为 [27]，如果综合考虑车辆行驶过程中会更加平稳、制动性能好、外视更加美观，根据GBT 2798-2008轮胎可取为65系列[28]，但根据试验要求需要选取最大轮胎，选取75系列轮胎规格为 。最大轮胎静负载半径 ；车轮内偏距 ；车轮额定负荷 ； 取强化试验系数 ，由公式 计算得试验弯 。 试验轴长可以取 ，由 得加载在试验轴末端的力： ， 力的方向这里同样选取两个不同的方向，在XY平面 Load1为力指向X轴方向，即指向轮辐，Load2为力的方向与XY轴成45°角，指向车轮通风孔方向，施加载荷后有限元模型如图4.4所示： Load1(力方向为X轴方向) Load2(力方向与X轴成45°角) 图4.4 车轮弯曲疲劳分析的有限元模型 Fig.4.4 FEA model of the wheel bending fatigue analysis 4.2.2 车轮弯曲疲劳有限元分析 进入UG求解模块(Solusion),求解器类型为NX-NASTRAN,位移变化情况如图4.5所示： Load1 变形图 Load2 变形图 图4.5 变形图 Fig.4.5 Deformation 在两种加载条件下(Load1、Load2)，试验轴末端力的方向指向X轴时(Load1)， 车轮轮辐部位的最大变形量为 。试验轴末端力的方向与XY轴成45°角时(Load2)，轮辐部位圆角处变形量最大，最大变形量为 ；在两种不同方向载荷条件下，轮辐部位的最大变形量有所不同，而且变形关于车轮轴心对称。综合分析、第一种情况下承受弯矩的主要是两根辐条，所以轮辐变形量会大一些。第二种情况四根辐条均有承受弯矩，力相对分散，所以变形量较小。总体看来，车轮结构的刚度性能满足要求。 两种不同方向载荷条件下，车轮应力分布如图4.6所示： Load1 应力分布 Load2 应力分布 图4.6 应力分布 Fig.4.6 The stress distribution 由变形图可以看出车轮在不同方向载荷条件下，应力集中部位都在通风孔靠近轮辋处的小圆角部位，第一种加载方式下应力分散在四根辐条上，小圆角处最大为 ，第二种加载方式下应力较为集中在两根辐条上，小圆角处应力最大为 。 综上所述，通风孔小圆角处为危险区域，而车轮材料A356.0对称弯曲疲劳极限 ，应力集中超过了疲劳极限要求，所以车轮还需要进行优化、对其进行参数化建模，适当调整圆角大小等。 4.3 本章小结 本章的主要内容是车轮出厂测试的另一项重要内容，车轮弯曲疲劳的有限元分析，主要阐述了弯曲疲劳测试的试验原理及方法，在UG装配环境中给车轮装配加载轴、螺栓、螺母，然后进入FEA（Finite Element Analysis）模块对车轮指派材料属性、网格划分、边界条件加载，来模拟车轮实际的工况，求解出车轮在模拟工况中轮辐部分的应力情况及最大变形量。 第五章 总结与展望 5.1 毕业设计总结 本课题基于UG的“火轮”汽车轮毂三维建模及有限元应力分析是从汽车铝合金车轮的基本特征入手，结合工业设计美学设计出具有中国元素的汽车车轮，用UG三维软件建立模型，并用有限元法分析结构，进行优化设计，得到理想模型。通过本课题的研究，熟练掌握并应用UG三维软件，综合应用专业知识能够独立解决实际的工程问题，作为我们从高校迈向社会重要的一步。 毕业设计是我作为一名学生即将完成学业的最后一份综和性作业，不仅是对前面所学知识的一种检验，而且也是对自己能力的一种提高。也是对所学专业学知识的全面总结和综合应用，通过这次毕业设计使我明白了光有我们这几年的理论学习是远远不够，更多的是要实践，要去面对实际的工程问题，提出解决方案，去查阅资料，向他人请教、并不断地优化改进，得出趋近理想的成果，才能使我们的理论知识得到进升华，同时学习新的知识、用知识武装自己，做一个真正意义上的技术人才 当时拿到这个课题，脑海中还没有一个思路，一切从零开始。先使用UG软件从建模开始，这个软件也不熟悉，建模过程中边学习软件边建立模型，严格按照设计指导书的要求，还有参考国家标准的汽车轮辋规格，在强度分析过程中查阅了大量资料，国内目前对这部分的研究也是没有绝对权威的资料，主要涉及到一些高校师生的期刊论文，也有基于实践的工程师指导资料，结合铝合金车轮的实际工况完成了强度分析，从中学到了不少知识，通过这次毕业设计，才清楚了自己要学习的东西还太多，学习是一个长期积累的过程，在以后的工作、生活中都应该不断的学习，努力提高自己知识和综合素质。 此次毕业设计，是我对所学知识理论的检验与总结，能够培养和提高设计者独立分析和解决问题的能力。我不仅把大学四年所学的专业知识融会贯通，又掌握了一门设计软件，同时在查找资料的过程中也了解了许多课外知识，开拓了视野，认识了汽车行业的发展方向，使自己在专业知识方面和动手能力方面有了质的飞跃。也为今后走向社会从事工作铸就了一个良好开端。 在此要感谢我们的指导老师郝老师对我悉心的指导，感谢老师们给我的帮助。在设计过程中，我通过查阅大量有关资料，与同学交流经验和自学，并向师兄、师姐请教等方式，使自己学到了不少知识，也经历了不少艰辛，但收获同样巨大。在整个设计中我懂得了许多东西，也培养了我独立工作的能力，树立了对自己工作能力的信心，相信会对今后的学习工作生活有非常重要的影响。而且大大提高了动手的能力，使我充分体会到了在创造过程中探索的艰难和成功时的喜悦。虽然这个设计做的也不太好，但是在设计过程中所学到的东西是这次毕业设计的最大收获和财富，使我终身受益。 5.2 未来工作展望 本课题着重对具有中国风元素的汽车铝合金车轮进行三维建模，及有限元应力分析，研究内容比较不够深入，只是基于软件建模并分析，缺乏相关实验论证，未免会出现不足。 首先是车轮模型的设计，根据设计指导书和国家标准建立了模型，并未参数化建立模型，没有进行优化设计，也是未来车轮研究方向的重点——车轮轻量化设计。 车轮结构静态有限元分析包括弯曲疲劳试验、径向疲劳试验、13°冲击试验，条件有限没能进行实验研究，用软件模拟实际工况进行分析，车轮实际运行中所受应力是很复杂的，如路面不平整引起的震动，路面石块的冲击，温度对车轮材料的影响，以及汽车急停时螺栓孔部位的变形等，所以在进行径向疲劳分析及弯曲疲劳分析过程中模拟的边界条件有所简化，与实际工况还有偏差是在所难免的。 有限元分析(Finite Element Analysis，FEA)是工程科学的重要工具，有限元技术的应用给企业带来了许多实际效益。然而作为一种近似的计算方法，必热存在误差，如何保证有限元分析的精度，也是有限元分析研究中的一个重点所在。比如本次研究中车轮网格划分疏密，有限元模型与车轮实体的一致程度。参数设定是否合理、软件正确使用与否等都有可能使铝合金车轮有限元分析结果产生偏差，因此，分析中还要注重这些，保证车轮有限元分析模型尽可能的接近实际，而且，采用有效的方法对分析结果的准确程度进行估测，也是十分必要的。 参考文献 [1] 陆斌, 程华, 柳琳. 汽车车轮进化与流行时尚[J]. 现代零部件, 2009,(03) [2] 钟翠霞. 铝合金车轮设计及结构分析[D]. 浙江大学硕士论文，2006 [3] 铝合金车轮制造技术[M] 机械工业出版社. [4] 孟模，张治民，李大旭，李慧娟. 高性能铝合金车轮制造技术[J]. 科学创新与生产力，2011(12) [5] 杜风义.汽车轮毂的发展展望[J] [6] 钟翠霞，童水光，闫胜昝，张响，徐立. 铝合金车轮的有限元分析与疲劳寿命预测[J].机械设计与制造，2006(12) [7] Zhihua Zhu，Jinhua Hu，Xiaoming Yuan，Huixue Sun.，Hongmei Sun. Research on structural optimization of the aluminum alloy wheel[R]. 2012 WASE International Conference on Information Engineering. [8] 王霄锋, 王波. 汽车车轮结构强度分析[J]. 机械强度, 2002, 24(1),66-99. [9] 臧孟炎，秦滔. 铝合金车轮强度有限元分析研究[J]. 科学技术与工程，2009(02). [10] 李世徳. 浅谈铝合金车轮疲劳监测与实验研究[J] [11] 刘刚，胡家琨. 铝合金车轮检验的研究现状[J]. 科技信息，2010(01) [12] 武海滨，闫绍峰，仪登利. 铝合金轮毂的有限元分析[J]. 辽宁工程技术大学学报，2010(04). [13] 李人宪. 有限元法基础[M]. 北京，国防工业出版社 [14] 张斌. 汽车车轮有限元分析与抗疲劳设计[D], 哈尔滨工业大学, 2003. [15] 闫胜昝，童水光，钟翠霞，朱训明，王树山. 基于FEA的车轮结构形状优化设计[J]. 机械设计,2009(01) . [16] 张永昌. MSC.Nastran有限元理论基础与应用[M]. 科学出版社, 2004. [17] 朱爱民，其木格，郑鑫. 基于UG的轮毂结构可靠性分析[J]. 机械工程及自动化，2010(01）. [18] 毛协民 汽车轮毂用A356铝合金特点 [19] GBT 3487-2005 汽车轮辋规格系列[S] [20] 颜伟泽，郝艳华，黄致建，蒋万标. 车轮径向疲劳试验有限元仿真及疲劳寿命估算[J]. 机械设计与制造，2011(06). [21] 韦东来，崔振山. 铝合金车轮径向疲劳试验的数值仿真[J]. 机械强度，2008(06). [22] 闫胜昝，童水光. 铝合金车轮结构设计有限元分析与实验研究[D]. 浙江大学博士论文，2008 [23] GBT 2978-2008 轿车轮胎规格、尺寸、气压与负荷[S]. [24] 汽车铝合金轮毂VIA测试标准(JWL). 自動車技術会. [25] 李平化，罗永新，李强，邓卫军. 铝合金车轮弯曲疲劳性能的有限元分析[J]. 现代机械，2008(03)，54-55. [26] 梁艳，池成忠，张伟，李晋敏. 基于Pro_E铝合金车轮弯曲疲劳试验的研究[J]. 锻压装备与制造技术，2011(02) [27] 周堃，赵玉涛，苏大为，汪煦. 轿车铝合金车轮弯曲疲劳性能的有限元分析[J]. 机械设计与制造，2009（05）. [28] 刘建军. 汽车轮胎的扁平率和速度等级[M].汽车与配件,，1999(44) 致谢 时光飞逝，一学期的毕业设计圆满结束了，在本次毕业设计期间，得到了精心指导和热情鼓励。给予了我学习上的鼓励、课题研究上的指导和生活上的关怀。老师渊博的知识、严谨的治学态度以及对教育事业的饱满热情和献身精神，给我留下了深刻的印象，并将在我今后的工作和学习中不断激励我前进。在此，我向敬爱的郝老师表示我真心的感谢、真诚的敬意和衷心的祝愿! 同时，还要感谢研究员对我的鼓励与帮助，老师治学严谨，为人谦和诚恳，在毕业设计过程中给予了我很大的帮助，在毕业论文上更是提出了许多宝贵的意见。要特别感谢康淑贤师姐，她在整个设计过程中给予我们莫大的关怀与帮助，感谢设计工程系实验室所有的师兄、师姐对我研究过程中悉心的指导，以及提供良好的研究环境。在此，我谨向他们表示我最衷心的感谢! 袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇羅膃蚈螂羁膂莈蚅袇膁蒀袀螃膀薂蚃肂腿节衿羈腿莄蚂袄芈蒇袇螀芇蕿蚀聿芆艿蒃肅芅蒁螈羁芄薃薁袆芃芃螆螂芃莅蕿肁节蒈螅羇莁薀薈袃莀艿螃蝿荿莂薆膈莈薄袁肄莇蚆蚄羀莇莆袀袆羃蒈蚂螂羂薁袈肀肁芀蚁羆肁莃袆袂肀薅虿袈聿蚇蒂膇肈莇螇肃肇葿薀罿肆薂螆袅肅芁薈螁膅莃螄聿膄蒆薇袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈 芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈葿螈聿蒄葿袁羁莀蒈羃膇芆蒇蚃羀膂蒆螅膅蒁薅袇羈莇薄罿膄芃薃虿羆艿薃袁节膅薂羄肅蒃薁蚃芀荿薀螆肃芅蕿袈芈膁蚈羀肁蒀蚇蚀袄莆蚇螂肀莂蚆羅袂芈蚅蚄膈膄蚄螇羁蒂蚃衿膆莈蚂羁罿芄螁蚁膄膀螁螃羇葿螀袅膃蒅蝿肈羆莁螈螇芁芇莄袀肄膃莄羂艿蒂莃蚂肂莈蒂螄芈芄蒁袆肀膀蒀罿袃薈 _1234567921.unknown _1234567937.unknown _1234567953.unknown _1234567961.unknown _1234567965.unknown _1234567967.unknown _1234567969.unknown _1234567970.unknown _1234567971.unknown _1234567968.unknown _1234567966.unknown _1234567963.unknown _1234567964.unknown _1234567962.unknown _1234567957.unknown _1234567959.unknown _1234567960.unknown _1234567958.unknown _1234567955.unknown _1234567956.unknown _1234567954.unknown _1234567945.unknown _1234567949.unknown _1234567951.unknown _1234567952.unknown _1234567950.unknown _1234567947.unknown _1234567948.unknown _1234567946.unknown _1234567941.unknown _1234567943.unknown _1234567944.unknown _1234567942.unknown _1234567939.unknown _1234567940.unknown _1234567938.unknown _1234567929.unknown _1234567933.unknown _1234567935.unknown _1234567936.unknown _1234567934.unknown _1234567931.unknown _1234567932.unknown _1234567930.unknown _1234567925.unknown _1234567927.unknown _1234567928.unknown _1234567926.unknown _1234567923.unknown _1234567924.unknown _1234567922.unknown _1234567905.unknown _1234567913.unknown _1234567917.unknown _1234567919.unknown _1234567920.unknown _1234567918.unknown _1234567915.unknown _1234567916.unknown _1234567914.unknown _1234567909.unknown _1234567911.unknown _1234567912.unknown _1234567910.unknown _1234567907.unknown _1234567908.unknown _1234567906.unknown _1234567897.unknown _1234567901.unknown _1234567903.unknown _1234567904.unknown _1234567902.unknown _1234567899.unknown _1234567900.unknown _1234567898.unknown _1234567893.unknown _1234567895.unknown _1234567896.unknown _1234567894.unknown _1234567891.unknown _1234567892.unknown _1234567890.unknown