谢谢各位帮助，但是我的问题还是没有解决。我还是先说第一个问题，以下是命令流文件，当r=5e-3时有就算的结果，但是当r=10e-3的时候就计算慢，且没有计算结果，实在不知道问题出在什么地方，
/BATCH
/TITLE,yuanhuan
! LGWRITE,'file','lgw','E:\pvdf\PVDFAC~1\',COMMENT
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@圆环材料 参数
mat1dens=8600
mat1yang=10.4e10
mat1person=0.37
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@质量块 参数
mat2dens=8600
mat2yang=10.4e10
mat2person=0.37
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@柱坐标底下的pvdf
mat3dens=1780
c33=8.1e9
c13=4.84e9
c23=4.84e9
c11=6.92e9
c12=4.38e9
c22=6.92e9
c66=1.38e9
c44=1.38e9
c55=1.38e9
e11=-6.5e-2
e21=-1.64e-2
e31=1.04e-2
e42=-3.8e-2
e63=-3.8e-2
fz=9.3e-11
fx=9.3e-11
fy=9.3e-11
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@结构参数
pi=3.1415926
tonghou=0.2e-3
pvdfhou=0.195e-3
r=5e-3
z=0.5e-3
t=6e-3
L=1e-3
r1=r-pvdfhou
r2=r+tonghou
r3=r+tonghou+pvdfhou
r4=3.5e-3
theta=pi/6
a=r-z-r4+sqrt(r4*r4-L*L)
!@@@@@@@@@@@@@@@@@@@@@@@网格划分的参数
tonghoufen=2
pvdfhoufen=2
hudu1fen=8
hudu2fen=8
hudu3fen=8
lfen=2
zfen=2
r4fen=8
!@@@@@@@@@@@@@@@@@@@@@@@@@@程序
/PREP7
ET,1,SOLID45
ET,2,SOLID45
ET,3,SOLID5
KEYOPT,3,1,3
KEYOPT,3,3,0
KEYOPT,3,5,0
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@ 定义第一种材料，圆环
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,mat1yang
MPDATA,PRXY,1,,mat1person
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,mat1dens
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@定义第二种材料，质量块
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,mat2yang
MPDATA,PRXY,2,,mat2person
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,mat2dens
!!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@定义第三种材料，pvdf
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,mat3dens
TB,ANEL,3,1,21,0
TBTEMP,0
TBDATA,,c11,c12,c13,,,
TBDATA,,c22,c23,,,,c33
TBDATA,,,,,c44,,
TBDATA,,c55,,c66,,,
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,PERX,3,,fx
MPDATA,PERY,3,,fy
MPDATA,PERZ,3,,fz
TB,PIEZ,3,,,0
TBMODIF,1,1,e11
TBMODIF,1,2,
TBMODIF,1,3,
TBMODIF,2,1,e21
TBMODIF,2,2,
TBMODIF,2,3,
TBMODIF,3,1,e31
TBMODIF,3,2,
TBMODIF,3,3,
TBMODIF,4,1,
TBMODIF,4,2,e42
TBMODIF,4,3,
TBMODIF,5,1,
TBMODIF,5,2,
TBMODIF,5,3,
TBMODIF,6,1,
TBMODIF,6,2,
TBMODIF,6,3,e63
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@生成关键点
K,1,0,0,0
K,2,r1*cos(theta),r1*sin(theta),0
K,3,r*cos(theta),r*sin(theta),0
k,4,r2*cos(theta),r2*sin(theta),0
k,5,r3*cos(theta),r3*sin(theta),0
k,6,r1*cos(theta),-r1*sin(theta),0
k,7,r*cos(theta),-r*sin(theta),0
k,8,r2*cos(theta),-r2*sin(theta),0
k,9,r3*cos(theta),-r3*sin(theta),0
k,10,L,a,0
k,11,L,sqrt(r*r-L*L),0
k,12,L,sqrt(r*r-L*L)+tonghou,0
K,13,-r1*cos(theta),r1*sin(theta),0
K,14,-r*cos(theta),r*sin(theta),0
k,15,-r2*cos(theta),r2*sin(theta),0
k,16,-r3*cos(theta),r3*sin(theta),0
k,17,-r1*cos(theta),-r1*sin(theta),0
k,18,-r*cos(theta),-r*sin(theta),0
k,19,-r2*cos(theta),-r2*sin(theta),0
k,20,-r3*cos(theta),-r3*sin(theta),0
k,21,-L,a,0
k,22,-L,sqrt(r*r-L*L),0
k,23,-L,sqrt(r*r-L*L)+tonghou,0
k,24,0,-r,0
k,25,0,-r2,0
k,26,0,r-z-r4,0
k,27,0,r-z-2*r4,0
!@@@@@@@@@@@@@@@@@@@@@@@@@@@@@点连成线
LSTR, 2, 3
LSTR, 3, 4
LSTR, 4, 5
LSTR, 6, 7
LSTR, 7, 8
LSTR, 8, 9
LSTR, 10, 11
LSTR, 11, 12
LSTR, 21, 22
LSTR, 22, 23
LSTR, 13, 14
LSTR, 14, 15
LSTR, 15, 16
LSTR, 17, 18
LSTR, 18, 19
LSTR, 19, 20
LSTR, 24, 25
LSTR, 26, 10
LSTR, 26, 21
LSTR, 26, 27
LARC,6,2,1,r1
!*
LARC,13,17,1,r1
!*
LARC,7,3,1,r
LARC,3,11,1,r
LARC,11,22,1,r
LARC,22,14,1,r
LARC,14,18,1,r
LARC,18,24,1,r
LARC,24,7,1,r
LARC,8,4,1,r2
LARC,4,12,1,r2
LARC,12,23,1,r2
LARC,23,15,1,r2
LARC,15,19,1,r2
LARC,19,25,1,r2
LARC,25,8,1,r2
!*
LARC,5,9,1,r3
!*
LARC,16,20,1,r3
LARC,21,10,26,r4
LARC,10,27,26,r4
LARC,27,21,26,r4
FLST,2,4,4
FITEM,2,38
FITEM,2,13
FITEM,2,34
FITEM,2,16
AL,P51X
FLST,2,4,4
FITEM,2,34
FITEM,2,12
FITEM,2,27
FITEM,2,15
AL,P51X
FLST,2,4,4
FITEM,2,27
FITEM,2,11
FITEM,2,22
FITEM,2,14
AL,P51X
FLST,2,4,4
FITEM,2,12
FITEM,2,33
FITEM,2,10
FITEM,2,26
AL,P51X
FLST,2,4,4
FITEM,2,15
FITEM,2,28
FITEM,2,17
FITEM,2,35
AL,P51X
FLST,2,3,4
FITEM,2,41
FITEM,2,19
FITEM,2,20
AL,P51X
FLST,2,3,4
FITEM,2,40
FITEM,2,18
FITEM,2,20
AL,P51X
FLST,2,3,4
FITEM,2,18
FITEM,2,39
FITEM,2,19
AL,P51X
FLST,2,4,4
FITEM,2,39
FITEM,2,7
FITEM,2,25
FITEM,2,9
AL,P51X
FLST,2,4,4
FITEM,2,25
FITEM,2,8
FITEM,2,32
FITEM,2,10
AL,P51X
FLST,2,4,4
FITEM,2,2
FITEM,2,31
FITEM,2,8
FITEM,2,24
AL,P51X
FLST,2,4,4
FITEM,2,37
FITEM,2,3
FITEM,2,30
FITEM,2,6
AL,P51X
FLST,2,4,4
FITEM,2,30
FITEM,2,2
FITEM,2,23
FITEM,2,5
AL,P51X
FLST,2,4,4
FITEM,2,1
FITEM,2,21
FITEM,2,4
FITEM,2,23
AL,P51X
FLST,2,4,4
FITEM,2,29
FITEM,2,17
FITEM,2,36
FITEM,2,5
AL,P51X
VOFFST,1,t, ,
VOFFST,2,t, ,
VOFFST,3,t, ,
VOFFST,4,-t, ,
VOFFST,5,t, ,
VOFFST,6,-t, ,
VOFFST,7,-t, ,
VOFFST,8,t, ,
VOFFST,9,t, ,
VOFFST,10,-t, ,
VOFFST,11,t, ,
VOFFST,12,-t, ,
VOFFST,13,t, ,
VOFFST,14,-t, ,
VOFFST,15,-t, ,
FLST,2,15,6,ORDE,2
FITEM,2,1
FITEM,2,-15
VGLUE,P51X
FLST,5,7,6,ORDE,7
FITEM,5,18
FITEM,5,20
FITEM,5,-21
FITEM,5,23
FITEM,5,-24
FITEM,5,26
FITEM,5,28
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 1, , 1, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,4,6,ORDE,4
FITEM,5,6
FITEM,5,16
FITEM,5,-17
FITEM,5,22
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 2, , 2, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FLST,5,4,6,ORDE,4
FITEM,5,1
FITEM,5,19
FITEM,5,25
FITEM,5,27
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT, 3, , 3, 0
CMSEL,S,_Y
CMDELE,_Y
CMDELE,_Y1
FLST,5,8,4,ORDE,8
FITEM,5,28
FITEM,5,-29
FITEM,5,35
FITEM,5,-36
FITEM,5,165
FITEM,5,-166
FITEM,5,180
FITEM,5,-181
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,hudu1fen, , , , ,1
!*
FLST,5,8,4,ORDE,8
FITEM,5,24
FITEM,5,26
FITEM,5,31
FITEM,5,33
FITEM,5,163
FITEM,5,-164
FITEM,5,172
FITEM,5,-173
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,hudu2fen, , , , ,1
!*
FLST,5,16,4,ORDE,15
FITEM,5,21
FITEM,5,-23
FITEM,5,27
FITEM,5,30
FITEM,5,34
FITEM,5,37
FITEM,5,-38
FITEM,5,42
FITEM,5,44
FITEM,5,52
FITEM,5,60
FITEM,5,124
FITEM,5,143
FITEM,5,175
FITEM,5,177
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,hudu3fen, , , , ,1
FLST,5,16,4,ORDE,16
FITEM,5,1
FITEM,5,3
FITEM,5,-4
FITEM,5,6
FITEM,5,11
FITEM,5,13
FITEM,5,-14
FITEM,5,16
FITEM,5,43
FITEM,5,45
FITEM,5,125
FITEM,5,161
FITEM,5,-162
FITEM,5,174
FITEM,5,178
FITEM,5,-179
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,pvdfhoufen, , , , ,1
!*
FLST,5,8,4,ORDE,8
FITEM,5,2
FITEM,5,5
FITEM,5,12
FITEM,5,15
FITEM,5,116
FITEM,5,159
FITEM,5,-160
FITEM,5,176
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,tonghoufen, , , , ,1
FLST,5,14,4,ORDE,12
FITEM,5,18
FITEM,5,-19
FITEM,5,25
FITEM,5,32
FITEM,5,39
FITEM,5,-41
FITEM,5,82
FITEM,5,-83
FITEM,5,156
FITEM,5,-158
FITEM,5,169
FITEM,5,171
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,r4fen, , , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,20
FITEM,5,84
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,r4fen, , , , ,1
!*
FLST,5,4,4,ORDE,4
FITEM,5,7
FITEM,5,9
FITEM,5,167
FITEM,5,-168
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,zfen, , , , ,1
FLST,5,4,4,ORDE,4
FITEM,5,8
FITEM,5,10
FITEM,5,68
FITEM,5,170
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,zfen, , , , ,1
FLST,5,26,4,ORDE,22
FITEM,5,46
FITEM,5,-49
FITEM,5,56
FITEM,5,-57
FITEM,5,64
FITEM,5,-65
FITEM,5,72
FITEM,5,-73
FITEM,5,80
FITEM,5,-81
FITEM,5,85
FITEM,5,-87
FITEM,5,91
FITEM,5,106
FITEM,5,115
FITEM,5,120
FITEM,5,-121
FITEM,5,128
FITEM,5,-130
FITEM,5,139
FITEM,5,144
FITEM,5,147
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,tfen, , , , ,1
!*
FLST,5,15,6,ORDE,4
FITEM,5,1
FITEM,5,6
FITEM,5,16
FITEM,5,-28
CM,_Y,VOLU
VSEL, , , ,P51X
CM,_Y1,VOLU
CHKMSH,'VOLU'
CMSEL,S,_Y
!*
MSHAPE,0,3d
MSHKEY,1
VMESH,_Y1
MSHKEY,0
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
! LGWRITE,'mesh','lgw','E:\JIASUD~1\060707\',COMMENT
CSYS,1
FLST,2,7392,1,ORDE,2
FITEM,2,1
FITEM,2,-7392
NROTAT,P51X
FINISH
/SOL
FLST,2,1,5,ORDE,1
FITEM,2,39
!*
/GO
DA,P51X,ALL,
!*
ANTYPE,2
!*
MSAVE,0
!*
MODOPT,LANB,5
EQSLV,SPAR
MXPAND,5, , ,0
LUMPM,0
PSTRES,0
!*
MODOPT,LANB,5,0,0, ,OFF
! /STATUS,SOLU
SOLVE
FINISH
/POST1
SET,LIST
! LGWRITE,'solve','lgw','E:\JIASUD~1\060707\',COMMENT |
发表于 2006-10-30 18:30
提升卡
变色卡
显身卡
发表于 2006-10-30 18:36
楼主