/filname.plane42/title,plane predictionAinits、si/PREP7!Defines a local element type fiom tlie element library. !3 - Plane stress with thickness input!单位为mmet,l,plane42,n31；1A2*dinKdens01600*doja, 1,1600!单位为MPaDENS0(ia=0.800 exll=DENS0(ia)*2*100 ?DENS0(ia 尸 1.0?exl l=dens0(ia)*2*5000 mpexja.exllMP,niiXYd久0.3*enddojj 1=20k、2,l k、3 丄 1 kA0,lU,21.2.31.3.41,4,1al, 1,2,3,4I*lesizeJ,40 lesize,2,40 lesize,3”,40 lesize,4,40 amesh.allFINISH/SOLU!FLST21 丄 ORDE,1!FITEM21/GO!D.P51X, ,UX,UY, D丄，,UX,UY,FLST,2,40 丄 ORDE,2 FITEM22FITEM.2,-41!dojl,l,jjl /postl !jja=jl/GOD.P51X,.UY, !D,2, ,UY,，八， !*SFL，3,PRES, 2 !SFLSPRES，2000000 !SFLSPRES，40000 SFTRAN!E100.FY0l*getaielenLeleniMnunLniax *dinKesednelem *dmKesedlMnelem *dmKesed2Mnelem *dniKesed 1 a.jielem *dmKdenslnelem *dim.detMnelem *dinKexlnelem *dmKdens2nelem *dmKStialvnelem !*dmLniaxstraJjl*dojill Jjielem!此处应该是将其初始化一下!DENS 1(1111 )=DENSO(iill) esed )=0 esedla(ii 11)=0*enddoeqsl;iter!开始J1循坏SOLVE fiiushetable,stre,s,eqvetable,stra,epto,eqvsmult,esed,stie,stra,0.5,lI*! * vget, STRA1 (1 ),elem, 1 ,etab,STRA严DO JI丄40! STRA 1 (1 +(JI-1) *40)=0.000001! STRA 1 (40+( JI J )*40尸0.000001 !*ENDDO元内的应变值设为0.00000b why?!MAXXJJ=STRA1(1)严DO, JU ,41,1560! STRA 1 (JI 1 )=ABS(STRA 1 (JI 1)! *IF,STRA 1 (JI 1 人GT.MAXXJ J.THEN!MAXXJJ=STRAl(ni)严ENDIF!*ENDDO! MAXSTRA(J1 )=MAXXJJt*vget,det( 1 ),elem,l ,etab.esed数组-1600维! * vget,esed2 (1 ),elem, 1 ,etab,esed*dojil Jaielem*SET,esed 1 (iil ),(det(iil)/dens 1 (ii1)-0.025)! *SET,esed 1 (iil),esed2(l )/dens l(iil)!赋1, 41, 81, 1561号单元应变。 !赋40, 80, 120, 1600 号单元应变。!以上两行即表示把平板左右两边单!取41到1560号单元的应变的绝对值。!找每一 J1循环的最人应变。!从单元1开始，把单元表中的应变能赋给det!esed2单位为MPa!此公式代表激励与密度的关系，esedl单位为J/gesed la(iil )=esed l(iil) *enddo!相邻两个单元中心的间距。mm!*影响系数Dl=l/40!DD=0.001*0.025DD=0.045!影响距离。!单位：nunDA=D1/DD*set,si,exp(-DA)*set,sj,exp(-DA*sqit(2)* set,sk,exp(-2 *DA)*set,sLexp(-DA*sqit(5)*set,sm,exp(-DA*sqrt(8)!垂直和水平方向相邻两个单元间的影响系数。!45角方向相邻两个单元间的影响系数。!垂直和水平方向相隔一个单元的两个单元间的影响系数。 !例如：82号和43号单元对1号单元的影响系数。!例如：83号单元对1号单元的影响系数。*do0丄40*ii0,eq,l Jhen 单元。*do,il?L40 iz=(i0-1)*40+11 *ifjl,eqj .then swa=esedl(iz)!开始10循坏!开始10, if语句，eq-Equal ,此处开始计算第一行的40个!开始11循坏!开始订，if语句，计算第一行第一个单元swb=(esed 1 (iz+40)+esed 1 (iz+ l)*si swc=(esed 1 (iz+41 )*sjswd=(esed 1 (iz+2)+esedl (iz+80)*sk swe=(esed 1 (iz+42)+esed 1 (iz+8 l)*sl swf=(es 亡 d 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf*elseiil,eq,2!第二行第二个单元swa=esedl(iz)swb=(esed 1 (iz-1 )+esed 1 (iz+1 )+esed 1 (iz+40)*siswc=(esed 1 (iz+41 )+esed 1 (iz+39)*sjswd=(esed 1 (iz+2)+esedl (iz+80)*skswe=(esed 1 (iz+42)+esed 1 (iz+79)+esed 1 (iz+81 )*sl swf=(es 亡 d 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf*elseif,il Je,38!边上中间的36个单元。le为小于等于swa=esedl(iz)swb=(esed 1 (iz-1 )+esed 1 (iz+1 )+esed 1 (iz+40)*siswc=(esed 1 (iz+39)+esed 1 (iz+4 l)*sjswd=(esed 1 (iz-2)+esed 1 (iz+2)+esed 1 (iz+80)*skswe=(esed 1 (iz+3 8)+esed 1 (iz+42)+esed 1 (iz+79)+esed 1 (iz+8 l)*sl swf=(esed 1 (iz+78)+esed 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf*elseifjl,eq,39!39号单元swa=esedl(iz)swb=(esed 1 (iz-1 )+esed 1 (iz+1 )+esed 1 (iz+40)*si swc=(esed 1 (iz+39)+esed 1 (iz+4 l)*sjswd=(esed 1 (iz-2)+esed 1 (iz+80)*skswe=(esed 1 (iz+3 8)+esed 1 (iz+79)+esed 1 (iz+81 )*sl swf=(esedl(iz+78)*smesed 1 a(iz)=swa+swb-rswc+swd+swe+swf*else!40号单元swa=esedl(iz)swb=(esed 1 (iz-1 )+esed 1 (iz+40) *siswc=(esed 1 (iz+39)*sjswd=(esed 1 (iz-2)+esed 1 (iz+80)*sk swe=(esed 1 (iz+38)+esed 1 (iz+79)*sl swf=(esed 1 (iz+78)*smesed 1 a(iz)=swa+swb-rswc+swd+swe+swf*endif!结束il, if语句。*enddo*elseifj0,eq.2 *do,il,L40 iz=(iO-1)*40+11 *ifjl.eq.ljhen swa=esedl(iz)!结束il循环。!计算第二行单元!41号单元swb=(esed 1 (iz-40)+esed 1 (iz+1 )+esedl (iz+40)*siswc=(esed 1 (iz-39)+esed 1 (iz+4 l)*sjswd=(esed 1 (iz+2)+esedl (iz+80)*slswe=(esed 1 (iz-3 8)+esed 1 (iz+42)+esed 1 (iz+8 l)*sl swf=(esed 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf *elseif,il,eq,2!42 号单元swa=esedl(iz)swb=(esed 1 (iz-40)+esed 1 (iz-1 )+esed 1 (iz+1 )+esed 1 (iz+40) *si swc=(esed 1 (iz-39)+esed 1 (iz-41 )+esed 1 (iz+39)+esed 1 (iz+4 l)*sj swd=(esed 1 (iz+2)+esedl (iz+80)*skswe=(esed 1 (iz-3 8)+esed 1 (iz+42)+esed 1 (iz+7 9)+esed 1 (iz+8 l)*sl swf=(esed 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf*elseif,ilJe,38!第二行中间36个单元。swa=esedl(iz)swb=(esed 1 (iz-40)+esed 1 (iz-1 )+esed 1 (iz+1 )+esed 1 (iz+40) *siswc=(esed 1 (iz-39)+esed 1 (iz-41 )+esed 1 (iz+39)+esed 1 (iz+4 l)*sjswd=(esed 1 (iz-2)+esed 1 (iz+2)+esed 1 (iz+80)*skswe=(esed 1 (iz-42)+esed 1 (iz-38)+esed 1 (iz+38)+esed 1 (iz+42)+es 亡 d 1 (iz+79)+esed 1 (iz+8 l)*sl swf=(esed 1 (iz+78)+esed 1 (iz+82)*smesed 1 a(iz)=swa+swb+swc+swd+sw 亡+swf*elseif,il,eq,39!79 号单元swa=esedl(iz)