变压器低压绕组加载10*SIN(100*PI*T)+5电流密度,出现警告,是BFE命令用错了吧?菜鸟,求高手指点

浏览:52470 回答:5

命令如下,大神帮忙啊!感激不尽
! transformer for 31500kVA/35kV/400V, developed by Ou
/FILNAME,Single
!Part 1: Initializing the parameter of geometry and properties.
Pi=3.1415926
Ri=94E-03 !Radius of iron core ;
G1=10E-03 !Gap between core and inner HV winding;
G2=5E-03 !Gap between LV and HV;
G3=5.5E-03 !Gap between winding and yoke;
Hhv=460E-03 !Height of HV winding;
Hlv=420E-03 !Height of LV winding;
Hi=Max(Hlv,Hhv)+2*G3
Whv=(12*2.6+11*5)*1E-03 !Width of inner HV winding;
Wlv=(3+5)*1E-03 !Width of LV winding;
Wt=450E-03 !Length between two core center lines;
Nlv=24 !Half turns of LV winding;
Slv=3*16E-06 !Area of core conductor about LV;
Flv=(Nlv*Slv)/(2*Wlv*Hlv) !Fill factor abou LV;
Shv=5.31E-06 !Areas of core conductor about HV;
Nhv=12*175 !Turns of HV winding;
Fhv=(Nhv*Shv)/(2*Whv*Hhv) !Fill factor abou HV;
!Part 2: Creating half FEM model;
/PREP7
R1=Ri+G1
R2=R1+Wlv+G2
Rectng,0,Ri,-Hhv,Hhv
Rectng,R1,R1+Wlv,-Hlv,Hlv !LV winding;
Rectng,R2,R2+Whv,-Hhv,Hhv !HV winding;
Rectng,0,Wt-Ri,-Hi,Hi !Air;
AOVLAP,ALL
NumCmp,All
!Part 3: Defining the properties.
ET,1,53,,,1 ! Air,iron
ET,2,53,3,,1 ! Primary (LV) winding;
ET,3,53,3,,1 ! Secondary (HV) winding;
MP,MURX,1,1 ! No.1= Air;
MP,MURX,3,2000 ! No.4= Iron;
Rsvx_Cu=1.24E-08 ! No_Load=1.24;
UIMP,2,MURX,RSVX,,1,Rsvx_cu !No.2=copper conductor;
R,1,Wlv*Hlv,Nlv,,1,Flv !LV windings;
R,2,Whv*Hhv,Nhv,,-1,Fhv
ASEL,S,,,1
AATT,3,,1 !iron core
ASEL,S,,,2 !LV winding;
AATT,2,1,2
Asel,S,,,3 !HV winding;
AATT,2,2,3
Asel,s,,,4 !Air;
AATT,1,,1
Allsel,All
!Part 4:Meshing model.
!Part 4.1: Meshing coil;
Lsel,S,,,6,8,2
Lesize,All,,,24
Lsel,S,,,5,7,2
Lesize,All,,,1
Lsel,S,,,10,12,2
Lesize,All,,,175
Lsel,S,,,9,11,2
Lesize,All,,,12
MShape,0,2D !Qualaterial elements;
Mshkey,1 !Mapped Meshing;
ASel,s,Type,,2,3
Amesh,All
!Part 4.2:Meshing air region;
Smrt,1
Mshape,1,2D
MshKey,0
ASel,S,mat,,1
ASel,a,mat,,3
Amesh,All
Allsel,All
!Part 6: Applying the boundary conditions.
ESel,S,TYPE,,2
NSLE,S
CP,1,VOLT,ALL ! Couple Current inside LV winding;
nsel,s,loc,x,0
D,all,az,0
Allsel,All
FINISH
/SOLU
EQSLV,JCG,1E-9
T=1E-8
C=0
N=80
PI=2*ASIN(1)
CON=1/4000
NEQIT,1
*CREATE,LOAD
TIME,T
I=10*SIN(100*PI*T)+5
ESel,S,TYPE,,2
!F,228,AMPS,I
BFE,ALL,JS,,I
T=T+CON
C=C+1
OUTRES,ALL,1
SOLVE
*END
*DO,I,1,81
*USE,LOAD
*ENDDO
FINISH
出现警告,是BFE命令用错了吧
1.PNG
邀请回答 我来回答

全部回答

(5)
默认 最新
晴博

我没运行你的命令,不过从字面警告上来看,你虽然在某个节点上加电流,但是这个节点貌似没选中,建议用一下全选命令。
另外你第一次提到的那个警告,应该是说你加载电流密度的单元不能加载电流密度,所以是你的单元类型的问题。线圈的单元类型第一个关键项设置为0就可以加载电流密度的。
2014年5月28日
评论 点赞
09wangzhen

试了很多次,还是有问题,我又换了个例子,加载用F命令加载的,先耦合volt自由度,然后将电流加载到界面上某一个节点上,会出现这种警告:specified forces from previous F command were applied only to defined and selected nodal,程序如下,能否帮看看呀
!/BACH,LIST
/TITLE,2D SOLENOID ACTUATOR STATIC ANALYSIS !定义工作标题
/FILNAME,EMAGE_2D,1 ! 定义工作文件名
KEYW,MAGNOD,1 ! 指定磁场分析
/PREP7
ET,1,PLANE53,,,1 ! 指定单元类型并定义单元选项
ET,2,PLANE53,2,,1 ! 指定单元类型并定义单元选项
MP,MURX,1,1 ! 定义空气材料属性
MP,MURX,2,1000 ! 定义铁芯材料属性
MP,MURX,3,1 ! 定义线圈材料属性
MP,MURX,4,2000 ! 定义电枢材料属性
MP,RSVX,3,3E-8 ! 定义线圈电阻率
MP,RSVX,4,70E-8 ! 定义电枢电阻率
N=650 ! 定义参数
TA=0.75
TB=0.75
TC=0.50
TD=0.75
WC=1
HC=2
GAP=0.25
SPACE=0.25
WS=WC+2*SPACE
HS=HC+0.75
W=TA+WS+TC
HB=TB+HS
H=HB+GAP+TD
ACOIL=WC*HC ! 线圈横截面积(CM2)
R,1,ACOIL*.01**2,N,,1,.95 ! 定义线圈实常数
/PNUM,AREA,1 ! 打开面区域编号
/PNUM,AREA,1
RECTNG,0,W,0,TB ! 生成几何模型
RECTNG,0,W,TB,HB
RECTNG,TA,TA+WS,0,H
RECTNG,TA+SPACE,TA+SPACE+WC,TB+SPACE,TB+SPACE+HC
AOVLAP,ALL
RECTNG,0,W,0,HB+GAP
RECTNG,0,W,0,H
AOVLAP,ALL
NUMCMP,AREA ! 压缩编号
APLOT
SAVE,Emage_2D_geom.db ! 保存几何模型到文件
ASEL,S,AREA,,2 ! 给线圈区域赋予材料特性
AATT,3,1,2,0
ASEL,S,AREA,,1 ! 给电枢区域赋予材料特性
ASEL,A,AREA,,12,13
AATT,4,1,1
ASEL,S,AREA,,3,5 ! 给铁芯区域赋予材料特性
ASEL,A,AREA,,7,8
AATT,2,1,1,0
/PNUM,MAT,1 ! 打开材料编号
ALLSEL,ALL
APLOT
SMRTSIZE,4 ! 设置智能化划分网格等级L
AMESH,ALL ! 划分自由网格
SAVE,Emage_2D_Mesh.db ! 保存网格单元数据到文件
ESEL,S,MAT,,4 ! 选择电枢上的所有单元
CM,ARM,ELEM ! 将电枢单元定义为一个组件
FMAGBC,'ARM' ! 给电枢加力边界条件
ESEL,S,MAT,,3
NSLE,S
CP,1,volt,ALL
CM,COIL,ELEM ! 将线圈单元定义为一个组件
ALLSEL,ALL
ARSCAL,ALL,,,.01,.01,1,,0,1 ! 改变单位制为MSK单位(米)
NSEL,EXT ! 选择外层节点
D,ALL,AZ,0 ! 施加磁力线平行边界条件
FINISH
/SOLU
EQSLV,JCG,1E-9 ! 选择雅克比求解器,误差为1×10-9
T=1E-8 ! 设置开始求解时间
C=0 ! 定义一个计数器
N=80 ! 设置一个时间步长
PI=2*ASIN(1) ! 获得π的值
CON=2*PI/N ! 设置时间增量大小
NEQIT,1 ! 设置时间步的迭代次数为1
*CREATE,LOAD ! 创建宏设置载荷步
TIME,T ! 设置载荷步终止时间
I=4*SIN(T) ! 计算激励电流
F,511,AMPS,I ! 把电流加到节点12上
T=T+CON ! 时间递增
C=C+1 ! 计数器递增
OUTRES,ALL,1 ! 设置输出控制
*IF,C,EQ,((N*.75)+1),THEN ! 条件判断
OUTPR,,1 ! 输出打印结果
*ELSEIF,C,EQ,(N+1),THEN ! 条件转移
OUTPR,,1
*ELSE
OUTPR,,0
*ENDIF
SOLVE
*END
*DO,I,1,81 ! 循环命令
*USE,LOAD ! 执行宏
*ENDDO
FINISH
2014年5月10日
评论 点赞
liusongoku

你看到bfe命令的用法了,val1代表jx,val2代表jy,val3代表jz,所以你只需要给val3赋值就行了。其它为空
2014年5月7日
评论 点赞
09wangzhen

谢谢你,那我想在z方向加电流为i,那该如何使用命令
2014年5月7日
评论 点赞
liusongoku
有一个问题,你用的53单元建模,那么肯定是二维模型了,那为什么
BFE,ALL,JS,,I
加载了x方向的电流?二维分析只能加载z方向电流啊
For Lab = JS and STLOC = 1, VAL1, VAL2 and VAL3 are the X, Y, and Z components of current density (in the element coordinate system), and VAL4 is the phase angle.
2014年4月28日
评论 点赞

没解决?试试专家一对一服务

换一批
    App下载
    技术邻APP
    工程师必备
    • 项目客服
    • 培训客服
    • 平台客服

    TOP