Ansys-Fatigue疲劳分析问题
浏览：44810
最近在做Anasys-Fatigue疲劳分析问题，结果做的不是很理想，结构允许的循环次数总是我输入S-N曲线20个点(S-N)中的一个点的N值，而不是处于两N值之间的某值，参考资料上输入的S-N曲线结果的循环次数就可以，请问问题出在哪？下面是疲劳分析的命令流，材料是黄铜。
dpzt=0.0064
hpzt=0.0001
dcu=0.008
hcu=0.0001
hpzthcu=hpzt+hcu
meshlinesize=0.0006
! 进入前处理器
/prep7
！设置单元类型
et,1,solid5,3
et,2,mesh200,6
et,3,solid45
! 设置材料参数
tb,anel,1,,,0
tbdata,1,1.2035e11,7.5179e10,7.509e10
tbdata,7,1.2035e11,7.509e10
tbdata,12,1.1087e11
tbdata,16,2.2584e10
tbdata,19,2.1053e10
tbdata,21,2.1053e10
mp,dens,1,7750
mp,perx,1,919
mp,pery,1,919
mp,perz,1,826
tb,piez,1
tbdata,3,-5.3512
tbdata,6,-5.3512
tbdata,9,15.784
tbdata,14,12.259
tbdata,16,12.259
mp,ex,3,118e9
mp,prxy,3,0.35
mp,dens,3,8920
!实体建模(基板要分为环形和圆形)
cylind,0,dpzt,0,hcu,0,90 !CYLIND，RAD1，RAD2，Z1，Z2，THETA1，THETA2（以工作平面圆点为中心生成圆柱体）
cylind,0,dpzt,0,hcu,90,180
cylind,0,dpzt,0,hcu,180,270
cylind,0,dpzt,0,hcu,270,360
cylind,dpzt,dcu,0,hcu,0,90
cylind,dpzt,dcu,0,hcu,90,180
cylind,dpzt,dcu,0,hcu,180,270
cylind,dpzt,dcu,0,hcu,270,360
cylind,0,dpzt,hcu,hpzthcu,0,90
cylind,0,dpzt,hcu,hpzthcu,90,180
cylind,0,dpzt,hcu,hpzthcu,180,270
cylind,0,dpzt,hcu,hpzthcu,270,360
cylind,0,dpzt,0,-hpzt,0,90
cylind,0,dpzt,0,-hpzt,90,180
cylind,0,dpzt,0,-hpzt,180,270
cylind,0,dpzt,0,-hpzt,270,360
allsel,all
vglue,all
!划分网格
!辅助网格划分
type,2
esize,meshlinesize
allsel,all
asel,s,loc,z,0
mshape,0,2d
mshkye,1
amesh,all
!基板网格划分
allsel,all
type,3
mat,3
vsel,s,loc,z,0,hcu
vsweep,all
!压电晶体网格划分
allsel,all
type,1
mat,1
vsel,s,loc,z,hcu,hpzthcu
vsweep,all
allsel,all
vsel,s,loc,z,0,-hpzt
vsweep,all
aclear,all
!施加载荷和边界条件
allsel,all
asel,s,,,127
asel,a,,,130
asel,a,,,123
asel,a,,,23
da,all,all
dtran
allsel,all
nsel,s,loc,z,hpzthcu
cp,1,volt,all
d,all,volt,100
nsel,s,loc,z,hcu
cp,2,volt,all
d,all,volt,0
nsel,s,loc,z,-hpzt
cp,4,volt,all
d,all,volt,-100
allsel,all
/solu
antype,static
EQSLV,SPAR
solve
FINISH
/POST1
FTSIZE,1,1,2
!设置S-N曲线
FP,1,10,1000,100000000, , ,
FP,7, , , , , ,
FP,13, , , , , ,
FP,19, ,
FP,21,325000000,292500000,100000000, , ,
FP,27, , , , , ,
FP,33, , , , , ,
FP,39, ,
FS,3207,1,1,1,0,0,0,0,0,0,
set,1,last
FSNODE,3207,1,2,
FE,1,72000000,1,
FTCALC,1
命令流中S-N曲线是理论曲线，主要是3个点，之间是线性关系，输入S-N点是，完全可以在3点之间再取其他17个点，可结果允许循环次数的还是最后一个点的N值：1E8,不知道问题出在哪了，希望哪位高手指教:handshake 。
dpzt=0.0064
hpzt=0.0001
dcu=0.008
hcu=0.0001
hpzthcu=hpzt+hcu
meshlinesize=0.0006
! 进入前处理器
/prep7
！设置单元类型
et,1,solid5,3
et,2,mesh200,6
et,3,solid45
! 设置材料参数
tb,anel,1,,,0
tbdata,1,1.2035e11,7.5179e10,7.509e10
tbdata,7,1.2035e11,7.509e10
tbdata,12,1.1087e11
tbdata,16,2.2584e10
tbdata,19,2.1053e10
tbdata,21,2.1053e10
mp,dens,1,7750
mp,perx,1,919
mp,pery,1,919
mp,perz,1,826
tb,piez,1
tbdata,3,-5.3512
tbdata,6,-5.3512
tbdata,9,15.784
tbdata,14,12.259
tbdata,16,12.259
mp,ex,3,118e9
mp,prxy,3,0.35
mp,dens,3,8920
!实体建模(基板要分为环形和圆形)
cylind,0,dpzt,0,hcu,0,90 !CYLIND，RAD1，RAD2，Z1，Z2，THETA1，THETA2（以工作平面圆点为中心生成圆柱体）
cylind,0,dpzt,0,hcu,90,180
cylind,0,dpzt,0,hcu,180,270
cylind,0,dpzt,0,hcu,270,360
cylind,dpzt,dcu,0,hcu,0,90
cylind,dpzt,dcu,0,hcu,90,180
cylind,dpzt,dcu,0,hcu,180,270
cylind,dpzt,dcu,0,hcu,270,360
cylind,0,dpzt,hcu,hpzthcu,0,90
cylind,0,dpzt,hcu,hpzthcu,90,180
cylind,0,dpzt,hcu,hpzthcu,180,270
cylind,0,dpzt,hcu,hpzthcu,270,360
cylind,0,dpzt,0,-hpzt,0,90
cylind,0,dpzt,0,-hpzt,90,180
cylind,0,dpzt,0,-hpzt,180,270
cylind,0,dpzt,0,-hpzt,270,360
allsel,all
vglue,all
!划分网格
!辅助网格划分
type,2
esize,meshlinesize
allsel,all
asel,s,loc,z,0
mshape,0,2d
mshkye,1
amesh,all
!基板网格划分
allsel,all
type,3
mat,3
vsel,s,loc,z,0,hcu
vsweep,all
!压电晶体网格划分
allsel,all
type,1
mat,1
vsel,s,loc,z,hcu,hpzthcu
vsweep,all
allsel,all
vsel,s,loc,z,0,-hpzt
vsweep,all
aclear,all
!施加载荷和边界条件
allsel,all
asel,s,,,127
asel,a,,,130
asel,a,,,123
asel,a,,,23
da,all,all
dtran
allsel,all
nsel,s,loc,z,hpzthcu
cp,1,volt,all
d,all,volt,100
nsel,s,loc,z,hcu
cp,2,volt,all
d,all,volt,0
nsel,s,loc,z,-hpzt
cp,4,volt,all
d,all,volt,-100
allsel,all
/solu
antype,static
EQSLV,SPAR
solve
FINISH
/POST1
FTSIZE,1,1,2
!设置S-N曲线
FP,1,10,1000,100000000, , ,
FP,7, , , , , ,
FP,13, , , , , ,
FP,19, ,
FP,21,325000000,292500000,100000000, , ,
FP,27, , , , , ,
FP,33, , , , , ,
FP,39, ,
FS,3207,1,1,1,0,0,0,0,0,0,
set,1,last
FSNODE,3207,1,2,
FE,1,72000000,1,
FTCALC,1
命令流中S-N曲线是理论曲线，主要是3个点，之间是线性关系，输入S-N点是，完全可以在3点之间再取其他17个点，可结果允许循环次数的还是最后一个点的N值：1E8,不知道问题出在哪了，希望哪位高手指教:handshake 。
