Gasket해석 방법과 bolt해석 방법

앞으로 좋은 자료가 있으면 많이 공개하도록 하겠습니다. .......................................................................... 이번에 올리는 파일은 좀 난해할 수도 있습니다. 내용을 간단히 설명드리면, 우선 볼트해석방법을 이해할 수 있을 것입니다. 반드시 ansys help를 참조해야만 내용이 이해될 것으로 생각됩니다. 그리고, ansys에서 제공하는 gasket요소를 사용하여, gasket을 구현하였습니다. 이 역시 ansys help를 참조하셔야 합니다. input파일이 다소 복잡한 경향이 있습니다만, 비선형해석이 다 그런 것이러니 하시고, 전체적인 프로세스를 중심으로 먼저 보신 다음, 상세한 적용방법을 보시면 무난하게 분석할 수 있을 것으로 봅니다. 난이도는 중상정도라 생각합니다. 생소한 내용이나 명령어가 들어 있어 제가 아래에 keyword를 정리해 놓겠습니다. help에서 keyword를 검색하여 공부하시면 도움이 많이 될 것입니다. 참, bolt 해석에서 bolt가 두개 있는데, 순차적으로 하중이 적용됩니다. 아주 재미있을 것입니다. bolt가 10개가 있는 경우도 순서를 정해주면 순서대로 bolt를 해석할 수 있습니다. 응용이 될 분야가 많겠지요. 예를들어서 bolt를 어떤 순서로 조이느냐에 따라 발생하는 잔류응력이 달라지겠으며, contact pressure 또한 달라지므로 누수문제를 해결하는 방법이 되기도 할 것입니다. 그럼 input파일을 볼까요.. 1) keyword:INTER195, GASKET, GASKET JOINTS SIMULATION, imesh, psmesh, sload 2) input file 아래의 input 파일은 5개의 load step으로 구성되어 있습니다 load step 1은 bolt를 조이는 해석이고, load step2 도 bolt를 조이는 해석입니다. load step 3은 온도를 180도로 올리고, 내압1.2MPadmf 적용하여 해석하는 단계이고, load step 4는 압력을 0으로 온도를 -10도로 내린 상태에서의 해석입니다. load step 5는 온도를 상온(20도)으로 올린 상태에서의 해석입니다. GASKET의 material을 정의할 때 조심해야합니다. input파일 참조 finish /clear /prep7 !!!!!!!!!!!!!!!!!!!!!!!!!!!! modeling k,1,93/2 k,2,229/2 k,3,229/2,23.9 k,4,93/2,23.9 a,1,2,3,4 LWPL,-1,1,0 wpoff,,,135/2-93/2 asbw,all wpoff,,,157/2-135/2 asbw,all wpoff,1.6 wprot,,,90 asbw,all kwpave, 4 wpoff,,,76.2 csys,4 k, wpoff,-6 k, l,6,14 l,14,13 l,13,4 al,9,2,4,7 wpoff,,,60 k, wpoff,6 k, l,14,15 l,15,16 l,16,13 al,4,10,22,23 adel,2,,,1 kwpave,1 wpoff,,,-3 arsym,z,all KWPAVE, 20 wpoff,,,-25 k, wpoff,-21 k, L, 20, 31 L, 31, 32 L, 32, 21 al,27,43,44,45 wpoff,,,-100 k, wpoff,21 k, L, 31, 34 L, 34, 33 L, 33, 32 al,44,46,47,48 ADELE, 3, , ,1 LWPL,-1,12,0.5 asbw,all LWPL,-1,9,0.15 asbw,all L, 38, 15 L, 27, 5 L, 46, 43 FLST,2,4,4 FITEM,2,57 FITEM,2,1 FITEM,2,51 FITEM,2,4 AL,P51X FLST,2,4,4 FITEM,2,75 FITEM,2,4 FITEM,2,67 FITEM,2,8 AL,P51X KWPAVE, 4 wpoff,,,48.75 asbw,all WPCSYS,-1,0 k, k,,,-100 VROTAT,all, , , , , ,55,54,22.5, , VSYMM,Z,all, , , ,0,0 nummrg,kp KWPAVE, 50 wprot,,90 CYL4, , ,14.5/2 FLST,2,2,4,ORDE,2 FITEM,2,225 FITEM,2,228 LCOMB,P51X, ,0 FLST,2,2,4,ORDE,2 FITEM,2,231 FITEM,2,233 LCOMB,P51X, ,0 VOFFST,153,70 FLST,2,5,6,ORDE,5 FITEM,2,4 FITEM,2,-5 FITEM,2,36 FITEM,2,-37 FITEM,2,65 VOVLAP,P51X FLST,2,5,6,ORDE,5 FITEM,2,6 FITEM,2,9 FITEM,2,38 FITEM,2,41 FITEM,2,70 VOVLAP,P51X FLST,2,6,6,ORDE,3 FITEM,2,65 FITEM,2,-69 FITEM,2,79 VDELE,P51X, , ,1 FLST,2,4,6,ORDE,4 FITEM,2,4 FITEM,2,-5 FITEM,2,36 FITEM,2,-37 VDELE,P51X, , ,1 FLST,5,32,6,ORDE,10 FITEM,5,33 FITEM,5,-35 FITEM,5,39 FITEM,5,-40 FITEM,5,42 FITEM,5,-64 FITEM,5,73 FITEM,5,-74 FITEM,5,77 FITEM,5,-78 VSEL,S, , ,P51X vdel,all,,,1 allsel WPCSYS,-1,0 wprot,,,22.5 csys,4 FLST,3,32,6,ORDE,10 FITEM,3,1 FITEM,3,-3 FITEM,3,7 FITEM,3,-8 FITEM,3,10 FITEM,3,-32 FITEM,3,71 FITEM,3,-72 FITEM,3,75 FITEM,3,-76 VSYMM,Z,P51X, , , ,0,0 nummrg,kp WPCSYS,-1,0 FLST,2,2,3,ORDE,2 FITEM,2,109 FITEM,2,114 KWPAVE,P51X wprot,,90 CYL4, , ,13/2 CYL4, , ,27/2 FLST,2,2,5,ORDE,2 FITEM,2,49 FITEM,2,59 AOVLAP,P51X LWPL,-1,222,0 FLST,2,2,5,ORDE,2 FITEM,2,49 FITEM,2,152 ASBW,P51X, , FLST,2,2,5,ORDE,2 FITEM,2,59 FITEM,2,164 ADELE,P51X, , ,1 csys,4 FLST,2,2,5,ORDE,2 FITEM,2,158 FITEM,2,169 VEXT,P51X, , ,0,-8,0,,,, FLST,2,2,5,ORDE,2 FITEM,2,158 FITEM,2,169 VEXT,P51X, , ,0,53.8,0,,,, FLST,2,2,5,ORDE,2 FITEM,2,200 FITEM,2,215 VEXT,P51X, , ,0,8,0,,,, vdel,64,,,1 FLST,5,5,6,ORDE,4 FITEM,5,61 FITEM,5,-63 FITEM,5,65 FITEM,5,-66 VSEL,S, , ,P51X cm,bolt,volu allsel cmsel,u,bolt WPCSYS,-1,0 wprot,,,10 vsbw,all wprot,,,25 vsbw,all allsel cmsel,s,bolt wprot,,,-35 wprot,,,22.5 csys,4 VSYMM,Z,all cm,bolt_all,volu allsel WPCSYS,-1,0 wprot,,90 wprot,,,90 csys,4 CSWPLA,11,1,1,1, VGEN, ,all, , , ,90, , , ,1 csys,0 /VIEW, 1 ,,-1 /rep !!!!!!!!!!!!!!!!!!!!!!!! Define Element Type et,1,185 !! flange et,2,185 !! bolt et,3,195 !! gasket et,4,170 et,5,173 !!!!!!!!!!!!!!!!!!!!!!!! Define real Constants r,(1:4), !!!!!!!!!!!!!!!!!!!!!!!! Material Properties ! material 1 - flange ! material 2 - bolts ! material 3 - gasket body mp,ex ,(1:2),210e3 mp,prxy,(1:2),0.3 mp,alpx,(1:2),15e-6 mp,reft,(1:2),20 mp,alpx,3,50e-6 mp,reft,3,20 mp,mu,5,0.2 !!!!!!!!!!!!!!!!!! material 3 (gasket body) tb,gasket,3,1,,para tbdata,1,1/11,5,1e-5 tb,gasket,3,1,6,comp tbpt,,2.700000000000e-002,2.080000000000e+000 tbpt,,5.400000000000e-002,8.320000000000e+000 tbpt,,8.100000000000e-002,1.872000000000e+001 tbpt,,1.080000000000e-001,3.328000000000e+001 tbpt,,1.350000000000e-001,5.200000000000e+001 tbpt,,1.750000000000e-001,5.600000000000e+001 tb,gasket,3,1,7,nunl tbpt,,1.750000000000e-001,5.600000000000e+001 tbpt,,1.675000000000e-001,4.536000000000e+001 tbpt,,1.600000000000e-001,3.584000000000e+001 tbpt,,1.525000000000e-001,2.744000000000e+001 tbpt,,1.375000000000e-001,1.400000000000e+001 tbpt,,1.225000000000e-001,5.040000000000e+000 tbpt,,1.000000000000e-001,0.000000000000e+000 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! FLST,5,76,6,ORDE,34 FITEM,5,9 FITEM,5,12 FITEM,5,14 FITEM,5,-15 FITEM,5,21 FITEM,5,-24 FITEM,5,28 FITEM,5,-33 FITEM,5,65 FITEM,5,-66 FITEM,5,70 FITEM,5,-72 FITEM,5,75 FITEM,5,-76 FITEM,5,81 FITEM,5,-84 FITEM,5,88 FITEM,5,-98 FITEM,5,103 FITEM,5,105 FITEM,5,-106 FITEM,5,112 FITEM,5,-115 FITEM,5,119 FITEM,5,-131 FITEM,5,134 FITEM,5,-135 FITEM,5,138 FITEM,5,144 FITEM,5,-147 FITEM,5,151 FITEM,5,-161 FITEM,5,164 FITEM,5,-165 VSEL,S, , ,P51X cmsel,u,bolt_all aslv,s,1 lsla,s,1 FLST,5,45,4,ORDE,45 FITEM,5,5 FITEM,5,25 FITEM,5,-26 FITEM,5,28 FITEM,5,42 FITEM,5,46 FITEM,5,48 FITEM,5,73 FITEM,5,85 FITEM,5,134 FITEM,5,136 FITEM,5,184 FITEM,5,189 FITEM,5,-190 FITEM,5,193 FITEM,5,214 FITEM,5,216 FITEM,5,219 FITEM,5,262 FITEM,5,267 FITEM,5,364 FITEM,5,372 FITEM,5,374 FITEM,5,377 FITEM,5,412 FITEM,5,417 FITEM,5,420 FITEM,5,435 FITEM,5,450 FITEM,5,510 FITEM,5,514 FITEM,5,517 FITEM,5,519 FITEM,5,540 FITEM,5,543 FITEM,5,545 FITEM,5,621 FITEM,5,624 FITEM,5,670 FITEM,5,675 FITEM,5,-676 FITEM,5,678 FITEM,5,699 FITEM,5,702 FITEM,5,704 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,15,4,ORDE,15 FITEM,5,43 FITEM,5,45 FITEM,5,82 FITEM,5,206 FITEM,5,208 FITEM,5,212 FITEM,5,399 FITEM,5,404 FITEM,5,409 FITEM,5,532 FITEM,5,535 FITEM,5,537 FITEM,5,691 FITEM,5,694 FITEM,5,697 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,2, , , , ,1 FLST,5,62,4,ORDE,43 FITEM,5,89 FITEM,5,-90 FITEM,5,104 FITEM,5,192 FITEM,5,195 FITEM,5,202 FITEM,5,209 FITEM,5,474 FITEM,5,-475 FITEM,5,547 FITEM,5,563 FITEM,5,-564 FITEM,5,566 FITEM,5,568 FITEM,5,575 FITEM,5,-579 FITEM,5,581 FITEM,5,590 FITEM,5,-607 FITEM,5,617 FITEM,5,-618 FITEM,5,620 FITEM,5,622 FITEM,5,625 FITEM,5,627 FITEM,5,650 FITEM,5,652 FITEM,5,655 FITEM,5,661 FITEM,5,667 FITEM,5,669 FITEM,5,672 FITEM,5,-673 FITEM,5,677 FITEM,5,680 FITEM,5,683 FITEM,5,685 FITEM,5,690 FITEM,5,692 FITEM,5,695 FITEM,5,698 FITEM,5,700 FITEM,5,703 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,2, , , , ,1 FLST,5,39,4,ORDE,39 FITEM,5,29 FITEM,5,31 FITEM,5,37 FITEM,5,-38 FITEM,5,56 FITEM,5,80 FITEM,5,117 FITEM,5,119 FITEM,5,142 FITEM,5,167 FITEM,5,175 FITEM,5,198 FITEM,5,204 FITEM,5,-205 FITEM,5,210 FITEM,5,225 FITEM,5,231 FITEM,5,319 FITEM,5,325 FITEM,5,338 FITEM,5,-339 FITEM,5,394 FITEM,5,396 FITEM,5,445 FITEM,5,447 FITEM,5,453 FITEM,5,461 FITEM,5,490 FITEM,5,495 FITEM,5,500 FITEM,5,529 FITEM,5,-530 FITEM,5,565 FITEM,5,569 FITEM,5,628 FITEM,5,653 FITEM,5,659 FITEM,5,688 FITEM,5,-689 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,2, , , , ,1 esiz,4 type,1 mat,1 vsweep,all !!!!!!!!!!!!!!!!!!!!!!!!!!!!! FLST,5,8,6,ORDE,8 FITEM,5,1 FITEM,5,-2 FITEM,5,7 FITEM,5,11 FITEM,5,64 FITEM,5,67 FITEM,5,77 FITEM,5,99 VSEL,s, , ,P51X aslv,s,1 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! Gasket Part Mesh type,3 mat,3 imesh,area,298,297 imesh,area,402,404 imesh,area,39,42 imesh,area,80,82 imesh,area,232,228 imesh,area,331,333 imesh,area,5,12 imesh,area,61,63 allsel !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! vsel,u,mat,,1 vsel,u,mat,,3 cmsel,u,bolt_all aslv,s,1 lsla,s,1 FLST,5,15,4,ORDE,15 FITEM,5,2 FITEM,5,7 FITEM,5,21 FITEM,5,126 FITEM,5,128 FITEM,5,147 FITEM,5,249 FITEM,5,254 FITEM,5,274 FITEM,5,313 FITEM,5,405 FITEM,5,418 FITEM,5,459 FITEM,5,616 FITEM,5,633 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,3, , , , ,1 !* FLST,5,36,4,ORDE,36 FITEM,5,14 FITEM,5,-15 FITEM,5,19 FITEM,5,-20 FITEM,5,50 FITEM,5,69 FITEM,5,107 FITEM,5,114 FITEM,5,159 FITEM,5,163 FITEM,5,165 FITEM,5,179 FITEM,5,181 FITEM,5,-182 FITEM,5,259 FITEM,5,281 FITEM,5,283 FITEM,5,294 FITEM,5,300 FITEM,5,312 FITEM,5,350 FITEM,5,353 FITEM,5,423 FITEM,5,425 FITEM,5,440 FITEM,5,479 FITEM,5,482 FITEM,5,486 FITEM,5,506 FITEM,5,-507 FITEM,5,554 FITEM,5,558 FITEM,5,643 FITEM,5,649 FITEM,5,665 FITEM,5,-666 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,2, , , , ,1 FLST,5,43,4,ORDE,33 FITEM,5,92 FITEM,5,97 FITEM,5,110 FITEM,5,122 FITEM,5,130 FITEM,5,132 FITEM,5,149 FITEM,5,-150 FITEM,5,156 FITEM,5,160 FITEM,5,170 FITEM,5,217 FITEM,5,-218 FITEM,5,310 FITEM,5,356 FITEM,5,476 FITEM,5,-477 FITEM,5,548 FITEM,5,-553 FITEM,5,555 FITEM,5,557 FITEM,5,560 FITEM,5,571 FITEM,5,-574 FITEM,5,584 FITEM,5,-589 FITEM,5,630 FITEM,5,632 FITEM,5,635 FITEM,5,637 FITEM,5,639 FITEM,5,642 FITEM,5,645 CM,_Y,LINE LSEL, , , ,P51X CM,_Y1,LINE CMSEL,,_Y !* LESIZE,_Y1, , ,2, , , , ,1 type,1 mat,1 vsweep,all allsel cmsel,s,bolt_all esiz,6 type,2 mat,2 vsweep,all allsel !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!save,gasket_2bolt,db !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! FLST,5,2,5,ORDE,2 FITEM,5,493 FITEM,5,496 ASEL,S, , ,P51X nsla,s,1 cm,con1_up,node asel,s,,,169 nsla,s,1 cm,tar1_up,node FLST,5,2,5,ORDE,2 FITEM,5,604 FITEM,5,606 ASEL,S, , ,P51X nsla,s,1 cm,con2_up,node asel,s,,,54 nsla,s,1 cm,tar2_up,node FLST,5,2,5,ORDE,2 FITEM,5,488 FITEM,5,490 ASEL,S, , ,P51X nsla,s,1 cm,con1_low,node asel,s,,,215 nsla,s,1 cm,tar1_low,node FLST,5,2,5,ORDE,2 FITEM,5,611 FITEM,5,-612 ASEL,S, , ,P51X nsla,s,1 cm,con2_low,node asel,s,,,165 nsla,s,1 cm,tar2_low,node allsel type,4 real,1 mat,5 cmsel,s,tar1_up esurf type,5 cmsel,s,con1_up esurf type,4 real,2 mat,5 cmsel,s,tar1_low esurf type,5 cmsel,s,con1_low esurf type,4 real,3 mat,5 cmsel,s,tar2_up esurf type,5 cmsel,s,con2_up esurf type,4 real,4 mat,5 cmsel,s,tar2_low esurf type,5 cmsel,s,con2_low esurf allsel !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !save,gasket_2bolt_2,db !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! WPCSYS,-1,0 csys,4 CSWPLA,12,1 asel,s,loc,x,93/2-0.5,93/2+0.5 cm,in_pres,area allsel vsel,s,,,63 eslv,s,1 cm,bolt1,elem vsel,s,,,6 eslv,s,1 cm,bolt2,elem allsel !!!!!!!!!!!!!! create pretension csys,0 cmsel,s,bolt1 psmesh,1,BOLT11,,all,,0,z,0.1 esel,r,ename,,179 PRETELEM=ELNEXT(0) *get,MYBOLT1,elem,pretelem,node,3 cmsel,s,bolt2 psmesh,2,BOLT22,,all,,0,z,0.1 esel,r,ename,,179 PRETELEM=ELNEXT(0) *get,MYBOLT2,elem,pretelem,node,3 cm,lines,line cmplot allsel !!!!!!!!!!!!!!!!!! apply b.c. WPCSYS,-1,0 csys,4 CSWPLA,11,1 asel,s,loc,y,-0.2,0.2 DA,all,SYMM asel,s,loc,y,44.8,45.2 DA,all,SYMM allsel csys,4 KWPAVE,33 asel,s,loc,z nsla,s,1 csys,11 nrot,all d,all,uy d,all,uz allsel csys,0 /solu nlgeom,on rescon,define,none pred,off tunif,20 allsel,all outres,all,all !!!!!!!!!!!!!!!!!!!!!!! load step 1: apply pretension nsubst,10,100,10 sload,1,9,TINY,forc,2500,1,2 sload,2,9,TINY,forc,2500,2,3 solve !!!!!!!!!!!!!!!!!!!!!!! load step 2: apply pretension solve !!!!!!!!!!!!!!!!!! load step 3: heat up and Internal Pressure nsubst,10,1e2,10 cmsel,s,in_pres sfa,all,,pres,1.2 allsel,all tunif,180 solve !!!!!!!!!!!!!!!!!!!!!!! load step 4: cool down tunif,-10 solve !!!!!!!!!!!!!!!!!!!!!! load step 5: ambient cmsel,s,in_pres sfa,all,,pres,0 allsel,all tunif,20 solve