[共建]＜ANSYS命令流实例库＞敬请您的加盟

chinatan · 发表于 2005-8-17 09:38:27

下面是自己完成的一个项目，我看还没有人贴过动力学方面的命令流，我就把自己写的这个瞬态动力学方面的命令流贴出来，希望对那些开始做动力学的朋友能有所帮助。命令流较长，要看朋友要放点耐心。建模的过程可能不是很关键，因为每次分析，模型都不一样，关键是看瞬态动力学的分析步骤以及相关的ansys设置。另外，循环语句、条件语句、复合材料中缠绕角的输入都是我们经常要用到的，还不清楚的朋友可以看看我的命令流相关的部分。

/prep7

!=======在本命令流中,均采用国际制单位
!=======长度单位：m
!=======密度单位：kg/m^3
!=======加速度的单位：m/s^2
!=======输出的应力单位为：pa

et,1,shell99
keyopt,1,8,1 !==存储并输出每一层的结果数据
et,2,shell93
et,3,surf154

!====内外蒙皮的模量、泊松比、密度
mp,ex,1,10.7e9
mp,ey,1,10.7e9
mp,ez,1,48e9
mp,prxy,1,0.36
mp,pryz,1,0.36
mp,prxz,1,0.36
mp,dens,1,2000

!====蜂窝
mp,ex,2,9e+007
mp,ey,2,9e7
mp,ez,2,9e7
mp,prxy,2,0.3
mp,pryz,2,0.3
mp,prxz,2,0.3
mp,gxy,2,60e6
mp,gyz,2,60e6
mp,gxz,2,60e6
mp,dens,2,48

!====frp
mp,ex,3,10e9
mp,prxy,3,0.3
mp,dens,3,1900

!====cfrp
mp,ex,4,65e9
mp,prxy,4,0.31
mp,dens,4,1780

!!!==============================下面是实常数的确定=========================!!!
!====内蒙皮：等角度变厚度；
!====外蒙皮：变角度等厚度；
!====桅杆的结构形式为：(外蒙皮)2h-3s-1h-夹层-1h-3s-1h(内蒙皮),总厚度：22.5mm
!====内蒙皮厚度的减少，我们只假设第二和第三层减少，第一层厚度不变

lay=0.0005    !==单层厚度(0.5mm)
wtk=lay*9       !==外蒙皮的总厚度
switch=0.014 !==夹层的厚度(15mm)
bthk=lay*9    !==底端内蒙皮的初始总厚度
tthk=bthk-0.0005 !==顶端内蒙皮厚度
dt=0.0005       !==内蒙皮厚度总减少量
dc=dt/2          !==每一层的总减少量，只假设内蒙皮第二和第三层减少，第一层厚度不变
dx=dc/29       !==每段内蒙皮厚度的减少量

lay1=lay    !==第1层厚度
lay2=lay*6  !==第2层初始厚度
lay3=lay*2  !==第3层厚度
lay4=0.014  !==第4层厚度(夹层厚度)
lay5=lay    !==第5层厚度
lay6=lay*6  !==第6层厚度
lay7=lay*2  !==第7层厚度

*dim,tk2,array,29  !定义一个一维数组,存储内蒙皮第二层每一段的厚度
*dim,tk3,array,29  !定义一个一维数组,存储内蒙皮第三层每一段的厚度
*do,i,1,29,1
tk2(i)=lay2-dx*(i-1)
tk3(i)=lay3-dx*(i-1)
*enddo

pi=3.1415926
jh=86*pi/180 !==所有环向缠绕角
j0=75*pi/180 !==内蒙皮的纵向缠绕角

!====下面是外蒙皮每段缠绕角的确定
j1=75*pi/180 !==第1段的缠绕角

*dim,h,,29  !==存储第1到第25段的高度
*dim,r,,29  !==每一段的半径
*dim,j,,29  !==每一段的缠绕角

j(1)=j1
*do,i,2,29,1
*if,i,gt,25,then
h(i)=h(25)+(i-25)*0.3
*else
h(i)=0.2+(i-2)*0.2
*endif
r(i)=0.25*(6.3-h(i))/6.3+0.125
j(i)=acos(0.375*cos(j1)/r(i))
*enddo

!h2=0.2
!r2=0.25*(6.3-h2)/6.3+0.125
!*set,j2,acos(0.375*cos(j1)/r2)

*do,i,1,29,1
r,i
rmodif,i,1,7

*if,i,lt,3,then !==第1段和第2段
rmodif,i,13,1,jh,lay1,1,j0,tk2(i),
rmodif,i,19,1,jh,tk3(i),4,0,switch,
rmodif,i,25,1,jh,lay5,1,j(i),lay6,
rmodif,i,31,1,jh,lay7

*elseif,i,eq,29 !==第29段
rmodif,i,13,1,jh,lay1,1,j0,tk2(i),
rmodif,i,19,1,jh,tk3(i),3,0,switch,
rmodif,i,25,1,jh,lay5,1,j(i),lay6,
rmodif,i,31,1,jh,lay7

*else !==第3段到第28段
rmodif,i,13,1,jh,lay1,1,j0,tk2(i),
rmodif,i,19,1,jh,tk3(i),2,0,switch,
rmodif,i,25,1,jh,lay5,1,j(i),lay6,
rmodif,i,31,1,jh,lay7

*endif
*enddo

r,30,0.006  !==支撑的厚度
r,31,0.04 !==法兰的厚度

!!!============================实常数输入完毕===============================!!!

!!!!===============================建模====================================!!!!

wpstyle,,,,,,,,1  !====显示工作平面
wpro,,,-7       !====将工作平面绕其y轴逆时针旋转7度

csys,4          !====将工作平面激活为当前的坐标系

!!!=================桅杆底端截面(包括法兰、金属护套)=================!!!
k,1,-0.3,-0.355,0
k,2,0,-0.375,0
k,3,0.3,-0.355
k,4,0.355,-0.3
k,5,0.375
k,6,0.355,0.3
k,7,0.3,0.355
k,8,0,0.375
k,9,-0.3,0.355
k,10,-0.355,0.3
k,11,-0.375
k,12,-0.355,-0.3 !====底端桅杆界面

bsplin,1,2,3
bsplin,4,5,6
bsplin,7,8,9
bsplin,10,11,12

l,3,4!l2tan,1,2
l,6,7!l2tan,2,3
l,9,10!l2tan,3,4
l,12,1!l2tan,4,1

k,13,-0.515,-0.515
k,14,0.765,-0.515
k,15,0.765,0.515
k,16,-0.515,0.515
k,17,0.625,-0.355
k,18,0.625,0.355
l,13,14
l,14,15
l,15,16
l,16,13
l,3,17
l,7,18
l,17,18

al,1,5,2,6,3,7,4,8
al,13,15,14,6,2,5
al,9,10,11,12  !====法兰

asel,s,area,,1,2,1
cm,sa,area
allsel,all
a***a,3,sa,,,keep

vext,1,,,0,0,6.347,0.333,0.333
vdele,1
adele,1
adele,3

wpro,,,7       !====将工作平面绕其y轴顺时针旋转7度，用于切割筒体

!========================灯具和天线的支撑==========================!

!====灯具1,2的支撑
!==计算得支撑所在桅杆的横界面的中心坐标为：(-0.658,0,5.4)

wpave,-1.220,0,5.400  !==把workplane的原点移到(-1.22,0,5.4)上
asel,s,area,,5,12,1
a***w,all
allsel,all

k,1000,0,0,0
circle,1000,0.161
k,1001,1.124
circle,1001,0.161 !====前后两个支撑

ldele,24
ldele,27
ldele,29
ldele,30

l,32,36
l,27,38
l,31,40
l,28,42

al,46,43,47,27,26,25,24
al,45,29,28,31,30,44,37

k,1002,0,0,0
k,1003,0,0,-0.050
l,1002,1003

adrag,25,26,,,,,56
adrag,28,31,,,,,56

wpoff,,,-0.3

asel,s,area,,1,3,2
asel,a,area,,13,18,1
a***w,all
allsel,all

l,45,49
l,43,54
l,46,53
l,48,50
al,27,83,32,61
al,24,88,33,58
al,30,66,36,84
al,29,87,35,63

!====天线2,3的支撑
!==计算得支撑所在桅杆的横界面的中心坐标为：(-0.7,0,5.7)
wpoff,,,0.600

asel,s,area,,19,26,1
a***w,all

csys,0
wpave,-0.700,0,5.700 !==把workplane的原点移到(-0.7,0,5.7)上
csys,4

wpoff,,0.800

k,2000,0,0,0
circle,2000,0.149
k,2001,0,-1.600,0
circle,2001,0.149
ldele,52,53,1,1
ldele,50,51,1,1
allsel,all !==画两个半圆

l,59,69
l,64,71
l,63,67
l,60,65

k,2002,0,0,0
k,2003,0,0,-0.050
l,2002,2003

adrag,48,49,,,,,111
adrag,54,55,,,,,111

a,29,76,69,59
a,34,64,71,74
a,33,73,67,63
a,30,60,65,68

al,48,53,95,97,99,52,49
al,50,55,54,51,102,41,91

!====顶部天线1、避雷针支座
!==计算得支撑所在桅杆的横界面的中心坐标为：(-0.737,0,6)
asel,s,area,,15,16,1
asel,a,area,,17,18,1
asel,a,area,,41,42,1
asel,a,area,,43,44,1
/replot

csys,0
wpave,-0.737,0,6.000 !==把workplane的原点移到(-0.737,0,6)上
csys,4
a***w,all

csys,0
wpave,-0.773,0,6.280
csys,4
!asel,s,area,,55,62,1
a***w,55
a***w,56
a***w,57
a***w,58
a***w,59
a***w,60
a***w,61
a***w,62

adele,15,17,2,1
adele,41,43,2,1
adele,55,57,2,1
adele,59,61,2,1

wpoff,0.680

k,3000,0,0,0
circle,3000,0.125
ldele,17,18,1,1
k,3001,0,0,-0.050
l,3000,3001
adrag,16,19,,,,,17

l,89,20
l,86,22

al,16,19,39,98,127,150,38
a,22,24,78,86
a,20,89,81,21
allsel,all

!====灯具3支座
!==计算得支撑所在桅杆的横界面的中心坐标为：(-0.466,0,3.8)

csys,0
wpave,-0.466,0,3.800  !==把workplane的原点移到(-0.466,0,3.8)上
asel,s,area,,11,12,1
asel,a,area,,27,32,1
a***w,all
wpoff,,,-0.200
asel,s,area,,57,59,2
asel,a,area,,61,62,1
asel,a,area,,72,75,1
a***w,all
csys,4

wpoff,,0.700
wpoff,,,0.200
k,4000,0,0,0
circle,4000,0.2258
ldele,137,158,21,1
k,4001,0,0,-0.050
l,4000,4001

adrag,131,133,,,,,137

l,94,107
l,97,109
l,102,110
l,105,112

al,131,194,167,162,168,195,133
al,160,196,189,194
al,164,197,192,195
allsel,all

!====灯具4支座
!==计算得支撑所在桅杆的横界面的中心坐标为：(-0.221,0,1.8)
csys,0
wpave,-0.221,0,1.800  !==把workplane的原点移到(-0.221,0,1.8)上
asel,s,area,,11,12,1
asel,a,area,,27,32,1
a***w,all
wpoff,,,-0.200

asel,s,area,,73,75,1
asel,a,area,,92,96,1
a***w,all

csys,4

wpoff,,0.700
wpoff,,,0.200
k,5000,0,0,0
circle,5000,0.3065

ldele,201,202,1,1
k,5001,0,0,-0.050
l,5000,5001
adrag,198,199,,,,,201

l,129,116
l,132,124
l,131,119
l,134,127

al,198,238,211,206,212,240,199
al,204,239,233,238
al,208,240,236,241

allsel,all
aplot
numcmp,all

csys,0
wpave,0,0,0 !==把workplane的原点移到全局坐标的原点上
cyl4,0,0,1.000
!==以上面的圆面进行划分
agen,8,110,,,,,0.200
agen,2,110,,,,,2.000
agen,8,118,,,,,0.200
agen,2,125,,,,,0.600
agen,5,126,,,,,0.200

asel,s,area,,1,109,1
cm,a1,area
allsel,all
asel,s,area,,111,130,1
cm,a2,area
allsel,all
a***a,a1,a2
adele,110
numcmp,all
wpstyle,,,,,,,,0
ldele,17,,,1
ldele,56,,,1
ldele,107,,,1
ldele,133,,,1
ldele,187,,,1
wpstyle,,,,,,,,0
!!!!=============================完成建模==================================!!!!

!!!==============划分网格=============!!!!
!==先对筒体进行划分
lsel,s,line,,5,8,1
lsel,a,line,,44,47,1
lsel,a,line,,75,78,1
lsel,a,line,,87,91,4
lsel,a,line,,94,95,1
lsel,a,line,,98
lsel,a,line,,126,134,4
lsel,a,line,,135,137,2
lsel,a,line,,146,150,4
lsel,a,line,,162,165,1
lsel,a,line,,168,171,1
lsel,a,line,,196,199,1
lsel,a,line,,202,205,1
lsel,a,line,,234,237,1
lsel,a,line,,262,265,1
lsel,a,line,,278,281,1
lsel,a,line,,288,291,1
lsel,a,line,,304,307,1
lsel,a,line,,314,317,1
lsel,a,line,,367,369,1
lsel,a,line,,382
lsel,a,line,,386,390,1
lsel,a,line,,439,461,1
lsel,a,line,,542,565,1

lesize,all,,,1
allsel,all

type,1
mshape,0,2d
mshkey,1

real,1
asel,s,area,,86,93,1
amesh,all

real,2
asel,s,area,,134,138,4
asel,a,area,,143,144,1
asel,a,area,,156,159,3
asel,a,area,,187,191,4
amesh,all

real,3
asel,s,area,,139,145,3
asel,a,area,,157,160,3
asel,a,area,,188,192,4
asel,a,area,,201
amesh,all

real,4
asel,s,area,,135,140,5
asel,a,area,,146,152,6
asel,a,area,,161,189,28
asel,a,area,,193,202,9
amesh,all

real,5
asel,s,area,,141,147,6
asel,a,area,,150,153,3
asel,a,area,,158,162,4
asel,a,area,,194,203,9
amesh,all

real,6
asel,s,area,,136,148,12
asel,a,area,,151,154,3
asel,a,area,,163,165,2
asel,a,area,,195,204,9
amesh,all

real,7
asel,s,area,,137,149,12
asel,a,area,,155,164,9
asel,a,area,,166,190,24
asel,a,area,,196,205,9
amesh,all

real,8
asel,s,area,,96,100,4
asel,a,area,,128,133,1
amesh,all

real,9
asel,s,area,,78,85,1
amesh,all

real,10
asel,s,area,,120,127,1
amesh,all

real,11
asel,s,area,,214,263,7
amesh,all

real,12
asel,s,area,,215,264,7
amesh,all

real,13
asel,s,area,,216,265,7
amesh,all

real,14
asel,s,area,,217,266,7
amesh,all

real,15
asel,s,area,,218,267,7
amesh,all

real,16
asel,s,area,,219,268,7
amesh,all

real,17
asel,s,area,,220,269,7
amesh,all

real,18
asel,s,area,,95,98,3
asel,a,area,,114,119,1
amesh,all

real,19
asel,s,area,,68,75,1
amesh,all

real,20
asel,s,area,,106,113,1
amesh,all

real,21
asel,s,area,,170,171,1
asel,a,area,,175,183,4
asel,a,area,,197,206,9
asel,a,area,,210
amesh,all

real,22
asel,s,area,,169,172,3
asel,a,area,,176,184,4
asel,a,area,,198,207,9
asel,a,area,,211
amesh,all

real,23
asel,s,area,,167,173,6
asel,a,area,,177,185,4
asel,a,area,,199,208,9
asel,a,area,,212
amesh,all

real,24
asel,s,area,,168,174,6
asel,a,area,,178,186,4
asel,a,area,,200,209,9
asel,a,area,,213
amesh,all

real,25
asel,s,area,,94,97,3
asel,a,area,,99
asel,a,area,,101,105,1
amesh,all

real,26
asel,s,area,,25,32,1
amesh,all

real,27
asel,s,area,,37,44,1
amesh,all

real,28
asel,s,area,,55,62,1
amesh,all

real,29
asel,s,area,,14,16,2
asel,a,area,,34,36,2
asel,a,area,,48,52,2
asel,a,area,,63
amesh,all

allsel,all

!!==再对支撑进行划分
asel,s,area,,1,3,2
asel,a,area,,7,13,1
asel,a,area,,17,24,1
asel,a,area,,15,35,20
asel,a,area,,47,51,2
asel,a,area,,54
asel,a,area,,64,66,1
asel,a,area,,76,77,1

type,2
mat,3
real,30
mshape,0,2d
mshkey,1
amesh,all

allsel,all

asel,s,area,,5,6,1
asel,a,area,,33
asel,a,area,,45,46,1
asel,a,area,,53,67,14

mshkey,0
amesh,all

asel,s,area,,3,11,8
asel,a,area,,5,8,1
asel,a,area,,17,24,1
asel,a,area,,35
asel,a,area,,46,47,1
asel,a,area,,49,51,2
asel,a,area,,54
asel,a,area,,65,66,1
asel,a,area,,76
areverse,all,0  !!==让这些面反向

!!==再对进行划分法兰进行划分
allsel,all
csys,wp
wpro,,,-7
asel,s,area,,2,4,2
pcirc,0.015,0,0,360 !==螺栓的直径为30mm
agen,2,270,,,-0.425,-0.425
agen,7,271,,,0,0.850/6
agen,9,271,,,1.100/8,0
agen,7,285,,,0,0.850/6
agen,8,277,,,1.100/8

adele,270,,,1
aovlap,all

adele,271,298,1,1

type,2
mat,4
real,31
smrt,6
asel,s,area,,2,270,268
amesh,all
allsel,all

!!!===========================瞬态动力学分析=============================!!!
/solu

/NERR,5,100000000 !===只输出5个warnings and errors，如果总的警告或错误数超过100000000则自动退出ansys

g=9.8
!==x负方向风载作用面
asel,s,area,,26
asel,a,area,,34
asel,a,area,,38
asel,a,area,,57
asel,a,area,,69
asel,a,area,,79
asel,a,area,,87
asel,a,area,,97
asel,a,area,,98
asel,a,area,,100
asel,a,area,,107
asel,a,area,,121
asel,a,area,,138,141,1
asel,a,area,,165,170,1
asel,a,area,,221,227,1
cm,xplane,area
allsel,all
wpstyle,,,,,,,,0

antype,trans

outpr,all,all
outres,all,all!===输出控制

!===第一步载荷计算（create initial condition）

lsel,s,line,,187
lsel,a,line,,566,676,1
dl,all,,all !===施加约束条件
time,0.01
autots,1    !===打开自动时间步长

kbc,1 !===stepped方式加载
acel,,,g
solve

!===第二步载荷计算
time,30
deltim,0.8
kbc,1
sfa,xplane,1,pres,-1962 !===风载
kbc,0
acel,36*g,,61*g!===纵向和垂直方向的惯性力
solve

!===第三步载荷计算
time,75
deltim,0.8
kbc,1
sfa,xplane,1,pres,-1962 !===风载
kbc,0
acel,0,,g!===纵向和垂直方向的惯性力
solve

/post26
/ANG, 1 ,-60,YS,1

总的应力图（时间为：30）：

MYP123 · 发表于 2005-9-28 15:35:33

一个比较典型的好例子，与大家共享
!进行初始地应力的计算
esel,s,type,,3    !选择支护壳单元为当前有效单元，然后将其杀死
ekill,all
esel,all
esel,s,live       !选择活的单元，即所有土体单元
nsle,s             !选择活单元上的节点
nsel,invert       !反向选择，即选择了死单元上的节点
d,all,all          !将死单元上的节点约束所有位移，使其不参与矩阵运算
nsel,all
esel,all
/PBC,ALL,,1
gplot
solve             !进行初始地应力的计算

!开挖过程的基本思路是：将挖去土体并将它杀死的同时，激活支护单元。可以将每天开挖后的计算状态保存为一个
载荷文件，然后统一计；也可以用循环语句来实现。以下施工模拟过程是基于简化假设；每天挖5m，分10天完成。因此，
整个计算过程应该包括1个初始地应力计算载荷步和10个开挖过程计算载荷步。

*do,ii,1,10,1
!以下步骤表示，先选择每天挖去的岩体单元为有效单元，然后将其杀死
esel,s,mat,,3
nsle,s
nsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)
esln,r,1
ekill,all
!以下步骤表示，先选择每天挖去的岩体单元对应的支护壳单元为当前有效单元，然后将其激活，同时将其节点上的约束
删除
esel,s,type,,3
nsle,s
nsel,r,loc,z,0.1-(ii-1)*5,-(5.1+(ii-1)*5)
esln,r,1
ealive,all
nsle,s
ddele,all,all
!选择活单元，此时应该包括两部分：刚被激活的壳单元和未挖去的岩体单元
esel,all
esel,s,live
nsle,s
nsel,invert !反向选择，将死单元上的节点约束所有自由度
d,all,all,
nsel,all
esel,all

solve
*enddo
fini

songfeng2002 · 发表于 2005-10-3 18:21:11

我看了上面大家的命令流库，感觉很多的命令流都是摘抄别人的，甚至是完全抄袭，我想给大家说的是，大家不要这样子，最好的是各位同仁能够把自己的项目或者手头上的工程所用的命令流贴出来供大家参考！为了表示诚意，我首先把我以前的一个项目贴出来给大家参考！
需要说明的是，我的这个分析结果和思路未必正确，因为这个只是我几个方案中的一个分析方案，最终方案并没有采用这个分析方案，但可能对于初级水平或者刚开始接触命令流的朋友具有一定的帮助！

显示全部楼层 · 发表于 2005-10-7 17:06:27

提示: 作者被禁止或删除内容自动屏蔽

Bob · 发表于 2005-10-8 13:57:59

碳纤维布加固钢筋混凝土
*SET,H,80
*SET,B,500/2
*SET,L,1800
*SET,LF,100
*SET,LS,50
*SET,A,20
*SET,PI,ACOS(-1)
*SET,SR,PI*(8/2)**2
*SET,RO,PI*(6/2)**2/200/H
*SET,CB,300
*SET,CH,0.165
*SET,CL,L-2*LS

*SET,F1,0.002
*SET,F2,30/LF

/PREP7
ET,1,SOLID65
ET,2,LINK8
ET,3,SHELL41
KEYOPT,3,1,2

ET,4,SOLID45
MP,EX,1,18000
MP,PRXY,1,0.2
TB,KINH,1,1,7
TBPT,,0.0001,1.8
TBPT,,0.0004,6.66
TBPT,,0.0008,11.84
TBPT,,0.0012,15.54
TBPT,,0.0016,17.76
TBPT,,0.0020,18.5
TBPT,,0.0033,18.5

!TB,CONCR,1,1,9
!TBDATA,,0.3,0.5,1.75,-1

TB,CONCR,1
TBDATA,1,0.5,0.95,1.75,-1

MP,EX,2,2.1E5
MP,PRXY,2,0.3
TB,BKIN,2,1,2,1
TBDATA,,235,0

MP,EX,3,2.35E5
MP,PRXY,3,0

MP,EX,4,2.1E5
MP,PRXY,4,0.3

R,1,2,RO
RMORE,

R,2,SR,
R,3,SR/2,
R,4,CH,
RMORE,

block,0,-b,0,h,-lf/2,-(l/2-ls),

lgen,2,8,,,50,
lgen,2,13,,,100,

adrag,13,,,,,,9
adrag,14,,,,,,9

v***a,1,7
v***a,3,8

lsel,s,loc,z,-lf/2
lsel,r,loc,y,0

lgen,2,all,,,,a,

lsel,s,loc,z,-lf/2
lsel,r,loc,y,a

adrag,all,,,,,,9

allsel,all

v***a,2,15
v***a,4,14
v***a,1,13

lsel,s,loc,y,0
lsel,a,loc,y,a
lsel,a,loc,y,h

lesize,all,50

lsel,s,loc,z,-lf/2
lsel,a,loc,z,-(l/2-ls)
lsel,u,loc,y,0
lsel,u,loc,y,a
lsel,u,loc,y,h

lesize,all,20

allsel,all

type,1
mat,1
real,1
esys,0
mshape,0,3d
mshkey,1

vmesh,all

asel,s,loc,z,-lf/2
type,1
extopt,esize,lf/50
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,lf,,,,

asel,s,loc,z,lf/2
type,1
extopt,esize,(l/2-lf/2-ls)/50
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,l/2-lf/2-ls,,,,

allsel,all
asel,s,loc,z,l/2-ls
type,1
extopt,esize,1
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,25

asel,s,loc,z,l/2-ls+25
type,1
extopt,esize,2
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,ls

asel,s,loc,z,-(l/2-ls)
type,1
extopt,esize,1
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,-25

asel,s,loc,z,-(l/2-ls+25)
type,1
extopt,esize,2
extopt,aclear,1
extopt,attr,0,0,0
MAT,1
REAL,1
ESYS,0

VEXT,all, , ,0,0,-ls

allsel,all

!
划分受力钢筋
lsel,s,loc,x,-200
lsel,a,loc,x,-100
lsel,r,loc,y,a

type,2
MAT,2
REAL,2
ESYS,0

lmesh,all

lsel,s,loc,x,0
lsel,r,loc,y,a

type,2
MAT,2
REAL,3
ESYS,0

lmesh,all

allsel,all

!炭纤维
ksel,s,loc,z,850
ksel,r,loc,y,0
ksel,r,loc,x,0

kgen,2,all,,,-cb/2,,,,0

lstr,41,97
allsel,all
lstr,41,4

adrag,221,,,,,,222
asel,s,,,168

lsla,s

lesize,all,50,,,,,,,1
type,3
mat,3
real,4
esys,0
mshape,0,2d
amesh,all
allsel,all
nummrg,node,,,,low

numcmp,node

!支座
asel,s,loc,y,0
asel,r,loc,z,l/2-25,l/2+25
type,4

extopt,esize,1
extopt,aclear,1
extopt,attr,0,0,0
MAT,4
REAL,4
ESYS,0

VEXT,all, , ,0,-20

asel,s,loc,y,0
asel,r,loc,z,-(l/2-25),-(l/2+25)
type,4

extopt,esize,1
extopt,aclear,1
extopt,attr,0,0,0
MAT,4
REAL,4
ESYS,0

VEXT,all, , ,0,-20

allsel,all
finish

/solu
nsel,s,loc,x,0
dsym,symm,x

nsel,s,loc,z,l/2
nsel,r,loc,y,-20
d,all,,,,,,uy,

nsel,s,loc,z,-l/2
nsel,r,loc,y,-20
d,all,,,,,,uy,uz,

antype,0
!nlgeom,1
nropt,full
eqslv,spar,,0
outres,all,all

time,1
autots,1
!nsubst,10,20,5,1
nsubst,10
cnvtol,f,,0.05,2
neqit,30
pred,on,,on

esel,mat,3
ekill,all

allsel,all

asel,s,loc,y,h
sfa,all,1,pres,f1
solve

time,2
autots,0
nsubst,25,,,1
pred,-1
arclen,1,10,0

asel,s,loc,y,h
asel,r,loc,z,-lf/2,lf/2
sfa,all,1,pres,f2*0.75

esel,mat,3
ealive,all
allsel,all
solve

*cfopen,uy,dat
*dim,uy,array,31
*do i,1,31
SET,,, ,,, ,i
*get,uy(i),node,540,u,y
*vwrite,uy(i)
*enddo

/post26
allsel,all
nsel,s,loc,x,0
nsel,r,loc,y,0
nsel,r,loc,z,0
*get,n1,node,0,num,max
nsol,2,n1,u,y
prod,3,2,,,,,,-1,1,1
xvar,3
plvar,1
!
/post1
etable,zxyl,ls,1
plls,zxyl,zxyl,1

fpemail · 发表于 2005-10-12 23:51:16

大家好！
三个多月没有上网了，因为个人家里住址变动，所以上直以来没有条件上网。现在本人安装了1M宽带，可以每天都来看看。
三个多月来，我没想到自已的帖子被大家如此关注！并且还被置顶！
感谢大家！谢谢！现在我被提为公司的强度组组长，日后还会有很多心得与大家一起分享与探讨，希望还能够得到大家的支持！同时我会尽我所能把这个帖子做得更加丰富！

显示全部楼层 · 发表于 2005-10-16 10:25:33

提示: 作者被禁止或删除内容自动屏蔽

fpemail · 发表于 2005-10-20 02:15:43

lsh0016 wrote:
好贴啊。顺便问一下楼主，命令流文件跟gui交互式处理有没有相互转化的命令

你的gui交互操作可以转化成命令流，不过有很多选择之类的命令需要转化！

wl2081 · 发表于 2005-10-21 20:36:42

用坐标系操作来绘制圆：

命令流:

CSWPLA,11,1,0.5,1
/PREP7
K,1,-5
K,2,5
L,1,2
CSYS,0
LSYMM,Y,1,,,,0,0
lglue,all
al,1,3

详细解释:

查阅帮助得知:

1.CSWPLA：
/*
CSWPLA, KCN, KCS, PAR1, PAR2

Defines a local coordinate system at the origin of the working plane.

KCN
Arbitrary reference number assigned to this coordinate system. Must be greater than 10. A coordinate system previously defined with this number will be redefined.

KCS
Coordinate system type:

0 or CART --
Cartesian

1 or CYLIN --
Cylindrical (circular or elliptical)

2 or SPHE --
Spherical (or spheroidal)

3 or TORO --
Toroidal

PAR1
Used for elliptical, spheroidal, or toroidal systems. If KCS = 1 or 2, PAR1 is the ratio of the ellipse Y-axis radius to X-axis radius (defaults to 1.0 (circle)). If KCS = 3, PAR1 is the major radius of the torus.

PAR2
Used for spheroidal systems. If KCS = 2, PAR2 = ratio of ellipse Z-axis radius to X-axis radius (defaults to 1.0 (circle)).

*/ 可以知道第一句是定义了一个圆柱坐标系
第三四句就创建了椭圆的两个端点.
第五句创建连接两个点构成一条线（因为是圆柱坐标系，所以绘制的线不是直线了：）
csys，0切换到整体坐标系.
最后"LSYMM,Y,1,,,,0,0"一句进行对称绘制.
然后lglue,all把两条线连接起来，最后一句绘制椭圆面.
这个里面比较难以理解的是用l,1,2命令连接两点得到的不是直线而是曲线.

yymaverick · 发表于 2005-10-25 15:58:11

不错。不过这样一个一个贴显得有点乱，是否可以只贴出问题的简单描述和结果图，命令流贴成附件。管理员可按专题分贴，这样大家也好找。

9921224 · 发表于 2005-10-27 09:22:37

刚看到这么好的贴,也跟一个
!曲轴结构分析（刘鸿文编《材料力学》第362页例9.6）
fini
/CLEAR,START
/prep7
/GRAPHICS,FULL
/COLOR,PBAK,OFF
*AFUN,DEG
!
*SET,d1,50
*SET,d,60
*SET,h,102
*SET,b,22
*SET,pi,3.1415926
*SET,beta,0.288
! ------------------------------------------------
! 定义连杆轴颈截面的几何参数
*SET,a1,0.25*pi*d1*d1
*SET,i1xx,pi*d1*d1*d1*d1/32
*SET,i1y,0.5*i1xx
*SET,i1z,0.5*i1xx
*SET,tkz1,d1
*SET,tky1,d1
!*
R,1,a1,i1z,i1y,tkz1,tky1, ,
RMORE, ,i1xx, , , , ,
! ------------------------------------------------
! 定义曲柄截面的几何参数
*SET,a2,b*h
*SET,i2xx,beta*b*b*b*h
*SET,i2y,b*b*b*h/12
*SET,i2z,b*h*h*h/12
*SET,tky2,h
*SET,tkz2,b
!*
R,2,a2,i2z,i2y,tkz2,tky2, ,
RMORE, ,i2xx, , , , ,
! -----------------------------------------------
! 定义主轴颈截面的几何参数
*SET,a3,0.25*pi*d*d
*SET,i3xx,pi*d*d*d*d/32
*SET,i3y,0.5*i3xx
*SET,i3z,0.5*i3xx
*SET,tky3,d
*SET,tkz3,d
!*
R,3,a3,i3z,i3y,tkz3,tky3, ,
RMORE, ,i3xx, , , , ,
! -----------------------------------------------
/PREP7
!*
ET,1,BEAM4
!*
MAT,1
MP,EX,1,2e5
MP,PRXY,1,.3
MP,DENS,1,7.8e-9
!
N,1,0,0,0,,,,
N,4,33,0,0,,,,
N,14,33,60,0,,,,
N,22,97,60,0,,,,
N,32,97,0,0,,,,
N,35,130,0,0,,,,
!
FILL,1,4,2, , ,1,1,1,
FILL,4,14,9, , ,1,1,1,
FILL,14,22,7, , ,1,1,1,
FILL,22,32,9, , ,1,1,1,
FILL,32,35,2, , ,1,1,1,
! --------------------------------------------------
! 定义主轴颈的单元
TYPE,1
MAT,1
REAL,3
ESYS,0
SECNUM,
TSHAP,LINE
E,1,2,5
E,2,3,6
E,3,4,7
E,32,33,22
E,33,34,23
E,34,35,24
! --------------------------------------------------
! 定义曲柄的单元
TYPE,1
MAT,1
REAL,2
ESYS,0
SECNUM,
TSHAP,LINE
!*
E,4,5,15
E,5,6,15
E,6,7,15
E,7,8,15
E,8,9,15
E,9,10,15
E,10,11,15
E,11,12,15
E,12,13,15
E,13,14,15
E,22,23,21
E,23,24,21
E,24,25,21
E,25,26,21
E,26,27,21
E,27,28,21
E,28,29,21
E,29,30,21
E,30,31,21
E,31,32,21
! -----------------------------------------------
! 定义主轴颈的单元
TYPE,1
MAT,1
REAL,1
ESYS,0
SECNUM,
TSHAP,LINE
!*
E,14,15,23
E,15,16,24
E,16,17,25
E,17,18,26
E,18,19,27
E,19,20,28
E,20,21,29
E,21,22,30
!*
! -------------------------------------------
EPLOT
!*
FINISH
/SOL
D,1, ,0, , , ,UX,UY,UZ,ROTX, ,
D,35, ,0, , , ,UY,UZ, , , ,
!*
F,18,FY,-32000
F,18,FZ,17000
F,35,MX,-1020000
F,14,FY,3000
F,22,FY,3000
F,4,FY,-7000
F,32,FY,-7000
!
/STATUS,SOLU
SOLVE
FINISH
!*
/POST1
ETABLE,momentxi,SMISC,4
ETABLE,momentxj,SMISC,10
ETABLE,momentyi,SMISC,5
ETABLE,momentyj,SMISC,11
ETABLE,momentzi,SMISC,6
ETABLE,momentzj,SMISC,12
ETABLE,forcexi,SMISC,1
ETABLE,forcexj,SMISC,7
ETABLE,forcezi,SMISC,3
ETABLE,forcezj,SMISC,9
ETABLE,forceyi,SMISC,2
ETABLE,forceyj,SMISC,8
!
ETABLE,stressxi,LS,1
ETABLE,stressxj,LS,6
ETABLE,stresypi,LS,2
ETABLE,stresypj,LS,7
ETABLE,stresyni,LS,3
ETABLE,stresynj,LS,8
ETABLE,streszpi,LS,4
ETABLE,streszpj,LS,9
ETABLE,streszni,LS,5
ETABLE,stresznj,LS,10
ETABLE,strsmaxi,NMISC,1
ETABLE,strsmaxj,NMISC,3
ETABLE,strsmini,NMISC,2
ETABLE,strsminj,NMISC,4
!
! 计算主轴颈截面2的等效应力
*GET,m2xi,ELEM,4,ETAB,MOMENTXI
*GET,m2xj,ELEM,4,ETAB,MOMENTXJ
*GET,m2yi,ELEM,4,ETAB,MOMENTYI
*GET,m2yj,ELEM,4,ETAB,MOMENTYJ
*GET,m2zi,ELEM,4,ETAB,MOMENTZI
*GET,m2zj,ELEM,4,ETAB,MOMENTZJ
*SET,M2,SQRT(M2Yj*M2Yj+M2Zj*M2Zj)
*SET,W3,2*I3Y/D
*SET,SIGMA2,SQRT(M2*M2+0.75*M2Xj*M2Xj)/W3
!* -------- -------------------------------------
! 计算连杆轴颈1的等效截面
*GET,m1XI,ELEM,30,ETAB,MOMENTXI
*GET,m1XJ,ELEM,30,ETAB,MOMENTXJ
*GET,m1YI,ELEM,30,ETAB,MOMENTYI
*GET,m1YJ,ELEM,30,ETAB,MOMENTYJ
*GET,m1ZI,ELEM,30,ETAB,MOMENTZI
*GET,m1ZJ,ELEM,30,ETAB,MOMENTZJ
*SET,M1,SQRT(M1YJ*M1YJ+M1ZJ*M1ZJ)
*SET,W1,2*I1Y/D1
*SET,SIGMA1,SQRT(M1*M1+0.75*M1ZJ*M1XJ)/W1
!*  ------------------------------------------------
! 计算曲柄截面3的应力
*GET,S3MAXI,ELEM,21,ETAB,STRSMAXI
*GET,S3MAXJ,ELEM,21,ETAB,STRSMAXJ
*GET,S3MINI,ELEM,21,ETAB,STRSMINI
*GET,S3MINJ,ELEM,21,ETAB,STRSMINJ
*SET,S3MAXI,ABS(S3MAXI)
*SET,S3MAXJ,ABS(S3MAXJ)
*SET,S3MINI,ABS(S3MINI)
*SET,S3MINJ,ABS(S3MINJ)
*SET,sigma3,max(s3maxj,s3minj)

fpemail · 发表于 2005-10-27 22:25:36

yymaverick wrote:
不错。不过这样一个一个贴显得有点乱，是否可以只贴出问题的简单描述和结果图，命令流贴成附件。管理员可按专题分贴，这样大家也好找。

这个主意好!不过我只是一名普通的会员，在论坛里没有什么权限!还希望斑竹多多给把关!谢谢!

fpemail · 发表于 2005-10-30 22:30:56

为了能够更好把这个主题帖做好!为大家切实的服务!请大家真对帖子的沾帖方法提出意见!因为现在的帖了格局不是太合理，看起来很乱!谢谢东FPEMAIL

happyhoney555 · 发表于 2005-11-9 22:36:07

发一个简单的混凝土初应力问题
/prep7
et,1,solid45
mp,ex,1,3e4
mp,alpx,1,1e-5
mp,nuxy,1,0.3
block,0,10,0,10,0,10
esize,1
vmesh,1
FINISH
/SOL
bf,all,temp,-5
asel,s,loc,y,0
da,all,all
solve
FINISH
/POST1
/view, 1 ,1,1,1
/ANG,1
/REP,FAST
pldisp,1
plnsol,s,x,0,1
！将x方向的应力写入x.txt文件
*get,n_num,node,,count
*cfopen,x,txt
*do,j,1,n_num
*get,sx,node,j,s,x
*vwrite,sx
(f25.15)
*enddo
*cfclos
fini

savagebob · 发表于 2005-11-19 19:16:43

转贴！！！
幕墙企业玻璃简化计算
/encrypt,yuhui4,Magic Glass,mac,d:\
/nopr
finish
/CLEAR,NOSTART
multipro,'start',10
*cset,1,3,H,'The height of glass(unit:m)',2
*cset,4,6,W,'The width of glass(unit:m)',2
*cset,7,9,a1,'Angle left side to down(unit:C)',90
*cset,10,12,a2,'Angle right side to down(unit:C)',90
*cset,13,15,hole,'the number of hole',4
*cset,16,18,type,'Layered Glass Type(!SeePeak!)',1
*cset,19,21,thi,'Thickness of inner Glass(unit:m)',0.01
*cset,22,24,tho,'Thickness of outer Glass(unit:m)',0.008
*cset,25,27,Wk,'Designed wind pressure(unit:N/m2)',1000
*cset,61,62,'Please Fill in related blanks to',' build FEA model And apply loads'
*cset,28,30,Seismic,'Designed seismic accelerate(m/s2)',0.08
*cset,63,64,'FOR LAYERED GLASS TYPE:1--[dan','pian],2--[jiajiao],3--[ganghua]'
multipro,'end'
*if,hole,eq,0,then
multipro,'start',8
*cset,1,3,doftop,'constrain type on the top side',2
*cset,4,6,dofdown,'constrain type on the down side',1
*cset,7,9,dofleft,'constrain type on the left side',2
*cset,10,12,dofright,'constrain type on the right side',1
*cset,13,15,Wdiv,'divide number on the width side',40
*cset,16,18,Hdiv,'divide number on the heigth side',40
*cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0
*cset,22,24,switch,'Choose solution type',1
*cset,61,62,'Please Fill in related blanks to',' Get FEA solution'
*cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
multipro,'end'
*elseif,hole,eq,2,then
multipro,'start',9
*cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175
*cset,4,6,kbw,'Distance to width side(unit:m)',0.106
*cset,7,9,kbh,'Distance to height side(unit:m)',0.106
*cset,10,12,dofdown,'constrain type on the down side',1
*cset,13,15,Wdiv,'divide number on the width side',40
*cset,16,18,Hdiv,'divide number on the heigth side',40
*cset,19,21,holediv,'divide number on the hole side',32
*cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0
*cset,25,27,switch,'Choose solution type',1
*cset,61,62,'Please Fill in related blanks to',' Get FEA solution'
*cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
multipro,'end'
*elseif,hole,eq,4,then
multipro,'start',8
*cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175
*cset,4,6,kbw,'Distance to width side(unit:m)',0.106
*cset,7,9,kbh,'Distance to height side(unit:m)',0.106
*cset,10,12,Wdiv,'divide number on the width side',40
*cset,13,15,Hdiv,'divide number on the heigth side',40
*cset,16,18,holediv,'divide number on the hole side',32
*cset,19,21,gravity,'The gravity Acceleration(unit:m/s2)',0
*cset,22,24,switch,'Choose solution type',1
*cset,61,62,'Please Fill in related blanks to',' Get FEA solution'
*cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
multipro,'end'
*elseif,hole,eq,6,then
multipro,'start',9
*cset,1,3,kr,'Radium of glass hole(unit:m)',0.0175
*cset,4,6,kbw,'Distance to width side(unit:m)',0.106
*cset,7,9,kbh,'Distance to height side(unit:m)',0.106
*cset,10,12,Dbottom,'hole central distance to bottom',1
*cset,13,15,Wdiv,'divide number on the width side',40
*cset,16,18,Hdiv,'divide number on the heigth side',40
*cset,19,21,holediv,'divide number on the hole side',32
*cset,22,24,gravity,'The gravity Acceleration(unit:m/s2)',0
*cset,25,27,switch,'Choose solution type',1
*cset,61,62,'Please Fill in related blanks to',' Get FEA solution'
*cset,63,64,'SWITCH:0--[By Youself]1--[Small',',Deformed] 2--[LargeDeformed]'
multipro,'end'
*endif
*if,hole,eq,0,then!!!!!!!!!!平板建模
/prep7
csys,0
k,1,1,1
k,2,1+w,1
*afun,deg
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
a,1,2,4,3
et,1,shell63
*if,type,eq,1,then
thforcal=1.1*thi
*elseif,type,eq,2,then
thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
thforcal=1.1*1.20*thi
*endif
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
lesize,all,,,hdiv
lsel,all
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
lesize,all,,,hdiv
lsel,all
amesh,1
*if,doftop,eq,0,then
*elseif,doftop,eq,1,then
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
dl,all,,uz
dl,all,,uX
allsel
*elseif,doftop,eq,2,then
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
dl,all,,uz
dl,all,,uy
allsel
*endif
*if,dofdown,eq,0,then
*elseif,dofdown,eq,1,then
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
dl,all,,uz
dl,all,,uX
allsel
*elseif,dofdown,eq,2,then
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
dl,all,,uz
dl,all,,uy
allsel
*endif
*if,dofleft,eq,0,then
*elseif,dofleft,eq,1,then
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
dl,all,,uz
dl,all,,uY
allsel
*elseif,dofleft,eq,2,then
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
dl,all,,uz
dl,all,,ux
allsel
*endif
*if,dofright,eq,0,then
*elseif,dofright,eq,1,then
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
dl,all,,uz
dl,all,,uY
allsel
*elseif,dofright,eq,2,then
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
dl,all,,uz
dl,all,,ux
allsel
*endif
finish
*elseif,hole,eq,2,then!!!!!!!!!!!!两点建模
/prep7
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
thforcal=1.1*thi
*elseif,type,eq,2,then
thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
thforcal=1.1*1.20*thi
*endif
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
csys,1
n,1,
n,2,kr,
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
et,1,shell63
real,2
mat,2
! :do loop1
*do,i,2,kld,1
e,1,i,i+1
*enddo
finish
/prep7
e,1,kld+1,2
*end
make
finish
/prep7
csys,0
k,1,1,1
k,2,1+w,1
*afun,deg
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,6,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
a,1,2,4,3
CYL4,kx(5),ky(5),kr
CYL4,kx(6),ky(6),kr
a***a,1,2
a***a,4,3
aplot
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
SMRT,6
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
lesize,all,,,hdiv
lsel,all
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
lesize,all,,,hdiv
lsel,all
LSEL,S,LENGTH,,2*kr,2*(h+w)
LSEL,INVE
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
mshkey,0
amesh,1
eplot
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6)
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
/auto,1
/rep
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
allsel
nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001
nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001
d,all,uy,0,,,,uz
allsel
CPINTF,ALL,0.0001,
*if,dofdown,eq,0,then
*elseif,dofdown,eq,1,then
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
dl,all,,uz
allsel
*elseif,dofdown,eq,2,then
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
dl,all,,uz
dl,all,,uy
allsel
*endif
finish
*elseif,hole,eq,4,then!!!!!!!!!!!!!!!!!四点建模
/prep7
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
thforcal=1.1*thi
*elseif,type,eq,2,then
thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
thforcal=1.1*1.20*thi
*endif
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
csys,1
n,1,
n,2,kr,
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
et,2,shell63
real,2
mat,2
! :do loop1
*do,i,2,kld,1
e,1,i,i+1
*enddo
finish
/prep7
e,1,kld+1,2
*end
make
finish
/prep7
csys,0
k,1,1,1
k,2,1+w,1
*afun,deg
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1)
k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2)
k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
a,1,2,4,3
CYL4,kx(5),ky(5),kr
CYL4,kx(6),ky(6),kr
CYL4,kx(7),ky(7),kr
CYL4,kx(8),ky(8),kr
a***a,1,2
a***a,6,3
a***a,1,4
a***a,2,5
aplot
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
SMRT,6
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
lesize,all,,,hdiv
lsel,all
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
lesize,all,,,hdiv
lsel,all
LSEL,S,LENGTH,,2*kr,l+w
LSEL,INVE
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
mshkey,0
amesh,1
eplot
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6)
*get,nnn3,node,0,count
egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7)
*get,nnn4,node,0,count
egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8)
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
/auto,1
/rep
finish
/solu
allsel
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
allsel
nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001
nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001
d,all,uy,0,,,,uz
allsel
nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001
nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001
d,all,ux,0,,,,uz
allsel
nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001
nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001
d,all,uz,0,,,,
allsel
CPINTF,ALL,0.0001,
finish
gplot
*elseif,hole,eq,6,then!!!!!!!!!!!!!!!!!六点建模
/prep7
/UIS,MSGPOP,3
et,1,shell63
*if,type,eq,1,then
thforcal=1.1*thi
*elseif,type,eq,2,then
thforcal=1.1*1.25*thi
*elseif,type,eq,3,then
thforcal=1.1*1.20*thi
*endif
r,2,thforcal
mp,ex,2,2.06e11
mp,prxy,2,0.3
csys,1
n,1,
n,2,kr,
kld=holediv
ngen,kld,1,2,,,,360/kld
*creat,make,mac
et,2,shell63
real,2
mat,2
! :do loop1
*do,i,2,kld,1
e,1,i,i+1
*enddo
finish
/prep7
e,1,kld+1,2
*end
make
finish
/prep7
csys,0
k,1,1,1
k,2,1+w,1
*afun,deg
k,3,1+H*cos(a1)/sin(a1),1+H
k,4,1+w-H*cos(a2)/sin(a2),1+H
k,5,kbw/sin(a1)+kbh*cos(a1)/sin(a1)+kx(1),kbh+ky(1)
k,6,kx(2)-kbw/sin(a2)+kbh*cos(a2)/sin(a2),kbh+ky(2)
k,7,kx(3)+kbw/sin(a1)-kbh*cos(a1)/sin(a1),ky(3)-kbh
k,8,kx(4)-kbw/sin(a2)+kbw*cos(a2)/sin(a2),ky(4)-kbh
hp=ky(7)-ky(5)
xx9=(hp-dbottom)*(kx(7)-kx(5))/hp
yy9=ky(5)+dbottom
xx10=(hp-dbottom)*(kx(6)-kx(8))/hp
yy10=ky(6)+dbottom
k,9,kx(7)-xx9,yy9
k,10,kx(8)+xx10,yy10
a,1,2,4,3
CYL4,kx(5),ky(5),kr
CYL4,kx(6),ky(6),kr
CYL4,kx(7),ky(7),kr
CYL4,kx(8),ky(8),kr
CYL4,kx(9),ky(9),kr
CYL4,kx(10),ky(10),kr
a***a,1,2
a***a,8,3
a***a,1,4
a***a,2,5
a***a,1,6
a***a,2,7
aplot
/UIS,MSGPOP,3
R,1,thforcal,
MP,EX,1,7.2e10
MP,PRXY,1,0.21
mp,dens,1,2560
SMRT,6
TYPE, 1
MAT, 1
REAL, 1
ESYS, 0
lsel,s,loc,y,ky(1)-0.01,ky(1)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,y,ky(3)-0.01,ky(3)+0.01
lesize,all,,,wdiv
lsel,all
lsel,s,loc,x,kx(1)-0.01,kx(3)+0.01
lesize,all,,,hdiv
lsel,all
lsel,s,loc,x,kx(4)-0.01,kx(2)+0.01
lesize,all,,,hdiv
lsel,all
LSEL,S,LENGTH,,2*kr,l+w
LSEL,INVE
lesize,all,,,kld/4
smrtsize,6
mshape,1,2d
mshkey,0
amesh,1
eplot
*get,nnn1,node,0,count
egen,2,nnn1+kld,1,kld,1,,,,,,kx(5),ky(5)
*get,nnn2,node,0,count
egen,2,nnn2+kld,1,kld,1,,,,,,kx(6),ky(6)
*get,nnn3,node,0,count
egen,2,nnn3+kld,1,kld,1,,,,,,kx(7),ky(7)
*get,nnn4,node,0,count
egen,2,nnn4+kld,1,kld,1,,,,,,kx(8),ky(8)
*get,nnn5,node,0,count
egen,2,nnn5+kld,1,kld,1,,,,,,kx(9),ky(9)
*get,nnn6,node,0,count
egen,2,nnn6+kld,1,kld,1,,,,,,kx(10),ky(10)
/UIS,MSGPOP,3
edele,1,kld,1
ndele,1,kld+1
/auto,1
/rep
finish
/solu
allsel
nsel,s,loc,x,kx(5)-0.001,kx(5)+0.001
nsel,r,loc,y,ky(5)-0.001,ky(5)+0.001
d,all,ux,0,,,,uy,uz
allsel
nsel,s,loc,x,kx(6)-0.001,kx(6)+0.001
nsel,r,loc,y,ky(6)-0.001,ky(6)+0.001
d,all,uy,0,,,,uz
allsel
nsel,s,loc,x,kx(7)-0.001,kx(7)+0.001
nsel,r,loc,y,ky(7)-0.001,ky(7)+0.001
d,all,ux,0,,,,uz
allsel
nsel,s,loc,x,kx(8)-0.001,kx(8)+0.001
nsel,r,loc,y,ky(8)-0.001,ky(8)+0.001
d,all,uz,0,,,,
allsel
nsel,s,loc,x,kx(9)-0.001,kx(9)+0.001
nsel,r,loc,y,ky(9)-0.001,ky(9)+0.001
d,all,ux,0,,,,uz
allsel
nsel,s,loc,x,kx(10)-0.001,kx(10)+0.001
nsel,r,loc,y,ky(10)-0.001,ky(10)+0.001
d,all,uz,0,,,,
allsel
CPINTF,ALL,0.0001,
finish
gplot
*endif
finish
/solu!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!求解部分
Ek=3*seismic*(thi+tho)*26.0*1000*9.8
w=1.4*wk
e=1.3*ek
q1=0.6*ek+wk
q2=w+0.6*e
finish
/solu
/UIS,MSGPOP,3
esel,s,real,,1,1,1
eplot
sfe,all,1,pres,,q1
acel,,gravity,
allsel
gplot
*if,switch,eq,1,then
nlgeom,off
kbc,1
*elseif,switch,eq,2,then
nlgeom,on
sstif,on
nsubst,10
kbc,1
*endif
solv
finish
/post1
esel,s,real,,1,1,1
AVPRIN
plnsol,u,sum,0,1
*get,numall,node,0,count
*get,numstart,node,0,num,min
umax=0
inuse=numstart
*do,i,1,numall,1
*get,usum,node,inuse,u,sum
*if,umax,ge,usum,then
umax=umax
*elseif,umax,lt,usum,then
umax=usum
*endif
inuse=ndnext(inuse)
*enddo
finish
/solu
/UIS,MSGPOP,3
esel,s,real,,1,1,1
eplot
sfe,all,1,pres,,q2
allsel
gplot
*if,switch,eq,1,then
nlgeom,off
kbc,1
solv
finish
/post1
esel,s,real,,1,1,1
AVPRIN
plnsol,s,eqv,0,1
allsel
esel,s,real,,1,1,1
NSLE,S
*get,numall,node,0,count
*get,numstart,node,0,num,min
smax=0
inuse=numstart
*do,i,1,numall,1
*get,ssum,node,inuse,s,eqv
*if,smax,ge,ssum,then
smax=smax
*elseif,smax,lt,ssum,then
smax=ssum
*endif
inuse=ndnext(inuse)
*enddo
allsel
*if,thforcal,le,0.012,then
dmq=84
kbq=58.8
*elseif,thforcal,gt,0.012,then
dmq=59
kbq=41.3
*endif
*if,h,lt,w,then
duanb=h
*elseif,h,ge,w,then
duanb=w
*endif
str=smax/1e6
esel,s,real,,1,1,1
/UIS,MSGPOP,1
*msg,note,str,umax*1000,dmq,10*duanb
Ultimate Limit States Glass STRESS=%gMa,%/&
Serviceabilitylimits DISPLACEMENT=%gmm,%/&
[Toughened Glass Allowble Stress]=%gMa,%/&
[Control Displacement]=%gmm,
/UIS,MSGPOP,3
*elseif,switch,eq,2,then
nlgeom,on
sstif,on
nsubst,10
kbc,1
solv
finish
/post1
esel,s,real,,1,1,1
AVPRIN
plnsol,s,eqv,0,1
allsel
esel,s,real,,1,1,1
NSLE,S
*get,numall,node,0,count
*get,numstart,node,0,num,min
smax=0
inuse=numstart
*do,i,1,numall,1
*get,ssum,node,inuse,s,eqv
*if,smax,ge,ssum,then
smax=smax
*elseif,smax,lt,ssum,then
smax=ssum
*endif
inuse=ndnext(inuse)
*enddo
allsel
*if,thforcal,le,0.012,then
dmq=84
kbq=58.8
*elseif,thforcal,gt,0.012,then
dmq=59
kbq=41.3
*endif
*if,h,lt,w,then
duanb=h
*elseif,h,ge,w,then
duanb=w
*endif
str=smax/1e6
esel,s,real,,1,1,1
/UIS,MSGPOP,1
*msg,note,str,umax*1000,dmq,10*duanb
Ultimate Limit States Glass STRESS=%gMa,%/&
Serviceabilitylimits DISPLACEMENT=%gmm,%/&
[Toughened Glass Allowble Stress]=%gMa,%/&
[Control Displacement]=%gmm,
/UIS,MSGPOP,3
*elseif,swtich,eq,0
allsel
gplot
*endif
/gopr
finish
/encrypt

那个红色的小鬼是比(5)大1个的（six）我输6就变成红色小鬼了。
那个音乐符号是（eight）

wangyongxin · 发表于 2005-11-29 17:16:36

大家贴的帖子都不错啊,受益菲浅!

显示全部楼层 · 发表于 2005-11-29 19:15:06

提示: 作者被禁止或删除内容自动屏蔽

luowl_hust · 发表于 2005-12-6 19:53:08

!混凝土板按照实体建模silod45，（先取为45单元，后期分析时采用65
!考虑预应力钢筋link8
!
!普通钢筋用8
/clear
/COM,ANSYS RELEASE  8.1 UP20040329    11:30:45 04/20/2005
/input,start81,ans,'d:\Program Files\Ansys Inc\v81\ANSYS\apdl\',,,,,,,,,,,,,,,,1
/filename, xxtbridge1
/title, xinxingtang bridge static analysis
/prep7
!
!混凝土65
et,1,solid65
!
et,2,45
!
ET,3,8
!
ET,4,10
!
r,
!
!直径为16
R,2,201.1,,
!直径为14
R,3,78.5, ,
!直径为8
R,4,50.3,,
!钢绞线，每束暂定为5－15.24
R,5,700,,
!
!混凝土材料特性及本构关系材料01
MP,EX,1,32845.5
MP,PRXY,1,0.167
MP,DENS,1,0.0025
TB,MELA,1,1,7,
TBPT,,0.000246,8.08
TBPT,,0.000595,17.727
TBPT,,0.000943,25.224
TBPT,,0.00130,30.712
TBPT,,0.0016,33.6
TBPT,,0.002,35
TBPT,,0.0033,29.75
TB,CONC,1,1,9,
TBTEMP,0
TBDATA,,0.25,1,3,-1,0,0
TBDATA,,0,0,0.6,,,
!材料2
MP,EX,2,33000
MP,PRXY,2,0.167
MP,DENS,2,0.0025
!材料3
!
mp,ex,3,200000
mp,prxy,3,0.3
mp,dens,3,0.00785
TB,BISO,3,1,2,
TBTEMP,0
TBDATA,,340,,,,,
!材料4
!
mp,ex,4,190909.1
mp,prxy,4,0.3
MP,ALPX,4,0.00001
mp,dens,4,0.00785
TB,MELA,4,1,3,
TBPT,,0.0066,1260
TBPT,,0.01,1570
TBPT,,0.035,1860
!材料5
!
mp,ex,5,200000
mp,prxy,5,0.3
mp,dens,5,0.00785
TB,BISO,5,1,2,
TBTEMP,0
TBDATA,,240,,,,,
!
!从8
!y
!!x
!z
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!00!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!顶板层结点1~41
!y=0以及y=16300
n,1,0-6400,0,0
n,2,42-6400,0,0
n,3,202-6400,0,0
n,19,2762-6400,0,0
fill
n,20,2960-6400,0,0
n,21,3059-6400,0,0
n,22,3158-6400,0,0
n,23,3386-6400,0,0
n,24,3580-6400,0,0
n,25,3790-6400,0,0
n,26,4000-6400,0,0
n,41,0,0,0
fill
!第二层结点42~82,z坐标减少42
!
ngen,2,41,1,41,,,,-42
!第三层83~123，z坐标减小58
n,83,0-6400,0,-100
n,84,42-6400,0,-101.9
n,85,202-6400,0,-109.3
n,86,362-6400,0,-116.6
n,87,522-6400,0,-124
n,88,682-6400,0,-131.3
n,89,842-6400,0,-138.7
n,90,1002-6400,0,-146
n,91,1162-6400,0,-153.3
n,92,1322-6400,0,-160.7
n,93,1482-6400,0,-168
n,94,1642-6400,0,-175.4
n,95,1802-6400,0,-182.7
n,96,1962-6400,0,-190.1
n,97,2122-6400,0,-198
n,98,2282-6400,0,-204.8
n,99,2442-6400,0,-212.1
n,100,2602-6400,0,-219.4
n,101,2762-6400,0,-226.8
n,102,2960-6400,0,-235.9
n,103,3120.4-6400,0,-173.8
n,104,3192.3-6400,0,-145
n,105,3420.3-6400,0,-145
n,106,3579.9-6400,0,-145
ngen,2,82,25,41,,,,-145
!第四层结点124~164，对应于下层钢筋，需要计算z
n,124,0-6400,0,-158
n,125,42-6400,0,-161.1
n,126,202-6400,0,-173.4
n,142,2762-6400,0,-371
fill
n,143,2960-6400,0,-386.2
n,144,3169.3-6400,0,-277.7
n,145,3226.7-6400,0,-248
n,146,3454.7-6400,0,-248
n,147,3579.8-6400,0,-248
n,148,3790-6400,0,-248
n,149,4000-6400,0,-248
n,150,4176.4-6400,0,-248
ngen,2,123,28,41,,,,-248
!第五层结点，对应于箱梁下表面，需要计算z
n,165,0-6400,0,-200
n,166,42-6400,0,-203.9
n,167,202-6400,0,-216.2
n,183,2762-6400,0,-412.5
fill
n,184,2960-6400,0,-427.7
n,185,3250-6400,0,-450
n,186,3295.2-6400,0,-453.7
n,187,3529.2-6400,0,-471.5
n,188,3579.6-6400,0,-490
n,189,3790-6400,0,-419.9
n,190,4000-6400,0,-349.9
n,191,4179.8-6400,0,-290
ngen,2,164,28,41,,,,-290
!!!
n,206,3276-6400,0,-528
n,207,3320-6400,0,-528
n,208,3548-6400,0,-528
n,209,3592.2-6400,0,-528
n,242,3733.3-6400,0,-1900
n,243,3777.3-6400,0,-1900
n,244,4005.3-6400,0,-1900
n,245,4049.6-6400,0,-1900
*do,i,1,4
fill,205+i,241+i,8,209+i,4
*enddo
n,246,3782.7-6400,0,-2048
n,247,3826.7-6400,0,-2048
n,248,4054.7-6400,0,-2048
n,249,4197.6-6400,0,-2048
n,250,3847.3-6400,0,-2242
n,251,3891.3-6400,0,-2242
n,252,4119.3-6400,0,-2242
n,253,4245.8-6400,0,-2242
n,254,3866.7-6400,0,-2300
n,255,3910.7-6400,0,-2300
n,256,4138.7-6400,0,-2300
n,257,4260.2-6400,0,-2300
n,258,3886.0-6400,0,-2358
n,259,3930.0-6400,0,-2358
n,260,4158.0-6400,0,-2358
n,261,4274.6-6400,0,-2358
n,262,3900.0-6400,0,-2400
n,263,3944.0-6400,0,-2400
n,264,4172.0-6400,0,-2400
n,265,4285.0-6400,0,-2400
!!!
n,266,4350.0-6400,0,-2200
n,279,6400-6400,0,-2200
fill
n,280,4349.6-6400,0,-2242
n,293,6400-6400,0,-2242
fill
n,294,4349.6-6400,0,-2300
n,307,6400-6400,0,-2300
fill
n,308,4349.6-6400,0,-2358
n,321,6400-6400,0,-2358
fill
n,322,4349.6-6400,0,-2400
n,335,6400-6400,0,-2400
fill
!
!nsym,,400,1,335
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!16300!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为16500!!!!!!!!!!!!!!!!!!!!!!!!!!!1
!插入0～16300
ngen,34,1000,1,999,,,500,
!y坐标为0时结点编号为1～335，，401，735
!y坐标为16500时结点编号为30001～30335，，30401，30735
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!20100!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为20000!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为20100时，增加6－6
!!
ngen,2,40000,1,22,,,20000
ngen,2,40000,30,63,,,20000
ngen,2,40000,71,104,,,20000
ngen,2,40000,112,145,,,20000
ngen,2,40000,153,186,,,20000
ngen,2,40000,194,205,,,20000
n,40023,3450-6400,20000,0
n,40024,3580-6400,20000,0
n,40025,3800-6400,20000,0
n,40026,4000-6400,20000,0
n,40027,4160-6400,20000,0
n,40028,4320-6400,20000,0
n,40029,4480-6400,20000,0
n,40023+41,3450-6400,20000,-42
n,40024+41,3580-6400,20000,-42
n,40025+41,3800-6400,20000,-42
n,40026+41,4000-6400,20000,-42
n,40027+41,4172-6400,20000,-42
n,40028+41,4320-6400,20000,-42
n,40029+41,4480-6400,20000,-42
n,40023+82,3484.3-6400,20000,-145
n,40024+82,3594.7-6400,20000,-145
n,40025+82,3800-6400,20000,-145
n,40026+82,4000-6400,20000,-145
n,40027+82,4201.9-6400,20000,-145
n,40028+82,4320-6400,20000,-145
n,40029+82,4480-6400,20000,-145
n,40023+123,3518.7-6400,20000,-248
n,40024+123,3609.4-6400,20000,-248
n,40025+123,3800-6400,20000,-248
n,40026+123,4000-6400,20000,-248
n,40027+123,4231.7-6400,20000,-248
n,40028+123,4320-6400,20000,-248
n,40029+123,4480-6400,20000,-248
n,40023+164,3597.8-6400,20000,-485.3
n,40024+164,3643.9-6400,20000,-490
n,40025+164,3800-6400,20000,-438
n,40026+164,4000-6400,20000,-371.3
n,40027+164,4243.9-6400,20000,-290
n,40028+164,4320-6400,20000,-290
n,40029+164,4480-6400,20000,-290
!!!
n,40206,3276-6400,20000,-528
n,40207,3320-6400,20000,-528
n,40208,3612-6400,20000,-528
n,40209,3656.5-6400,20000,-528
n,40242,3715.3-6400,20000,-1846
n,40243,3759.3-6400,20000,-1846
n,40244,4051.3-6400,20000,-1846
n,40245,4095.9-6400,20000,-1846
*do,i,1,4
fill,40205+i,40241+i,8,40209+i,4
*enddo
n,40246,3782.7-6400,20000,-2048
n,40247,3826.7-6400,20000,-2048
n,40248,4118.7-6400,20000,-2048
n,40249,4297.8-6400,20000,-2048
n,40250,3829.3-6400,20000,-2188
n,40251,3873.3-6400,20000,-2188
n,40252,4165.3-6400,20000,-2188
n,40253,4297.8-6400,20000,-2188
n,40254,3857.7-6400,20000,-2273
n,40255,3901.7-6400,20000,-2273
n,40256,4193.7-6400,20000,-2273
n,40257,4297.8-6400,20000,-2273
n,40258,3886.0-6400,20000,-2358
n,40259,3930.0-6400,20000,-2358
n,40260,4222.0-6400,20000,-2358
n,40261,4297.8-6400,20000,-2358
n,40262,3900.0-6400,20000,-2400
n,40263,3944.0-6400,20000,-2400
n,40264,4236.0-6400,20000,-2400
n,40265,4297.8-6400,20000,-2400
!!!
n,40266,4395.8-6400,20000,-2146
n,40279,6400-6400,20000,-2146
fill
n,40280,4395.8-6400,20000,-2188
n,40293,6400-6400,20000,-2188
fill
n,40294,4395.8-6400,20000,-2273
n,40307,6400-6400,20000,-2273
fill
n,40308,4395.8-6400,20000,-2358
n,40321,6400-6400,20000,-2358
fill
n,40322,4395.8-6400,20000,-2400
n,40335,6400-6400,20000,-2400
fill
!nsym,,400,40001,40335
!插入16300到20000
*do,i,1,335
fill,33000+i,40000+i,6,33000+1000+i,1000
*enddo
!*do,i,401,735
!fill,33000+i,40000+i,6,33000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!23600!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为23500!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为23500时，增加7－7
!!!
ngen,2,47000,1,22,,,23500
ngen,2,47000,30,63,,,23500
ngen,2,47000,71,104,,,23500
ngen,2,47000,112,145,,,23500
ngen,2,47000,153,186,,,23500
ngen,2,47000,194,205,,,23500
n,47023,3509-6400,23500,0
n,47024,3660-6400,23500,0
n,47025,3830-6400,23500,0
n,47026,4000-6400,23500,0
n,47027,4129-6400,23500,0
n,47028,4379-6400,23500,0
n,47029,4480-6400,23500,0
n,47023+41,3509-6400,23500,-42
n,47024+41,3660-6400,23500,-42
n,47025+41,3830-6400,23500,-42
n,47026+41,4000-6400,23500,-42
n,47027+41,4231-6400,23500,-42
n,47028+41,4379-6400,23500,-42
n,47029+41,4480-6400,23500,-42
n,47023+82,3543.3-6400,23500,-145
n,47024+82,3669.9-6400,23500,-145
n,47025+82,3830-6400,23500,-145
n,47026+82,4000-6400,23500,-145
n,47027+82,4260.9-6400,23500,-145
n,47028+82,4379-6400,23500,-145
n,47029+82,4480-6400,23500,-145
n,47023+123,3577.7-6400,23500,-248
n,47024+123,3679.7-6400,23500,-248
n,47025+123,3830-6400,23500,-248
n,47026+123,4000-6400,23500,-248
n,47027+123,4290.8-6400,23500,-248
n,47028+123,4379-6400,23500,-248
n,47029+123,4480-6400,23500,-248
n,47023+164,3656.7-6400,23500,-485.1
n,47024+164,3702.9-6400,23500,-490
n,47025+164,3830-6400,23500,-447.6
n,47026+164,4000-6400,23500,-391.0
n,47027+164,4302.9-6400,23500,-290
n,47028+164,4379-6400,23500,-290
n,47029+164,4480-6400,23500,-290
!!!
n,47206,3276-6400,23500,-528
n,47207,3320-6400,23500,-528
n,47208,3671-6400,23500,-528
n,47209,3715.5-6400,23500,-528
n,47242,3699-6400,23500,-1797
n,47243,3743-6400,23500,-1797
n,47244,4094-6400,23500,-1797
n,47245,4138.6-6400,23500,-1797
*do,i,1,4
fill,47205+i,47241+i,8,47209+i,4
*enddo
n,47246,3782.7-6400,23500,-2048
n,47247,3826.7-6400,23500,-2048
n,47248,4177.7-6400,23500,-2048
n,47249,4389.6-6400,23500,-2048
n,47250,3813-6400,23500,-2139
n,47251,3857-6400,23500,-2139
n,47252,4208-6400,23500,-2139
n,47253,4380.5-6400,23500,-2139
n,47254,3849.5-6400,23500,-2248.5
n,47255,3893.5-6400,23500,-2248.5
n,47256,4244.5-6400,23500,-2248.5
n,47257,4370.0-6400,23500,-2248.5
n,47258,3886.0-6400,23500,-2358
n,47259,3930.0-6400,23500,-2358
n,47260,4281.0-6400,23500,-2358
n,47261,4358.9-6400,23500,-2358
n,47262,3900.0-6400,23500,-2400
n,47263,3944.0-6400,23500,-2400
n,47264,4295.0-6400,23500,-2400
n,47265,4355.0-6400,23500,-2400
!!!
n,47266,4438.6-6400,23500,-2097
n,47279,6400-6400,23500,-2097
fill
n,47280,4438.6-6400,23500,-2139
n,47293,6400-6400,23500,-2139
fill
n,47294,4438.6-6400,23500,-2248.5
n,47307,6400-6400,23500,-2248.5
fill
n,47308,4438.6-6400,23500,-2358
n,47321,6400-6400,23500,-2358
fill
n,47322,4438.6-6400,23500,-2400
n,47335,6400-6400,23500,-2400
fill
!nsym,,400,47001,47335
!插入20000到23500
*do,i,1,335
fill,40000+i,47000+i,6,40000+1000+i,1000
*enddo
!*do,i,401,735
!fill,40000+i,47000+i,6,40000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!27100!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为27000!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为27100时，增加8－8
!!!
ngen,2,54000,1,22,,,27000
ngen,2,54000,30,63,,,27000
ngen,2,54000,71,104,,,27000
ngen,2,54000,112,145,,,27000
ngen,2,54000,153,186,,,27000
ngen,2,54000,194,205,,,27000
n,54023,3568-6400,27000,0
n,54024,3660-6400,27000,0
n,54025,3830-6400,27000,0
n,54026,4030-6400,27000,0
n,54027,4278-6400,27000,0
n,54028,4408-6400,27000,0
n,54029,4480-6400,27000,0
n,54023+41,3568-6400,27000,-42
n,54024+41,3660-6400,27000,-42
n,54025+41,3830-6400,27000,-42
n,54026+41,4030-6400,27000,-42
n,54027+41,4290.2-6400,27000,-42
n,54028+41,4408-6400,27000,-42
n,54029+41,4480-6400,27000,-42
n,54023+82,3602.3-6400,27000,-145
n,54024+82,3683.4-6400,27000,-145
n,54025+82,3830-6400,27000,-145
n,54026+82,4030-6400,27000,-145
n,54027+82,4320-6400,27000,-145
n,54028+82,4408-6400,27000,-145
n,54029+82,4480-6400,27000,-145
n,54023+123,3636.7-6400,27000,-248
n,54024+123,3706.9-6400,27000,-248
n,54025+123,3830-6400,27000,-248
n,54026+123,4030-6400,27000,-248
n,54027+123,4349.8-6400,27000,-248
n,54028+123,4408-6400,27000,-248
n,54029+123,4480-6400,27000,-248
n,54023+164,3716.1-6400,27000,-486.3
n,54024+164,3761.9-6400,27000,-490
n,54025+164,3830-6400,27000,-467.3
n,54026+164,4030-6400,27000,-400.6
n,54027+164,4361.9-6400,27000,-290
n,54028+164,4408-6400,27000,-290
n,54029+164,4480-6400,27000,-290
!!!
n,54206,3286-6400,27000,-558
n,54207,3330-6400,27000,-558
n,54208,3740-6400,27000,-558
n,54209,3784.6-6400,27000,-558
n,54242,3682.7-6400,27000,-1748
n,54243,3726.7-6400,27000,-1748
n,54244,4136.7-6400,27000,-1748
n,54245,4181.3-6400,27000,-1748
*do,i,1,4
fill,54205+i,54241+i,8,54209+i,4
*enddo
n,54246,3752.0-6400,27000,-1956.2
n,54247,3796.1-6400,27000,-1956.2
n,54248,4206.1-6400,27000,-1956.2
n,54249,4389.4-6400,27000,-1956.2
n,54250,3813-6400,27000,-2139
n,54251,3857-6400,27000,-2139
n,54252,4267-6400,27000,-2139
n,54253,4400.8-6400,27000,-2139
n,54254,3849.5-6400,27000,-2248.5
n,54255,3893.5-6400,27000,-2248.5
n,54256,4303.5-6400,27000,-2248.5
n,54257,4407.7-6400,27000,-2248.5
n,54258,3886.0-6400,27000,-2358
n,54259,3930.0-6400,27000,-2358
n,54260,4340-6400,27000,-2358
n,54261,4414.5-6400,27000,-2358
n,54262,3900.0-6400,27000,-2400
n,54263,3944.0-6400,27000,-2400
n,54264,4354-6400,27000,-2400
n,54265,4417.1-6400,27000,-2400
!!!
n,54266,4481.3-6400,27000,-2048
n,54279,6400-6400,27000,-2048
fill
n,54280,4481.3-6400,27000,-2139
n,54293,6400-6400,27000,-2139
fill
n,54294,4481.3-6400,27000,-2248.5
n,54307,6400-6400,27000,-2248.5
fill
n,54308,4481.3-6400,27000,-2358
n,54321,6400-6400,27000,-2358
fill
n,54322,4481.3-6400,27000,-2400
n,54335,6400-6400,27000,-2400
fill
!nsym,,400,54001,54335
!插入23500到27000
*do,i,1,335
fill,47000+i,54000+i,6,47000+1000+i,1000
*enddo
!*do,i,401,735
!fill,47000+i,54000+i,6,47000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!30500!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为30500时，增加9－9
!!!
ngen,2,61000,1,22,,,30500
ngen,2,61000,30,63,,,30500
ngen,2,61000,71,104,,,30500
ngen,2,61000,112,145,,,30500
ngen,2,61000,153,186,,,30500
ngen,2,61000,194,205,,,30500
n,61023,3626-6400,30500,0
n,61024,3699.3-6400,30500,0
n,61025,3927.1-6400,30500,0
n,61026,4127.1-6400,30500,0
n,61027,4336-6400,30500,0
n,61028,4466-6400,30500,0
n,61029,4538-6400,30500,0
n,61023+41,3626-6400,30500,-42
n,61024+41,3699.3-6400,30500,-42
n,61025+41,3927.1-6400,30500,-42
n,61026+41,4127.1-6400,30500,-42
n,61027+41,4348.2-6400,30500,-42
n,61028+41,4466-6400,30500,-42
n,61029+41,4538-6400,30500,-42
n,61023+82,3660.3-6400,30500,-145
n,61024+82,3727-6400,30500,-145
n,61025+82,3927.1-6400,30500,-145
n,61026+82,4127.1-6400,30500,-145
n,61027+82,4377.9-6400,30500,-145
n,61028+82,4466-6400,30500,-145
n,61029+82,4538-6400,30500,-145
n,61023+123,3694.7-6400,30500,-248
n,61024+123,3754.8-6400,30500,-248
n,61025+123,3927.1-6400,30500,-248
n,61026+123,4127.1-6400,30500,-248
n,61027+123,4407.8-6400,30500,-248
n,61028+123,4466-6400,30500,-248
n,61029+123,4538-6400,30500,-248
n,61023+164,3774.4-6400,30500,-487
n,61024+164,3819.9-6400,30500,-490
n,61025+164,3927.1-6400,30500,-454.3
n,61026+164,4127.1-6400,30500,-387.6
n,61027+164,4419.9-6400,30500,-290
n,61028+164,4466-6400,30500,-290
n,61029+164,4538-6400,30500,-290
!!!
n,61206,3286-6400,30500,-558
n,61207,3330-6400,30500,-558
n,61208,3798-6400,30500,-558
n,61209,3842.6-6400,30500,-558
n,61242,3666.7-6400,30500,-1700
n,61243,3710.7-6400,30500,-1700
n,61244,4178.7-6400,30500,-1700
n,61245,4223.2-6400,30500,-1700
*do,i,1,4
fill,61205+i,61241+i,8,61209+i,4
*enddo
n,61246,3727.5-6400,30500,-1882.6
n,61247,3771.5-6400,30500,-1882.6
n,61248,4239.5-6400,30500,-1882.6
n,61249,4405.8-6400,30500,-1882.6
n,61250,3780.7-6400,30500,-2042
n,61251,3824.7-6400,30500,-2042
n,61252,4292.7-6400,30500,-2042
n,61253,4425.2-6400,30500,-2042
n,61254,3833.2-6400,30500,-2199.5
n,61255,3877.2-6400,30500,-2199.5
n,61256,4345.2-6400,30500,-2199.5
n,61257,4444.4-6400,30500,-2199.5
n,61258,3886.0-6400,30500,-2358
n,61259,3930.0-6400,30500,-2358
n,61260,4398-6400,30500,-2358
n,61261,4468.9-6400,30500,-2358
n,61262,3900.0-6400,30500,-2400
n,61263,3944.0-6400,30500,-2400
n,61264,4412-6400,30500,-2400
n,61265,4468.9-6400,30500,-2400
!!!
n,61266,4523.2-6400,30500,-2000
n,61279,6400-6400,30500,-2000
fill
n,61280,4523.2-6400,30500,-2042
n,61293,6400-6400,30500,-2042
fill
n,61294,4523.2-6400,30500,-2200
n,61307,6400-6400,30500,-2200
fill
n,61308,4523.2-6400,30500,-2358
n,61321,6400-6400,30500,-2358
fill
n,61322,4523.2-6400,30500,-2400
n,61335,6400-6400,30500,-2400
fill
!nsym,,400,61001,61335
!插入27000到30500
*do,i,1,335
fill,54000+i,61000+i,6,54000+1000+i,1000
*enddo
!*do,i,401,735
!fill,54000+i,61000+i,6,54000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!32500!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为32500时，增加10－10
!!!
ngen,2,65000,1,22,,,32500
ngen,2,65000,30,63,,,32500
ngen,2,65000,71,104,,,32500
ngen,2,65000,112,145,,,32500
ngen,2,65000,153,186,,,32500
ngen,2,65000,194,205,,,32500
n,65023,3626-6400,32500,0
n,65024,3699.3-6400,32500,0
n,65025,3927.1-6400,32500,0
n,65026,4127.1-6400,32500,0
n,65027,4336-6400,32500,0
n,65028,4466-6400,32500,0
n,65029,4538-6400,32500,0
n,65023+41,3626-6400,32500,-42
n,65024+41,3699.3-6400,32500,-42
n,65025+41,3927.1-6400,32500,-42
n,65026+41,4127.1-6400,32500,-42
n,65027+41,4348.2-6400,32500,-42
n,65028+41,4466-6400,32500,-42
n,65029+41,4538-6400,32500,-42
n,65023+82,3660.3-6400,32500,-145
n,65024+82,3727-6400,32500,-145
n,65025+82,3927.1-6400,32500,-145
n,65026+82,4127.1-6400,32500,-145
n,65027+82,4377.9-6400,32500,-145
n,65028+82,4466-6400,32500,-145
n,65029+82,4538-6400,32500,-145
n,65023+123,3694.7-6400,32500,-248
n,65024+123,3754.8-6400,32500,-248
n,65025+123,3927.1-6400,32500,-248
n,65026+123,4127.1-6400,32500,-248
n,65027+123,4407.8-6400,32500,-248
n,65028+123,4466-6400,32500,-248
n,65029+123,4538-6400,32500,-248
n,65023+164,3774.4-6400,32500,-487
n,65024+164,3819.9-6400,32500,-490
n,65025+164,3927.1-6400,32500,-454.3
n,65026+164,4127.1-6400,32500,-387.6
n,65027+164,4419.9-6400,32500,-290
n,65028+164,4466-6400,32500,-290
n,65029+164,4538-6400,32500,-290
!!!
n,65206,3286-6400,32500,-558
n,65207,3330-6400,32500,-558
n,65208,3798-6400,32500,-558
n,65209,3842.6-6400,32500,-558
n,65242,3666.7-6400,32500,-1700
n,65243,3710.7-6400,32500,-1700
n,65244,4178.7-6400,32500,-1700
n,65245,4223.2-6400,32500,-1700
*do,i,1,4
fill,65205+i,65241+i,8,65209+i,4
*enddo
n,65246,3727.5-6400,32500,-1882.6
n,65247,3771.5-6400,32500,-1882.6
n,65248,4239.5-6400,32500,-1882.6
n,65249,4405.8-6400,32500,-1882.6
n,65250,3780.7-6400,32500,-2042
n,65251,3824.7-6400,32500,-2042
n,65252,4292.7-6400,32500,-2042
n,65253,4425.2-6400,32500,-2042
n,65254,3833.2-6400,32500,-2199.5
n,65255,3877.2-6400,32500,-2199.5
n,65256,4345.2-6400,32500,-2199.5
n,65257,4444.4-6400,32500,-2199.5
n,65258,3886.0-6400,32500,-2358
n,65259,3930.0-6400,32500,-2358
n,65260,4398-6400,32500,-2358
n,65261,4468.9-6400,32500,-2358
n,65262,3900.0-6400,32500,-2400
n,65263,3944.0-6400,32500,-2400
n,65264,4412-6400,32500,-2400
n,65265,4468.9-6400,32500,-2400
!!!
n,65266,4523.2-6400,32500,-2000
n,65279,6400-6400,32500,-2000
fill
n,65280,4523.2-6400,32500,-2042
n,65293,6400-6400,32500,-2042
fill
n,65294,4523.2-6400,32500,-2200
n,65307,6400-6400,32500,-2200
fill
n,65308,4523.2-6400,32500,-2358
n,65321,6400-6400,32500,-2358
fill
n,65322,4523.2-6400,32500,-2400
n,65335,6400-6400,32500,-2400
fill
!nsym,,400,65001,65335
!插入30500到32500
*do,i,1,335
fill,61000+i,65000+i,3,61000+1000+i,1000
*enddo
!*do,i,401,735
!fill,61000+i,65000+i,3,61000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!34500!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为34500时，增加11－11截面结点，与9－9截面一样,由9截面生成11
ngen,2,8000,61001,61999,,,4000,,
!
!插入32500到34500
*do,i,1,335
fill,65000+i,69000+i,3,65000+1000+i,1000
*enddo
!*do,i,401,735
!fill,65000+i,69000+i,3,65000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!37900!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为38000!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为37900时，增加12－12截面结点，与8－8截面一样,由8截面生成12
ngen,2,22000,54001,54999,,,11000,,
!
!插入34500到37900
*do,i,1,335
fill,69000+i,76000+i,6,69000+1000+i,1000
*enddo
!*do,i,401,735
!fill,69000+i,76000+i,6,69000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!41400!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为41400时，增加13－13截面结点，与7－7截面一样,由7截面生成11
ngen,2,36000,47001,47999,,,18000,,
!
!插入37900到41400
*do,i,1,335
fill,76000+i,83000+i,6,76000+1000+i,1000
*enddo
!*do,i,401,735
!fill,76000+i,83000+i,6,76000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!44900!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!调整为45000!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为44900时，增加14－14截面结点，与6－6截面一样,由6截面生成14
ngen,2,50000,40001,40999,,,25000,,
!
!插入41400到44900
*do,i,1,335
fill,83000+i,90000+i,6,83000+1000+i,1000
*enddo
!*do,i,401,735
!fill,83000+i,90000+i,6,83000+1000+i,1000
!*enddo
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!57500!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!y坐标为57500时，增加18－18截面结点，与1－1截面一样,由1截面生成18
ngen,2,115000,1,999,,,57500,,
!
!插入44900到57500
*do,i,1,335
fill,90000+i,115000+i,24,90000+1000+i,1000
*enddo
!*do,i,401,735
!fill,90000+i,115000+i,24,90000+1000+i,1000
!*enddo
cs,100,cart,115041,115001,115082
nsym,z,1000000,1,115999
!murge,
!!!
!!
*get,e_num,elem,,num,max
type,2
mat,2
real,1
e,1,2,43,42,1001,1002,1043,1042
egen,40,1,e_num+1,e_num+1
egen,4,41,e_num+1,e_num+40
*get,e_num,elem,,num,max
e,185,186,207,206,1185,1186,1207,1206
egen,3,1,e_num+1,e_num+1
*get,e_num,elem,,num,max
e,206,207,211,210,1206,1207,1211,1210
egen,3,1,e_num+1,e_num+1
egen,14,4,e_num+1,e_num+3
*get,e_num,elem,,num,max
e,249,266,280,253,1249,1266,1280,1253
e,253,280,294,257,1253,1280,1294,1257
e,257,294,308,261,1257,1294,1308,1261
e,261,308,322,265,1261,1308,1322,1265
*get,e_num,elem,,num,max
e,266,267,281,280,1266,1267,1281,1280
egen,13,1,e_num+1,e_num+1
egen,4,14,e_num+1,e_num+13
egen,2,1000,1,261
egen,2,64000,1,522
egen,2,111000,1,522
egen,2,113000,1,522
*get,e_num,elem,,num,max
type,1
mat,1
real,1
e,2001,2002,2043,2042,2000+1001,2000+1002,2000+1043,2000+1042
egen,40,1,e_num+1,e_num+1
egen,4,41,e_num+1,e_num+40
*get,e_num,elem,,num,max
e,2000+185,2000+186,2000+207,2000+206,2000+1185,2000+1186,2000+1207,2000+1206
egen,3,1,e_num+1,e_num+1
*get,e_num,elem,,num,max
e,2000+206,2000+207,2000+211,2000+210,2000+1206,2000+1207,2000+1211,2000+1210
egen,3,1,e_num+1,e_num+1
egen,14,4,e_num+1,e_num+3
*get,e_num,elem,,num,max
!type,2
!mat,2
e,2000+249,2000+266,2000+280,2000+253,2000+1249,2000+1266,2000+1280,2000+1253
e,2000+253,2000+280,2000+294,2000+257,2000+1253,2000+1280,2000+1294,2000+1257
e,2000+257,2000+294,2000+308,2000+261,2000+1257,2000+1294,2000+1308,2000+1261
e,2000+261,2000+308,2000+322,2000+265,2000+1261,2000+1308,2000+1322,2000+1265
*get,e_num,elem,,num,max
!type,1
!mat,1
e,2000+266,2000+267,2000+281,2000+280,2000+1266,2000+1267,2000+1281,2000+1280
egen,13,1,e_num+1,e_num+1
egen,4,14,e_num+1,e_num+13
*get,e_num,elem,,num,max
egen,62,1000,e_num-260,e_num
egen,2,3000,18270-260,18270
egen,45,1000,18271,18531
!单元反射！！！！！！！！！！!!!!!!!!
!!!ENSYM, , ,400,all
!
!直径为16
!1#
type,3
mat,3
real,2
*get,e_num,elem,,num,max
e,43,44
egen,39,1,e_num+1,e_num+1
egen,115,1000,e_num+1,e_num+39
!2#
*get,e_num,elem,,num,max
e,84,85
egen,39,1,e_num+1,e_num+1
egen,115,1000,e_num+1,e_num+39
!3#
*get,e_num,elem,,num,max
e,125,126
egen,39,1,e_num+1,e_num+1
egen,115,1000,e_num+1,e_num+39
!5#
*get,e_num,elem,,num,max
e,63,104
egen,3,41,e_num+1,e_num+1
e,186,207
e,207,211
egen,13,4,e_num+5,e_num+5
egen,2,1,e_num+1,e_num+17
egen,115,1000,e_num+1,e_num+34
!6－1#
*get,e_num,elem,,num,max
e,251,252
e,252,253
e,253,280
e,255,256
e,256,257
e,257,294
e,259,260
e,260,261
e,261,308
e,280,281
egen,13,1,e_num+10,e_num+10
egen,3,14,e_num+10,e_num+22
egen,115,1000,e_num+1,e_num+48
!直径为10
!
!1#
type,3
mat,5
real,3
!
*get,e_num,elem,,num,max
e,43,1043
egen,115,1000,e_num+1,e_num+1
egen,40,1,e_num+1,e_num+115
!
*get,e_num,elem,,num,max
e,125,1125
egen,115,1000,e_num+1,e_num+1
egen,40,1,e_num+1,e_num+115
!
!!!
*get,e_num,elem,,num,max
e,251,1251
e,252,1252
e,253,1253
e,259,1259
e,260,1260
e,261,1261
egen,115,1000,e_num+1,e_num+6
!!!
*get,e_num,elem,,num,max
e,280,1280
e,308,1308
egen,115,1000,e_num+1,e_num+2
egen,13,1,e_num+1,e_num+230
SAVE
!直径为8
type,3
mat,3
real,4
*get,e_num,elem,,num,max
*do,i,0,10
e,1207+4*i,2207+4*i
e,1208+4*i,2208+4*i
*enddo
egen,114,1000,e_num+1,e_num+22
!
!1#
type,3
mat,3
real,4
*get,e_num,elem,,num,max
e,207,208
e,215,216
e,223,224
e,231,232
e,239,240
e,243,244
e,244,245
egen,115,1000,e_num+1,e_num+7
!!!!!!!!!!!!!!!!!!!!!!!生成钢绞线!!!!!!!!!!!!11111
type,4
mat,4
real,5
!
!B1钢绞线，76～54
*get,e_num,elem,,num,max
e,54123,55123
e,54118,55118
e,54109,55109
e,54101,55101
e,54096,55096
e,54091,55091
egen,22,1000,e_num+1,e_num+6
!B2钢绞线，47～83
*get,e_num,elem,,num,max
e,47122,48122
e,47116,48116
e,47098,48098
e,47092,48092
egen,36,1000,e_num+1,e_num+4
!B3钢绞线，40～90
*get,e_num,elem,,num,max
e,40120,41120
e,40115,41115
e,40107,41107
e,40093,41093
egen,50,1000,e_num+1,e_num+4
!B4钢绞线，33～97
*get,e_num,elem,,num,max
e,33119,34119
e,33113,34113
e,33103,34103
e,33095,34095
egen,64,1000,e_num+1,e_num+4
!B5钢绞线，25～105
*get,e_num,elem,,num,max
e,25112,26112
e,25108,26108
e,25099,26099
egen,80,1000,e_num+1,e_num+3
!B6钢绞线，17～111
*get,e_num,elem,,num,max
e,17114,18114
e,17110,18110
e,17106,18106
e,17102,18102
egen,94,1000,e_num+1,e_num+4
!C钢绞线，1～115
*get,e_num,elem,,num,max
e,1105,2105
egen,114,1000,e_num+1,e_num+1
!
!D1钢绞线，17～33
*get,e_num,elem,,num,max
e,17307,18307
egen,16,1000,e_num+1,e_num+1
!D2钢绞线，9～33
*get,e_num,elem,,num,max
e,9294,10294
egen,24,1000,e_num+1,e_num+1
!D3钢绞线，9～40
*get,e_num,elem,,num,max
e,9305,10305
egen,31,1000,e_num+1,e_num+1
!D4钢绞线，5～40
*get,e_num,elem,,num,max
e,5304,6304
egen,35,1000,e_num+1,e_num+1
!D5钢绞线，1～47
*get,e_num,elem,,num,max
e,300,1300
e,297,1297
egen,47,1000,e_num+1,e_num+2
!D6钢绞线，1～54
*get,e_num,elem,,num,max
e,302,1302
e,299,1299
e,295,1295
egen,54,1000,e_num+1,e_num+3
!E1钢绞线，105～115
*get,e_num,elem,,num,max
e,105307,106307
egen,10,1000,e_num+1,e_num+1
!E2钢绞线，97～115
*get,e_num,elem,,num,max
e,97305,98305
e,97300,98300
e,97294,98294
egen,18,1000,e_num+1,e_num+3
!E3钢绞线，90～115
*get,e_num,elem,,num,max
e,90306,91306
e,90303,91303
e,90301,91301
e,90299,91299
egen,25,1000,e_num+1,e_num+4
!E4钢绞线，83～115
*get,e_num,elem,,num,max
e,83304,84304
e,83302,84302
e,83298,84298
e,83296,84296
egen,32,1000,e_num+1,e_num+4
!E5钢绞线，76～115
*get,e_num,elem,,num,max
e,76297,77297
e,76295,77295
egen,39,1000,e_num+1,e_num+2
!!!!!!A类别束曲线束等效为直线束。。。。!!!!
!A1钢绞线，54～76
*get,e_num,elem,,num,max
e,54187,55187
egen,22,1000,e_num+1,e_num+1
!A2钢绞线，54～76
*get,e_num,elem,,num,max
e,54186,55186
egen,22,1000,e_num+1,e_num+1
!A3钢绞线，47～83
*get,e_num,elem,,num,max
e,47146,48146
egen,36,1000,e_num+1,e_num+1
!A4钢绞线，47～83
*get,e_num,elem,,num,max
e,47145,48145
egen,36,1000,e_num+1,e_num+1
!A5钢绞线，40～90
*get,e_num,elem,,num,max
e,40105,41105
egen,50,1000,e_num+1,e_num+1
!A6钢绞线，33～97
*get,e_num,elem,,num,max
e,33104,34104
egen,64,1000,e_num+1,e_num+1
ALLSEL
save
ESYM, ,1000000,all
NUMMRG,NODE, , , ,LOW
csdele,100
SAVE
!
!3月20
!
!
!第三步    3m
!第四步    3.5m
!第五步    4.0m
!第六步    4.5m
!第七～九步    5.5m～7.5m
!第十步    8.0m
!第十一步    8.2m
!第十二步       8.4m
!第十三步       8.6m
!第十四步       8.8m
!第十五步       9.0m
!第十六步       9.2m
!第十七步       9.4m
!第十八步       9.6m
!
FINISH
/SOL
CNVTOL,F,8000000000,0.001,2, ,
NSUBST,1,1,1
!
!
NSEL,S,LOC,Y,500
NSEL,R,LOC,Z,-2400
D,all, , , , , ,UZ, , , , ,
ALLSEL,ALL
!
NSEL,S,LOC,Y,32500
NSEL,R,LOC,Z,-2400
D,all, , , , , ,UY,UZ, , , ,
ALLSEL,ALL
!
NSEL,S,LOC,Y,82500
NSEL,R,LOC,Z,-2400
D,all, , , , , ,UZ, , , , ,
ALLSEL,ALL
!
NSEL,S,LOC,Y,114500
NSEL,R,LOC,Z,-2400
D,all, , , , , ,UZ, , , , ,
ALLSEL,ALL
!
NSEL,S,LOC,X,0
D,all, , , , , ,UX, , , , ,
!第1
ACEL,0,0,0.01,
ALLSEL
!第1荷载步  预应力，控制应力为1000
TREF,0.0
TUNIF,0.0
ESEL,s,TYPE,,4
BFE,all,TEMP,1,-1000/(196000*0.00001), , ,
ALLSEL
!第1荷载步  二期恒载 6.75T/m
!67500/2/82
nsel,s,loc,z,0
F,all,FZ,-67500/164
ALLSEL
lswrite,1
lsread,1
!
!
!加载范围为10.6～24.6 （10.5＋0.5－－14.5＋0.5）（15.5＋0.5－－19.5＋0.5）（20.5＋0.5－－24.5＋0.5
!加载范围为32.5＋10.25＝42.75～56.75 （43－－47）（48－－52）（53－－57
!加载范围为32.5＋10.25＝42.75～56.75 （58－－62）（63－－67）（68－－72
!钢锭7.85T/M3=785KN/M=785000N/M
!每10cm钢锭重量为：!785000*4*4*0.1
!第二荷载步开始10cm
!结点数量为25×9＝225
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-4*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-4*78500*4*4*0.1/351
ALLSEL
lswrite,2
lsread,2
*do,i,3,4
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-2*i*78500*4*4*0.1/351
ALLSEL
lswrite,i
lsread,i
*enddo
*do,i,5,7
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-(i+4)*78500*4*4*0.1/351
ALLSEL
lswrite,i
lsread,i
*enddo
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-13*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-13*78500*4*4*0.1/351
ALLSEL
lswrite,8
lsread,8
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-15*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-15*78500*4*4*0.1/351
ALLSEL
lswrite,9
lsread,9
*do,i,10,17
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,11000,15000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/225
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,16000,20000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,21000,25000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,43000,47000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,48000,52000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,57000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,58000,62000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,53000,67000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
NSEL,S,LOC,Z,0
NSEL,R,LOC,X,-4000,0
NSEL,R,LOC,y,68000,72000
F,all,FZ,-67500/164-(16+0.5*(i-10))*78500*4*4*0.1/351
ALLSEL
lswrite,i
lsread,i
*enddo
lssolve,1,17
/REPLOT,RESIZE
WPSTYLE,,,,,,,,0
/USER,  1
/VIEW,  1, -0.107982741112E-01, -0.278950934295 ,  0.960244642543
/ANG, 1,  0.232639069558
/REPLO
FINISH
/POST1
FINISH
/POST26
FINISH
/CLEAR
/COM,ANSYS RELEASE  8.1 UP20040329    08:52:01 04/21/2005
/input,start81,ans,'d:\Program Files\Ansys Inc\v81\ANSYS\apdl\',,,,,,,,,,,,,,,,1
/POST1
INRES,BASIC
FILE,'xxtbridge1','rst','.'
SET,LAST
FINISH
/POST26
FILE,'xxtbridge1','rst','.'
/UI,COLL,1
NUMVAR,200
SOLU,191,NCMIT
STORE,MERGE
FILLDATA,191,,,,1,1
REALVAR,191,191
!*
!*
!*
NSOL,2,115041,U,Z,
!*
!*
!*
NSOL,3,33041,U,Z,
!*
!*
!*
NSOL,4,33023,U,Z,
!*
!*
!*
NSOL,5,115023,U,Z,
!*
/USER,  1
/REPLO
NUMVAR,200
FILLDATA,191,,,,1,1
REALVAR,191,191
PRVAR,5,
PRVAR,2,
PRVAR,5,
PRVAR,3,
PRVAR,4,
NSEL,S,LOC,X,98500
/REPLOT
NSEL,S,LOC,X,98500
NLIST,ALL, , , ,NODE,NODE,NODE
NSEL,S,LOC,Y,98500
/REPLOT
NLIST,ALL, , , ,NODE,NODE,NODE
!*
!*
!*
NSOL,6,1033041,U,Z,
!*
!*
!*
NSOL,7,1033023,U,Z,
!*
NUMVAR,200
FILLDATA,191,,,,1,1
REALVAR,191,191
PRVAR,6,
PRVAR,7,
ALLSEL,ALL
ALLSEL,ALL
/REPLOT
FINISH
/POST1
SET,LIST
SET,LIST,999
SET,,, ,,, ,1
PRRSOL,FZ
FLST,2,24,1
FITEM,2,1023
FITEM,2,11023
FITEM,2,21023
FITEM,2,31023
FITEM,2,41023
FITEM,2,51023
FITEM,2,61023
FITEM,2,71023
FITEM,2,81023
FITEM,2,91023
FITEM,2,101023
FITEM,2,111023
FITEM,2,1109023
FITEM,2,1099023
FITEM,2,1089023
FITEM,2,1079023
FITEM,2,1069023
FITEM,2,1059023
FITEM,2,1049023
FITEM,2,1039023
FITEM,2,1029023
FITEM,2,1019023
FITEM,2,1009023
FITEM,2,1000023
!*
PATH,wjq22,24,30,1,
PPATH,P51X,1
PATH,STAT
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,YG,UZ
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,1.08222638492,1
/REP,FAST
SET,LIST,999
SET,,, ,,, ,17
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
PRRSOL,FZ
SET,LIST,999
SET,,, ,,, ,16
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,2
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,15
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,3
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,14
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,4
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,13
PRRSOL,FZ
/DIST,1,0.924021086472,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,5
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,12
PRRSOL,FZ
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,6
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,11
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,7
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,10
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,8
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
SET,LIST,999
SET,,, ,,, ,9
PRRSOL,FZ
AVPRIN,0, ,
!*
PDEF, ,U,Z,AVG
/PBC,PATH, ,0
!*
PRPATH,UZ
FINISH
! /EXIT,NOSAV

显示全部楼层 · 发表于 2005-12-6 20:29:35

提示: 作者被禁止或删除内容自动屏蔽

阿波罗 · 发表于 2005-12-7 18:04:38

本算例特色：
(1) 荷载是在APDL中用程序计算的；
(2) 用Link10单元模拟地层抗力作用，Link10的弹性模量的确定比较有技巧；
(3) 考虑了地震荷载，用静力法。
                                    浙江省交通规划设计研究院郑国平
                                    tjzhenggp@sohu.com

fini
/clear
/filename,MingDong
/triad,off
/prep7
*SET,gamma1,22      !拱圈填土的重度
*SET,gamma2,22      !边墙填土的重度
*SET,phi1,40      !拱圈填土的计算内摩擦角
*SET,phi2,40      !边墙填土的计算内摩擦角
*SET,mu,0.4      !回填土石与开挖边坡面间的摩擦系数
*SET,PI,3.1415926
*SET,up_shirtsegment_len,26.85/27  !拱圈,边墙的分段长度
*SET,down_shirtsegment_len,26.85/27  !仰拱的分段长度
*SET,Left_k,300000      !左侧边墙的浆砌片石的弹性抗力(单位：kPa)
*SET,Mid_k,500000      !中间边墙的浆砌片石的弹性抗力(单位：kPa)

*SET,link_area,0.01
*SET,link_len,1.0
*SET,soilline_k,0      !填土线的斜率 y=kx+b
*SET,soilline_b,19.5      !填土线的截距
*SET,wall_height,5.34
*Set,alpha,atan(soilline_k)
*SET,rho,alpha      !侧压力作用方向与水平线夹角
*SET,concrete_gamma,25    !仰拱填充材料的重度
*SET,roadsurf_height,2.25  !路面离计算基点的高度

ET,1,LINK10
KEYOPT,1,2,0
KEYOPT,1,3,1
ET,2,BEAM3
R,1,link_area, ,  !拱圈链杆
R,2,0.85,0.0512,0.85, , , ,

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,Left_k*up_shirtsegment_len*1*link_len/link_area     !拱部链杆的折算弹性模量
MPDATA,PRXY,1,,0.2

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,Mid_k*up_shirtsegment_len*1*link_len/link_area     !拱部链杆的折算弹性模量
MPDATA,PRXY,2,,0.2

MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,29500000           !C25混凝土
MPDATA,PRXY,3,,0.2
MPDATA,DENS,3,,2500

K,1,0,0
K,2,-.99782359912183,1.79778906053798E-02
K,3,-1.99435198066192,7.18882263520584E-02
K,4,-2.98829160828598,.16166102932533
K,5,-3.97835230597286,.287179770598769
K,6,-4.96324893271787,.448281521493712
K,7,-5.94170305070065,.644757165067631
K,8,-6.91244458475171,.876351667556889
K,9,-7.87421347096472,1.14276440942135
K,10,-8.81337745272924,1.47574415383816
K,11,-9.60299972317716,2.07747331407765
K,12,-10.1330205166551,2.91691316942824
K,13,-10.362486416782,3.88539361509889
K,14,-10.4299999901402,4.88038074860016
K,15,-10.3623834715027,5.87536089153166
K,16,-10.1608783046729,6.85206614193302
K,17,-9.82918413798916,7.7925641260928
K,18,-9.37339089853265,8.679587238129
K,19,-8.8018669806595,9.49684967438779
K,20,-8.12586539875008,10.2301299330995
K,21,-7.37904001648599,10.8917423991212
K,22,-6.57856844012281,11.4873298924619
K,23,-5.73027016019853,12.0125624457816
K,24,-4.84031237110128,12.4636215776195
K,25,-3.91516513508668,12.8372280530502
K,26,-2.96155434442289,13.130665723977
K,27,-1.98641282363144,13.3418012757388
K,28,-.996829927314252,13.4690997364742
K,29,5.11590769747272E-12,13.5116356364854
K,30,.996829927324484,13.4690997364734
K,31,1.98641282364156,13.3418012757373
K,32,2.9615543444329,13.1306657239747
K,33,3.91516513509646,12.8372280530472
K,34,4.84031237111071,12.4636215776158
K,35,5.73027016020762,12.0125624457772
K,36,6.57856844013145,11.4873298924569
K,37,7.37904001649417,10.8917423991156
K,38,8.12586539875781,10.2301299330933
K,39,8.80186698066666,9.49684967438104
K,40,9.37339089853913,8.67958723812181
K,41,9.82918413799496,7.79256412608525
K,42,10.160878304678,6.85206614192522
K,43,10.3623834715071,5.87536089152371
K,44,10.4299999901439,4.88038074859216
K,45,10.3624864167848,3.88539361509095
K,46,10.1330205166572,2.91691316942047
K,47,9.60299972317864,2.07747331407029
K,48,8.81337745273026,1.4757441538314
K,49,7.87421347096551,1.14276440941531
K,50,6.91244458475228,.876351667551589
K,51,5.94170305070099,.644757165063069
K,52,4.9632489327181,.448281521489903
K,53,3.97835230597298,.287179770595714
K,54,2.98829160828598,.161661029323042
K,55,1.9943519806618,7.18882263505378E-02
K,56,.997823599121716,1.79778906046124E-02
K,57,-4.54747350886464E-13,-1
K,58,-1.03384611172589,-.981373088073127
K,59,-2.06635024711227,-.925516530819365
K,60,-3.09617217176242,-.832502832431899
K,61,-4.12197513290369,-.702452728657619
K,62,-5.1424275945493,-.535535030076943
K,63,-6.15620496588883,-.331966402980484
K,64,-7.16199132066379,-9.20110881270233E-02
K,65,-8.158481105296,.184019442252406
K,66,-9.26888767884464,.585513602854746
K,67,-10.3638868367301,1.4279026717532
K,68,-11.0598318370295,2.5413858382985
K,69,-11.3533133755853,3.75025687801995
K,70,-11.4299999888021,4.88043248074692
K,71,-11.3531964431756,6.01060014309844
K,72,-11.1243128569398,7.12000991121739
K,73,-10.7475515480436,8.18829294757279
K,74,-10.2298298793131,9.19583550417914
K,75,-9.58065325520727,10.124139032321
K,76,-8.82028673019408,10.9496985639023
K,77,-8.0096390409069,11.667851199741
K,78,-7.14076065363031,12.3143364916324
K,79,-6.21996838173118,12.8844544412878
K,80,-5.25395645655033,13.3740602469284
K,81,-4.24974785982056,13.7795944363161
K,82,-3.21464326604337,14.0981087443937
K,83,-2.15616796601932,14.3272875474001
K,84,-1.08201715739835,14.4654646976358
K,85,4.54747350886464E-12,14.5116356364854
K,86,1.08201715740938,14.465464697635
K,87,2.15616796603024,14.3272875473985
K,88,3.21464326605405,14.0981087443912
K,89,4.24974785983102,13.7795944363129
K,90,5.25395645656045,13.3740602469244
K,91,6.21996838174096,12.8844544412831
K,92,7.14076065363963,12.3143364916269
K,93,8.00963904091577,11.6678511997348
K,94,8.82028673020238,10.9496985638956
K,95,9.58065325521488,10.1241390323136
K,96,10.2298298793201,9.1958355041713
K,97,10.7475515480497,8.18829294756455
K,98,11.1243128569452,7.12000991120885
K,99,11.3531964431801,6.01060014308973
K,100,11.4299999888058,4.88043248073815
K,101,11.3533133755881,3.75025687801124
K,102,11.0598318370313,2.54138583829001
K,103,10.3638868367311,1.42790267174524
K,104,9.26888767884498,.585513602847641
K,105,8.158481105296,.184019442246154
K,106,7.16199132066356,-9.20110881325229E-02
K,107,6.15620496588849,-.331966402985202
K,108,5.14242759454874,-.535535030080879
K,109,4.121975132903,-.702452728660774
K,110,3.09617217176162,-.832502832434272
K,111,2.06635024711147,-.925516530820943
K,112,1.03384611172498,-.981373088073909

*do,i,1,55
  l,i,i+1
*enddo

!以下建立链杆
*do,i,1,56
  l,i,i+56
*enddo
l,56,1

!建立选择集
ksel,s,kp,,1,56
cm,kWall,KP    
lslk,s,1
cm,lShirt,LINE

allsel
ksel,s,kp,,57,112
cm,kLink,kp    !弹性链杆的关键点

lslk,s,0
cm,lLink,LINE  !左侧边墙的弹性链杆

allsel
!赋材料
cmsel,s,lshirt,LINE
latt,3,2,2

allsel
ksel,s,kp,,85,103
lslk,s,0
latt,2,1,1

ksel,s,kp,,57,112
ksel,u,kp,,85,103
lslk,s,0
latt,1,1,1

allsel
lsel,all
LESIZE,ALL, , ,1, , , , ,1
LMESH,ALL
/PBC,ALL, ,1
/PBF,DEFA, ,1
/PIC,DEFA, ,1

/solu
allsel
*AFUN,DEG    !设置三角计算的角度和弧度模式:角度
nsel,all
DOFSEL,all  !设置荷载的添加模式为叠加,即原先有荷载的话,则增加新增的荷载
FCUM,ADD,1,
      !添加桩底位移约束
!acel,,9.8
cmsel,s,klink,kp
DK,all, , , ,0,ALL, , , , , ,

allsel
ksel,s,kp,,1,56
DK,all, , , ,0,UZ, , , , , ,

!施加土压力荷载
*Set,total_weight,0
allsel
*do,i,13,44
  *get,this_x,KP,i,LOC,x
  *get,this_y,KP,i,LOC,y
  *get,this_z,KP,i,LOC,z
    *get,prev_x,KP,i-1,LOC,x
    *get,prev_y,KP,i-1,LOC,y
    *get,prev_z,KP,i-1,LOC,z
    *get,next_x,KP,i+1,LOC,x
    *get,next_y,KP,i+1,LOC,y
    *get,next_z,KP,i+1,LOC,z
  
  dDepth = soilline_k*this_x+soilline_b-this_y
  *if,dDepth,le,0,then
    dDepth=0
  *endif
  *if,this_y,gt,wall_height,then    !拱圈的土压力
    lambda=cos(alpha)*(cos(alpha)-sqrt(cos(alpha)**2-cos(phi1)**2))/(cos(alpha)+sqrt(cos(alpha)**2-cos(phi1)**2))
    nodefy=gamma1*dDepth*(next_x-prev_x)*0.5
    nodefx=nodefy*lambda*(next_y-prev_y)*0.5
  *else          !边墙的土压力
    *if,soilline_k,gt,0,then  !填土斜面向上倾斜
      alpha_p=atan(gamma1*tan(alpha)/gamma2)
      *if,this_x,gt,0,then
        aa=cos(phi2)**2
        bb=1+sqrt(sin(phi2)*sin(phi2-alpha_p)/cos(alpha_p))
        lambda=aa/(bb**2)
      *else
        aa=-tan(phi2)+sqrt((1+tan(phi2)**2)*(1+tan(alpha_p)/tan(phi2)))
        bb=1+(1+tan(phi2)**2)*tan(alpha_p)/tan(phi2)
        tan_theta0=aa/bb
        theta0=atan(tan_theta0)
        lambda=tan_theta0/(tan(theta0+phi2)*(1+tan(alpha_p)*tan(theta0)))
      *endif
    *endif
    *if,soilline_k,lt,0,then
      alpha_p=atan(gamma1*tan(alpha)/gamma2)
      *if,this_x,lt,0,then
        aa=cos(phi2)**2
        bb=1+sqrt(sin(phi2)*sin(phi2-alpha_p)/cos(alpha_p))
        lambda=aa/(bb**2)
      *else
        aa=-tan(phi2)+sqrt((1+tan(phi2)**2)*(1+tan(alpha_p)/tan(phi2)))
        bb=1+(1+tan(phi2)**2)*tan(alpha_p)/tan(phi2)
        tan_theta0=aa/bb
        theta0=atan(tan_theta0)
        lambda=tan_theta0/(tan(theta0+phi2)*(1+tan(alpha_p)*tan(theta0)))
      *endif
    *endif
    *if,soilline_k,eq,0,then
      lambda=(tan(45-phi2*0.5))**2
    *endif
    nodefx=gamma2*(wall_height-this_y+gamma1*(dDepth-wall_height+this_y)/gamma2)*lambda*(next_y-prev_y)*0.5
    nodefy=0
  *endif
  nodefy=-nodefy
  fk,i,FX,nodefx + 0.85*25*1.7*0.2*0.1    !包含了地震荷载
  fk,i,FY,nodefy
*enddo

allsel
time,1
AUTOTS,1
NSUBST,20, , ,1
KBC,1
!输出每个子步的结果
OUTPR,BASIC,ALL,
OUTRES,ALL,ALL,
!设定大变形,Newton-Raphson选项为full N-R,应力刚化
NLGEOM,ON
NROPT,FULL, ,OFF
SSTIF,ON
/STATUS,SOLU
allsel
SOLVE
FINISH

/post1
!!!!!!!!!!!!!!Beam3梁的单元表!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
etable,beam_mz1,SMISC, 6  !I节点绕z轴的弯矩
etable,beam_mz2,SMISC, 12   !J节点绕z轴的弯矩
etable,beam_fx1,SMISC, 1  !轴力MFORX
etable,beam_fx2,SMISC, 7
etable,beam_fy1,SMISC, 2  !剪力MFORY
etable,beam_fy2,SMISC, 8
etable,beam_smax1,NMISC, 1  !大应力
etable,beam_smax2,NMISC, 3
etable,beam_smin1,NMISC, 2  !小应力
etable,beam_smin2,NMISC, 4
etable,beam_sdir1,LS, 1    !中心应力
etable,beam_sdir2,LS, 4
ETABLE,Beam_uy,U,Y
ETABLE,Beam_ux,U,X
!!!!!!!!!!!!!!LINK10的单元表!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
etable,link_saxl,LS,1
etable,link_mforx,SMISC,1
cmsel,s,lshirt,line
esll,s
plls,beam_MZ1,beam_MZ2,1,0
finish

账号		自动登录	找回密码
密码			注册

账号		自动登录	找回密码
密码			立即注册

yudou007 该用户已被删除	发表于 2005-10-7 17:06:27 \| 显示全部楼层来自甘肃兰州提示: 作者被禁止或删除内容自动屏蔽
yudou007 该用户已被删除
	回复赞不支持使用道具举报显身卡

lsh0016 该用户已被删除	发表于 2005-10-16 10:25:33 \| 显示全部楼层来自上海大学提示: 作者被禁止或删除内容自动屏蔽
lsh0016 该用户已被删除
	回复赞不支持使用道具举报显身卡

chxyz_china 该用户已被删除	发表于 2005-11-29 19:15:06 \| 显示全部楼层来自江苏南京提示: 作者被禁止或删除内容自动屏蔽
chxyz_china 该用户已被删除
	回复赞不支持使用道具举报显身卡

[共建]＜ANSYS命令流实例库＞敬请您的加盟

结构瞬态动力学分析（原创）

本帖子中包含更多资源

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

本帖子中包含更多资源

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

结构

本帖子中包含更多资源

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

楼主请大家提意见!

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

本帖子中包含更多资源

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

评分

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

评分

隧道衬砌荷载结构法算例

评分

浏览过的版块

zoudongping 该用户已被删除	发表于 2005-12-6 20:29:35 \| 显示全部楼层来自广东深圳提示: 作者被禁止或删除内容自动屏蔽
zoudongping 该用户已被删除
	回复赞不支持使用道具举报显身卡

[共建]＜ANSYS命令流实例库＞敬请您的加盟

结构瞬态动力学分析（原创）

本帖子中包含更多资源

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

本帖子中包含更多资源

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

结构

本帖子中包含更多资源

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

楼主请大家提意见!

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

本帖子中包含更多资源

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

评分

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

Re:[共建] &lt;&lt;ANSYS命令流实例库&gt;&gt; 敬请您的加盟

评分

隧道衬砌荷载结构法算例

评分

浏览过的版块

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟

Re:[共建] <<ANSYS命令流实例库>> 敬请您的加盟