1、6-16 应用梁单元应用梁单元可以模拟结构,该结构一个方向的尺度(长度)是明显地大于其它两个方向的尺度,并且沿长度方向的应力是最重要的。梁理论的基本假设是由变量可以完全地确定结构的变形,这些变量是沿着结构长度方向位置的函数。为了应用梁理论产生可接受的结果,横截面的尺度必须小于结构典型轴向尺度的 1/10。下面是典型轴向尺度的例子: 支承点之间的距离, 横截面发生显著变化部分之间的距离, 所关注的最高阶振型的波长。ABAQUS 梁单元的假设是在变形中垂直于梁轴线的平截面保持为平面。不要误解所谓横截面的尺度必须小于典型 单元 长度的 1/10 的提法。高度精细的网格中可能包含梁单元,其长度小于其横
2、截面尺寸,尽管一般不建议这样做,在这种情况下实体单元可能更适合。6.1 梁横截面几何形状定义梁横截面的轮廓可以应用三种方法之一:从 ABAQUS 提供的横截面库中选择和指定梁横截面的形状和尺度;应用截面工程性质来定义一个一般性的梁轮廓,诸如面积和惯性矩;或者利用特殊二维单元的一个网格,由数值计算得到它的几何量,称为划分网格的梁横截面(meshed beam cross-section) 。ABAQUS 提供了各种常用的横截面形状,如图 6-1 所示,你必须选择并定义梁的几何轮廓。利用任意横截面(arbitrary cross-section)的定义,也可以定义几乎任意的薄壁横截面。关于在 AB
3、AQUS 中应用梁横截面的详细讨论,请参阅 ABAQUS 分析用户手册的第 15.3.9 节“Beam cross-section library” 。如果应用了在 ABAQUS 库中建立横截面的方法之一, ABAQUS /CAE 会提示所需要的横截面尺寸,每一种类型的横截面是不同的。当梁的轮廓与梁的截面特性相关时,可以指定或者是在分析过程中计算截面的工程性质,或者是让 ABAQUS 预先计算它们(在分析开始时) 。选择前者可以应用于线性或者非线性的材料行为(例如,截面刚度6-2因非弹性屈服而改变) ;选择后者虽然是计算效率高,但是只适用于线弹性材料行为。图 6-1 梁横截面相应的,你可以提供
4、截面的工程性质(面积、惯性矩和扭转常数) ,以代替提供横截面尺寸。材料行为可以是线性或者非线性。这样,你可以将梁的几何和材料特性组合起来定义它对载荷的响应,这些响应可能是线性或者是非线性。关于进一步详细的内容,请参阅 ABAQUS 分析用户手册的第 15.3.7 节“Using the *BEAM GENERAL SECTION option to define the section behavior”。划分网格的梁横截面允许表现包括多种材料和复杂几何形状的梁横截面。在ABAQUS 分析用户手册的第 7.13.1 节“Meshed beam cross-section”中进一步讨论了这种类型
5、的梁轮廓。6.1.1 截面点(section points)当应用在 ABAQUS 横截面库的建立梁轮廓的方式来定义梁横截面,并要求在分析过程中计算横截面的工程性质时,在通过分布于梁横截面上的一组截面点上,ABAQUS 计算梁单元的响应。截面点的数目及其位置,详见 ABAQUS 分析用户手册第 15.3.9 节“Beam cross-section library” 。单元的输出变量可在任何一个截面点上输出,如应力和应变;然而,仅在几个选定的截面点上提供了默认的输出。如在 ABAQUS 分6-3析用户手册的第 15.3.9 节“Beam cross-section library”所描述的。对
6、于矩形横截面,所有的截面点如图 6-2 所示。对于该横截面,在点 1、5、21 和 25 上提供了默认的输出。在图 6-2 中所示的梁单元中总共使用了 50 个截面点(两个积分点,每个积分点上有 25个截面点)来计算刚度。图 6-2 在 B32 矩形梁单元中的积分点和默认截面点当要求在分析前计算梁截面的性质时,ABAQUS 就不在截面点上计算梁的响应,而是应用截面的工程性质确定截面的响应。因此,ABAQUS 仅应用截面点作为输出的位置,而你需要指定出理想输出所在的截面点。6.1.2 横截面方向用户必须在整体笛卡儿空间中定义梁横截面的方向。从单元的第一节点到下一个节点的矢量被定义为沿着梁单元的局
7、部切线 t,梁的横截面垂直于这个局部切线矢量。矢量 n1 和 n2 代表了局部(1-2)梁截面轴。这三个矢量 t、n 1、n 2 构成了局部、右手法则的坐标系(见图 6-3) 。6-4图 6-3 梁单元切线矢量 t ,梁截面轴 n1和 n2的取向对于二维梁单元,n 1 的方向总是(0.0, 0.0, -1.0 ) 。对于三维梁单元,有几种方法来定义局部梁截面轴的方向。第一种方法是在定义单元的数据行中指定一个附加的节点(这种方法需要对由 ABAQUS/CAE 产生的输入文件(.inp )进行人工编辑) 。应用从梁单元的第一个节点到这个附加节点的矢量v(见图 6-3) ,作为初始的近似 n1 方向
8、。然后,ABAQUS 定义梁的 n2 方向为 tv。在n2 确定后,ABAQUS 定义实际的 n1 方向为 n2t。上述过程确保了局部切线与局部梁截面轴构成了一个正交系。另一种方法是当在 ABAQUS/CAE 中定义梁截面特性时,可以给定一个近似的 n1方向,然后 ABAQUS 应用上面描述的过程计算实际的梁截面轴。如果你不但指定了一个附加的节点,而且又给出一个近似的 n1 方向,将优先采用附加节点的方法。如果你没有提供近似的 n1 方向,ABAQUS 将从原点到点(0.0,0.0,-1.0 )的矢量作为默认的n1 方向。这里有两种办法可以用来覆盖由 ABAQUS 定义的 n2 方向;这两种办
9、法都要求人工编辑输入文件。一种是给出 n2 矢量的分量作为节点坐标的第 4、第 5 和第 6 个数据值;另一种是使用*NORMAL 选项直接地指定法线方向(添加该选项可以通过ABAQUS/CAE 中的 Keywords Editor(关键词编辑器) ) 。如果同时采用了这两种办法,后者优先。ABAQUS 再由 n2t 定义 n1 方向。用户提供的 n2 方向不必垂直于梁单元的切线 t。当你提供了 n2 方向,局部梁单元切线 t 将被重新定义作为 n1n2 的叉积。在这种情况下,很有可能重新定义的局部梁切线 t 与从第一节点到第二节点的矢量所定义的梁轴线不一致。如果 n2 方向与垂直于单元轴线的
10、平面的夹角超过了 20,ABAQUS 将在数据文件(.dat)中发出一个警告6-5信息。在第 6.4 节“Example: cargo crane”展示的例子中,将说明如何用 ABAQUS/CAE赋值梁横截面的方向。6.1.3 梁单元曲率梁单元的曲率是基于梁的 n2 方向相对于梁轴的取向。如果 n2 方向不与梁轴正交(即,梁轴的方向不与切向 t 一致) ,则认为梁单元有初始弯曲。由于曲梁的行为与直梁的行为不同,用户必须经常检查模型以确保应用了正确的法线,进而正确的曲率。对于梁和壳体,ABAQUS 使用同样的算法来确定几个单元共享节点的法线。在ABAQUS 分析用户手册第 15.3.4 节“Be
11、am element cross-section orientation”中给出了描述。如果你打算模拟曲梁结构,你可能需要使用在前面描述的直接定义 n2 方向的两种方法之一,它允许你更好地控制对曲率进行模拟。即使你打算模拟由直梁组成的结构,由于在共享节点处被平均化的法线,也可能要引入曲率。如前面所解释的,通过直接定义梁的法线,你可以矫正这个问题。6.1.4 梁截面的节点偏置当应用梁单元作为壳模型的加强件时,使梁和壳单元应用相同的节点是很方便的。壳单元的节点是位于壳的中面上,而梁单元的节点是位于梁的横截面上某点。因此,如果壳和梁单元使用相同的节点,壳与梁加强件将会重叠,除非梁横截面是偏置于节点位
12、置(见图 6-4) 。图 6-4 梁作为壳单元的加强部件:( a)梁截面无偏置 (b)梁截面有偏置6-6采用工字型、梯型和任意多边形的梁截面形式,可能要将该截面几何形状定位在与截面的局部坐标系的原点(原点位于单元节点处)具有一定距离的位置上。由于很容易使采用这些横截面的梁偏离它们的节点,可以应用它们作为如图 6-4(b)所示的加强件(如果加强件的翼缘或腹板的屈曲是很重要的,则应该采用壳单元来模拟加强件) 。图 6-5 所示的工字型梁附着在一个 1.2 单位厚的壳上。通过定义梁的节点从 I-截面的底部的偏移量,梁截面的定位可以如图所示。在这种情况下,偏移量为 0.6,亦即壳厚度的一半。壳截面厚度
13、 1 . 2图 6-5 工字型梁用作壳单元的加强件你也可以指定形心和剪切中心的位置;这些位置也可以从梁的节点偏置,从而使你很容易地模拟加强件。另外也可以分别定义梁节点和壳节点,并在两个节点之间采用一个刚性梁的约束连接梁和壳。关于进一步详细的内容,请参阅 ABAQUS 分析用户手册第 20.2.1 节“Linear Constraint equations”。6.2 计算公式和积分在 ABAQUS 中的所有梁单元都是梁柱类单元,这意味着它们可以产生轴向、弯曲和扭转变形。Timoshenko 梁单元还考虑了横向剪切变形的影响。6-76.2.1 剪切变形线性单元(B21 和 B31)和二次单元(B2
14、2 和 B32)是考虑剪切变形的Timoshenko 梁单元;因此,它们既适用于模拟剪切变形起重要作用的深梁又适用于模拟剪切变形不太重要的细长梁。这些单元横截面的特性与厚壳单元横截面的特性相同,如图 6-6(b) 所示,在第 5.2 节 “壳体公式厚壳或薄壳” 给出了讨论。图 6-6 梁的横截面特性 (a)细长梁 ( b)深梁ABAQUS 假设这些梁单元的横向剪切刚度为线弹性和常数。另外,建立了这些梁的数学公式,因此它们的横截面面积可以做为轴向变形的函数而变化,仅在几何非线性模拟中考虑它的影响,此时截面的泊松比具有非零值(详见第 8 章“非线性” ) 。只要梁的横截面尺寸是小于结构的典型轴向尺
15、寸的 1/10(这通常被考虑是梁理论的适用性的界限) ,这些单元就可以提供有用的结果;如果梁的横截面在弯曲变形时不能保持为平面,梁理论是不适合模拟这种变形的。在 ABAQUS/Standard 中的三次单元,称为 Euler-Bernoulli 梁单元(B23 和 B33) ,它们不能模拟剪切变形。这些单元的横截面在变形过程中与梁的轴线保持垂直(见图6-6(a)) 。因此,应用三次梁单元模拟相对细长构件的结构更为有效。由于三次单元可以模拟沿单元长度方向位移的三阶变量,所以对于静态分析,常常可用一个三次单元模拟一个结构构件,而对于动态分析,也只采用很少数量的单元。这些单元假设剪切变形是可以忽略的
16、。一般情况下,如果横截面尺寸是小于结构的典型轴向尺寸的6-81/15,这个假设就是有效的。6.2.2 扭转响应翘曲结构构件经常承受扭矩,几乎所有的三维框架结构都会发生这种情况。在一个构件中引起弯曲的载荷,可能在另一个构件中引起了扭转,如图 6-7 所示。图 6-7 框架结构中的扭转梁对扭转的响应依赖于它的横截面形状。一般说来,梁的扭转会使横截面产生翘曲或非均匀的离面位移。ABAQUS 仅对三维单元考虑扭转和翘曲的影响。在翘曲计算中假设翘曲位移是小量。在扭转时,以下横截面的行为是不同的:实心横截面;闭口薄壁横截面;和开口薄壁横截面。实心横截面在扭转作用下,非圆型的实心横截面不再保持平面;而是发生
17、翘曲。ABAQUS 应用 St.Venant 翘曲理论在横截面上每一个截面点处计算由翘曲引起的剪切应变的分量。这种实心横截面的翘曲被认为是无约束的,所产生的轴向应力可以忽略不计(翘曲约束仅仅影响非常靠近约束端处的结果) 。实心横截面梁的扭转刚度取决于材料的剪切模量 G 和梁截面的扭转常数 J。扭转常数取决于梁横截面的形状和翘曲特征。对于在横截面上产生较大非弹性变形的扭转载荷,应用这种方法不能够得到精确的模拟。6-9闭口薄壁横截面闭口薄壁非圆型横截面(箱型或六边型)的梁具有明显的抗扭刚度,因此,其性质与实心横截面梁类似。ABAQUS 假设在这些横截面上的翘曲也是无约束的。横截面的薄壁性质允许 A
18、BAQUS 考虑剪应变沿壁厚是一个常数。当壁厚是典型梁横截面尺寸的 1/10 时,一般的薄壁假设是有效的。关于薄壁横截面的典型横截面尺寸的例子包括: 管截面的尺寸。 箱型截面的边长。 任意形状截面的典型边长。开口薄壁横截面当翘曲是无约束时,开口薄壁横截面在扭转中是非常柔性的,而这种结构抗扭刚度的主要来源是对于轴向翘曲应变的约束。约束开口薄壁梁的翘曲会引起轴向应力,该应力又会影响梁对其它类型载荷的响应。ABAQUS/Standard 具有剪切变形梁单元,B31OS 和 B32OS,它们包括了在开口薄壁横截面中的翘曲影响。当模拟采用开口薄壁横截面的结构在承受显著的扭转载荷时,例如管道(定义为任意多
19、边形截面)或者工字型截面,必须使用这些单元。翘曲函数翘曲引起了整个梁横截面的轴向变形,由截面的翘曲函数定义了翘曲的变化。在开口截面梁单元中,采用一个附加的自由度 7 来处理这个函数的量值。约束住这个自由度可以使被约束的节点不会发生翘曲。在每个构件分支上的翘曲量可以不同,在框架结构中开口截面梁之间的连接点处,一般每个构件分支应该使用各自不同的节点(见图 6-8) 。然而,如果连接方式的设计已经防止了翘曲,则所有的构件应该共享同一个节点,并必须约束住翘曲的自由度。6-10图 6-8 连接开口截面梁当剪力没有通过梁的剪切中心作用时会产生扭转,扭转力矩等于剪力乘以它到剪切中心的偏心距。对于开口薄壁梁截
20、面,其形心和剪切中心常常并不重合(见图 6-9) 。如果节点不是位于横截面的剪切中心,在载荷作用下截面可能扭转。图 6-9 关于一些梁横截面的剪切中心 s 和形心 c 的近似位置6.3 选择梁单元 在任何包含接触的模拟中,应该使用一阶、剪切变形梁单元(B21,B31 )。 如果横向剪切变形是非常重要的,采用 Timoshenko(二阶)梁单元(B22,B32 )。 如果结构是或者非常刚硬或者非常柔软,在几何非线性模拟中,应当使用ABAQUS/Standard 中的杂交梁单元(B21H ,B32H 等)。6-11 在 ABAQUS/Standard 中的 Euler-Bernoulli 三次梁单
21、元(B23,B33 )模拟承受分布载荷作用的梁有很高的精度,例如动态振动分析。 在 ABAQUS/Standard 中,模拟开口薄壁横截面的结构应该采用那些应用了开口横截面翘曲理论的梁单元(B31OS, B32OS)。6.4 例题:货物吊车一个轻型的货物吊车,如图 6-10 所示,要求你确定当它承受 10kN 载荷时的静挠度。你也要标识结构中的关键部件和节点,即它们有最大的应力和载荷。由于这是一个静态分析,你将应用 ABAQUS/Standard 分析这个货物吊车。图 6-10 轻型货物吊车的草图吊车由两榀桁架结构组成,通过交叉支撑连接在一起。每榀桁架结构的两个主要构件是箱型截面钢梁(箱型横截
22、面) 。每榀桁架结构由内部支撑加固,内部支撑焊接在主要构件上。连接两榀桁架结构的交叉支撑通过螺栓连接在桁架结构上,这些连接不能传递弯矩(如果存在弯矩的话) ,因此,将它们作为铰节点处理。内部支撑和交叉支撑均采用箱型横截面钢梁,采用远小于桁架结构主要构件的横截面。两榀桁架结构在它们的端点(在点 E)连接,这种连接方式允许它们各自独立地沿 3 方向移动和所有的转动,而约束它们在 1 方向和 2 方向的位移相等。吊车在点 A、B、C 和 D 牢固地焊6-12接在巨大的结构上。吊车的尺寸如图 6-11 所示。在下面的图中,桁架 A 是包括构件AE、BE 和它们的内部支撑的结构;桁架 B 是包括构件 C
23、E、DE 和它们的内部支撑的结构。图 6-11 货物吊车的尺寸(m)在吊车的主要构件中,典型横截面的尺寸与总长度的比值远小于 1/15。在内部支撑应用的最短构件中,这个比值近似为 1/15。因此,应用梁单元模拟吊车是合理的。6.4.1 前处理应用 ABAQUS/CAE 创建模型现在讨论如何应用 ABAQUS/CAE 创建这个模型。在本手册的在线文档第 A.4 节“Cargo crane”提供了输入文件。当通过 ABAQUS/CAE 运行这个输入文件时,将创建关于该问题的完整的分析模型。根据下面给出的指导如果你遇到困难,或者如果你希望检查你的工作,则可以运行这个输入文件。在附录 A“Exampl
24、e Files”中,给出了如何提取和运行输入文件的指导。如果你没有进入 ABAQUS/CAE 或者其它的前处理器,可以人工创建关于这个问6-13题的输入文件,关于这方面的讨论,见 Getting Started with ABAQUS/Standard:Keywords Version,第 6.4 节“Example:Cargo crane” 。创建部件在吊车的主要构件和它的内部支撑之间是焊接节点,从模型的一个区域到邻近的区域,焊接节点提供了平移和转动的完全的连续性。因此,在模型中每一个焊接节点仅需要一个几何实体(即顶点) ,可应用单一部件代表主要构件和内部支撑。为了方便,两个桁架结构将作为一
25、个部件处理。将交叉支撑连接到桁架结构是螺栓节点,它们在桁架结构端点的连接区别于焊接节点的连接。由于这些节点对于所有的自由度不提供完全的连续性,所以在连接处需要分别给出各自的节点。由于需要明确的几何实体模拟螺钉,这样,交叉支撑需要作为独立的实体。在分离的节点之间需要定义适当的约束。我们从讨论如何定义桁架几何形体的技术开始。由于两个桁架结构完全一致,仅应用单一桁架结构的几何形状来定义部件的基本特征是足够的。该桁架的几何构图可以被存储,然后应用它再增加第二个桁架结构到部件的定义中。图 6-11 显示的尺寸是相对于在图中的笛卡儿坐标系给出的,然而,基本特征是画在部件的局部 1-2 平面上;因此,当绘制
26、桁架时所指定的尺寸需要做相应的调整。一旦所有的部件装配在一个公共坐标系中,它们可以根据需要进行旋转和重新定位,这样,结构是与整体坐标系一致的,如图 6-11 所示。定义单一桁架的几何形状:1 首先创建一个三维、可变形的平面线框,设置近似的部件尺寸为 15.0,并命名部件为 Truss。2 应用 Create Lines(创建线): Connected(连接)工具,创建两条几何线代表桁架的主要构件。在图中标注尺寸;并用 Edit Dimension Value(编辑尺寸值)工具 编辑尺寸以给出桁架的准确水平跨度,如图 6-12 所示。提示:应用 Sketcher Options 修改尺寸文本的格
27、式。3生成五个独立点,如图 6-13 所示。对每个点创建和编辑尺寸标注,如图所示,然后通过每个点创建一条竖直辅助线。在主要构件上确定了在辅助线与两个主要构件之间的交叉点,在这些点上将内部支撑焊接到桁架上。6-14图 6-12 桁架的主要构件图 6-13 用于定位辅助几何的点6-154采用在辅助线和几何线(即代表桁架主要构件的线)之间预选的点,在焊接的位置创建独立点(isolated point)。此外,在两条几何线的端点创建独立点。5删除几何线,并应用一系列的连接线重新定义桁架的几何形状。例如,从位于结构左下角的点开始,以逆时针的方式依次连接相邻点,可以定义整个桁架的几何形状。最终的图形如图
28、6-14 所示。图 6-14 单榀桁架结构的最终图形6应用 Save Sketch As(保存草图)工具 ,将图形保存为 Truss。7点击 Done 退出绘图环境并保存部件的基本特征。另一个桁架将作为一个平面线框特性加入。当加入一个平面特性时,不仅需要指定一个构图平面,而且还要指定它的方位。我们将应用基准面(datum plane)定义该平面;应用一个基准轴(datum axis)定义平面的方位。然后,将桁架草图投影到这个平面上。定义第二个桁架结构的几何形体:1从桁架的端点应用偏置(offset)定义三个数据点,如图 6-15 所示。6-16图 6-15 基准点、基准面和基准轴从母体端点的偏
29、置已在图中标出。此外,如图中标示的定义了第四个基准点。可以根据需要旋转视图以观察基准点。前三个基准点用来定义基准面;第四个基准点用来定义基准轴。2应用 Create Datum Plane(创建基准面): 3 Points(3 点法)工具来定义一个基准面,和 Create Datum Axis(创建基准轴): 2 points(2 点法)工具来定义一个基准轴,如图 6-15 所示。3应用 Create Wire(创建线) : Planar(平面法)工具给部件增加一个特性。选取基准面作为绘图平面,和选取基准轴作为边界,该边界为竖向显示在图的右侧。4应用 Add Sketch(添加草图)工具重新获
30、得桁架草图。通过选择新桁架端部的顶点作为平移矢量的起点,和标记为 E的基准点作为平移矢量的终点,平移(translate)草图,如图 6-15 所示。如果需要,可以放大和旋转图形以便于选取。5在提示区,点击 Done 退出绘图环境。最终的桁架部件如图 6-16 所示。6-17图 6-16 最终的桁架结构的几何图形;放大标识的点表示铰接点的位置回顾交叉支撑必须作为分离的部件处理,才能正确地代表在支撑与桁架之间的铰接关系。然而,绘制交叉支撑最简单的办法是直接在桁架连接点的位置之间创建线框特征。因此,我们将采用如下的方法创建交叉支撑部件:首先,创建一个桁架部件的复制件,在其上添加代表交叉支撑的线段(
31、我们不能采用这一新的部件,这是因为在铰接点处是共享的,这样不能代表一个铰接点) ;然后,我们将使用在装配(Assembly)模块中的切割性能,在含有交叉支撑的桁架和没有交叉支撑的桁架之间进行布尔(Boolean)切割,保留交叉支撑几何部分作为一个明确的部件。其过程详细地描述如下。创建交叉支撑的几何形体:1 从主菜单栏中,选择 PartCopyTruss;在 Part Copy(复制部件)对话框中,命名新部件为 Truss-all,并点击 OK。2 铰接位置在图 6-16 中以高亮度显示。利用 Creat Wire: Poly Line(多线)工具添加交叉支撑几何到新部件中,如图 6-17 所示
32、(在该图中的点对应于在图 6-16 中标记的点;在图 6-17 中隐去了桁架的可视性)。采用如下的坐标指定类似的视图:Viewpoint(1.19, 5.18, 7.89),Up Vector(-0.40, 0.76, -0.51)。提示:如果在连接交叉支撑几何时出现错误,你可以应用 Delete Feature(删6-18除特征)工具 删除线段;你不能恢复被删除的特征。图 6-17 交叉支撑几何形体3进入 Assembly 模块,创建每个部件的实体(Truss 和 Truss-all)。4从主菜单栏中,选择 InstanceMerge/Cut。在 Merge/Cut Instance(合并/
33、切割实体)对话框中,命名新部件为 Cross brace,在 Operations(操作)域中选择 Cut geometry(切割几何体),并点击 Continue。5从 Instance List(实体列表)中,选择 Truss-all-1 作为被切割的实体和Truss-1 作为将用于切割的实体。在切割完成后,创建了一个名为 Cross brace,且仅包含交叉支撑几何的新部件。当前的装配模型只包含这个部件的一个实体;而原来的实体被默认删除了。由于在模型的装配中,我们需要使用原来的桁架,打开 Feature Manager(特征管理器)恢复名为 Truss-1 的部件实体。现在,我们回到 P
34、roperty 模块,定义梁截面性质。定义梁截面性质在该分析中,由于假设材料行为是线弹性,从计算的观点考虑采用预先计算梁截面性质的方法将大幅提高计算效率。假定桁架和支撑都是由中等强度的钢材制造,E = 200.0 109 Pa,v = 0.25,G = 80.0 109 Pa。在该结构中所有的梁都是箱型横截面。箱型截面如图 6-18 所示。在吊车中两榀桁架的主要构件的尺寸如图 6-18 所示。6-19支撑构件的梁截面的尺寸如图 6-19 所示。图 6-18 主要构件的横截面几何形状和尺寸(m)图 6-19 内部和交叉支撑构件的横截面几何形状和尺寸(m)6-20定义梁截面性质:1 在 Prope
35、rty 模块中,创建两个箱型轮廓(Profile):一个是桁架结构的主要构件,另一个是内部和交叉支撑;将两个轮廓分别命名为 MainBoxProfile 和BraceBoxProfile。采用在图 6-18 和 6-19 所示的尺寸完成轮廓的定义。2为桁架结构的主要构件和内部和交叉支撑各创建一个梁截面,并分别命名截面为MainMemberSection 和 BracingSection。a. 对于两个截面的定义,在分析前指定截面的积分方式。当选择了这种类型的截面积分,材料性质定义作为截面定义的组成部分,而不需要另外给出材料的定义。b. 选择 MainBoxProfile 作为主要构件的截面定义
36、,和选择BraceBoxProfile 作为支撑截面的定义。c. 点击 Linear properties(线性性质),在 Beam Linear Behavior(梁线性行为)对话框中的相应文本区域中(如前面描述的)输入杨氏模量和剪切模量。d. 在 Edit Beam Section(编辑梁截面)对话框的相应文本区域中输入泊松比。3将 MainMemberSection 赋予几何区域代表桁架的主要构件,并将BracingSection 赋予区域代表内部和交叉支撑构件。应用位于工具栏下面的Part 列表恢复每一个部件。由于不再需要 Truss-all 部件,可以忽略它。定义梁截面方向主要构件的
37、梁截面轴定位为:梁的 1 轴是正交于桁架结构的平面,如在正视图(图 6-11)中所示,而梁的 2 轴是正交于该平面中的单元。对于内部桁架支撑和与之相应的桁架结构的主要构件,其近似的 n1 矢量是相同的。在它的局部坐标系中,Truss 部件的方向如图 6-20 所示。6-21图 6-20 桁架在它的局部坐标系中的方向从主菜单栏中,选择 AssignBeam Section Orientation 为每个桁架结构指定一个近似的 n1 矢量。如前面所述,该矢量的方向必须正交于桁架的平面。因此,对于平行于部件局部 1-2 平面的桁架(桁架 B) ,近似的 n1 =(0.0,0.0,1.0) ;而另一个
38、桁架结构(桁架 A) ,其近似的 n1 =(0.2222, 0.0, 0.975) 。从主菜单栏中,选择 AssignTangent 指定梁的切线方向。当必要时翻转切线方向,这样,它们显示的结果如图 6-21 所示。图 6-21 梁的切线方向所有交叉支撑和每榀桁架结构内部的支撑都具有相同的梁截面几何,而它们的梁截面轴的方向却各不相同。由于方形交叉支承构件主要是承受轴向载荷,它们的变形对横截面的取向并不敏感;因此,我们可以采用某些假定,以便使指定6-22交叉支撑的方向多少更容易些。所有的梁法线(n 2 矢量)必须是近似地位于货物吊车俯视图的平面内(见图 6-19) ,这个平面相对于整体 1-3
39、平面只有轻微的倾斜。定义这个方向的一个简单办法就是提供一个正交于这个平面的近似的 n1 矢量。该矢量应该几乎是平行于整体的 2 方向。因此,对于交叉支撑,指定n1(0.0,1.0,0.0) ,以使它与部件的(后面也将看到,整体的)y 轴一致。梁的法线在这个模型中,如果你提供的数据仅定义了近似的 n1 矢量的方向,将会引起模拟的误差。如果不覆盖的话,梁法线的平均化方式(见第 6.1.3 节“梁截面曲率”)将引起 ABAQUS 对于货物吊车模型使用不正确的几何形状。为了看到这点,你可以应用 Visualization 模块显示梁截面轴和梁切线矢量(见第 6.4.2 节“后处理” ) 。如果没有对梁
40、法线的方向给出进一步的修正,在吊车模型中的法线将在Visualization 模块中显示是正确的;可是事实上,它们将有轻微的偏差。图 6-22 显示了桁架结构的几何形状。从该图可以看出,对于吊车模型,其正确的几何形状要求在顶点 V1 有三个独立的梁法线:R1 区及 R2 区各一个,R3 区和 R4 区共用一个。应用 ABAQUS 关于平均法线的逻辑,显然在 R2 区中顶点 V1的梁法线将与在该点的相邻区的法线进行平均。在这个例子中,平均化逻辑的重要性在于当法线对参考法线的夹角小于 20时,将该法线与参考法线平均化从而定义一条新的参考法线。假设在该点的初始参考法线是对于 R3 区和 R4 区的法
41、线,由于在顶点 V1 处 R2 区的法线与初始参考法线的夹角小于 20,所以它与初始参考法线取平均化,在该点处定义新的参考法线。另一方面,由于在顶点 V1 处 R1区的法线与初始参考法线的夹角大约是 30,因此它有一个独立的法线。图 6-22 需要指定梁法线的位置6-23这个不正确的平均法线意味着在 R2、R3 和 R4 区中创建的单元,并共享在V1 处创建的节点,这些单元具有围绕从单元的一端到另一端的梁轴线产生扭曲的截面几何形状,而这并不是所希望的几何形体。所以在相邻区域法线的夹角小于20的位置,你必须明确地指定法线方向,这样才能避免 ABAQUS 应用它的平均算法。在这个例子中,在吊车的两
42、侧桁架中的相应区域都需要采取这种办法。 在桁架结构的顶点 V2 处的法线也存在这个问题,因为相连这个顶点的两个区域之间的夹角也小于 20。由于我们模拟的是直梁,在每根梁的两个端点的法线是常数,通过直接地指定梁的法线方向,可以修正这一问题。如前所述,在吊车的两侧桁架中的相应区域都需要采取这种办法。目前,在 ABAQUS/CAE 中指定梁法线方向的唯一方法是应用 Keywords Editor(关键词编辑器) 。Keywords Editor 是一个专门的文本编辑器,在提交输入文件进行分析前,它允许你修改由 ABAQUS/CAE 生成的 ABAQUS 输入文件。这样,它允许你添加 ABAQUS/S
43、tandard 或者 ABAQUS/Explicit 的功能,这些功能在目前的 ABAQUS/CAE 版本还不能够得到支持。关于 Keywords Editor 的更多信息,请参阅 ABAQUS/CAE 用户手册的第 13.8.1 节, “Adding unsupported keywords to your ABAQUS/CAE model”。由于 Keywords Editor 只能在 Job 模块中调用,所以你只能在此指定梁的法线方向。创建一个装配件我们现在关注如何组装模型。这包括旋转和平移前面已经创建的部件实体,使它们在整体笛卡儿坐标下组装成一个整体。如图 6-11 所示。定位吊车装配
44、件:1 从主菜单栏中,选择 InstanceRotate,将桁架部件实体绕着由 C 点和 D 点定义的轴旋转 6.4188(定义在图 6-10 中)。该轴平行于整体 y 轴。判断旋转的方位取决于轴的方向,而轴的方向取决于首选了哪个顶点;例如,如果选择 C 点作为定义旋转轴的起点,转动必须指定作为正值(即对应于 y 轴负方向)。ABAQUS/CAE 显示一个临时的效果图,演示对所选择部件实体将要施加的旋转;在提示区点击 OK,以接受部件实体的新位置。2 对交叉支撑部件的实体重复上述步骤,确保旋转该实体的轴与对桁架所用的旋转轴一致(即再次应用 C 点和 D 点)。6-243 在 B 点和 D 点之
45、间的中点处创建一个基准点(见图 6-10),然后,从主菜单栏中,选择 InstanceTranslate 平移桁架部件的实体。指定这个基准点作为平移矢量的起点和点(0.0, 0.0, 0.0)作为矢量的终点。你也必须将交叉支撑部件的实体平移相同的量。最简单办法是在交叉支撑上任选一点作为平移矢量的起点,并在桁架中相应的点作为矢量的终点。为了方便,此时将定义一些以后会用到的几何集合。定义一个包含从点 A 到点 D 的几何集合(参见图 6-10 的准确位置) ,并命名集合为 Attach。另外创建位于桁架顶端处的点的集合(在图 6-10 中的位置 E) ,分别命名为 Tip-a 及Tip-b,用 T
46、ip-a 表示与桁架 A 相关的几何集合(见图 6-16) 。最后,参照图 6-16 和图 6-22,为每一个需要指定梁法线的区域创建一个集合。对于桁架 A,为由R2 表示的区域创建一个命名为 Inner-a 的集合,和为由 R5 表示的区域创建一个命名为 Leg-a 集合;为桁架 B,创建相应的集合 Inner-b 和 Leg-b。创建分析步定义和指定输出在 Step 模块,创建一个一般静态( static,general )分析步,命名该步骤为Tip load,并输入下面的分析步描述:Static tip load on crane。为了应用 ABAQUS/CAE 作后处理,将在节点处的位
47、移()和约束反力(RF)以及在单元中的截面力(SF)作为场变量写入输出数据库。定义约束方程在 Interaction(相互作用)模块中指定在节点自由度之间的约束。每个方程的形式是 120,nAuAu式中 Ai 是与自由度 ui 相关的系数。在吊车模型中,将两榀桁架的顶端连接在一起,这样,自由度 1 和 2(在 1和 2 方向的平移)是相等的,而其它的自由度(36)是独立的。我们需要两个线性约束,一个是在两个顶点的自由度 1 相等,而另一个是在两个顶点的自由度2 相等。创建线性方程:1 切换到 Interaction 模块,从主菜单栏中,选择 ConstraintCreate,命名约束为TipC
48、onstraint-1,并指定为方程约束。6-252 在 Edit Constraint(编辑约束)对话框中,在第一行中输入系数(Coefficient )1.0,集合名(Set Name)Tip-a 和自由度(DOF)1。在第二行中输入系数-1.0,集合名 Tip-b 和自由度 1。点击 OK。这样就生定义了自由度 1 的约束方程。注意:在 ABAQUS/CAE 中的文本输入是区分大小写的。3 从主菜单栏中,选择 ConstraintCopy,将 TipConstraint-1 复制到TipConstraint-2。4 选择 ConstraintEditTipConstraint-2,将两行
49、的自由度改为 2。在刚度矩阵中,将消去定义在约束方程中与第一个集合相关的自由度。因此,这个集合将不能在其它的约束方程中出现,而且边界条件也不能施加在消去的自由度上。模拟在桁架和交叉支撑之间的铰接与桁架内部的支撑不同,交叉支撑是用螺栓连接在桁架构件上。你可以假设这些螺栓连接处是不能传递转动和扭转。在这些需要定义约束的位置处定义了两个相同的节点。在 ABAQUS 中可以用多点约束、约束方程或者连接件来定义这样的约束。在本例中,将采用最后一种方法。连接件允许模拟在模型装配件中任意两点之间的连接(或者在装配件中的任意一点与地面之间) 。在 ABAQUS 中包含一个庞大的连接件库。关于每种连接件类型的描述和全部列表,请参阅 ABAQUS 分析用户手册的第 17.1.5 节“Connector element library”。应用 JOIN 连接件模拟螺栓连接。由这种连接件建立的铰连接约束了相等的位移,而转动(如果它们存在)则保持独立。每一个连接件必须提供一个连接件特性以定义它的类型(类似于在单元与截面特性之间的关系) 。因此,首先要定义特性,然后是各个连接件。定义连接件特性:1 从主菜单栏中,选择 ConnectorPropertyCreate。在 Creat