ansys 蠕变例子.doc

ansys 蠕变例子绷紧的螺钉的应力松弛现象/COM,ANSYS MEDIA REL. 8.1 (03-15-2004) REF. VERIF. MANUAL: REL. 8.1/VERIFY,VM132/PREP7/TITLE, VM132, STRESS RELAXATION OF A BOLT DUE TO CREEPC* STR. OF MATL., TIMOSHENKO, PART 2, 3RD ED., PAGE 531ANTYPE,STATICET,1,LINK1 ! SPAR ELEMENTR,1,1,(1/30000) ! INITIAL STRAINMP,EX,1,30E6TB,CREEP,1TBDATA,1,4.8E-30,7 ! CREEP PROPERTIESN,1N,2,10E,1,2BFUNIF,TEMP,900 ! UNIFORM TEMPERATURETIME,1000KBC,1D,ALL,ALL ! FIX ALL DOFSFINISH/SOLU SOLCONTROL,0NSUBST,100OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEPOUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEPSOLVEFINISH/POST26ESOL,2,1,LS,1,SIG ! STORE AXIAL STRESSPRVAR,2 ! PRINT AXIAL STRESS VS TIME*GET,T190,VARI,2,RTIME,190*GET,T420,VARI,2,RTIME,420*GET,T690,VARI,2,RTIME,690*GET,T880,VARI,2,RTIME,880*GET,T950,VARI,2,RTIME,950*status,parm*DIM,LABEL,CHAR,5,2*DIM,VALUE,5,3LABEL(1,1) = SIG ,SIG ,SIG ,SIG ,SIG LABEL(1,2) = 190 hr,420 hr,690 hr,880 hr,950 hr*VFILL,VALUE(1,1),DATA,975,950,925,910,905*VFILL,VALUE(1,2),DATA,T190,T420,T690,T880,T950V1 = ABS(T190/975) V2 = ABS(T420/950) V3 = ABS(T690/925)V4 = ABS(T880/910) V5 = ABS(T950/905) *VFILL,VALUE(1,3),DATA,V1,V2,V3,V4,V5/COM/OUT,vm132,vrt/COM,- VM132 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,A8, ,F10.0, ,F10.0, ,1F5.2)/COM,-/OUTFINISH*LIST,vm132,vrt=/COM,ANSYS MEDIA REL. 8.1 (03-15-2004) REF. VERIF. MANUAL: REL. 8.1/VERIFY,VM133/PREP7MP,PRXY,0.3/TITLE, VM133, MOTION OF A ROD DUE TO IRRADIATION INDUCED CREEPC* REFERENCE - ANY BASIC CALCULUS BOOKANTYPE,STATICET,1,BEAM23 ! PLASTIC BEAMR,1,.25,.0052083,.5 ! AREA, IZZ, HEIGHTMP,EX,1,300TB,CREEP,1TBDATA,55,0.5E-12,1E10 ! CREEP EQUATION CONSTANTS K1 AND K2TBDATA,66,5 ! SELECT IRRADIATION INDUCED CREEP EQUATIONN,1N,2,1E,1,2D,1,ALL ! FIX ONE ENDF,2,FX,.25 ! FORCE INDUCING CONSTANT STRESSFINISH/SOLU !SOLCONTROL,0BFE,1,TEMP,1,1000,1000,1000,1000 ! APPLY CONSTANT TEMPERATUREBFE,1,FLUE,1,0,0,0,0 ! APPLY ZERO FLUENCETIME,1E-8 ! NEAR ZERO TIME FOR FIRST LOAD STEPOUTPR,BASIC,1 ! PRINT BASIC ELEMENT SOLUTION OUTRES,EPCR,1 ! STORE CREEP STRAIN RESULTS FOR EVERY SUBSTEPCNVTOL,F,1E-6 ! NEAR ZERO VALUE FOR MINREF FIELDCNVTOL,M,-1 ! CONVERGENCE CRITERION BASED UPON MOMENTS IS ! REMOVED AS IT IS NOT NEEDED FOR THIS TESTSOLVE ! LOAD STEP 1NSUBST,50,500,50TIME,5OUTPR,BASIC,5BFE,1,FLUE,1,5E10,5E10,5E10,5E10 ! FINAL FLUENCES (RAMPED)SOLVE ! LOAD STEP 2FINISH/POST26ESOL,2,1,LEPCR,1,EPCR ! STORE CREEP STRAINPRVAR,2 ! PRINT STRAIN VARIATION WITH TIME*GET,T1,VARI,2,RTIME,0*GET,T2,VARI,2,RTIME,.5*GET,T3,VARI,2,RTIME,1*GET,T4,VARI,2,RTIME,5*status,parm*DIM,LABEL,CHAR,4,2*DIM,VALUE,4,3LABEL(1,1) = CRP STR ,CRP STR ,CRP STR ,CRP STR LABEL(1,2) = 0 hr,.5 hr,1 hr,5 hr*VFILL,VALUE(1,1),DATA,0,.00197,.00316,.00497*VFILL,VALUE(1,2),DATA,T1,T2,T3,T4*VFILL,VALUE(1,3),DATA,000,ABS(T2/.00197),ABS(T3/.00316),ABS(T4/.00497)/COM/OUT,vm133,vrt/COM,- VM133 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,A8, ,F10.5, ,F10.5, ,1F5.3)/COM,-/OUTFINISH*LIST,vm133,vrt=/COM,ANSYS MEDIA REL. 8.1 (03-15-2004) REF. VERIF. MANUAL: REL. 8.1/VERIFY,vm224JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS/SHOW,JPEG/TITLE,VM224, Implicit Creep Under Biaxial Load/COM, NAFEMS Fundamental Tests of Creep Behavior, Becker and Hyde/NOPR/COM,/COM, 2D CREEP TESTS WITH BIAXIAL CONSTANT LOAD, /COM, REFERENCE: TEST 10(A) FROM NAFEMS R0027. /COM,/COM, EXPECTED RESULTS:/COM, TIME | EPCRX | EPCRY /COM, - /COM, 0.0 | 0.0 | -0.0 /COM, 0.1 | 0.0427 | -0.0427 /COM, 1.0 | 0.135 | -0.135 /COM, 5.0 | 0.3019 | -0.3019 /COM, 10.0 | 0.4269 | -0.4269 /COM, 50.0 | 0.9546 | -0.9546 /COM, 100.0 | 1.35 | -1.35 /COM, 500.0 | 3.019 | -3.019 /COM, 1000.0 | 4.2691 | -4.2691 /COM,/PREP7C* PARAMETRIC INPUT FOR CREEP CONSTANTS (PRIMARY CREEP)!*SET,C1,1.5625E-14 !ASSIGN VALUE*SET,C2,5.0 !ASSIGN VALUE*SET,C3,-0.5 !ASSIGN VALUE*SET,C4,0 !ASSIGN VALUEC* TIME PARAMETER*SET,HOUR,1000 !ASSIGN VALUEC* ELASTIC CONSTANTMP,EX,1,200E3 !DEFINE YOUNGS MODULUSMP,NUXY,1,0.3 !DEFINE POISONS RATIOTUNIF,HOT !ASSIGN TEMP TO NODESTOFF,OFFS !SPECIFY TEMP RELATIVE TO ABSOLUTE VALUESTB,CREEP,1,6 !ACTIVATE DATA TABLETBDATA,1,C1,C2,C3,C4 !DEFINE DATA FOR TABLESAVE !SAVE /PREP7N,1,0,0,0N,2,100,0N,3,100,100N,4,0,100N,5,50,0N,6,100,50N,7,50,100N,8,0,50N,9,50,50ET,1,PLANE182KEYOPT,1,1,1KEYOPT,1,3,0E,1,5,9,8E,5,2,6,9E,9,6,3,7E,8,9,7,4NALLNSEL,S,LOC,XD,ALL,UXNSEL,S,LOC,YD,ALL,UYNALL/SOLUNSEL,S,LOC,X,100SF,ALL,PRES,-200NALLNSEL,S,LOC,Y,100SF,ALL,PRES,200NALLRATE, OFFDELT,1.0E-8,1.0E-9,1.0E-8TIME, 1.0E-8/OUT,SCRATCHOUTRES,ESOL,ALLSOLVE/OUT RATE, ON, ONDELT,1E-5,1E-5,100TIME,1000/OUT,SCRATCHSOLVE/OUTFINISH/POST26ESOL,2,1,EPCR,XESOL,3,1,EPCR,YPRVAR,2,3PLVAR,2,3*GET,RES1X,VARI,2,RTIME,1000*GET,RES1Y,VARI,3,RTIME,1000FINISHPARSAV,ALLRESUMEPARRES,CHANGE/PREP7N,1,0,0,0N,2,100,0N,3,100,100N,4,0,100N,5,50,0N,6,100,50N,7,50,100N,8,0,50ET,1,PLANE183KEYOPT,1,3,0E,1,2,3,4,5,6,7,8NALLNSEL,S,LOC,XD,ALL,UXNSEL,S,LOC,YD,ALL,UYNALL/SOLUNSEL,S,LOC,X,100SF,ALL,PRES,-200NALLNSEL,S,LOC,Y,100SF,ALL,PRES,200NALLRATE, OFFDELT,1.0E-8,1.0E-9,1.0E-8TIME, 1.0E-8/OUT,SCRATCHSOLVE/OUT RATE, ON, ONDELT,1E-5,1E-5,100TIME,1000/OUT,SCRATCHOUTRES,ESOL,ALLSOLVE/OUTFINISH/POST26ESOL,2,1,EPCR,XESOL,3,1,EPCR,YPRVAR,2,3PLVAR,2,3*GET,RES2X,VARI,2,RTIME,1000*GET,RES2Y,VARI,3,RTIME,1000*DIM,LABEL1,CHAR,2*DIM,VALUE1,2,3LABEL1(1) = ECRXX , ECRYY *VFILL,VALUE1(1,1),DATA,4.2691,-4.2691*VFILL,VALUE1(1,2),DATA,RES1X,RES1Y*VFILL,VALUE1(1,3),DATA,ABS(RES1X/4.2691),ABS(RES1Y/(-4.2691)*DIM,LABEL2,CHAR,2*DIM,VALUE2,2,3LABEL2(1) = ECRXX , ECRYY *VFILL,VALUE2(1,1),DATA,4.2691,-4.2691*VFILL,VALUE2(1,2),DATA,RES2X,RES2Y*VFILL,VALUE2(1,3),DATA,ABS(RES2X/4.2691),ABS(RES2Y/(-4.2691)/OUT,vm224,vrt/COM/COM,- VM224 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,/COM, PLANE182/COM,*VWRITE,LABEL1(1),VALUE1(1,1),VALUE1(1,2),VALUE1(1,3)(1X,A8, ,F7.4, ,F7.4, ,1F5.3)/COM,/COM, PLANE183/COM,*VWRITE,LABEL2(1),VALUE2(1,1),VALUE2(1,2),VALUE2(1,3)(1X,A8, ,F7.4, ,F7.4, ,1F5.3)/COM,-/OUTFINISH*LIST,vm224,vrt=/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0/VERIFY,VM200JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS/SHOW,JPEG/PREP7/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS/COM, - 2-D ANALYSIS -/COM, FE CALCULATIONS OF RESIDUAL STRESSES .,SOULES ET AL.ET,1,VISCO88,1 ! AXISYMMETRIC 2-D VISCOELASTIC ELEMENT*CREATE,MAC1 ! MACRO FOR MATERIAL PROPERTIES/COM, MATERIAL ONE IS G-11 GLASS AND MATERIAL TWO IS ALUMINA/COM, NOTE THAT ALUMINA IS AN ELASTIC MATERIAL THEREFORE IT/COM, DOES NOT HAVE VISCOELASTICITY AND STRUCTURAL RELAXATION/COM, MATERIAL PROPERTIES OF ALUMINATB,EVISC,2TBDATA,31,52.6E-7 ! THERMAL EXPANSION COEFFICENTS FOR ALUMINATBDATA,32,.119E-7 ! LOCATIONS 31-33TBDATA,33,-1.0E-11TBDATA,46,1.435E5 ! SHEAR MODULUS AT TIME = 0TBDATA,47,1.435E5 ! SHEAR MODULUS AT TIME = INFINITYTBDATA,48,3.11E5 ! BULK MODULUS AT TIME = 0TBDATA,49,3.11E5 ! BULK MODULUS AT TIME = INFINITY/COM, MATERIAL PROPERTIES OF G-11 GLASSTB,EVISC,1TBDATA,1,6.45D4 ! H/RTBDATA,2,0.53 ! VALUE OF XTBDATA,3,6 ! NO. OF MAXWELL ELEMENTS FOR STRUCTURAL RELAXATIONTBDATA,6,0.108 ! COEFFICIENTS OF THE MAXWELL ELEMENTS FORTBDATA,7,0.443 ! VOLUME DECAY (STRUCTURAL RELAXATION)TBDATA,8,0.166 ! LOCATIONS 6-11TBDATA,9,0.161TBDATA,10,0.046TBDATA,11,0.076TBDATA,16,3.00 ! RELAXATION TIMES FOR VOLUME DECAY FUNCTIONTBDATA,17,0.671 ! LOCATIONS 16-21TBDATA,18,0.247TBDATA,19,0.091TBDATA,20,0.033TBDATA,21,0.008TBDATA,26,3.43E-5 ! THERMAL EXPANSION COEFF. IN LIQUID STATETBDATA,31,64.7E-7 ! THERMAL EXPANSION COEFF. IN GLASS (SOLID) STATETBDATA,32,.02E-7 ! LOCATIONS 31-32TBDATA,36,618 ! FICTIVE TEMPERATURE LOCATIONS 36-41TBDATA,37,618TBDATA,38,618TBDATA,39,618TBDATA,40,618TBDATA,41,618TBDATA,46,2.79D4 ! SHEAR MODULUS AT TIME = 0TBDATA,47,0.0D0 ! SHEAR MODULUS AT TIME = INFINITYTBDATA,48,6.05D4 ! BULK MODULUS AT TIME = 0TBDATA,49,6.05D4 ! BULK MODULUS AT TIME = INFINITYTBDATA,50,3 ! THREE MAXWELL ELEMENTS FOR SHEAR RELAXATIONTBDATA,51,0.422 ! COEFF. FOR SHEAR RELAXATION LOCATIONS 51-53TBDATA,52,0.423TBDATA,53,0.155TBDATA,61,0.0689 ! RELAXATION TIMES FOR SHEAR RELAXATIONTBDATA,62,0.0065 ! LOCATIONS 61-63TBDATA,63,0.0001TBDATA,71,0 ! NO BULK MODULUS RELAXATION*END*USE,MAC1 ! EXECUTE MACRO FOR MATERIAL PROPERTIES:COM, CREATE FINITE ELEMENT MODELN,1,N,3,0.00025FILLN,5,0,(0.00025+0.00325)FILLNGEN,3,10,1,5,1,.001MAT,2E,1,21,23,3,11,22,13,2MAT,1E,3,23,25,5,13,24,15,4/COM, APPLY BOUNDARY CONDITIONS AND COUPLINGNSEL,S,LOC,YDSYM,SYMM,YNSEL,S,LOC,XDSYM,SYMM,XNSEL,ALLD,1,ALLCP,1,UX,21,22,23,24,25CP,2,UY,2,22CP,3,UY,3,13,23CPSGEN,2,2,2,3,1FINISH/COM SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A/COM SCRATCH FILE/OUTPUT,SCRATCH*CREATE,MAC2 ! CREATE MACRO FOR ANALYSIS TYPE AND LOADING/SOLUSOLCONTROL,0ANTYPE,STATIC/COM, TEMPERATURE SET UPTREF,618TOFFST,273TUNIF,618TIME,1E-5CNVTOL,F,.00001 ! VERY SMALL MINIMUM ENFORCED! FOR CONVERGENCESOLVEOUTRES,ESOL,1 ! STORE RESULTS FOR EVERY SUBSTEPNSUBST,200TUNIF,460 ! COOLINGTIME,3160SOLVETIME,(14400+3160) ! ISOTHERMAL HOLDSOLVETUNIF,18 ! FURTHER COOLINGTIME,(14400+12000)SOLVE*ENDFINISH*USE,MAC2 ! EXECUTE ANALYSIS AND LOADING MACRO/OUTPUT/POST26ESOL,2,2,BFE,TEMPESOL,3,2,3,S,X,STRESS*CREATE,MAC3 ! MACRO FOR PROCESSING RESULTSXVAR,2/GRID,1/AXLAB,X,TEMPERATURE/AXLAB,Y,IN-PLANE STRESS (MPA)PLVAR,3*GET,MXSX,VARI,3,EXTREM,VMAX ! MAXIMUM IN-PLANE STRESSNSTORE,20 ! STORE EVERY 20TH TIME POINT RESULTSPRVAR,2,3*END*USE,MAC3 ! EXECUTE POSTPROCESSING MACRO*SET,P1,(MXSX)*GET,T1,VARI,3,EXTREM,TMAX*GET,TE,VARI,2,RTIME,T1*DIM,LABEL,CHAR,2,2*DIM,VALUE,2,3LABEL(1,1) = PRES MX ,TEMP LABEL(1,2) = MPa,DEG C*VFILL,VALUE(1,1),DATA,12.5,460*VFILL,VALUE(1,2),DATA,P1,TE*VFILL,VALUE(1,3),DATA,ABS(P1/12.5),ABS(TE/460)SAVE,TABLE_1FINISH/CLEAR,NOSTART ! CLEAR THE DATABASEJPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS/PREP7/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS/COM, - 3-D ANALYSIS -ET,1,89 ! 3-D VISCOELASTIC ELEMENT*USE,MAC1 ! EXECUTE MACRO FOR MATERIAL PROPERTIES/COM, CREATE FINITE ELEMENT MODELN,1,N,3,0.00025FILLN,5,(0.00025+0.00325)FILLNGEN,3,10,1,5,1,.001NGEN,3,100,1,25,1,0.001MAT,2E,1,3,23,21,201,203,223,221EMORE,2,13,22,11,202,213,222,211EMORE,101,103,123,121EGEN,2,2,1,1,1,-1NSLE,SNSEL,INVENDELE,ALLNSLE,S/COM, APPLY BOUNDARY CONDITIONS AND COUPLINGNSEL,S,LOC,YDSYM,SYMM,YNSEL,S,LOC,XDSYM,SYMM,XNSEL,S,LOC,ZDSYM,SYMM,ZNSEL,S,LOC,Y,0.002CP,1,UY,ALLNSEL,S,LOC,Z,0.002CP,2,UZ,ALLNSEL,S,LOC,X,0.00025CP,3,UX,ALLNSEL,S,LOC,X,0.0035CP,4,UX,ALLNSEL,ALLFINISH/COM, SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A/COM, SCRATCH FILE/OUTPUT,SCRATCH*USE,MAC2 ! EXECUTE ANALYSIS AND LOADING MACRO/OUTPUT/POST26ESOL,2,2,BFE,TEMPESOL,3,2,3,S,Y,STRESS*USE,MAC3 ! EXECUTE POSTPROCESSING MACRO*SET,P2,(MXSX)*GET,T2,VARI,3,EXTREM,TMAX*GET,TE2,VARI,2,RTIME,T2*DIM,LABEL,CHAR,2,2*DIM,VALUE,2,3LABEL(1,1) = PRES MX ,TEMP LABEL(1,2) = MPa,DEG C*VFILL,VALUE(1,1),DATA,12.5,460*VFILL,VALUE(1,2),DATA,P2,TE2*VFILL,VALUE(1,3),DATA,ABS(P2/12.5),ABS(TE2/460)SAVE,TABLE_2RESUME,TABLE_1/COM/OUT,vm200,vrt/COM,- VM200 RESULTS COMPARISON -/COM,/COM, | TARGET | ANSYS | RATIO/COM,/COM,RESULTS USING VISCO88*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,A8, ,F10.1, ,F10.1, ,1F5.3)/NOPRRESUME,TABLE_2/GOPR/COM,/COM,RESULTS USING VISCO89*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)(1X,A8,A8, ,F10.1, ,F10.1, ,1F5.3)/COM,-/OUTFINISH/NOPR/DELETE,TABLE_1/DELETE,TABLE_2/DELETE,MAC1/DELETE,MAC2/DELETE,MAC3FINISH*LIST,vm200,vrt=然后再贴几个其他的例子FINI /CLE /PREP7 !建立有限元模型ET,1,VISCO107N,1 N,2,1E-2 N,3,1