声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2372|回复: 2

再次请教ansys分析粘弹性问题

[复制链接]
发表于 2006-3-18 19:39 | 显示全部楼层 |阅读模式

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

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

x
<P>我现在做标准三参数固体的模型而ansys给的好象是maxwell模型。能否再次给一些详细的提示,或者我举个例子:矩形板一边固支,对边受常力,采用三参数固体模型。用ansys具体怎么做呀?<BR>谢谢</P>
回复
分享到:

使用道具 举报

发表于 2006-3-19 11:47 | 显示全部楼层

回复:(wenyue)再次请教ansys分析粘弹性问题

<P>ansys自带的vm200不就是相关的例子吗?<BR><BR>/COM,ANSYS MEDIA REL. 10.0 (05/31/2005) REF. VERIF. MANUAL: REL. 10.0<BR>/VERIFY,VM200<BR>JPGPRF,500,100,1             ! MACRO TO SET PREFS FOR JPEG PLOTS<BR>/SHOW,JPEG<BR>/PREP7<BR>/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS<BR>/COM,    ----- 2-D ANALYSIS -----<BR>/COM,    "FE CALCULATIONS OF RESIDUAL STRESSES ....",SOULES ET AL.<BR>ET,1,VISCO88,,,1            ! AXISYMMETRIC 2-D VISCOELASTIC ELEMENT<BR>*CREATE,MAC1                ! MACRO FOR MATERIAL PROPERTIES<BR>/COM,    MATERIAL ONE IS G-11 GLASS AND MATERIAL TWO IS ALUMINA<BR>/COM,    NOTE THAT ALUMINA IS AN ELASTIC MATERIAL THEREFORE IT<BR>/COM,     DOES NOT HAVE VISCOELASTICITY AND STRUCTURAL RELAXATION<BR>/COM,    MATERIAL PROPERTIES OF ALUMINA<BR>TB,EVISC,2<BR>TBDATA,31,52.6E-7           ! THERMAL EXPANSION COEFFICENTS FOR ALUMINA<BR>TBDATA,32,.119E-7           !    LOCATIONS 31-33<BR>TBDATA,33,-1.0E-11<BR>TBDATA,46,1.435E5           ! SHEAR MODULUS AT TIME = 0<BR>TBDATA,47,1.435E5           ! SHEAR MODULUS AT TIME = INFINITY<BR>TBDATA,48,3.11E5            ! BULK MODULUS AT TIME  = 0<BR>TBDATA,49,3.11E5            ! BULK MODULUS AT TIME  = INFINITY<BR>/COM,    MATERIAL PROPERTIES OF G-11 GLASS<BR>TB,EVISC,1<BR>TBDATA,1,6.45D4             ! H/R<BR>TBDATA,2,0.53               ! VALUE OF X<BR>TBDATA,3,6                  ! NO. OF MAXWELL ELEMENTS FOR STRUCTURAL RELAXATION<BR>TBDATA,6,0.108              ! COEFFICIENTS OF THE MAXWELL ELEMENTS FOR<BR>TBDATA,7,0.443              !    VOLUME DECAY (STRUCTURAL RELAXATION)<BR>TBDATA,8,0.166              !    LOCATIONS 6-11<BR>TBDATA,9,0.161<BR>TBDATA,10,0.046<BR>TBDATA,11,0.076<BR>TBDATA,16,3.00              ! RELAXATION TIMES FOR VOLUME DECAY FUNCTION<BR>TBDATA,17,0.671             !   LOCATIONS 16-21<BR>TBDATA,18,0.247<BR>TBDATA,19,0.091<BR>TBDATA,20,0.033<BR>TBDATA,21,0.008<BR>TBDATA,26,3.43E-5           ! THERMAL EXPANSION COEFF. IN LIQUID STATE<BR>TBDATA,31,64.7E-7           ! THERMAL EXPANSION COEFF. IN GLASS (SOLID) STATE<BR>TBDATA,32,.02E-7            !      LOCATIONS 31-32<BR>TBDATA,36,618               ! FICTIVE TEMPERATURE LOCATIONS 36-41<BR>TBDATA,37,618<BR>TBDATA,38,618<BR>TBDATA,39,618<BR>TBDATA,40,618<BR>TBDATA,41,618<BR>TBDATA,46,2.79D4            ! SHEAR MODULUS AT TIME = 0<BR>TBDATA,47,0.0D0             ! SHEAR MODULUS AT TIME = INFINITY<BR>TBDATA,48,6.05D4            ! BULK MODULUS AT TIME = 0<BR>TBDATA,49,6.05D4            ! BULK MODULUS AT TIME = INFINITY<BR>TBDATA,50,3                 ! THREE MAXWELL ELEMENTS FOR SHEAR RELAXATION<BR>TBDATA,51,0.422             ! COEFF. FOR SHEAR RELAXATION LOCATIONS 51-53<BR>TBDATA,52,0.423<BR>TBDATA,53,0.155<BR>TBDATA,61,0.0689            ! RELAXATION TIMES FOR SHEAR RELAXATION<BR>TBDATA,62,0.0065            !    LOCATIONS 61-63<BR>TBDATA,63,0.0001<BR>TBDATA,71,0                 ! NO BULK MODULUS RELAXATION<BR>*END<BR>*USE,MAC1                   ! EXECUTE MACRO FOR MATERIAL PROPERTIES<BR>:COM,   CREATE FINITE ELEMENT MODEL<BR>N,1,<BR>N,3,,0.00025<BR>FILL<BR>N,5,0,(0.00025+0.00325)<BR>FILL<BR>NGEN,3,10,1,5,1,.001<BR>MAT,2<BR>E,1,21,23,3,11,22,13,2<BR>MAT,1<BR>E,3,23,25,5,13,24,15,4<BR>/COM,    APPLY BOUNDARY CONDITIONS AND COUPLING<BR>NSEL,S,LOC,Y<BR>DSYM,SYMM,Y<BR>NSEL,S,LOC,X<BR>DSYM,SYMM,X<BR>NSEL,ALL<BR>D,1,ALL<BR>CP,1,UX,21,22,23,24,25<BR>CP,2,UY,2,22<BR>CP,3,UY,3,13,23<BR>CPSGEN,2,2,2,3,1<BR>FINISH<BR>/COM    SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A<BR>/COM     SCRATCH FILE<BR>/OUTPUT,SCRATCH<BR>*CREATE,MAC2                   ! CREATE MACRO FOR ANALYSIS TYPE AND LOADING<BR>/SOLU<BR>SOLCONTROL,0<BR>ANTYPE,STATIC<BR>/COM,    TEMPERATURE SET UP<BR>TREF,618<BR>TOFFST,273<BR>TUNIF,618<BR>TIME,1E-5<BR>CNVTOL,F,,,,.00001             ! VERY SMALL MINIMUM ENFORCED<BR>                               !  FOR CONVERGENCE<BR>SOLVE<BR>OUTRES,ESOL,1                  ! STORE RESULTS FOR EVERY SUBSTEP<BR>NSUBST,200<BR>TUNIF,460                      ! COOLING<BR>TIME,3160<BR>SOLVE<BR>TIME,(14400+3160)              ! ISOTHERMAL HOLD<BR>SOLVE<BR>TUNIF,18                       ! FURTHER COOLING<BR>TIME,(14400+12000)<BR>SOLVE<BR>*END<BR>FINISH<BR>*USE,MAC2                      ! EXECUTE ANALYSIS AND LOADING MACRO<BR>/OUTPUT<BR>/POST26<BR>ESOL,2,2,,BFE,TEMP<BR>ESOL,3,2,3,S,X,STRESS<BR>*CREATE,MAC3                   ! MACRO FOR PROCESSING RESULTS<BR>XVAR,2<BR>/GRID,1<BR>/AXLAB,X,TEMPERATURE<BR>/AXLAB,Y,IN-PLANE STRESS (MPA)<BR>PLVAR,3<BR>*GET,MXSX,VARI,3,EXTREM,VMAX   ! MAXIMUM IN-PLANE STRESS<BR>NSTORE,20                      ! STORE EVERY 20TH TIME POINT RESULTS<BR>PRVAR,2,3<BR>*END<BR>*USE,MAC3                      ! EXECUTE POSTPROCESSING MACRO<BR>*SET,P1,(MXSX)<BR>*GET,T1,VARI,3,EXTREM,TMAX<BR>*GET,TE,VARI,2,RTIME,T1<BR>*DIM,LABEL,CHAR,2,2<BR>*DIM,VALUE,,2,3<BR>LABEL(1,1) = 'PRES MX ','TEMP '<BR>LABEL(1,2) = 'MPa','DEG C'<BR>*VFILL,VALUE(1,1),DATA,12.5,460<BR>*VFILL,VALUE(1,2),DATA,P1,TE<BR>*VFILL,VALUE(1,3),DATA,ABS(P1/12.5),ABS(TE/460)<BR>SAVE,TABLE_1<BR>FINISH<BR>/CLEAR,NOSTART                 ! CLEAR THE DATABASE<BR>JPGPRF,500,100,1             ! MACRO TO SET PREFS FOR JPEG PLOTS<BR>/PREP7<BR>/TITLE, VM200, VISCOELASTIC SANDWICH SEAL ANALYSIS<BR>/COM,    ----- 3-D ANALYSIS -----<BR>ET,1,89                        ! 3-D VISCOELASTIC ELEMENT<BR>*USE,MAC1                      ! EXECUTE MACRO FOR MATERIAL PROPERTIES<BR>/COM,   CREATE FINITE ELEMENT MODEL<BR>N,1,<BR>N,3,0.00025<BR>FILL<BR>N,5,(0.00025+0.00325)<BR>FILL<BR>NGEN,3,10,1,5,1,,.001<BR>NGEN,3,100,1,25,1,,,0.001<BR>MAT,2<BR>E,1,3,23,21,201,203,223,221<BR>EMORE,2,13,22,11,202,213,222,211<BR>EMORE,101,103,123,121<BR>EGEN,2,2,1,1,1,-1<BR>NSLE,S<BR>NSEL,INVE<BR>NDELE,ALL<BR>NSLE,S<BR>/COM,    APPLY BOUNDARY CONDITIONS AND COUPLING<BR>NSEL,S,LOC,Y<BR>DSYM,SYMM,Y<BR>NSEL,S,LOC,X<BR>DSYM,SYMM,X<BR>NSEL,S,LOC,Z<BR>DSYM,SYMM,Z<BR>NSEL,S,LOC,Y,0.002<BR>CP,1,UY,ALL<BR>NSEL,S,LOC,Z,0.002<BR>CP,2,UZ,ALL<BR>NSEL,S,LOC,X,0.00025<BR>CP,3,UX,ALL<BR>NSEL,S,LOC,X,0.0035<BR>CP,4,UX,ALL<BR>NSEL,ALL<BR>FINISH<BR>/COM,    SINCE THE SOLUTION OUTPUT IS VOLUMINOUS IT IS DIVERTED TO A<BR>/COM,     SCRATCH FILE<BR>/OUTPUT,SCRATCH<BR>*USE,MAC2                       ! EXECUTE ANALYSIS AND LOADING MACRO<BR>/OUTPUT<BR>/POST26<BR>ESOL,2,2,,BFE,TEMP<BR>ESOL,3,2,3,S,Y,STRESS<BR>*USE,MAC3                       ! EXECUTE POSTPROCESSING MACRO<BR>*SET,P2,(MXSX)<BR>*GET,T2,VARI,3,EXTREM,TMAX<BR>*GET,TE2,VARI,2,RTIME,T2<BR>*DIM,LABEL,CHAR,2,2<BR>*DIM,VALUE,,2,3<BR>LABEL(1,1) = 'PRES MX ','TEMP '<BR>LABEL(1,2) = 'MPa','DEG C'<BR>*VFILL,VALUE(1,1),DATA,12.5,460<BR>*VFILL,VALUE(1,2),DATA,P2,TE2<BR>*VFILL,VALUE(1,3),DATA,ABS(P2/12.5),ABS(TE2/460)<BR>SAVE,TABLE_2<BR>RESUME,TABLE_1<BR>/COM<BR>/OUT,vm200,vrt<BR>/COM,------------------- VM200 RESULTS COMPARISON --------------<BR>/COM,<BR>/COM,                 |   TARGET   |   ANSYS   |   RATIO<BR>/COM,<BR>/COM,RESULTS USING VISCO88<BR>*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)<BR>(1X,A8,A8,'   ',F10.1,'  ',F10.1,'   ',1F5.3)<BR>/NOPR<BR>RESUME,TABLE_2<BR>/GOPR<BR>/COM,<BR>/COM,RESULTS USING VISCO89<BR>*VWRITE,LABEL(1,1),LABEL(1,2),VALUE(1,1),VALUE(1,2),VALUE(1,3)<BR>(1X,A8,A8,'   ',F10.1,'  ',F10.1,'   ',1F5.3)<BR>/COM,-----------------------------------------------------------<BR>/OUT<BR>FINISH<BR>/NOPR<BR>/DELETE,TABLE_1<BR>/DELETE,TABLE_2<BR>/DELETE,MAC1<BR>/DELETE,MAC2<BR>/DELETE,MAC3<BR>FINISH<BR>*LIST,vm200,vrt<BR></P>
发表于 2006-3-19 14:45 | 显示全部楼层
<P>     哇,高手,这些看着好像有点头痛,我也有相似的问题想问您,我也是想用三元件模型来定义模型的材料属性,但是不知如何ansys输入,因为材料的粘弹性模型它只有Maxwell模型呀.请高手指点!!!!</P>
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-4-29 19:41 , Processed in 0.102772 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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