VibInfo 发表于 2005-11-17 08:12

三杆桁架优化分析APDL命令

三杆桁架优化分析APDL命令

命令流如下:
/filnam,truss
/title, optimization of a three-bar truss
!初始化设计变量参数
B=1000 !基本尺寸
A1=1000 !第一个面积
A2=1000 !第二个面积
A3=1000 !第三个面积

!进入PREP7并建模
/prep
et,1,link1 !二维单元
r,1,A1 !以参数形式的实参
r,2,A2
r,3,A3
mp,ex,1,2.1E6 !杨氏模量
n,1,-B,0,0
n,2,0,0,0
n,3,B,0,0
n,4,0,-1000,0
e,1,4
real,2
e,2,4
real,3
e,3,4
finish
!
!进入求解器,定义载荷和求解
/solu
d,1,all,0,,3
f,4,fx,200000
f,4,fy,-20000
solve
finish
!
!进入POST1并读出状态变量数值
/post1
set,last
etable,evol,volu !将每个单元的体积放入ETABLE
ssum !将单元表格内数据求和
*get,vtot,ssum,,item,evol !VTOT=总体积
rho=2.85e-4
wt=tho*vtot !计算总质量
etable,sig,ls,1 !将轴向应力放入ETABLE
!
*get,sig1,elem,1,etab,sig !SIG1=第一个单元的轴向应力
*get,sig2,elem,2,etab,sig !SIG2=二单元的轴向应力
*get,sig3,elem,3,etab,sig !SIG3=三单元的轴向应力

sig1=abs(sig1) !计算轴向应力的绝值
sig2=abs(sig2)
sig3=abs(sig3)
!
/eshape,2 !以实体单元模式显示壳单元
/view,1,1,1,1 !轴测视图
eplot
!
/opt !进入优化处理器
opanl,truss,lgw !指定分析文件(批处理方式中不用这个命令)

opvar,B,dv,400,2000 !定义设计变量
opvar,A1,dv,1,1000
opvar,A2,dv,1,1000
opvar,A3,dv,1,1000
opvar,sig1,sv,,400 !定义状态变量
opvar,sig2,sv,,400
opvar,sig3,sv,,400

opsave,trussvar,opt !存储数据
!
opvar,wt,obj,,,2, !定义目标函数
!
optype,first !定义一阶方法
opfrst,45 !最大45次迭代
opexs !开始优化分析
!
oplist,16 !列出最佳设计序列,号为16
oplist,all
!
/view,1,,,1 !前视图
!
/axlab,x,iteration number !画重量对迭代数图形
/axlab,y,structure weight
plvaropt,wt
!
/axlab,y,base dimension !画B对迭代数图形
plvaropt,B

/axlab,y,max stress !画最大应力对迭代图形
plvaropt,sig1,sig2,sig3
!
/axlab,y,cross-sectional area !画面积对迭代图形
plavaropt,A1,A2,A3
!
finish
/exit

[ 本帖最后由 rodge 于 2008-3-24 03:36 编辑 ]

dlt2000 发表于 2008-3-23 19:02

回复 楼主 的帖子

错误
1.rho=2.85e-4
wt=tho*vtot !计算总体积
2.*get,sig,elem,1,etab,sig !SIG1=第一个单元的轴向应力
*get,sig,elem,2,etab,sig !SIG2=二单元的轴向应力
*get,sig,elem,3,etab,sig !SIG3=三单元的轴向应力

sig1=abs(sig1) !计算轴向应力的绝值
sig2=abs(sig2)
sig3=abs(sig3)

主题贴已经编辑过

[ 本帖最后由 rodge 于 2008-3-24 03:36 编辑 ]
页: [1]
查看完整版本: 三杆桁架优化分析APDL命令