十杆桁架结构优化设计.docVIP

题目：十杆桁架结构优化设计

日期：2023.09.16

目录

TOC\o1-3\h\z\u1 设计题目 1

2 设计过程 2

2.1 一、运用Abaqus求解各杆轴力应力 2

2.1.1 Abaqus计算流程 2

2.1.2 结果 3

2.2 二、运用材料力学知识求解 4

2.2.1 基本思绪 4

2.2.2 解题过程 4

2.2.3 结果 5

2.3 三、编写有限元程序求解 6

2.3.1 程序基本环节 6

2.3.2 Vs2023中重要的程序段 6

2.3.3 程序输出文献 9

2.3.4 材料力学、有限元程序、Abaqus结果比较 10

2.4 四、装配应力计算 11

2.4.1 解决技巧 11

2.4.2 Abaqus解决技巧 11

2.4.3 不加外力（P1,P2,P3）时材力，Ansys与Abaqus结果 12

2.4.4 不加外力（P1,P2,P3）时材力，Ansys与Abaqus误差分析 12

2.4.5 加外力（P1,P2,P3）时Ansys与Abaqus结果 12

2.5 五、优化设计 14

2.5.1 设计中变量的概念 14

2.5.2 优化环节运用VS2023编写复合形法进行约束优化。 14

2.5.3 VS2023优化程序 16

2.5.4 优化结果 19

2.5.5 结果说明 19

3 设计感想 20

4 备注 20

4.1 参考书目 20

4.2 说明 20

设计题目

十字桁架结构优化设计

现有十字桁架结构见图1，材料泊松比为0.3，E=2.1e11，密度为7.8×103kg/m3,许用应力为160Mpa，P1=600kN，P2=900kN，P3=600kN，杆1-6面积为A1=0.03m2,杆7-10面积为A2=0.02m。

1、运用Ａｂａｑｕｓ计算各杆的应力；

2、运用材料力学的知识求解，并与1计算出的结果做比较；

3、编写有限元程序求解，与1和2计算结果进行比较；

4、若杆5制作时短了0.001m，试求各杆的应力；

5、若令2节点的位移小于0.005m，A1、A2为0.005～0.05m2，试对结构进行优化，使其重量最小。（同材料力学优化结果比较）。

图

图SEQ图\*ARABIC1十杆桁架

设计过程

一、运用Abaqus求解各杆轴力应力

运用Abaqus求解，十字桁架结构可用2Dtruss单元模拟。单元参数为：弹性模量，1-6杆截面面积，7-10杆截面面积。加载求解输出各杆应力，结点位移。

Abaqus计算流程

Part:创建parttrussmain,part45l,part45r,parttrussmain,包涵除8,10杆外的所有杆，part45l包涵8杆，part45r包涵10杆。

Property：createMaterial:Elastic:弹性模量,泊松比0.3

CreateSection：beam\Truss：SectionA1,截面面积30000。SectionA2,截面面积20230。并给各杆赋材料属性。

Assembly：组装parttrussmain,part45l,part45r。

Step:创建一个分析步，step1。

Interaction：用Tie把parttrussmain,part45l,part45r,绑定。

Load：createload：5,6点加铰接约束，固定x,y方向位移。

CreateboundaryCondition：2,4点加相应力。

Mesh：划分网格，一个杆为一个单元。Elementtape，选truss

Job：创建一个job，WriteInput,DataCheck,Submit,通过Result来查看应力云图。

结果

图

图SEQ图\*ARABIC2Abaqus各杆应力云图

二、运用材料力学知识求解

基本思绪

显然题目中的十字桁架结构是两次静不定问题。对于一次静不定问题，材料力学给出了两类解法：①去掉约束加力，找位移协调关系解题；②力法正则方程求解。对于多次静不定，特别是上述桁架问题，找出其协调关系基本上是不也许的，而力法正则方程更适合于解这种结构。如图3所示，去掉多余约束，建立力法正则方程：

图

图SEQ图\*ARABIC3去多余约束

解题过程

分别求出外力作用下各杆内力和单位力作用下的各杆内力，为计算方便，将其结果列入下表1中。应用莫尔积分