费花花 发表于 2013-4-17 16:15:52

求助，三维块导体加环形电流

一个圆盘状的导体线圈


导体线圈内要加载环形电流，在开缝处，一个面耦合volt自由度，把电流加在面上的一个点上，另一个面设置volt=0，外面加球形的空气场（这个图里外面的空气场还没画）
可是求解完以后，没有结果显示。
我想问一下，这样加载电流对不对？加载荷时要修改单元局部坐标系吗？平行边界条件应该如何加？
求指教，非常感谢！

费花花 发表于 2013-4-17 16:16:40

我的命令流

/PREP7
!*
ET,1,SOLID97
!*
ET,2,SOLID97
!*
KEYOPT,1,1,1
KEYOPT,1,2,0
KEYOPT,1,5,0
!*
KEYOPT,1,1,1
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,1,,1
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,RSVX,1,,1.7e-8   
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,MURX,2,,1
RECTNG,0.02,0.06,0.005,-0.005,   
CYL4,0.02,0,0.005   
FLST,2,2,5,ORDE,2   
FITEM,2,1   
FITEM,2,-2
AADD,P51X   
K,1,0,0.005,0,
K,100,0,0.005,0,
K,101,0,-0.005,0,      
FLST,2,1,5,ORDE,1   
FITEM,2,3   
FLST,8,2,3
FITEM,8,100
FITEM,8,101
VROTAT,P51X, , , , , ,P51X, ,360,2,
/VIEW,1,1,1,1   
/ANG,1
/REP,FAST      
BLOCK,-0.001,0.001,0.006,-0.006,0.01,0.07,
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,-2
VSBV,P51X,       3
FLST,2,3,6,ORDE,3   
FITEM,2,2   
FITEM,2,4   
FITEM,2,-5
VADD,P51X   
SPH4,0,0,0.15   
FLST,2,2,6,ORDE,2   
FITEM,2,1   
FITEM,2,-2
VOVLAP,P51X
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       1, ,   1,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
CM,_Y,VOLU
VSEL, , , ,       3
CM,_Y1,VOLU
CMSEL,S,_Y
!*
CMSEL,S,_Y1
VATT,       2, ,   2,       0   
CMSEL,S,_Y
CMDELE,_Y   
CMDELE,_Y1
!*
ESIZE,0.005,0,
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU
VSEL, , , ,       1
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
ESIZE,0.01,0,   
MSHKEY,0
MSHAPE,1,3d
CM,_Y,VOLU
VSEL, , , ,       3
CM,_Y1,VOLU
CHKMSH,'VOLU'   
CMSEL,S,_Y
!*
VMESH,_Y1   
!*
CMDELE,_Y   
CMDELE,_Y1
CMDELE,_Y2
!*
! SAVE, mesh,db,D:\Users\   
LPLOT   
ASEL,S, , ,      24
NSLA,S,1
FINISH
/PREP7
FLST,4,35,1,ORDE,12
FITEM,4,2   
FITEM,4,9   
FITEM,4,-17
FITEM,4,145
FITEM,4,150
FITEM,4,159
FITEM,4,166
FITEM,4,176
FITEM,4,202
FITEM,4,-209
FITEM,4,312
FITEM,4,-323
CP,1,VOLT,P51X   
ASEL,S, , ,      30
NSLA,S,1
FINISH
/SOL
FLST,2,35,1,ORDE,10
FITEM,2,895
FITEM,2,902
FITEM,2,921
FITEM,2,-928
FITEM,2,1029
FITEM,2,1034
FITEM,2,1044
FITEM,2,1051
FITEM,2,1058
FITEM,2,-1078   
/GO
!*
D,P51X,VOLT,0   
ALLSEL,ALL   
!*
ANTYPE,3
HARFRQ,2.8e6,0,
NSUBST, ,   
KBC,0   
!*
FLST,2,1,1,ORDE,1   
FITEM,2,23089   
/GO
!*
F,P51X,AMPS,672,1164
/STATUS,SOLU
SOLVE   
FINISH

tsfly 发表于 2013-5-17 14:10:35

电流是无旋的，这样直接加载电流是不行，可以尝试加载电压，或者建立1/4或1/2模型，使求解域边界与导体横截面边界（对称边界）相切，在一端耦合电压自由度吗，加载总电流，另一端设基准电压为零即可。

perch88 发表于 2013-5-17 14:21:17

可以尝试36+96单元配合
另外可以简化为2d，能简化的尽量简化，不然3d网格划分增加很多麻烦，尤其新手不一定能做好

kkyukk 发表于 2013-5-23 11:03:13

我感觉得将单元坐标系改成柱坐标系吧。我也是新手。。。
页: [1]
查看完整版本: 求助，三维块导体加环形电流