求助--关于压电陶瓷分析

eaglerna

   请问谁有关于压电陶瓷分析的实例，查遍了很多书都没有找到，请高手不吝赐教。
  还有请问压电陶瓷的参数怎么设置？极化方向怎么设置呢？

NASA

[ 本帖最后由 supervb 于 2006-12-20 19:23 编辑 ]

雪缘

1. /TITLE,VM237, PIEZOELECTRIC-CIRCUIT ANALYSIS
   
2. JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
   
3. /SHOW,JPEG
   
4. /COM, ---------------------------------------------------------------
   
5. /COM, FINITE ELEMENT MODEL OF A PIEZOELECTRIC CIRCULAR PLATE
   
6. /COM, ---------------------------------------------------------------
   
7. A=1E-3 ! RADIUS, M
   
8. T=0.1E-3 ! THICKNESS, M
   
9. /COM,
   
10. /COM, MATERIAL PROPERTIES OF LEAD ZIRCONATE TITANATE (PZT-5A)
    
11. /COM,
    
12. /COM, -- MATERIAL MATRICES (POLAR AXIS ALONG Y-AXIS): IEEE INPUT
    
13. /COM,
    
14. /COM, [s11 s13 s12 0 0 0 ] [ 0 d31 0 ] [ep11 0 0 ]
    
15. /COM, [s13 s33 s13 0 0 0 ] [ 0 d33 0 ] [ 0 ep33 0 ]
    
16. /COM, [s12 s13 s11 0 0 0 ] [ 0 d31 0 ] [ 0 0 ep11]
    
17. /COM, [ 0 0 0 s44 0 0 ] [ 0 0 d15]
    
18. /COM, [ 0 0 0 0 s66 0 ] [ 0 0 0 ]
    
19. /COM, [ 0 0 0 0 0 s44] [d15 0 0 ]
    
20. /COM,
    
21. /COM, - COMPLIANCE COEFFICIENTS, M2/N
    
22. S11=16.4E-12
    
23. S12=-5.74E-12
    
24. S13=-7.22E-12
    
25. S33=18.8E-12
    
26. S44=47.5E-12
    
27. S66=44.3E-12
    
28. /COM, - PIEZOELECTRIC STRAIN COEFFICIENTS, C/N
    
29. D15=5.84E-10
    
30. D31=-1.71E-10
    
31. D33=3.74E-10
    
32. /COM, - RELATIVE PERMITTIVITY AT CONSTANT STRESS
    
33. EP11=1730
    
34. EP33=1700
    
35. /COM, - DENSITY, KG/M3
    
36. RHO=7750
    
37. /NOPR
    
38. 
    
39. /PREP7
    
40. /COM,
    
41. /COM, -- MATERIAL MATRICES (POLAR AXIS ALONG Y-AXIS): ANSYS INPUT
    
42. /COM,
    
43. /COM, [s11 s13 s12 0 0 0 ] [ 0 d31 0 ] [ep11 0 0 ]
    
44. /COM, [s13 s33 s13 0 0 0 ] [ 0 d33 0 ] [ 0 ep33 0 ]
    
45. /COM, [s12 s13 s11 0 0 0 ] [ 0 d31 0 ] [ 0 0 ep11]
    
46. /COM, [ 0 0 0 s44 0 0 ] [d15 0 0 ]
    
47. /COM, [ 0 0 0 0 s44 0 ] [ 0 0 d15]
    
48. /COM, [ 0 0 0 0 0 s66] [ 0 0 0 ]
    
49. 
    
50. TB,ANEL,1,,,1 ! ANISOTROPIC ELASTIC COMPLIANCE MATRIX
    
51. TBDA,1,S11,S13,S12
    
52. TBDA,7,S33,S13
    
53. TBDA,12,S11
    
54. TBDA,16,S44
    
55. 
    
56. TB,PIEZ,1,,,1 ! PIEZOELECTRIC STRAIN MATRIX
    
57. TBDA,2,D31
    
58. TBDA,5,D33
    
59. TBDA,8,D31
    
60. TBDA,10,D15
    
61. 
    
62. TB,DPER,1,,,1 ! DIELECTRIC PERMITTIVITY AT CONSTANT STRESS
    
63. TBDA,1,EP11,EP33
    
64. 
    
65. TBLIS,ALL ! LIST INPUT AND CONVERTED MATRICES
    
66. 
    
67. MP,DENS,1,RHO ! DENSITY
    
68. 
    
69. ET,1,PLANE223,1001,,1 ! PIEZOELECTRIC AXISYMMETRIC ELEMENT TYPE
    
70. RECT,,A,,T
    
71. ESIZE,T ! MESH SOLID MODEL
    
72. NUMSTR,NODE,10
    
73. AMESH,1
    
74. ! APPLY STRUCTURAL BC, SIMPLY SUPPORTED PLATE
    
75. NSEL,S,LOC,X
    
76. D,ALL,UX,0
    
77. NSEL,S,LOC,X,A
    
78. NSEL,R,LOC,Y,0
    
79. D,ALL,UY,0
    
80. NSEL,ALL
    
81. ! APPLY ELECTRIC BC
    
82. NSEL,S,LOC,Y,0
    
83. CP,1,VOLT,ALL ! COUPLE BOTTOM ELECTRODE
    
84. *GET,NBOT,NODE,0,NUM,MIN ! GET MASTER NODE ON BOTTOM ELECTRODE
    
85. NSEL,S,LOC,Y,T
    
86. CP,2,VOLT,ALL ! COUPLE TOP ELECTRODE
    
87. *GET,NTOP,NODE,0,NUM,MIN ! GET MASTER NODE ON TOP ELECTRODE
    
88. NSEL,ALL
    
89. D,NBOT,VOLT,0 ! GROUND BOTTOM ELECTRODE
    
90. D,NTOP,VOLT,1 ! APPLY 1V LOAD ON TOP ELECTRODE
    
91. 
    
92. FINI
    
93. /SOLU
    
94. ANTYPE,STATIC ! STATIC ANALYSIS
    
95. SOLVE
    
96. FINI
    
97. 
    
98. /POST1
    
99. *GET,QT,NODE,NTOP,RF,CHRG ! GET TOTAL CHARGE ON TOP ELECTRODE
    
100. CP=ABS(QT) ! CAPACITANCE CP=Q/V, WHERE V=1V
     
101. EPZ0=8.854E-12 ! FREE SPACE PERMITTIVITY
     
102. PI=3.1415 ! PI CONSTANT
     
103. C=EP33*EPZ0*PI*A**2/T ! ANALYTICAL CAPACITANCE
     
104. /COM, 2-D CAPACITANCE (ANALYTICAL) =%C%, F
     
105. /COM, 2-D CAPACITANCE (ANSYS) = %CP%, F
     
106. FINI
     
107. 
     
108. /COM, ---------------------------------------------------------------
     
109. /COM, FINITE ELEMENT MODEL OF RLC-CIRCUIT
     
110. /COM, ---------------------------------------------------------------
     
111. /PREP7
     
112. DDELE,NTOP,VOLT ! DELETE VOLTAGE LOAD ON TOP ELECTRODE
     
113. ET,2,CIRCU94,0 ! DEFINE A RESISTOR
     
114. R=3000 ! RESISTANCE, OHM
     
115. R,1,R
     
116. N,1
     
117. TYPE,2
     
118. REAL,1
     
119. E,1,NTOP
     
120. 
     
121. ET,3,CIRCU94,1 ! DEFINE AN INDUCTOR
     
122. L=15 ! INDUCTANCE, H
     
123. R,2,L
     
124. N,2
     
125. TYPE,3
     
126. REAL,2
     
127. E,2,1
     
128. 
     
129. ET,4,CIRCU94,4 ! DEFINE A VOLTAGE SOURCE
     
130. V=1 ! VOLTAGE LOAD, V
     
131. R,3,V
     
132. N,3
     
133. TYPE,4
     
134. REAL,3
     
135. E,2,NBOT,3
     
136. FINI
     
137. 
     
138. /SOLU
     
139. ANTYPE,TRANS ! TRANSIENT ANALYSIS
     
140. NSUB,100 ! NUMBER OF TIME STEPS
     
141. TIME,2E-3 ! ANALYSIS TIME, S
     
142. TINTP,,0.25,0.5,0.5 ! INTEGRATION PARAMETERS FOR A PIEZOELECTRIC ANALYSIS
     
143. OUTRES,ALL,ALL
     
144. SOLVE
     
145. FINI
     
146. 
     
147. /GOPR
     
148. /COM, ANALYTICAL SOLUTION:
     
149. /COM, V_C = 1-EXP(-D*T)*COS(B*T)-D/B*EXP(-D*T)*SIN(B*T)
     
150. /COM, WHERE:
     
151. D=R/(2*L)
     
152. B=SQRT(1/(L*C)-D**2)
     
153. /NOPR
     
154. 
     
155. /POST26
     
156. NUMVAR,20
     
157. NSOL,2,NTOP,VOLT,,V_C_ANSYS
     
158. 
     
159. ! DERIVE EXACT SOLUTION
     
160. *DIM,WORK1,ARRAY,100
     
161. *DIM,WORK2,ARRAY,100
     
162. FILLDATA,3,,,,1 ! 1
     
163. EXP,4,1,,,,,,-D,-1 ! -EXP(-D*T)
     
164. PROD,5,1,,,,,,B ! B*T
     
165. VGET,WORK1(1),5
     
166. *VFUN,WORK2(1),COS,WORK1(1)
     
167. VPUT,WORK2(1),6 ! COS(B*T)
     
168. *VFUN,WORK2(1),SIN,WORK1(1)
     
169. VPUT,WORK2(1),7 ! SIN(B*T)
     
170. ADD,8,6,7,,,,,,D/B ! COS(B*T) + D/B*SIN(B*T)
     
171. PROD,9,4,8 ! -EXP(-D*T)*[COS(B*T) + D/B*SIN(B*T)]
     
172. ADD,10,3,9,,V_C_EXACT ! 1-EXP(-D*T)*[COS(B*T) + D/B*SIN(B*T)]
     
173. 
     
174. ! PRINT AND PLOT ANSYS AND EXACT VOLTAGE ACROSS THE PZT CAPACITOR
     
175. PRVAR,2,10
     
176. PLVAR,2,10
     
177. 
     
178. /NOPR
     
179. *DIM,VCE,ARRAY,5 ! EXACT SOLUTION FOR RESULTS TABLE
     
180. *DIM,VCA,ARRAY,5 ! ANSYS SOLUTION FOR RESULTS TABLE
     
181. 
     
182. *GET,VCE(1),VARI,10,RTIME,0.18E-3
     
183. *GET,VCE(2),VARI,10,RTIME,0.40E-3
     
184. *GET,VCE(3),VARI,10,RTIME,0.88E-3
     
185. *GET,VCE(4),VARI,10,RTIME,0.13E-2
     
186. *GET,VCE(5),VARI,10,RTIME,0.186E-2
     
187. 
     
188. *GET,VCA(1),VARI,2,RTIME,0.18E-3
     
189. *GET,VCA(2),VARI,2,RTIME,0.40E-3
     
190. *GET,VCA(3),VARI,2,RTIME,0.88E-3
     
191. *GET,VCA(4),VARI,2,RTIME,0.13E-2
     
192. *GET,VCA(5),VARI,2,RTIME,0.186E-2
     
193. 
     
194. *DIM,LABEL,CHAR,5
     
195. *DIM,VALUE,,5,3
     
196. LABEL(1) = '1.8E-2s','4.0E-2s','8.8E-2s','1.3E-1s','1.86E-1s'
     
197. *VFILL,VALUE(1,1),DATA,VCE(1),VCE(2),VCE(3),VCE(4),VCE(5)
     
198. *VFILL,VALUE(1,2),DATA,VCA(1),VCA(2),VCA(3),VCA(4),VCA(5)
     
199. *VFILL,VALUE(1,3),DATA,ABS(VCA(1)/VCE(1)),ABS(VCA(2)/VCE(2)),ABS(VCA(3)/VCE(3)),
     
200. ABS(VCA(4)/VCE(4)),ABS(VCA(5)/VCE(5))
     
201. /OUT,vm237,vrt
     
202. /COM
     
203. /COM,------------------- VM237 RESULTS COMPARISON ---------------------
     
204. /COM,
     
205. /COM, | TARGET | ANSYS | RATIO
     
206. /COM,
     
207. /COM, VC for t @ ...
     
208. /COM
     
209. *VWRITE,LABEL(1),VALUE(1,1),VALUE(1,2),VALUE(1,3)
     
210. (1X,A8,' ',F12.4,' ',F12.4,' ',1F5.3)
     
211. /COM,----------------------------------------------------------------
     
212. /OUT
     
213. FINI

复制代码

雪缘

另外可以help VM175 VM176看一下

glaciery

系统帮助文件里面有，最好学习命令流，这样跟方便交流，网上很多地方有这方面的例子！上网多找找吧

lixinxin5833