第一步：建立框架结构模型
第二步：框架结构施加重力荷载
第三步：框架结构施加活荷载
第四步：框架结构施加风荷载
第五步：框架结构荷载组合运算
第六步：框架结构模态分析
第七步：框架结构屈曲特征值分析
第八步：框架结构地震时程弹性分析
第九步：结束
采用地震波时间文件为
:TIME.TXT
!
采用地震波文件分别为
:AC_X.TXT/AC_Y.TXT
FINISH
/CLEAR
/FILENAME,FRAME_SHEAR_WALL
/TITLE, FRAME_SHEAR_WALL
第一步：框架结构建模型
/PREP7      进入前处理器（默认单位为m，kg,s，所以建模时要特别注意） 
ET,1,BEAM4      定义BEAM4梁单元，编号为1
ET,2,SHELL63    定义Shell63壳单元，编号为2
R,1,0.25,0.0052,0.0052,0.5,0.5  定义实常数组号为1的各个量
R,2,0.2025,0.0034,0.0034,0.45,0.45
R,3,0.06,0.0002,0.00045,0.3,0.2
R,4,0.25,0.25,0.25,0.25
R,5,0.1,0.1,0.1,0.1
MP,EX,1,3.0E10    定义编号为1的单元的弹性模量为3.0E10
MP,NUXY,1,0.2    定义编号为1的单元的次泊松比为0.2
MP,DENS,1,2500    定义编号为1的密度为2500
*DO,II,1,1      利用循环命令，循环变量II从1~1
K,1+(II-1)*20,0,0,(II-1)*6      即：K,1
K,2+(II-1)*20,6,0,(II-1)*6        K,2,6         
K,3+(II-1)*20,12,0,(II-1)*6        K,3,12
K,4+(II-1)*20,17,0,(II-1)*6        K,4,17
K,5+(II-1)*20,0,5,(II-1)*6        K,5,0,5
K,6+(II-1)*20,6,5,(II-1)*6        K,6,6,5
K,7+(II-1)*20,12,5,(II-1)*6        K,7,12,5
K,8+(II-1)*20,17,5,(II-1)*6        K,8,17,5
K,9+(II-1)*20,0,10,(II-1)*6        K,9,0,10
K,10+(II-1)*20,6,10,(II-1)*6      K,10,6,10
K,11+(II-1)*20,12,10,(II-1)*6      K,11,12,10
K,12+(II-1)*20,17,10,(II-1)*6      K,12,17,10
K,14+(II-1)*20,6,1.875,(II-1)*6    K,14,6,1.875
K,16+(II-1)*20,6,3.125,(II-1)*6    K,16,6,3.125
K,18+(II-1)*20,6,1.875,3.9        K,18,6,1.875,3.9           
K,20+(II-1)*20,6,3.125,3.9        K,20,6,3.125,3.9
*ENDDO    结束循环
*DO,II,1,10        利用循环命令，循环变量II从1~10                                     
K,1+II*20,0,0,(II-1)*3+6  创建关键点
K,2+II*20,6,0,(II-1)*3+6
K,3+II*20,12,0,(II-1)*3+6
K,4+II*20,17,0,(II-1)*3+6
K,5+II*20,0,5,(II-1)*3+6
K,6+II*20,6,5,(II-1)*3+6
K,7+II*20,12,5,(II-1)*3+6
K,8+II*20,17,5,(II-1)*3+6
K,9+II*20,0,10,(II-1)*3+6
K,10+II*20,6,10,(II-1)*3+6
K,11+II*20,12,10,(II-1)*3+6
K,12+II*20,17,10,(II-1)*3+6
K,13+II*20,0,1.875,(II-1)*3+6
K,14+II*20,6,1.875,(II-1)*3+6
K,15+II*20,0,3.125,(II-1)*3+6
K,16+II*20,6,3.125,(II-1)*3+6
K,17+II*20,0,1.875,(II-1)*3+6+1.5     
K,18+II*20,6,1.875,(II-1)*3+6+2     
K,19+II*20,0,3.125,(II-1)*3+6+1.5     
K,20+II*20,6,3.125,(II-1)*3+6+2     
*ENDDO                                       
*DO,II,11,11                                     
K,1+II*20,0,0,(II-1)*3+6     
K,2+II*20,6,0,(II-1)*3+6     
K,3+II*20,12,0,(II-1)*3+6     
K,4+II*20,17,0,(II-1)*3+6     
K,5+II*20,0,5,(II-1)*3+6     
K,6+II*20,6,5,(II-1)*3+6     
K,7+II*20,12,5,(II-1)*3+6     
K,8+II*20,17,5,(II-1)*3+6     
K,9+II*20,0,10,(II-1)*3+6     
K,10+II*20,6,10,(II-1)*3+6     
K,11+II*20,12,10,(II-1)*3+6     
K,12+II*20,17,10,(II-1)*3+6     
K,13+II*20,0,1.875,(II-1)*3+6     
K,14+II*20,6,1.875,(II-1)*3+6     
K,15+II*20,0,3.125,(II-1)*3+6     
K,16+II*20,6,3.125,(II-1)*3+6     
K,17+II*20,0,1.875,(II-1)*3+6+1.8   
K,18+II*20,6,1.875,(II-1)*3+6+2.5 
K,19+II*20,0,3.125,(II-1)*3+6+1.8   
K,20+II*20,6,3.125,(II-1)*3+6+2.5 
*ENDDO                                                 
*DO,II,12,12                               
K,1+II*20,0,0,39.5
K,2+II*20,6,0,39.5
K,3+II*20,12,0,39.5
K,4+II*20,17,0,39.5
K,5+II*20,0,5,39.5
K,6+II*20,6,5,39.5
K,7+II*20,12,5,39.5
K,8+II*20,17,5,39.5
K,9+II*20,0,10,39.5
K,10+II*20,6,10,39.5
K,11+II*20,12,10,39.5
K,12+II*20,17,10,39.5
*ENDDO
*DO,II,2,13                         
L,(II-1)*20+1,(II-1)*20+2    通过两关键点创建线
L,(II-1)*20+2,(II-1)*20+3 
L,(II-1)*20+3,(II-1)*20+4 
L,(II-1)*20+5,(II-1)*20+6 
L,(II-1)*20+6,(II-1)*20+7 
L,(II-1)*20+7,(II-1)*20+8 
L,(II-1)*20+9,(II-1)*20+10 
L,(II-1)*20+10,(II-1)*20+11 
L,(II-1)*20+11,(II-1)*20+12 
L,(II-1)*20+5,(II-1)*20+9 
L,(II-1)*20+6,(II-1)*20+10 
L,(II-1)*20+3,(II-1)*20+7 
L,(II-1)*20+7,(II-1)*20+11 
L,(II-1)*20+4,(II-1)*20+8 
L,(II-1)*20+8,(II-1)*20+12 
*ENDDO                             
*DO,II,1,12                     
L,(II-1)*20+1,II*20+1 
L,(II-1)*20+2,II*20+2 
L,(II-1)*20+3,II*20+3 
L,(II-1)*20+4,II*20+4 
L,(II-1)*20+5,II*20+5 
L,(II-1)*20+6,II*20+6 
L,(II-1)*20+7,II*20+7 
L,(II-1)*20+8,II*20+8 
L,(II-1)*20+9,II*20+9 
L,(II-1)*20+10,II*20+10
L,(II-1)*20+11,II*20+11
L,(II-1)*20+12,II*20+12
*ENDDO
*DO,II,1,1                     
L,(II-1)*20+14,(II-1)*20+18
L,(II-1)*20+16,(II-1)*20+20
*ENDDO                         
*DO,II,2,12
L,(II-1)*20+13,(II-1)*20+17
L,(II-1)*20+15,(II-1)*20+19
L,(II-1)*20+14,(II-1)*20+18
L,(II-1)*20+16,(II-1)*20+20
*ENDDO
A,1,21,25,5    通过四个关键点创建面
*DO,II,2,12
A,(II-1)*20+1,II*20+1,II*20+5,(II-1)*20+5,(II-1)*20+15,(II-1)*20+19,(II-1)*20+17,(II-1)*20+13    通过8个关键点创建面
*ENDDO
*DO,II,1,12                                                                                                 
array工艺详解
A,(II-1)*20+2,II*20+2,II*20+6,(II-1)*20+6,(II-1)*20+16,(II-1)*20+20,(II-1)*20+18,(II-1)*20+14
*ENDDO
*DO,II,1,12
A,II*20+1,II*20+2,II*20+6,II*20+5
A,II*20+2,II*2
0+3,II*20+7,II*20+6
A,II*20+3,II*20+4,II*20+8,II*20+7
A,II*20+5,II*20+6,II*20+10,II*20+9
A,II*20+6,II*20+7,II*20+11,II*20+10
A,II*20+7,II*20+8,II*20+12,II*20+11
*ENDDO
立柱网格划分
底层立柱网格划分
LSEL,S,,,181,192    从所有线中选择最小线号为181，最大线号为192的一组线   
LATT,1,1,1          设置线的单元属性，材料号为1，实常数号为1，单元类型号为1     
LESIZE,ALL,,,10      对所选的每条线的分段数为10
LMESH,ALL            在所有的线上生成线单元
2~12层立柱网格划分
LSEL,S,,,192,324    从所有线中选择最小线号为192，最大线号为324的一组线   
LATT,1,2,1          设置线的单元属性，材料号为1，实常数号为2，单元类型号为1
LESIZE,ALL,,,5      对所选的每条线的分段数为5
LMESH,ALL            在所有的线上生成线单元
梁网格划分
LSEL,S,,,1,180      从所有线中选择最小线号为1，最大线号为180的一组线
LATT,1,3,1,        设置线的单元属性，材料号为1，实常数号为3，单元类型号为1
LESIZE,ALL,,,8      对所选的每条线的分段数为8
LMESH,ALL          在所有的线上生成线单元
剪力墙网格划分
外剪力墙两端、内剪力墙顶端
LSEL,S,,,372          从所有线中选择编号为372的线为当前子集   
LSEL,A,,,413          另选线编号为413的线添加到当前子集中
LSEL,A,,,461          另选线编号为461的线添加到当前子集中
LESIZE,ALL,,,8        对所选的每条线的分段数为8
内剪力墙底端
LSEL,S,,,418          从所有线中选择编号为418的线为当前子集   
LSEL,A,,,420          另选线编号为420的线添加到当前子集中
LESIZE,ALL,,,3        对所选的每条线的分段数为3
内剪力墙底端内侧
LSEL,S,,,325          从所有线中选择编号为325的线为当前子集 
LSEL,A,,,326          另选线编号为326的线添加到当前子集中               
LESIZE,ALL,,,6        对所选的每条线的分段数为6
剪力墙空洞顶部
LSEL,S,,,375          从所有线中选择编号为375的线为当前子集
LSEL,A,,,379          另选线编号为379的线添加到当前子集中
LSEL,A,,,383
LSEL,A,,,387
LSEL,A,,,391
LSEL,A,,,395
LSEL,A,,,399
LSEL,A,,,403
LSEL,A,,,407
LSEL,A,,,411
LSEL,A,,,415
LSEL,A,,,419
LSEL,A,,,423
LSEL,A,,,427
LSEL,A,,,431
LSEL,A,,,435
LSEL,A,,,439
LSEL,A,,,443
LSEL,A,,,447
LSEL,A,,,451
LSEL,A,,,455
LSEL,A,,,459
LSEL,A,,,463
LESIZE,ALL,,,2          对所选的每条线的分段数为2
ASEL,S,,,1,24            从所有面中选择最小编号为1，最大编号为24的面为当前子集
AATT,1,4,2              设置面的单元属性，材料号为1，实常数号为4，单元类型号为2
AMESH,ALL                在所有的面上生成面单元
楼板网格划分
ASEL,S,,,25,96        从所有面中选择最小编号为25，最大编号为96的面为当前子集
AATT,1,5,2            设置面的单元属性，材料号为1，实常数号为5，单元类型号为2
AMESH,ALL              在所有的面上生成面单元
NSEL,S,LOC,Z,0        选择Z向为0的节点（即选取模型底端节点）         
D,ALL,ALL            对所选节点施加自由度约束           
ALLSEL,ALL            选中所有项目         
/ESHAPE,1            显示由实常数或截面定义的单元形状（p95)
/VIEW,1,1,1,1        设置视图方向，对窗口号1进行视图设置，视图方向为（1,1,1）与原点（0,0,0）连线的方向
/ANGLE,1,270,XM,0    对窗口号1以旋转轴(总体直角坐标系的XM)旋转270度（顺时针为正）,0为采用绝对角度旋转（可缺省）
/REPLOT              显示各实常数号             
FINISH                退出前处理器
第二步：框架结构施加重力荷载
/SOLU                进入求解器
ANTYPE,STATIC        （P262)指定分析类型为静态分析
NSEL,S,LOC,Z,0        选择Z向为0的节点（即选取模型底端节点）       
D,ALL,ALL            对所选节点施加自由度约束           
ALLSEL,ALL            选中所有项目         
ACEL,0,0,9.8        （P248)施加的重力加速度为9.8
SOLVE                      求解
FINISH                      退出求解器
/
POST1                进入后处理器
SET,FIRST            (P295）从结果文件中读出指定的数据组，此为读出第1个载荷步的数据
PLNSOL,U,Z,0,1      （P300)显示节点结果，仅显示结构z方向的结构变形图，缩放因子为1（缺省）         
FINISH                    退出后处理器
第三步：框架结构施加楼面活荷载D=3KN/M^2
FINISH
/SOLU
ANTYPE,STATIC
NSEL,S,LOC,Z,0        选取底层所有节点，选择Z=0的所有节点       
D,ALL,ALL            约束所有自由度             
SFA,ALL,1,PRES,3000        （P244)在面上施加面荷载，对所选的面施加的面荷载为3000
ALLSEL,ALL                  （P122)选择所有图素           
SOLVE
FINISH
/POST1
SET,FIRST          (P295）从结果文件中读出指定的数据组，此为读出第1个载荷步的数据
PLNSOL,U,Z,0,1    （P300)显示节点结果，仅显示结构z方向的结构变形图，缩放因子为1（缺省）
FINISH
第四步：框架结构施加风荷载(基本风压=0.25KN/M^2，体形系数为1.0）
框架结构施加风荷载(先Y方向，后X方向）分两次分别施加
FINISH
/SOLU
ANTYPE,STATIC        （P262)指定分析类型为静态分析
*DIM,LOAD_1,ARRAY,12  定义数值数组名为LOAD_1的数组为12行，1列
LOAD_1(1)=3.78,2.16,2.39,2.57,2.72,2.84,2.95,3.17,3.20,3.29,3.39,3.51
*DIM,LOAD_2,ARRAY,12    定义数值数组名为LOAD_2的
数组为12行，1列
LOAD_2(1)=7.56,4.32,4.78,5.14,5.44,5.68,5.90,6.34,6.40,6.58,6.76,7.02
*DIM,LOAD_3,ARRAY,12    定义数值数组名为LOAD_3的数组为12行，1列
LOAD_3(1)=6.93,3.96,4.38,4.71,4.99,5.21,5.41,5.81,5.87,6.03,6.20,6.44
*DIM,LOAD_4,ARRAY,12    定义数值数组名为LOAD_4的数组为12行，1列
LOAD_4(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93
*DIM,LOAD_A,ARRAY,12    定义数值数组名为LOAD_A的数组为12行，1列
LOAD_A(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93
*DIM,LOAD_B,ARRAY,12    定义数值数组名为LOAD_B的数组为12行，1列
LOAD_B(1)=6.30,3.60,3.98,4.28,4.54,4.74,4.92,5.28,5.34,5.48,5.64,5.86   
*DIM,LOAD_C,ARRAY,12    定义数值数组名为LOAD_C的数组为12行，1列
LOAD_C(1)=3.15,1.80,1.99,2.14,2.27,2.37,2.46,2.64,2.67,2.74,2.82,2.93 
*DO,II,1,12                  利用循环，循环变量从1~12
FK,1+20*II,FY,LOAD_1(II)    （P236）施加关键点集中荷载，即对关键点1+20*II施加Y方向的集中荷载值为LOAD_1(II)
*ENDDO                        结束循环
*DO,II,1,12
FK,2+20*II,FY,LOAD_2(II)
*ENDDO
*DO,II,1,12
FK,3+20*II,FY,LOAD_3(II)
*ENDDO
*DO,II,1,12
FK,4+20*II,FY,LOAD_4(II)
*ENDDO 
SOLVE
FKDELE,ALL,ALL  （P236)删除所有的关键点上的集中荷载
第二次施加风荷载
*DO,II,1,12
FK,1+20*II,FX,LOAD_A(II)
*ENDDO 
*DO,II,1,12
FK,5+20*II,FX,LOAD_B(II)
*ENDDO       
*DO,II,1,12
FK,9+20*II,FX,LOAD_C(II)
*ENDDO 
SOLVE
FINISH
/POST1                  进入后处理器
SET,FIRST            （P295)从结果文件中读出第1个载荷步的数据
PLNSOL,U,Y,0,1        （P300)显示节点结果，即显示Y方向的位移云图
SET,NEXT                    从结果文件中读出下一个个载荷步的数据
PLNSOL,U,X,0,1              显示节点结果，即显示x方向的位移云图
FINISH                  退出后处理器
第五步：框架结构荷载组合计算，框架结构荷载组合计算：(1.D+L+WIND_X;2.D+L+WIND_Y)
第六步：框架结构模态
FINISH
/SOLU                  进入求解器
ANTYPE,MODAL          （P262)定义分析类型，定义分析类型为模态分析         
MODOPT,LANB,6          （P273)定义模态分析方法，用分块兰索斯法提取的模态数目为6
MODOPT,LANB,6,0,0,,OFF                      用分块兰索斯法提取的模态数目为6，程序从0开始计算，相对质量矩阵进行归一化处理
SOLVE
FINISH
/POST1                进入后处理器
SET,LIST              (P295)从结果文件中读出指定的数据组，列表显示荷载步的汇总
SET,FIRST                                                  读出第1个载荷步的数据
PLDISP,0              （P300)显示结构变形图，0为仅显示结构变形图
SET,NEXT                                                  读出第1个载荷

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。