声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 4147|回复: 9

[电磁场分析] 【转帖】Piezoelectric - 压电材料振动能量采集的有限元分析

[复制链接]
发表于 2012-9-11 10:18 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
Piezoelectric - 压电材料振动能量采集的有限元分析
200552317113389216.gif
回复
分享到:

使用道具 举报

 楼主| 发表于 2012-9-11 10:20 | 显示全部楼层


  1. !1) Initialisation of elements and material properties

  2. /TITLE, Piezoelectric Vibration Energy Harvester
  3. /UNITS,MKS

  4. accelx=0     !Acceleration in X
  5. accely=0     !Acceleration in Y
  6. accelz=2.5     !Acceleration in Z

  7. /COM,  
  8. /COM,   Finite Element Model of a Piezoelectric Energy Harvester
  9. /COM,  

  10. /PREP7

  11. /COM,  Define Element Types Used During FEA Analysis

  12. ET,1,SOLID5,0    !8-Node Coupled Field Piezoelectric Element For PZT
  13. ET,2,SOLID45    !8-Node Solid Element For Modelling Tungsten Mass and
  14.       !steel Shim
  15. ET,3,SHELL41    !4-Node Membrane Shell Element (Electrode), Not used in
  16.       !this project
  17. ET,4,TARGE170,0    !3D Target surface for modelling bonding contact
  18.       !between tungsten mass and beam
  19. ET,5,CONTA174    !Coupled field contact surface for modelling contact
  20.       !between tungsten mass and beam

  21. /COM,  Define Element Type Used During Electrical Analysis

  22. ET,6,CIRCU94,0    !Define Resistor

  23. /COM,      !Define Keyopts for Finite Elements

  24. KEYOPT,4,5,1    !Set TARGE 170 to solid-solid constraint
  25. KEYOPT,5,1,0    !UX,UY,UZ
  26. KEYOPT,5,2,2    !Contact algorithm: MPC
  27. KEYOPT,5,4,1    !Contact detection on nodal point
  28. KEYOPT,5,5,0    !No automated contact adjustment
  29. KEYOPT,5,7,0    !No element level time increment control
  30. KEYOPT,5,8,0    !No asymmetric contact selection
  31. KEYOPT,5,9,1    !Exclude both initial geometrical penetration or gap
  32.       !and offset
  33. KEYOPT,5,10,1    !Update contact stiffness each substep based on mean stress
  34.       !of underlying elements
  35. KEYOPT,5,12,5    !Bonded always contact behaviour

  36. /COM,  ZT-SH4 (Navy Type VI) Material Properties
  37. /COM,  iezo Systems Inc.

  38. EMUNIT, MKS    ! Free space permittivity

  39. MP,DENS,1,7500    !Material density

  40. MP,PERX,1,3130    !Relative permittivity MP,PERY,1,3130
  41. MP,PERZ,1,3400

  42. TB,ANEL,1,1,,1    !Flexibility matrix - inverted by Ansys
  43. TBDATA,1,16.5e-12,-4.78e-12,-8.45e-12
  44. TBDATA,7,16.5e-12,-8.45e-12
  45. TBDATA,12,20.7e-12
  46. TBDATA,16,43.5e-12
  47. TBDATA,19,43.5e-12
  48. TBDATA,21,42.6e-12

  49. TB,PIEZ,1    ! Piezoelectric constant matrix
  50. TBDATA,3,-6.622
  51. TBDATA,6,-6.622
  52. TBDATA,9,23.24
  53. TBDATA,11,17.034
  54. TBDATA,13,17.394

  55. MP,MURX,1,0    !False material properties to suppress error messages
  56. MP,KXX,1,0

  57. /COM,  Steel Shim Material Properties

  58. MP,DENS,2,7700    !Density
  59. MP,EX,2,207e9    !Young's Modulus
  60. MP,PRXY,2,0.3    !Poisson's ratio
  61. !MP,RSVX,2,72e-8  !Resistivity (Volt DOF will be coupled later)

  62. /COM,  Tungsten Mass Material Properties
  63. MP,DENS,3,17000    !Density
  64. MP,EX,3,540e9    !Mass deformation will be negligible
  65. MP,PRXY,3,0.28    !Poisson's ratio

  66. /COM,  Nickel Electrode Material Properties, not used in this
  67. /COM,  project

  68. R,1,0.2e-6    !Thickness of top nickel electrode
  69. R,2,0.2e-6    !Thickness of bottom nickel electrode

  70. MP,EX,4,207e9    !Young's Modulus
  71. MP,PRXY,4,0.3    !Poisson's ratio

  72. /COM,  Contact material properties

  73. REAL,3      !Set element real constant attribute pointer
  74. R,3,,,,,,
  75. RMORE,,,,,,
  76. RMORE,,0,,,,
  77. RMORE,0      !Set electrical contact conductance
  78. ALLSEL,ALL,ALL    !Select all entities

  79. /COM,  Define real constants for circuit elements

  80. R,4,180e3    !Define Resistor Properties (Defaults)

  81. !Application of material properties to solid geometry

  82. /COM,  Apply material properties to volumes of imported
  83. /COM,  ro/Engineer Wildfire 2 IGES model

  84. VSEL,S,VOLU,,2,,,0  !Select volume 2 (Bottom Piezo layer)
  85. VSEL,A,VOLU,,3,,,0  !Additionally select volume 3 (Top Piezo layer)
  86. VATT,1,1,1,0    !Apply PZT-5144 material properties

  87. VSEL,S,VOLU,,1,,,0  !Select volume 1(Steel centre shim)
  88. VATT,2,1,2,0    !Apply steel material properties

  89. VSEL,S,VOLU,,4,,,0  !Select volume 4 (Tungsten mass)
  90. VATT,3,1,2,0    !Apply tungsten material properties

  91. /COM,  Set mesh densities

  92. /COM,  ZT thickness density

  93. LSEL,S,LINE,,38,,,0  !Select bottom PZT layer thickness
  94. LSEL,A,LINE,,84,,,0  !Select top PZT layer thickness
  95. LESIZE,ALL,,,3    !Specify 3 divisions mesh density

  96. /COM,  Steel shim thickness density

  97. LSEL,S,LINE,,12,,,0
  98. LESIZE,Al1,,,2    !Specify 2 divisions mesh density

  99. /COM  Beam width density

  100. LSEL,S,LINE,,4,,,0
  101. LSEL,A,LINE,,10,,,0
  102. LSEL,A,LINE,,41,,,0
  103. LESIZE,ALL,,,4    !4 divisions mesh density

  104. /COM,  Beam length density

  105. LSEL,S,LINE,,2,,,0
  106. LSEL,A,LINE,,9,,,O
  107. LSEL,A,LINE,,81,,,0
  108. LESIZE,ALL,,,120  ! 120 divisions mesh density

  109. /COM,  Tungsten mass mesh density

  110. LSEL,S,LINE,,109
  111. LESIZE,ALL,,,3
  112. LSEL,S,LINE,,111,,,0
  113. LESIZE,ALL,,,2
  114. LSEL,S,LINE,,110,,,0
  115. LESIZE,ALL,,,4

  116. /COM,  Mesh all volumes and areas

  117. MSHKEY,1    !Set mapped mesh

  118. ALLSEL,ALL,VOLU   !Select all volumes to be meshed
  119. VMESH,ALL    !Mesh all volumes

  120. !Coupling DOF
  121. !The VOLT DOF is coupled between the dielectric and
  122. !conductive layers in the model. This procedure uses the CP command to couple the
  123. !DOF set of the nodes in the selected region, hence, the procedure must be performed after
  124. !meshing.

  125. /COM,  Couple voltage degree of freedom between layers

  126. ASEL,S,AREA,,41,,,0  !Select top electrode
  127. NSLA,S,1    !Select nodes related to top electrode
  128. NSEL,U,LOC,X,10.7e-3,30e-3
  129. CM,top_electrode,NODE   !Create component from nodes

  130. ASEL,S,AREA,,34,,,0  !Select bottom electrode
  131. NSLA,S,1    !Select nodes related to top electrode
  132. NSEL,U,LOC,X,10.7e-3,30e-3
  133. CM,bottom_electrode,NODE!Create component from nodes

  134. NSEL,S,NODE,,top_electrode
  135. NSEL,A,NODE,,bottom_electrode
  136. CP,I,VOLT,ALL    !Couple voltage degree of freedom for all nodes
  137. *GET,OUT ELECT,NODE,O,NUM,MIN !Get master node on top electrode

  138. VSEL,S,VOLU,,1,,,1  !Select centre shim
  139. NSLV,S,1    !Select nodes associated with volume
  140. CM,centre shim,NODE   !Create component from nodes
  141. CP,NEXT,VOLT,ALL   !Couple voltage degree of freedom for all nodes
  142. *GET,NCENTRE,NODE,O,NUM,MIN !Get master centre electrode node

  143. ALLSEL

  144. !D.1.5 Modelling contact of proof mass
  145. !In this section of the analysis, a target element is defined on the base of the tungsten
  146. !mass and a contactelement is defined on the top surface of the beam as described in ANSYS
  147. !Inc. (2004d).

  148. ASEL,S,AREA,,48,,,0   !Select bottom of tungsten mass
  149. NSLA,S,1    !Select nodes related to top electrode
  150. CM,target_surface,NODE   !Create component from nodes

  151. ALLSEL,ALL,ALL    !Select all entities

  152. /COM,  Define contact pair's

  153. NSEL,S,,,target_surface !Select nodes on base of tungsten mass
  154. TYPE,4      ! Set target element
  155. ESLN,S,0    !Select elements attached to nodes
  156. ESURF,ALL    !Create target elements

  157. ASEL,S,AREA,,41,,,0  !Select top electrode
  158. NSLA,S,1    !Select nodes related to top electrode
  159. TYPE,5      !Set contact element type
  160. ESLN,S,0    !Select elements attached to nodes
  161. ESURF,ALL    !Create contact elements
  162. ALLSEL  !Select all

  163. FINISH

  164. !2) Analysis

  165. !2.1) Static analysis

  166. /COM,  Static analysis to test regime

  167. /SOLU      !Enter solution pre-processor

  168. /COM,  Apply boundary conditions
  169. /COM  Completely rigid fixture, constrained Y for all

  170. NSEL,S,LOC,X,0    !Select nodes at fixture location
  171. D,ALL,UX,0,,,,UY,UZ,   !Constrain all DOF

  172. /COM,  Define Model symmetry

  173. NSEL,S,LOC,Y,0    !Select nodes at model symmetry line
  174. DSYM,SYMM,Y    !Apply symmetry boundary conditions
  175. ALLSELL,ALL,ALL    !Select all entities

  176. /COM,  Apply 5V to top and bottom electrode, OV to shim
  177. D,top_electrode,VOLT,5
  178. D,bottom_electrode,VOLT,5
  179. D,centre_shim,VOLT,0
  180. ALLSELL,ALL,ALL    !Select all entities
  181. ANTYPE,STATIC
  182. SOLVE

  183. FINISH

  184. !2.2) Modal analysis

  185. /COM  Modal analysis for determining resonance response of
  186. /COM,  harvester
  187. /SOLU
  188. /COM,  Delete DOF from static analysis
  189. DDELE,top_electrode,VOLT
  190. DDELE,bottom_electrode,VOLT

  191. ANTYPE,MODAL    !Modal analysis
  192. MODOPT,REDUC,8    !4 modes using reduced method
  193. MXPAND,8    !Expand a114 modes
  194. TOTAL,10,1
  195. ALLSEL
  196. D,NCENTRE,VOLT,0.0   !Ground centre electrode
  197. DMPRAT,0.031    !Constant damping ratio
  198. NSEL,ALL    !Select all nodes

  199. SOLVE      !Solve current load step

  200. FINISH

  201. !2.3) Harmonic analysis without resistor

  202. /COM,  Haitnonic analysis around first resonant frequency

  203. /SOLU
  204. ANTYP,HARM    !Set harmonic analysis
  205. HARFRQ,30,120    !Set frequency range
  206. NSUBST,36    !Set 2.5 Hz increments
  207. KBC,1      !Loads step changed
  208. D,NCENTRE,VOLT,0.0   !Ground centre electrode
  209. DMPRAT,0.031    !Constant damping ratio
  210. ALLSEL

  211. ACEL,accelx,accely,accelz !Set accelerations

  212. SOLVE
  213. FINISH

  214. /POST26
  215. NSOL,9, OUT_ELECT, VOLT,,output_voltage ! Store output power

  216. FINISH

  217. !2.4) Harmonic analysis with resistor

  218. /COM,  Harmonic analysis around first frequency with
  219. /COM,  resistance load

  220. /COM,  Create circuit geometry

  221. /PREP7

  222. N,,0,-5e-3,0,    !Defines node 1 for circuit
  223. N,,5e-3,-5e-3,0,  !Defines node 2 for circuit

  224. /COM,  Create circuit elements on geometry
  225. TYPE,6      !Set circuit element type (Resistor)
  226. REAL,4      !Set circuit element constants (Resistance)
  227. E,5505,5506    !Plot Resistor

  228. /ICSCALE,1,0.1    !Scale circuit elements to correct size

  229. /COM,  Couple voltage degrees of freedom between circuit
  230. /COM,  and FE model. Places resistor between top and bottom
  231. /COM   electrode

  232. NSEL,S,NODE,,5505  !Select input node in circuit
  233. CP,1,VOLT,ALL    !Couple voltage Degree of Freedom on top elect

  234. NSEL,S,NODE,,5506  !Select input node in circuit
  235. CP,4,VOLT,ALL    !Couple voltage Degree of Freedom on centre elect

  236. ALLSELL,ALL,ALL    !Select all entities

  237. FINISH

  238. /COM,  Harmonic analysis with resistor coupled to outer electrodes

  239. /SOLU

  240. ANTYP,HARM    !Set harmonic analysis

  241. /COM  Set constraints on DOF

  242. NSEL,S,LOC,X,0    !Select nodes at fixture location
  243. D,ALL,UX,O,,,,UY,UZ,   !Constrain all DOF

  244. /COM  Define Model symmetry

  245. NSEL,S,LOC,Y,0     !Select nodes at model symmetry line
  246. DSYM,SYMM,Y     !Apply symmetry boundary conditions
  247. ALLSELL,ALL,ALL   !Select all entities

  248. HARFRQ,30,120    !Set frequency range
  249. NSUBST,36    !Set 2.5 Hz increments
  250. KBC,1      !Loads step changed
  251. D,NCENTRE,VOLT,0.0   !Ground centre electrode
  252. DMPRAT,0.031    !Constant damping ratio
  253. ALLSEL

  254. ACEL,accelx,accely,accelz !Set accelerations

  255. SOLVE
  256. FINISH

  257. !2.5) Spectral Analysis

  258. /COM,  Spectral analysis

  259. /SOLU

  260. ALLSEL
  261. DDELE,ALL,ALL    !Delete all previous constraints
  262. ACEL,0,0,0    !Set accelerations to zero (zero gravity)

  263. D,NCENTRE,VOLT,0.0   !Ground centre electrode

  264. NSEL,S,LOC,X,0    !Select nodes at fixture location
  265. D,ALL,UX,0,,,,UY,UZ,   !Constrain all DOF

  266. /COM,  Define Model symmetry

  267. NSEL,S,LOC,Y,0    !Select nodes at model symmetry line
  268. DSYM,SYMM,Y    !Apply symmetry boundary conditions
  269. ALLSELL,ALL,ALL    !Select all entities

  270. /COM,  Modal analysis
  271. /SOLU
  272. ANTYPE,MODAL    ! Mode-frequency analysis
  273. MODOPT,REDUC,,,,4   ! Householder, print first 3 reduced mode shapes
  274. MXPAND,4,    ! Expand first mode shape
  275. TOTAL,10,1
  276. OUTPR,BASIC,1
  277. SOLVE

  278. FINISH

  279. /COM,  Spectrum Analysis

  280. /SOLU
  281. ANTYPE,SPECTR    ! Spectrum analysis
  282. SPOPT,SPRS    ! Single point spectrum
  283. DMPRAT,0.031    !Constant damping ratio
  284. SED,0,0,1    ! Global Z-axis as spectrum direction
  285. SVTYPE,2    ! Seismic acceleration spectrum
  286. MCOMB,SRSS,0.001,

  287. /COM  Define frequency table (from table 4.1)

  288. FREQ,12,31,32,34,35,57,59,67,68 !Frequency points 1-9 (Hz)
  289. FREQ,69,70,72,73  ! Frequency 10-13 (Hz)

  290. /COM,  Define acceleration table to match frequency table (from
  291. /COM,  table 4.1)

  292. SV,0.031,0.167,0.052,0.067,0.067,0.059,0.168,0.054,0.072,0.284
  293. SV,0.031,0.078,0.078,0.068,0.07 ! Acceleration magnitudes (m/s2)

  294. OUTRES,ALL,ALL,    !Output all results to file
  295. ALLSEL

  296. SOLVE
  297. FINISH

  298. /Post26      !Enter post processing

  299. /INPUT,,mcom,,1,0   !Input combined modes results file

  300. FINISH
复制代码
 楼主| 发表于 2012-9-11 10:20 | 显示全部楼层
一秒时间历程的瞬态分析:
200552317363712545.gif
 楼主| 发表于 2012-9-11 10:22 | 显示全部楼层
电压:; \ 200552317405184815.gif
 楼主| 发表于 2012-9-11 10:24 | 显示全部楼层

  1. 瞬态分析的独立程序如下:

  2. /PREP7
  3. RESUME,PZT_geo,db,,0,0

  4. !D.1.1 Initialisation of elements and material properties

  5. /UNITS,MKS

  6. /COM,  
  7. /COM,   Finite Element Model of a Piezoelectric Energy Harvester
  8. /COM,  

  9. /PREP7

  10. /COM,    Define Element Types Used During FEA Analysis

  11. ET,1,SOLID5,0    !8-Node Coupled Field Piezoelectric Element For PZT
  12. ET,2,SOLID45    !8-Node Solid Element For Modelling Tungsten Mass and
  13.       !steel Shim
  14. ET,3,SHELL41    !4-Node Membrane Shell Element (Electrode), Not used in
  15.       !this project
  16. ET,4,TARGE170,0    !3D Target surface for modelling bonding contact
  17.       !between tungsten mass and beam
  18. ET,5,CONTA174    !Coupled field contact surface for modelling contact
  19.       !between tungsten mass and beam

  20. /COM,  Define Element Type Used During Electrical Analysis

  21. ET,6,CIRCU94,0    !Define Resistor

  22. /COM,      !Define Keyopts for Finite Elements

  23. KEYOPT,4,5,1    !Set TARGE 170 to solid-solid constraint
  24. KEYOPT,5,1,0    !UX,UY,UZ
  25. KEYOPT,5,2,2    !Contact algorithm: MPC
  26. KEYOPT,5,4,1    !Contact detection on nodal point
  27. KEYOPT,5,5,0    !No automated contact adjustment
  28. KEYOPT,5,7,0    !No element level time increment control
  29. KEYOPT,5,8,0    !No asymmetric contact selection
  30. KEYOPT,5,9,1    !Exclude both initial geometrical penetration or gap
  31.       !and offset
  32. KEYOPT,5,10,1    !Update contact stiffness each substep based on mean stress
  33.       !of underlying elements
  34. KEYOPT,5,12,5    !Bonded always contact behaviour

  35. /COM  ZT-SH4 (Navy Type VI) Material Properties
  36. /COM,  iezo Systems Inc.

  37. EMUNIT, MKS    ! Free space permittivity

  38. MP,DENS,1,7500    !Material density

  39. MP,PERX,1,3130    !Relative permittivity MP,PERY,1,3130
  40. MP,PERZ,1,3400

  41. TB,ANEL,1,1,,1    !Flexibility matrix - inverted by Ansys
  42. TBDATA,1,16.5e-12,-4.78e-12,-8.45e-12
  43. TBDATA,7,16.5e-12,-8.45e-12
  44. TBDATA,12,20.7e-12
  45. TBDATA,16,43.5e-12
  46. TBDATA,19,43.5e-12
  47. TBDATA,21,42.6e-12

  48. TB,PIEZ,1    ! Piezoelectric constant matrix
  49. TBDATA,3,-6.622
  50. TBDATA,6,-6.622
  51. TBDATA,9,23.24
  52. TBDATA,11,17.034
  53. TBDATA,13,17.394

  54. MP,MURX,1,0    !False material properties to suppress error messages
  55. MP,KXX,1,0

  56. /COM,  Steel Shim Material Properties

  57. MP,DENS,2,7700    !Density
  58. MP,EX,2,207e9    !Young's Modulus
  59. MP,PRXY,2,0.3    !Poisson's ratio
  60. !MP,RSVX,2,72e-8  !Resistivity (Volt DOF will be coupled later)

  61. /COM,  Tungsten Mass Material Properties
  62. MP,DENS,3,17000    !Density
  63. MP,EX,3,540e9    !Mass deformation will be negligible
  64. MP,PRXY,3,0.28    !Poisson's ratio

  65. /COM  Nickel Electrode Material Properties, not used in this
  66. /COM,  project

  67. R,1,0.2e-6    !Thickness of top nickel electrode
  68. R,2,0.2e-6    !Thickness of bottom nickel electrode

  69. MP,EX,4,207e9    !Young's Modulus
  70. MP,PRXY,4,0.3    !Poisson's ratio

  71. /COM,  Contact material properties

  72. REAL,3      !Set element real constant attribute pointer
  73. R,3,,,,,,
  74. RMORE,,,,,,
  75. RMORE,,0,,,,
  76. RMORE,0      !Set electrical contact conductance
  77. ALLSEL,ALL,ALL    !Select all entities

  78. /COM,  Define real constants for circuit elements

  79. R,4,180e3    !Define Resistor Properties (Defaults)

  80. /COM,  Apply material properties to volumes of imported
  81. /COM,  ro/Engineer Wildfire 2 IGES model

  82. VSEL,S,VOLU,,1,,,0  !Select volume 2 (Bottom Piezo layer)
  83. VSEL,A,VOLU,,3,,,0  !Additionally select volume 3 (Top Piezo layer)
  84. VATT,1,1,1,0    !Apply PZT-5144 material properties

  85. VSEL,S,VOLU,,2,,,0  !Select volume 1(Steel centre shim)
  86. VATT,2,1,2,0    !Apply steel material properties

  87. VSEL,S,VOLU,,4,,,0  !Select volume 4 (Tungsten mass)
  88. VATT,3,1,2,0    !Apply tungsten material properties

  89. /COM,  Set mesh densities

  90. /COM,  ZT thickness density

  91. LSEL,S,LINE,,10,,,0  !Select bottom PZT layer thickness
  92. LSEL,A,LINE,,26,,,0  !Select top PZT layer thickness
  93. LESIZE,ALL,,,3    !Specify 3 divisions mesh density

  94. /COM,  Steel shim thickness density

  95. LSEL,S,LINE,,18,,,0
  96. LESIZE,Al1,,,2    !Specify 2 divisions mesh density

  97. /COM  Beam width density

  98. LSEL,S,LINE,,1,,,0
  99. LSEL,A,LINE,,5,,,0
  100. LSEL,A,LINE,,13,,,0
  101. LESIZE,ALL,,,4    !4 divisions mesh density

  102. /COM,  Beam length density

  103. LSEL,S,LINE,,2,,,0
  104. LSEL,A,LINE,,6,,,O
  105. LSEL,A,LINE,,14,,,0
  106. LESIZE,ALL,,,120  ! 120 divisions mesh density

  107. /COM,  Tungsten mass mesh density

  108. LSEL,S,LINE,,38
  109. LESIZE,ALL,,,3
  110. LSEL,S,LINE,,33,,,0
  111. LESIZE,ALL,,,2
  112. LSEL,S,LINE,,34,,,0
  113. LESIZE,ALL,,,4

  114. /COM,  Mesh all volumes and areas

  115. MSHKEY,1    !Set mapped mesh

  116. ALLSEL,ALL,VOLU   !Select all volumes to be meshed
  117. VMESH,ALL    !Mesh all volumes

  118. /COM,  Couple voltage degree of freedom between layers

  119. ASEL,S,AREA,,12,,,0  !Select top electrode
  120. NSLA,S,1    !Select nodes related to top_electrode
  121. NSEL,U,LOC,X,10.7e-3,30e-3
  122. CM,top_electrode,NODE   !Create component from nodes

  123. ASEL,S,AREA,,1,,,0  !Select bottom electrode
  124. NSLA,S,1    !Select nodes related to top_electrode
  125. NSEL,U,LOC,X,10.7e-3,30e-3
  126. CM,bottom_electrode,NODE!Create component from nodes

  127. NSEL,S,NODE,,top_electrode
  128. NSEL,A,NODE,,bottom_electrode
  129. CP,1,VOLT,ALL    !Couple voltage degree of freedom for all nodes
  130. *GET,OUT_ELECT,NODE,0,NUM,MIN !Get master node on top_electrode

  131. VSEL,S, , ,       2
  132. !VSEL,S,VOLU,,2,,,1  !Select centre shim
  133. NSLV,S,1    !Select nodes associated with volume
  134. CM,centre_shim,NODE   !Create component from nodes
  135. CP,NEXT,VOLT,ALL   !Couple voltage degree of freedom for all nodes
  136. *GET,NCENTRE,NODE,0,NUM,MIN !Get master centre electrode node

  137. ALLSEL

  138. ASEL,S,AREA,,17,,,0   !Select bottom of tungsten mass
  139. NSLA,S,1    !Select nodes related to top electrode
  140. CM,target_surface,NODE   !Create component from nodes

  141. ALLSEL,ALL,ALL    !Select all entities

  142. /COM,  Define contact pair's

  143. NSEL,S,,,target_surface !Select nodes on base of tungsten mass
  144. TYPE,4      ! Set target element
  145. ESLN,S,0    !Select elements attached to nodes
  146. ESURF,ALL    !Create target elements

  147. ASEL,S,AREA,,12,,,0  !Select top electrode
  148. NSLA,S,1    !Select nodes related to top electrode
  149. TYPE,5      !Set contact element type
  150. ESLN,S,0    !Select elements attached to nodes
  151. ESURF,ALL    !Create contact elements
  152. ALLSEL  !Select all

  153. FINISH

  154. CSYS,0

  155. /SOLU      !Enter solution pre-processor

  156. NSEL,S,LOC,X,0    !Select nodes at fixture location
  157. D,ALL,UX,0,,,,UY,UZ,   !Constrain all DOF

  158. ALLSELL,ALL,ALL    !Select all entities

  159. !*
  160. NT=100        
  161. DT=0.01           !time increment
  162. *dim,ac,,NT
  163. /input,vibration_in,txt
  164. !*
  165. /SOLU
  166. NSUBST,1, , ,1         
  167. OUTRES,ALL,1         !output result of each substep
  168. ANTYPE,TRANS         !transient
  169. TINTP,,0.25,0.5,0.5,0.5,0,,,,,,
  170. *do,i,1,NT
  171. ACEL,0,0,ac(i)
  172. TIME,i*DT
  173. solve
  174. *enddo
  175. !*
  176. /POST26
  177. NSOL,2,6069,U,Z,
  178. /GRID,1
  179. PLVAR,2
复制代码
 楼主| 发表于 2012-9-11 10:25 | 显示全部楼层
震动波

vibration_in.rar

825 Bytes, 下载次数: 11

 楼主| 发表于 2012-9-11 10:27 | 显示全部楼层
部分GUI结果:
2005523185546497.gif
20055231811174981.gif
20055231883830821.gif
发表于 2012-9-12 08:58 | 显示全部楼层
谢谢楼主
发表于 2012-9-16 23:20 | 显示全部楼层
请问这个材料是从学术会议的报告中摘出来的么?
发表于 2014-5-30 10:03 | 显示全部楼层
几何模型建立的程序没有
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 00:40 , Processed in 0.075032 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表