马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
命令流如下,我看优化的结果,设计变量改变了,但是目标函数和状态变量没有改变。不知道是不是跟我的施工阶段有关系呢?
/prep7
/title,the anslyse of small arch bridge
!**************************************************************
!**************************************************************
!以下为建立各个截面的命令
!**************************************************************
et,1,82
!********截面9-9****************
!为什么只有一个82。作用何在?
rectng,-0.1,0.1,-0.1,0.1 !1, 与图纸编号不同,先做框框
esize,0.1 !自动划分网格,5为划分的细度,第一个数字为 1-10
amesh,all
secwrite,xxx11,sect,,1 !建立一个包含单元划分信息的ASC_2码文件,第一个为文件名,和路径,路径省略则为文件默认路径
SECTYPE,1,BEAM,MESH, !截面代号1
SECOFFSET,CENT,,,
SECREAD,'xxx11','SECT',,MESH
asel,all
ACLEAR,ALL
ADELE,ALL,,,1
!*************9-9**********************
etdele,1 !删除代理单元
!定义钢材特性
et,1,beam188 !不知道此处的截面和上面的截面有什么区别?
mp,ex,1,2.1e11
mp,dens,1,7800
mp,prxy,1,0.3
mp,alpx,1,1e-5
et,2,link10
mp,ex,2,2.1e11
mp,dens,2,7800
mp,prxy,2,0.3
mp,alpx,2,1e-5
r,2,0.1,0.0001
!************拱肋计算************
!**********计算拱肋中心线坐标**********
*dim,cs,array,48
*dim,y,array,48
*dim,z,array,48
m=1.9
k=1.257195827
f=8.5
y(1)=0
*do,i,1,46,1
y(i+1)=y(i)+0.5
*enddo
y(48)=y(47)+0.349166446
*do,i,1,48,1
cs(i)=(23.34916645-y(i))*2/23.34916645/2
z(i)=f-f*(cosh(k*cs(i))-1)/(m-1)
*enddo
*do,i,1,47,1
n,i,y(i)-23.34916645,z(i)-8.5
*enddo
n,48,y(48)-23.34916645,z(48)-8.5
*do,i,47,1,-1
n,96-i,23.34916645-y(i),z(i)-8.5
*enddo
!********************************************type,1
real,1
mat,1
secnum,1
*dim,kzd,array,5
kzd(1)=16
kzd(2)=32
kzd(3)=48
kzd(4)=64
kzd(5)=80
*dim,yyc,array,5
*do,i,1,5
*get,yyc(i),node,kzd(i),loc,y
*enddo
*do,i,1,20,1
e,i,i+1
*enddo
n,150,-24,5
type,2
mat,2
real,2
e,150,18
cm,jd1,elem
esel,none
type,1
real,1
mat,1
*do,i,21,40,1
e,i,i+1
*enddo
type,2
mat,2
real,2
n,151,-24,6.5
e,151,38
cm,jd2,elem
esel,none
type,1
real,1
mat,1
*do,i,95,76,-1
e,i,i-1
*enddo
type,2
mat,2
real,2
n,153,24,5
e,153,78
cm,jd3,elem
esel,none
type,1
real,1
mat,1
*do,i,75,56,-1
e,i,i-1
*enddo
type,2
mat,2
real,2
n,154,24,6.5
e,154,58
cm,jd4,elem
esel,none
type,1
real,1
mat,1
*do,i,41,54
e,i,i+1
*enddo
cm,helong,elem
esel,none
allsel
/solu
h1=200
h2=200
bfe,21,temp,,-538.02
bfe,42,temp,,-618.30
bfe,63,temp,,-618.3
bfe,84,temp,,-538.02
acel,,10
time,1
nlgeom,on
nropt,full
estif,1e-10
allsel
ekill,all
esel,s,live
esel,all
ealive,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
nsel,s,,,150,154
d,all,all
nsel,s,,,1
nsel,a,,,95
d,all,all
allsel
solve
save
time,2
allsel
ddele,all,all
ekill,all
esel,s,live
esel,s,mat,,1
ealive,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
nsel,s,,,1
nsel,a,,,95
d,all,all
allsel
solve
save
!************优化成桥索力过程
/post1
set,last
!esort,etab,uy,0,0
*get,uyy1,node,16,u,y
*get,uyy2,node,32,u,y
*get,uyy3,node,48,u,y
*get,uyy4,node,64,u,y
*get,uyy5,node,80,u,y
uyyma=max(uyy1,uyy2,uyy3,uyy4,uyy5)
uyymi=min(uyy1,uyy2,uyy3,uyy4,uyy5)
uyyx=max(uyyma,-uyymi)
esel,s,mat,,1
etable,wanjs,smisc,32
etable,wanjx,smisc,33
plls,wanjs,wanjs
plls,wanjx,wanjx
!*dim,xigms,array,94
!*dim,xigmx,array,94
esort,etab,wanjs,0,0
*get,sminT,sort,,max
esort,etab,wanjx,0,0
PRETAB,wanjx
*get,sminB,sort,,max
smin=mAX(sminT,sminB)
!/eof
lgwrite,arch,lgw
/input,arch.lgw
/opt !进入优化处理器
OPANL,arch,lgw
opvar,h1,dv,0,1000
opvar,h2,dv,0,1000 !设计变量
!opvar,h3,dv,0,1000
!opvar,h4,dv,0,1000
opvar,smin,sv,-2e7,2e7
OPSAVE,trussvar,opt
opvar,uyyx,obj,,,1e-3 !目标函数
optype,subp
opsubp,40
opexe
/opt
oplist,all,,1
finish |