! Build basic 2x2x60 tower finish /clear /prep7 blc4,-1,-1,2,2,60 et,1,185 mptemp,1,70 mpdata,ex,1,1,20e6 mpdata,nuxy,1,1,.23 mpdata,dens,1,1,.001 alphad,5 !Mesh esize,.9 vmesh,all /view,1,1,1,1 /vup,1,z eplot !create component of bottom nodes nsel,s,loc,z,0 cm,nbt,node nsel,all !Fix in UY and UZ, move UX 1" in first time step d,nbt,uy,0 d,nbt,uz,0 d,nbt,ux,1 time,.1 ! Set up Transient, Large Deflection Run antype,4 trnopt,full lumpm,0 nlgeom,1 kbc,0 timint,1 nsubst,40,40,40 outres,all,1 !Save first Load Step lswrite,1 !LS 2, UX = 5 inches d,nbt,ux,5 time,.2 lswrite,2 !LS 3, UX = 0 inches d,nbt,ux,0 time,.3 lswrite,3 !LS 3, no change, let it damp out time,2 lswrite,4 ! Solve all four load steps /solu lssolve ,1,4 finish ! Enter POST26 /POST26 NUMVAR,200 ! Max out number of allowed variables !SOLU,191,NCMIT !STORE,MERGE !FILLDATA,191,,,,1,1 !REALVAR,191,191 !!* ! Get UX on node 18 NSOL,2,18,U,X, UX_2 STORE,MERGE !plot UX results XVAR,1 PLVAR,2 ! Calculate total force on node 8 (sqrt-sum-of-squares) ! Put reaction forces in each direction in a variable (11-13) rforce,11,8,f,x,FX_8 rforce,12,8,f,y,FY_8 rforce,13,8,f,z,FZ_8 ! Square the variables prod,14,11,11 prod,15,12,12 prod,16,13,13 ! Add them together add,17,14,15,16 ! Square Root sqrt,18,17 ! Plot the result plvar,18 ! No find max stress on "edge" at 1,1 ! Grab the nodes nsel,s,loc,x,1 nsel,r,loc,y,-1 ! Get the node count *get,nnd,node,,count ! Get first node in selected set nn = ndnext(-1) ! Store its SIG1 values in var 20 ansol,20,nn,s,1 ! Get next node nn = ndnext(nn) !loop on remaining nodes *do,i,1,nnd-1 !Put SIG1 of current node in var 21 ansol,21,nn,s,1 ! Compare 20 and 21, and at each time step, but the larger value in 20 large,20,20,21 ! Move on to the next node nn = ndnext(nn) *enddo ! Plot the max values plvar,20 ! Use *get to get the maximum value over time, and when that time is *get,mxs,vari,20,extrem,vmax *get,mxt,vari,20,extrem,tmax ! Select all the nodes nsel,all /eof PLNSOL, U,SUM, 2,1.0 ANTIME,80,0.25, ,1,2,0,1 data,20,,,,ld_1,0 (F,F)