钢筋混凝土有限元分析的几个问题
1、钢筋混凝土有限元模型（引自江见鲸等《混凝土结构有限元分析》）
钢筋混凝土有限元模型包括实体单元模型和杆系单元模型。实体单元模型多用于分析独立的构件，如一个梁或一个节点，以及一个连续大型结构（水坝、核反应堆安全壳等）；而杆系单元模型多用于常见的建筑工程，如框架结构、框剪结构等，这种整体结构往往有成千上万个构件组成，用实体单元对整个结构建模分析几乎是不可能的。
实体单元模型包括分离式、组合式和整体式。三种模型，只有分离式能考虑粘结和滑移。
分离式模型：将混凝土和钢筋作为不同的单元来处理，可以插入联结单元模拟钢筋和混凝土之间的粘结和滑移，也可以不用联结单元，直接将钢筋和混凝土单元共节点或进行耦合。根据需要决定。
组合式模型：包括分层组合式、带钢筋的四边形单元、带钢筋膜的8节点六面体单元。这种模型运用较少。
整体式模型：将钢筋分布于整个单元中，并把单元视为连续均匀材料，单元刚度矩阵综合了混凝土与钢筋单元的刚度矩阵，将钢筋化为等效的混凝土，然后计算单元刚度。一般实际工程均采用这种模型，单元划分少，计算量小，可适应复杂配筋。缺点是只能求得钢筋的平均应力，不能计算粘结应力。
钢筋混凝土杆系单元模型：通常用梁单元模拟结构中的梁柱构件，用带刚域的梁单元模拟剪力墙和连梁，用桁架单元模拟支撑。根据各个杆件的类型、尺寸和材料组成（混凝土强度、钢筋强度、配筋率、钢筋布置方式）、受力工况，设定相应的非线性恢复力关系（滞回曲线），然后输入荷载工况，得到整个结构的非线性行为。杆系有限元分析要解决两个核心问题：选择合适的杆件刚度模型，选择合适的恢复力关系。
现阶段对钢筋混凝土杆系结构的刚度计算，主要有三种方法：
A、平均刚度模拟
B、分布刚度模拟
C、集中塑性铰模拟
对于框架结构中众多的杆系构件，在地震作用下往往有这样的破坏机制：反弯点在中间，塑性铰在两端。因此，可以认为塑性变形都集中在两端塑性铰区，中段构件保持弹性。
2、引自https://www.doczj.com/doc/7810617774.html,/cgi-bin/ut/topic_show.cgi?id=160085&h=1#711329混凝土本构关系应如何设置
Q：我在进行的钢管混凝土滞回曲线分析，见过有些文章中在对混凝土本构关系的定义只定义concrete（即五参数破坏准则），而不定义等向或随动的强化模型的本构关系。我不知道这样不定义屈服准则计算混凝土对不对？
我在滞回分析中，当不定义屈服准则时，能够得到滞回曲线，只是形状不好，没有强度和刚度的退化，骨架曲线没有下降段。
而当定义了屈服准则过后，整个模型就没法进行计算。
想请问的是，如果我不定义混凝土的屈服准则，而只定义破坏准则，这样做对不对？恳请高手指点，谢过！！
A1:：solid65号单元本身就有默认的本构关系，如果你不定义本构关系，可能就用本身的本构关系，这与你实际的本构关系有差别，所以得到的滞回曲线有问题，建议最好还是自己定义
对于下降段的问题的确很难处理，很多时候本构关系中定义下降段的时候，最后计算不能收敛，遇到这种情况的时候，首选还是将下降段拉直吧，因为一般的下降段在分析的时候用不到，没有什么实际意义（少数情况例外，比如分析混凝土压碎的问题），还有一个就是考虑好收敛准则，最好是多种准则结合使用，使用单一的收敛准则很容易不收敛，当然有的时候本构关系中有下降段又使用单一的收敛准则的时候还是能收敛的。
MATLAB: MATLAB是美国MathWorks公司出品的商业数学软件，用于数据分析、无线通信、深度学习、图像处理与计算机视觉、信号处理、量化金融与风险管理、机器人，控制系统等领域。 MATLAB是matrix&laboratory两个词的组合，意为矩阵工厂（矩阵实验室），软件主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连接其他编程语言的程序等。MATLAB的基本数据单位是矩阵，它的指令表达式与数学、工程中常用的形式十分相似，故用MATLAB来解算问题要比用C，FORTRAN等语言完成相同的事情简捷得多，并且MATLAB也吸收了像Maple等软件的优点，使MATLAB成为一个强大的数学软件。在新的版本中也加入了对C，FORTRAN，C++，JAVA的支持。 MATLAB有限元分析与应用:
《MATLAB有限元分析与应用》是2004年4月清华大学出版社出版的图书，作者是卡坦，译者是韩来彬。 内容简介： 《MATLAB有限元分析与应用》特别强调对MATLAB的交互应用，书中的每个示例都以交互的方式求解，使读者很容易就能把MATLAB用于有限分析和应用。另外，《MATLAB有限元分析与应用》还提供了大量免费资源。 《MATLAB有限元分析与应用》采用当今在工程和工程教育方面非常流行的数学软件MATLAB来进行有限元的分析和应用。《MATLAB有限元分析与应用》由简单到复杂，循序渐进地介绍了各种有限元及其分析与应用方法。书中提供了大量取自机械工程、土木工程、航空航天工程和材料科学的示例和习题，具有很高的工程应用价值。
Ansys有限元分析实例[教学] 有限元分析案例:打点喷枪模组(用于手机平板电脑等电子元件粘接)，该产品主要是使用压缩空气推动模组内的顶针作高频上下往复运动，从而将高粘度的胶水从喷嘴中打出(喷嘴尺寸,0.007”)。顶针是这个产品中的核心零件，设计使用材料是:AISI 4140 最高工作频率是160HZ(一个周期中3ms开3ms关)，压缩空气压力3-8bar, 直接作用在顶针活塞面上，用Ansys仿真模拟分析零件的强度是否符合要求。 1. 零件外形设计图:
2. 简化模型特征后在Ansys14.0 中完成有限元几何模型创建:
3. 选择有限元实体单元并设定，单元类型是SOILD185，由于几何建模时使用的长度单位是mm, Ansys采用单位是长度:mm 压强: 3Mpa 密度:Ton/M。根据题目中的材料特性设置该计算模型使用的材料属性:杨氏模量 2.1E5; 泊松比:0.29; 4. 几何模型进行切割分成可以进行六面体网格划分的规则几何形状后对各个实体进行六面体网格划分，网格结果: 5. 依据使用工况条件要求对有限元单元元素施加约束和作用载荷:
说明: 约束在顶针底端球面位移全约束; 分别模拟当滑块顶断面分别以8Bar，5Bar，4Bar和3Bar时分析顶针的内应力分布，根据计算结果确定该产品允许最大工作压力范围。 6. 分析结果及讨论: 当压缩空气压力是8Bar时: 当压缩空气压力是5Bar时:
当压缩空气压力是4Bar时: 结论: 通过比较在不同压力载荷下最大内应力的变化发现，顶针工作在8Bar时最大应力达到250Mpa，考虑到零件是在160HZ高频率在做往返运动，疲劳寿命要求50百万次以上，因此采用允许其最大工作压力在5Mpa，此时内应力为156Mpa，按线性累积损伤理论[3 ]进行疲劳寿命L-N疲劳计算，进一部验证产品的设计寿命和可靠性。
目录 摘要 ............................................................................................ 错误！未定义书签。Abstract (2) 第一章分析方法和研究对象 ........................................... 错误！未定义书签。 1.1 有限单元法的概述....................................................... 错误！未定义书签。 1.1.1 有限单元法的历史 (4) 1.1.2 有限单元法的基本概念 (4) 1.2 ANSYS软件简介 (4) 1.2.1 ANSYS主要应用领域 (4) 1.2.2 ANSYS操作界面 (5) 1.2.3 ANSYS的主要功能 (6) 1.2.4 ANSYS主要特点 (7) 1.3 曲柄滑块机构简介 (7) 1.3.1 曲柄滑块定义 (8) 1.3.2 曲柄滑块机构特性应用以及分类 (8) 第二章曲柄滑块机构的求解 (10) 2.1 曲柄滑块机构的问题描述 (10) 2.2 曲柄滑块机构问题的图解法 (10) 2.2.1 图解法准备工作 (11) 2.2.2 图解法操作步骤 (11) 第三章有限元瞬态动力学概述 (14) 3.1 有限元瞬态动力学定义 (14) 3.2 瞬态动力学问题求解方法........................................... 错误！未定义书签。 3.2.1 完全法 (14) 3.2.2 模态分析法 (14) 3.2.2 缩减法 (15) 3.1 有限元结构静力学分析基本概念 (15) 3.1 有限元结构静力学分析步骤 (16) 第四章曲柄滑块的有限元瞬态动力学分析 (17) 4.1 曲柄滑块机构瞬态简要概述 (17) 4.2曲柄滑块有限元瞬态动力学分析步骤 (18)
论文题目：钢筋混凝土有限元分析技术在结构工程中的应用 学生姓名：刘畅 学号：2014105110 学院：建筑与工程学院 2015年06月30日
有限元分析在钢筋混凝土结构中的应用【摘要】在国内外的土木工程中，钢筋混凝土结构因具有普遍性、可靠性良好、操作简单等优点，而得到了广泛的应用。钢筋混凝土结构是钢筋与混凝土两种性质截然不同的材料组合而成，由于其组合材料的性质较为复杂，同时存在非线性与几何线形的特征，应用传统的解析方法进行材料的分析与描述在受力复杂、外形复杂等情况下较为困难，往往不能得到准确的数据，给工程安全带来隐患。而有限元分析方法则充分利用现代电子计算机技术，借助有限元模型有效解决了各种实际问题。 【关键词】有限元分析；钢筋混凝土结构；应用 随着计算机在工程设计领域中的广泛应用，以及非线性有限元理论研究的不断深入，有限元作为一个具有较强能力的专业数据分析工具，在钢筋混凝土结构中得到了广泛的应用。在现代建筑钢筋混凝土结构的分析中，有限元分析方法展现了较强的可行性、实用性与精确性。例如：在计算机上应用有限元分析法，对形状复杂、柱网复杂的基础筏板，转换厚板，体型复杂高层建筑侧向构件、楼盖，钢-混凝土组合构件等进行应力，应变分析，使设计人员更准确的掌握构件各部分内力与变形，进而进行设计，有效解决传统分析方法的不足，满足当前建筑体型日益复杂，工程材料多样化的实际情况。但是在有限元分析方法的应用中，必须结合钢筋混凝土结构工程的实际情况，选取作为合理的有限元模型，才能保证模拟与分析结果的真实性、精确性与可靠性。 在钢筋混凝土结构工程中，非线性有限元分析的基本理论可以概括为：1）通过分离钢筋混凝土结构中的钢筋、混凝土，使其成为有限单位、二维三角形单元，钢箍离散为一维杆单元，以利于分析模型的构建；2）为了合理模拟钢筋、混凝土之间的粘结滑移关系，以及
有限元分析案例 图1 钢铸件及其砂模的横截面尺寸 砂模的热物理性能如下表所示: 铸钢的热物理性能如下表所示: 一、初始条件:铸钢的温度为2875o F，砂模的温度为80o F；砂模外边界的对流边界条件：对流系数0.014Btu/hr.in2.o F，空气温度80o F；求3个小时后铸钢及砂模的温度分布。 二、菜单操作: 1.Utility Menu>File>Change Title, 输入Casting Solidification; 2.定义单元类型:Main Menu>Preprocessor>Element Type>Add/Edit/Delete, Add, Quad 4node 55； 3.定义砂模热性能:Main Menu>Preprocessor>Material Props>Isotropic,默认材料编号1, 在Density(DENS)框中输入0.054,在Thermal conductivity (KXX)框中输入0.025,在S pecific heat(C)框中输入0.28; 4.定义铸钢热性能温度表:Main Menu>Preprocessor>Material Props>-Temp Dependent->Temp Table,输入T1=0,T2=2643, T3=2750, T4=2875; 5.定义铸钢热性能:Main Menu>Preprocessor>Material Props>-Temp Dependent ->Prop Table, 选择Th Conductivity,选择KXX, 输入材料编号2,输入C1=1.44, C2=1.54, C3=1.22, C4=1.22,选择Apply,选择Enthalpy,输入C1=0, C2=128.1, C3=163.8, C4=174.2; 6.创建关键点:Main Menu>Preprocessor>-Modeling->Create>Keypoints>In Active
Abaqus分析报告 （齿轮轴） 名称：Abaqus齿轮轴 姓名： 班级： 学号： 指导教师：
一、简介 所分析齿轮轴来自一种齿轮泵，通过用abaqus软件对齿轮轴进行有限元分析和优化。齿轮轴装配结构图如图1，分析图1中较长的齿轮轴。 图1.齿轮轴装配结构图 二、模型建立与分析 通过part、property、Assembly、step、Load、Mesh、Job等步骤建立齿轮轴模型，并对其进行分析。 1.part 针对该齿轮轴，拟定使用可变型的3D实体单元，挤压成型方式。 2.材料属性 材料为钢材，弹性模量210Gpa，泊松比0.3。
3.截面属性 截面类型定义为solid，homogeneous。 4.组装 组装时选择dependent方式。 5.建立分析步 本例用通用分析中的静态通用分析（Static，General）。 6.施加边界条件与载荷 对于齿轮轴，因为采用静力学分析，考虑到前端盖、轴套约束，而且根据理论，对受力部分和轴径突变的部分进行重点分析。 边界条件：分别在三个轴径突变处采用固定约束，如图2。 载荷：在Abaqus中约束类型为pressure，载荷类型为均布载荷,分别施加到齿轮接触面和键槽面，根据实际平衡情况，两力所产生的绕轴线的力矩方向相反，大小按比例分配。 均布载荷比计算： 矩形键槽数据： 长度：8mm、宽度：5mm、高度：3mm、键槽所在轴半径：7mm 键槽压力面积：S1 = 8x3=24mm2 平均受力半径：R1=6.5mm 齿轮数据：= 齿轮分度圆半径：R2 =14.7mm、压力角：20°、 单个齿轮受力面积：S2 ≈72mm2 通过理论计算分析，S1xR1xP1=S2xR2xP2，其中，P1为键槽均布载荷
ABAQUS钢筋混凝土损饬塑性模型有限元分析 发表时间：2009-10-12 刘劲松刘红军来源：万方数据 钢筋混凝土材料，是一种非匀质的力学性能复杂的建筑材料。随着计算机和有限元方法的发展，有限元法已经成为研究混凝土结构的一个重要的手段。由于数值计算具有快速、代价低和易于实现等诸多优点，这种分析方法已经广泛用于实际工程中。然而，要在有限元软件中尽可能准确地模拟混凝土这种材料，是不容易的，国内外学者提出了基于各种理论的混凝土本构模型。但是迄今为止，还没有一种理论被公认为可以完全描述混凝土的本构关系。 ABAQUS是大型通用的有限元分析软件，其在非线性分析方面的巨大优势，获得了广大用户的认可，在结构分析领域的应用趋于广泛。本文把规范建议的混凝土本构关系，应用到损伤塑性模型，对一悬臂梁进行了精细的有限元建模计算和探讨。 1 混凝土损伤塑性模型 ABAQUS在钢筋混凝土分析上有很强的能力。它提供了三种混凝土本构模型：混凝土损伤塑性模型，混凝土弥散裂缝模型和ABAQUS/Explicit中的混凝土开裂模型。其中混凝土损伤塑性模型可以用于单向加载、循环加载以及动态加载等场合，它使用非关联多硬化塑性和各向同性损伤弹性相结合的方式描述了混凝土破碎过程中发生的不可恢复的损伤。这一特性使得损伤塑性模型具有更好的收敛性。 2 模型材料的定义 2.1 混凝土的单轴拉压应力-应变曲线 本模型中选用的混凝土本构关系是《混凝土结构设计规范》所建议的曲线，其应力应变关系可由函数表达式定义。 2.2 钢筋的本构关系 钢筋采用本构关系为强化的二折线模型，无刚度退化。折线第一上升段的斜率，为钢筋本身的弹性模量，第二上升段为钢筋强化段，此时的斜率大致可取为第一段的1/100。 2.3 损伤的定义 损伤是指在单调加载或重复加载下，材料性质所产生的一种劣化现象，损伤在宏观方面的表现就是(微)裂纹的产生。材料的损伤状态，可以用损伤因子来描述。根据前面确定的混凝土非弹性阶段的应力一应变关系。可求得损伤因子的数值。 2.4混凝土塑性数值的计算 混凝土在单向拉伸，压缩试验中得到的数据，通常是以名义应变和名义应力表示的，为了准确地描述大变形过程中截面积的改变，需要使用真实应变和真实应力，可通过它们之间的换算公式计算。真实应变是由塑性应变和弹性应变两部分构成的。在ABAQUS中定义塑性材料参数时，需要使用塑性应变。 3 钢筋混凝土悬臂梁实例分析 3.1 模型设计 该悬臂梁的具体情况如图1所示，梁截面尺寸为200mm×300mm，梁长1500mm；纵筋为HRB335钢筋，箍筋为HPB235钢筋，混凝土强度等级为C30。混凝土和钢筋的各力学参数均取自《混凝土结构设计规范》的标准值。
用ANSYS进行四连杆机构的有限元分析 作者：谭辉 日期：08年3月6日 分析目的 1、利用ANSYS对典型的四连杆机构进行分析，主要包含各点的轨迹分 析，例如X和Y方向的位移等。 2、为五连杆和六连杆机构的分析提供可行的分析方法以及原型代码。 问题简述 分析主动杆1绕节点1旋转一周时节点4的运动轨迹，杆2和杆3为从动杆，具体问题见下图：
分析思路 1、根据分析目的，在ANSYS选用link1单元进行单元建模，主要考虑 是link1单元具有X和Y方向的自由度，可以获得各个节点的位移轨迹。 之后可以用梁单元等实现更高级的分析目的，例如获得杆上的力，位移， 加速度等相关信息。 2、该模型结构简单，可以利用直接建模方法进行有限元系统建模，主 要命令：N，E。 3、利用自由度耦合对重合节点进行建模，例如节点2和节点3、节点4 和节点5进行建模，主要命令：cpintf，利用该命令可以一次性将重合节 点生成自由度耦合。 4、利用表数组对于杆1（主动杆）的节点2进行瞬态边界条件的载荷施 加，分析类型为瞬态分析，主要命令：*dim，d等。 5、生成节点位移的对应变量，从而获得节点4的随时间的位移曲线， 主要命令：nsol，plvar等。 命令流如下 行号命令符号注释 结束上一次的分析 1finish ！ 清除数据库，并读取启动配置文件2/clear,start ！ 3 ！ 设置图形显示的背景颜色 4/color,pbak,on,1,5 ！ 5 ！
6/units,si ！ 设置单位制：国际单位制 7*afun,deg ！ 设置三角函数运算采用度为单位 8 ！ 9/prep7 ！ 进入前处理模块 10et,1,link1 ！ 设置单元类型：link1 11mp,ex,1,2.07e11 ！ 设置材料的弹性模量 12r,1,1 ！ 设置单元的实常数，面积为1 13n,1,0,0,0 ！ 在（0，0，0）处建立节点1 14n,2,3,0,0 ！ 在（3，0，0）处建立节点2 15n,3,3,0,0 ！在（3，0，0）处建立节点3，和节点2重合 16n,4,8,7,0 ！ 在（8，7，0）处建立节点4 17n,5,8,7,0 ！在（8，7，0）处建立节点4，和节点4重合 18n,6,10,0,0 ！ 在（10，0，0）处建立节点6 19e,1,2 ！ 建立单元1（连接节点1和2） 20e,3,4 ！ 建立单元2（连接节点3和4） 21e,5,6 ！ 建立单元3（连接节点5和6） 22 ！ 23cpintf,all,1e-3 ！对于重合节点一次性的建立耦合自由度，容差1e－3 24 ！ 25/pnum,node,1 ！ 显示节点编号 26/pnum,elem,1 ！ 显示单元编号 27eplot ！ 显示单元
板结构有限元分析实例详解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对话框
ABAQUS中的钢筋混凝土剪力墙建模 曲哲 2006-5-29 一、试验标定 选用ABAQUS中的塑性损伤混凝土本构模型，分离式钢筋建模，建立平面应力模型模拟钢筋混凝土剪力墙的单调受力行为。李宏男（2004）本可以提供比较理想的基准试验。然而计算发现，该文中试验记录的初始刚度普遍偏小，仅为弹性分析结果的1/5~1/8，原因不明，故此处不予采用。左晓宝（2001）研究了小剪跨比开缝墙的低周滞回性能，其中有一片整体墙作为对照试件，本文仅以这片墙为基准标定有限元模型。 图1：剪力墙尺寸与配筋 该试件尺寸及配筋如图1所示。墙全高750mm，宽800mm，厚75mm，墙内布有间距φ6@100的分布钢筋，墙两端设有暗柱。混凝土立方体抗压强度为54.9MPa，钢筋均为一级光圆筋。 （a）墙体分区及网格（b）钢筋网 图2：ABAQUS中的有限元模型 剪力墙采用平面应力八节点全积分单元，墙上下两端各加设100mm高的弹性梁。钢筋采用两节点梁单元，通过Embed方式内嵌于墙体内。模型网格及外观如图2所示。墙下弹性梁底面嵌固。分析中，先在墙顶施加160kN均布轴压力，再在墙上方弹性梁的左端缓缓施加位移荷载。 ABAQUS中损伤模型各参数取值如表1、图3所示。未说明的参数均使用ABAQUS默认值。
表1：有限元模型材料属性 混凝土 钢筋 材料非线性模型 Damaged Plasticity Plasticity 初始弹性模量（GPa ） 38.1 210 泊松比 0.2 0.3 膨胀角（deg ） 50 初始屈服应力（MPa ） 13 235 峰值压应力（MPa ） 44 峰值压应变（με） 2000 峰值拉应力（MPa ） 3.65 注：其中混凝土弹性模量为文献中提供的试验值，其余均为估计值。 （a ）压应力-塑性应变曲线 （b ）拉应力-非弹性应变曲线 （c ）受拉损伤指标-开裂应变曲线 图3：混凝土塑性硬化及损伤参数 ABAQUS 的混凝土塑性损伤模型用两个硬化参数分别控制混凝土的拉压行为，同时可以分别引入受压和受拉损伤指标。本文受压硬化曲线采用Saenz 曲线（式1），可用表1中列出的初始弹性模量、峰值应力和峰值应变唯一确定。受拉软化曲线采用Gopalaratnam 和Shah （1985）曲线（式2），并采取江见鲸建议参数k =63，λ=1.01，如图3（b ）所示。本文模型只定义受拉损伤指标，损伤指标随开裂应变的变化如图3（c ）所示，当开裂应变小于0.0014时，损伤指标线性增大，开裂应变超过0.0014后，损伤指标保持固定值0.6。 02 0000012c c c c E E εσεεεσεε= ??????+?+???????????? （1） e k t t f λ ωσ?= （2） 图4比较了采用4节点单元和8节点单元得到的剪力墙荷载-位移曲线，并同时画出了 文献中提供的荷载-位移骨架线。可见8节点单元模型的计算结果较4节点单元模型更加平滑顺畅，下降段也比较稳定。二者在达到峰值之前差别不大，但软化行为则相差较多。这可能与基于开裂应变定义的损伤指标引入的网格依赖性有关，本文对此不做深入讨论。 与试验曲线相比，有限元分析得到的荷载-位移曲线初始刚度略大，且墙底开裂（图中1点）时刚度退化不如试验中显著，导致之后的分析结果位移偏小。受拉侧钢筋屈服后计算得到的刚度与试验曲线比较接近，不久主斜裂缝的出现使墙的承载力进入软化段，被主要裂缝穿过的钢筋均进行屈服段。软化过程中墙体形成了新的主斜裂缝并最终沿这条主斜裂缝破坏。图5、6分别展示了剪力墙在受力全过程中关键点处的混凝土主拉应变和钢筋大主应力。 与试验曲线相比，计算结果刚度偏差较大，承载力基本一致。
基于ANSYS Workbench的定位卡锁机构有限元分析 摘要本文首先在Pro/E中建立了定位卡锁机构受最大外力时的简化模型，然后将该模型导入到ANSYS Workbench 13平台中进行了有限元模型的分析求解，最后结合求解结果用第四强度理论对定位卡锁机构各零件进行了强度校核，同时对该定位卡锁机构的改进提出了建议。 关键词定位卡锁机构；有限元分析 在某工程项目中应用的定位卡锁机构承担着为某输送设备准确定位的作用。由于该输送设备运行一个周期位就要启停一次，启停工作由定位卡锁机构配合实现。定位卡锁机构收回，输送设备开始运转，一个周期位后电机停转，定位卡锁机构伸出，进入与之配合的凹槽使输送设备完全停位。因此，定位卡锁机构成为该输送设备的关键部件，是保证输送设备正常工作的必备条件。所以，对定位卡锁机构的研究与分析有着重要的意义。 定位卡锁机构在伸出状态受最大外力时，其所受最大应力不应超过材料的许用应力是保证定位卡锁机构实现其功能的充分条件。为了保证定位卡锁机构的工作可靠性，本文利用ANSYS Workbench对该机构进行有限元分析，研究在定位卡锁机构受最大外力时的受力及变形情况，并依据理论知识对其强度进行校核。 1 定位卡锁机构模型的建立与导入 在对定位卡锁机构进行有限元分析之前，首先应建好定位卡锁机构的三维模型。一般在整个有限元分析的过程中，几何建模的工作量占据了非常多的时间，同时也是非常重要的过程[2]。ANSYS Workbench 13中，建模工作主要由ANSYS Workbench 自带的几何建模工具Design Modeler模块完成。对于小型或简单模型的建立可以直接在Design Modeler模块中建模，这样避免了从CAD系统中导入ANSYS的模型可能不能直接进行网格划分，需进行大量修补完善工作的麻烦。对于零部件较多的装配体的建模，通常先利用专业的三维建模软件完成模型的建立，然后再把它导入到ANSYS中进行分析。这样，工程技术人员就可以使用自己擅长的CAD软件建好模型，从而避免了重复现有CAD模型的劳动。 本文采用PTC公司的Pro/Engineer对定位卡锁机构进行三维建模。定位卡锁机构简化模型由液压缸、卡锁活塞杆、端盖、螺塞、螺钉组成，建好的三维模型如图1所示。建好后的三维模型可以在Pro/E中直接导入到ANSYS Workbench 13 中进行有限元分析。 图1 定位卡锁机构的三维模型 2 定位卡锁机构的有限元分析 2.1 定义模型材料属性
有限元分析 一个厚度为20mm的带孔矩形板受平面内张力，如下图所示。左边固定，右边受载荷p=20N/mm作用，求其变形情况 P 一个典型的ANSYS分析过程可分为以下6个步骤： ①定义参数 ②创建几何模型 ③划分网格 ④加载数据 ⑤求解 ⑥结果分析 1定义参数 1.1指定工程名和分析标题 (1)启动ANSYS软件，选择File→Change Jobname命令，弹出如图所示的[Change Jobname]对话框。 (2)在[Enter new jobname]文本框中输入“plane”，同时把[New log and error files]中的复选框选为Yes，单击确定 (3)选择File→Change Title菜单命令，弹出如图所示的[Change Title]对话框。 (4)在[Enter new title]文本框中输入“2D Plane Stress Bracket”，单击确定。 1.2定义单位
在ANSYS软件操作主界面的输入窗口中输入“/UNIT,SI” 1.3定义单元类型 (1)选择Main Menu→Preprocessor→Element Type→Add/Edit/Delete命令，弹出如图所示[Element Types]对话框。 (2)单击[Element Types]对话框中的[Add]按钮，在弹出的如下所示[Library of Element Types]对话框。 (3)选择左边文本框中的[Solid]选项，右边文本框中的[8node 82]选项，单击确定，。 (4)返回[Element Types]对话框，如下所示 (5)单击[Options]按钮，弹出如下所示[PLANE82 element type options]对话框。
基于ANSYS的齿轮传动有限元分析和优化 摘要 ANSYS是随着电子计算机的发展而迅速发展起来的一种在计算数学，计算力学和计算工程科学领域最有效的通用有限元分析软件。它是融结构，热，流体，电磁，声学于一体的大型通用有限元商用分析软件。利用ANSYS有限元分析，可以对各种机械零件，构件进行应力，应变，变形，疲劳分析，并对某些复杂系统进行仿真，实现虚拟的设计，从而大大节省人力，财力和物力。由于其方便性、实用性和有效性，ANSYS软件在各个领域，特别是机械工程当中得到了广泛的应用。 齿轮是机械中常用的一种零件，其在工作的过程中会产生应力，应变和变形，为保证其正常工作需要对齿轮的轮齿和整体受力进行分析，保证其刚度和强度的要求。本论文采用ANSYS软件对齿轮进行静力学分析和优化实现对齿轮的虚拟设计。 齿轮是最重要的零件之一。它具有功率范围大，传动效率高，传动比正确，使用寿命长等特点，但从零件失效的情况来看，齿轮也是最容易出故障的零件之一。据统计，在各种机械故障中，齿轮失效就占故障总数的60%以上。其中轮齿的折断又是齿轮失效的主要原因之一。 齿轮啮合过程作为一种接触行为, 因涉及接触状态的改变而成为一个复杂的非线性问题。传统的齿轮理论分析是建立在弹性力学基础上的, 对于齿轮的接触强度计算均以两平行圆柱体对压的赫兹公式为基础,在计算过程中存在许多假设,不能准确反映齿轮啮合过程中的应力以及应变分布与变化。相对于理论分析,有限元法则具有直观、准确、快速方便等优点。 齿廓曲面是渐开线曲面，所以建模的难点和关键在于如何确定精确的渐开线。通过PDL命令流直接在ANSYS中创建标准直齿圆柱齿轮，学习应用ANSYS软件进行零件的几何建模和网格划分，并进行静力加载和求解，对求解的结果进行查看，分析和优化。 关键词：ANSYS；有限元；齿轮；CAE
混凝土有限元分析 廖奕全 （06级防灾减灾工程及防护工程，06114249） 摘要：用传统的理论解析方法分析钢筋混凝土结构，只能解决一些非常简单的构件或结构的非线性问题，对大量的钢筋混凝土结构的非线性分析问题只能用数值方法解决，因此，有限元方法作为一个强有力的数值分析工具，在钢筋混凝土结构的非线性分析中得到了广泛地应用。随着有限元理论和计算机技术的进步，钢筋混凝土非线性有限元分析方法也得以迅速的发展并发挥出巨大的作用。 关键词：钢筋混凝土有限元分析有限元模型 钢筋混凝土结构是土木工程中应用最广泛的一种建筑结构。相比其它材料结构，钢筋混凝土结构有以下特点：①造价低，往往是建筑结构的首选材料；②易于浇注成各种形状，满足建筑功能及各种工艺的要求；⑧充分发挥钢筋和混凝土的作用，结构受力合理：④材料的重度与强度之比不大；⑤材料性能复杂，一般的计算模型难与实际结构的受力情况相符。正因为钢筋混凝土材料的这些优缺点，长期以来，钢筋混凝土在工程中的应用如此广泛；为了满足工程需要所建立的反映混凝土材料性能的计算模型也不断完善。然而，混凝土是一种由水泥、水、砂、石及各种掺合料、外加剂混合而成的成分复杂、性能多样的材料。到目前为止，还没有一种公认的、能全面反映混凝土的力学行为和性质的计算模型或本构关系。因此，对钢筋混凝土的力学性能研究还需要学术界和工程人员继续努力。长期以来，人们用线弹性理论来分析钢筋混凝土结构的受力和变形，以极限状态的设计方法来确定构件的承载能力。这种设计方法在一定程度上能满足工程的要求。随着国民经济的发展，越来越多大型、复杂的钢筋混凝土结构需要修建，而且对设计周期和工程质量也提出了更高的要求。这样一来，常规的线弹性理论分析方法用于钢筋混凝土结构和构件的设计就力不从心。设计人员常有“算不清楚”以及“到底会不会倒”的困惑。为此，钢筋混凝土非线性有限元分析方法开始受到重视。同时，随着有限元理论和计算机技术的进步，钢筋混凝土非线性有限元分析方法也得以迅速的发展并发挥出巨大的作用。 一、钢筋混凝土结构有限元分析的意义 钢筋混凝土结构是目前各种建筑结构物的主要结构形式，由于钢筋混凝土结构受到较大的荷载(如地震荷载)作用时其非线性特性对结构的性能影响很大，所以钢筋混凝土结构的非线性分析在结构抗震工程领域中十分重要并成为一个研究热点。用传统的理论解析方法分析钢筋混凝土结构，只能解决一些非常简单的构件或结构的非线性问题，对大量的钢筋混凝土结构的非线性分析问题只能用数值方法解决，因此，有限元方法作为一个强有力的数值分析工具，在钢筋混凝土结构的非线性分析中得到了广泛地应用。由于钢筋混凝土是由两种性质不同的材料——混凝土和钢筋组合而成的，它的性能明显地依赖于这两种材料的性能以及它们的相互作用，特别是在非线性阶段，混凝土钢筋本身的各种非线性性能，都不同程度地在这种组合材料中反映出来。以下是与钢筋混凝土结构计算分析有关的一些非线性问题： 1)由于钢筋和混凝土的抗拉强度相差很大，钢筋混凝土结构在正常使用状态下，大部分受弯构件都已经开裂而进入非线性状态。2)混凝土和钢筋在一个结构中共同工作的条件是两者之间的变形协调而且没有相对的滑移，但实际上，这种条件并不能完全满足，特别是在反
有 限 元 分 析 作 业 作业名称 输气管道有限元建模分析 姓 名 陈腾飞 学 号 3070611062 班 级 07机制（2）班 宁波理工学院
题目描述： 输气管道的有限元建模与分析 计算分析模型如图1所示 承受内压：1.0e8 Pa R1=0.3 R2=0.5 管道材料参数：弹性模量E=200Gpa；泊松比v=0.26。 图1受均匀内压的输气管道计算分析模型（截面图） 题目分析： 由于管道沿长度方向的尺寸远远大于管道的直径，在计算过程中忽略管道的断面效应，认为在其方向上无应变产生。然后根据结构的对称性，只要分析其中1/4即可。此外，需注意分析过程中的单位统一。 操作步骤 1.定义工作文件名和工作标题 1.定义工作文件名。执行Utility Menu-File→Chang Jobname-3070611062，单击OK按钮。 2.定义工作标题。执行Utility Menu-File→Change Tile-chentengfei3070611062，单击OK按钮。 3.更改目录。执行Utility Menu-File→change the working directory –D/chen 2.定义单元类型和材料属性 1.设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK
2.选择单元类型。执行ANSYS Main Menu→Preprocessor →Element Type→Add/Edit/Delete →Add →select Solid Quad 8node 82 →apply Add/Edit/Delete →Add →select Solid Brick 8node 185 →OK Options…→select K3: Plane strain →OK→Close如图2所示，选择OK接受单元类型并关闭对话框。 图2 3.设置材料属性。执行Main Menu→Preprocessor →Material Props →Material Models →Structural →Linear →Elastic →Isotropic，在EX框中输入2e11，在PRXY框中输入0.26，如图3所示，选择OK并关闭对话框。 图3 3.创建几何模型 1. 选择ANSYS Main Menu: Preprocessor →Modeling →Create →Keypoints →In Active CS →依次输入四个点的坐标：input:1(0.3,0),2(0.5,0),3(0,0.5),4(0,0.3) →OK
ABAQUS钢筋混凝土有限元分析 发表时间：2009-10-12 刘劲松刘红军来源：万方数据 钢筋混凝土材料，是一种非匀质的力学性能复杂的建筑材料。随着计算机和有限元方法的发展，有限元法已经成为研究混凝土结构的一个重要的手段。由于数值计算具有快速、代价低和易于实现等诸多优点，这种分析方法已经广泛用于实际工程中。然而，要在有限元软件中尽可能准确地模拟混凝土这种材料，是不容易的，国内外学者提出了基于各种理论的混凝土本构模型。但是迄今为止，还没有一种理论被公认为可以完全描述混凝土的本构关系。 ABAQUS是大型通用的有限元分析软件，其在非线性分析方面的巨大优势，获得了广大用户的认可，在结构分析领域的应用趋于广泛。本文把规范建议的混凝土本构关系，应用到损伤塑性模型，对一悬臂梁进行了精细的有限元建模计算和探讨。 1 混凝土损伤塑性模型 ABAQUS在钢筋混凝土分析上有很强的能力。它提供了三种混凝土本构模型：混凝土损伤塑性模型，混凝土弥散裂缝模型和ABAQUS/Explic it中的混凝土开裂模型。其中混凝土损伤塑性模型可以用于单向加载、循环加载以及动态加载等场合，它使用非关联多硬化塑性和各向同性损伤弹性相结合的方式描述了混凝土破碎过程中发生的不可恢复的损伤。这一特性使得损伤塑性模型具有更好的收敛性。 2 模型材料的定义 2.1 混凝土的单轴拉压应力-应变曲线 本模型中选用的混凝土本构关系是《混凝土结构设计规范》所建议的曲线，其应力应变关系可由函数表达式定义。 2.2 钢筋的本构关系 钢筋采用本构关系为强化的二折线模型，无刚度退化。折线第一上升段的斜率，为钢筋本身的弹性模量，第二上升段为钢筋强化段，此时的斜率大致可取为第一段的1/100。 2.3 损伤的定义 损伤是指在单调加载或重复加载下，材料性质所产生的一种劣化现象，损伤在宏观方面的表现就是(微)裂纹的产生。材料的损伤状态，可以用损伤因子来描述。根据前面确定的混凝土非弹性阶段的应力一应变关系。可求得损伤因子的数值。 2.4混凝土塑性数值的计算 混凝土在单向拉伸，压缩试验中得到的数据，通常是以名义应变和名义应力表示的，为了准确地描述大变形过程中截面积的改变，需要使用真实应变和真实应力，可通过它们之间的换算公式计算。真实应变是由塑性应变和弹性应变两部分构成的。在ABAQUS中定义塑性材料参数时，需要使用塑性应变。 3 钢筋混凝土悬臂梁实例分析 3.1 模型设计 该悬臂梁的具体情况如图1所示，梁截面尺寸为200mm×300mm，梁长1500mm；纵筋为HRB335钢筋，箍筋为HPB235钢筋，混凝土强度等级为C30。混凝土和钢筋的各力学参数均取自《混凝土结构设计规范》的标准值。
2009 年 9 月 第 6 卷 第 3 期 深圳土木与 与建筑 VOL.6?NO O.3?SEP2009 9
有限元分 有 分析中钢筋 筋混凝土 土梁的预 预应力模拟 拟?
陈宜 宜言 林松 尧国皇 尧
（深圳市市 市政设计研究院 院有限公司）
【摘要】通 ABAQUS软 通过 软件对钢筋混 混凝土梁的预 预应力进行模 模拟，通过计 计算结果进行 行比较分析知 知：用降温度 度 法和初始应 应力法施加预 预应力时，所 所得到的预应 应力效果相同 同，相当于先 先张法，而reb 施加初始 bar 始应力法相当 于后张法。 。 【关键词】有限元 ABA AQUS 预应力 力模拟
1 引言?
当梁的 的跨度或荷载较大时，其变 变形和裂缝宽 宽度 可能无法满 满足正常使用要求，研究和 和工程经验表 表明， 使用预应力 力技术能较好 好地解决上述问题，同时还 还可 增加梁的弹 弹性工作范围，提高承载力 力，充分利用 用材 料性能，从 从而降低结构高度、减轻自重、减小地 地震 作用，增加 加强度储备，延长使用期限 1 ，因此预 预应 力钢筋混凝 凝土梁得到广泛应用，本文 文将采用通用 用有 限元软件 ABAQUS 中不同的方法对钢 钢筋混凝土梁 梁的 行模拟，并探讨不同计算结 结果之间的区 区别。 预应力进行
（a） 整体 体有限元?
2 有限元模 模型?
对 一 跨 度 为 8000 ， 截 面 高 800 ，宽 400 的预 预应力钢筋混 混凝土简支梁 梁进行弹性分 分析， 其目的主要 要是探讨 ABA AQUS 中不同预 预应力模拟方 方法 之间的区别 别。混凝土取为 C30，按规 规范取其弹性 性模 量 3.0 10 /mm ；梁顶缘和 m 和底缘分别设 设置 四根直径为 为20 的二级 级钢筋，并在 在沿梁高方向 向左 右两边各设 设置两根直径 径12 的二级 级钢筋，箍筋 筋采 用直径10 的一级钢 钢筋，取钢筋 筋的弹性模量 量为 2.1 10 / ，预应力钢筋 筋以直线的形 形式 布置在梁底 底缘，共设置 4 根，直径 径为15.2 。在 。 梁的跨中施 施加1200 的 的集中荷载， ，并对预应力 力钢 1 筋施加600 的预应力 均施加在 力， 在荷载步 中，图 1 示 示处了简支梁 梁的有限元模型。 对于混 混凝土采用实体单元 C3D D20R，普通钢 钢筋 采用桁架单 单元 T3D2；边 边界条件为： ：约束梁一端 端的 平动自由度 度，而另一端只约束竖向自由度。
（b） 钢筋 筋骨架 图 1 有限元模型?
降温法其工 工作机理就是 是通过设置 材料的线膨 膨 胀系 系数，并对索 索单元进行降 降温，从而达 达到施加预应 应 力的 的目的，所施 施加的温度荷 荷载可以通过 过下式求得： 1 式中： — —施加的温度； —预应力钢绞 — 绞线的弹性摸 摸量； —材料的线膨胀系数； — —预应力钢绞 — 绞线的面积； —预加力的大 — 大小。 （b）初始 始应力法 初始应力法 法与降温法相 相类似，只是 是形式上有所 所 不同 同，其实质是 是一样的，初 初始应力法即 即在预应力钢 钢 筋上 上施加初始预 预应力。 （c）rebar 施加初始应 r 应力法 rebar 施加 加初始应力法 法与上述两种 种方法不同， 钢筋 筋不是由杆单 单元模拟，而 而是通过创建 建一个具有钢 钢 筋属 属性的几何面 面，网格划分 分时选取几何 何面的单元种 种 类为 surface， 为 然后通过关键 Initial?co 然 键词 onditions 及
3?ABAQUS 中常用的预应力模拟 S 拟方法?
ABAQU US通用有限元 元软件中，常 常用的有限元 元模 拟方法有： 降温法、 始应力法以及 初始 及rebar施加初 初始 2 。 应力法 降温法 （a）降
陈宜言，深圳 圳市市政设计研 研究院有限公司 司，教授级高工 工 地址： 深圳市 市笋岗西路 300 号市政设计大 07 大厦附楼， 8029 518 电话：83324 4956 22
The finite element analysis Finite element method, the solving area is regarded as made up of many small in the node connected unit (a domain), the model gives the fundamental equation of sharding (sub-domain) approximation solution, due to the unit (a domain) can be divided into various shapes and sizes of different size, so it can well adapt to the complex geometry, complex material properties and complicated boundary conditions Finite element model: is it real system idealized mathematical abstractions. Is composed of some simple shapes of unit, unit connection through the node, and under a certain load. Finite element analysis: is the use of mathematical approximation method for real physical systems (geometry and loading conditions were simulated. And by using simple and interacting elements, namely unit, can use a limited number of unknown variables to approaching infinite unknown quantity of the real system. Linear elastic finite element method is a ideal elastic body as the research object, considering the deformation based on small deformation assumption of. In this kind of problem, the stress and strain of the material is linear relationship, meet the generalized hooke's law; Stress and strain is linear, linear elastic problem boils down to solving linear equations, so only need less computation time. If the efficient method of solving algebraic equations can also help reduce the duration of finite element analysis. Linear elastic finite element generally includes linear elastic statics analysis and linear elastic dynamics analysis from two aspects. The difference between the nonlinear problem and linear elastic problems: 1) nonlinear equation is nonlinear, and iteratively solving of general; 2) the nonlinear problem can't use superposition principle; 3) nonlinear problem is not there is always solution, sometimes even no solution. Finite element to solve the nonlinear problem can be divided into the following three categories: 1) material nonlinear problems of stress and strain is nonlinear, but the stress and strain is very small, a linear relationship between strain and displacement at this time, this kind of problem belongs to the material nonlinear problems. Due to theoretically also cannot provide the constitutive relation can be accepted, so, general nonlinear relations between stress and strain of the material based on the test data, sometimes, to simulate the nonlinear material properties available mathematical model though these models always have their limitations. More important material nonlinear problems in engineering practice are: nonlinear elastic (including piecewise linear elastic, elastic-plastic and viscoplastic, creep, etc. 2) geometric nonlinear geometric nonlinear problems are caused due to the nonlinear relationship between displacement. When the object the displacement is larger, the strain and displacement relationship is nonlinear relationship. Research on this kind of problem Is assumes that the material of stress and strain is linear relationship. It consists