初始应力添加与多工况任意组合实例.doc

一） 初始应力场施加实例实例1要在分层的SHELL181单元中的一指定层上施加初始应力,可以使用指定层的材料ID号 　　本文介绍了ANSYS施加初始应力场的相关方法　　在用ANSYS做结构分析时,可以把初始应力指定为一项载荷，只在静态分析和全瞬态分析中被承认(分析可以是线性,也可以是非线性), 初始应力载荷只能在分析的第一载荷步中施加.恒应力可以用ISTRESS命令指定.运用ISFILE命令从输入文件中获取初始应力值,为了列表或是删除初始应力,命令变元允许初始应力被限制为具体的材料类型.　　要在分层的SHELL181单元中的一指定层上施加初始应力,可以使用指定层的材料ID号在初始应力命令中使用MATx 字段可以根据层的材料ID号施加需要的应力，或者可以在SHELL181单元的所有层上施加初始应力，然后删除那些不希望有初始应力的层上的初始应力用delete, elementID, 层号命令，还可以用初始应力命令列表或删除分层SHELL181单元上指定层上的初始应力值　　初始应力功能只支持以下单元类型:PLANE2, PLANE42, SOLID45, PLANE82, SOLID92, SOLID95, LINK180, SHELL181, PLANE182, PLANE183, SOLID185, SOLID186, SOLID187, BEAM188,还有BEAM189.ISFILE命令可以将合成初始应力写入文件。

这些命令可以在/SOLUTION处理器中运用，在ANSYS Commands Reference(ANSYS命令参考手册)中可以看到所有初始应力命令的描述　　【1】首先将初始应力数据用一个外部的ASII文件给出，假设文明名为istress.ist，内容如下　　/CSYS,0　　! ELEM ID ELEM INTG LAY/CELL SECT INTG SX SY SZ SXY SYZ SX　　1,all,all,all,100,0,0,0,0,0　　2,all,all,all,100,0,0,40,0,0　　3,all,all,all,100,0,0,0,0,0　　4,all,all,all,100,0,0,0,0,40　　5,all,all,all,100,20,0,0,0,0　　【2】读取应力场文件　　/prep7　　/title, Example of Initial stress import into ANSYS　　et,1,182　　! Plane stress PLANE182 element　　mp,ex,1,1.0e9　　mp,nuxy,1,0.3　　!　　! Define the nodes　　!　　n,1　　n,2,2.0　　n,3,4.0　　n,4,6.0　　n,5,8.0　　n,6,10.0　　n,7,,1.0　　n,8,2.0,1.0　　n,9,4.0,1.0　　n,10,6.0,1.0　　n,11,8.0,1.0　　n,12,10.0,1.0　　!　　! Define the 5 elements　　!　　e,1,2,8,7　　e,2,3,9,8　　e,3,4,10,9　　e,4,5,11,10　　e,5,6,12,11　　! Constrain all dofs on all nodes at x=0 to be zero　　nsel,s,loc,x,　　d,all,all　　nall　　finish　　!　　/solu　　!读应力文件　　inis,read,istress,ist　　! List the initial stresses　　inis,list　　outres,all,all　　solve　　finish　　!　　/post1　　set,last　　PLNSOL,S,EQV,0,1.0　　finish实例2 一个岩土工程中常用的岩层简化平面模型同时受重力和地表均压作用，比较以下两种方法所获的结果。

1.常规方法，同时施加重力加速度和均压 2.初应力载荷法，先计算模型在重力作用下的结果，同时输出初应力文件，然后在全新的模型中调入此初应力文件，计算等效节点力反力；施加等效节点力，进行求解，得到与只有重力场作用下的模型同样的受力和变形;最后施加均压，求解finish /filename,try01,1/prep7et,1,42mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800blc4,,,1,2esize,0.1amesh,all /solu nsel,s,loc,y,0d,all,uy,0nsel,s,loc,x,0nsel,a,loc,x,1d,all,ux,0allsel acel,,10,inistate,write,1,,ssolve !得到只有重力作用下的初应力文件finish /clear,nostart /filename,try02,1/prep7et,1,42mp,ex,1,2e11mp,nuxy,1,0.3mp,dens,1,7800blc4,,,1,2esize,0.1amesh,all /solu d,all,all isfile,read,try01,ist,,2isfile,list outres,all,all allsel solve !计算等效节点力反力finish/filename,try03,1/solu lsclear,all nsel,s,loc,y,0d,all,uy,0nsel,s,loc,x,0nsel,a,loc,x,1d,all,ux,0allsel ldread,reac,1,1,,,'try02','rst'solvesfl,3,pres,1e5allsel solve/post1plnsol,u,y 实例3ANSYS中关于初应力和荷载步设置的算例算例为一个地基+地基上面的一块方块墙吧，先通过一次计算仅加边界条件和自重，计算得到自重应力场，并输出初应力文件用来模拟初应力场。

然后施加应力和自重并进行计算，此时的位移基本为0，即消除了初位移，所得应力场即为自重应力场在些基础上进一步施加墙上法向面荷载，并进行第二步计算，得到的位移应该是仅有面荷载引起的位移关于初应力想说明的几点：1、初应力只能加在第一个荷载步，用命令流和GUI方式均可但在求解前不能退出求解器，如果加了初应力后退出求解器到前处理器或者后处理器后再回到前处理器，刚才施加的初应力就没有了，必须再次施加2、也可以不用初应力而直接分两个荷载步进行计算，如第一步仅计算自重，第二步再加面荷载后进行计算，在后处理中用工况组合来得到“净位移”，但工况组合中的应力结果似乎是不正确的3、如果用LSSOLVE从荷载文件进行求解，在写荷载文件时初应力的设置并不会写入荷载文件，所以，在命令流或者GUI方式下在求解前必须显示指定加载初应力关于荷载步设置的几点建议：1、在GUI方式下，每次进入求解器进行求解似乎都是开始一个新的分析（这一点偶也不是很明白）如果不退出求解器，即便不改变约束和荷载，只要求解一次，就会多一个荷载步结果，但所有结果是一样的；如果退出求解器后再进来，求解就重新开始（根据时间值）2、对点、线、面、体荷载都有替换和叠加两种方式，在替换方式下，在同一位置重复加荷载，只有最后一次加的荷载有效；在叠加方式下，在同一位置重复加荷载，所有荷载会叠加后共同作用在结构上。

/BATCH  KEYW,PR_SET,1   KEYW,PR_STRUC,1KEYW,PR_THERM,0KEYW,PR_FLUID,0KEYW,PR_ELMAG,0KEYW,MAGNOD,0   KEYW,MAGEDG,0   KEYW,MAGHFE,0   KEYW,MAGELC,0   KEYW,PR_MULTI,0KEYW,PR_CFD,0   /GO!*  /prep7!* define material protertiesmp,dens,1,2490mp,ex,1,20e9mp,prxy,1,.167  mp,dens,2,2660  mp,ex,2,10e9mp,prxy,2,.28!* define element typeet,1,solid45!* define geometryblock,-5,-0.5,-5,5,-5,0block,-0.5,0.5,-5,5,-5,0block,0.5,5,-5,5,-5,0block,-0.5,0.5,-5,5,0,5vglue,allvsel,s,loc,z,-6,0vatt,1,,1,0vsel,invevatt,2,,1,0!* mesh geometryESIZE,0,10lsel,s,loc,x,-0.4,0.4lesize,all,,,4allsel,allvmesh,allASEL,S,LOC,Z,-6,-0.1ASEL,U,LOC,X,-0.6,-0.4  ASEL,U,LOC,X,0.4,0.6cm,uxyz,areaasel,s,loc,z,0,5asel,r,loc,x,0.4,0.6cm,sfload,areaALLSEL,ALL  FINISH  /SOLU   !*  ANTYPE,0!* the zero load step only for get initial stress filelsclear,allTIME,1  AUTOTS,-1   NSUBST,1, , ,1  KBC,0!* define constrainscmsel,s,uxyzda,all,all,0.!* deifine gravityacel,0.,0.,10.allsel,alliswrite,on   !* open option for write initial stress filesolveFINISH  /SOLU   !*  ANTYPE,0!* the first load steplsclear,allTIME,1  AUTOTS,-1   NSUBST,1, , ,1  KBC,0!* define constrainscmsel,s,uxyzda,all,all,0.!* deifine gravityacel,0.,0.,10.allsel,allISFILE,READ,AINIStress,ist, ,1  !*初应力仅在第一个荷载步施加iswrite,offsolve!* the second load steplsclear,allTIME,2  AUTOTS,-1   NSUBST,1, , ,1  KBC,0!* define constrainscmsel,s,uxyz。