! Test für KONTAKT ! =================================================================================== ! Beschreibung: Blech auf Blech, dazwischen Kontakt; MAT: Linear ! D_0: Abstand zwischen Blechen ! oberes Blech an weichen federn aufgehängt, unteres fest ! Belastung: Grav. um Kontakt zu schließen ! Belastung mit Flaechendruck Druck + ZUG ! Einheiten: m, kN, t, s ! =================================================================================== finish $ /clear,nostart /PREP7 ! ----- Parameter ------------------------------------------------------------------- b1_x=2 $ b1_y=1.0 $ b1_z=0.010 ! Blech 1: Abmessungen e1_x=50 $ e1_y=4 ! Blech 1: NDIVs b2_x=2.4 $ b2_y=1.2 $ b2_z=0.001 ! Blech 2: Abmessungen e2_x=50 $ e2_y=6 ! Blech 2: NDIVs REIB=0.99 ! Reibung in Kontaktfuge a_1=0 !10e-10 ! Reals fuer Kontakt D_0=0.001 ! Abstand zw. Blechen k_fed=1e-5 ! Steifigkeit der Federn SWITCH_1= 0 ! 0: Randlinienlast at TIME=2 1: Trapezdruckverteilung at TIME=2 ! ----- KO-Systeme ------------------------------------------------------------------ LOCAL, 99, 0, 0, 0, D_0, , , , , ! ----- Material -------------------------------------------------------------------- ! ---------- Blech 1 --------------------------------------( Festigkeit entsp. Alu) - MP,EX ,1,0.7e8 $ MP,DENS,1,1 $ MP,NUXY,1,0.34 ! DENS NICHT Aendern ! ---------- Blech 2 ---------------------------------------------------------------- MP,EX ,2,1e20 $ MP,DENS,2,1e-10 $ MP,NUXY,2,0.3 ! ---------- Reibungskontakt -------------------------------------------------------- MP,MU ,10,REIB ! ----- Querschnitte ---------------------------------------------------------------- SECTYPE,1,SHELL $ SECDATA,b1_z ! Blech 1 SECTYPE,2,SHELL $ SECDATA,b2_z ! Blech 2 ! ----- Elemente -------------------------------------------------------------------- ET,1,181 $ KEYOPT,1,3,2 $ KEYOPT,1,8,2 ! Blech 1 ET,2,181 $ KEYOPT,2,3,2 $ KEYOPT,2,8,2 ! Blech 2 ET,5,14 $ KEYO,5,3,0 $ r,5,k_FED ! Federn ET,10,173 ! 3D-surface-to-surface CONTA173 KEYOPT,10,1,0 ! DoF ux,uy,uz KEYOPT,10,2,0 ! augment.LaGrange KEYOPT,10,4,0 ! Gausspunkt KEYOPT,10,5,4 ! AUTO ICONT KEYOPT,10,6,0 ! default stiffness updating KEYOPT,10,10,1 ! contact stiffnes update KEYOPT,10,11,0 ! exclude shell THK KEYOPT,10,12,0 ! 0: Standard, 1: rough et,11,170 ! Targetelement TARGE170 KEYOPT,11,2,0 ! KEYOPT,11,4,000000 r,10,a_1,a_1 ! ===== Modellierung ================================================================ ! ----- GEOMETRIE ------------------------------------------------------------------- ! ---------- Blech 1 CSYS,99 $ WPCSYS,1,99 BLC4, -b1_x/2, -b1_y/2, b1_x, b1_y, $ CM,BL_1,AREA lsel,s,loc,x,-b1_x/2 $ LSEL,a,loc,x,b1_x/2 $ LESI, all, , , e1_y LSEL,INVE $ LESI, all, , , e1_x $ alls AATT,1,,1,,1 $ ASEL,s,area,,BL_1 $ AMESH,all ! ---------- Blech 2 CSYS,0 $ WPCSYS,1,0 BLC4, -b2_x/2, -b2_y/2, b2_x, b2_y, $ asel,u,area,, BL_1 $ CM,BL_2,AREA lsel,s,loc,x,-b2_x/2 $ LSEL,a,loc,x,b2_x/2 $ Lsel,r,line,,BL_2 $ LESI, all, , , e2_y LSEL,INVE $ Lsel,r,line,,BL_2 $ LESI, all, , , e2_x alls ASEL,s,area,,BL_2 $ AATT,2,,2,,2 ASEL,s,area,,BL_2 $ AMESH,all ! ----- Auflager: ------------------------------------------------------------------- ! ---------- Volleinspannung Blech 2 an allen Knoten (Rigid) ------------------------ ASEL,s,area,,BL_2 $ NSLA,S,1 $ d,all,all $ alls ! ---------- Federlager ------------------------------------------------------------- CSYS,0 $ WPCSYS,1,0 *DO,i_nx,-1,1,2 *DO,i_ny,-1,1,2 n,,(i_nx)*b1_x/2,(i_ny)*b1_y ,D_0+1 $ FED_LAG_temp_1=node((i_nx)*b1_x/2,(i_ny)*b1_y ,D_0+1) n,,(i_nx)*b1_x ,(i_ny)*b1_y/2,D_0+1 $ FED_LAG_temp_2=node((i_nx)*b1_x ,(i_ny)*b1_y/2,D_0+1) ESEL,s,TYPE,,1 $ NSLE,s NSEL,r,loc,x,(i_nx)*b1_x/2 $ NSEL,r,loc,y,(i_ny)*b1_y/2 *GET,ECK_N_TEMP,NODE,,num,min TYPE,5 $ REAL,5 e,FED_LAG_temp_1,ECK_N_TEMP $ e,FED_LAG_temp_2,ECK_N_TEMP alls *ENDDO *ENDDO NSEL,s,loc,z,D_0+1 $ d,ALL,ALL $ ALLS ! ----- Kontaktdefinition ----------------------------------------------------------- MAT,10 $ REAL,10 $ TYPE,10 $ TSHAP,QUAD ESEL,s,elem,,BL_1 $ ESURF,,BOTTOM $ ALLS MAT,10 $ REAL,10 $ TYPE,11 $ TSHAP,QUAD ESEL,s,elem,,BL_2 $ ESURF,,TOP $ ALLS ! ---------- Lastordintaten G_TOTAL=b1_x*b1_y*b1_z*10 ! Gewichtskraft Blech 1 p_G=G_TOTAL/(b1_x*b1_y) ! Anpressdruck inf Gewicht ! ===== SOLU ======================================================================== ! --- A ---- Set-up /solu $ antype,static,new $ csys,0 $ alls AUTOTS,on $ TIMINT,off $ solcon,on,on ! $ NLGEOM,on ! --- B ---- Vorbelastung vertikal via acel+dens von Blech_1, um kontakt zu finden TIME,1 ACEL,0,0,10 solve ! Anmerkung: post1 -> contact pressure ~= p_G -> ok ! --- C ---- Abhebe-Drucklast via tab. BC -p_r at -b1_x/2 bis p_r at ---------------- /SOLU Time,2 FAKT=1.4 ! ---------- Lastdefinition *DIM,DRUCK,TABLE,2,2,1,y,x,,0 *VFILL,DRUCK(1,0),dATA,-b1_y/2,b1_y/2 ! INDEXSPALTE DRUCK(0,1)=-b1_x/2 $ DRUCK(0,2)=b1_x/2 ! INDEXZEILE DRUCK(1,1 )=+p_G*FAKT $ DRUCK(1,2)=-p_G*FAKT ! DATEN DRUCK(2,1 )=+p_G*FAKT $ DRUCK(2,2)=-p_G*FAKT ! DATEN /PSF,PRES,NORM,2,1,1 sfe,BL_1,1,pres,1,%DRUCK% $ alls /EOF SOLVE ! ===== POST ======================================================================== /post1 plnsol,u,z