想问一个post26后处理时的报错（瞬态分析）

wuming1112

我加载的条件是同幅值的两个正弦信号。但是相位差90度。一个sinwt一个coswt
算了40个周期后，后处理一点就报错。错误请见图。各位大侠帮我分析下

wuming1112

有人知道吗？

yueshenabc

我也出现这样的问题，正在解决中，知道了告诉我。

wuming1112

兄弟你输出设置是全部吗？我的设成了只输出节点解。和这个有关系吗？

16443

把你选用的单元、求解时候的警告或者错误信息也说出来才能解决。

wuming1112

单元的话结构中的金属基体我是选用solid95 单元
压电陶瓷模块选用的是solid5单元

16443

你这样说看不出哪里错误来，把完整过程命令贴出来

wuming1112

*SET,l1,
*SET,b1,
*SET,h1,
*set,l2,
*set,b2,
*set,h2,
*set,vol,0
*set,voh,200
*SET,l3,
*SET,l4,
*SET,h3,
*SET,long,
/config,nres,50000   !扩展结果输出
/PREP7  
ET,1,SOLID95
et,2,solid5,3
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,7.3e10
MPDATA,PRXY,1,,.34  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2700   

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,3,7500       !+     正向极化压电陶瓷      
MP,PERX,3,7.124e-09   
MP,PERY,3,7.124e-09   
MP,PERZ,3,5.841e-09   
   
TB,ANEL,3                    
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10                 

TB,PIEZ,3                 
tbdata,3,-4.1   
tbdata,6,-4.1   
tbdata,9,14.1   
tbdata,14,10.5   
tbdata,16,10.5           

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,4,7500            !-     反向极化压电陶瓷
MP,PERX,4,7.124e-09  
MP,PERY,4,7.124e-09   
MP,PERZ,4,5.841e-09
         
TB,ANEL,4        
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10  

TB,PIEZ,4           
tbdata,3,4.1           
tbdata,6,4.1   
tbdata,9,-14.1         
tbdata,14,-10.5      
tbdata,16,-10.5      

!~~~~~~~~~~~~

K,1,0,0,l1,
K,2,h1,0,l1,
K, 3,h1,0,-l1,   
K, 4,0,0,-l1,
*do,x,1,3
  lstr,x,x+1
*enddo
lstr,4,1
al,1,2,3,4
wpoff,0,long,0
CSWPLA,11,0,1,1,
CSYS,11,
K,5 ,h1,0,l3,
K, 6,h1+h3,0,l3,
K, 7,h1+h3,0,l4,
K, 8,h1+h3+h3,0,l4,  
K, 9,h1+h3+h3,0,-l4,
K, 10,h1+h3,0,-l4,
K, 11,h1+h3,0,-l3,
K, 12,h1,0,-l3,  
*do,x,5,11
  lstr,x,x+1
*enddo
lstr,12,5
al,5,6,7,8,9,10,11,12
VOFFST,2,b1-long, ,
csys,0
VOFFST,1,b1, ,
VSYMM,Y,1,2 , , ,0,0  
wpoff,0,-long,0
BLOCK,-h2,0,-b2,b2,-l2,l2,
VGEN,2,5, , , h1+h2, ,, ,1            
VGEN,2,5, , , h1+h2, ,l1/2, ,1            
VGEN,2,5, , ,h1+h2, ,-l1/2, ,1  
VADD,1,2,3,4   
NUMCMP,ALL
VGLUE,all
NUMCMP,ALL
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.004,0,
VSEL, , , ,1,3,1
VSWEEP,all
TYPE,   2   
MAT,       3
REAL,   
ESYS,       0   
SECNUM,   
VSEL, , , ,       4
VSWEEP,all
SMRT,6  
ESIZE,0.012,0,  
MSHAPE,1,3D
MSHKEY,0
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
VSEL, , , ,       5
VMESH,all


   
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~模态
ASEL,S,,,39     
NSLA,S,1
CM,n_39,NODE
cp,39,volt,all
allsel,all
ASEL,S,,,40   
NSLA,S,1
CM,n_40,NODE
cp,40,volt,all
allsel,all

ASEL,S,,,43     
NSLA,S,1
CM,n_43,NODE
cp,43,volt,all
allsel,all

ASEL,S,,,44     
NSLA,S,1
CM,n_44,NODE
cp,44,volt,all
allsel,all

ASEL,S,,,51     
NSLA,S,1
CM,n_51,NODE
cp,51,volt,all
allsel,all

ASEL,S,,,52     
NSLA,S,1
CM,n_52,NODE
cp,52,volt,all
allsel,all

ASEL,S,,,47     
NSLA,S,1
CM,n_47,NODE
cp,47,volt,all
allsel,all

ASEL,S,,,48   
NSLA,S,1
CM,n_48,NODE
cp,48,volt,all
allsel,all

d,n_52,volt,VOL
d,n_51,volt,VOL
d,n_39,volt,VOL
d,n_40,volt,VOL
d,n_43,volt,VOL
d,n_44,volt,VOL
d,n_47,volt,VOL
d,n_48,volt,VOL


ASEL,S, , ,6
NSLA,S,1
D,all, , , , , ,UX,UY,UZ, , ,  
ALLSEL,ALL
ASEL,S, , ,17
NSLA,S,1
D,all, , , , , ,UX,UY,UZ, , ,  
ALLSEL,ALL
fini
save


!~~~~~~~~~~~~~~~~~~~~~~~~瞬态
*SET,f,20000        !激振频率               
*SET,number,800      !总的载荷步数
*dim,v,array,number  !定义电压加载时输入的参数数组
*dim,t,array,number  !定义载荷步每步的时间
*do,i,1,number
  v(i) =voh*sin(2*i*acos(-1)/5) !将每个周期划分为8份。vohsin（2πft/8）
  t(i)=i/(5*f)                !载荷步结束时间
*enddo

/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
TINTP, ,0.25,0.5,0.5, , ,   !压电瞬态分析必须定义的
OUTPR,NSOL,ALL,              !只输出节点解
!~~~~~~~~~~~~~~~~~~~~~~加载
*do,j,1,number
TIME,t(j)  !载荷步结束时间         
AUTOTS,1  !自动步长开启（我不知道要不要开）
DELTIM,2e-6,0.5e-6,4e-6,1 !有关积分步长的设置
KBC,0                   !连续加载
TSRES,ERASE


!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~纵振压电陶瓷（中间上下对称两个）

d,n_39,volt,v(j)     !加正弦信号
d,n_40,volt,VOl
d,n_43,volt,v(j)      !加正弦信号
d,n_44,volt,VOl
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~弯曲激振压电陶瓷 （下面左右各一个）
d,n_51,volt,,v(j)!加余弦信号
d,n_52,volt,,VOl
d,n_47,volt,,v(j)!加余弦信号
d,n_48,volt,,VOl   

ALLSEL,ALL  
LSWRITE,j, !写入载荷步
*enddo
ALLSEL,ALL  
LSSOLVE,1,number,1,!求解
fini
/post26

wuming1112

这是我的命令流

16443

你的命令留没法运行，检查你计算时候的提示。

allenfieldin

可以直接点“确定”，我试过，没问题的

wuming1112

*SET,l1,64e-3   
*SET,b1,25e-3  
*SET,h1,33.5e-3
*set,l2,10e-3
*set,b2,4e-3
*set,h2,4e-3
*set,vol,0
*set,voh,200
*SET,l3,2.5e-3
*SET,l4,9e-3
*SET,h3,8e-3
*SET,long,8e-3
/config,nres,50000   !扩展结果输出
/PREP7  
ET,1,SOLID95
et,2,solid5,3
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,EX,1,,7.3e10
MPDATA,PRXY,1,,.34  
MPTEMP,,,,,,,,  
MPTEMP,1,0  
MPDATA,DENS,1,,2700   

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,3,7500       !+     正向极化压电陶瓷      
MP,PERX,3,7.124e-09   
MP,PERY,3,7.124e-09   
MP,PERZ,3,5.841e-09   
   
TB,ANEL,3                    
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10                 

TB,PIEZ,3                 
tbdata,3,-4.1   
tbdata,6,-4.1   
tbdata,9,14.1   
tbdata,14,10.5   
tbdata,16,10.5           

!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
MP,DENS,4,7500            !-     反向极化压电陶瓷
MP,PERX,4,7.124e-09  
MP,PERY,4,7.124e-09   
MP,PERZ,4,5.841e-09
         
TB,ANEL,4        
tbdata,1,1.32e11,7.1e10,7.3e10,0,0,0   
tbdata,7,1.32e11,7.3e10,0,0,0        
tbdata,12,11.5E10,0,0,0                 
tbdata,16,3.06E10,0,0                 
tbdata,19,2.56E10,0               
tbdata,21,2.56E10  

TB,PIEZ,4           
tbdata,3,4.1           
tbdata,6,4.1   
tbdata,9,-14.1         
tbdata,14,-10.5      
tbdata,16,-10.5      

!~~~~~~~~~~~~

K,1,0,0,l1,
K,2,h1,0,l1,
K, 3,h1,0,-l1,   
K, 4,0,0,-l1,
*do,x,1,3
  lstr,x,x+1
*enddo
lstr,4,1
al,1,2,3,4
wpoff,0,long,0
CSWPLA,11,0,1,1,
CSYS,11,
K,5 ,h1,0,l3,
K, 6,h1+h3,0,l3,
K, 7,h1+h3,0,l4,
K, 8,h1+h3+h3,0,l4,  
K, 9,h1+h3+h3,0,-l4,
K, 10,h1+h3,0,-l4,
K, 11,h1+h3,0,-l3,
K, 12,h1,0,-l3,  
*do,x,5,11
  lstr,x,x+1
*enddo
lstr,12,5
al,5,6,7,8,9,10,11,12
VOFFST,2,b1-long, ,
csys,0
VOFFST,1,b1, ,
VSYMM,Y,1,2 , , ,0,0  
wpoff,0,-long,0
BLOCK,-h2,0,-b2,b2,-l2,l2,
VGEN,2,5, , , h1+h2, ,, ,1            
VGEN,2,5, , , h1+h2, ,l1/2, ,1            
VGEN,2,5, , ,h1+h2, ,-l1/2, ,1  
VADD,1,2,3,4   
NUMCMP,ALL
VGLUE,all
NUMCMP,ALL
TYPE,   2   
MAT,       2
REAL,   
ESYS,       0   
SECNUM,
ESIZE,0.004,0,
VSEL, , , ,1,3,1
VSWEEP,all
TYPE,   2   
MAT,       3
REAL,   
ESYS,       0   
SECNUM,   
VSEL, , , ,       4
VSWEEP,all
SMRT,6  
ESIZE,0.012,0,  
MSHAPE,1,3D
MSHKEY,0
TYPE,   1   
MAT,       1
REAL,   
ESYS,       0   
SECNUM,
VSEL, , , ,       5
VMESH,all


   
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~模态
ASEL,S,,,39     
NSLA,S,1
CM,n_39,NODE
cp,39,volt,all
allsel,all
ASEL,S,,,40   
NSLA,S,1
CM,n_40,NODE
cp,40,volt,all
allsel,all

ASEL,S,,,43     
NSLA,S,1
CM,n_43,NODE
cp,43,volt,all
allsel,all

ASEL,S,,,44     
NSLA,S,1
CM,n_44,NODE
cp,44,volt,all
allsel,all

ASEL,S,,,51     
NSLA,S,1
CM,n_51,NODE
cp,51,volt,all
allsel,all

ASEL,S,,,52     
NSLA,S,1
CM,n_52,NODE
cp,52,volt,all
allsel,all

ASEL,S,,,47     
NSLA,S,1
CM,n_47,NODE
cp,47,volt,all
allsel,all

ASEL,S,,,48   
NSLA,S,1
CM,n_48,NODE
cp,48,volt,all
allsel,all

d,n_52,volt,VOL
d,n_51,volt,VOL
d,n_39,volt,VOL
d,n_40,volt,VOL
d,n_43,volt,VOL
d,n_44,volt,VOL
d,n_47,volt,VOL
d,n_48,volt,VOL


ASEL,S, , ,6
NSLA,S,1
D,all, , , , , ,UX,UY,UZ, , ,  
ALLSEL,ALL
ASEL,S, , ,17
NSLA,S,1
D,all, , , , , ,UX,UY,UZ, , ,  
ALLSEL,ALL
fini
save


!~~~~~~~~~~~~~~~~~~~~~~~~瞬态
*SET,f,20000        !激振频率               
*SET,number,800      !总的载荷步数
*dim,v,array,number  !定义电压加载时输入的参数数组
*dim,t,array,number  !定义载荷步每步的时间
*do,i,1,number
  v(i) =voh*sin(2*i*acos(-1)/5) !将每个周期划分为8份。vohsin（2πft/8）
  t(i)=i/(5*f)                !载荷步结束时间
*enddo

/SOL
ANTYPE,4
TRNOPT,FULL
LUMPM,0
TINTP, ,0.25,0.5,0.5, , ,   !压电瞬态分析必须定义的
OUTPR,NSOL,ALL,              !只输出节点解
!~~~~~~~~~~~~~~~~~~~~~~加载
*do,j,1,number
TIME,t(j)  !载荷步结束时间         
AUTOTS,1  !自动步长开启（我不知道要不要开）
DELTIM,2e-6,0.5e-6,4e-6,1 !有关积分步长的设置
KBC,0                   !连续加载
TSRES,ERASE


!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~纵振压电陶瓷（中间上下对称两个）

d,n_39,volt,v(j)     !加正弦信号
d,n_40,volt,VOl
d,n_43,volt,v(j)      !加正弦信号
d,n_44,volt,VOl
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~弯曲激振压电陶瓷 （下面左右各一个）
d,n_51,volt,,v(j)!加余弦信号
d,n_52,volt,,VOl
d,n_47,volt,,v(j)!加余弦信号
d,n_48,volt,,VOl   

ALLSEL,ALL  
LSWRITE,j, !写入载荷步
*enddo
ALLSEL,ALL  
LSSOLVE,1,number,1,!求解
fini
/post26

wuming1112

原帖由 16443 于 2009-5-27 08:42 发表

                               
登录/注册后可看大图

你的命令留没法运行，检查你计算时候的提示。

刚刚上传的命令流应该是能运行的。

wuming1112

/COM,ANSYS RELEASE 10.0    UP20050718       08:47:43    05/27/2009

*** WARNING ***                         CP =       5.188   TIME= 08:47:54
Structural elements without mid nodes usually produce much more         
  accurate results in quad or brick shape.                              

*** WARNING ***                         CP =       6.000   TIME= 08:47:57
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =       6.000   TIME= 08:47:57
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =      45.484   TIME= 08:48:46
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =      45.484   TIME= 08:48:46
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =      72.734   TIME= 08:49:18
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =      72.734   TIME= 08:49:18
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     101.125   TIME= 08:49:48
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     101.125   TIME= 08:49:48
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     132.453   TIME= 08:50:21
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     132.453   TIME= 08:50:21
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     164.016   TIME= 08:50:54
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     164.016   TIME= 08:50:54
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     200.969   TIME= 08:51:32
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     200.969   TIME= 08:51:33
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     237.906   TIME= 08:52:11
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     237.906   TIME= 08:52:11
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     274.812   TIME= 08:52:49
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     274.812   TIME= 08:52:49
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     311.734   TIME= 08:53:28
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     311.734   TIME= 08:53:28
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     348.578   TIME= 08:54:06
For material 2, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.                                                     

*** WARNING ***                         CP =     348.578   TIME= 08:54:06
For material 3, the value of permittivity (PERX) is less than 1.0, and  
will be interpreted as absolute permittivity (i.e., as input) by the   
elements of type 5.

wuming1112

错误里面都是这些类似的警告