中间包结构有限元分析
摘要介绍了某钢厂中间包结构存在的问题，简要论述了中间包产生变形和裂纹的机理。利用数值模拟的方法对中间包结构强度和刚度进行有限元分析，通过计算所得的中间包温度场和应力场，显示中间包结构高应力区和强度的薄弱位置，提出改进方案。此外还改变中间包的耐火材料层的厚度和综合导热系数，分析这些因素对中间包温度场和应力场的影响，为中间包结构的优化提供理论支持。
关键词中间包结构强度刚度有限元分析
Finite Element Analysis of the Tundish Structure
NI Sai-zhen, LI Fu-shuai, TAO Jin-ming
(Metallurgical technology research institute of Beijing in CCTEC , Beijing 100028, China) Abstract In this paper, we introduced the problems of the tundish structure in a steel work at first, and briefly discussed the mechanism for the occurring of the tundish deformation and crack. Using the finite element method to analyze strength and stiffness for the tundish structure, according to the temperature and stress field, we can find hot point and high stress region. The effects of the fireproof material thickness and the total conductive coefficient on the temperature and stress field are also studied.
Key words tundish structure strength stiffness finite element analysis
1 前言
一直以来对中间包的研究都侧重于中间包内流场的分析，有关中间包流场研究的文章很多[1-5]，而对于中间包包体本身结构分析，研究者却很少关注，在这方面的文献也较少。中间包的强度以及结构的稳定性对于浇注的顺利进行以及保证铸坯质量方面同样起到很重要的作用。在热应力以及外载荷作用下中间包会产生变形，改变水口间的相对位置。如果变形过大的话，会影响到水口对中操作。包体的变形还可能使其产生裂纹，严重破坏包体结构，从而发生事故，不利于安全生产的进行。
某钢厂中间包为七机七流，铸机断面尺寸为150mm×150mm，流间距为1250mm，浇注周期约为36min，主要生产碳素结构钢Q235B，优质碳素结构钢45#，低合金结构钢
25MnSiV、Q345B等钢种，该中间包为T形结构，容量为40吨。中间包内衬耐火材料由外向内依次为工作层、永久层、保温层。该中间包存在以下问题：
（1）现场反应变形比较严重，而相应结构的六机六流的中间包变形问题不明显；
（2）新的中间包在开始浇铸时，靠四个耳轴支撑，中间底部与中间罐车横梁不接触，但随着浇铸时间的不断增加，中间就会慢慢凹陷，浇铸大约5-6小时后，中间
底部就会与横梁接触；
（3）旧中间包或多或少都存在中间凹陷的永久变形，有些变形较大，在浇铸前中间
底部就已经与横梁接触；
（4） 出现过包壳发红的现象，这说明现场出现包壳温度过高的状况，对浇注不利。 文献[6]中提到中间包耳轴产生裂纹的力学原因是由于局部温度梯度较大，产生热应力集中。中间包在热应力的作用下发生了扭曲变形，使得水口相对位置发生改变。文献[7]同样认为中间包产生裂纹的主要原因是由热应力引起的，重力的影响较小。文献[8]指出中间包起吊和工作时，最大应力分别出现在吊耳和平板上， 其中需要重点关注起吊状态下的吊耳处。
2 有限元模型
2.1 几何模型的建立
利用商业软件ANSYS 根据几何平面图建立三维立体几何模型，考虑到中间包几何模型的对称性，建立1/2实体模型图。
2.2 有限元模型的建立
中间包的耐火材料分成3层，紧靠钢结构的一层约为10mm ，为保温层，导热系数为0.37 /W m K ?； 接下来是永久层，侧边厚度为135mm ，底面厚度为175mm ，导热系数为1.4 /W m K ?；靠近钢水的一层为工作层，厚度为35mm ，导热系数为0.7 /W m K ?。为了模型计算方便，将它们综合为一层，侧边厚度为180mm ，底面厚度为220mm ，综合导热系数为
1.2 /W m K ?。
中间包大部分材料为Q235，采用SOLID70单元对模型进行网格划分，划分后的总单元数为80125，节点数21535，划分后的网格图如图1所示。
图1 中间包单元网格图
根据中间包的实际工作状况，本文采用热-结构分步耦合计算，即先根据热边界条件计算出温度场，然后将其结果作为热载荷施加到结构计算中去。
2.2.1 温度场计算载荷与边界条件
根据中间包的载荷描述，中间包工作液面800mm ，溢流高度为900mm ，最大承载钢水42吨。因此模型内衬面与钢水接触位置节点施加固定温度载荷，本模型中选择高温1560℃。中间包外表面施加对流换热系数为2
10/W m K ，环境温度30℃。中间包初始温度30℃。计算浇铸24小时后中间包温度场的分布。
2.2.2 结构计算载荷与边界条件
因为中间包是对称模型，在对称面上施加对称位移边界条件；考虑到内衬在铺设时存在间隙，因此在结构计算时去掉内衬模型；选择结构内腔表面0-900mm 的面单元，根据不同高度施加铁水静压力，即
P=ρgh
其中h 为单元中心距离底面高度，g 为重力加速度，ρ为铁水密度，取37000/kg m 。 分别在耳轴位置的轴面上施加Z 向位移约束，并约束其中一个耳轴的横向，即Y 向位移。其中Z 向竖直向上，Y 向为中间包宽度方向，X 向为中间包长度方向，如图[2]所示：
图2 结构计算载荷与边界条件示意图
3 有限元分析结果
3.1 温度场分布
图3 整个中间包温度场的分布图4 中间罐钢结构温度场的分布如图3 所示，为经过24小时后整个中间包温度场的分布，可以看出，由于耐火材料具有良好的隔热效果，热扩散速度比较慢，大的温度梯度主要在耐火材料层内部。
图4为经过24小时后中间包钢结构温度场的分布，可以看出，最高温度位于靠近耐火材料一侧，中间包的侧面，拐角位置，其最高温度为461.1℃，其它位置温度如云图所示。需要指出的是，该温度是在24小时后的温度场分布，随着浇铸时间的延长，钢结构的温度
会逐渐升高，最后会达到一个热平衡温度场，此时中间包钢结构的最高温度为597℃。
3.2 浇注状态下应力分布
浇铸状态下，承载耳轴为左右两侧，即图2中所示位置约束。中间包钢结构的应力场分布云图如图5所示。中间包最大应力为478MPa，位于最高温度对应位置，大于Q235的屈服极限。图6为大于250 MPa应力的单元显示，可以看出大于材料屈服极限的区域较大面积的出现，主要位于溢流高度的钢结构内面位置，分析其原因主要是该位置温度高，温度梯度大。因此可以认为如果中间包一次性长期工作达到24小时以上，中间包钢结构本体可能会发生塑性变形，即便冷却后也不能完全恢复。
图5 中间包钢结构的应力场分布图6 中间钢结构大应力单元显示图
3.3 净载荷状态下应力分布
如果不考虑热应力的影响，只考虑重力载荷，此时中间包最大应力为110MPa，位于顶板拐角处，小于Q235的屈服极限。由于顶板不构成中间包整体结构强度的破坏，因此如果去掉顶板，中间包最大应力为64.2MPa，位于耳轴的根部，远小于Q235的屈服极限。因此可以说中间包净载荷强度满足够，满足工程要求。
3.4 位移场分布
中间包钢结构Z向位移场分布如图7所示，可以看出，由于钢水对其向下的压力作用，尤其是钢结构受热后引起的温度不均匀分布，中间包钢结构中部会向下凹沉，对于浇铸状态，最大向下位移为31.35 mm。
中间包钢结构X向位移场分布如图8所示，可以看出，由于钢结构随着时间，温度逐渐升高，热膨胀的作用使得中间包钢结构沿X向延展，边部最大延展为20.0mm。而对于水口位置，各水口（从中心到外侧的顺序）与冷态相比，分别向外延展1mm，7.45mm，13.9mm，20mm。设计者可以将此计算作为参考，合理布置水口位置，以便使其浇铸过程中偏流情况最小。
图7中间包钢结构Z 向位移场分布 图8 中间包钢结构X 向位移场分布
3.5 中间包耐火材料增厚温度场及应力场分析
考虑到钢结构在长时间使用时出现的热变形问题，假设将长时间使用的中间罐耐火材料增加厚度，侧面厚度为250mm ，底面厚度为270mm ，综合导热系数为1.2 W.m/k ， 计算其温度场及应力场分布情况。
浇铸状态下，经过24小时后，中间包钢结构的最高温度仍然位于靠近耐火材料一侧，中间包的侧面，拐角位置，最高温度为305.4℃，比现有结构温度降低了155℃。 经过无限长时间后， 最高温度为473.5℃，也低于现有结构的597.4℃。
连续工作24小时后中间包最大应力为291MPa ，远小于现有结构的478MPa ，位于最高温度对应位置，虽然也大于Q235的屈服极限，但是通过大250MPa 应力的单元显示，发现：大于材料屈服极限的区域很小，只是以点状零星分布在溢流高度位置，及个别拐点位置。因此可以认为耐火材料增厚的中间包一次性长期工作达到24小时以上，中间包钢结构本体虽然也可能会发生局部塑性变形，但变形量较小。
耐材增厚后中间包钢结构中部向下凹沉20.3mm ，小于现有结构的31.4 mm 。 对于相应的X 向位移场分布：由于钢结构随着时间，温度逐渐升高，热膨胀的作用使得中间包钢结构沿X 向延展，边部最大延展为13.6mm 。而对于水口位置，各水口（从中心到外侧的顺序）与冷态相比，分别向外延展0.64mm ，4.9mm ，9.01mm ，12.8mm 。
3.6 中间罐耐火材料增厚并改变导热系数温度场及应力场分析
如果在上述增厚的基础上，改变耐火材料的导热性能，取综合导热系数为1/W m K ， 计算其温度场及应力场分布情况。
浇铸状态下，经过24小时后，中间包钢结构的最高温度也是位于靠近耐火材料一侧，中间罐的侧面，拐角位置，最高温度为253.8 ℃，比现有结构温度降低了207℃。 经过无限长时间后，最高温度为416.7℃，也低于现有结构的597.4℃。
连续工作24小时中间包最大应力为248MPa ，远小于现有结构的 478MPa ，位于最高温度对应位置，虽然也大于Q235的屈服极限，但是通过大于250 MPa 应力的单元显示，发现：大于材料屈服极限的区域很小，存在于个别拐点位置，并少于导热系数为1.2 /W m K 的结构。
改变综合导热系数后中间包钢结构中部向下凹沉18.1mm ，小于现有结构的31.4 mm 。 对于相应的X 向位移场分布：由于钢结构随着时间，温度逐渐升高，热膨胀的作用使得中间罐车钢结构沿X 向延展，边部最大延展为11.45mm 。而对于水口位置，各水口（从中心到外侧的顺序）与冷态相比，分别向外延展0.54mm ，4.12mm ，7.43mm ，10.51mm 。
通过以上数据可以发现，采用阻热效果好的保温材料，能有效地减少热量向中间包钢结构传递，从而减小钢结构的热应力，保证钢结构的强度和刚度。由于本文中的中间包变形过大，建议增加中间包耐火材料的厚度或者采用保温效果更好的材料，另外在中间包中间底部加支撑结构也能起到减小变形的效果。
4 结语
通过对某钢厂中间包存在问题的分析，发现热应力集中是使中间包结构产生变形和裂纹的主要原因。对现有的结构进行有限元分析，得到中间包结构最大应力超过材料的屈服极限，会产生结构的永久变形，而且结构的最高温度较高，对中间包的安全性产生重要影响。通过增加耐火材料的厚度和减小综合导热系数进行进一步计算，可以看出能有效减少热量向中间包钢结构传递，降低结构的最高温度和最大变形，减少集中应力的产生。
文章来源：中国冶金装备网出师表
两汉：诸葛亮
先帝创业未半而中道崩殂，今天下三分，益州疲弊，此诚危急存亡之秋也。然侍卫之臣不懈于内，忠志之士忘身于外者，盖追先帝之殊遇，欲报之于陛下也。诚宜开张圣听，以光先帝遗德，恢弘志士之气，不宜妄自菲薄，引喻失义，以塞忠谏之路也。
宫中府中，俱为一体；陟罚臧否，不宜异同。若有作奸犯科及为忠善者，宜付有司论其刑赏，以昭陛下平明之理；不宜偏私，使内外异法也。
侍中、侍郎郭攸之、费祎、董允等，此皆良实，志虑忠纯，是以先帝简拔以遗陛下：愚以为宫中之事，事无大小，悉以咨之，然后施行，必能裨补阙漏，有所广益。
将军向宠，性行淑均，晓畅军事，试用于昔日，先帝称之曰“能”，是以众议举宠为督：愚以为营中之事，悉以咨之，必能使行阵和睦，优劣得所。
亲贤臣，远小人，此先汉所以兴隆也；亲小人，远贤臣，此后汉所以倾颓也。先帝在时，每与臣论此事，未尝不叹息痛恨于桓、灵也。侍中、尚书、长史、参军，此悉贞良死节之臣，愿陛下亲之、信之，则汉室之隆，可计日而待也。
臣本布衣，躬耕于南阳，苟全性命于乱世，不求闻达于诸侯。先帝不以臣卑鄙，猥自枉屈，三顾臣于草庐之中，咨臣以当世之事，由是感激，遂许先帝以驱驰。后值倾覆，受任于败军之际，奉命于危难之间，尔来二十有一年矣。
先帝知臣谨慎，故临崩寄臣以大事也。受命以来，夙夜忧叹，恐托付不效，以伤先帝之明；故五月渡泸，深入不毛。今南方已定，兵甲已足，当奖率三军，北定中原，庶竭驽钝，攘除奸凶，兴复汉室，还于旧都。此臣所以报先帝而忠陛下之职分也。至于斟酌损益，进尽忠言，则攸之、祎、允之任也。
愿陛下托臣以讨贼兴复之效，不效，则治臣之罪，以告先帝之灵。若无兴德之言，则责攸之、祎、允等之慢，以彰其咎；陛下亦宜自谋，以咨诹善道，察纳雅言，深追先帝遗诏。臣不胜受恩感激。
今当远离，临表涕零，不知所言。
第1次作业（结构力学二） 一、单项选择题（本大题共40分，共 20 小题，每小题 2 分） 1. 位移法的基本结构是( ) A. 静定刚架； B. 单跨静定梁的组合体； C. 单跨超静定梁的组合体 D. 铰结体系 2. ：以下关于影响线的说法不正确的一项为( ) A. 影响线指的是单位力在结构上移动时所引起的结构的某一内力(或反力)变化规律的图形 B. 利用影响线可以求结构在固定荷载作用下某个截面的内力 C. 利用影响线可以求结构某个截面内力的最不利荷载位置 D. 影响线的横坐标是截面位置，纵坐标为此截面位置处的截面内力值 3. A. B. C. D. 仅由平衡条件不能确定 4. 不计杆的分布质量，图示体系的动力自由度为( ) A. 1； B. 2； C. 3； D. 4 5. 用力法计算超静定结构时，其基本未知量为 A. 杆端弯矩； B. 结构角位移； C. 结点线位移； D. 多余未知力 6. 单元坐标转换矩阵是（） A. 奇异矩阵 B. 对称三对角矩阵 C. 对称非奇异矩阵 D. 正交矩阵 7. 位移法的基本未知量包括（） A. 独立的角位移 B. 独立的线位移 C. 独立未知的结点角位移和线位移 D. 结点位移 8. 图乘法计算位移的公式中( ) A. A和y C 可取自任何图形B. A和y C 必须取自直线图形 C. 仅要求A必须取自直线图形 D. 仅要求y C 必须取自直线图形 9. 已知材料屈服极限 =300MPa,结构截面形状如图所示，则极限弯矩Mu=（）
10. 整体坐标系下单元刚度矩阵与下面的哪一个因素无关 A. 局部坐标与整体坐标的选取 B. 结构的约束信息 C. 单元的几何参数 D. 杆端位移与杆端力之间的变换关系 11. 欲减小图示结构的自振频率，可采取的措施有（） A. 减小质量m B. 增大刚度EI C. 将B支座改为固定端 D. 去掉B支座 12. 图（b）为图（a）所示结构MK影响线，利用该影响线求得图（a）所示固定荷载作用下的MK值为（） A. 4kN?m B. 2kN?m C. -2kN?m D. -4kN?m 13. 图示为三自由度体系的振型，其相应的频率是ω a 、ω b 、ω c ，它们之间的大小关系应是( ) A. B. C. D. 14. 图（a）所示一组移动荷载作用在图（b）所示的梁上，则C截面弯矩的最不利位置为（） A. P 1作用在C点上 B. P 2 作用在C点上 C. P 3 作用在C点上 D. P 3 作用在B点上 15. 平面杆件自由单元(一般单元)的单元刚(劲)度矩阵是( ) A. 非对称、奇异矩阵 B. 对称、奇异矩阵 C. 对称、非奇异矩阵 D. 非对称、非奇异矩阵 16. 对称结构在反对称荷载作用下，内力图中为正对称的是( ) A. 弯矩图 B. 剪力图 C. 轴力图 D. 弯矩图、剪力图和轴力图 17. 由于温度改变，静定结构（） A. 会产生内力，也会产生位移； B. 不产生内力，会产生位
有限元法的基本思想是将连续的结构离散成有限个单元，并在每一个单元中设定有限个节点，将连续体看做是只在节点处相连接的一组单元的集合体；同时选定场函数的节点值作为基本未知量，并在每一个单元中假设一个近似插值函数表示单元中场函数的分布规律；然后利用力学中的变分原理建立求解节点未知量的有限元方程，这样就将一个连续域中的无限自由度的问题转化为离散域的自由度问题。求解后可以利用已知的节点值和插值函数确定单元以及整个集合体上场函数。 ANSYS结构有限元分析流程 1．前处理 前处理的目的是建立一个符合实际情况的结构有限元模型。在Preprocessor 处理器中进行。包括：分析环境设置（指定分析工作名称、分析标题）、定义单元类型、定义实常数、定义材料属性（如线弹性材料的弹性模量、泊松比、密度）、建立几何模型（一般用自底向上建模：先定义关键点，由这些点连成线，由线组成面，再由线形成体）、对几何模型进行网格划分（分为三个步骤：赋予单元属性、指定网格划分密度、网格划分） 2．施加载荷、设置求解选项并求解 这些工作通过SOLUTION 处理器实现。 指定分析类型（静力分析、模态分析、谐响应分析、瞬态动力分析、谱分析等）、设置分析选项（不同分析类型设置不同选项，有非线性选项设置、线性设置和求解器设置）、设置载荷步选项(包括时间、
子步数、载荷步、平衡迭代次数和输出控制)、加载（ANSYS结构分析的载荷包括位移约束、集中力、面载荷、体载荷、惯性力、耦合场载荷，将其施加于几何模型的关键点、线、面、体上）然后求解。3．后处理 当完成计算以后，通过后处理模块查看结果。ANSYS软件的后处理模块包括通用后处理模块（POST1）和时间历程后处理模块(POST26)。可以轻松获得求解计算结果，包括位移、温度、应变、热流等，还可以对结果进行数学运算，然后以图形或者数据列表的形式输出。结构的变形图、内力图（轴力图、弯矩图、剪力图），各节点的位移、应力、应变，还有位移应力应变云图都可以得出，为我们分析问题提供重要依据。 ANSYS软件提供了100种以上的单元类型，用来模拟工程中的各种材料和结构，各种不同单元组合在一起，成为具体物理问题的抽象模型。如在隧道工程中衬砌用beam3梁单元模拟，弹簧单元COMBIN14模拟围岩与结构的相互作用。边坡工程中边坡土体用平面单元来模拟。水利工程中对大坝进行三维模拟分析时用实体单元，二维分析时用平面单元；水库闸门用壳单元模拟。桥梁结构模拟分析中，用梁单元模拟不同截面的钢梁、混凝土梁，壳单元模拟桥面板箱梁等薄壁结构，杆单元可以模拟预应力钢筋和桁架。房屋建筑结构中，梁单元模拟框架柱，壳单元模拟屋面板，实体单元模拟大体积混凝土，杆单元模拟预应力钢筋等。 一般都要对结构进行静力分析，结果必须满足设计要求。当动荷
1、结构的动力特性一般指什么？ 答：结构的动力特性是指：频率（周期）、振型和阻尼。动力特性是结构固有的，这是因为它们是由体系的基本参数（质量、刚度）所确定的、表征结构动力响应特性的量。动力特性不同，在振动中的响应特点亦不同。 2、什么是阻尼、阻尼力，产生阻尼的原因一般有哪些？什么是等效粘滞阻尼？ 答：振动过程的能量耗散称为阻尼。 产生阻尼的原因主要有：材料的内摩擦、构件间接触面的摩擦、介质的阻力等等。当然，也包括结构中安装的各种阻尼器、耗能器。 阻尼力是根据所假设的阻尼理论作用于质量上用于代替能量耗散的一种假想力。粘滞阻尼理论假定阻尼力与质量的速度成比例。 粘滞阻尼理论的优点是便于求解，但其缺点是与往往实际不符，为扬长避短，按能量等效原则将实际的阻尼耗能换算成粘滞阻尼理论的相关参数，这种阻尼假设称为等效粘滞阻尼。 3、采用集中质量法、广义位移法（坐标法）和有限元法都可使无限自由度体系简化为有限自由度体系，它们采用的手法有何不同？ 答：集中质量法：将结构的分布质量按一定规则集中到结构的某个或某些位置上，认为其他地方没有质量。质量集中后，结构杆件仍具有可变形性质，称为“无重杆”。 广义坐标法：在数学中常采用级数展开法求解微分方程，在结构动力分析中，也可采用相同的方法求解，这就是广义坐标法的理论依据。所假设的形状曲线数目代表在这个理想化形式中所考虑的自由度个数。考虑了质点间均匀分布质量的影响（形状函数），一般来说，对于一个给定自由度数目的动力分析，用理想化的形状函数法比用集中质量法更为精确。 有限元法：有限元法可以看成是广义坐标法的一种特殊的应用。一般的广义坐标中，广义坐标是形函数的幅值，有时没有明确的物理意义，并且在广义坐标中，形状函数是针对整个结构定义的。而有限元法则采用具有明确物理意义的参数作为广义坐标，且形函数是定义在分片区域的。在有限元分析中，形函数被称为插值函数。 综上所述，有限元法综合了集中质量法和广义坐标法的特点： (l) 与广义坐标法相似，有限元法采用了形函数的概念。但不同于广义坐标法在整体结构上插值（即定义形函数），而是采用了分片的插值，因此形函数的表达式（形状）可以相对简单。 (2) 与集中质量法相比，有限元法中的广义坐标也采用了真实的物理量，具有直接、直观的优点，这与集中质量法相同。 4、直接动力平衡法中常用的有哪些具体方法？它们所建立的方程各代表什么条件？ 答：常用方法有两种：刚度法和柔度法。刚度法方程代表的是体系在满足变形协调条件下所应满足的动平衡条件；而柔度法方程则代表体系在满足动平衡条件下所应满足的变形协调条件。 5、刚度法与柔度法所建立的体系运动方程间有何联系？各在什么情况下使用方便？ 答：刚度法与柔度法建立的运动方程在所反映的各量值之间的关系上是完全一致的。由于刚度矩阵与柔度矩阵互逆，刚度法建立的运动方程可转化为柔度法建立的方程。一般来，对于单自由度体系，求[δ]和求[k]的难易程度是相同的，因为它们互为倒数，都可以用同一方法求得，不同的是一个已知力求位移，一个已知位移求力。对于多自由度体系，若是静定结构，一般情况下求柔度系数容易些，但对于超静定结构就要根据具体情况而定。若仅从建立运动方程来看，当刚度系数容易求时用刚度法，柔度系数容易求时用柔度法。 6、计重力与不计重力所得到的运动方程是一样的吗？ 答：如果计与不计重力时都相对于无位移的位置来建立运动方程，则两者是不一样的。但如果计重力时相对静力平衡位置来建立运动方程，不计重力仍相对于无位移位置来建立，
有限元分析软件 有限元分析是对于结构力学分析迅速发展起来的一种现代计算方法。它是50 年代首先在连续体力学领域--飞机结构静、动态特性分析中应用的一种有效的数值分析方法，随后很快广泛的应用于求解热传导、电磁场、流体力学等连续性问题。 有限元分析软件目前最流行的有：ANSYS、ADINA、ABAQUS、MSC 四个比较知名比较大的公司，其中ADINA、ABAQUS 在非线性分析方面有较强的能力目前是业内最认可的两款有限元分析软件，ANSYS、MSC 进入中国比较早所以在国内知名度高应用广泛。目前在多物理场耦合方面几大公司都可以做到结构、流体、热的耦合分析，但是除ADINA 以外其它三个必须与别的软件搭配进行迭代分析，唯一能做到真正流固耦合的软件只有ADINA。ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 结构分析能力排名：ABAQUS、ADINA、MSC、ANSYS 流体分析能力排名：ANSYS、ADINA、MSC、ABAQUS 耦合分析能力排名：ADINA、ANSYS、MSC、ABAQUS 性价比排名：最好的是ADINA，其次ABAQUS、再次ANSYS、最后MSC ABAQUS 软件与ANSYS 软件的对比分析： 1．在世界范围内的知名度：两种软件同为国际知名的有限元分析软件，在世界范围内具有各自广泛的用户群。ANSYS 软件在致力于线性分析的用户中具有很好的声誉，它在计算机资源的利用，用户界面开发等方面也做出了较大的贡献。ABAQUS软件则致力于更复杂和深入的工程问题，其强大的非线性分析功能在设计和研究的高端用户群中得到了广泛的认可。由于ANSYS 产品进入中国市场早于ABAQUS，并且在五年前ANSYS 的界面是当时最好的界面之一，所以在中国，ANSYS 软件在用户数量和市场推广度方面要高于ABAQUS。但随着ABAQUS北京办事处的成立，ABAQUS软件的用户数目和市场占有率正在大幅度和稳步提高，并可望在今后的几年内赶上和超过ANSYS。 2．应用领域：ANSYS 软件注重应用领域的拓展，目前已覆盖流体、电磁场和多物理场耦合等十分广泛的研究领域。ABAQUS 则集中于结构力学和相关领域研究，致力于解决该领域的深层次实际问题。 3．性价比：ANSYS 软件由于价格政策灵活，具有多种销售方案，在解决常规的
有限元分析软件的比较（购买必看）－转贴 随着现代科学技术的发展，人们正在不断建造更为快速的交通工具、更大规模的建筑物、更大跨度的桥梁、更大功率的发电机组和更为精密的机械设备。这一切都要求工程师在设计阶段就能精确地预测出产品和工程的技术性能，需要对结构的静、动力强度以及温度场、流场、电磁场和渗流等技术参数进行分析计算。例如分析计算高层建筑和大跨度桥梁在地震时所受到的影响，看看是否会发生破坏性事故；分析计算核反应堆的温度场，确定传热和冷却系统是否合理；分析涡轮机叶片内的流体动力学参数，以提高其运转效率。这些都可归结为求解物理问题的控制偏微分方程式，这些问题的解析计算往往是不现实的。近年来在计算机技术和数值分析方法支持下发展起来的有限元分析（FEA，Finite Element A nalysis）方法则为解决这些复杂的工程分析计算问题提供了有效的途径。在工程实践中，有限元分析软件与CAD系统的集成应用使设计水平发生了质的飞跃，主要表现在以下几个方面： 增加设计功能，减少设计成本； 缩短设计和分析的循环周期； 增加产品和工程的可靠性； 采用优化设计，降低材料的消耗或成本； 在产品制造或工程施工前预先发现潜在的问题； 模拟各种试验方案，减少试验时间和经费； 进行机械事故分析，查找事故原因。 在大力推广CAD技术的今天，从自行车到航天飞机，所有的设计制造都离不开有限元分析计算，FEA在工程设计和分析中将得到越来越广泛的重视。国际上早20世纪在50年代末、60年代初就投入大量的人力和物力开发具有强大功能的有限元分析程序。其中最为著名的是由美国国家宇航局（NASA）在1965年委托美国计算科学公司和贝尔航空系统公司开发的NASTRAN有限元分析系统。该系统发展至今已有几十个版本，是目前世界上规模最大、功能最强的有限元分析系统。从那时到现在，世界各地的研究机构和大学也发展了一批规模较小但使用灵活、价格较低的专用或通用有限元分析软件，主要有德国的ASKA、英国的PA FEC、法国的SYSTUS、美国的ABQUS、ADINA、ANSYS、BERSAFE、BOSOR、COSMOS、ELAS、MARC和STARDYNE等公司的产品。 以下对一些常用的软件进行一些比较分析： 1. LSTC公司的LS-DYNA系列软件
板结构有限元分析实例详解1：带孔平板结构静力分析本节介绍带孔平板结构静力分析问题，同时介绍布尔操作的基本用法。 8.3.1 问题描述与分析 有孔的矩形平板，左侧边缘固定，长400mm，宽200 mm，厚度为10 mm，圆孔在板的正中心，半径为40 mm，左侧全约束，右侧边缘均布应力1MPa，如图8.7所示。求板的变形、位移及应力变化情况。（材料的材料属性为：弹性模量为300000 MPa，剪切模量为0.31。） 图8.7 带孔的矩形平板 由于小孔处边缘不规则，本文采用PLANE82高阶平面单元进行分析。 8.3.2 求解过程 8.3.2.1 定义工作目录及文件名 启动ANSYS Mechanical APDL Product Launcher窗口，如图8.8所示。在License下 拉选框中选择ANSYS Multiphysics产品，在Working Directory输入栏中输入工作目 录：C:\ANSYS12.0 Structural Finite Elements Analysis and Practice\Chapter 8\8-1，在Job Name一栏中输入工作文件名：Chapter8-1。以上参数设置完毕后，单 击Run按钮运行ANSYS。
图8.8 ANSYS设置窗口菜单 可以先在目标文件位置建立工作目录，然后单击Browse按钮选择工作目录；也 可以通过单击Browse按钮选择工作文件名。 8.3.2.2 定义单元类型和材料属性 选择Main Menu>Preferences命令，出现Preferences for GUI Filtering对话框， 如图8.9所示，在Individual discipline(s) to show in the GUI中勾选Structural，过滤掉ANSYS GUI菜单中与结构分析无关的选项，单击OK按钮关闭该对话框。 图8.9 Preferences for GUI Filtering对话框
概念题 1.1 结构动力计算与静力计算的主要区别是什么？ 答：主要区别表现在：(1) 在动力分析中要计入惯性力，静力分析中无惯性力；(2) 在动力分析中，结构的内力、位移等是时间的函数，静力分析中则是不随时间变化的量；(3) 动力分析方法常与荷载类型有关，而静力分析方法一般与荷载类型无关。 1.2 什么是动力自由度，确定体系动力自由度的目的是什么？ 答：确定体系在振动过程中任一时刻体系全部质量位置或变形形态所需要的独立参数的个数，称为体系的动力自由度（质点处的基本位移未知量）。确定动力自由度的目的是：(1) 根据自由度的数目确定所需建立的方程个数（运动方程数=自由度数），自由度不同所用的分析方法也不同；(2) 因为结构的动力响应（动力内力和动位移）与结构的动力特性有密切关系，而动力特性又与质量的可能位置有关。 1.3 结构动力自由度与体系几何分析中的自由度有何区别？ 答：二者的区别是：几何组成分析中的自由度是确定刚体系位置所需独立参数的数目，分析的目的是要确定体系能否发生刚体运动。结构动力分析自由度是确定结构上各质量位置所需的独立参数数目，分析的目的是要确定结构振动形状。 1.4 结构的动力特性一般指什么？ 答：结构的动力特性是指：频率（周期）、振型和阻尼。动力特性是结构固有的，这是因为它们是由体系的基本参数（质量、刚度）
所确定的、表征结构动力响应特性的量。动力特性不同，在振动中的响应特点亦不同。 1.5 什么是阻尼、阻尼力，产生阻尼的原因一般有哪些？什么是等效粘滞阻尼？ 答：振动过程的能量耗散称为阻尼。 产生阻尼的原因主要有：材料的内摩擦、构件间接触面的摩擦、介质的阻力等等。当然，也包括结构中安装的各种阻尼器、耗能器。阻尼力是根据所假设的阻尼理论作用于质量上用于代替能量耗散的一种假想力。粘滞阻尼理论假定阻尼力与质量的速度成比例。粘滞阻尼理论的优点是便于求解，但其缺点是与往往实际不符，为扬长避短，按能量等效原则将实际的阻尼耗能换算成粘滞阻尼理论的相关参数，这种阻尼假设称为等效粘滞阻尼。 1.6 采用集中质量法、广义位移法（坐标法）和有限元法都可使无限自由度体系简化为有限自由度体系，它们采用的手法有何不同？ 答：集中质量法：将结构的分布质量按一定规则集中到结构的某个或某些位置上，认为其他地方没有质量。质量集中后，结构杆件仍具有可变形性质，称为“无重杆”。 广义坐标法：在数学中常采用级数展开法求解微分方程，在结构动力分析中，也可采用相同的方法求解，这就是广义坐标法的理论依据。所假设的形状曲线数目代表在这个理想化形式中所考虑的自由度个数。考虑了质点间均匀分布质量的影响（形状函数），一般来说，
1 ： 图 a 桁 架, 力 法 基 本 结 构 如 图 b ，力 法 典 型 方 程 中 的 系 数 为 ：( ) 3. 2：图示结构用力矩分配法计算时，结点A 的约束力矩（不平衡 力矩）为（以顺时针转为正） ( ) 4.3Pl/16 3：图示桁架1,2杆内力为： 4. 4：连续梁和 M 图如图所示，则支座B 的竖向反力 F By 是：
4.17.07(↑) 5：用常应变三角形单元分析平面问题时，单元之间（）。 3.应变、位移均不连续； 6：图示体系的几何组成为 1.几何不变,无多余联系； 7：超静定结构在荷载作用下的内力和位移计算中，各杆的刚度为（） 4.内力计算可用相对值，位移计算须用绝对值 8：图示结构用力矩分配法计算时，结点A之杆AB的分配系数
μAB 为（各杆 EI= 常数）( ) 4.1/7 9：有限元分析中的应力矩阵是两组量之间的变换矩阵，这两组量是（ ）。 4.单元结点位移与单元应力 10：图示结构用位移法计算时，其基本未知量数目为( ) 4.角位移=3，线位移=2 11：图示结构，各柱EI=常数，用位移法计算时，基本未知量数 目是（ ） 3.6 12：图示结构两杆长均为d,EI=常数。则A 点的垂直位移为（ ） 4.qd 4/6EI （↓） 13：图示桁架,各杆EA 为常数,除支座链杆外,零杆数为：
1.四 根 ； 14：图示结构，各杆线刚度均为i,用力矩分配法计算时，分配 系数μAB 为（ ） 2. 15：在位移法中，将铰接端的角位移，滑动支撑端的线位移作为基本未知量： 3.可以，但不必； 1：用图乘法求位移的必要条件之一是：( ) 2.结构可分为等截面直杆段； 2：由于静定结构内力仅由平衡条件决定,故在温度改变作用下静定结构将（ ） 2.不产生内力 3：图示结构，各杆EI=常数，欲使结点B 的转角为零，比值P1/P2应 为（ ） 2.1
岩土工程有限元分析软件 PLAXIS 2D 2015? 案例教程 北京市古城西街19号研发主楼4层，100043
目录 新奥法（NATM）隧道开挖 (1) 1.1 输入 (2) 1.2生成网格 (5) 1.3计算 (6) 1.4 结果 (8)
新奥法（NATM）隧道开挖 本例利用PLAXIS分析NATM隧道施工过程。NATM是在地下开挖时，利用喷射混凝土作为临时支护，保证开挖稳定性的一种施工方法。 图1.1 项目几何尺寸 目标： ●模拟NATM隧道施工（β法）。 ●用重力加载生成初始应力。
PLAXIS 2D AE案例教程：新奥法（NATM）隧道开挖 1.1 输入 1.1.1一般设置 ●打开PLAXIS 2D AE软件，在出现的快速选择对话框中选择一个新的项目。 ●在工程属性窗口的工程标签下，键入一个合适标题。 ●在模型标签下，模型（平面应变）和单元（15-Node）保持默认选项。 ●保持单位和一般设置框为默认值。 ●在几何形状设定框中设定土层模型尺寸xmin=-50，xmax=50，ymin=0，ymax=35。 ●点击OK即关闭工程属性窗口，完成设定。 1.1.2土层定义 利用钻孔生成土层，模型中考虑11m厚的泥灰岩，这层的底部y min=0作为参考点，定义土层： 在x=-22处创建第一个钻孔。 ●修改土层窗口将出现。为钻孔添加三层土。钻孔Borehole_1第一层的深度为0.指 定第一层土的顶部和底部值为24。第二层土层的顶部=24和底部=11.第三层土层的 顶部=11和底部=0。 ●单击在修改土层窗口的底部钻孔按钮。 ●在出现的菜单中选择添加选项。添加钻孔窗口出现。 ●指定第二个钻孔的位置为x=-14. ●注意：钻孔Borehole_1的特性复制给了Borehole_2。 ●Borehole_2第一层的深度也是0。修改土层的顶部=30和底部=30.第二层土顶部=30 和底部=11。第三层土顶部=11和底部=0. ●指定第三个钻孔的位置为x=-7. ●Borehole_3第一层土顶部=35和底部=30.第二层土顶部=30和底部=11,。第三层土 的顶部=11和底部=0. ●所有钻孔设置水头高度为y=0m。土层分布如图1.2。 ●根据表1.1定义土层材料属性，并分别指定给相应土层（图1.2）. ●关闭修改土层窗口，切换到结构模式定义结构单元。 图1.2 土层分布
有限元分析及应用大作业 作业要求： 1）个人按上机指南步骤至少选择习题中3个习题独立完成，并将计算结果上交； 也可根据自己科研工作给出计算实例。 2）以小组为单位完成有限元分析计算； 3）以小组为单位编写计算分析报告； 4）计算分析报告应包括以下部分： A、问题描述及数学建模； B、有限元建模（单元选择、结点布置及规模、网格划分方案、载荷及边界 条件处理、求解控制） C、计算结果及结果分析（位移分析、应力分析、正确性分析评判） D、多方案计算比较（结点规模增减对精度的影响分析、单元改变对精度的 影响分析、不同网格划分方案对结果的影响分析等） 题一：图示无限长刚性地基上的三角形大坝，受齐顶的水压力作用，试用三节点常应变单元和六节点三角形单元对坝体进行有限元分析，并对以下几种计算方案进行比较： 1）分别采用相同单元数目的三节点常应变单元和六节点三角形单元计算；（注意ANSYS中用四边形单元退化为三节点三角形单元） 2）分别采用不同数量的三节点常应变单元计算； 3）当选常应变三角单元时，分别采用不同划分方案计算。 解：1.建模： 由于大坝长度>>横截面尺寸，且横截面沿长度方向保持不变，因此可将大坝看作无限长的实体模型，满足平面应变问题的几何条件；对截面进行受力分析，作
用于大坝上的载荷平行于横截面且沿纵向方向均匀分布，两端面不受力，满足平面应变问题的载荷条件。因此该问题属于平面应变问题，大坝所受的载荷为面载荷，分布情况P=98000-9800*Y；建立几何模型，进行求解；假设大坝的材料为钢，则其材料参数：弹性模量E=2.1e11，泊松比σ=0.3； 2：有限元建模过程： 2.1 进入ANSYS : 程序→ANSYS APDL 15.0 2.2设置计算类型: ANSYS Main Menu: Preferences →select Structural →OK 2.3选择单元类型: ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 4node 182(三节点常应变单元选择Solid Quad 4node 182，六节点三角形单元选择Solid Quad 8node 183)→OK (back to Element Types window) →Option →select K3: Plane Strain →OK→Close (the Element Type window) 2.4定义材料参数： ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic →input EX:2.1e11, PRXY:0.3 →OK 2.5生成几何模型： 生成特征点： ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints→In Active CS →依次输入四个点的坐标：input:1(0,0),2(10,0),3(1,5),4(0.45,5) →OK 生成坝体截面： ANSYS Main Menu: Preprocessor →Modeling →Create →Areas →Arbitrary →Through KPS →依次连接四个特征点，1(0,0),2(6,0),3(0,10) →OK 2.6 网格划分： ANSYS Main Menu: Preprocessor →Meshing →Mesh Tool→(Size Controls) lines: Set →依次拾取两条直角边:OK→input NDIV: 15 →Apply→依次拾取斜边:OK →input NDIV: 20 →OK →(back to the mesh tool window)Mesh:Areas, Shape: tri, Mapped →Mesh →Pick All (in Picking Menu) →Close( the Mesh Tool window) 2.7 模型施加约束： 给底边施加x和y方向的约束： ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Displacement →On lines →pick the lines →OK →select Lab2:UX, UY →OK 给竖直边施加y方向的分布载荷： ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result窗口中出现{X},写入所施加的载荷函数： 98000-9800*{Y}；3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file：将需要的.func文件打开，参数名取meng，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取竖直边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷为meng参数名→OK 2.8 分析计算： ANSYS Main Menu: Solution →Solve →Current LS →OK(to close the solve Current Load
浅谈结构力学在结构设计中的体现 摘要： 随着计算在工程上应用的日益广泛，结构设计是把数学上最优化理论结合计算机技术应用于结构设计。结构计算简图的选择经历一个复杂的过程，需要各种力学知识并结合工程实践经验，经过科学抽象、实验论证，根据实际受力、变形规律等主要因素，对结构进行合理简化。 关键词： 结构力学结构设计应用 1前言 结构力学是固体力学的一个分支，它主要研究工程结构受力和传力的规律，以及如何进行结构优化的学科。所谓工程结构是指能够承受和传递外载荷的系统，包括杆、板、壳以及它们的组合体，如桥梁、屋架和承重墙等。 随着现代经济的发展，高层建筑及各种地下复杂结构也逐步增多，结构力学的在工程上应用也越来越广泛，当然这也促进了结构理论的发展。特别是20世纪中叶，随着电子计算机和有限元法的问世使得大型结构的复杂计算成为可能，从而将结构力学的研究和应用水平提到了一个新的高度。结构力学是一门古老的学科，又是一门迅速发展的学科。随着新型工程材料和新型工程结构的大量出现，向结构力学提供了新的研究内容并提出新的要求。计算机的发展，为结构力学提供了有力的计算工具，另一方面，结构力学对数学及其他学科的发展也起了推动作用。有限元法这一数学方法的出现和发展就与结构力学的研究有密切关系。 2结构力学的重要性 实际结构是很复杂的，在对实际结构（如高层建筑、大跨度桥梁、大型水工结构）进行力学分析和计算之前必须加以简化，用一个简化图形（结构计算简图）来代替实际结构，略其次要细节，显示其基本特点，作为力学计算的基础，这一过程通常称为力学建模，用于结构计算的称为计算简图。
计算简图由实际结构简化抽象而成，取杆件轴线，或板壳中面，或块体轮廓加上结构内部的结点、结线联系，或外部的支杆、支座等边界约束，并考虑简化或分配的荷载，构成力学计算模型。 结构计算简图的选择经历一个复杂的过程，需要力学知识、结构知识、工程实践经验和洞察力，经过科学抽象、实验论证，根据实际受力、变形规律等主要因素，对结构进行合理简化。它不仅与结构的种类、功能有关，而且与作用在结构上的荷载、计算精度要求、结构构件的刚度比、安装顺序、实际运营状态及其它指标有关。计算简图的选择可能因计算状态（是考虑强度或刚度，计算稳定或振动，还是钢筋混凝土抗裂验算）而异，也依赖于所要采用的计算理论和计算方法，方能完成结构构件线性或非线性的应力和应变状态分析。实用上可以参考同类工程实例。 结构设计是先有“设想”后有“计算”，“设想”是建立在定性分析的基础上。力学始于定性分析，终于定性分析；定性分析在先，定量分析在后；定性失准，定量准偏。在进行工程设计和处理工程实际问题时，需要设计人员对结构的合理形式以及相应的结构变形和内力等具有总体概念和定性分析能力，还需要具有对工程中计算的数据、发生的现象和出现的问题能够做出迅速科学判断的能力，这就是所谓概念设计和概念分析理念。 结构力学是一切工程进行设计的基础。实际工程中都是将工程实践中的实际问题抽象为相应的力学计算公式进行求解；作为工程技术设计人员应该掌握工程结构的基本理论和实用设计方法，具备根据建筑工程项目的特点、性质、功能和业主的要求正确、合理地进行工程结构设计的基本能力。 2在xx中的应用 中国以木结构为主体的古建筑，在世界建筑之林中独树一帜。木结构它以木构为骨、砖石为体、结瓦为盖、油饰彩绘为衣，经历代能工巧匠精心设计，巧妙施工，潜心装饰，付诸心血和智慧建造而成，体现出东方古典建筑独有的艺术魅力和中国古建筑木结构的历史性、艺术性和科学性。 巧妙而科学的框架式结构是中国古代建筑在建筑结构上最重要的一个特征。因为中国古代建筑主要是木构架结构，即采用木柱、木梁构成房屋的框
有限元软件与多体动力学软件 数值分析技术与传统力学的结合在结构力学领域取得了辉煌的成就，出现了以ANSYS 、NASTRAN 等为代表的应用极为广泛的结构有限元分析软件。计算机技术在机构的静力学分析、运动学分析、动力学分析以及控制系统分析上的应用，则在二十世纪八十年代形成了计算多体系统动力学，并产生了以ADAMS 和DADS 为代表的动力学分析软件。两者共同构成计算机辅助工程（CAE ）技术的重要内容。 商业通用软件的广泛应用给我们工程师带来了极大的便利，很多时候我们不需要精通工程问题中的力学原理，依然可以通过商业软件来解决问题，不过理论基础的缺失还是会给我们带来不少的困扰。随着动力有限元与柔性多体系统分析方法的成熟，有时候正确区分两者并不是很容易。 机械领域应用比较广泛的有两类软件，一类是有限元软件，代表的有：ANSYS, NASTRAN, ABAQUS, LS-DYNA, Dytran 等；另一类是多体动力学软件，代表的有ADAMS, Recurdyn ， Simpack 等。在使用时，如何选用这两类软件并不难，但是如果深究这两类软件根本区别并不容易。例如，有限元软件可以分析静力学问题，也可以分析“动力学”问题，这里的“动力学”与多体动力学软件里面的动力学一样吗？有限元软件在分析动力学问题时，可以模拟物体的运动，它与多体动力学软件中模拟物体运动相同吗？多体动力学软件也可以分析柔性体的应力、应变等，这与有限元软件分析等价吗？ 1 有限元软件 有限单元法是一种数学方法，不仅可以计算力学问题，还可以计算声学，热，磁等多种问题，我们这里只探讨有限元法在机械领域的应用。 计算结构应力、应变等的力学基础是弹性力学，弹性力学亦称为弹性理论，主要研究弹性体在外力作用或温度变化等外界因素下所产生的应力、应变和位移，从而为工程结构或构件的强度、刚度设计提供理论依据和计算方法。也就是说用有限元软件分析力学问题时，是用有限元法计算依据弹性力学列出的方程。 考虑下面这个问题，在()0t ， 时间内给一个结构施加一个随时间变化的载荷()P t ，我们希望得到结构的应力分布，在刚刚施加载荷的时候，结构中的应力会有波动，应力场是变化的，但很久以后，应力场趋于稳定。 如果我们想得到载荷施加很久以后，稳定的应力场分布，那么应该用静力学分析方法分析
各种有限元分析软件比较 有限元分析（FEA，Finite Element Analysis）利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。还利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 有限元分析具有确保产品设计的安全合理性，同时采用优化设计，找出产品设计最佳方案，降低材料的消耗或成本; 在产品制造或工程施工前预先发现潜在的问题; 模拟各种试验方案，减少试验时间和经费等作用，越来越被应用，越来越的有限元分析也不断被开发出来，当我们在做有限元分析时，我们该选择什么样的软件？或者我们该学习什么软件？成了大多数人困惑的问题。看板网根据自己超过十年的有限元分析项目经验和培训经验，对各种有限元分析软件进行了一些比较，希望大家在选择时能够大家做参考。 有限元分析常用软件 国外软件 大型通用有限元商业软件：如ANSYS可以分析多学科的问题，例如：机械、电磁、热力学等；电机有限元分析软件NASTRAN等。还有三维结构设计方面的UG,CATIA,Proe等都是比较强大的。 国内软件 国产有限元软件：FEPG,SciFEA,JiFEX,KMAS等。 当然首先要明确你要用这个软件进行什么分析，一般会用到有限元分析的地方有：1.模流分析；2.结构强度分析；3.电磁场分析；4.谐响应分析（比如查找共振频率）；5. 铸造分析。等等 ANSYS是商业化比较早的一个软件，目前公司收购了很多其他软件在旗下。ABAQUS专注结构分析目前没有流体模块。MSC是比较老的一款软件目前更新速度比较慢。ADINA是在同一体系下开发有结构、流体、热分析的一款软件，功能强大但进入中国时间比较晚市场还没有完全铺开。 workbench是一个综合性的有限元分析软件，几乎囊括了所有有限元分析领域，传统的优势领域有强度分析、谐响应分析和电磁分析。workbench是ansys
对称结构有限元分析 ----3节点三角形单元的分析 一问题分析(对称框架线弹性实体的静力平衡问题) 图是一个方形弹性实体，单位边长、单位厚度、承受等效竖向压力2 1m，其中边界条 KN 件暗示着存在两组相对称的平面，因此现考虑的仅是问题的。每个节点上的自由度号码代表了各自在x和y方向上可能的位移。 结构和单元信息NELS NCE NN NIP 8 2 9 1 AA BB E V
.5 .55 1.E6 .3 约束节点自由度信息NR 5 K , NF(:,K), I=1，NR 10 1 4 0 1 7 0 0 8 1 9 1 0 载荷信息LOADED_NODES 3 (K, LOADS(NF(:,K)), I=1 , LOADED_NODES) 1 .0 -.25 2 .0 -.5 3 .0 -.25 333 3节点三角形单元网络的总体节点和单元编号 3节三角形单元局部坐标系中节点和自由度编号
二理论基础（有限元方法原理） 通过弹性力学变分原理建立弹性力学问题有限元方法表达格式的基本步骤。最小位能原理的未知场变量是位移，以结点位移为基本未知量，并以最小位能原理为基础建立的有限元为位移元。它是有限元方法中应用最为普遍的单元，也是本书主要讨论的单元。 对于一个力学或无力问题，在建立其数学模型以后，用有限元方法对它进行分析的首要步骤是选择单元形式。平面问题3结点三角形单元是有限元方法最早采用，而且至今仍经常采用的单元形式。我们将以它作为典型，讨论如何应用广义坐标建立单元位移模式与位移插值函数，以及如何根据最小位能原理建立有限元求解方程的原理、方法与步骤，并进而引出弹性力学问题有限元方法的一般表达格式。对于前一问题，着重讨论选择广义坐标和有限元位移模式的一般原则和建立其位移插值函数的一般步骤。对于后一问题，着重讨论单元刚度矩阵和单元载荷向量的形式，总体刚度矩阵和总体载荷向量集成的原理和方法，以及它们各自的特性。 作为一种数值方法，有限元解的收敛性无疑是十分重要的问题，以后将讨论解的收敛准则及其物理意义，所阐明的原则在以后还将得到进一步的应用和具体化。 在建立了有限元的一般表达格式以后，原则上可以将它推广到平面问题以外的其他弹性力学问题和采用任何形式的单元。轴对称问题具有很广泛的应用领域，轴对称问题3结点三角形 单元的表达格式可以看作是平面问题此种单元表达格式的直接推广。 一）弹性力学平面问题的有限元格式 结点三角形单元是有限元方法中最早提出，并且至今仍广泛应用的单元，由于三角形单元对复杂边界有较强的适应能力，因此很容易将一个二维离散成有限个三角形单元，如图1所示。在边界上以若干段直线近似原来的曲线边界，随着单元增多，这种拟合将趋于精确。我们在讨论如何应用有限元方法分析各类具体问题的开始，将以平面问题3结点三角形单元 为例来阐明弹性力学问题有限元分析的表达格式和一般步 1.1）单元位移模式及插值函数的构造 典型的3节点三角形单元节点编码i,j,m ，以逆时针方向编码为正向。每个节点有位移分量如图所示。 ?? ? ???=i i v u i a (i,j,m) 每个单元有6个节点位移即6个节点自由度，亦即 [ ] T m m j j i i m j i e v u v u v u a a a =??? ? ??????=a 1.2） 单元的位移模式和广义坐标 在有限元方法中单元的位移模式或称位移函数一般采用多项式作为近似函数，因为 多项式运算简便，并且随着项数的增多，可以逼近任何一段光滑的函数曲线。多项式的选取由低次到高次。
2 受料仓与给料机的钢结构有限元分析 2.1建立有限元模型 如图2.1破碎站主视图和图2.2破碎机布置图，它的工作过程是：卸料卡车间歇把最大入料粒度为1500mm的煤块倒入受料仓，受料仓存储大粒度煤块。刮板给料机把受料仓的大粒度的煤块连续的刮给破碎平台的破碎机。破碎机把最大入料粒度为1500mm的煤块破碎成最大排料粒度为300mm的煤块，煤块由底部的传送带传出。 图2.1 破碎站主视图
图2.2 破碎机布置图 破碎站钢结构的弹性模量E＝200000MPa，泊松比，质量密度 ×10-3kg/cm3。破碎站由支撑件型钢和斜支撑角钢组成。在结构离散化时，由于角钢和其它部位铰接，铰接是具有相同的线位移，而其角位移不同。承受轴向力，不承受在其它方向的弯矩，相当于二力杆，所以型钢用梁单元模拟，角钢用杆单元模拟。破碎站是由受料仓与给料机和破碎平台与控制室两部分组成，故计算时是分别对这两部分进行的。离散后，受料仓和给料机共个单元，其中梁单元个，杆单元个，节点总数为个，有限元模型如图和图所示。
图2.3 受料仓与给料机有限元模型 图2.4 受料仓与给料机有限元模型俯视图 2.2载荷等效计算 2.2.1主要结构截面几何参数 破碎站主要结构采用H型钢梁，截面尺寸如图2.5所示，各截面横截面积A，截面惯性矩I y，I z和极惯性矩I如下。
图2.5 截面尺寸 料仓及给料机支撑结构 料仓及给料机六根支撑立柱(H500×400×12×20) A= 215.2mm2，I y＝101947×104mm4，I z＝21340×104mm4，I＝240×104mm4料仓B－B面横梁和给料机E－E、F－F面横梁(H400×300×12×20) A=16320mm2，I y＝48026×104mm4，I z＝9005×104mm4，I＝181×104mm4料仓C－C面和D－D面横梁(H400×400×12×20) A=20320mm2，I y＝62479×104mm4，I z＝21339×104mm4，I＝234×104mm4给料机两根纵梁(H550×400×12×20) A=22120mm2，I y＝125678×104mm4，I z＝21341×104mm4，I＝243×104mm4给料机六根横梁(H400×400×12×20) A=20320mm2，I y＝62479×104mm4，I z＝21339×104mm4，I＝234×104mm4其它横梁(H400×300×12×20) A=16320mm2，I y＝48026×104mm4，I z＝9005×104mm4，I＝181×104mm4 斜支撑的横截面积 ∠125×12：A＝2856mm2 ∠75×6：A＝864mm2 2.2.1实际载荷情况