Maxwell 3D v11 VBScript Example

hottomson

参考官方的文档An M3DFS Macro Example编写了Maxwell 3D v11的脚本，与原文档格式基本一致，方便对照，主要的变动是原mac采用重新建立模型，这里采用修改模型的方法。没有采用中文编写，我英语水平超水（四级都没有过） ，语法、词汇错误就别指正了 。

这个问题不用宏，直接采用参数化分析就可以完成。该脚本仅做参考，不保证正常运行 。

遗留问题：关于Force的输出存在问题，以下语句不对：
oDesign.AddOutputVariable "ForceVal", "Force_magnitude"

在Maxwell 12中则可以使用如下语句：
Set oModule = oDesign.GetModule("OutputVariable")
oModule.CreateOutputVariable "ForceVal", "group1.Force_mag",  _
   "SolnOptnSetup : LastAdaptive", "Magnetostatic", Array()

注：初步发现v12的脚本和v11脚本有一些不同，使用12的朋友注意了，具体差异可以参考软件自带的online help。

[ 本帖最后由 hottomson 于 2008-6-20 11:58 编辑 ]

hottomson

1. 
   
2. Dim oAnsoftApp
   
3. Dim oDesktop
   
4. Dim oProject
   
5. Dim oDesign
   
6. Dim oEditor
   
7. Dim oModule
   
8. 'Define
   
9. Dim u_plate_length
   
10. Dim u_plate_width
    
11. Dim u_plate_thick
    
12. Dim u_plate_high
    
13. Dim armature_length
    
14. Dim armature_width
    
15. Dim armature_thick
    
16. Dim core_radius
    
17. Dim core_height
    
18. Dim coil_1_height
    
19. Dim coil_1_inner_rad
    
20. Dim coil_1_outer_rad
    
21. Dim coil_1_position
    
22. Dim no_current_steps
    
23. Dim no_turns_coil1
    
24. Dim current_min
    
25. Dim current_max
    
26. Dim no_gaps
    
27. Dim gap_min
    
28. Dim gap_max
    
29. u_plate_length = 100
    
30. u_plate_width = 60
    
31. u_plate_thick = 5
    
32. u_plate_high = 75
    
33. armature_length = 120
    
34. armature_width = 80
    
35. armature_thick = 5
    
36. core_radius = 10
    
37. core_height = 70
    
38. coil_1_height = 5
    
39. coil_1_inner_rad = 15
    
40. coil_1_outer_rad = 25
    
41. coil_1_position = 25
    
42. no_current_steps = 2
    
43. no_turns_coil1 = 1000
    
44. current_min = 0.5
    
45. current_max = 1
    
46. no_gaps = 2
    
47. gap_min = 0.01
    
48. gap_max = 1
    
49. Set oAnsoftApp = CreateObject("AnsoftMaxwell.MaxwellScriptInterface.11")
    
50. Set oDesktop = oAnsoftApp.GetAppDesktop()
    
51. oDesktop.RestoreWindow
    
52. Set oProject = oDesktop.NewProject
    
53. oProject.InsertDesign "Maxwell", "MaxwellDesign1", "Magnetostatic", ""
    
54. Set oDesign = oProject.SetActiveDesign("MaxwellDesign1")
    
55. Set oEditor = oDesign.SetActiveEditor("3D Modeler")
    
56. Dim temp,temp1,temp2,temp3,gap,gap_loop_index
    
57. gap_loop_index=0
    
58. temp = gap_max - gap_min
    
59. temp3 = no_gaps - 1
    
60. temp1 = temp/temp3
    
61. temp2 = temp1*gap_loop_index
    
62. gap = gap_min + temp2
    
63. Dim u_plate_x,u_plate_y
    
64. u_plate_x=u_plate_width/2
    
65. u_plate_y=u_plate_length/2
    
66. oEditor.CreateBox Array("NAME:BoxParameters", "CoordinateSystemID:=", -1, "XPosition:=",  _
    
67.   "0mm", "YPosition:=", "0mm", "ZPosition:=", "0mm", "XSize:=", CStrmm(u_plate_x), "YSize:=",  _
    
68.   CStrmm(u_plate_y), "ZSize:=", CStrmm(u_plate_thick)), Array("NAME:Attributes", "Name:=", "u_plate", _
    
69.   "Flags:=", "", "Color:=",  "(132 132 193)", "Transparency:=", 0, "PartCoordinateSystem:=", "Global", "MaterialName:=",  _
    
70.   "vacuum", "SolveInside:=", true)
    
71. '
    
72. 'Calculate dimensions for quarter model U Plate
    
73. '  
    
74. temp1 = u_plate_x*(-1)
    
75. temp2 = u_plate_thick*(-1)
    
76. oEditor.CreateBox Array("NAME:BoxParameters", "CoordinateSystemID:=", -1, "XPosition:=",  _
    
77.   CStrmm(u_plate_x), "YPosition:=", CStrmm(u_plate_y), "ZPosition:=", CStrmm(u_plate_thick), "XSize:=", CStrmm(temp1), "YSize:=",  _
    
78.   CStrmm(temp2), "ZSize:=", CStrmm(u_plate_high)), Array("NAME:Attributes", "Name:=", "u1", _
    
79.   "Flags:=", "", "Color:=",  "(132 132 193)", "Transparency:=", 0, "PartCoordinateSystem:=", "Global", "MaterialName:=",  _
    
80.   "vacuum", "SolveInside:=", true)
    
81. oEditor.Unite Array("NAME:Selections", "Selections:=", "u_plate,u1"), Array("NAME:UniteParameters", "CoordinateSystemID:=",  _
    
82.   -1, "KeepOriginals:=", false)
    
83. '
    
84. 'Calculate dimensions for armature
    
85. '  
    
86. temp = u_plate_thick
    
87. temp = temp + u_plate_high
    
88. temp = temp + gap
    
89. arm_y = armature_length/2
    
90. arm_x = armature_width/2
    
91. oEditor.CreateBox Array("NAME:BoxParameters", "CoordinateSystemID:=", -1, "XPosition:=",  _
    
92.   "0mm", "YPosition:=", "0mm", "ZPosition:=", CStrmm(temp), "XSize:=", CStrmm(arm_x), "YSize:=",  _
    
93.   CStrmm(arm_y), "ZSize:=", CStrmm(armature_thick)), Array("NAME:Attributes", "Name:=", "armature", "Flags:=", _
    
94.   "", "Color:=", "(132 132 193)", "Transparency:=", 0, "PartCoordinateSystem:=",  _
    
95.   "Global", "MaterialName:=", "vacuum", "SolveInside:=", true)
    
96. '
    
97. 'Calculate dimensions for the core
    
98. '
    
99. temp = u_plate_thick
    
100. oEditor.CreateRectangle Array("NAME:RectangleParameters", "CoordinateSystemID:=",  _
     
101.   -1, "IsCovered:=", true, "XStart:=", "0mm", "YStart:=", "0mm", "ZStart:=", CStrmm(temp), "Width:=",  _
     
102.   CStrmm(core_height), "Height:=", CStrmm(core_radius), "WhichAxis:=", "Y"), Array("NAME:Attributes", "Name:=",  _
     
103.   "core", "Flags:=", "", "Color:=", "(132 132 193)", "Transparency:=", 0, "PartCoordinateSystem:=",  _
     
104.   "Global", "MaterialName:=", "vacuum", "SolveInside:=", true)
     
105. oEditor.SweepAroundAxis Array("NAME:Selections", "Selections:=", "core"), Array("NAME:AxisSweepParameters", "CoordinateSystemID:=",  _
     
106.   -1, "DraftAngle:=", "0deg", "DraftType:=", "Round", "CheckFaceFaceIntersection:=",  _
     
107.   false, "SweepAxis:=", "Z", "SweepAngle:=", "90deg", "NumOfSegments:=", "8")
     
108. '
     
109. 'Calculate dimensions for coil_1
     
110. '
     
111. temp = u_plate_thick
     
112. temp = temp + coil_1_position
     
113. temp1 = coil_1_outer_rad - coil_1_inner_rad
     
114. oEditor.CreateRectangle Array("NAME:RectangleParameters", "CoordinateSystemID:=",  _
     
115.   -1, "IsCovered:=", true, "XStart:=", CStrmm(coil_1_inner_rad), "YStart:=", "0mm", "ZStart:=",  _
     
116.   CStrmm(temp), "Width:=", CStrmm(temp1), "Height:=", CStrmm(coil_1_height), "WhichAxis:=", "Y"), Array("NAME:Attributes", "Name:=",  _
     
117.   "coil1", "Flags:=", "", "Color:=", "(132 132 193)", "Transparency:=", 0, "PartCoordinateSystem:=",  _
     
118.   "Global", "MaterialName:=", "vacuum", "SolveInside:=", true)
     
119. oEditor.SweepAroundAxis Array("NAME:Selections", "Selections:=", "coil1"), Array("NAME:AxisSweepParameters", "CoordinateSystemID:=",  _
     
120.   -1, "DraftAngle:=", "0deg", "DraftType:=", "Round", "CheckFaceFaceIntersection:=",  _
     
121.   false, "SweepAxis:=", "Z", "SweepAngle:=", "90deg", "NumOfSegments:=", "8")
     
122.   
     
123. oEditor.Section Array("NAME:Selections", "Selections:=", "coil1"), Array("NAME:SectionToParameters", "CoordinateSystemID:=",  _
     
124.   -1, "SectionPlane:=", "ZX")
     
125. oEditor.ChangeProperty Array("NAME:AllTabs", Array("NAME:Geometry3DAttributeTab", Array("NAME:PropServers",  _
     
126. "Section1"), Array("NAME:ChangedProps", Array("NAME:Name", "Value:=", "term1"))))
     
127. oEditor.Section Array("NAME:Selections", "Selections:=", "coil1"), Array("NAME:SectionToParameters", "CoordinateSystemID:=",  _
     
128. -1, "SectionPlane:=", "YZ")
     
129. oEditor.ChangeProperty Array("NAME:AllTabs", Array("NAME:Geometry3DAttributeTab", Array("NAME:PropServers",  _
     
130.   "Section1"), Array("NAME:ChangedProps", Array("NAME:Name", "Value:=", "term2"))))
     
131. '
     
132. 'the new region defination
     
133. '
     
134. oEditor.CreateRegion Array("NAME:RegionParameters", "CoordinateSystemID:=", -1, "+XPadding:=", _
     
135.   "300", "-XPadding:=", "0", "+YPadding:=", "300", "-YPadding:=", "0", "+ZPadding:=",  _
     
136.   "100", "-ZPadding:=", "50"), Array("NAME:Attributes", "Name:=", "Region", "Flags:=", _
     
137.   "Wireframe#", "Color:=", "(255 0 0)", "Transparency:=", 0, "PartCoordinateSystem:=", _
     
138.   "Global", "MaterialName:=", "vacuum", "SolveInside:=", true)
     
139. oEditor.ChangeProperty Array("NAME:AllTabs", Array("NAME:Geometry3DAttributeTab", Array("NAME:PropServers",  _
     
140.   "armature", "u_plate", "core"), Array("NAME:ChangedProps", Array("NAME:Material", "Material:=",  _
     
141.   "steel_1010"))))
     
142. oEditor.ChangeProperty Array("NAME:AllTabs", Array("NAME:Geometry3DAttributeTab", Array("NAME:PropServers",  _
     
143.   "coil1"), Array("NAME:ChangedProps", Array("NAME:Material", "Material:=", "copper"))))
     
144. Set oModule = oDesign.GetModule("BoundarySetup")
     
145. '
     
146. 'Loop for current steps
     
147. '
     
148. loop_index = 0
     
149. loop_limit = no_current_steps
     
150. temp = current_max - current_min
     
151. temp3 =no_current_steps -1
     
152. temp1 = temp/temp3
     
153. temp2 = temp1*loop_index
     
154. curr = current_min + temp2
     
155. curr = curr * no_turns_coil1
     
156. oModule.AssignCurrent Array("NAME:curr1_in", "Current:=", CStrA(curr), "IsSolid:=",  _
     
157.   true, "Point out of terminal:=", false, "Objects:=", Array("term1"))
     
158. oModule.AssignCurrent Array("NAME:curr1_out", "Current:=", CStrA(curr), "IsSolid:=",  _
     
159.   true, "Point out of terminal:=", true, "Objects:=", Array("term2"))
     
160.   
     
161. Set oModule = oDesign.GetModule("MaxwellParameterSetup")
     
162. oModule.AssignForce Array("NAME:group1", "Is Virtual:=", true, "Reference CS:=",  _
     
163.   "Global", "Objects:=", Array("armature"))
     
164.   
     
165. oModule.AssignMatrix Array("NAME:Matrix1", Array("NAME:MatrixEntry", Array("NAME:MatrixEntry", "Source:=", _
     
166.   "curr1_in", "NumberOfTurns:=", "1")), Array("NAME:MatrixGroup"), "Objects:=", Array( _
     
167.   "coil1"))
     
168.   
     
169. Set oModule = oDesign.GetModule("AnalysisSetup")
     
170. oModule.InsertSetup "Magnetostatic", Array("NAME:SolnOptnSetup", "MaximumPasses:=", 10, "MinimumPasses:=",  _
     
171.   2, "MinimumConvergedPasses:=", 1, "PercentRefinement:=", 50, "LinearResidual:=", _
     
172.   1E-006, "ConvergenceDispParamId:=", -1, "SolveFieldOnly:=", false, "PercentError:=", _
     
173.   1, "SolveMatrixAtLast:=", true, "UseOutputVariable:=", false, "PreAdaptMesh:=",  _
     
174.   false, "NonLinearResidual:=", 0.01, "SolverType:=", 0, "MuNonLinearBH:=", true, "ComputeHc:=",  _
     
175.   false, "HcNonLinearBH:=", true)
     
176.   
     
177. oDesign.AddOutputVariable "ForceVal", "Force_magnitude" '这里不正确
     
178. oDesign.AddOutputVariable "MatrixVal", "L(curr1_in,curr1_in)"
     
179. Dim ForceVal,MatrixVal
     
180. '  
     
181. 'Create the file, open it for data entry, and output the column labels.
     
182. '
     
183. Set oFS = CreateObject("Scripting.FileSystemObject")
     
184. Set ofile1 = oFS.CreateTextFile ("force.txt")
     
185. line = chr(9) & "Force vs gap"
     
186. ofile1.WriteLine line
     
187. line = "GAP [ mm ]" & chr(9)& "FORCE [ Newtons ]"
     
188. ofile1.WriteLine line
     
189. Set ofile2 = oFS.CreateTextFile ("L11.txt")
     
190. line = chr(9) & "L_11 vs gap"
     
191. ofile2.WriteLine line
     
192. line = "GAP [ mm ]" & chr(9)& "INDUCTANCE [ Henries ]"
     
193. ofile2.WriteLine line
     
194. While gap_loop_index < no_gaps
     
195. temp = gap_max - gap_min
     
196. temp3 = no_gaps - 1
     
197. temp1 = temp/temp3
     
198. temp2 = temp1*gap_loop_index
     
199. gap = gap_min + temp2
     
200. temp = u_plate_thick
     
201. temp = temp + u_plate_high
     
202. temp = temp + gap
     
203. oEditor.ChangeProperty Array("NAME:AllTabs", Array("NAME:Geometry3DCmdTab", Array("NAME:PropServers",  _
     
204.   "armature:CreateBox:1"), Array("NAME:ChangedProps", Array("NAME:Position", "X:=", "0mm", "Y:=",  _
     
205.   "0mm", "Z:=", CStrmm(temp)))))
     
206.   
     
207. loop_limit = no_current_steps
     
208. While loop_index < loop_limit
     
209.   temp = current_max - current_min
     
210.   temp3 =no_current_steps -1
     
211.   temp1 = temp/temp3
     
212.   temp2 = temp1*loop_index
     
213.   curr = current_min + temp2
     
214.   curr = curr * no_turns_coil1
     
215.   
     
216.   line = "Current: "  & curr
     
217.   ofile1.WriteLine line
     
218.   ofile2.WriteLine line
     
219.   
     
220.   oDesign.ChangeProperty Array("NAME:AllTabs", Array("NAME:MaxwellTab", Array("NAME:PropServers",  _
     
221.    "BoundarySetup:curr1_in"), Array("NAME:ChangedProps", Array("NAME:Current", "Value:=",  _
     
222.    CStrA(curr)))))
     
223.   oDesign.ChangeProperty Array("NAME:AllTabs", Array("NAME:MaxwellTab", Array("NAME:PropServers",  _
     
224.    "BoundarySetup:curr1_out"), Array("NAME:ChangedProps", Array("NAME:Current", "Value:=",  _
     
225.    CStrA(curr)))))
     
226.   
     
227.   oDesign.Analyze "SolnOptnSetup"
     
228.   
     
229.   ForceVal = oDesign.GetOutputVariableValue("ForceVal", "SolnOptnSetup : LastAdaptive","", "")
     
230.   MatrixVal = oDesign.GetOutputVariableValue("MatrixVal", "SolnOptnSetup : LastAdaptive","", "")
     
231.   
     
232.   ForceVal = ForceVal*4
     
233.   line = gap & chr(9) & ForceVal
     
234.   ofile1.WriteLine line
     
235.   
     
236.   MatrixVal = MatrixVal*4*no_turns_coil1*no_turns_coil1
     
237.   line = gap & chr(9) & MatrixVal
     
238.   ofile2.WriteLine line
     
239.   loop_index = loop_index +1
     
240. Wend
     
241. gap_loop_index = gap_loop_index+1
     
242. Wend
     
243. ofile1.close
     
244. ofile2.close
     
245.   
     
246. Function CStrmm(nDig)
     
247. CStrmm = CStr(nDig)+"mm"
     
248. End Function
     
249. Function CStrA(nDig)
     
250. CStrA = CStr(nDig)+"A"
     
251. End Function
     

复制代码

[ 本帖最后由 hottomson 于 2008-6-20 17:19 编辑 ]

moon_night

v11 2D里面只在后处理里面用到的宏和3D里面的Script不一样吧
mac里面连Dim命令都不认

hottomson

你说vbscript，我以为是v12的2D了，v11 2D用的是macro，那你可以参考Maxwell 2D参数化后处理宏编写小技巧http://www.simwe.com/forum/thread-790983-1-1.html

Simol_wayne

Hottomson大哥编程相当厉害了！这看这满篇的程序我头就疼，至少我要先学下VB了（笨的可以了，VB俺也不懂）

samuelch

有maxwell v12 了？

zhaobo_ren

今年3月底公布的

samuelch

:lol :lol :lol

Simol_wayne

号称是Maxwell的最终完美版 绝唱了都