24小时热门版块排行榜

返回列表

清光彩纱帐

木虫 (小有名气)

应助: 0 (幼儿园)
金币: 4630.7
帖子: 90
在线: 43.6小时
虫号: 1725790
注册: 2012-03-29
性别: GG
专业: 有色金属冶金

[求助] Ansys网格、载荷都是对称的，可结果不对称已有1人参与

采用的旋转复制的方法生成旋转对称的体，载荷是120度相位差的电流（我做过尝试，不管电极1上的电流初始相位是0、120或240结果都完全一样，所以可以认为载荷是旋转对称）加在三根电极上，但结果却不对称，复制得到的两个1/3模型结果对称，与被复制的1/3模型的结果不对称。。下面是命令流。
/TITLE,THREE PHASE KRL
/FILNAME,THREE PHASE KRL             !定义工作文件名
KEYW,MAGEDG,1                      !指定磁场分析
/UNITS,MKS
/SHAPP,OFF
/PREP7
!建模和划分网格
!指定单元类型
ET,1,SOLID117,1
!设置材料参数
MP,RSVX,1,0.89          !料层比电阻
MP,RSVX,2,1.93E-4       !渣层电阻率
MP,RSVX,3,1.82E-4       !硅层电阻率???????
MP,RSVX,4,6E-4          !电极电阻率
MP,RSVX,5,4.35E13       !SiO+CO电阻率
MP,RSVX,6,1E-5          !电弧电阻率(自己估的)????????
MP,RSVX,7,4.351E13       !空气电阻率
MP,MURX,1,1       !料层磁导率
MP,MURX,2,1       !渣层磁导率
MP,MURX,3,1       !镍铁层磁导率
MP,MURX,4,1       !电极磁导率
MP,MURX,5,1       !SiO+CO磁导率
MP,MURX,6,1       !电弧磁导率
MP,MURX,7,1       !空气磁导率
I=46562          !操作电流
!设置结构参数
He=3.000 !电极高度
De=0.96    !电极直径
Dm=6.00 !熔体直径
Dp=2.45 !极心圆直径
Da=1.800 !电弧区直径
Th1=1.500 !料层厚度
Th2=0.200 !渣层厚度
Th3=0.400 !金属液层进度
Dep=1.000 !电极插入深度
Hc=0.5    !铜瓦高度
L=16.000 !向右复制的距离
!开始建模
/PNUM,VOLU,1
/PNUM,AREA,1
/PNUM,LINE,1
WPOFFS,Dp/2
SPHERE,0,Da/2,0,180
VSBW,1
VDELE,2,,,1
CYLIND,0,De/2,0,Th1-Dep+He,0,180
VOVLAP,ALL
VADD,2,4
AADD,16,14,8,17
L,12,14
AL,2,20,22
VSBA,1,4
WPCSYS,-1
CYLIND,0,Dm/2,0,Th1,0,60
CYLIND,0,Dm/2,0,-1*Th2,0,60
CYLIND,0,Dm/2,-1*Th2,-1*(Th2+Th3),0,60
WPOFFS,,,Th1+Hc
ASBW,18
VOVLAP,ALL
VADD,9,7
WPCSYS,-1
WPOFFS,Dp/2
WPROTA,,,90
VSBW,1
VSBW,2
WPCSYS,-1
WPOFFS,,,Th1
SPHERE,0,Dm,0,60
VOVLAP,ALL
NUMCMP,ALL
/TRLCY,VOLU,1,7

!给图形区域指定单元和材料属性
VSEL,S,VOLU,,5             !选中料层
VATT,1,,1                   !将单元属性1、材料属性1赋予所选中的体
VSEL,S,VOLU,,9             !选中渣层
VATT,2,,1                   !将单元属性1、材料属性2赋予所选中的体
VSEL,S,VOLU,,6          !选中金属层
VATT,3,,1                   !将单元属性1、材料属性3赋予所选中的体
VSEL,S,VOLU,,2,4,2          !选中电极
VATT,4,,1                   !将单元属性1、材料属性4赋予所选中的体
VSEL,S,VOLU,,8             !选中SiO+CO
VATT,5,,1                   !将单元属性1、材料属性5赋予所选中的体
VSEL,S,VOLU,,1,3,2          !选中电弧
VATT,6,,1                   !将单元属性1、材料属性6赋予所选中的体
VSEL,S,VOLU,,7             !选中空气
VATT,7,,1                   !将单元属性1、材料属性7赋予所选中的体

!开始划分网格
ALLSEL,ALL
MSHAPE,0,3D
MSHKEY,1
!划分料层
LESIZE,65,,,6
LESIZE,34,,,8
LESIZE,69,,,8
LESIZE,17,,,9
LCCAT,32,52
LCCAT,30,51
VSWEEP,5,39,40
!划分电弧
LESIZE,24,,,6
LESIZE,25,,,6
VMESH,1
VMESH,3
!划分电极
LESIZE,46,,,8
LESIZE,56,,,3
LESIZE,47,,,8
LESIZE,57,,,3
VSWEEP,2,5,3
VSWEEP,4,17,16
!划分Si层
LESIZE,19,,,12
LESIZE,54,,,4
LESIZE,23,,,8
VSWEEP,6,27,28
!划分SiO+CO
MSHAPE,1,3D
MSHKEY,0
SMRTSIZE,3
VMESH,8
!划分渣
VMESH,9
!划分空气
SMRTSIZE,3
VMESH,7

!镜像复制
VSYMM,Y,ALL

!旋转复制
CSYS,1
NROTAT,ALL
ALLSEL,ALL
VGEN,3,ALL,,,,120,,,,
CSYS,0
NUMMRG,ALL

!透明化
/TRLCY,VOLU,1,7
/TRLCY,VOLU,1,16
/TRLCY,VOLU,1,25
/TRLCY,VOLU,1,34
/TRLCY,VOLU,1,43
/TRLCY,VOLU,1,52

!开始求解
/SOLU
!求解设置
ANTYPE,HARM
HARFRQ,60
HROPT,FULL
HROUT,OFF
OUTPR,BASIC,ALL,
!开始加载边界条件
!加载交流电
ASEL,R,AREA,,3
ASEL,A,AREA,,16
ASEL,A,AREA,,55
ASEL,A,AREA,,61
NSLA,S,1
CP,1,VOLT,ALL
*GET,P1,NODE,,NUM,MIN
F,P1,AMPS,I,0
ASEL,S,AREA,,99
ASEL,A,AREA,,105
ASEL,A,AREA,,143
ASEL,A,AREA,,149
NSLA,S,1
CP,2,VOLT,ALL
*GET,P2,NODE,,NUM,MIN
F,P2,AMPS,-0.5*I,-I*SQRT(3)/2
ASEL,S,AREA,,187
ASEL,A,AREA,,193
ASEL,A,AREA,,231
ASEL,A,AREA,,237
NSLA,S,1
CP,3,VOLT,ALL
*GET,P3,NODE,,NUM,MIN
F,P3,AMPS,-1*I/2,I*SQRT(3)/2
!加载磁平行边界条件
ALLSEL,ALL
NSEL,EXT
D,ALL,AZ,0
ALLSEL,ALL
SOLVE
FINISH

模型图3.PNG

回复此楼