求助hy1966老师关于circu94压电电路分析【附命令流】

steven710

想要在一个压电陶瓷圆片的两级施加电压激励，观察压电陶瓷片两端的电压波形，应用circu94单元，出现了关于circu94单元定义节点node的错误！请hy1966老师指教！！！
/PREP7
/VIEW,1,,-1
/TITLE, RADIUS PIEZO TRANSIENT ANALYSIS
ET,1,SOLID5,0
ET,2,MESH200
KEYOPT,2,1,6
KEYOPT,2,2,0
MP,DENS,1,7500
MP,DMPR,1,0.005
MP,PERX,1,804.6
MP,PERY,1,804.6
MP,PERZ,1,659.7
TB,PIEZ,1
TBDATA,16,10.5
TBDATA,14,10.5
TBDATA,3,-4.1
TBDATA,6.-4.1
TBDATA,9,14.1
TB,ANEL,1
TBDATA,1,13.2E10,7.1E10,7.3E10
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA,21,2.6E10
R=15E-3
T=15.5E-3
K,1,0,0,0
K,2,0,0,T
K,3,R,0,0
K,4,R,0,T
A,1,2,4,3
AATT,1
ESIZE,5E-3
AMESH,ALL
EXTOPT,ESIZE,16,0,
VROTAT,ALL,,,,,,1,2,360
EPLOT

NSEL,S,LOC,Z,T
CP,1,VOLT,ALL
*GET,N_TOP,NODE,O,NUM,MIN
D,ALL,VOLT,1
NSEL,S,LOC,Z,0
CP,2,VOLT,ALL
*GET,N_BOT,NODE,0,NUM,MIN
D,ALL,VOLT,0
NSEL,ALL
FINISH
/SOLU
ANTYPE,STATIC
SOLVE
FINI
/PREP7
DDELE,N_TOP,VOLT
ET,3,CIRCU94,0
R=3000
R,1,R
N,1,0,0
TYPE,3
REAL,1
E,,N_TOP
ET,4,CIRCU94,4
V=1
R,2,V
N,2,0,1,0
TYPE,4
REAL,2
E,N_TOP,N_BOT,2
FINI
/SOLU
ANTYPE,TRANS
NSUB,100
TIME,2E-3
TINTP,,0.25,0.5,0.5
OUTERS,ALL,ALL
SOLVE
FINI
/POST26
NUMVAR,20
NSOL,2,N_TOP,VOLT,,V_C
PRVAR,2
PLVAR,2

steven710

我觉得是我在circu94单元内定义节点的问题，在circu94单元之前的命令应该没有什么问题，其实我对下面的这段代码不是很理解，求详细的解释！
ET,3,CIRCU94,0
R=3000
R,1,R
N,1,0,0                   ？这个节点是用来做什么的？如果直接写成N，1同样出错
TYPE,3
REAL,1
E,1,N_TOP               
ET,4,CIRCU94,4
V=1
R,2,V
N,2,0,1,0                 ？这句的问题与上面相同！
TYPE,4
REAL,2
E,N_TOP,N_BOT,2
FINI 1# steven710

steven710

应用了numstr,将起始节点数设为10，更新程序，出现新的问题：
/PREP7
/VIEW,1,,-1
/TITLE, RADIUS PIEZO TRANSIENT ANALYSIS
ET,1,SOLID5,0
ET,2,MESH200
KEYOPT,2,1,6
KEYOPT,2,2,0
MP,DENS,1,7500
MP,DMPR,1,0.005
MP,PERX,1,804.6
MP,PERY,1,804.6
MP,PERZ,1,659.7
TB,PIEZ,1
TBDATA,16,10.5
TBDATA,14,10.5
TBDATA,3,-4.1
TBDATA,6.-4.1
TBDATA,9,14.1
TB,ANEL,1
TBDATA,1,13.2E10,7.1E10,7.3E10
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA,21,2.6E10
R=2E-3
T=1E-3
K,1,0,0,0
K,2,0,0,T
K,3,R,0,0
K,4,R,0,T
A,1,2,4,3
AATT,1
ESIZE,0.5E-3
NUMSTR,NODE,10
AMESH,ALL
EXTOPT,ESIZE,16,0,
VROTAT,ALL,,,,,,1,2,360
EPLOT

NSEL,S,LOC,Z,T
CP,1,VOLT,ALL
*GET,N_TOP,NODE,O,NUM,MIN
D,ALL,VOLT,1
NSEL,S,LOC,Z,0
CP,2,VOLT,ALL
*GET,N_BOT,NODE,0,NUM,MIN
D,ALL,VOLT,0
NSEL,ALL
FINISH
/SOLU
ANTYPE,STATIC
SOLVE
FINI
/PREP7
DDELE,N_TOP,VOLT
ET,3,CIRCU94,0
R=3000
R,1,R
N,1
TYPE,3
REAL,1
E,1,N_TOP
ET,4,CIRCU94,4
V=1
R,2,V
N,2
TYPE,4
REAL,2
E,N_TOP,N_BOT,2
FINI
/SOLU
DDELE,N_TOP,VOLT
D,1,VOLT,0
ANTYPE,TRANS
NSUB,100
TIME,2E-3
TINTP,,0.25,0.5,0.5
KBC,1
OUTRES,ALL,ALL
SOLVE
FINI
/POST26
NUMVAR,20
NSOL,2,N_TOP,VOLT,,V_C
PRVAR,2
PLVAR,2

hy1966

我觉得是我在circu94单元内定义节点的问题，在circu94单元之前的命令应该没有什么问题，其实我对下面的这段代码不是很理解，求详细的解释！
ET,3,CIRCU94,0
R=3000
R,1,R
N,1,0,0                   ？这个节点是 ...
steven710 发表于 2010-9-5 16:40

看ANSYS HELP帮助中定义节点的命令，不要用你已经在使用的节点编号，不在实体模型上定义节点，也不要在同一位置定义两个节点
N, NODE, X, Y, Z, THXY, THYZ, THZX
再看一下circu94单元如何定义电阻、电感、电容、电压源的节点，这些节点表示你定义的单元是如何连接的

steven710

如何做到不在实体上定义节点？如果不在实体上定义节点，那么N命令后面的X,Y,Z又都表示什么含义呢？ 4# hy1966

hy1966

如何做到不在实体上定义节点？如果不在实体上定义节点，那么N命令后面的X,Y,Z又都表示什么含义呢？ 4# hy1966
steven710 发表于 2010-9-6 10:05

呵呵，你去认真看了ANSYS HELP没有？
坐标值

steven710

我认真看了呀，我的意思是，如果指定坐标值，那不就是定义在实体模型上了么？ 6# hy1966

steven710

hy1966老师，下面是我修改后的程序，解决了上面提到的问题，可是，出现了新的问题，大概意思是说我在节点上施加了多余的约束？这个应该怎么解决呀，我程序的目的是希望在压电陶瓷片两端加脉冲激励信号，观察陶瓷片两端电压的变化，请多多指教，十分感谢！
/PREP7
/VIEW,1,,-1
/TITLE, RADIUS PIEZO TRANSIENT ANALYSIS


ET,1,SOLID5,0
ET,2,MESH200
KEYOPT,2,1,6
KEYOPT,2,2,0


MP,DENS,1,7500
MP,DMPR,1,0.005
MP,PERX,1,804.6
MP,PERY,1,804.6
MP,PERZ,1,659.7


TB,PIEZ,1
TBDATA,16,10.5
TBDATA,14,10.5
TBDATA,3,-4.1
TBDATA,6.-4.1
TBDATA,9,14.1


TB,ANEL,1
TBDATA,1,13.2E10,7.1E10,7.3E10
TBDATA,7,13.2E10,7.3E10
TBDATA,12,11.5E10
TBDATA,16,3.0E10
TBDATA,19,2.6E10
TBDATA,21,2.6E10


R=2E-3
T=1E-3
K,1,0,0,0
K,2,0,0,T
K,3,R,0,0
K,4,R,0,T
A,1,2,4,3
AATT,1
ESIZE,0.5E-3
NUMSTR,NODE,10
AMESH,ALL


EXTOPT,ESIZE,16,0,
VROTAT,ALL,,,,,,1,2,360
EPLOT




NSEL,S,LOC,Z,T
CP,1,VOLT,ALL
*GET,N_TOP,NODE,O,NUM,MIN
D,ALL,VOLT,1


NSEL,S,LOC,Z,0
CP,2,VOLT,ALL
*GET,N_BOT,NODE,0,NUM,MIN
D,ALL,VOLT,0
NSEL,ALL
FINISH


/SOLU
ANTYPE,STATIC
SOLVE
FINI


/PREP7
DDELE,N_TOP,VOLT


ET,3,CIRCU94,4,3
R,1,0,50,10E-9,5E-9,20E-9,5E-6
N,1
TYPE,3 $ REAL,1
E,N_TOP,N_BOT,1
FINI


/SOLU
ANTYPE,TRANS
NSUB,100
TIME,2E-3
TINTP,,0.25,0.5,0.5
KBC,1
OUTRES,ALL,ALL
SOLVE
FINI


/POST26
NUMVAR,20
NSOL,2,N_TOP,VOLT,,V_C
PRVAR,2
PLVAR,2 6# hy1966

hy1966

删除多余的约束
定义了耦合组只要在耦合组代表节点上加电压就可以了。你在耦合组的所有节点上加负载（这不对），而删除时只删除了主节点上的电压负载，其它的仍然在

steven710

十分感谢hy1966老师，您的建议真是一针见血呀！学了很多东西，再次感谢！ 9# hy1966

qoohotboy

可以請問一下壓電材料加入circu94元素的電阻，在諧波外力所產生的壓電電壓是"峰值"還是"峰對峰值"?