ansys桩土分析命令流
/config,nres,2000000
/units,si!国际单位制
fk1=1!法向接触刚度
ftoln=0.1!初始渗透因子
/prep7
et,1,solid45
A=3.14*0.3*0.3
cylind,0,0.192,0,11,0,90!旋喷桩上部内侧v1
cylind,0.192,0.2,0,11,0,90!钢管v2
cylind,0.2,0.3,0,11,0,90!旋喷桩上部外侧v3
cylind,0,0.3,11,22,0,90!旋喷桩下部v4
cylind,0.3,1.5,0,2.6,0,90!2.6m淤泥质粉质粘土v5
cylind,0.3,1.5,2.6,7.92,0,90!5.32m粉土v6
cylind,0.3,1.5,7.92,16.3,0,90!8.38m淤泥v7
cylind,0.3,1.5,16.3,18.67,0,90!2.37m粉质粘土v8
cylind,0.3,1.5,18.67,20.47,0,90!1.8m粉质粘土v9
cylind,0.3,1.5,20.47,25.3,0,90!4.83m粉土
v10cylind,0,0.3,22,42,0,90!桩底内侧20m粉质粘土
v11cylind,0.3,1.5,25.3,42,0,90!桩底外侧20m粉质粘土v12
allsel,al
lvsel,s,,
,2
mp,ex,2,2e11!钢管
mp,nuxy,2,0.3
mp,dens,2,7850
allsel,all
vsel,s,,,1,4
vsel,u,,,2
mp,ex,3,1.5e8
mp,nuxy,3,0.2mp,
dens,3,2131
!tb,concr,3!水泥土破坏模型，材料属性
!tbdata,,0.7,0.9,0.3e6,2e6!张开裂缝的剪力传递系数,闭合裂缝的剪力传递系数, 单轴抗拉强度,单轴抗压强度
allsel,all
vsel,s,,,5
vplot
mp,ex,5,3e6!淤泥质粉质黏土
mp,nuxy,5,0.35
mp,dens,5,1950
!mp,mu,5,0.633!设置摩擦系数为0.633
tb,dp,5!土层5的D-P参数
tbdata,1,13.2e3,10.1,0!土的粘聚力，摩擦角，膨胀角
allsel,al
lvsel,s,,
,6
vplot
mp,ex,6,9e6!粉土
mp,nuxy,6,0.3
mp,dens,6,1970
!mp,mu,6,0.633!设置摩擦系数为0.633
tb,dp,6!土层6的D-P参数
tbdata,1,15.2e3,21.2,30!土的粘聚力，摩擦角，膨胀角
allsel,all
vsel,s,,,7
vplot
mp,ex,7,1.5e6!淤泥
mp,nuxy,7,0.4
mp,dens,7,1670
!mp,mu,6,0.633!设置摩擦系数为0.633
tb,dp,7!土层7的D-P参数
tbdata,1,12.9e3,9,0!土的粘聚力，摩擦角，膨胀角
allsel,all
vsel,s,,,8
vplot
mp,ex,8,5e6!粉质粘土
mp,nuxy,8,0.35
mp,dens,8,1880
!mp,mu,8,0.633!设置摩擦系数为0.633
tb,dp,8!土层8的D-P参数
tbdata,1,15.8e3,15.2,0!土的粘聚力，摩擦角，膨胀角
allsel,all
vsel,s,,,9
vplot
mp,ex,9,8e6!粉质粘土
mp,nuxy,9,0.35
mp,dens,9,2000
!mp,mu,9,0.633!设置摩擦系数为0.633
tb,dp,9!土层9的D-P参数
tbdata,1,37.1e3,19.8,0!土的粘聚力，摩擦角，膨胀角
allsel,all
vsel,s,,,10
vplot
mp,ex,10,10e6!粉土
mp,nuxy,10,0.3
mp,dens,10,1980
!mp,mu,10,0.633!设置摩擦系数为0.633
tb,dp,10!土层10的D-P参数
tbdata,1,14.3e3,20.1,10!土的粘聚力，摩擦角，膨胀角
allsel,all
vsel,s,,,11,12
vplot
mp,ex,11,6e6!粉质粘土
mp,nuxy,11,0.35
mp,dens,11,1910
!mp,mu,10,0.633!设置摩擦系数为0.633
tb,dp,11!土层11的D-P参数
tbdata,1,23.4e3,15.9,0!土的粘聚力，摩擦角，膨胀角
!布尔运算粘结体所有的除桩之外的部分
ALLSEL,ALLvsel
,s,,,5,12,1
vplotvglue,all
vsel,s,,,1,4vp
lotvglue,allal
lsel,all/pnum,
line,1
saveallsel,all
allsel,all
lsel,s,,,3,4
lsel,a,,,13,14
lesize,all,,,6
lsel,s,,,7,8
lsel,a,,,18,19
lesize,all,0.1
lsel,s,,,55,58
lesize,all,,,1
type,1
mat,2!钢管
vmesh,6
vsel,s,,,1,7,6 vsel,a,,,8
vplott
ype,1
mat,3!水泥土
lsel,s,,,1,2 lesize,all,,,6 lsel,s,,,55,56 lesize,all,,,1 lsel,s,,,22,24,2le size,all,,,10,0.8, 1
lsel,s,,,3 lesize,all,,,6 lsel,s,,,13 lesize,13,,,6 lsel,s,,,25 lesize,25,,,6 lsel,s,,,7,9 lsel,a,,,18,1
9lsel,a,,,30,
31
lsel,a,,,67,6
9lesize,all,0
.1
lsel,s,,,4,6!体8 lsel,a,,,14,26,12 lesize,all,,,6 lsel,s,,,57,58 lesize,all,,,1 lsel,s,,,65,66 lesize,all,,,10,0. 8,1aadd,35,31,2 aslv,slsla,s FLST,2,3,4,ORDE,3 FITEM,2,6 FITEM,2,57 FITEM,2,65 LCCAT,P51X FLST,2,3,4,ORDE,3 FITEM,2,5FITE
M,2,58
FITEM,2,66LCC AT,P51Xvsweep
,1
vsweep,8
allsel,all
vsel,u,,,1,6,5 vsel,u,,,7,8
vsel,u,,,18
vplot
lsel,s,,,44,46,2 lesize,all,,,6
lsel,s,,,43,45,2l esize,all,,,5,0.8
,1
lsel,s,,,51,54, lesize,all,0.1
lsel,s,,,47,49,2 lesize,all,,,6
lsel,s,,,48,50,2 lesize,all,,,5,0.8,1 type,1
mat,5
vsweep,5
vsel,s,,,13
lsel,s,,,136,139 lesize,all,0.1
lsel,s,,,59,61,2 lesize,all,,,6
lsel,s,,,60,62,2l esize,all,,,5,0.8
,1
type,1
mat,6vsweep
,13
vsel,s,,,14
lsel,s,,,140,143 lesize,all,0.1
lsel,s,,,71,73,2 lesize,all,,,6
lsel,s,,,72,74
,2
lesize,all,,,5,0.8,1 type,1
mat,7
vsweep,14
lsel,s,,,144,147 lesize,all,0.1lse
l,s,,,83,85,2 lesize,all,,,6lse
l,s,,,84,86,2lesi ze,all,,,5,0.8,1 type,1
mat,8
vsweep,15
vsel,s,,,16
lsel,s,,,148,151 lesize,all,0.1
lsel,s,,,95,97,2 lesize,all,,,6
lsel,s,,,96,98,2 lesize,all,,,5,0.8,1 type,1
mat,9
vsweep,16
vsel,s,,,19
lsel,s,,,158,163 lesize,all,0.1
lsel,s,,,107,109,2 lesize,all,,,6
lsel,s,,,108,110,2 lesize,all,,,5,0.8,1 type,1
mat,10vswee
p,19
vsel,s,,,17
lsel,s,,,152,156 lesize,all,0.1
lsel,s,,,118,128,10 lesize,all,,,6
lsel,s,,,155,157,2 lesize,all,,,5,0.8,1 type,1m
at,11
vsweep,17
lsel,s,,,154,158,2
lsel,a,,,159
lsel,a,,,123
lesize,all,0
.1
lsel,s,,,118,120,
lesize,all,,,6
lsel,s,,,115,117
lesize,all,,,6
type,1
mat,11vswee
p,18
!用接触向导建立接触对
et,2,170
et,3,173
keyopt,3,9,0!包括由几何形状和接触表面偏移所引起的穿透keyopt,3,12,2!不分开接触，允许滑动，法向不分开
!桩建立目标面
r,20,,,fk1,ftoln,,,
mp,mu,100,0.85
allsel,all
vsel,s,,,1,6,5
vsel,a,,,7,8
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20
aplot
cm,target,area
type,2
real,20
mat,100
nsla,s,1
esurf,all
esel,s,type,,2
eplot
r,21,,,fk1,ftoln,,,
mp,mu,101,0.3
allsel,all
vsel,s,,,1,6,5
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20 aplot
cm,target,area type,2
real,21
mat,101
nsla,s,1
esurf,all
esel,s,type,,2
eplot
r,22,,,fk1,ftoln,,, mp,mu,102,0.15 allsel,all
vsel,s,,,1,6,5 vsel,a,,,7,8
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20 aplot
cm,target,area type,2
real,22m
at,102
nsla,s,1
esurf,all
esel,s,type,,2
eplot
r,23,,,fk1,ftoln,,, mp,mu,103,0.3 allsel,all
vsel,s,,,1,6,5 vsel,a,,,7,8
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20 aplot
WORD格式
cm,target,area
type,2
real,23m
at,103
nsla,s,1
esurf,al
l
esel,s,type,,2
eplot
r,24,,,fk1,ftoln,,,
mp,mu,104,0.35
allsel,all
vsel,s,,,1,6,5
vsel,a,,,7,8
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20
aplot
cm,target,area
type,2r
eal,24
mat,104
nsla,s,1
esurf,all
esel,s,type,,2
eplot
r,25,,,fk1,ftoln,,,
mp,mu,105,0.3
allsel,all
vsel,s,,,1,6,5
vsel,a,,,7,8
aslv,s
aplot
asel,s,,,14
asel,a,,,19,39,20
aplot
cm,target,area
type,2
real,25m
at,105
nsla,s,1
esurf,all
WORD格式
esel,s,type,,2
eplot
!桩建立接触面
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,26
asel,a,,,72
asel,a,,,76
asel,a,,,80
asel,a,,,84
asel,a,,,95
asel,a,,,59
aplot
type,3
real,20
mat,100
nsla,s,1
esurf,all
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,72
aplott
ype,3
real,21
mp,mu,101,0.3
mat,101
nsla,s,1
esurf,all
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,76
aplott
ype,3
real,22
mp,mu,102,0.15
WORD格式
mat,102
nsla,s,1
esurf,all
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,80
aplot
type,3r
eal,23
mp,mu,103,0.3
mat,103
nsla,s,1
esurf,all
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,84
aplot
type,3
real,24
mp,mu,104,0.35
mat,104
nsla,s,1
esurf,all
allsel,all
vsel,u,,,1,6,5
vsel,u,,,7,8
vplot
asel,s,,,59,95,36
aplot
type,3
real,25
mp,mu,105,0.3
mat,105n
sla,s,1
esurf,al
l
/psymb,esys,1
allsel
gplot
savef
inish
!求解过程
/solu
!定义面约束条件
asel,s,loc,x,0!侧面施加对称约束
da,all,symm
asel,s,loc,y,0
da,all,symm
asel,s,loc,z,42底!面试加约束
da,all,all
allsel,all!土外层施加约束
asel,s,,,86asel,a,,,94
asel,a,,,8
2asel,a,,,
78
asel,a,,,7
4asel,a,,,
70asel,a,,
,25aplot
da,all,all
ALLSEL,ALL
antype,static
nlgeom,onsolc
ontrol,on
autot,onpred,
on
neqit,50!每个子步允许的最大迭代次数outpr,all,all
OUTRES,all,all
FirstF=840e3/A
LastF=840e3/A!考虑注浆会增加30％的部分Increment=60e3/A
*do,Force,FirstF,LastF,Increment allsel,al
lvsel,s,,
,1
vsel,a,,,6
vsel,a,,,7
aslv,s,
asel,s,,,1,
asel,a,,,2
9asel,a,,,
34
aplot
sfa,all,,pres,Force allsel,all
nsubst,10,50,5 solve
*enddo
finish
ansys桩土分析命令流
/config,nres,2000000 /units,si !国际单位制 fk1=1!法向接触刚度 ftoln=0.1!初始渗透因子 /prep7 et,1,solid45 A=3.14*0.3*0.3 cylind,0,0.192,0,11,0,90 !旋喷桩上部内侧v1 cylind,0.192,0. 2,0,11,0,90!钢管v2 cylind,0.2,0.3,0,11,0,90 !旋喷桩上部外侧v3 cylind,0,0.3,11,22,0,90 !旋喷桩下部v4 cylind,0.3,1.5,0,2.6,0,90 !2.6m淤泥质粉质粘土v5 cylind,0.3,1.5,2.6,7.92,0,90 !5.32m粉土v6 cylind,0.3,1.5,7.92,16.3,0,90 !8.38m淤泥v7 cylind,0.3,1.5,16.3,18.67,0,90!2.37m粉质粘土v8 cylind,0.3,1.5, 18.67,20.47 ,0,90!1.8m粉质粘土v9 cylind,0.3,1.5, 20.47,25.3 ,0,90!4.83m粉土v10 cylind,0,0.3,22,42,0,90 !桩底内侧20m粉质粘土v11 cylind,0.3,1.5,25.3,42,0,90 !桩底外侧20m粉质粘土v12 allsel,all vsel,s,,,2 mp,ex,2,2e11 !钢管 mp,nuxy,2,0.3 mp,dens,2,7850 allsel,all vsel,s,,,1,4 vsel,u,,,2 mp,ex,3,1.5e8 mp,nuxy,3,0.2 mp,dens,3,2131 !tb,concr,3 !水泥土破坏模型，材料属性 !tbdata,,0.7,0.9,0.3e6,2e6!张开裂缝的剪力传递系数 ,闭合裂缝的剪力传递系数 ,单轴抗拉强度, 单轴抗压强度 allsel,all vsel,s,,,5 vplot mp,ex,5,3e6 !淤泥质粉质黏土 mp,nuxy,5,0.35 mp,dens,5,1950 !mp,mu,5,0.633 !设置摩擦系数为0.633 tb,dp,5 !土层5的D-P参数
ANSYS命令流实例
/PREP7 !进入前处理 ANTYPE,STATIC !设置分析类型为静力结构分析 PSTRES,ON !用于后面的模态分析中考虑预应力(该开关不影响静力分析) ET,1,LINK10 !选取单元类型1(单向杆单元) KEYOPT,1,3,0 !设置仅承受拉应力，KEYOPT(3)=0 R,1,306796E-8,543248E-8 !设置实常数，包括绳索截面积(306796E-8)，初始应变(543248E-8) MP,EX,1,30E6 !定义材料的弹性模量(1号材料) MP,DENS,1,73E-5 !定义材料的密度(1号材料) N,1 ! 定义第1号节点 N,14,100 ! 定义第14号节点 FILL ! 均分填满第2号至第13号节点 E,1,2 !由节点1及节点2生成单元 EGEN,13,1,1 !依序复制生成13个单元 D,ALL,ALL ! 对所有节点施加固定约束 FINISH ! 前处理结束 /SOLU ! 进入求解模块，求解预应力引起的应力状态 SOLVE ! 求解 FINISH ! 退出求解模块 /POST1 ! 进入一般的后处理 ETABLE,STRS,LS,1 !针对LINK10单元,建立单元列表STRS,通过LS及特征号1来获得单元的轴向应力 *GET,STRSS,ELEM,13,ETAB,STRS !针对单元列表STRS, 提取13号单元的应力 FINISH ! 后处理结束 /POST26 ! 进入时间历程后处理，处理支反力 RFORCE,2,1,F,X !将1号节点上的x方向支反力提取，并存储到2号变量中 STORE ! 存储 *GET,FORCE,V ARI,2,EXTREM,VMAX !将2号变量的最大值赋给参数FORCE /SOLU ! 再次进入求解模块，模态分析 ANTYPE,MODAL ! 模态分析 MODOPT,SUBSP,3 ! 选择子空间迭代法，求3阶模态 MXPAND,3 ! 设定3阶模态扩展 PSTRES,ON ! 用于在模态分析中考虑预应力(还需在前面的静力分析中也同时打开) DDELE,2,UX,13 ! 删除从2号节点到13号节点上的UX约束 DDELE,2,UY,13 !删除从2号节点到13号节点上的UY约束 SOLVE !求解 *GET,FREQ1,MODE,1,FREQ ! 提取第1阶模态共振频率，并赋值给参数FREQ1 *GET,FREQ2,MODE,2,FREQ ! 提取第2阶模态共振频率，并赋值给参数FREQ2 *GET,FREQ3,MODE,3,FREQ ! 提取第3阶模态共振频率，并赋值给参数FREQ3 *STATUS !列出所有参数的实际内容
ansys旋转经典命令流
1 旋转摩擦 (1) 2. 电磁三d命令流实例（论坛看到） (11) 3. 帮助感应加热例子induction heating of a solid cylinder billet (15) 4. 感应加热温度场的数值模拟（论文）inducheat30命令流 (19) 5. 如何施加恒定的角速度？Simwe仿真论坛 (24) 6. 旋转一个已经生成好的物体 (27) 7. 产生这样的磁力线 (28) 8. 旋转摩擦生热简单例子（二维旋转） (32) 8.1. 原版 (32) 8.2. 部分gui操作 (35) 9. VM229 Input Listing (39) 10 轴承---耦合+接触分析 (47) 11. 板的冲压仿真 (52) 1 旋转摩擦 FINISH /FILNAME,Exercise24 !定义隐式热分析文件名 /PREP7 !进入前处理器 ET,1,SOLID5 !选择单元类型 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7800 !定义材料1的密度 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,C,1,,460 !定义材料1的比热 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,KXX,1,,66.6 !定义材料1的热传导系数 MPTEMP,,,,,,,, MPTEMP,1,0 UIMP,1,REFT,,,30 !定义材料1的热膨胀系数的参考温度 MPDATA,ALPX,1,,1.06e-5 !定义材料1的热膨胀系数MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,206e9 !定义材料1的弹性模量 MPDATA,PRXY,1,,0.3 !定义材料1的泊松比 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,2,,8900 !定义材料2的密度 MPTEMP,,,,,,,, MPTEMP,1,0
ansys实例命令流-弹塑性分析命令流
/FILNAME,Elastic-Plasitc,1 /TITLE, Elastic-Plasitc Analysis !前处理。 /PREP7 !**定义梁单元189。 ET,1,BEAM189 !定义单元。 !**梁截面1。 SECTYPE, 1, BEAM, HREC, , 0 !定义梁截面。SECOFFSET, CENT SECDATA,50,100,6,6,6,6,0,0,0,0 !定义梁截面完成。 !**定义材料。 MPTEMP,,,,,,,, !定义弹塑性材料模型。MPTEMP,1,0 MPDATA,EX,1,,2.05e5 MPDATA,PRXY,1,,0.3 TB,BISO,1,1,2, TBTEMP,0 TBDATA,,150,18600,,,, !定义弹塑性材料模型。!**建立几何模型。 K,1, , , , K,2 ,900, K,3 ,,50 LSTR, 1, 2 !**网格划分。 FLST,5,1,4,ORDE,1 !定义网格密度。FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义网格密度完成。CM,_Y,LINE !网格划分。 LSEL, , , , 1 CM,_Y1,LINE CMSEL,S,_Y CMSEL,S,_Y1 LATT,1, ,1, , 3, ,1 CMSEL,S,_Y CMDELE,_Y CMDELE,_Y1 LMESH, 1 !网格划分完成。 !施加载荷及求解。 FINISH /SOL
!**施加约束。 FLST,2,1,3,ORDE,1 !施加约束。FITEM,2,1 /GO DK,P51X, , , ,0,UX,UY,UZ,ROTX, , , FLST,2,1,3,ORDE,1 FITEM,2,2 /GO DK,P51X, , , ,0,UY,UZ,ROTX, , , , !施加约束完成。 !**加载。 FLST,2,50,2,ORDE,2 FITEM,2,1 FITEM,2,-50 SFBEAM,P51X,1,PRES,100, , , , , , LSWRITE,1, !定义载荷步1完成。FLST,2,50,2,ORDE,2 !定义载荷步2。FITEM,2,1 FITEM,2,-50 SFEDELE,P51X,1,PRES LSWRITE,2, !定义载荷步2完成。!设定求解步并求解。 LSSOLVE,1,2,1,
Ansys常见命令流
Ansys命令流 第一天 目标:熟悉ANSYS基本关键字的含义 k --> Keypoints 关键点 l --> Lines 线 a --> Area 面 v --> Volumes 体 e --> Elements 单元 n --> Nodes 节点 cm --> component 组元 et --> element type 单元类型 mp --> material property 材料属性 r --> real constant 实常数 d --> DOF constraint 约束 f --> Force Load 集中力 sf --> Surface load on nodes 表面载荷 bf --> Body Force on Nodes 体载荷 ic --> Initial Conditions 初始条件 第二天 目标:了解命令流的整体结构,掌握每个模块的标识 !文件说明段 /BATCH /TITILE,test analysis !定义工作标题 /FILENAME,test !定义工作文件名 /PREP7 !进入前处理模块标识 !定义单元,材料属性,实常数段 ET,1,SHELL63 !指定单元类型 ET,2,SOLID45 !指定体单元 MP,EX,1,2E8 !指定弹性模量 MP,PRXY,1,0.3 !输入泊松比 MP,DENS,1,7.8E3 !输入材料密度 R,1,0.001 !指定壳单元实常数-厚度...... !建立模型 K,1,0,0,, !定义关键点 K,2,50,0,,
K,3,50,10,, K,4,10,10,, K,5,10,50,, K,6,0,50,, A,1,2,3,4,5,6, !由关键点生成面 ...... !划分网格 ESIZE,1,0, AMESH,1 ...... FINISH !前处理结束标识 /SOLU !进入求解模块标识 !施加约束和载荷 DL,5,,ALL SFL,3,PRES,1000 SFL,2,PRES,1000 ...... SOLVE !求解标识 FINISH !求解模块结束标识 /POST1 !进入通用后处理器标识 ...... /POST26 !进入时间历程后处理器 …… /EXIT,SAVE !退出并存盘 以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助 /ANGLE !指定绕轴旋转视图 /DIST !说明对视图进行缩放 /DEVICE !设置图例的显示,如:风格,字体等 /REPLOT !重新显示当前图例 /RESET !恢复缺省的图形设置 /VIEW !设置观察方向 /ZOOM !对图形显示窗口的某一区域进行缩放
ANSYS应用——m法计算桩基时模拟桩与土作用的命令流
ANSYS应用——m法计算桩基时模拟桩与土作用的命令流 !z1spring1.mac /PREP7 WPCSYS,-1 !工作平面恢复到默认状态 CSYS,0 dpile1=1.5!桩直径***** allsel,all cmsel,s,z nsll,s,1 NSEL,r,LOC,Z,-5-5-5.36,-6-5-5.36!建立泥面以下桩的节点集合 CM,SPNODE1,NODE allsel,all,all cmsel,s,z nsll,s,1 NSEL,r,LOC,Z,-6-5-5.36,-20-5-5.36 CM,SPNODE2,NODE allsel,all,all cmsel,s,z nsll,s,1 NSEL,r,LOC,Z,-20-5-5.36,-25-5-5.36 CM,SPNODE3,NODE allsel,all,all cmsel,s,z nsll,s,1 NSEL,r,LOC,Z,-25-5-5.36,-27-5-5.36 CM,SPNODE4,NODE allsel,all,all cmsel,s,z nsll,s,1 NSEL,r,LOC,Z,-27-5-5.36,-31-5-5.36 CM,SPNODE5,NODE allsel,all,all !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!spring1 /PREP7 CMSEL,S,SPNODE1 !选择泥面以下,桩的节点集合 *GET,enum,ELEM,,COUNT, , , , *GET, NNUM, NODE, 0, COUNT, *GET, nmin, NODE, 0, num, min,
ANSYS-结构稳态(静力)分析之经典实例-命令流格式
ANSYS 结构稳态(静力)分析之经典实例-命令流格式.txt两人之间的感情就像织毛衣，建立 的时候一针一线，小心而漫长，拆除的时候只要轻轻一拉。。。。/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生 一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT
ansys命令流
第一天目标： 熟悉ANSYS基本关键字的含义k --> Keypoints关键点l --> Lines线a --> Area 面v --> Volumes体e --> Elements单元n --> Nodes节点cm --> component组元et --> element type单元类型mp --> material property材料属性r --> real constant实常数d --> DOF constraint约束f --> Force Load集中力sf --> Surface load on nodes 表面载荷bf --> Body Force on Nodes体载荷ic --> Initial Conditions初始条件第二天目标： 了解命令流的整体结构,掌握每个模块的标识!文件说明段/BATCH/TILE,test analysis!定义工作标题/FILENAME,test!定义工作文件名/PREP7!进入前处理模块标识!定义单元,材料属性,实常数段ET,1,SHELL63!指定单元类型ET,2,SOLID45!指定体单元MP,EX,1,2E8!指定弹性模量MP,PRXY,1, 0.3!输入泊松比MP,DENS,1, 7.8E3!输入材料密度R,1, 0.001!指定壳单元实常数-厚度......!建立模型K,1,0,0,,!定义关键点 K,2,50,0,,K,3,50,10,,K,4,10,10,,K,5,10,50,,K,6,0,50,,A,1,2,3,4,5,6,!由关键点生成面......!划分网格ESIZE,1,0,AMESH, 1......FINISH!前处理结束标识/SOLU!进入求解模块标识!施加约束和载荷DL,5,,ALLSFL,3,PRES,1000SFL,2,PRES, 1000......SOLVE!求解标识FINISH!求解模块结束标识/POST1!进入通用后处理器标识....../POST26!进入时间历程后处理器……/EXIT,SAVE!退出并存盘以下是日志文件中常出现的一些命令的标识说明,希望能给大家在整理LOG文件时有所帮助/ANGLE!指定绕轴旋转视图/DIST!说明对视图进行缩放/DEVICE!设置图例的显示,如： 风格,字体等/REPLOT!重新显示当前图例/RESET!恢复缺省的图形设置/VIEW!设置观察方向/ZOOM!对图形显示窗口的某一区域进行缩放第三天生成关键点和线部分 1.生成关键点K,关键点编号,X坐标,Y坐标,Z坐标例：
ansys实例命令流-谱分析命令流
/FILNAME, Beam,1 ！定义工作文件名。 /TITLE, Beam Analysis ！定义工作标题。/PREP7 !定义单元。 ET,1,BEAM188 !定义材料属性。 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,2.1e5 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,7.9e-6 ! 定义杆件截面■200。 SECTYPE, 1, BEAM, RECT, , 0 SECOFFSET, CENT SECDATA,10,10,0,0,0,0,0,0,0,0 !建立几何模型。 K,1, ,, , K,2,350,, , !生成立柱。 LSTR, 1, 2 !以上完成几何模型。 !以下进行网格划分。 FLST,5,1,4,ORDE,1 FITEM,5,1 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,50, , , , ,1 !定义单元大小。!分配、划分平板结构。 LMESH, 1 !分析类型施加载荷并求解。 FINISH /SOLUTION ANTYPE,2 ！定义分析类型及求解设置。MSAVE,0 ！模态提取方法。
MODOPT,LANB,10 EQSLV,SPAR MXPAND,10, , ,0 ！模态扩展设置。 LUMPM,0 PSTRES,0 MODOPT,LANB,10,0,0, ,OFF MXPAND,10,0,0,1,0.001, !施加约束。 FLST,2,2,3,ORDE,2 FITEM,2,1 FITEM,2,-2 /GO DK,P51X, , , ,0,ALL, , , , , , !求解。 FINISH /SOL /STATUS,SOLU SOLVE !定义谱分析。 FINISH /SOLUTION ANTYPE,8 SPOPT,PSD,10,1 PSDUNIT,1,DISP,386.4, PSDFRQ,1, ,13.8,40,50.6,73,120 !定义谱—频率表。PSDFRQ,1, ,134,178,233, , PSDV AL,1,1,4,0.6,3,5 PSDV AL,1,6,2,6, , FLST,2,2,1,ORDE,2 !施加谱。 FITEM,2,1 FITEM,2,-2 D,P51X,UX,1.0 PFACT,1,BASE, !计算PSD激励参与系数。PSDRES,DISP,REL !设置输出选项。PSDRES,VELO,OFF PSDRES,ACEL,OFF
ansys桩土相互作用例子复习课程
a n s y s 桩土相互 作用 例子 /prep7 et,1,pla ne42 et,2,solid45 mp,ex,1,2.5e10 mp, nu xy,1,0.2 mp,de ns,1,2500 mp,ex,2,2.5e8 mp, nu xy,2,0.4 mp,de ns,2,2000 tb,dp,2 tbdata,1,19,32,30 桩的弹性模量桩的泊松比 桩的密度 土的弹性模量土的泊松比 土的密度 粘聚力c为19摩擦角为32度，膨胀角为30 RECTNG,0,1,0,8, RECTNG,1,5,0,6, 面1 面2 RECTNG,0,5,0,-16,面3 /pnu m,area,1 /pnu m,li ne,1 asel,s,,,2,3,1 aglue,all nu mcmp,all allsel aplot lsel,s,,,1,3,2 划分面1 lesize,all,,,2 lsel,s,,,2,4,2 lesize,all,,,16 amesh,1
lsel,s,,,6,8,2 划分面2 lesize,all,0.5 lsel,s,,,5 lesize,all,0.5 lsel,s,,,7 lesize,all,0.5 amesh,2 lsel,s,,,11 划分面3 lesize,all,0.5 lsel,s,,,12 lesize,all,0.5 lsel,s,,,10 lesize,all,0.5
lsel,s,,,9 lesize,all,0.5 amesh,3 EXTOPT,ESIZE,2,0, EXTOPT,ACLEAR,1 type,2 mat,1 VEXT,1,,,,,-1 ! 拉伸成体1 EXTOPT,ESIZE,2,0, ! 拉伸成体2，3 EXTOPT,ACLEAR,1 mat,2 VEXT,2,3,1,,,-1 allsel /view,1,1,1,1 eplot k,1001,5,6,-5 k,1002,5,0,-5 k,1003,5,-16,-5 l,18,1001 l,1001,1002 l,1002,17 l,1002,1003 l,1003,21 al,22,36,37,38 al,30,38,39,40 面! 4 和面5 EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1 VEXT,19,20,1,-5 !拉伸成体4 和5 vplot vsel,s,,,4,5,1 vplot lsel,s,,,36,52,1 lesize,all,0.5 type,2 mat,2 vmesh,all allsel !划分体4和 5 vsel,s,,,2,5,1 VPLOT !合并体2，3，4，5 上的重复单元及节点
几个ansys经典实例(长见识)
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
ansys命令流最全详细介绍二
三 生成关键点和线部分 1.生成关键点 K,关键点编号,X坐标,Y坐标,Z坐标 例:K,1,0,0,0 2.在激活坐标系生成直线 LSTR,关键点P1,关键点P2 例LSTR,1,2 3.在两个关键点之间连线 L,关键点P1,关键点P2 例L,1,2 注:此命令会随当前的激活坐标系不同而生成直线或弧线 4.由三个关键点生成弧线 LARC,关键点P1,关键点P2,关键点PC,半径RAD 例LARC,1,3,2,0.05 注:关键点PC是用来控制弧线的凹向 5.通过圆心半径生成圆弧
CIRCLE,关键点圆心,半径RAD,,,,圆弧段数NSEG 例:CIRCLE,1,0.05,,,,4 6.通过关键点生成样条线 BSPLIN,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6 例:BSPLIN,1,2,3,4,5,6 7.生成倒角线 LFILLT,线NL1,线NL2,倒角半径RAD 例LFILLT,1,2,0.005 8.通过关键点生成面 A,关键点P1,关键点P2,关键点P3,关键点P4,关键点P5,关键点P6,P7,P8... 例:A,1,2,3,4 9.通过线生成面 AL,线L1,线L2,线L3,线L4,线L5,线L6,线L7,线L8,线L9,线L10 例:AL,5,6,7,8 10.通过线的滑移生成面
ASKIN,线NL1,线NL2,线NL3,线NL4,线NL5,线NL6,线NL7,线NL8,线NL9 例:ASKIN,1,4,5,6,7,8 注:线1为滑移的导向线 四 目标:掌握常用的实体-面的生成 生成矩形面 1.通过矩形角上定位点生成面 BLC4,定位点X方向坐标XCORNER,定位点Y方向坐标YCORNER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 例:BLC4,0,0,5,3,0 2.通过矩形中心定位点生成面 BLC5,定位点X方向坐标XCENTER,定位点Y方向坐标YCENTER,矩形宽度WIDTH,矩形高度HEIGHT,矩形深度DEPTH 注:与上条命令的不同就在于矩形的定位点不一样 例:BLC5,2.5,1.5,5,3,0 3.通过在工作平面定义矩形X.Y坐标生成面 RECTNG,矩形左边界X坐标X1,矩形右边界X坐标X2,矩形下边界Y
ANSYSAPDL命令流悬臂梁分析教程
ANSYS APDL命令流悬臂梁分析教程 本文通过分析悬臂梁介绍了ANSYS APDL相关命令流方法。 考虑悬臂梁如图2-2，求x=L变形量。已知条件：杨氏系数E=200E9;截面参数：t=0.01m, w=0.03m, A=3E-4,I=2.5E-9;几何参数：L=4m, a=2m, b=2m;边界外力F=2N,q=0.05N/m. 使用ANSYS解决该问题的命令如下： /FILNAM,EX2-1? ! 定义文件名 /TITLE,CANTILEVER BEAM DEFLECTION? !定义分析的标题 /UNITS,SI !定义单位制(注意观察输出窗口的单位) /PREP7 !进入前置处理 ET,1,3? !定义元素类型为beam3 MP,EX,1,200E9 ! 定义杨氏模量 R,1,3E-4,2.5E-9,0.01 !定义实常数(要严格根据该元素类型的说明文档所给出的实常数格式) N,1,0,0!定义第1号节点X坐标为0，Y坐标为0 N,2,1,0!定义第2号节点X坐标为1，Y坐标为0 N,3,2,0 !定义第3号节点X坐标为2，Y坐标为0 N,4,3,0 !定义第4号节点X坐标为3，Y坐标为0 N,5,4,0!定义第5号节点X坐标为4，Y坐标为0 E,1,2!把1、2号节点相连构成单元，系统将自定义为1号单元 E,2,3!把2、3号节点相连构成单元，系统将自定义为2号单元 E,3,4!把3、4号节点相连构成单元，系统将自定义为3号单元 E,4,5!把4、5号节点相连构成单元，系统将自定义为4号单元 FINISH? !退出该处理层 /SOLU!进入求解处理器 D,1,ALL,0 !对1节点施加约束使它X，Y向位移都为0
ANSYS建模apdl命令流实例应用
大桥全长2996.8m，其中主桥采用跨度为101.5+188.5+580+217.5+159.5+116m的钢桁梁斜拉桥；非通航孔正桥采用6孔跨径64m预应力混凝土简支箱梁；东引桥采用16孔梁长32.6m预应力混凝土简支箱梁；跨大堤桥采用48.9+86+48.8m预应力混凝土连续箱梁；西引桥采用15孔梁长32.6m预应力混凝土简支梁及2孔梁长24.6m预应力混凝土简支梁，其中宁安线采用箱梁，阜景线采用T梁。 主桥采用103+188.5+580+217.5+159.5+117.5m两塔钢桁斜拉桥方案，全长1366m。主梁为三片主桁钢桁梁，桁间距2x14m，节间长14.5m，桁高15m。主塔为钢筋混凝土结构，塔顶高程+204.00m，塔底高程-6.00m，斜拉索为空间三索面，立面上每塔两侧共18对索，全桥216根斜拉索。所有桥墩上均设竖向和横向约束，4#塔与主梁之间设纵向水平约束，3#塔与梁间使用带限位功能的粘滞阻尼器。主梁为”N”字型桁式，横向采用三片桁结构，主桁的横向中心距各为14m，桁高15m，节间距14.5m[2]。 结构构造 主桥采用两塔钢桁斜拉桥方案，主梁为三片主桁钢桁梁，主桁上下弦杆均为箱型截面，上弦杆内高1000mm，内宽1200mm，板厚20～48mm。下弦杆内高1400mm，宽1200mm，板厚20～56mm。下弦杆顶板向桁内侧加宽700mm与整体桥面板焊接。腹杆主要采用H型截面。H型杆件宽1200mm，高720和760mm，板厚20～48mm。根据不同的受力区段选用不同的杆件截面，在辅助墩附近的压重区梁段，腹杆采用箱型截面杆件。主桁采用焊接杆件，整体节点。在节点外以高强度螺栓拼接的结构形式，上下弦杆四面等强对接拼装。H型腹杆采用插入式连接。箱型腹杆采用四面与主桁节点对拼的连接形式。主桁拼接采用M30高强螺栓。
ansys实例命令流-实体梁分析命令流
/FILNAME,SolidBeam ,1 !定义工作文件名。/TITLE,SolidBeam Analysis !定义工作标题。/PREP7 ET,1,SOLID95 !定义材料属性。 MPTEMP,,,,,,,, !定义材料属性。MPTEMP,1,0 MPDATA,EX,1,,2.06e5 MPDATA,PRXY,1,,0.3 !建立几何模型 K,1,,,, K,2,450,,, K,3,450,-55,, K,4,,-100,, FLST,2,4,3 FITEM,2,1 FITEM,2,2 FITEM,2,3 FITEM,2,4 A,P51X VOFFST,1,45, , !网格划分。 FLST,5,4,4,ORDE,2 FITEM,5,9 FITEM,5,-12 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,4,4,ORDE,4
FITEM,5,2 FITEM,5,4 FITEM,5,6 FITEM,5,8 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,6, , , , ,1 !* FLST,5,4,4,ORDE,4 FITEM,5,1 FITEM,5,3 FITEM,5,5 FITEM,5,7 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y LESIZE,_Y1, , ,30, , , , ,1 !* CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,'VOLU' CMSEL,S,_Y !* VSWEEP,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加载。
ANSYS_结构稳态(静力)分析之经典实例-命令流格式
/FILNAME,Allen-wrench,1 ! Jobname to use for all subsequent files /TITLE,Static analysis of an Allen wrench /UNITS,SI ! Reminder that the SI system of units is used /SHOW ! Specify graphics driver for interactive run; for batch ! run plots are written to pm02.grph ! Define parameters for future use EXX=2.07E11 ! Young's modulus (2.07E11 Pa = 30E6 psi) W_HEX=.01 ! Width of hex across flats (.01m=.39in) *AFUN,DEG ! Units for angular parametric functions定义弧度单位 W_FLAT=W_HEX*TAN(30) ! Width of flat L_SHANK=.075 ! Length of shank (short end) (.075m=3.0in) L_HANDLE=.2 ! Length of handle (long end) (.2m=7.9 in) BENDRAD=.01 ! Bend radius of Allen wrench (.01m=.39 in) L_ELEM=.0075 ! Element length (.0075 m = .30 in) NO_D_HEX=2 ! Number of divisions on hex flat TOL=25E-6 ! Tolerance for selecting nodes (25e-6 m = .001 in) /PREP7 ET,1,SOLID45 ! 3维实体结构单元；Eight-node brick element ET,2,PLANE42 ! 2维平面结构；Four-node quadrilateral (for area mesh) MP,EX,1,EXX ! Young's modulus for material 1；杨氏模量 MP,PRXY,1,0.3 ! Poisson's ratio for material 1；泊松比 RPOLY,6,W_FLAT ! Hexagonal area创建规则的多边形 K,7 ! Keypoint at (0,0,0) K,8,,,-L_SHANK ! Keypoint at shank-handle intersection K,9,,L_HANDLE,-L_SHANK ! Keypoint at end of handle L,4,1 ! Line through middle of hex shape L,7,8 ! Line along middle of shank L,8,9 ! Line along handle LFILLT,8,9,BENDRAD ! Line along bend radius between shank and handle! 产生一个倒角圆，并生成三个点 /VIEW,,1,1,1 ! Isometric view in window 1 /ANGLE,,90,XM ! Rotates model 90 degrees about X! 不用累积的旋转 /TRIAD,ltop /PNUM,LINE,1 ! Line numbers turned on LPLOT ! Line numbers off !
ansys桩土相互作用例子
et,1,plane42 et,2,solid45 mp,ex,1,2.5e10 !桩的弹性模量 mp,nuxy,1,0.2 !桩的泊松比 mp,dens,1,2500 !桩的密度 mp,ex,2,2.5e8 !土的弹性模量 mp,nuxy,2,0.4 !土的泊松比 mp,dens,2,2000 !土的密度 tb,dp,2 tbdata,1,19,32,30 !粘聚力c为19,摩擦角为32度,膨胀角为30 RECTNG,0,1,0,8, !面1 RECTNG,1,5,0,6, !面2 RECTNG,0,5,0,-16, !面3 /pnum,area,1 /pnum,line,1 asel,s,,,2,3,1 aglue,all numcmp,all allsel aplot lsel,s,,,1,3,2 !划分面1 lesize,all,,,2 lsel,s,,,2,4,2 lesize,all,,,16 amesh,1 lsel,s,,,6,8,2 !划分面2 lesize,all,0.5 lsel,s,,,5 lesize,all,0.5 lsel,s,,,7 lesize,all,0.5 amesh,2 lsel,s,,,11 !划分面3 lesize,all,0.5
lesize,all,0.5 lsel,s,,,10 lesize,all,0.5 lsel,s,,,9 lesize,all,0.5 amesh,3 EXTOPT,ESIZE,2,0, !拉伸成体1 EXTOPT,ACLEAR,1 type,2 mat,1 VEXT,1,,,,,-1 EXTOPT,ESIZE,2,0, !拉伸成体2，3 EXTOPT,ACLEAR,1 mat,2 VEXT,2,3,1,,,-1 allsel /view,1,1,1,1 eplot k,1001,5,6,-5 k,1002,5,0,-5 k,1003,5,-16,-5 l,18,1001 l,1001,1002 l,1002,17 l,1002,1003 l,1003,21 al,22,36,37,38 !面4和面5 al,30,38,39,40 EXTOPT,ESIZE,10,0, EXTOPT,ACLEAR,1 VEXT,19,20,1,-5 !拉伸成体4和5 vplot vsel,s,,,4,5,1 vplot
命令流ansys经典实例
ansys钢筋混凝土建模实例 finish /clear /prep7 et,1,solid65,,,,,,,1 et,2,link8 et,3,185 et,4,solid45 !************************定义材料属性***************************** !混凝土材料属性 mp,ex,1,1.596e10 mp,prxy,1,0.2 mp,dens,1,2400 fc=1.68e7 !c30混凝土轴心抗压强度设计值 ft=1.86e6 tb,concr,1 tbdata,,0.5,0.95,ft,-1 tb,miso,1,,11 tbpt,,0.0002,fc*0.19 tbpt,,0.0004,fc*0.36 tbpt,,0.0006,fc*0.51 tbpt,,0.0008,fc*0.64 tbpt,,0.0010,fc*0.75 tbpt,,0.0012,fc*0.84 tbpt,,0.0014,fc*0.91 tbpt,,0.0016,fc*0.96 tbpt,,0.0018,fc*0.99 tbpt,,0.002,fc tbpt,,0.0033,fc !57到68 v,6,7,56,54,10,11,60,58 vgen,4,57,57,1,0,0.38,0 vgen,2,57,60,1,0.25,0,0 vgen,2,57,60,1,-0.25,0,0 !69到80 v,189,190,234,233,191,192,238,237 vgen,4,69,69,1,0,0.38,0
vgen,3,69,72,1,0,0,0.25 !121到136 vgen,2,105,120,1,0,0,0.107 !193到208 k,1178,0.81,0,0.107 k,1179,0.823,0,0.107 k,1180,0.81,0,0.417 k,1181,0.81,0.16,0.107 k,1182,0.823,0.16,0.107 k,1183,0.81,0.16,0.417 k,1184,0.81,0.38,0.107 k,1185,0.823,0.38,0.107 k,1186,0.81,0.38,0.417 v,1178,1179,187,1180,1181,1182,189,1183 v,1181,1182,189,1183,1184,1185,191,1186 vgen,4,193,194,1,0,0.38,0 vgen,2,193,200,1,0.107,0,0 !209到232 k,1271,-0.02,0,-0.013 k,1272,-0.02,0,0 k,1273,-0.02,0.16,-0.013 k,1274,-0.02,0.16,0 k,1275,-0.02,0.38,-0.013 k,1276,-0.02,0.38,0 v,677,1,1272,1271,679,5,1274,1273 v,679,5,1274,1273,681,9,1276,1275 vgen,4,209,210,1,0,0.38,0 vgen,3,209,216,1,0.25,0,0 !233到256 vgen,2,209,232,1,0,0,0.107 vgen,2,257,280,1,0.107,0,0 k,5001,-0.1,1.3,0 k,5002,-0.02,1.3,0 k,5003,-0.02,1.3,0.094