Документ взят из кэша поисковой машины. Адрес оригинального документа : http://www.apo.nmsu.edu/Telescopes/SDSS/eng.papers/19990112_PlugPlateDistortion/pl124.fea Дата изменения: Wed Oct 31 21:02:15 2001 Дата индексирования: Sun Apr 10 08:48:59 2016 Кодировка: Поисковые слова: п п п п п п п п п п п п п п п п п р п

c***c*** Initialize variablesc*** SLOP = 0.37254E-01 OUTY = 0.35241E-02 INY = 0.10442E-02c*** INY: 0.001 added to INY since design variable must be greater than 0.c***/PREP7C***C*** The next line is the only one which needs the run number changed./TITLE, pl124c***C*** AN ELEMENT ALWAYS HAS IDENTICAL MAT, TYPE, & REAL TABLE NOS.:C***C***C***C*** TABLE DESCRIPTIONC***C*** 1 plug-plateC***C*** 2-d isoparametric solid, large deflection, axisymmetricET,1,plane42,1,0,1,0,0*repeat,3,1,0C***C***mp,NUXY,1,.334 ! Poisson's ratio - Aluminum*repeat,3,,1,0mp,ALPX,1,23.2 ! (deg C)^-1 - Aluminum*repeat,3,,1,0mp,DENS,1,2770 ! kg/m^3 - Aluminum*repeat,3,,1,0mp,EX,1,71.0E9 ! N/m^2 - Aluminum*repeat,3,,1,0C***thic = 0.00318c*** radi is actually 0.3975 (31.3 inches).radi = 0.3937ndbg = 11ndb1 = ndbg+1ndif = 100ndum = 90elum = ndum-1nlay = 6elay = nlay-1ntbg = (elay*ndif)+ndbgorin = 1orad = 0.39243irin = 14irad = 0.330945delx = orad-iraddely = slop*delxC***C*** ***************** NODES **************************C***csysn,ndbg,radi,-thicn,ndbg+orin,orad,-thicfill,,,,,,,,1n,ndbg+irin,irad,-thicfill,,,,,,,,1n,ndbg+elum,0,-thicfill,,,,,,,,1n,ntbg,radi,0n,ntbg+orin,orad,0fill,,,,,,,,1n,ntbg+irin,irad,0fill,,,,,,,,1n,ntbg+elum,0,0fill,,,,,,,,1fill,ndbg,ntbg,nlay-2,,,ndum,1c***c*** ******************** Elements **********************************c***type,1mat,1e,ndbg,ndb1,ndb1+ndif,ndbg+ndifegen,elum,1,-1egen,elay,ndif,-elumcsyssavefinishc***C*** ********************* Solution *********************************C***/solutionantype,staticnlgeom,onc*** sstif,onc*** deltim,.5,.05c*** autots,onoutpr,all,lastcnvtol,u,,.001c***C*** *************************** LOADS *******************************C***C*** Axial constraintscsysd,ndbg+elum,uy,iny-0.001d,ndbg+orin,uy,outy+delyd,ntbg+irin,uy,outyc***d,ndbg+elum,ux,0*repeat,6,ndif !rp nlayc***solvefinishc***c*** Enter post processor and calculate the variance of the error weightedc*** with area.c***/post1setnsel,s,node,,44fsum*get,f01,fsum,fynsel,s,node,,77fsum*get,f02,fsum,fynallf01=-f01f02=-f02c***aa00 = -1.468e-003 aa02 = -3.104e-006 aa04 = -4.111e-010 aa06 = 1.377e-015 c***ab00 = 35.881ab02 = -0.00073173ab04 = 1.3682e-08ab06 = -4.0372e-14c***vsum=0.0wsum=0.0*get,y1,node,528,u,yy1=(1000.0*y1)*get,t1,node,528,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,532,u,yy1=(1000.0*y1)*get,t1,node,532,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,536,u,yy1=(1000.0*y1)*get,t1,node,536,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,540,u,yy1=(1000.0*y1)*get,t1,node,540,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,545,u,yy1=(1000.0*y1)*get,t1,node,545,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,550,u,yy1=(1000.0*y1)*get,t1,node,550,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,556,u,yy1=(1000.0*y1)*get,t1,node,556,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,563,u,yy1=(1000.0*y1)*get,t1,node,563,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,571,u,yy1=(1000.0*y1)*get,t1,node,571,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempc****get,y1,node,583,u,yy1=(1000.0*y1)*get,t1,node,583,loc,xt1=(1000.0*t1)t2=t1*t1y2=t2*(aa04+(aa06*t2))y2=-(aa00+(t2*(aa02+y2)))w1=ab00+(ab02+(ab04+ab06*t2)*t2)*t2w1=w1*w1wsum=wsum+w1temp=w1*(y1-y2)*(y1-y2)vsum=vsum+tempvsum=vsum/wsumrms=sqrt(vsum)*statusfinish/eofc***c*** Enter optimizer and minimize vsum, the variance of the errorc***/optopvar,slop,dv,0.035,0.045,0.000000001opvar,outy,dv,0.0031,0.0044,0.000000001opvar,iny,dv,0.0006,0.0014,0.000000001c***opvar,vsum,obj,,,1.0e-12oplistopcopyoprun,49oplist,allfinish