急！急！急！有哪位大虾能指点一下 我的ANSYS程序
浏览：60561
我正在做毕业论文 关于桥梁横向地震反应分析 主梁和桥墩用BEAM189 水平向支座（用铅芯橡胶支座）用COMBIN40单元 竖向支座用COMBIN14单元 但是 不管我怎么改变COMBIN40单元的参数 即K C 值 桥梁的地震反应都不会变化 这是怎么回事 ？怎么处理此问题？？？？？？？？？？？
fini
/clear
/prep7
k,1
K,2,3
K,3,3,3.95
K,4,6.4,3.95
K,5,6.4,4.2
K,6,0,4.2
a,1,2,3,4,5,6
blc4,0,0.4,2.6,3.55
asba,1,2
wprot,0,0,90
wpoff,0,0,2.6
asbw,all
wpoff,0,0,0.4
asbw,all
wprot,0,90,0
wpoff,0,0,-0.4
asbw,all
wpoff,0,0,-3.55
asbw,all
wprot,0,-90,0
arsym,x,all
wpcsys
aglue,all
et,1,plane82
lesize,all,,,1
amesh,all
secwrite,mybeam!主梁截面定义 单箱单室 高4。2米，顶板厚0.25米，底板厚0。4米，肋宽0.4米
fini
/clear
/prep7
ET,1,beam189
MP,EX,1,3.0e10
MP,prxy,1,0.1667
MP,DAMP,1,0.05
MP,dens,1,2.549e3
sectype,1,beam,mesh
secread,mybeam,,,mesh
k,1,-32,25.4
k,33,0,25.4
k,73,40,25.4
k,105,72,25.4!定义生成主梁关键点 32米+40米+32米
k,1000,50,50
k,34,0,25
k,35,0,0
k,74,40,25
k,75,40,0 !定义生成桥墩关键点 墩高25米
n,36,-0.25,25
n,77,39.75,25!定义节点 以和桥顶节点 生成COMBIN40单元
l,1,33
l,33,73
l,73,105
lesize,all,0.5
latt,1,,1,,1000,,1
lmesh,all!主梁网格划分
SECTYPE,2,BEAM,RECT
SECDATA,3,6,6,12
l,34,35,25
l,74,75,25
latt,1,,1,,,,2
lmesh,all!桥墩网格划分
/eshape,1
eplot
!ksel,s,loc,x,-32
!dk,all,uy,uz,ux,rotx
!ksel,s,loc,x,72
!dk,all,uy,uz,rotx
nsel,all
d,78,uy,uz,rotx
d,511,uy,uz,rotx!约束梁端处 78号节点为主梁最左端点 511号节点为主梁最右端点
k1=6.3e5
k2=7e4
c=2e5
fd=4e4
ET,3,COMBIN40 !用COMBIN40单元模拟铅芯支座
KEYOPT,3,1,1
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,6,0
R,1,k1,c,1,,fd,k2
type,3
real,1
e,36,703
e,77,754!生成COMBIN40单元 即水平向的铅芯橡胶支座
et,4,combin14
keyopt,4,1,0
keyopt,4,2,0
keyopt,4,3,2
R,2,5000000,5000000,,
type,4
real,2
e,79,703
e,271,754!生成COMBIN14单元
!CPINTF,all
cp,1,ux,36,79
cp,2,ux,77,271 !顺桥向耦合
!定义和读入时程曲线
NT=1000!时程曲线有NT个点
DT=0.01!时间间隔
*dim,ac,,NT
!读入数据
/input,tianjin,txt !天津波共10 秒
/solu
NSUBST,1, , ,1 !1 个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
ksel,s,loc,y,0
dk,all,all!约束墩底
*do,i,1,NT
ACEL,0,0,ac(i)!在Z方向加载 即在横桥向加载地震波
TIME,i*DT
solve
*enddo
!--------------------------------------------------
/post26
nsol,2,79,U,Z!定义主梁Z方向的位移
plvar,2,3,,,!绘制Z方向位移-时间曲线图
fini
/clear
/prep7
k,1
K,2,3
K,3,3,3.95
K,4,6.4,3.95
K,5,6.4,4.2
K,6,0,4.2
a,1,2,3,4,5,6
blc4,0,0.4,2.6,3.55
asba,1,2
wprot,0,0,90
wpoff,0,0,2.6
asbw,all
wpoff,0,0,0.4
asbw,all
wprot,0,90,0
wpoff,0,0,-0.4
asbw,all
wpoff,0,0,-3.55
asbw,all
wprot,0,-90,0
arsym,x,all
wpcsys
aglue,all
et,1,plane82
lesize,all,,,1
amesh,all
secwrite,mybeam!主梁截面定义 单箱单室 高4。2米，顶板厚0.25米，底板厚0。4米，肋宽0.4米
fini
/clear
/prep7
ET,1,beam189
MP,EX,1,3.0e10
MP,prxy,1,0.1667
MP,DAMP,1,0.05
MP,dens,1,2.549e3
sectype,1,beam,mesh
secread,mybeam,,,mesh
k,1,-32,25.4
k,33,0,25.4
k,73,40,25.4
k,105,72,25.4!定义生成主梁关键点 32米+40米+32米
k,1000,50,50
k,34,0,25
k,35,0,0
k,74,40,25
k,75,40,0 !定义生成桥墩关键点 墩高25米
n,36,-0.25,25
n,77,39.75,25!定义节点 以和桥顶节点 生成COMBIN40单元
l,1,33
l,33,73
l,73,105
lesize,all,0.5
latt,1,,1,,1000,,1
lmesh,all!主梁网格划分
SECTYPE,2,BEAM,RECT
SECDATA,3,6,6,12
l,34,35,25
l,74,75,25
latt,1,,1,,,,2
lmesh,all!桥墩网格划分
/eshape,1
eplot
!ksel,s,loc,x,-32
!dk,all,uy,uz,ux,rotx
!ksel,s,loc,x,72
!dk,all,uy,uz,rotx
nsel,all
d,78,uy,uz,rotx
d,511,uy,uz,rotx!约束梁端处 78号节点为主梁最左端点 511号节点为主梁最右端点
k1=6.3e5
k2=7e4
c=2e5
fd=4e4
ET,3,COMBIN40 !用COMBIN40单元模拟铅芯支座
KEYOPT,3,1,1
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,6,0
R,1,k1,c,1,,fd,k2
type,3
real,1
e,36,703
e,77,754!生成COMBIN40单元 即水平向的铅芯橡胶支座
et,4,combin14
keyopt,4,1,0
keyopt,4,2,0
keyopt,4,3,2
R,2,5000000,5000000,,
type,4
real,2
e,79,703
e,271,754!生成COMBIN14单元
!CPINTF,all
cp,1,ux,36,79
cp,2,ux,77,271 !顺桥向耦合
!定义和读入时程曲线
NT=1000!时程曲线有NT个点
DT=0.01!时间间隔
*dim,ac,,NT
!读入数据
/input,tianjin,txt !天津波共10 秒
/solu
NSUBST,1, , ,1 !1 个子步
OUTRES,ALL,1 !输出每个子步的结果
ANTYPE,TRANS !时程分析
ksel,s,loc,y,0
dk,all,all!约束墩底
*do,i,1,NT
ACEL,0,0,ac(i)!在Z方向加载 即在横桥向加载地震波
TIME,i*DT
solve
*enddo
!--------------------------------------------------
/post26
nsol,2,79,U,Z!定义主梁Z方向的位移
plvar,2,3,,,!绘制Z方向位移-时间曲线图
