有限元三杆桁架优化分析报告

利用 ANSYS 有限元分析软件对三杆组成的桁架结构进行数值模拟,并根据计算结果,建立优化设计数学模型,在优化处理器指定分析文件, 对三根横截面积为
A1A2A3基本尺寸 B 为变量进行分析对比, 通过数值迭代模拟主要的到如下结论
(1横截面积迭代进行 ANSYS 优化分析时,在分析得到的重量,应力,横截面,三个图中当寻优迭代进行到第 16次主动变量被调整到相同的优化效率时 A1为 1
10
7056
. 4-
⨯A2为 4
10
0000
. 6-
⨯A3为 2
10
3055
. 3-
⨯, 桁架重量取得最小值 130370kg 与初始设计重量 481520.422kg 相比,得到了很大程度的减轻。

符合最优化准则 (2根据计算结果,改进的桁架明显好于其他情况, ansys 软件数值模拟得到最优解,其计算误差很小,完全能满足工程精度要求
ANSYS 程序中进行优化的方法是成功的 , 方法本身收敛速度快 , 精度高 , 稳定性强。

本文使用迭代法得到的最优解都非常接近于或优于所求问题的最优解 , 这表明将迭代法一类的高效优化方法用 APDL 语言嵌套到 AnSYS 程序中来求解优化问题的方法既可行又简便 , 结构优化设计领域具有很好的应用前景。

第9章桁架和梁的有限元分析第1节基本知识一、桁架和梁的有限元分析概要1．桁架杆系的有限元分析概要桁架杆系系统的有限元分析问题是工程中最常见的结构形式之一，常用在建筑的屋顶、机械的机架及各类空间网架结构等多种场合。

桁架结构的特点是，所有杆件仅承受轴向力，所有载荷集中作用于节点上。

由于桁架结构具有自然离散的特点，因此可以将其每一根杆件视为一个单元，各杆件之间的交点视为一个节点。

2．梁的有限元分析概要梁的有限元分析问题也是是工程中最常见的结构形式之一，常用在建筑、机械、汽车、工程机械、冶金等多种场合。

梁结构的特点是，梁的横截面均一致，可承受轴向、切向、弯矩等载荷。

根据梁的特点，等截面的梁在进行有限元分析时，需要定义梁的截面形状和尺寸，用创建的直线代替梁，在划分网格结束后，可以显示其实际形状。

二、桁架和梁的常用单元桁架和梁常用的单元类型和用途见表９-1。

通过对桁架和梁进行有限元分析，可得到其在各个方向的位移、应力并可得到应力、位移动画等结果。

第2节 桁架的有限元分析实例一、案例1——2D 桁架的有限元分析图9-1 人字形屋架的示意图 问题人字形屋架的几何尺寸如图9-1所示。

杆件截面尺寸为0.01m 2，试进行静力分析，对人字形屋架进行静力分析，给出变形图和各点的位移及轴向力、轴力图。

条件人字形屋架两端固定，弹性模量为2.0×1011 N/m 2，泊松比为0.3。

解题过程制定分析方案。

材料弹性材料，结构静力分析，属2D 桁架的静力分析问题，选用Link1单元。

建立坐标系及各节点定义如图9-1所示，边界条件为1点和5点固定，6、7、8点各受1000 N 的力作用。

1．ANSYS 分析开始准备工作（1）清空数据库并开始一个新的分析 选取Utility>Menu>File>Clear & Start New ，弹出Clears database and Start New 对话框，单击OK 按钮，弹出Verify 对话框，单击OK 按钮完成清空数据库。

x
由单元①的刚度方程：
Fj
①
k
① ji
i
①
k
① jj
j
①
k
① ji
2
k
① jj
1
由单元③的刚度方程：
Fj
③
k
③ ji
i
③
k
③ jj
j
③
k
③ ji
3
k
③ jj
1
§2.3 结点平衡与整体刚度矩阵的集成
代入结点1的平衡条件:
k
l
xi
)
(dx j
dxi
)
(
yj
l
yi )
(dy j
dyi )
(dx j dxi ) (dy j dyi )
cos sin
由于杆件的变形产生位移：
ui dxi vi dyi
u j dxj v j dy j
因此，杆件应变为：
dl l
l
(ui
uj)
l
(vi
vj)
杆件轴力为：
(2k1 k2 )v4 P
结构的整体刚度系数
v4
P 2k1
k2
12 3
l2 l1 l1
4 P
N1
N1y
cos
k1v4
cos
k1P
(2k1 k2 ) cos
N2
k2v4
k2P 2k1 k2
位移法求解超静定结构。
§2.1 平面桁架单元的离散
结构的离散化：尽量将结构离散成数量最少的等截面直 杆单元
kki③ ③jii
ki③j
k
③ jj
3 3 3 3
§2.3 结点平衡与整体刚度矩阵的集成

【 关键词】 平面圆钢管桁架； 极限温度； 有限元； 抗火性能
【 中图分类号】 Ｔ ３２３ Ｕ ９．
【 文献标识码】 Ｂ
【 文章编号】 １ １ ６６ （０２ ０ — ０７ ０ ０ — ８４ ２ １）２ ０２ — ３ ０
ＡＮＴ ． ＲＥ ＥＸＰ Ｉ ＦＩ ＥＲⅡ ＥＮＴ ＡＮＤ Ｉ ＴＥ ＥＭ ＥＮＴ Ｆ ＮＩ ＥＬ ＡＮＡＬＹＳ Ｓ Ｉ
Ａｂ ｔ ａ ｔ Ｅ ｐ ｒ ｎ ａ ｅ ｅ ｒ ｈ ｈ ｓ ｂ ｅ ｏ ｄ ｃ ｅ ｏ ｓｕ ｙ ｔ ｅ ｍｅ ｈ ｎ ｃ ｌｂ ｈ ｖ ｏ ｆｔ ｕ ｌ ｓ ｒ ｃ ： ｘ ｅ ｍｅ ｔｌｒ ｓ ａｃ ａ ｅ ｎ ｃ ｎ ｕ ｔ ｄ ｔ ｔ ｄ ｈ ｃ ａ ｉａ ｅ ａ ｉｒ ｏ ｉ ｗｏ ｆ ｌ －
本次火灾试验 的试验装置 详见 图 １ 。反力架 为一 平面钢架 ， 由一个工字钢梁和两个工字钢柱组成 ， 工字
钢梁截面尺寸 为 １０ ｍ ０ ｍ Ｘ ０ ｍ Ｘ ０ ｍ， ０ ０ ｍ ｘ ０ ｍ ２ ｍ ２ ｍ 工 ５
材 的屈服强度和弹性模 量会 降： 温时 的一半 ， 勾室 因此 没 有采取任何防火措施 的钢结 构 ， 在火 灾发 生时极 易 引起建筑物 的倒塌 。大量钢结构建筑 由于火灾造成破
字钢柱 截 面尺 寸为 ５ ０ ｍ×５０ ｍ Ｘ２ ｍ ｍ ０ｍ ０ ｍ ｍ Ｘ２ ｍ。 ０ ０
为保证反力架可 以承受试 验时产 生的水 平反力 ， 工字
【 摘
要】 进行了２ 个足尺钢管桁架的抗火试验， 建立了桁架的有限元模型， 经比较， 有限元与试验结果吻
合 良好 。火灾试验和有 限元分析结果 表明 ： 灾作 用下平 面 圆钢 管桁架 的破 坏主要是 由于 腹杆管壁 局部屈 服引 火 起的 ； 并且 随着荷载 比的增加 ， 面圆钢管Hale Waihona Puke 架的耐火极 限温度逐渐 降低 。 平

Ｃｙ 。 Ｃｙ ＣＺ
ｃｘ Ｚ ｍ
Ｃｙ ＣＹ ＣＺ
ＣＺ ＣＺ
ＣＸ Ｃｙ
ＣＸ ＣＺ
Ｃ ＣＺ
ＣＺ ６６ ×
—
ＣＸ ＣＺ
一
一
其 中 ，足 ］为二 力杆单 元 轴 向 刚度 矩阵 ， ［ 愚
ＥＡ
Ｌ 。
２ 编 程 思 路
２ １ 工 程 实 例 ．
如 图 ２所示 ， 一机 架 由空间桁 架杆组 成 ， 其杆 件单 元 的横 截面 积为 １ ｍ。 由钢制 成 （ ５ｃ ， Ｅ一 ２ ０ＧＰ ） 试 ０ ａ ，
Ｃ ０ ＯＳ ｙ
＿—
—
［ 卅 尼］
—
Ｃ ｘＣ ＳＯ ＯＳＯ Ｏ ｚ
—
Ｃ Ｓ０ ｏ ０ Ｏ ｙｃ ｓ ｚ
。— —
ＣＳ ｚ Ｏ Ｏ
—
—
Ｃ Ｓ０ Ｃ ＳＯ — ＯＳ０ ＯＳＯ Ｏ ｘ Ｏ ｚ —Ｃ ｙＣ ｚ
ＣＳ ｘ Ｏ ０ Ｃ ｘＣ ｒ ＯＳ０ ＯＳＯ Ｃ Ｓ０ Ｏ Ｏ ｘＣ Ｓ
Ｃ Ｓ０ Ｏ ｚ Ｏ ｘＣ Ｓ ０
Ｃ ｙ Ｏ ｚ ＯＳ０ Ｃ Ｓ Ｏ
—
一
ＣＳ ｘ Ｏ Ｏ
。 Ｏ ｘＣ ｙ — —Ｃ Ｓ Ｏ ＯＳ０
—
・—
—
Ｃ ＳＯ Ｏ ｚ Ｏ ｘＣ Ｓ Ｏ Ｃ Ｓ ０ Ｃ ＳＯ Ｏ ｙ Ｏ ｚ
—
Ｃ Ｓｅ Ｏ ｙ Ｏ ｘＣ Ｓ ０
—
—
ＣＸ ｚ
一
ＣＸ， ＣＹ
ＣＸ ＣＺ Ｉ
ＣＸ ＣＹｍ
ＣＸ ＣＺ ＣＸ 。
一
Ｃ ＣＺ
ＣＺ
ＣＺ
。
［ 。 ｒ Ｋ ］＝ － ］ ｋ
一

有限元上机分析报告学院：机械工程专业及班级：机械设计及其自动化08级7班姓名：王浩煜学号：20082798题目编号： 21.题目概况1.1 结构组成和基本数据结构：该结构为一个六根杆组成的桁架结构，其中四根杆组成了直径为800cm的正方形，其他两根杆的两节点为四边形的四个角。

材料：该六根杆截面面积均为100cm2，材料均为Q235，弹性模量为200GPa，对于直径或厚度大于100mm的截面其强度设计值为190Mpa。

载荷：结构的左上和左下角被铰接固定，限制了其在平面内x和y方向的位移，右上角受到大小为2000KN的集中载荷。

结构的整体状况如下图所示：1.2 分析任务该分析的任务是对该结构的静强度进行校核分析以验算该结构否满足强度要求。

2.模型建立2.1 物理模型简化及其分析由于该结构为桁架结构，故认为每根杆件只会沿着轴线进行拉压，而不会发生弯曲和扭转等变形。

结构中每根杆为铰接连接，有集中载荷作用于最上方的杆和最右方杆的铰接点。

2.2单元选择及其分析由于该结构的杆可以认为是只受拉压的杆件，故可以使用LINK180单元，该单元是有着广泛工程应用的杆单元，它可以用来模拟桁架、缆索、连杆、弹簧等等。

这种三维杆单元是杆轴方向的拉压单元，每个节点具有三个自由度：沿节点坐标系X、Y、Z方向的平动。

就像铰接结构一样，不承受弯矩。

输入的数据有：两个节点、横截面面积（AREA）、单位长度的质量(ADDMAS)及材料属性。

输出有：单元节点位移、节点的应力应变等等。

由此可见，LINK180单元适用于该结构的分析。

3.3 模型建立及网格划分（1）启动Ansys软件，选择Preferences→Structural，即将其他非结构菜单过滤掉。

（2）选择单元类型：选择Preprocessor→Element Type→Add/Edit/Delete→Add，在出现的对话框中选择Link→3d finit stn 180，即LINK180，点击“OK”（3）选择实常数：选择Preprocessor→Real Constants→Add/Edit/Delete→Add，在出现的对话框中的Cross-sectional area中输入100，点击“OK”。

平面桁架结构的有限元分析平面桁架结构是一种经常在建筑和工程领域中使用的结构形式。

它由直杆组成，连接在节点上，形成一个稳定的平面结构。

平面桁架结构的设计和分析需要使用有限元分析方法来确定结构的受力状态和稳定性。

本文将介绍平面桁架结构的有限元分析方法，包括模型建立、加载条件、应力和变形分析等。

首先，建立平面桁架结构的有限元模型。

模型应包括杆件和节点两个基本元素。

杆件是结构的主要受力元素，节点是杆件的连接点。

通过连接节点和杆件，可以构建起整个桁架结构。

在有限元模型中，每个节点被赋予一个坐标，每个杆件的长度和截面积也需要定义。

通过这些信息，可以建立结构的有限元模型。

加载条件是进行有限元分析的第二个关键步骤。

加载条件包括结构所承受的外部力和约束条件。

外部力是指作用于结构上的力，包括重力、风力、地震力等。

约束条件是指限制结构自由运动的条件，例如固定节点或滑动支座等。

在有限元分析中，将这些加载条件应用到有限元模型中，以模拟真实结构的受力情况。

然后进行应力和变形分析。

在有限元分析中，结构的应力分布和变形情况可以通过求解有限元方程来得到。

有限元方程是由结构的力平衡和材料的应力-应变关系所组成的方程组。

通过求解有限元方程，可以计算出结构中每个节点的应力和变形情况。

这些结果可以用来评估结构的安全性和稳定性。

在进行有限元分析时，需要注意一些细节。

首先，选择合适的材料模型和参数。

不同的材料具有不同的力学特性，例如弹性模量、屈服强度等。

选择适当的材料模型和参数，以获得准确的分析结果。

其次，进行网格划分和单元类型选择。

将结构划分为小单元，并选择适当的单元类型，以确保每个单元的形状和大小适合结构的几何形状。

最后，进行后处理和结果分析。

得到应力和变形结果后，可以进行结果的可视化和分析，以评估结构的性能。

总之，平面桁架结构的有限元分析是一种有效的工具，可以用于评估结构的受力状态和稳定性。

通过合适的模型建立、加载条件选择以及应力和变形分析等步骤，可以得到准确的分析结果，为结构的设计和优化提供有力支持。

本 次计 算 工 程 分 析 的扶 梯 分 段 位 置 螺 栓 的最 大 应 力 为 ２２ ５ ｐ ．而 该 扶 梯 ７ ．Ｍ ａ 分 段 处 所 用 的 螺 栓 为 ｌ． Ｏ ９级 螺 栓 ，破
断 应 力 为 １４￣ ａ ００Ｐ． 安鲧 一
≈ ３． ８２
从 材料 特 性 来 说 ， 螺 栓 是 安全 的 。 该
本 工程 桁架所 用 材料 为 Ｑ３Ｂ ２ ５，各
参 数 如 下表 １ 示 ： 所
３ 载 荷处 理 、
自动 扶 梯 桁 架 结 构 计 中 的乘 客 载
荷 以梯 级 水 平投 影 面 积 为 基 准 的载 荷 ，
国标 乘 客 载 荷 为 ５ ０Ｎ ｍ 。扶 手 带 、梯 ００／ 级 及 外 装 饰 等 重 量 则 根据 具 体 梯 种 进 行 计 算 并 转 化 成 线 载 荷 加 载 到桁 架 上 。具
图 ５ 总位 移分布 云图
度 不 得 超 过 两 支 点之 间 距 离 的 １７ ｏ ／５。 对 于 此 扶 梯 ， 允 许 挠 度 为 ６ ４ １ ００ ．／５０＝ １ ２ ４ ＜ ／ ５ ，所 以挠 度 满 足 要 求 。 ／ ３ ３ １７ ０
２ 、桁 架 连 接 螺栓 强度 分 析
Ｐ ＝ ０ ０Ｘ２ ． ９ ＝ ４ ９ ５ （） Ａ ５ ０ ８ １ ７ １ ０ ８ Ｎ
所 以乘 客载 荷 引起 的挠度 为：
１ ０ ８ ４ ９５
２４ ． ５ 、
６ ・ ４
实测挠度 为 ６ｍ ｍ ，所 以计 算 与 实 测
值相 当吻合 ，验证 了该计算的可靠性。
根 据 Ｅ １５ １ ０８规 定 ： 自动 扶 梯 ， Ｎ １— ２ ０ 以乘 客 重 量 为 基 础 计 算 或 测 得 的 最 大挠

有限元分析实验报告有限元分析实验报告一、实验基本要求根据实验指导书的要求能够独立的使用ANSYS 软件操作并在计算机上运行，学会判断结果及结构的分析，学会建立机械优化设计的数学模型，合理选用优化方法，独立的解决机械优化设计的实际问题。

二、实验目的1. 加深对机械优化设计方法的理解2. 掌握几种常用的最优化设计方法3. 能够熟练使用ANSYS 软件操作，培养学生解决案例的能力4. 培养学生灵活运用优化设计方法解决机械工程中的具体实例三、实验软件及设备计算机一台、一种应用软件如ANSYS四、实验内容实验报告例题实训1——衍架的结构静力分析图2-2所示为由9个杆件组成的衍架结构，两端分别在1，4点用铰链支承，3点受到一个方向向下的力F y , 衍架的尺寸已在图中标出，单位： m。

试计算各杆件的受力。

其他已知参数如下: 弹性模量（也称扬式模量）E=206GPa；泊松比μ=0.3；作用力F y =-1000N；杆件的2横截面积A=0.125m.一、 ANSYS8.0的启动与设置图2-2 衍架结构简图1．启动。

点击：开始>所有程序> ANSYS8.0> ANSYS ，即可进入ANSYS 图形用户主界面。

图2-4 Preference 参数设置对话框2．功能设置。

电击主菜单中的“Preference ”菜单，弹出“参数设置”对话框，选中“Structural ”复选框，点击“OK ”按钮，关闭对话框，如图2-4所示。

本步骤的目的是为了仅使用该软件的结构分析功能，以简化主菜单中各级子菜单的结构。

3．系统单位设置。

由于ANSYS 软件系统默认的单位为英制，因此，在分析之前，应将其设置成国际公制单位。

在命令输入栏中键入“/UNITS，SI ”，然后回车即可。

（注：SI 表示国际公制单位）二单元类型，几何特性及材料特性定义1．定义单元类型。

2．定义几何特性。

3．定义材料特性。

三衍架分析模型的建立1．生成节点。

图 2 梁单元臂架有限元模型
机 械振 动理 论, 各 阶 固有 振 型的线性叠加 即为结构 表现
出来的振动。其 中低阶 振型
比高阶振型对 结构的振 动贡
献 大, 基 本决 定了 机 械结 构
的动态特性。
三维结构在 无约束 边界
条 件下 的模 态 分析, 计算 出 图 3 壳单元臂架有限元模型 来的前 6阶模态接近于 0, 是 所谓的刚体模态。因此, 真正有意义的模 态应该是从 第 7阶 开始的模态。故臂架有限元模态 分析采用 Lanczos法 来求解
相吻合, 频率误差均 在 10% 以 内, 验证 了 有限 元模 型 的正 确 性和 可 靠性, 确保 了 臂架 系 统 有限 元 分析 的 准
确性。
关键词: 桁架; 有限元模态; 试验模态
中图分类号: TU 322
文献标识码: A
F in ite E lem ent and Experim entalM odal Analysis of the T russ
采用有限元分 析方法 进行 履带起 重机 臂架 的结构 设计 计算将会大大的提高设计效率 [ 1, 2]。但如何保证 有限元模型 建立的正确性及结果的可靠性 和准确性, 很多学者 都做了大 量的研究工 作 [ 3- 6] 。下 面通过 建立 梁单 元及 壳单 元两 种有 限元模型, 进行有限元模态分析, 并与试验模 态分析做 比较, 来修正和验证有限元模型的正 确性, 确保 有限元分 析结果的 可靠性。
2. 北京邮电大学 自动化学院, 北京 100876)
摘 要: 在 ABAQU S软件中分别用 梁单元和壳单元建立了桁架 有限元模 型, 并利 用有限元 模态分析 得到
其低阶固有频率和振型。同时采用锤激法对中间节臂架进行了试验模态 分析, 利用 LM S P olyM AX 分析方 法对

考试题目：下图所示为一个有3根杆组成的桁架，承受纵向和横向载荷，试对该结构进行优化设计，使得桁架重量最少。

系学数K=班号（为5,6,7,8之一）×100+学号最后两位数，如7班同学，号最后两位为20号，那么K=720已知桁架的材料特性为：弹性模量E=0.5K×103MPa；泊松比：0.5K×10-3；密度ρ= K×10 kg/m3许用应力：σ=0.5K×10-2MPa几何属性如下所示：横截面面积变化范围：0.6×10-3～0.645m2）基本尺寸B变化范围：10～0.5K×10-1m集中载荷为：Fx= 2K×103N, Fy= -2K×103N要求：写出操作步骤和命令流，定义工作文件名和工作标题为你的姓名拼音+学号。

GUI操作方式（1）定义工作文件名及工作标题1）定义工作文件名2）定义工作标题（2）定义参数和材料属性定义参数的初始值2)设置材料属性(3)定义单元类型及属性1）定义单元类型定义实常数A2 A3同A1做法(4)建立有限元模型1）生成有限元节点（节点1 2 3 4做法雷同）2）关闭坐标符号的显示3）打开节点编号显示4）生成第一个单元5）改变第二个单元的属性6）生成第二个单元7）改变第三个单元的属性8）生成第三个单元（5）施加约束和载荷1）施加边界约束2)施加集中载荷3）保存数据4）求解运算结果如下：S O L U T I O N O P T I O N SPROBLEM DIMENSIONALITY. . . . . . . . . . . . .2-D DEGREES OF FREEDOM. . . . . . UX UYANALYSIS TYPE . . . . . . . . . . . . . . . . .STATIC (STEADY-STATE) GLOBALLY ASSEMBLED MATRIX . . . . . . . . . . .SYMMETRICL O A D S T E P O P T I O N SLOAD STEP NUMBER. . . . . . . . . . . . . . . . 1TIME AT END OF THE LOAD STEP. . . . . . . . . . 1.0000NUMBER OF SUBSTEPS. . . . . . . . . . . . . . . 1STEP CHANGE BOUNDARY CONDITIONS . . . . . . . . NOPRINT OUTPUT CONTROLS . . . . . . . . . . . . .NO PRINTOUTDATABASE OUTPUT CONTROLS. . . . . . . . . . . .ALL DATA WRITTENFOR THE LAST SUBSTEP 5）保存优化结果到文件（6）进入后处理，得到状态变量和目标函数的值1）定义单元表2)计算单元体积的总和结果如下：SUM ALL THE ACTIVE ENTRIES IN THE ELEMENT TABLETABLE LABEL TOTALVOLU 65.43743)取出体积的值4)计算初始重量5)设置单元表6)得到第一杆的轴向应力7）得到第二杆的轴向应力8）得到第三杆的轴向应力9）计算轴向力的绝对值（7）显示当前设计并生成分析文件1）显示杆的当前设计2）改变视图方向3)生成优化分析文件（8）进入处理器并分析文件1）指定分析文件2）指定设计变量（A1 A2 A3 B做法雷同）3）设置状态变量3）保存优化数据库4）设置目标函数5）指定一阶优化方法6）保存数据7）运行优化8）保存优化结果到文件（9）查看优化结果1）查看最佳设计序列如下：LIST OPTIMIZATION SETS FROM SET 10 TO SET 10 AND SHOW ONLY OPTIMIZATION PARAMETERSSET 10(FEASIBLE)SIG1 (SV) 0.26212E+07SIG2 (SV) 65340.SIG3 (SV) 0.25491E+07A1 (DV) 0.57422A2 (DV) 0.60000E-03A3 (DV) 0.30810E-01B (DV) 23.873WT (OBJ) 114.462）列出所有序列的结果如下：LIST OPTIMIZATION SETS FROM SET 1 TO SET 10 AND SHOWONLY OPTIMIZATION PARAMETERS. (A "*" SYMBOL IS USED TOINDICATE THE BEST LISTED SET)SET 1 SET 2 SET 3 SET 4(FEASIBLE) (FEASIBLE) (FEASIBLE) (FEASIBLE) SIG1 (SV) 0.16434E+07 0.24929E+07 0.25179E+07 0.24641E+07 SIG2 (SV) 0.96269E+06 0.12070E+07 0.20381E+06 50744. SIG3 (SV) 0.68072E+06 0.11304E+07 0.22835E+07 0.24060E+07 A1 (DV) 0.64500 0.44325 0.58658 0.61001 A2 (DV) 0.64500 0.50234 0.26012 0.23868 A3 (DV) 0.64500 0.44325 0.65182E-01 0.49212E-01 B (DV) 26.500 23.282 22.787 22.891 WT (OBJ) 346.82 236.29 157.26 155.87SET 5 SET 6 SET 7 SET 8(FEASIBLE) (FEASIBLE) (FEASIBLE) (FEASIBLE) SIG1 (SV) 0.24181E+07 0.25032E+07 0.24888E+07 0.25939E+07 SIG2 (SV) 0.15564E+06 0.10283E+07 0.11761E+07 0.11863E+06 SIG3 (SV) 0.22407E+07 0.13731E+07 0.12008E+07 0.24629E+07 A1 (DV) 0.61469 0.59326 0.60418 0.58028 A2 (DV) 0.23115 0.28213E-01 0.60000E-03 0.60000E-03 A3 (DV) 0.58005E-01 0.68621E-01 0.60115E-01 0.32016E-01 B (DV) 23.020 23.995 24.088 23.863 WT (OBJ) 157.61 129.37 126.17 115.81SET 9 *SET 10*(FEASIBLE) (FEASIBLE)SIG1 (SV) 0.26140E+07 0.26212E+07SIG2 (SV) 0.19452E+06 65340.SIG3 (SV) 0.23994E+07 0.25491E+07A1 (DV) 0.57577 0.57422A2 (DV) 0.60000E-03 0.60000E-03A3 (DV) 0.32593E-01 0.30810E-01B (DV) 23.885 23.873WT (OBJ) 115.11 114.463）显示目标函数的变化规律a设置坐标轴标题b显示目标函数的变化规律4)显示基本尺寸B的变化规律a设置坐标轴标题b显示基本尺寸的变化规律5)显示杆面积的变化规律a 设置坐标标题b 显示杆横截面的变化规律6）显示杆中应力的变化规律a设置坐标轴标题b显示杆中应力的变化规律（10）退出ANSYS命令流方式：/BATCH/FILNAME,zhangliwen+2,1 /TITLE,zhangliwen+2*SET,B,26.5*SET,A1,0.645*SET,A2,0.645*SET,A3,0.645/PREP7 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.65e11 MPDATA,PRXY,1,,0.265 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,5.3 ET,1,LINK1R,1,A1, ,R,2,A2, ,R,3,A3, ,N,1,-B,0,0,,,,N,2,0,0,0,,,,N,3,B,0,0,,,,N,4,0,-26.5,0,,,,/PLOPTS,INFO,3/PLOPTS,LEG1,1/PLOPTS,LEG2,1/PLOPTS,LEG3,1/PLOPTS,FRAME,1/PLOPTS,TITLE,1/PLOPTS,MINM,1/PLOPTS,FILE,0/PLOPTS,LOGO,1/PLOPTS,WINS,1/PLOPTS,WP,0/PLOPTS,DATE,2/TRIAD,OFF/REPLOT/PNUM,KP,0/PNUM,LINE,0/PNUM,AREA,0/PNUM,VOLU,0/PNUM,NODE,1/PNUM,TABN,0/NUMBER,0/PNUM,ELEM,0/REPLOTFLST,2,2,1FITEM,2,1FITEM,2,4E,P51XTYPE,1MAT,1REAL,2ESYS,0SECNUM,TSHAP,LINEFLST,2,2,1FITEM,2,2FITEM,2,4E,P51XTYPE,1MAT,1REAL,3ESYS,0SECNUM,TSHAP,LINEFLST,2,2,1FITEM,2,3FITEM,2,4E,P51XFINISH/SOLFLST,2,3,1,ORDE,2FITEM,2,1FITEM,2,-3D,P51X,ALL,FLST,2,1,1,ORDE,1FITEM,2,4F,P51X,FX,1.06e6FLST,2,1,1,ORDE,1FITEM,2,4F,P51X,FY,-1.06e6SAVE/STATUS,SOLUSOLVESAVE,'zhangliwen_2','db','C:\DOCUME~1\ADMINI~1\' FINISH/POST1ETABLE,EVOLUME,VOLU,SSUM*GET,VTOT,SSUM, ,ITEM,EVOLUME*SET,DENS,5.3*SET,WT,DENS*VTOTAVPRIN,0, ,ETABLE,SIGMA,LS,1*GET,sig1,ELEM,1,ETAB,SIGMA*GET,sig2,ELEM,2,ETAB,SIGMA*GET,sig3,ELEM,3,ETAB,SIGMA*SET,sig1,abs(sig1)*SET,sig2,abs(sig2)*SET,sig3,abs(sig3)/SHRINK,0/ESHAPE,2/EFACET,1/RATIO,1,1,1/CFORMAT,32,0/REPLOT/VIEW,1,1,1,1/ANG, 1/REP,FASTEPLOTLGWRITE,'zhangliwen_2','lgw','C:\DOCUME~1\ADMINI~1\',COMMENT FINISH/OPTOPANL,'zhangliwen_2','lgw','OPVAR,A1,DV,0.6e-3,0.645, ,OPVAR,A2,DV,0.6e-3,0.645, ,OPVAR,A3,DV,0.6e-3,0.645, ,OPVAR,B,DV,10,26.5, ,OPVAR,SIG1,SV, ,2.65E6, ,OPVAR,SIG2,SV, ,2.65E6, ,OPVAR,SIG3,SV, ,2.65E6, ,OPSAVE,'zhangliwen_var','opt',' 'OPVAR,WT,OBJ, , ,1,SAVEOPTYPE,FIRSOPFRST,15, , ,SAVEOPEXE! OPTIMIZATION LOOPING HAS CLEARED THE INTERNAL LOG KEYW,BETA,0SAVEOPLIST,16, ,0OPLIST,ALL, ,0/VIEW,1,,,1/AXLAB,X,Iteration Number /AXLAB,Y,Structural Weight /GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1 XVAROPT,' 'PLVAROPT,WT/AXLAB,X,Iteration Number /AXLAB,Y,Base Dimension/GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/YRANGE,DEFAULT,,1 XVAROPT,' 'PLVAROPT,B/AXLAB,X,Iteration Number/AXLAB,Y,Cross_Sec-tional Area /GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1 XVAROPT,' 'PLVAROPT,A1,A2,A3/AXLAB,X,Iteration Number/AXLAB,Y,Maximum Stress/GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/YRANGE,DEFAULT,,1 XVAROPT,' ' PLVAROPT,SIG1,SIG2,SIG3 SAVE/DIST,1,1.,1/REP,FAST/DIST,1,1.,1/REP,FASTSAVEFINISH。

・新材料新工艺・航天器发动机推力支架桁架结构的有限元分析与优化设计3彭超义 曾竟成 肖加余 杜　刚 杨孚标(国防科技大学航天与材料工程学院,长沙　410073)文　摘　对某航天器上的推力支架结构和承载状态进行了分析,设计了一系列桁架结构,以碳纤维/环氧复合材料管件作为桁架结构的组成构件,采用有限元分析软件Ansys7.0对桁架结构的尺寸进行了优化,并优选出桁架质量和载荷相同条件下承载性能最好的桁架结构形式。

关键词　推力支架,结构分析,桁架结构,优化设计,有限元Finite Element Analysis and Optimum Design forThrust T russ Structure of S pacecraftPeng Chaoyi Z eng Jingcheng X iao Jiayu Du G ang Y ang Fubiao(Astronautics and Materials Institute,National University of Defense and T echnology,Changsha　410073)Abstract　A series of thrust truss structures are designed by analyzing mechanical properties of the thrust truss structure and its loading status on s ome spacecraft.Carbon/epoxy com posite tubes are used as the truss spars.Finite ele2 ment analysis s oftware Ansys7.0is used to optimize dimensions of the truss structures,and the best truss structure is se2 lected under same truss mass and load conditions.K ey w ords　Thrust,Structural analysis,T russ structure,Optimum design,Finite1　引言航天器发动机产生的巨大推力经由推力支架传递到液氢和液氧储罐上,从而推动航天器飞行前进。

三杆桁架的优化设计本文旨在介绍三杆桁架优化设计的背景和目的。

三杆桁架是一种常用的结构形式，具有高强度、轻量化和刚性好的特点，广泛应用于建筑、航空航天等领域。

然而，在实际应用中，三杆桁架结构的设计效果往往会受到诸多因素的制约，包括材料选择、结构形式、荷载条件等。

因此，对三杆桁架进行优化设计，不仅可以提高结构的性能和稳定性，还可以降低材料的使用量，减少成本。

本文将分析三杆桁架优化设计的背景和目的，探讨优化设计的方法和策略，以期为相关领域的研究者和设计师提供参考和指导。

三杆桁架的结构分析三杆桁架是一种常见的结构形式，由三根杆件和若干个节点组成。

它具有简单的结构和良好的稳定性，在工程领域得到广泛应用。

三杆桁架的基本结构是由三根杆件连接而成的三角形，每个顶点都是一个节点，杆件在节点处连接。

三杆桁架的性质取决于杆件的材料特性和连接方式。

三杆桁架的优化设计在设计三杆桁架时，可以采用优化设计的方法来提高其性能和效率。

优化设计的目标是使得三杆桁架在给定约束条件下，达到最佳的结构性能。

优化设计中的关键是确定合适的优化目标和设计变量。

优化目标可以包括最小化杆件的重量、最大化桁架的刚度或最小化应力集中等。

而设计变量可以包括杆件的截面积、材料的选择等。

进行三杆桁架的优化设计时，可以采用数值计算方法，如有限元分析和遗传算法等。

通过建立数学模型和进行参数优化，可以找到最优的设计方案。

总之，三杆桁架的优化设计是一项复杂而重要的工作。

通过合理的优化设计，可以提升三杆桁架的性能，实现结构的优化和效率的提高。

本文将介绍三杆桁架的优化设计方法，包括有限元分析等相关内容。

通过优化设计，我们可以改善三杆桁架的性能和结构强度，以满足特定的工程需求。

三杆桁架的优化设计可以采用以下方法：1.确定设计目标在开始优化设计之前，需要明确设计目标。

这可以包括改善结构强度、减小重量、降低成本等方面。

明确设计目标可以帮助我们选择适当的优化方法和评估指标。

2.建立数学模型根据设计目标，我们需要建立三杆桁架的数学模型。

（ ）方案 一 ａ
’ —— 一 一 ｉ ——
１ ． ＋１０风 ＋Ｏ７活 ０恒 ． ．
桁架的上下弦杆是主要的受力构件 ，其轴力随 Ｘ方向坐
标 的变 化 如下 图 ３ 图 ４所示 。 、
从 图 中可 以看 出 ，该 三角 型空 问桁 架 呈现 出与 普通 桁架 相 似 的受 力特 征 。三角 形桁 架 结构 在 对称 荷 载作 用 下 ，杆件 的轴 力分 布 也是 对 称 的 。屋盖 上 弦杆 受压 ，下 弦杆 受拉 。且
第 １ 卷 第 ６ １ 期
２ １ 每 ０ １
中
国
水
运
Ｖｏ１ １ ．１
Ｊ ｅ ｕｎ
Ｎｏ ６ ．
２ Ｏ１１
６月
Ｏ ｎ Ｗ ｔｅ Ｔ ａｎ ｐ ｔ ｈ ｉ ａ ａ ｒ ｒ ｓ ｏｒ
某仓库钢桁架结构静动力有限元分析
杜双庆
（中交 第 二航 务工 程勘察 设 计 院有 限公 司 ，湖 北 武 汉 ４ ０７ ３ ０ １）
一
利 用 ｓ ｐ ０ ０建立 仓库 整体 模 型 ， ａ２ ０ 屋盖 各杆 件之 间均 为
铰接 。 钢柱和地面为固接。结构整体有限元模型如图 ２所对 于 本工 程 的三 角形 截 面桁 架 结构 在横 向 （ 度 方 向） 跨
是 稳定 的 ，可 以等效 为一 根 梁 ，而 在纵 向是不 稳 定 的 ，要靠 上 弦平 面 内水 平支 撑 才 能保 证 结构 的纵 向稳 定 性 ，而 檩 条不 足 以做 为 结构 纵 向支 撑 ，只 是做 为 屋 面板 的支 承 构件 而 已。 因此 ，在 屋盖 的 两端 和 中部 设置 水 平支 撑 ，其 它榀 桁 架 用刚 性 系杆 与之 连接 。 在 不 影响 使 用 的条 件 下 ，对 主桁 架 的 两种 结构 方案进 行 了比较 ，如图 １所 示 ，方案 二与 方 案一 的不 同为 加 设一 斜撑 杆 。结 果 表 明 ：方案 一主 桁 架 的跨 中挠 度 和 上、 下弦 杆 内力 均 较大 ，用钢 量 高 ，主 桁架 对 柱 的水 平推 力较 小 ，柱 顶 承受

三杆桁架的优化设计班级：机自Y107姓名：罗展雄学号：201000104091考试题目：下图所示为一个有3根杆组成的桁架，承受纵向和横向载荷，试对该结构进行优化设计，使得桁架重量最少。

系数K=班号（为5,6,7,8之一）×100+学号最后两位数，如7班同学，学号最后两位为20号，那么K=720K=791已知桁架的材料特性为：弹性模量E=0.5K ×103MPa=395500MPa泊松比：0.5K ×10-3=0.3955 密度ρ= K ×10 kg/m3=7910kg/m 3许用应力：σ=0.5K ×10-2MPa=3.955MPa 几何属性如下所示：横截面面积变化范围：0.6×10-3～0.645m 2） 基本尺寸B 变化范围：10～0.5K ×10-1m集中载荷为：Fx= 2K ×103N=11582×103N, Fy= -2K×103N=-11582×103N1、定义工作文件名及工作工作标题1）定义工作文件名：2）定义工作标题2、定义参数和材料属性1）定义参数的初始值：2）设计材料属性：3、定义单元类型及属性1）定义单元类型：2）定义实常数：3）打开节点编号显示：4）生成第一单元：6）生成第二个单元：8）生成第三个单元：1）施加边界载荷：2）施加集中载荷：3）求解运算：4）保存优化结果到文件：6、进入后处理，得到状态变量和目标函数的值1）定义单元表：2）计算单元体积的总和：3）取出体积的值：4）计算初始重量：）设计单元表：6）得到第一杆的轴向应力：7）得到第二杆的轴向应力：8）得到第三杆的轴向应力：9）计算轴向力的绝对应力：7、显示当前设计并生成分析文件1）显示杆的当前设计：2）改变视图方向：3）生成优化分析文件：8、进入优化处理器并指定分析文件1）指定分析文件：2）指定设计变量：3）设置状态变量：5）设计目标函数：6）指定一阶优化方法：7）运行优化：8）保存优化结果到文件：9、查看优化结果1）查看最佳设计序列：2）列出所有序列的结果：3）显示目标函数的变化规律：4）显示基本尺寸B的变化规律：5）显示杆面积的变化规律：6）显示杆中应力的变化规律：命令流方式：/filename,luozhanxiong/TITLE,201000104091*SET,b,25*SET,a1,0.645*SET,a2,0.645*SET,a3,0.645/PREP7 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,3.955e11 MPDATA,PRXY,1,,0.3955 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7910 ET,1,LINK1R,1,a1, ,R,2,a2, ,R,3,a3, ,N,1,-b,0,0,,,,N,2,0,0,0,,,,N,3,b,0,0,,,,N,4,0,-25,0,,,,/PLOPTS,INFO,3/PLOPTS,LEG1,1/PLOPTS,LEG2,1/PLOPTS,LEG3,1/PLOPTS,FRAME,1/PLOPTS,TITLE,1/PLOPTS,MINM,1/PLOPTS,FILE,0/PLOPTS,LOGO,1/PLOPTS,WINS,1/PLOPTS,WP,0/PLOPTS,DATE,2/TRIAD,OFF/REPLOT/PNUM,KP,0/PNUM,LINE,0/PNUM,AREA,0/PNUM,VOLU,0/PNUM,NODE,1/PNUM,TABN,0/PNUM,SV AL,0/NUMBER,0/PNUM,ELEM,0/REPLOTFLST,2,2,1FITEM,2,1FITEM,2,4E,P51XTYPE, 1MAT, 1REAL, 2ESYS, 0SECNUM,TSHAP,LINEFLST,2,2,1FITEM,2,2FITEM,2,4E,P51XTYPE, 1MAT, 1REAL, 3ESYS, 0SECNUM,TSHAP,LINEFLST,2,2,1FITEM,2,3FITEM,2,4E,P51XFINISH/SOLFLST,2,3,1,ORDE,2FITEM,2,1FITEM,2,-3D,P51X, , , , , ,ALL, , , , ,FLST,2,1,1,ORDE,1FITEM,2,4FLST,2,1,1,ORDE,1FITEM,2,4F,P51X,FX,1582000FLST,2,1,1,ORDE,1FITEM,2,4F,P51X,FY,-1582000SAVE/STATUS,SOLUSOLVESAVE,'truss_resu','db','C:\Users\luo\'FINISH/POST1A VPRIN,0, ,ETABLE,evolumb,VOLU,SSUM*GET,vtot,SSUM, ,ITEM,EVOLUMB*SET,dens,7910*SET,wt,dens*vtotA VPRIN,0, ,ETABLE,sigma,LS, 1*GET,sig1,ELEM,1,ETAB,SIGMA*GET,sig2,ELEM,2,ETAB,SIGMA*GET,sig3,ELEM,3,ETAB,SIGMA*SET,sig1,abs(sig1)*SET,sig2,abs(sig2)*SET,sig3,abs(sig3)/SHRINK,0/ESHAPE,2/EFACET,1/RATIO,1,1,1/CFORMAT,32,0/REPLOT/VIEW, 1 ,1,1,1/ANG, 1/REP,FASTEPLOT/REPLOT,RESIZE! LGWRITE,'truss_opt','lgw','C:\Users\luo\',COMMENT FINISH/OPTOPANL,'truss_opt','lgw',' 'OPV AR,A1,DV,0.6e-3,0.645, ,OPV AR,A2,DV,0.6e-3,0.645, ,OPV AR,A3,DV,0.6e-3,0.645, ,OPV AR,B,DV,10,39.55, ,OPV AR,DENS,SV, ,3.955e6, ,OPV AR,SIG2,SV, ,3.955e6, ,OPV AR,,DEL,2OPV AR,,DEL,1OPV AR,SIG1,SV, ,3.955e6, ,OPV AR,SIG2,SV, ,3.955e6, ,OPV AR,SIG3,SV, ,3.955e6, ,OPSAVE,'truss_var','opt',' 'OPV AR,WT,OBJ, , ,1,OPTYPE,FIRSOPFRST,15, , ,SAVEOPEXEKEYW,BETA,0SAVE,'truss_opt_resu','db','C:\Users\luo\' OPLIST,11, ,0OPLIST,ALL, ,0/AXLAB,X,iteration number/AXLAB,Y,structural weight/GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1XV AROPT,' 'PLV AROPT,WT/AXLAB,X,iteration number/AXLAB,Y,base dimension/GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1XV AROPT,' 'PLV AROPT,B/AXLAB,X,iteration number/AXLAB,Y,cross_sec_tional area /GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,DIG2,3,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1XV AROPT,' 'PLV AROPT,A1,A2,A3/AXLAB,X,iteration number/AXLAB,Y,maximun stress/GTHK,AXIS,2/GRTYP,0/GROPT,ASCAL,ON/GROPT,LOGX,OFF/GROPT,LOGY,OFF/GROPT,AXDV,1/GROPT,AXNM,ON/GROPT,AXNSC,1,/GROPT,DIG1,4,/GROPT,XAXO,0,/GROPT,YAXO,0,/GROPT,DIVX,/GROPT,DIVY,/GROPT,REVX,0/GROPT,REVY,0/GROPT,LTYP,0/XRANGE,DEFAULT/YRANGE,DEFAULT,,1XV AROPT,' 'PLV AROPT,SIG1,SIG2,SIG3 finish。

《有限元技术基础》考试题目：下图所示为一个有3根杆组成的桁架，承受纵向和横向载荷，杆件的横截面面积和基本尺寸B在指定范围内变化，要求桁架的每根杆件承受的最大应力小于（800+学号最后两位数）MPa，试对该结构进行优化设计，使得桁架重量最少。

弹性模量E=220GPa；泊松比：0.3；密度ρ=7800kg/m3材料最大许用应力：σ=890MPa横截面面积变化范围：0.01～10cm2（初始值为10）基本尺寸B变化范围：1～2m（初始值为2）要求：写出操作步骤和命令流，定义工作文件名和工作标题为你的姓名拼音。

GUI操作方式：(1) 定义工作文件名和工作标题：1)定义工作文件名：Utility Menu- File-Change Jobname,输入文件名“litao”,单击“OK”。

2)定义工作标题：Utility Menu- File-Change Title,输入工作标题“litao”,单击“ＯＫ”。

(2) 定义参数和材料属性：１）定义参数初始值：Utility Menu-Parameters-Scaler Parameters, 分别在“Selection”’下面的输入栏中输入：B=2，A1=0.001，A2=0.001，A3=0.001。

”所得结果如图所示，单击”close”。

２）设置材料属性：Main Menu-Preprocessor-Material Props-Material Models,设置材料属性“EX=2.2e11,PRXY=0.3”设置单元密度为“DENS=7800”，单击“OK”,设置如图所示，完成对材料属性的设置。

(3)定义单元类型及属性1)定义单元类型：Main Menu-Preprocessor-Element Type-Add/Edit/Delete,在“Library of Element Type”左面的列表栏选择Structural Link”,右边的为“2D spar 1”如图所示，完成单元类型设置。