ansys 蠕变例子
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 CREEP
C*** STR. OF MATL., TIMOSHENKO, PART 2, 3RD ED., PAGE 531 ANTYPE,STATIC
ET,1,LINK1 ! SPAR ELEMENT
R,1,1,(1/30000) ! INITIAL STRAIN
MP,EX,1,30E6
TB,CREEP,1
TBDATA,1,4.8E-30,7 ! CREEP PROPERTIES
N,1
N,2,10
E,1,2
BFUNIF,TEMP,900 ! UNIFORM TEMPERATURE
TIME,1000
KBC,1
D,ALL,ALL ! FIX ALL DOFS
FINISH
/SOLU
SOLCONTROL,0
NSUBST,100
OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEP OUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEP
SOLVE
FINISH
/POST26
ESOL,2,1,,LS,1,SIG ! STORE AXIAL STRESS
PRVAR,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,3
LABEL(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,T950
V1 = 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,----------------------------------------------------------
/OUT
FINISH
*LIST,vm132,vrt
==============================
/COM,ANSYS MEDIA REL. 8.1 (03-15-2004) REF. VERIF. MANUAL: REL. 8.1 /VERIFY,VM133
/PREP7
MP,PRXY,,0.3
/TITLE, VM133, MOTION OF A ROD DUE TO IRRADIATION INDUCED CREEP
C*** REFERENCE - ANY BASIC CALCULUS BOOK
ANTYPE,STATIC
ET,1,BEAM23 ! PLASTIC BEAM
R,1,.25,.0052083,.5 ! AREA, IZZ, HEIGHT
MP,EX,1,300
TB,CREEP,1
TBDATA,55,0.5E-12,1E10 ! CREEP EQUATION CONSTANTS K1 AND K2 TBDATA,66,5 ! SELECT IRRADIATION INDUCED CREEP EQUATION
N,1
N,2,1
E,1,2
D,1,ALL ! FIX ONE END
F,2,FX,.25 ! FORCE INDUCING CONSTANT STRESS
FINISH
/SOLU
!SOLCONTROL,0
BFE,1,TEMP,1,1000,1000,1000,1000 ! APPLY CONSTANT TEMPERATURE BFE,1,FLUE,1,0,0,0,0 ! APPLY ZERO FLUENCE
TIME,1E-8 ! NEAR ZERO TIME FOR FIRST LOAD STEP
OUTPR,BASIC,1 ! PRINT BASIC ELEMENT SOLUTION
OUTRES,EPCR,1 ! STORE CREEP STRAIN RESULTS FOR EVERY SUBSTEP CNVTOL,F,,,,1E-6 ! NEAR ZERO VALUE FOR MINREF FIELD
CNVTOL,M,-1 ! CONVERGENCE CRITERION BASED UPON MOMENTS IS
! REMOVED AS IT IS NOT NEEDED FOR THIS TEST
SOLVE ! LOAD STEP 1
NSUBST,50,500,50
TIME,5
OUTPR,BASIC,5
BFE,1,FLUE,1,5E10,5E10,5E10,5E10 ! FINAL FLUENCES (RAMPED)
SOLVE ! LOAD STEP 2
FINISH
/POST26
ESOL,2,1,,LEPCR,1,EPCR ! STORE CREEP STRAIN
PRVAR,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,3
LABEL(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/.0049 7)
/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,----------------------------------------------------------
/OUT
FINISH
*LIST,vm133,vrt
==============================
/COM,ANSYS MEDIA REL. 8.1 (03-15-2004) REF. VERIF. MANUAL: REL. 8.1 /VERIFY,vm224
JPGPRF,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,
/PREP7
C*** 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 VALUE
C*** TIME PARAMETER
*SET,HOUR,1000 !ASSIGN VALUE
C*** ELASTIC CONSTANT
MP,EX,1,200E3 !DEFINE YOUNG'S MODULUS
MP,NUXY,1,0.3 !DEFINE POISON'S RATIO
TUNIF,HOT !ASSIGN TEMP TO NODES
TOFF,OFFS !SPECIFY TEMP RELATIVE TO ABSOLUTE VALUES
TB,CREEP,1,,,6 !ACTIVATE DATA TABLE
TBDATA,1,C1,C2,C3,C4 !DEFINE DATA FOR TABLE
SAVE !SAVE
/PREP7
N,1,0,0,0
N,2,100,0
N,3,100,100
N,4,0,100
N,5,50,0
N,6,100,50
N,7,50,100
N,8,0,50
N,9,50,50
ET,1,PLANE182
KEYOPT,1,1,1
KEYOPT,1,3,0
E,1,5,9,8
E,5,2,6,9
E,9,6,3,7
E,8,9,7,4
NALL
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NALL
/SOLU
NSEL,S,LOC,X,100
SF,ALL,PRES,-200
NALL
NSEL,S,LOC,Y,100
SF,ALL,PRES,200
NALL
RATE, OFF
DELT,1.0E-8,1.0E-9,1.0E-8 TIME, 1.0E-8
/OUT,SCRATCH
OUTRES,ESOL,ALL
SOLVE
/OUT
RATE, ON, ON
DELT,1E-5,1E-5,100
TIME,1000
/OUT,SCRATCH
SOLVE
/OUT
FINISH
/POST26
ESOL,2,1,,EPCR,X
ESOL,3,1,,EPCR,Y
PRVAR,2,3
PLVAR,2,3
*GET,RES1X,VARI,2,RTIME,1000
*GET,RES1Y,VARI,3,RTIME,1000 FINISH
PARSAV,ALL
RESUME
PARRES,CHANGE
/PREP7
N,1,0,0,0
N,2,100,0
N,3,100,100
N,4,0,100
N,5,50,0
N,6,100,50
N,7,50,100
N,8,0,50
ET,1,PLANE183
KEYOPT,1,3,0
E,1,2,3,4,5,6,7,8
NALL
NSEL,S,LOC,X
D,ALL,UX
NSEL,S,LOC,Y
D,ALL,UY
NALL
/SOLU
NSEL,S,LOC,X,100
SF,ALL,PRES,-200
NALL
NSEL,S,LOC,Y,100
SF,ALL,PRES,200
NALL
RATE, OFF
DELT,1.0E-8,1.0E-9,1.0E-8 TIME, 1.0E-8
/OUT,SCRATCH
SOLVE
/OUT
RATE, ON, ON
DELT,1E-5,1E-5,100
TIME,1000
/OUT,SCRATCH
OUTRES,ESOL,ALL
SOLVE
/OUT
FINISH
ESOL,2,1,,EPCR,X
ESOL,3,1,,EPCR,Y
PRVAR,2,3
PLVAR,2,3
*GET,RES2X,VARI,2,RTIME,1000
*GET,RES2Y,VARI,3,RTIME,1000
*DIM,LABEL1,CHAR,2
*DIM,VALUE1,,2,3
LABEL1(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,3
LABEL2(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,---------------------------------------------------------------- /OUT
FINISH
*LIST,vm224,vrt
==============================
/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0 /VERIFY,VM200
JPGPRF,500,100,1 ! MACRO TO SET PREFS FOR JPEG PLOTS
/SHOW,JPEG
/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 ALUMINA
TB,EVISC,2
TBDATA,31,52.6E-7 ! THERMAL EXPANSION COEFFICENTS FOR ALUMINA TBDATA,32,.119E-7 ! LOCATIONS 31-33
TBDATA,33,-1.0E-11
TBDATA,46,1.435E5 ! SHEAR MODULUS AT TIME = 0
TBDATA,47,1.435E5 ! SHEAR MODULUS AT TIME = INFINITY
TBDATA,48,3.11E5 ! BULK MODULUS AT TIME = 0
TBDATA,49,3.11E5 ! BULK MODULUS AT TIME = INFINITY
/COM, MATERIAL PROPERTIES OF G-11 GLASS
TB,EVISC,1
TBDATA,1,6.45D4 ! H/R
TBDATA,2,0.53 ! VALUE OF X
TBDATA,3,6 ! NO. OF MAXWELL ELEMENTS FOR STRUCTURAL RELAXATION TBDATA,6,0.108 ! COEFFICIENTS OF THE MAXWELL ELEMENTS FOR TBDATA,7,0.443 ! VOLUME DECAY (STRUCTURAL RELAXATION)
TBDATA,8,0.166 ! LOCATIONS 6-11
TBDATA,9,0.161
TBDATA,10,0.046
TBDATA,11,0.076
TBDATA,16,3.00 ! RELAXATION TIMES FOR VOLUME DECAY FUNCTION TBDATA,17,0.671 ! LOCATIONS 16-21
TBDATA,18,0.247
TBDATA,19,0.091
TBDATA,20,0.033
TBDATA,21,0.008
TBDATA,26,3.43E-5 ! THERMAL EXPANSION COEFF. IN LIQUID STATE TBDATA,31,64.7E-7 ! THERMAL EXPANSION COEFF. IN GLASS (SOLID) STATE TBDATA,32,.02E-7 ! LOCATIONS 31-32
TBDATA,36,618 ! FICTIVE TEMPERATURE LOCATIONS 36-41
TBDATA,37,618
TBDATA,38,618
TBDATA,39,618
TBDATA,40,618
TBDATA,41,618
TBDATA,46,2.79D4 ! SHEAR MODULUS AT TIME = 0
TBDATA,47,0.0D0 ! SHEAR MODULUS AT TIME = INFINITY
TBDATA,48,6.05D4 ! BULK MODULUS AT TIME = 0
TBDATA,49,6.05D4 ! BULK MODULUS AT TIME = INFINITY
TBDATA,50,3 ! THREE MAXWELL ELEMENTS FOR SHEAR RELAXATION TBDATA,51,0.422 ! COEFF. FOR SHEAR RELAXATION LOCATIONS 51-53 TBDATA,52,0.423
TBDATA,53,0.155
TBDATA,61,0.0689 ! RELAXATION TIMES FOR SHEAR RELAXATION TBDATA,62,0.0065 ! LOCATIONS 61-63
TBDATA,63,0.0001
TBDATA,71,0 ! NO BULK MODULUS RELAXATION
*END
*USE,MAC1 ! EXECUTE MACRO FOR MATERIAL PROPERTIES
:COM, CREATE FINITE ELEMENT MODEL
N,1,
N,3,,0.00025
FILL
N,5,0,(0.00025+0.00325)
FILL
NGEN,3,10,1,5,1,.001
MAT,2
E,1,21,23,3,11,22,13,2
MAT,1
E,3,23,25,5,13,24,15,4
/COM, APPLY BOUNDARY CONDITIONS AND COUPLING
NSEL,S,LOC,Y
DSYM,SYMM,Y
NSEL,S,LOC,X
DSYM,SYMM,X
NSEL,ALL
D,1,ALL
CP,1,UX,21,22,23,24,25
CP,2,UY,2,22
CP,3,UY,3,13,23
CPSGEN,2,2,2,3,1
FINISH
/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
/SOLU
SOLCONTROL,0
ANTYPE,STATIC
/COM, TEMPERATURE SET UP
TREF,618
TOFFST,273
TUNIF,618
TIME,1E-5
CNVTOL,F,,,,.00001 ! VERY SMALL MINIMUM ENFORCED
! FOR CONVERGENCE
SOLVE
OUTRES,ESOL,1 ! STORE RESULTS FOR EVERY SUBSTEP NSUBST,200
TUNIF,460 ! COOLING
TIME,3160
SOLVE
TIME,(14400+3160) ! ISOTHERMAL HOLD
SOLVE
TUNIF,18 ! FURTHER COOLING
TIME,(14400+12000)
SOLVE
*END
FINISH
*USE,MAC2 ! EXECUTE ANALYSIS AND LOADING MACRO
/OUTPUT
/POST26
ESOL,2,2,,BFE,TEMP
ESOL,3,2,3,S,X,STRESS
*CREATE,MAC3 ! MACRO FOR PROCESSING RESULTS
XVAR,2
/GRID,1
/AXLAB,X,TEMPERATURE
/AXLAB,Y,IN-PLANE STRESS (MPA)
PLVAR,3
*GET,MXSX,VARI,3,EXTREM,VMAX ! MAXIMUM IN-PLANE STRESS NSTORE,20 ! STORE EVERY 20TH TIME POINT RESULTS PRVAR,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,3
LABEL(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_1
FINISH
/CLEAR,NOSTART ! CLEAR THE DATABASE
JPGPRF,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 MODEL
N,1,
N,3,0.00025
FILL
N,5,(0.00025+0.00325)
FILL
NGEN,3,10,1,5,1,,.001
NGEN,3,100,1,25,1,,,0.001
MAT,2
E,1,3,23,21,201,203,223,221
EMORE,2,13,22,11,202,213,222,211
EMORE,101,103,123,121
EGEN,2,2,1,1,1,-1
NSLE,S
NSEL,INVE
NDELE,ALL
NSLE,S
/COM, APPLY BOUNDARY CONDITIONS AND COUPLING
NSEL,S,LOC,Y
DSYM,SYMM,Y
NSEL,S,LOC,X
DSYM,SYMM,X
NSEL,S,LOC,Z
DSYM,SYMM,Z
NSEL,S,LOC,Y,0.002
CP,1,UY,ALL
NSEL,S,LOC,Z,0.002
CP,2,UZ,ALL
NSEL,S,LOC,X,0.00025
CP,3,UX,ALL
NSEL,S,LOC,X,0.0035
CP,4,UX,ALL
NSEL,ALL
FINISH
/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
/POST26
ESOL,2,2,,BFE,TEMP
ESOL,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,3
LABEL(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_2
RESUME,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)
/NOPR
RESUME,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,----------------------------------------------------------- /OUT
FINISH
/NOPR
/DELETE,TABLE_1
/DELETE,TABLE_2
/DELETE,MAC1
/DELETE,MAC2
/DELETE,MAC3
FINISH
*LIST,vm200,vrt
============================== 然后再贴几个其他的例子
FINI
/CLE
/PREP7
!建立有限元模型
ET,1,VISCO107
N,1
N,2,1E-2
N,3,1E-2,1E-2
N,4,,1E-2
NGEN,2,4,1,4,1,,,1E-2
E,1,2,3,4,5,6,7,8
D,1,ALL,,,2
D,5,ALL,,,6
Save,visoplas,db
!定义粘弹性材料参数
Resume,visoplas,db
MP,EX,1,60.6E9
MP,NUXY,1,0.4999
TB,ANAND,1 !!!!关键,定义材料内型TBDATA,1,29.7E6
TBDATA,2,21.08999E3
TBDATA,3,1.91E7
TBDATA,4,7.0
TBDATA,5,0.23348
TBDATA,6,1115.6E6
TBDATA,7,18.92E6
TBDATA,8,0.07049
TBDATA,9,1.3
FINISH
!求解(静力分析)
/SOLU
ANTYPE,STATIC
NLGEOM,ON ! 打开几何大变形OUTRES,RSOL,ALL
BFUNIF,TEMP,673 !施加初始温度D,3,ALL,0.0,,8
TIME,0.000001
SOLVE
NSUBST,20 !定义该载荷步的子步数D,3,UX,0.2E-2,,4
D,7,UX,0.2E-2,,8
TIME,20
SOLVE
FINISH
!后处理
/POST26
RFORCE,2,3,F,X
RFORCE,3,4,F,X
RFORCE,4,7,F,X
RFORCE,5,8,F,X
ADD,6,2,3,4
ADD,7,6,5,,LOAD
PRVAR,7
*GET,F1,VARI,7,RTIME,20 FINISH
============================== fini
/cle
/fil,exam0610
/COM, Structural
/config,NRES,6000
/prep7
et,1,beam23
r,1,10*9,10*9**3/12,9,,,, mp,ex,1,3.5e10
mp,prxy,1,0.3
mp,dens,1,7800
acel,9.8,
k,1
k,2,10
l,1,2
lsel,all
latt,1,1,1
esize,5
lmesh,all
nsel,s,loc,x,0,40
esln,s,1,all
cm,kuai0,elem
allsel,all
finish
/solu
eqslv,sparse
time,7
NLGEON,ON
NROPT,FULL
ESTIF,1e-5
ALLSEL
EKILL,ALL
allsel
nsel,s,loc,x,0
d,all,all
allsel
bfunif,temp,300
ncnv,2
kbc,1
autots,on
tb,creep,1
tbdata,1,0,0,1
nsubst,40
outpr,basic,1
outres,epcr,1
solve
save
rate,1
esel,s,,,kuai0
ealive,all
esel,s,live
nsle,s
NSEL,INVE
D,ALL,ALL,0
allsel
nsel,s,loc,x,0
d,all,all
time,100
tb,creep,1
tbdata,1,1e-5,0,1,0,0,
allsel
solve
finish
save
==============================
有限元ansys静力分析的一个小例子
有限元 学院：机电学院 专业： 姓名： 学号：
一、问题描述 如图所示的平面，板厚为0.01m，左端固定，右端作用50kg的均布载荷，对其进行静力分析。弹性模量为210GPa，泊松比为0.25. 二、分析步骤 1．启动ansys，进入ansys界面。 2．定义工作文件名 进入ANSYS/Multiphsics的的程序界面后，单击Utility Menu菜单下File中Change Jobname的按钮，会弹出Change Jobname对话框，输入gangban为工作文件名，点击ok。 3．定义分析标题 选择菜单File-Change Title在弹出的对话框中，输入Plane Model作为分析标题，单击ok。 4．重新显示 选择菜单Plot-Replot单击该按钮后，所命令的分析标题工作文件名出现在ANSYS 中。 5．选择分析类型 在弹出的对话框中，选择分析类型，由于此例属于结构分析，选择菜单Main Menu：Preferences，故选择Structural这一项，单击ok。 6．定义单元类型 选择菜单Main Menu-Preprocessor-Element Type-Add/Edit/Delete单击弹出对话框中的Add按钮，弹出单元库对话框，在材料的单元库中选Plane82单元。即在左侧的窗口中选取Solid单元，在右侧选择8节点的82单元。然后单击ok。 7．选择分析类型 定义完单元类型后，Element Type对话框中的Option按钮被激活，单击后弹出一个对话框，在Elenment behavior中选择Plane strs w/ thk，在Extra Element output 中，选择Nodal stress，单击close，关闭单元类型对话框。 8．定义实常数 选择菜单Main Menu-Preprocessor-Real Constants Add/Edit/Delete执行该命令后，在弹出Real Constants对话框中单击Add按钮，确认单元无误后，单击ok，弹出Real Constants Set Number 1，for Plane 82对话框，在thickness后面输入板的厚度0.01单击ok，单击close。 9．定义力学参数 选择菜单Main Menu-Preprocessor-Material Props-Material Model 在弹出的对
七个ansys经典入门实例
“有限元分析及应用”课程有限元分析软件ANSYS6.xed 上机指南 清华大学机械工程系 2002年9月
说明 本《有限元分析软件ANSYS6.1ed:上机指南》由清华大学机械工程系石伟老师组织编写，由助教博士生孔劲执笔, 于2002年9月完成，基本操作指南中的所有算例都在相应的软件系统中进行了实际调试和通过。 本上机指南的版权归清华大学机械工程系所有，未经同意，任何单位和个人不得翻印。
目录 Project1 简支梁的变形分析 (1) Project2 坝体的有限元建模与受力分析 (3) Project3 受内压作用的球体的应力与变形分析 (5) Project4 受热载荷作用的厚壁圆筒的有限元建模与温度场求解 (7) Project5 超静定桁架的有限元求解 (9) Project6 超静定梁的有限元求解 (11) Project7 平板的有限元建模与变形分析 (13)
Project1 梁的有限元建模与变形分析 计算分析模型如图1-1 所示, 习题文件名: beam。 NOTE:要求选择不同形状的截面分别进行计算。 梁承受均布载荷：1.0e5 Pa 图1-1梁的计算分析模型 梁截面分别采用以下三种截面（单位：m）： 矩形截面：圆截面：工字形截面： B=0.1, H=0.15 R=0.1 w1=0.1，w2=0.1，w3=0.2, t1=0.0114，t2=0.0114，t3=0.007 1.1进入ANSYS 程序→ANSYSED 6.1 →Interactive →change the working directory into yours →input Initial jobname: beam→Run 1.2设置计算类型 ANSYS Main Menu: Preferences →select Structural →OK 1.3选择单元类型 ANSYS Main Menu: Preprocessor →Element Type→Add/Edit/Delete… →Add… →select Beam 2 node 188 →OK (back to Element Types window)→Close (the Element Type window) 1.4定义材料参数 ANSYS Main Menu: Preprocessor →Material Props →Material Models →Structural→Linear→Elastic→Isotropic→input EX:2.1e11, PRXY:0.3→OK 1.5定义截面 ANSYS Main Menu: Preprocessor →Sections →Beam →Common Sectns→分别定义矩形截面、圆截面和工字形截面：矩形截面:ID=1,B=0.1，H=0.15 →Apply →圆截面：ID=2,R=0.1 →Apply →工字形截面：ID=3,w1=0.1，w2=0.1，w3=0.2，t1=0.0114，t2=0.0114，t3=0.007→OK
ansys工程实例(4经典例子)解析
输气管道受力分析（ANSYS建模） 任务和要求： 按照输气管道的尺寸及载荷情况,要求在ANSYS中建模，完成整个静力学分析过程。求出管壁的静力场分布。要求完成问题分析、求解步骤、程序代码、结果描述和总结五部分。所给的参数如下： 材料参数：弹性模量E=200Gpa; 泊松比0.26；外径R?=0.6m；内径R?=0.4m；壁厚t=0.2m。输气管体内表面的最大冲击载荷P为1Mpa。 四.问题求解 （一）．问题分析 由于管道沿长度方向的尺寸远大于管道的直径，在计算过程中忽略管道的端面效应，认为在其长度方向无应变产生，即可将该问题简化为平面应变问题，选取管道横截面建立几何模型进行求解。 （二）．求解步骤 定义工作文件名 选择Utility Menu→File→Chang Jobname 出现Change Jobname对话框，在[/FILNAM] Enter new jobname 输入栏中输入工作名LEILIN10074723，并将New log and eror file 设置为YES，单击[OK]按钮关闭对话框 定义单元类型 1)选择Main Meun→Preprocessor→Element Type→Add/Edit/Delte命令，出现Element Type 对话框，单击[Add]按钮，出现Library of Element types对话框。 2）在Library of Element types复选框选择Strctural、Solid、 Quad 8node 82,在Element type reference number输入栏中出入1，单击[OK]按钮关闭该对话框。 3. 定义材料性能参数 1）单击Main Meun→Preprocessor→Material Props→Material models出现Define Material Behavion 对话框。选择依次选择Structural、Linear、Elastic、Isotropic选项，出现Linear Isotropic Material Properties For Material Number 1对话框。 2）在EX输入2e11，在Prxy输入栏中输入0.26，单击OK按钮关闭该对话框。 3）在Define Material Model Behavion 对话框中选择Material→Exit命令关闭该对话框。 4.生成几何模型、划分网格 1）选择Main Meun→Preprocessor→Modeling→Create→Areas→Circle→Partail→Annulus出现Part Annulus Circ Area对话框，在WP X文本框中输入0，在WP Y文本框中输入0，在Rad1文本框中输入0.4，在Theate-1文本框中输入0，在Rad2文本框中输入0.6，在Theate-2文本框中输入90，单击OK按钮关闭该对话框。 2）选择Utility Menu→Plotctrls→Style→Colors→Reverse Video，设置显示颜色。 3）选择Utility Menu→Plot→Areas，显示所有面。 4) 选择Main Menu→Preprocessor→Modeling→Reflect→Areas，出现Reflect Areas拾取菜
ANSYS瞬态分析实例
例题：一根钢梁支撑着集中质量并承受一个动态载荷（如图1所示）。钢梁长为L，支撑着一个集中质量M。这根梁承受着一个上升时间为t1的值为F1 的动态载荷F（t）。梁的质量可以忽略，确定产生最大位移响应时的时间t max 和响应y max。 图1 钢梁支撑集中质量的几何模型 材料特性：弹性模量为2e5MPa，质量为M＝0.0215t，质量阻尼为8； 几何尺寸为：L＝450mm，I＝800.6mm4，h＝18mm； 载荷为：F1＝20N，t1＝0.075s GUI操作方式： 1．定义单元类型：Main Menu>Preprocessor>Element Type>Add/Edit/Delete，出现一个对话框，单击“Add”，又出现一个对话框，在对话框左面的列表栏中选择“Structural Beam”，在右面的列表栏中选择“2D elastic 3”，单击“Apply”，在对话框左面的列表栏中选择“Structural Mass”，在右边选择“3D mass 21”，单击“OK”，在单击“Options”，弹出对话框，设置K3为“2-D W/O rot iner”，单击“OK”，再单击“Close”。 2．设置实常数：Main Menu>Preprocessor>Real Constants> Add/Edit/Delete，出现对话框，单击“Add”，又弹出对话框，选择“Type1 BEAM3”，单击“OK”，
又弹出对话框，输入AREA为1，IZZ＝800.6，HEIGHT＝18，单击“OK”，在单击“Add”，选择Type 2 MASS21，单击“OK”，设置MASS为0.0215，单击“OK”，再单击“Close”。 3．定义材料属性：Main Menu>Preprocessor>Material Props>Material Modls，出现对话框，在“Material Models Available”下面的对话框中，双击打开“Structural>Linear>Elastic>Isotropic”，又出现一个对话框，输入弹性模量EX＝2e5，泊松比PRXY＝0，单击“OK”，单击“Materal>Exit”。 4．建立模型： 1）创建节点：依次单击Main Menu>Preprocessor>Modeling>Create>Nodes>In Active CS，在弹出对话框中，依次输入节点的编号1，节点坐标x＝0，y ＝0，然后单击“Apply”，输入节点编号2，节点坐标x＝450/2，y＝0，然后单击“Apply”，输入节点编号3，节点坐标x＝450，y＝0。单击“OK”。2）创建单元：依次单击Main Menu>Preprocessor>Modeling>Create>Elements >Auto Numbered>Thru Nodes，弹出拾取框，拾取节点1和2，2和3，单击“OK”。 3）指定单元实常数：Main Menu>Preprocessor>Modeling>Create>Elements> Elem Attributes，弹出对话框，设置TYPE为2，REAL为2，单击“OK”。4）创建单元：依次单击Main Menu>Preprocessor>Modeling>Create>Elements >Auto Numbered>Thru Nodes，弹出拾取框，拾取节点2，单击“OK”。5．定义分析类型：Main Menu>Solution>Analysis Type>New Analysis，弹出对话框，选择Trasiernt，单击“OK”，又弹出对话框，选择Reduced，单击“OK”。6．设置分析选项：Main Menu>Solution>Analysis Type>Analysis Options，弹出对话框，单击“OK”。
ANSYS 中使用函数加载的一个简单例子
ANSYS 中使用函数加载的一个简单例子 本文将通过一个具体实例说明在ANSYS 中如何使用函数加载，后续将通过该实例在分析过程中遇到的一个问题提出自己的一点看法。 实例的具体说明： 一个1/4 圆柱，内半径30 mm，外半径42 mm，长度100mm，如图1 所示： 所用材料为双线性弹塑性材料，其机械性能为： 弹性模量 E = 201000 Mpa；泊松比μ=0.3
屈服应力σ= 200 Mpa；切线模量Et = 2010使用单元类型solid185 (8 节点六面体单元)。 取整体单元边长4 mm，然后可以直接对该几何模型划分MAP 网格，划分网格结果如图2： 约束条件为： 轴向两个截面为对称边界条件；一个端面约束轴向位移Uz。 载荷条件为： 在外表面施加变化的压力载荷，载荷函数为： P (y) = 8e7 + 7E7 * (Y/42)
即： X = 0 ，Y = 42 (最高点) 时，P = 15E7； X = 42，Y = 0 （最低点）时，P = 8E7。 我们采用函数方式来施加这一压力载荷，首先定义函数： 在Solution 模块中，点击菜单路径： Solution > Define Loads > Apply > Functions > Define/Edit 将会弹出一个函数编辑器，可以在其中定义所需的函数。 在函数编辑器中，函数类型选择为Single equation，即单值函数；计算函数值时使用的插值坐标系( (x,y,z) interpreted in CSYS) 选择0，即总体直角坐标系，如图3 所示：
然后，在函数编辑器中间位置的“Result = “ 小窗口中输入要定义的函数表达式，如果表达式中有x, y, z, time 等变量(供定义函数时使用的“自变量”)，可以用{X},{Y},{Z},{TIME} 等的形式输入；或者点击下面一个小窗口右边的小箭头，会出现一个下拉列表，列出可以选择的变量，然后从该列表中选择某个自变量，则该自变量会按照上述格式写入函数中，如图5 所示：
ansys结构瞬态分析实例
第二日 练习主题：各种网格划分方法 输入实体模型尝试用映射、自由网格划分，并综合利用多种网格划分控制方法 一个瞬态分析的例子 练习目的：熟悉瞬态分析过程 练习过程：瞬态(FULL)完全法分析板-梁结构实例 如图1所示板-梁结构，板件上表面施加随时间变化的均布压力，计算在下列已知条件下结构的瞬态响应情况。 全部采用A3钢材料，特性： 杨氏模量=2e112 /m N 泊松比=0.3 密度=7.8e33/m Kg 板壳： 厚度=0.02m 四条腿（梁）的几何特性： 截面面积=2e-42m 惯性矩=2e-84m 宽度=0.01m 高度=0.02m 压力载荷与时间的关系曲线如图2所示。 图1 质量梁-板结构及载荷示意图 压力（N/m 2） 10000 5000 0 1 2 4 6 时间（s ） 图 2 板上压力-时间关系 分析过程 第1步：设置分析标题 1. 选取菜单途径Utility Menu>File>Change Title 。 2. 输入“ The Transient Analysis of the structure ”，然后单击OK 。 第2步：定义单元类型 单元类型1为SHELL63，单元类型2为BEAM4 第3步：定义单元实常数 实常数1为壳单元的实常数1，输入厚度为0.02（只需输入第一个值，即等78厚度壳）
实常数2为梁单元的实常数，输入AREA 为2e-4惯性矩IZZ=2e-8，IYY ＝2e-8宽度TKZ=0.01，高度TKY=0.02。 第5步：杨氏模量EX=2e112 /m N 泊松比NUXY=0.3 密度DENS=7.8e33/m Kg 第6步：建立有限元分析模型 1． 创建矩形，x1=0,x2=2,y1=0,y2=1 2． 将所有关键点沿Z 方向拷贝，输入DZ ＝－1 3． 连线。将关键点1，5；2，6；3，7；4，8分别连成直线。 4． 设置线的分割尺寸为0.1，首先给面划分网格；然后设置单元类型为2，实常数为2， 对线5到8划分网格。 第7步：瞬态动力分析 1. 选取菜单途径Main Menu>Solution>-Analysis Type-New Analysis ，弹出New Analysis 对话框。 2. 选择Transient ，然后单击OK ，在接下来的界面仍然单击OK 。 3. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Damping ，弹出 Damping Specifications 窗口。 4. 在Mass matrix multiplier 处输入5。单击OK 。 5. 选取菜单途径Main Menu > Solution > -Loads-Apply > -Structural- Displacement>On Nodes 。弹出拾取（Pick ）窗口，在有限元模型上点取节点232、242、252和262，单击OK ，弹出Apply U,ROT on Nodes 对话框。 6. 在DOFS to be constrained 滚动框中，选种“All DOF ”(单击一次使其高亮度显示， 确保其它选项未被高亮度显示）。单击OK 。 7. 选取菜单途径Utility Menu>Select>Everything 。 8. 选取菜单途径Main Menu>Solution>-Load Step Opts-Output Ctrls>DB/Results File ，弹 出Controls for Database and Results File Writing 窗口。 9. 在Item to be controlled 滚动窗中选择All items ，下面的File write frequency 中选择 Every substep 。单击OK 。 10. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Time – Time Step ，弹出Time – Time Step Options 窗口。 11. 在Time at end of load step 处输入1；在Time step size 处输入0.2；在Stepped or ramped b.c 处单击ramped ；单击Automatic time stepping 为on ；在Minimum time step size 处输入0.05；在Maximum time step size 处输入0.5。单击OK 。 12. 选取菜单途径Main Menu>Solution>-Loads-Apply>-Structure-Pressure>On Areas 。弹 出Apply PRES on Areas 拾取窗口。 13. 单击Pick All ，弹出Apply PRES on Areas 对话框。 14. 在pressure value 处输入10000。单击OK 15. 选取菜单途径Main menu>Solution>Write LS File ，弹出Write Load Step File 对话框。 16. 在Load step file number n 处输入1，单击OK 。 17. 选取菜单途径Main Menu>Solution>-Load Step Opts-Time/Frequenc> Time – Time Step ，弹出Time – Time Step Options 窗口。 18. 在Time at end of load step 处输入2。单击单击OK 。
第23例 材料蠕变分析实例
第23例材料蠕变分析实例—受拉平板本例简单地介绍了蠕变的概念及蠕变材料模型的创建方法，简单地介绍了结构蠕变分析的方法、步骤及要点。 23.1蠕变简介 蠕变是指金属材料在长时间的恒温、恒载作用下，持续发生缓慢塑性变形的行为，大多数金属材料在高温下都会表现出蠕变行为。 如果材料发生了蠕变，在恒载作用下结构会发生持续变形；如果结构承受恒位移，则应力会随时间而减小，即产生应力松弛。 图23-1 蠕变曲线 蠕变一般分为蠕变初始阶段（Primary）、蠕变稳定阶段（Secondary）和蠕变加速阶段(Tertiary)三个阶段，如图23-1所示。蠕变初始阶段时间很短，应变率随时间而减小；在蠕变稳定阶段，应变以常速率发展；在蠕变加速阶段，应变率急剧增大直至材料失效。研究蠕变行为，主要针对蠕变初始阶段和蠕变稳定阶段。 研究问题时一般以蠕变方程（又称本构关系）来表征蠕变行为，蠕变方程以蠕应变率的，形式表示dεcr/dt =AσBεC t P式中，εcr为蠕应变。A、B、C、D是由实验得到的材料特性参数。当D<0时，蠕应变率随时间减小，材料处于蠕变初始阶段；当D=0时，蠕应变率不随时间变化，材料处于蠕变稳定阶段。
在ANSYS中，有一个蠕应变率库供选择。 23.2问题描述 一矩形平板，左端固定，右端作用有恒定压力p=100MPa，矩形平板尺寸如图23-2所示，材料的弹性模量为2xl05MPa，泊松比为0.3，蠕变稳定阶段蠕变方程dεcr/dt =C1σC2。C2，式中，C1=3.125 x10-14，C2=5。试分析平板右端的位移随时间的变化情况。 提示：为避免出现较小值，力单位用N，长度单位用mm，时间单位为h。 图23-2受拉矩形平板 23.3分析步骤 23.3.1改变任务名 拾取菜单Utility Menu→File→Change Jobname，弹出如图23-3所示的对话框，在“[/FJLNAM]”文本框中输入EXAMPLE23，单击“OK”按钮。 图23-3改变任务名对话框 23.3.2选择单元类型 拾取菜单Main Menu→Preprocessor→Element Type→Add/Edit／Delete，弹出如图23-4所示的对话框，单击“Add…”按钮，弹出如图23-5所示的对话框，
ANSYS动力学瞬态分析完全法
完全法 完全法采用完整的系统矩阵计算瞬态响应（没有矩阵缩减）。它是三种方法中功能最强的，允许包括各类非线性特性（塑性、大变形、大应变等）。 注─如果并不想包括任何非线性，应当考虑使用另外两种方法中的一种。这是因为完全法是三种方法中开销最大的一种。 完全法的优点是： ·容易使用，不必关心选择主自由度或振型。 ·允许各种类型的非线性特性。 ·采用完整矩阵，不涉及质量矩阵近似。 ·在一次分析就能得到所有的位移和应力。 ·允许施加所有类型的载荷：节点力、外加的（非零）位移（不建议采用）和单元载荷（压力和温度），还允许通过TABLE数组参数指定表边界条件。 ·允许在实体模型上施加的载荷。 完全法的主要缺点是它比其它方法开销大。 §3.4 完全法瞬态动力学分析 首先，讲述完全法瞬态动力学分析过程，然后分别介绍模态叠加法和缩减法与完全法不相同的计算步骤。完全法瞬态动力分析（在ANSYS/Multiphsics、ANSYS/Mechauioal及ANSYS/Structural中可用）由以下步骤组成： 1.建造模型 2.建立初始条件 3.设置求解控制 4.设置其他求解选项 5.施加载荷 6.存储当前载荷步的载荷设置 7.重复步骤3-6定义其他每个载荷步
8.备份数据库 9.开始瞬态分析 10.退出求解器 11.观察结果 § 型 在这一步中，首先要指定文件名和分析标题，然后用PREP7定义单元类型，单元实常数，材料性质及几何模型。这些工作在大多数分析中是相似的。<
数值模拟分析实例
华中科技大学体育馆数值模拟分析 6.1分析模型的建立 采用有限元软件ANSYS建立该网壳结构有限元分析模型。整体屋盖结构共计1481个节点，4430个单元，16种截面类型。建模时，网壳结构主体结构部分（包括主桁架、次桁架、水平支撑和檩条）采用ANSYS的LINK8杆单元建模，两侧翼的主梁、次梁和支承钢管柱均采用BEAM4梁单元，网壳结构屋面下部混凝土支承结构亦采用BEAM4梁单元。分析时，屋面板、设备管线等荷载等效为节点荷载，施加在结构节点上。 在网壳结构有限元分析中，对于杆件采用的LINK8 3-D Spar单元为三维单元，假设材料为均质等直杆，且在轴向上施加载荷，可以承受单向的拉伸或者压缩，每个节点上具有三个自由度，即沿X、Y和Z坐标轴方向。该单元具有塑性、蠕变、应力硬化和大变形等功能，能较好的模拟三维空间桁架单元。 对于两侧翼结构和下部支撑体系的柱、梁等结构采用的BEAM4单元是一个轴向拉压、扭转和弯曲单元，每个节点有三个平动自由度和三个转动自由度，具有应力刚化和大变形功能。 施工过程模拟分析时考虑时，同时考虑温度效应影响，计算时材料假定为理想弹塑性材料。
图6-1 有限元分析模型 6.2分析工况选取 按照实际施工顺序，将网壳结构屋盖施工过程划分为5个工况进行施工数值模拟，计算温度取为该阶段施工完成时的环境温度。 工况1: 7榀拱形主桁架安装完毕，但临时支撑未撤除，计算温度为温度15℃； （a）短轴立面
（b）长轴立面 图6-2 工况1中屋盖结构平面图图6-3 工况1中屋盖结构立面图工况2: 两侧翼结构安装完毕，完成后拆除其临时支撑，计算温度为8℃； （a）短轴立面 （b）长轴立面 图6-4 工况2中屋盖结构平面图图6-5 工况2中屋盖结构立面图工况3: 次桁架、水平支撑及楼梯安装完毕，临时支撑拆除，计算温度为29℃； （a）短轴立面 （b）长轴立面 图6-6 工况2中屋盖结构平面图图6-7 工况2中屋盖结构立面图工况4: 檩条及设备管线安装完毕，计算温度为41℃；
几个ansys经典实例(长见识)
平面问题斜支座的处理 如图5-7所示，为一个带斜支座的平面应力结构，其中位置2及3处为固定约束，位置4处为一个45o的斜支座，试用一个4节点矩形单元分析该结构的位移场。 (a)平面结构(b)有限元分析模型 图5-7 带斜支座的平面结构 基于ANSYS平台，分别采用约束方程以及局部坐标系的斜支座约束这两种方式来进行处理。 (7) 模型加约束 左边施加X,Y方向的位移约束 ANSYS Main Menu: Solution →Define Loads →Apply →-Structural→Displacement On Nodes →选取2，3号节点→OK →Lab2: All DOF(施加X,Y方向的位移约束) →OK 以下提供两种方法处理斜支座问题，使用时选择一种方法。 ?采用约束方程来处理斜支座 ANSYS Main Menu：Preprocessor →Coupling/ Ceqn →Constraint Eqn ：Const :0, NODE1:4, Lab1: UX,C1:1,NODE2:4,Lab2:UY,C2:1→OK 或者?采用斜支座的局部坐标来施加位移约束 ANSYS Utility Menu：WorkPlane →Local Coordinate System →Create local system →At specified LOC + →单击图形中的任意一点→OK →XC、YC、ZC分别设定为2，0，0,THXY：45 →OK ANSYS Main Menu：Preprocessor →modeling →Move / Modify →Rotate Node CS →To active CS → 选择4号节点 ANSYS Main Menu：Solution →Define Loads →Apply →Structural →Displacement On Nodes →选取4号节点→OK →选择Lab2:UY(施加Y方向的位移约束) →OK 命令流； !---方法1 begin----以下的一条命令为采用约束方程的方式对斜支座进行处理 CE,1,0,4,UX,1,4,UY,-1 !建立约束方程(No.1): 0=node4_UX*1+node_UY*(-1) !---方法1 end --- !--- 方法2 begin --以下三条命令为定义局部坐标系，进行旋转，施加位移约束 !local,11,0,2,0,0,45 !在4号节点建立局部坐标系 !nrotat, 4 !将4号节点坐标系旋转为与局部坐标系相同 !D,4,UY !在局部坐标下添加位移约束 !--- 方法2 end
ansys动力学瞬态分析详解
§3.1瞬态动力学分析的定义 瞬态动力学分析（亦称时间历程分析）是用于确定承受任意的随时间变化载荷结构的动力学响应的一种方法。可以用瞬态动力学分析确定结构在稳态载荷、瞬态载荷和简谐载荷的随意组合作用下的随时间变化的位移、应变、应力及力。载荷和时间的相关性使得惯性力和阻尼作用比较重要。如果惯性力和阻尼作用不重要，就可以用静力学分析代替瞬态分析。 瞬态动力学的基本运动方程是： 其中： [M] =质量矩阵 [C] =阻尼矩阵 [K] =刚度矩阵 {}=节点加速度向量 {}=节点速度向量 {u} =节点位移向量 在任意给定的时间，这些方程可看作是一系列考虑了惯性力（[M]{}）和 阻尼力（[C]{}）的静力学平衡方程。ANSYS程序使用Newmark时间积分方法在离散的时间点上求解这些方程。两个连续时间点间的时间增量称为积分时间步长（integration time step）。 §3.2学习瞬态动力学的预备工作 瞬态动力学分析比静力学分析更复杂，因为按“工程”时间计算，瞬态动力学分析通常要占用更多的计算机资源和更多的人力。可以先做一些预备工作以理解问题的物理意义，从而节省大量资源。例如，可以做以下预备工作：
1.首先分析一个较简单模型。创建梁、质量体和弹簧组成的模型，以最小的代价深入的理解动力学认识，简单模型更有利于全面了解所有的动力学响应所需要的。 2.如果分析包括非线性特性，建议首先利用静力学分析掌握非线性特性对结构响应的影响规律。在某些场合，动力学分析中是没必要包括非线性特性的。 3.掌握结构动力学特性。通过做模态分析计算结构的固有频率和振型，了解这些模态被激活时结构的响应状态。同时，固有频率对计算正确的积分时间步长十分有用。 4.对于非线性问题，考虑将模型的线性部分子结构化以降低分析代价。<<高级技术分指南>>中将讲述子结构。 §3.3三种求解方法 瞬态动力学分析可采用三种方法：完全（Full）法、缩减（Reduced）法及模态叠加法。ANSYS/Professional产品中只允许用模态叠加法。在研究如何实现这些方法之前，让我们先探讨一下各种方法的优点和缺点。 §3.3.1完全法 完全法采用完整的系统矩阵计算瞬态响应（没有矩阵缩减）。它是三种方法中功能最强的，允许包括各类非线性特性（塑性、大变形、大应变等）。 注─如果并不想包括任何非线性，应当考虑使用另外两种方法中的一种。这是因为完全法是三种方法中开销最大的一种。 完全法的优点是： ·容易使用，不必关心选择主自由度或振型。 ·允许各种类型的非线性特性。 ·采用完整矩阵，不涉及质量矩阵近似。 ·在一次分析就能得到所有的位移和应力。 ·允许施加所有类型的载荷：节点力、外加的（非零）位移（不建议采用）和单元载荷（压力和温度），还允许通过TABLE数组参数指定表边界条件。 ·允许在实体模型上施加的载荷。 完全法的主要缺点是它比其它方法开销大。
一个经典的ansys热分析实例(流程序)
/PREP7 /TITLE,Steady-state thermal analysis of pipe junction /UNITS,BIN ! 英制单位；Use U. S. Customary system of units (inches) ! /SHOW, ! Specify graphics driver for interactive run ET,1,90 ! Define 20-node, 3-D thermal solid element MP,DENS,1,.285 ! Density = .285 lbf/in^3 MPTEMP,,70,200,300,400,500 ! Create temperature table MPDATA,KXX,1,,8.35/12,8.90/12,9.35/12,9.80/12,10.23/12 ! 指定与温度相对应的数据材料属性；导热系数；Define conductivity values MPDATA,C,1,,.113,.117,.119,.122,.125 ! Define specific heat values（比热） MPDATA,HF,2,,426/144,405/144,352/144,275/144,221/144 ! Define film coefficient;除144是单位问题，上面的除12也是单元问题 ! Define parameters for model generation RI1=1.3 ! Inside radius of cylindrical tank RO1=1.5 ! Outside radius Z1=2 ! Length RI2=.4 ! Inside radius of pipe RO2=.5 ! Outside pipe radius Z2=2 ! Pipe length CYLIND,RI1,RO1,,Z1,,90 ! 90 degree cylindrical volume for tank WPROTA,0,-90 ! 旋转当前工作的平面；从Y到Z旋转-90度；；Rotate working plane to pipe axis CYLIND,RI2,RO2,,Z2,-90 ! 角度选择在了第四象限；90 degree cylindrical volume for pipe WPSTYL,DEFA ! 重新安排工作平面的设置；另外WPSTYL,STAT to list the status of the working plane；；Return working plane to default setting BOPT,NUMB,OFF ! 关掉布尔操作的数字警告信息；Turn off Boolean numbering warning VOVLAP,1,2 ! 交迭体；Overlap the two cylinders /PNUM,VOLU,1 ! 体编号打开；Turn volume numbers on /VIEW,,-3,-1,1
ANSYSMaxwell瞬态分析案例解析
1.Maxwell 2D: 金属块涡流损耗 （一）启动W o r k b e n c h并保存 1.在windows系统下执行“开始”→“所有程序”→ANSYS 15.0→Workbench 15.0命令， 启动ANSYS Workbench 15.0，进入主界面。 2.进入Workbench后，单击工具栏中的 按钮，将文件保存。 （二）建立电磁分析 1.双击Workbench平台左侧的Toolbox→Analysis Systems→Maxwell 2D此时在Project Schematic中出现电磁分析流程图。 2.双击表A中的A2，进入Maxwell软件界面。在Maxwell软件界面可以完成有限元分析 的流程操作。 3.选择菜单栏中Maxwell 2D→Solution Type命令，弹出Solution Type对话框 (1)Geometry Mode:Cylinder about Z (2)Magnetic:Transient (3)单击OK按钮 4.依次单击Modeler→Units选项，弹出Set Model Units对话框，将单位设置成mm，并单 击OK按钮。 （三）建立几何模型和设置材料 1.选择菜单栏中Draw→Rectangle 命令，创建长方形 在绝对坐标栏中输入：X=500，Y=0，Z=0，并按Enter键 在相对坐标栏中输入：dX=20，dY=0，dZ=500，并按Enter键 2.选中长方形，选择菜单栏中Edit→Duplicate along line命令 在绝对坐标栏中输入：X=0，Y=0，Z=0，并按Enter键 在相对坐标栏中输入：dX=50，dY=0，dZ=0，并按Enter键 弹出Duplicate along line对话框，在对话框中Total Number:3，然后单击OK按钮。 3.选中3个长方形右击，在快捷菜单中选择Assign Material命令，在材料库中选择 Aluminum，然后单击OK按钮。 （四）设置求解域 选择菜单栏中Draw→Region命令，在弹出的Region对话框中输入Value=500，并单击OK按钮。 （五）添加激励 1.选中Rectangle1右击，在快捷菜单中选择Assign Excitations→Coil命令，弹出Coil Excitations对话框，在对话框中填入以下内容： (1)Name：CoilTerminal1 (2)Number of Conductors：100 (3)单击OK按钮 2.选中Rectangle1右击，在快捷菜单中选择Assign Excitations→Add Winding命令，弹出 Winding对话框，在对话框中填入以下内容： (1)Name:Winding_A (2)Type:Current (3)Stranded:?Checked (4)Current:50*sin(2*PI*50*Time)
ANSYS优化设计含几个实例
ANSYS 优化设计 1.认识ANSYS优化模块 1.1 什么时候我需要它的帮忙? 什么是ANSYS优化？我想说明一个例子要比我在这里对你絮叨半天容易理解的多。 注意过普通的水杯吗？底面圆圆的，上面加盖的哪一种。仔细观察一下，你会发现比较老式的此类水杯有一个共同特点：底面直径＝水杯高度。 图1 水杯的简化模型 为什么是这样呢？因为只有满足这个条件，才能在原料耗费最少的情况下使杯子的容积最大。在材料一定的情况下，如果水杯的底面积大，其高度必然就要小；如果高度变大了，底面积又大不了，如何调和这两者之间的矛盾？其实这恰恰就反应了一个完整的优化过程。 在这里，一个水杯的材料是一定的，所要优化的变量就是杯子底面的半径r和杯子的高度h，在ANSYS的优化模块里面把这些需要优化的变量叫做设计变量（DV）；优化的目标是要使整个水杯的容积最大，这个目标在ANSYS的优化过程里叫目标函数（OBJ）；再者，对设计变量的优化有一定的限制条件，比如说整个杯子的材料不变，这些限制条件在ANSYS 的优化模块中用状态变量（SV）来控制。下面我们就来看看ANSYS中怎么通过设定DV、SV、OBJ，利用优化模块求解以上问题。 首先参数化的建立一个分析文件（假设叫volu.inp），水杯初始半径为R＝1，高度为H ＝1（DV），由于水杯材料直接喝水杯的表面积有关系，这里假设水杯表面积不能大于100，这样就有S＝2πRH＋2πR2<100（SV），水杯的容积为V＝πR2H（OBJ）。 (用参数直接定义也可或者在命令栏内直接写) R=1 H=1 S=2*3.14*R*H+2*3.14*R*R V=10000/(3.14*R*R*H) 然后再建一个优化分析文件（假设叫optvolu.inp）,设定优化变量，并求解。 /clear,nostart /input,volu,inp /opt opanl,volu,inp opvar,R,dv,1,10,1e-2 opvar,H,dv,1,10,1e-2 opvar,S,sv,,100,1e-2 opvar,V,obj,,,1e-2 opkeep,on optype,subp opsave,optvolu,opt0 opexec 最后，打开Ansys6.1，在命令输入框中键入“/input,optvolu,inp”，整个优化过程就开始了。
ansys 蠕变例子
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 CREEP C*** STR. OF MATL., TIMOSHENKO, PART 2, 3RD ED., PAGE 531 ANTYPE,STATIC ET,1,LINK1 ! SPAR ELEMENT R,1,1,(1/30000) ! INITIAL STRAIN MP,EX,1,30E6 TB,CREEP,1 TBDATA,1,4.8E-30,7 ! CREEP PROPERTIES N,1 N,2,10 E,1,2 BFUNIF,TEMP,900 ! UNIFORM TEMPERATURE TIME,1000 KBC,1 D,ALL,ALL ! FIX ALL DOFS FINISH /SOLU SOLCONTROL,0 NSUBST,100 OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEP OUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEP SOLVE FINISH /POST26 ESOL,2,1,,LS,1,SIG ! STORE AXIAL STRESS PRVAR,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,3 LABEL(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,T950