声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1599|回复: 3

[综合讨论] 模型中的裂纹在画网格后为什么消失了

[复制链接]
发表于 2010-11-21 13:39 | 显示全部楼层 |阅读模式

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

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

x
画了一个波形板,有一些裂纹,但在画完网格后这些裂纹消失了,请问哪位大侠知道这是什么情况
原模型
C:\Documents and Settings\Administrator\5.jpg
画网格后
C:\Documents and Settings\Administrator\7.jpg
回复
分享到:

使用道具 举报

发表于 2010-11-22 10:19 | 显示全部楼层
裂纹问题没有做过,不太清楚,不过你可以参考相应的例子比对一下,比如VM143就是裂纹的例子

另外在附一个裂纹的分析实例
  1. !12.2  均布拉力作用下含裂纹板的应力强度因子计算
  2. ! 本程序来源于邢静忠等编著的《ANSYS分析实例与工程应用》,机械工业出版社(2004年)
  3. ! (1) 定义将二维SOLID45模型通过添加1/4位置的中间节点转化成SOLID95裂纹尖端单元的宏程序
  4. *CREATE,FRACT,MAC
  5. !  MACRO TO CREATE 3D SOLID95 CRACK TIP ELEMENTS FROM 3D SOLID45 ELEMENTS
  6. !     MAKE A COMPONENT CONTAINING THE CRACK TIP NODES (CRACKTIP)
  7. !     THE CRACK TIP IS BETWEEN NODES K AND O
  8. !     SET ELEMENT TYPE TO POINT TO SOLID95
  9. !     SET ARG1 TO  N  (THE TYPE OF THE ELEMENTS AROUND THE CRACK TIP)
  10. /NOPR
  11. NSEL,ALL
  12. *GET,N,NODE,,NUM,MAX
  13. CMSEL,S,CRACKTIP
  14. ESLN
  15. *GET,ELMAX,ELEM,,NUM,MAX
  16. *DO,IEL,1,ELMAX
  17.    ELMI=IEL
  18.    *IF,ELMI,LE,0,EXIT
  19.    *GET,ELTYPE,ELEM,ELMI,ATTR,TYPE
  20.    *IF,ELTYPE,NE,ARG1,CYCLE
  21.    N3 = NELEM(ELMI,3)
  22.    *IF,NSEL(N3),LE,0,CYCLE
  23.    N7 = NELEM(ELMI,7)
  24.    *IF,NSEL(N7),LE,0,CYCLE
  25.    N1 = NELEM(ELMI,1)
  26.    N2 = NELEM(ELMI,2)
  27.    N5 = NELEM(ELMI,5)
  28.    N6 = NELEM(ELMI,6)
  29.    X3 = 0.75*NX(N3)
  30.    Y3 = 0.75*NY(N3)
  31.    Z3 = 0.75*NZ(N3)
  32.    X  = 0.25*NX(N2) + X3
  33.    Y  = 0.25*NY(N2) + Y3
  34.    Z  = 0.25*NZ(N2) + Z3
  35.    N  = N + 1
  36.    N10 = N
  37.    N,N10,X,Y,Z
  38.    X  = 0.25*NX(N1) + X3
  39.    Y  = 0.25*NY(N1) + Y3
  40.    Z  = 0.25*NZ(N1) + Z3
  41.    N  = N + 1
  42.    N12= N
  43.    N,N12,X,Y,Z
  44.    X7 = 0.75*NX(N7)
  45.    Y7 = 0.75*NY(N7)
  46.    Z7 = 0.75*NZ(N7)
  47.    X  = 0.25*NX(N6) + X7
  48.    Y  = 0.25*NY(N6) + Y7
  49.    Z  = 0.25*NZ(N6) + Z7
  50.    N  = N + 1
  51.    N14 = N
  52.    N,N14,X,Y,Z
  53.    X  = 0.25*NX(N5) + X7
  54.    Y  = 0.25*NY(N5) + Y7
  55.    Z  = 0.25*NZ(N5) + Z7
  56.    N  = N + 1
  57.    N16 = N
  58.    N,N16,X,Y,Z
  59.    N4=N3
  60.    N8=N7
  61.    NSEL,ALL
  62.    TYPE,3
  63.    EN,ELMI,N1,N2,N3,N4,N5,N6,N7,N8
  64.    EMORE,0,N10,0,N12,0,N14,0,N16
  65.    EMORE,
  66. *ENDDO
  67. CMSEL,U,CRACKTIP
  68. NUMMRG,NODE
  69. NSEL,ALL
  70. ESEL,ALL
  71. /GOPR
  72. *END
  73. ! (2) 在前处理模块设置工程选项,分析类型,单元类型和材料参数
  74. FINISH
  75. /CLEAR,NOSTART
  76. /FILENAME,EX12.2
  77. /PREP7
  78. SMRT,OFF
  79. /TITLE, EX12.2, STRESS INTENSITY - CRACK IN FINITE WIDTH PLATE(SOLID95)
  80. C***    BROWN AND SRAWLEY, ASTM SPECIAL TECHNICAL PUBLICATION NO. 410.
  81. /COM,     ****** CRACK IN 3-DIMENSIONS USING SOLID45 AND SOLID95
  82. ANTYPE,STATIC
  83. ET,1,SOLID45
  84. ET,2,SOLID45
  85. ET,3,SOLID95
  86. MP,EX,1,207E3
  87. MP,NUXY,1,0.3
  88. ! (3) 在柱坐标系下定义节点
  89. CSYS,1
  90. N,1
  91. NGEN,9,20,1
  92. N,11,20.32
  93. N,171,20.32,180
  94. FILL,11,171,7,31,20
  95. CSYS,0
  96. FILL,1,11,9,2,1,9,20,3     为3
  97. N,15,101.6
  98. N,75,101.6,127
  99. FILL,15,75,2,35,20
  100. N,155,-25.4,127
  101. FILL,75,155,3,95,20
  102. N,172,-25.4
  103. FILL,155,172,5,177,-1,,,0.15
  104. FILL,11,15,3,,,7,20,3
  105. NGEN,2,200,1,177,,,,6.35
  106. ! (4) 定义单元连接
  107. /View,1,1,1,1
  108. /PNUM,NODE,1
  109. E,2,22,1,1,202,222,201,201
  110. EGEN,8,20,-1
  111. E,2,3,23,22,202,203,223,222
  112. EGEN,8,20,-1
  113. EGEN,9,1,-8
  114. EGEN,5,1,73,78
  115. E,171,151,173,172,371,351,373,372
  116. E,151,131,174,173,351,331,374,373
  117. E,131,132,175,174,331,332,375,374
  118. EGEN,3,1,-1
  119. E,134,135,155,177,334,335,355,377
  120. TYPE,2
  121. EMODIF,1
  122. *REPEAT,8,1
  123. NUMMRG,NODE
  124. ! (5) 定义位移约束
  125. NSEL,S,LOC,X,0
  126. NSEL,R,LOC,Y,0
  127. CM,CRACKTIP,NODE
  128. /NERR,0
  129. FRACT,2
  130. /NERR,DEFA
  131. /OUTPUT
  132. OUTPR,,ALL
  133. OUTPR,VENG,ALL
  134. NSEL,S,LOC,X,-25.4
  135. DSYM,SYMM,X
  136. NSEL,S,LOC,X,0,101.6
  137. NSEL,R,LOC,Y,0
  138. DSYM,SYMM,Y
  139. NSEL,ALL
  140. D,ALL,UZ
  141. NSEL,S,LOC,Y,127
  142. SF,ALL,PRES,-0.0038902
  143. NSEL,ALL
  144. ESEL,ALL
  145. /VIEW,1,0,0,1
  146. EPLOT
  147. FINISH
  148. ! (6) 在求解模块开始求解
  149. /OUTPUT,SCRATCH
  150. /SOLU
  151. SOLVE
  152. FINISH
  153. /OUTPUT
  154. ! (7) 在后处理模块中,定义路径后计算应力强度因子
  155. /POST1
  156. ETABLE,SENE,SENE
  157. ETABLE,VOLU,VOLU
  158. C*** IN POST1 DETERMINE KI (STRESS INTENSITY FACTOR) USING KCALC
  159. PATH,KI1,3,,48
  160. PPATH,1,1
  161. PPATH,2,406
  162. PPATH,3,162
  163. KCALC,,,1
  164. *GET,KI1,KCALC,,K,1
  165. ! (8) 定义完成J积分计算的宏程序JIN1
  166. ! *********************************************************************
  167. ! **************************  J-INTEGRAL USER FILE ******************
  168. ! *******************************************************************
  169. ! ***NOTE:- IN GENERAL USAGE, THE USER FILE WOULD BE AVAILABLE IN THE
  170. !            LOCAL DIRECTORY RATHER THAN BEING CREATED IN THE INPUT
  171. ! *******************************************************************
  172. *CREATE,JIN1
  173. STINFC
  174. SEXP,W,SENE,VOLU,1,-1
  175. PATH,JINT,4,50,48
  176. PPATH,1,ARG1
  177. PPATH,2,ARG2
  178. PPATH,3,ARG3
  179. PPATH,4,ARG4
  180. PDEF,W,ETAB,W
  181. PCALC,INTG,J,W,YG
  182. *GET,JA,PATH,,LAST,J
  183. PDEF,CLEAR
  184. PVECT,NORM,NX,NY,NZ
  185. PDEF,INTR,SX,SX
  186. PDEF,INTR,SY,SY
  187. PDEF,INTR,SXY,SXY
  188. PCALC,MULT,TX,SX,NX
  189. PCALC,MULT,C1,SXY,NY
  190. PCALC,ADD,TX,TX,C1
  191. PCALC,MULT,TY,SXY,NX
  192. PCALC,MULT,C1,SY,NY
  193. PCALC,ADD,TY,TY,C1
  194. *GET,DX,PATH,,LAST,S
  195. DX=DX/100
  196. PCALC,ADD,XG,XG,,,,-DX/2
  197. PDEF,INTR,UX1,UX
  198. PDEF,INTR,UY1,UY
  199. PCALC,ADD,XG,XG,,,,DX
  200. PDEF,INTR,UX2,UX
  201. PDEF,INTR,UY2,UY
  202. PCALC,ADD,XG,XG,,,,-DX/2
  203. C=(1/DX)
  204. PCALC,ADD,C1,UX2,UX1,C,-C
  205. PCALC,ADD,C2,UY2,UY1,C,-C
  206. PCALC,MULT,C1,TX,C1
  207. PCALC,MULT,C2,TY,C2
  208. PCALC,ADD,C1,C1,C2
  209. PCALC,INTG,J,C1,S
  210. *GET,JB,PATH,,LAST,J
  211. JINT=2*(JA-JB)
  212. PDEF,CLEAR
  213. *END
  214. ! (9) 在后处理模块用J积分确定应力强度因子KI
  215. C**************** IN POST1 DETERMINE KI FROM J-INTEGRAL
  216. CON1=207E3/(1-(0.3*0.3))
  217. *ULIB,JIN1
  218. *USE,STINFC,4,45,125,164
  219. KI2=SQRT(CON1*JINT)
  220. *STATUS,KI1
  221. *STATUS,KI2
  222. ! (10) 定义数组,将结果保存到这些数组
  223. *DIM,LABEL,CHAR,2,2
  224. *DIM,VALUE,,2,3
  225. LABEL(1,1) = 'BY DISP ','BY J-'
  226. LABEL(1,2) = 'EXTRP ','INT'
  227. *VFILL,VALUE(1,1),DATA,0.03561,0.03561
  228. *VFILL,VALUE(1,2),DATA,KI1,KI2
  229. *VFILL,VALUE(1,3),DATA,ABS(KI1/0.03561),ABS(KI2/0.03561)
  230. SAVE,TABLE_1
  231. FINISH
  232. ! (11) 改用PLANE82单元重新计算
  233. /CLEAR, NOSTART
  234. /FILENAME,EX12.2B
  235. /PREP7
  236. SMRT,OFF
  237. /TITLE, EX12.2B, FRACTURE MECHANIC - CRACK IN A FINITE WIDTH PLATE(PLANE82)
  238. /COM,   ****** CRACK IN 2-DIMENSIONS USING 2-D PLANE82 ******
  239. ET,1,PLANE82,,,2
  240. MP,EX,1,207E3
  241. MP,NUXY,1,0.3
  242. ! (12) 定义关键点和线,对称位移条件后将面划分为面单元
  243. K,1
  244. K,2,101.5
  245. K,3,101.6,127
  246. K,4,-25.4,127
  247. K,5,-25.4
  248. L,1,2
  249. L,2,3
  250. LESIZE,2,,,4
  251. L,3,4
  252. LESIZE,3,,,4
  253. L,4,5,
  254. LESIZE,4,,,6,.2
  255. L,5,1
  256. ESIZE,,5
  257. KSCON,1,3.81,1,8
  258. AL,1,2,3,4,5
  259. DL,1,1,SYMM
  260. DL,4,1,SYMM
  261. SFL,3,PRES,-0.0038902
  262. AMESH,1
  263. OUTPR,ALL
  264. FINISH
  265. ! (13) 进入求解模块求解
  266. /OUTPUT,SCRATCH
  267. /SOLU
  268. SOLVE
  269. FINISH
  270. /OUTPUT
  271. ! (14) 在后处理模块显示计算结果
  272. /POST1
  273. ETABLE,SENE,SENE
  274. ETABLE,VOLU,VOLU
  275. C*** IN POST1 DETERMINE KI (STRESS INTENSITY FACTOR) USING KCALC
  276. NSEL,S,LOC,Y,0
  277. NSEL,R,LOC,X,0
  278. *GET,NOD1,NODE,,NUM,MIN
  279. NSEL,A,LOC,Y
  280. NSEL,R,LOC,X,-0.127,-3.683
  281. *GET,NOD2,NODE,,NUM,MIN
  282. NSEL,A,LOC,Y
  283. NSEL,R,LOC,X,-3.683,-3.937
  284. *GET,NOD3,NODE,,NUM,MIN
  285. NSEL,ALL
  286. PATH,KI2,3,,48
  287. PPATH,1,NOD1
  288. PPATH,2,NOD2
  289. PPATH,3,NOD3
  290. KCALC,,,1
  291. *GET,KI1,KCALC,,K,1
  292. ! (15) 由J积分确定应力强度因子KI
  293. C**************** IN POST1 DETERMINE KI FROM J-INTEGRAL
  294. CSYS,1
  295. NSEL,S,LOC,X,12.7,20.32
  296. NSEL,R,LOC,Y,-1,1
  297. *GET,NOD4,NODE,,NUM,MAX
  298. NSEL,S,LOC,X,12.7,20.32
  299. NSEL,R,LOC,Y,35,55
  300. *GET,NOD5,NODE,,NUM,MAX
  301. NSEL,S,LOC,X,12.7,20.32
  302. NSEL,R,LOC,Y,120,145
  303. *GET,NOD6,NODE,,NUM,MAX
  304. NSEL,S,LOC,X,12.7,20.32
  305. NSEL,R,LOC,Y,179,181
  306. *GET,NOD7,NODE,,NUM,MIN
  307. NSEL,ALL
  308. CSYS,0
  309. *USE,STINFC,NOD4,NOD5,NOD6,NOD7
  310. CON1=207E3/(1-(0.3*0.3))
  311. KI2=SQRT(CON1*JINT)
  312. *STATUS,KI1
  313. *STATUS,KI2
  314. ! (16) 定义数组,将结果保存到这些数组
  315. *DIM,LABEL,CHAR,2,2
  316. *DIM,VALUE,,2,3
  317. LABEL(1,1) = 'BY DISP ','BY J-'
  318. LABEL(1,2) = 'EXTRP ','INT'
  319. *VFILL,VALUE(1,1),DATA,0.03561,0.03561
  320. *VFILL,VALUE(1,2),DATA,KI1,KI2
  321. *VFILL,VALUE(1,3),DATA,ABS(KI1/0.03561),ABS(KI2/0.03561)
  322. SAVE,TABLE_2
  323. ! (17) 从数据库文件中恢复计算结果后,将其输出到结果文件EX12_2
  324. RESUME,TABLE_1
  325. /COM
  326. /OUT,EX12_2,out
  327. /COM,------------------- EX12.2 RESULTS COMPARISON --------------
  328. /COM,
  329. /COM,                 |   TARGET   |   ANSYS   |   RATIO
  330. /COM,
  331. /COM,USING SOLID95 AND SOLID45 (3-D ANALYSIS)
  332. *VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
  333. (1X,A8,A8,'   ',F10.5,'  ',F10.5,'   ',1F5.3)
  334. /NOPR
  335. RESUME,TABLE_2
  336. /GOPR
  337. /COM,
  338. /COM,USING PLANE82 (2-D ANALYSIS)
  339. *VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)
  340. (1X,A8,A8,'   ',F10.5,'  ',F10.5,'   ',1F5.3)
  341. /COM,-----------------------------------------------------------
  342. /OUT
  343. FINISH
  344. *LIST,EX12_2,out
复制代码


 楼主| 发表于 2010-11-22 17:06 | 显示全部楼层
画网格前

                               
登录/注册后可看大图

画网格后

                               
登录/注册后可看大图


发表于 2010-11-22 18:48 | 显示全部楼层
本帖最后由 Chelsea 于 2010-11-22 18:48 编辑

请上传图片,粘贴的话只有你自己能够看到,别人是看不到的

我发现好多人都不会上传图片,给大家介绍一下,方法如下
上传图片.png
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-11 18:23 , Processed in 0.073113 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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