ansys的问题

发布网友 发布时间：2022-04-26 17:41

我来回答

共3个回答

热心网友 时间：2023-10-18 08:38

*CREATE,SCRATCH

/FILNAME,ZUOYE
/TITLE,NO6

A1=0.002
A2=0.002
disp=0.003

/PREP7
mp,prxy,1,0.3
ET,1,LINK1
R,1,A1
R,2,A2
MP,EX,1,2E11Y

N,1,1.44,0.72
N,2,1.44
N,3,0.72,0.72
N,4,0.72
N,5,,0.72
N,6

REAL,1
E,1,2
E,1,3
E,6,4
E,2,4
E,3,4
E,3,5
REAL,2
E,1,4
E,3,2
E,4,5
E,3,6

FINISH
/SOLU
D,5,UX
D,5,UY
D,6,UX
D,6,UY

F,2,FX,6E5
F,2,FY,-9E5
F,4,FY,-6E5
SOLVE

FINISH
/POST1
SET,LAST
NSORT,U,Y,,1
*GET,DISP,SORT,,MAX
disp=abs(disp)

ETABLE,EVOL,VOLU
SSUM
*GET,VTOT,SSUM,,ITEM,EVOL
RHO=7.8E3
WT=RHO*VTOT

ETABLE,SIG,LS,1
*GET,SIG1,ELEM,1,ETAB,SIG
*GET,SIG2,ELEM,2,ETAB,SIG
*GET,SIG3,ELEM,3,ETAB,SIG
*GET,SIG4,ELEM,4,ETAB,SIG
*GET,SIG5,ELEM,5,ETAB,SIG
*GET,SIG6,ELEM,6,ETAB,SIG
*GET,SIG7,ELEM,7,ETAB,SIG
*GET,SIG8,ELEM,8,ETAB,SIG
*GET,SIG9,ELEM,9,ETAB,SIG
*GET,SIG10,ELEM,10,ETAB,SIG
SIG1=ABS(SIG1)
SIG2=ABS(SIG2)
SIG3=ABS(SIG3)
SIG4=ABS(SIG4)
SIG5=ABS(SIG5)
SIG6=ABS(SIG6)
SIG7=ABS(SIG7)
SIG8=ABS(SIG8)
SIG9=ABS(SIG9)
SIG10=ABS(SIG10)

*END
*USE,SCRATCH

finish
/OPT
OPANL,SCRATCH

OPVAR,A2,DV,0.003,0.01
OPVAR,A1,DV,0.003,0.01

OPVAR,disp,SV,,0.005
OPVAR,SIG1,SV,,125E6
OPVAR,SIG2,SV,,125E6
OPVAR,SIG3,SV,,125E6
OPVAR,SIG4,SV,,125E6
OPVAR,SIG5,SV,,125E6
OPVAR,SIG6,SV,,125E6
OPVAR,SIG7,SV,,125E6
OPVAR,SIG8,SV,,125E6
OPVAR,SIG9,SV,,125E6
OPVAR,SIG10,SV,,125E6
OPVAR,WT,OBJ,,,2 !定义目标函数

OPSAVE,trussvar,opt !存储数据

OPTYPE,first !定义1阶方法
OPFRST,30 !最大30次迭代
OPEXE

STATUS
OPLIST,ALL

/AXLAB,X,ITERATION NUMBER
/AXLAB,Y,STRUCTURE WEIGHT
PLVAROPT,WT
/AXLAB,Y,DEFELECTION
/AXLAB,Y,CROSS-SECTIONAL AREA
PLVAROPT,A1,A2

这个是我改后的程序，是可以运行的！
运行方法是清除掉你的工作目录里面的所有不必要的文件，因为宏会冲突，你的优化数据库也会冲突，也就是说保证你的工作目录没有跟这个程序有关的任何东西！这样就不会出现警告！

但是我建议你在做优化的时候一定要主要你的优化变量的上下限以及公差！我觉得你的DV和SV都设置的不怎么样，这样的优化结果是不正确的：
SET 1 SET 2 *SET 3* SET 4
(INFEASIBLE) (INFEASIBLE) (INFEASIBLE) (INFEASIBLE)
DISP (SV) > 0.13957E-01 0.28094E-02 0.27914E-02 0.27914E-02
SIG1 (SV) > 0.15970E+09 0.32227E+08 0.31941E+08 0.31941E+08
SIG2 (SV) > 0.15970E+09 0.32227E+08 0.31941E+08 0.31941E+08
SIG3 (SV) > 0.56944E+09 0.11526E+09 0.11389E+09 0.11389E+09
SIG4 (SV) 0.97033E+07 0.18801E+07 0.19407E+07 0.19407E+07
SIG5 (SV) 0.40264E+08 0.80049E+07 0.80528E+07 0.80528E+07
SIG6 (SV) > 0.78056E+09 > 0.15786E+09 > 0.15611E+09 > 0.15611E+09
SIG7 (SV) > 0.22585E+09 0.45055E+08 0.45171E+08 0.45171E+08
SIG8 (SV) > 0.41054E+09 0.82224E+08 0.82108E+08 0.82108E+08
SIG9 (SV) > 0.59318E+09 0.11872E+09 0.11864E+09 0.11864E+09
SIG10 (SV) > 0.46748E+09 0.93416E+08 0.93497E+08 0.93497E+08
A2 (DV) > 0.20000E-02 0.10000E-01 0.10000E-01 0.10000E-01
A1 (DV) > 0.20000E-02 0.98856E-02 0.10000E-01 0.10000E-01
WT (OBJ) 130.93 650.80 654.65 654.65

你可以看下你的最优解发生在序列3，但是结果并不正确！

热心网友 时间：2023-10-18 08:39

1.估计你约束的命令不对，你要在solution-defind
loads-apply-structure-displacement-on
nodes(或者是其他的)里加
UY,UX.UZ是个复选框，可以多选
2.是你划分的网格超出了它的*
比如用SHELL63单元划分时，遇到小的多边形就会出现这种情况。
你可以用check
mesh命令找出来对应单元来，从而确定出问题的面，然后在模型中将面变为规整的
或者减小网格尺寸
出问题的地方一般是细小的面等
或者是进行布尔操作时留下来的
不过还要说句，这个问题不是特别厉害，只要能往下进行下去并且没有错误的话可以忽略。不过这种面能减少就减少点。

热心网友 时间：2023-10-18 08:39

如果直接定义截面惯性矩和面积，则直接通过定义实常数，而且section并不支持所有的beam单元，beam3和beam4只支持矩形截面，而且也是通过实常数定义，beam44可以定义截面，但是计算参数仍取实常数定义值，beam188/189只支持截面定义，不支持实常数。