计算结构力学 安建工 08土木完整版.docx
《计算结构力学 安建工 08土木完整版.docx》由会员分享，可在线阅读，更多相关《计算结构力学 安建工 08土木完整版.docx（30页珍藏版）》请在冰豆网上搜索。
计算结构力学安建工08土木完整版
《计算结构力学》
课程设计
姓名:
0000
学号:
班级:
指导:
苏少卿晏燕
安徽建筑工业学院土木工程学院
矩形悬臂柱长L=4m，高h=2m，横截面为矩形，宽t=15mm。
悬臂端受均匀分布压力6KN/m，材料弹性模量E=210GPa，泊松比ν=0.28。
讨论：
1、网格大小对有限元计算结果精度的影响（在模板的基础上进行分析）。
2、距固定端
的截面上各点应力状态，检验它是否为单向应力状态。
3、距固定端
的截面上各点变形，检验平面假定。
4、假设材料泊松比可以变化，讨论它的影响。
5、增大或者减小高度，对结果精度有什么影响？
（主要考虑减小高度方面分析）
【解】建立结构坐标系，把板放在第一象限
1、网格大小对有限元计算结果精度的影响（在模板基础上进行分析）。
方案1：
按16X8剖分，如图，板中点
是结点77（2,1），考察它的位移：
输入数据如下：
5
悬臂柱长4m,高2m,自由端载荷6kN/m均布压力
杨氏模量210GPa,泊桑比0.28,板厚15mm
自动划分为256个单元,结点编号顺序:
从左到右,自下向上
10181
171
80002
80.2500.252
80.500.52
80.7500.752
81012
81.2501.252
81.501.52
81.7501.752
82022
82.2502.252
82.502.52
82.7502.752
83032
83.2503.252
83.503.52
83.7503.752
84042
000
1.102.103.104.105.106.107.108.109.10
1.202.203.204.205.206.207.208.209.20
2100.280.015
0
8
1451461-6-6
1461471-6-6
1471481-6-6
1481491-6-6
1491501-6-6
1501511-6-6
1511521-6-6
1521531-6-6
计算结果：
Ua=-3.754μm,Va=0.036μm
方案2：
按32X8剖分，如图，板中点
是结点149（2,1），考察它的位移：
输入数据如下：
5
悬臂柱长4m,高2m,自由端载荷6kN/m均布压力
杨氏模量210GPa,泊桑比0.28,板厚15mm
自动划分为512个单元,结点编号顺序:
从左到右,自下向上
10181
331
80002
80.12500.1252
80.2500.252
80.37500.3752
80.500.52
80.62500.6252
80.7500.752
80.87500.8752
81012
81.12501.1252
81.2501.252
81.37501.3752
81.501.52
81.62501.6252
81.7501.752
81.87501.8752
82022
82.12502.1252
82.2502.252
82.37502.3752
82.502.52
82.62502.6252
82.7502.752
82.87502.8752
83032
83.12503.1252
83.2503.252
83.37503.3752
83.503.52
83.62503.6252
83.7503.752
83.87503.8752
84042
000
1.102.103.104.105.106.107.108.109.10
1.202.203.204.205.206.207.208.209.20
2100.280.015
0
8
2892901-6-6
2902911-6-6
2912921-6-6
2922931-6-6
2932941-6-6
2942951-6-6
2952961-6-6
2962971-6-6
计算结果：
Ua=-3.757μm,Va=0.020μm
方案3：
按32X16剖分，如图，板中点
是结点281（2,1）：
输入数据如下：
5
悬臂柱长4m,高2m,自由端载荷6kN/m均布压力
杨氏模量210GPa,泊桑比0.28,板厚15mm
自动划分为1024个单元,结点编号顺序:
从左到右,自下向上
10341
331
160002
160.12500.1252
160.2500.252
160.37500.3752
160.500.52
160.62500.6252
160.7500.752
160.87500.8752
161012
161.12501.1252
161.2501.252
161.37501.3752
161.501.52
161.62501.6252
161.7501.752
161.87501.8752
162022
162.12502.1252
162.2502.252
162.37502.3752
162.502.52
162.62502.6252
162.7502.752
162.87502.8752
163032
163.12503.1252
163.2503.252
163.37503.3752
163.503.52
163.62503.6252
163.7503.752
163.87503.8752
164042
000
1.102.103.104.105.106.107.108.109.1010.1011.1012.1013.1014.1015.1016.1017.10
1.202.203.204.205.206.207.208.209.2010.2011.2012.2013.2014.2015.2016.2017.20
2100.280.015
0
16
5455461-6-6
5465471-6-6
5475481-6-6
5485491-6-6
5495501-6-6
5505511-6-6
5515521-6-6
5525531-6-6
5535541-6-6
5545551-6-6
5555561-6-6
5565571-6-6
5575581-6-6
5585591-6-6
5595601-6-6
5605611-6-6
计算结果：
Ua=-3.760μm,Va=0.012μm
方案4：
按32X32剖分，如图，板中点
是结点545（2,1）：
输入数据如下：
5
悬臂柱长4m,高2m,自由端载荷6kN/m均布压力
杨氏模量210GPa,泊桑比0.28,板厚15mm
自动划分为2048个单元,结点编号顺序:
从左到右,自下向上
10661
331
320002
320.12500.1252
320.2500.252
320.37500.3752
320.500.52
320.62500.6252
320.7500.752
320.87500.8752
321012
321.12501.1252
321.2501.252
321.37501.3752
321.501.52
321.62501.6252
321.7501.752
321.87501.8752
322022
322.12502.1252
322.2502.252
322.37502.3752
322.502.52
322.62502.6252
322.7502.752
322.87502.8752
323032
323.12503.1252
323.2503.252
323.37503.3752
323.503.52
323.62503.6252
323.7503.752
323.87503.8752
324042
000
1.102.103.104.105.106.107.108.109.1010.1011.1012.1013.1014.1015.1016.1017.1018.1019.1020.1021.1022.1023.1024.1025.1026.1027.1028.1029.1030.1031.1032.1033.10
1.202.203.204.205.206.207.208.209.2010.2011.2012.2013.2014.2015.2016.2017.2018.2019.2020.2021.2022.2023.2024.2025.2026.2027.2028.2029.2030.2031.2032.2033.20
2100.280.015
0
32
105710581-6-6
105810591-6-6
105910601-6-6
106010611-6-6
106110621-6-6
106210631-6-6
106310641-6-6
106410651-6-6
106510661-6-6
106610671-6-6
106710681-6-6
106810691-6-6
106910701-6-6
107010711-6-6
107110721-6-6
107210731-6-6
107310741-6-6
107410751-6-6
107510761-6-6
107610771-6-6
107710781-6-6
107810791-6-6
107910801-6-6
108010811-6-6
108110821-6-6
108210831-6-6
108310841-6-6
108410851-6-6
108510861-6-6
108610871-6-6
108710881-6-6
108810891-6-6
计算结果：
Ua=-3.761μm,Va=0.007μm
方案5：
按64X32剖分，如图，板中点
是结点1073（2,1）：
输入数据如下：
5
悬臂柱长4m,高2m,自由端载荷6kN/m均布压力
杨氏模量210GPa,泊桑比0.28,板厚15mm
自动划分为4096个单元,结点编号顺序:
从左到右,自下向上
10661
651
320002
320.062500.06252
320.12500.1252
320.187500.18752
320.2500.252
320.312500.31252
320.37500.3752
320.437500.43752
320.500.52
320.562600.56252
320.62500.6252
320.687500.68752
320.7500.752
320.812500.81252
320.87500.8752
320.937500.93752
321012
321.062501.06252
321.12501.1252
321.187501.18752
321.2501.252
321.312501.31252
321.37501.3752
321.437501.43752
321.501.52
321.562601.56252
321.62501.6252
321.687501.68752
321.7501.752
321.812501.81252
321.87501.8752
321.937501.93752
322022
322.062502.06252
322.12502.1252
322.187502.18752
322.2502.252
322.312502.31252
322.37502.3752
322.437502.43752
322.502.52
322.562602.56252
322.62502.6252
322.687502.68752
322.7502.752
322.812502.81252
322.87502.8752
322.937502.93752
323032
323.062503.06252
323.12503.1252
323.187503.18752
323.2503.252
323.312503.31252
323.37503.3752
323.437503.43752
323.503.52
323.562603.56252
323.62503.6252
323.687503.68752
323.7503.752
323.812503.81252
323.87503.8752
323.937503.93752
324042
000
1.102.103.104.105.106.107.108.109.1010.1011.1012.1013.1014.1015.1016.1017.1018.1019.1020.1021.1022.1023.1024.1025.1026.1027.1028.1029.1030.1031.1032.1033.10
1.202.203.204.205.206.207.208.209.2010.2011.2012.2013.2014.2015.2016.2017.2018.2019.2020.2021.2022.2023.2024.2025.2026.2027.2028.2029.2030.2031.2032.2033.20
2100.280.015
0
32
211321141-6-6
211421151-6-6
211521161-6-6
211621171-6-6
211721181-6-6
211821191-6-6
211921201-6-6
212021211-6-6
212121221-6-6
212221231-6-6
212321241-6-6
212421251-6-6
212521261-6-6
212621271-6-6
212721281-6-6
212821291-6-6
212921301-6-6
213021311-6-6
213121321-6-6
213221331-6-6
213321341-6-6
213421351-6-6
213521361-6-6
213621371-6-6
213721381-6-6
213821391-6-6
213921401-6-6
214021411-6-6
214121421-6-6
214221431-6-6
214321441-6-6
214421451-6-6
计算结果：
Ua=-3.763μm,Va=0.004μm
从水平位移看，已经划分“够细”了，但从竖向位移看，还要再细分。
其分析结果如下所示：
A点位移随网格细分的变化：
a/h
u
u误差
v
v误差
0.125
-3.754
0.036
0.125
-3.757
0.08%
0.02
44.44%
0.0625
-3.76
0.08%
0.012
40.00%
0.03125
-3.761
0.03%
0.007
41.67%
0.03125
-3.763
0.05%
0.004
42.86%
A点位移散点图：
结论：
本题网格长度小于矩形板边长的1/8时，主要位移的计算已经基本收敛，但纵向位移计算的收敛速度较慢，需要划分到网格尺度小于边长的1/XXX
2、距固定端
的截面上各点应力状态，检验它是否为单向应力状态。
一、由方案一可知：
距固定端
的截面上各点应力状态为
结点应力（千帕）
结点
Sigmax
Sigmay
Tauxy
91
-399.191
0.179
0.418
92
-399.554
0.36
0.563
93
-400.106
0.674
0.655
94
-400.572
0.906
0.396
95
-400.717
0.954
-0.046
96
-400.48
0.793
-0.465
97
-399.962
0.482
-0.665
98
-399.406
0.162
-0.509
99
-399.232
0.08
-0.294
结论：
距固定端
的截面上各点应力状态为复杂应力状态，不是单向应力状态。
二、由方案二可知：
距固定端
的截面上各点应力状态为
结点应力（千帕）
结点
Sigmax
Sigmay
Tauxy
172
-399.039
0.169
0.446
173
-399.418
0.353
0.638
174
-400.085
0.703
0.769
175
-400.658
0.986
0.484
176
-400.858
1.071
-0.027
177
-400.602
0.921
-0.525
178
-399.998
0.592
-0.775
179
-399.33
0.238
-0.606
180
-399.066
0.111
-0.373
结论：
距固定端
的截面上各点应力状态为复杂应力状态，不是单向应力状态。
同理可知：
方案三、四和方案五中，距固定端
的截面上各点应力状态都为复杂应力状态。
3、距固定端
的截面上各点变形，检验平面假定。
一、由方案一可知：
距固定端
的截面上各点变形为
结构位移（微米）:
结点
U
V
91
-4.694
-0.49
92
-4.698
-0.357
93
-4.702
-0.223
94
-4.706
-0.089
95
-4.71
0.046
96
-4.715
0.181
97
-4.721
0.315
98
-4.727
0.448
99
-4.732
0.582
结论：
距固定端
的截面上各点变形不一致，但在一定的允许误差范围内符合平面假定。
二、由方案二可知：
距固定端
的截面上各点变形为
结构位移（微米）:
结点
U
V
172
-4.468
-0.513
173
-4.47
-0.38
174
-4.472
-0.246
175
-4.473
-0.111
176
-4.475
0.024
177
-4.478
0.159
178
-4.482
0.293
179
-4.486
0.427
180
-4.488
0.56
结论：
距固定端
的截面上各点变形不一致，但在一定的允许误差范围内符合平面假定。
同理可知：
方案三、四和方案五中，距固定端
的截面上各点变形不一致，但在一定的允许误差范围内符合平面假定。
4、假设材料泊松比可以变化，讨论它的影响。
当泊松比为0.1、0.25、0.28、0.32、0.45时，可从结构位移、支座反力、单元应力、结点应力中讨论它改变后的影响。
其结果如下：
①结构位移（微米）:
泊松比ν
结点
U
V
0.1
77
-3.801
0.009
0.25
77
-3.764
0.031
0.28
77
-3.754
0.036
0.32
77
-3.738
0.045
0.45
77
-3.673
0.079
②支座
升级会员 