Part 21 Report - 1996-040

ACCESSION #: 9601190322 2060 Fairfax Avenue, Cherry Hill, NJ 08003-1666 Telephone: (609) 424-0999 H O L T E C Telex: 910-240-6663 INTERNATIONAL Telefax: (609) 424-1710 October 16, 1995 Nuclear Regulatory Commission Office of Inspection and Enforcement 475 Allendale Road King of Prussia, Pennsylvania 19406 NOV 16 1995 Subject: 10CFR21 Notification Gentlemen: We herein enclose internal memorandums issued by Dr. Alan Soler, Executive Vice President, and Dr. James Bosnik, Principal Engineer. These memorandums document discrepant results obtained during nonlinear transient analyses using the finite element code ANSYS, versions 5.0A and 5.1. Please be advised that two distinct problems have been identified: 1. Disparate results for calculated displacements are obtained depending on the computer employed (described in attached memorandum from A. Soler to V. Gupta). 2. Disparate results for forces and displacements are obtained depending upon the technique employed to launch the analysis runs (described in attached memorandum from J. Bosnik to V. Gupta). The magnitude of the discrepancies in both cases is considered significant and should not occur. The originator of the finite element code, ANSYS Inc., has been advised of these disparities and has stated that the matter is under investigation. It is indefinite whether ANSYS plans to issue a general notice regarding this issue. The differences documented in this notification have no affect with respect to safety-related equipment or services previously supplied to licensees by Holtec International. Sincerely, * Received response from ANSYS Inc. dated 10/18/95. This letter, and A. Soler's counterresponse, is attached. V. Gupta 10/23/95 Vik Gupta, Quality Assurance Manager H O L T E C INTERNATIONAL MEMORANDUM TO: Vik Gupta, Quality Assurance Manager FROM: Alan Soler, Executive Vice President RE: Non-Linear Time History Analysis DATE: October 12, 1995 1. In the process of performing a non-linear transient analysis using ANSYS 5.0A, using the recommended default values for convergence, solution methodology, two distinct sets of displacement results occurred, depending on the computer used. Rerunning the problem using ANSYS 5.1 produced a third set of results. Attachment A, which was sent to ANSYS, Inc., outlines the problem. 2. ANSYS will issue an answer to us shortly which claims roundoff error as the problem. I am, however, not convinced and believe that the problem may lie with the default iteration scheme used. If different computers were truly the answer, then all systems would give different results. 3. In our application, forces were insensitive to the solution so that the analysis results could be used. However, in a problem where displacements are important, then the lack of convergence is a problem. 4. Please issue a 10CFR21 notice. Alan Soler Executive Vice President AS:mp Document ID: 40254AW 2060 Fairfax Avenue, Cherry Hill, NJ 08003-1666 Telephone: (609) 424-0999 H O L T E C Telex: 910-240-6663 INTERNATIONAL Telefax: (609) 424-1710 August 4, 1995 QA Manager, SASI Swanson Analysis Systems, Inc. P.O. Box 65 Houston, PA 15324-9909 Dear Sir, The attached information concerns a problem we are having with ANSYS which gives two distinct sets of answers on DOS platforms and a third set under Windows NT. We need to know whether we have simply a convergence problem (and that all solutions would be the same on all platforms if we reduce tolerances to the minim:n-rn) or whether there is a more fundamental problem. Sincerely, Alan I. Soler Exercutive Vice President AIS:mp Attachment Document ID: 40254AP QA NOTIFICATION TO SASI H-40254 To: QA Manager, SASI From Dr. Alan L Soler, Executive Vice President, Holtec International Subject: Inconsistent ANSYS results on a Nuclear Safety Related Analysis Date:August 4, 1995 By this letter, we are informing you of a possible problem with ANSYS dealing with the transient analysis of a structure involving substructuring, linear, and non-linear elements. The problem we are having involves the use of a substructure generation pass, a time history analysis using seismic acceleration files input as a gravity vector at each step in the 20 sec. solution, and a subsequent post-processing to find maximum displacements, and spring forces. Attachment 1 is a listing of the generation pass and the use pass which contains both the solve block and the post block. Also enclosed is a disc containing the script files, the input seismic loadings, and certain macros called during the post processing. The problem encountered is that the script has been tried on 5 different commuters at two different locations using ANSYS 5.0A and 5.0A-95 under DOS and results in two distinct sets of answers on the 5 platforms. In addition, the problem has been run on ANSYS5.1 on one of the computers and gives a third set of results. Mallet Technology in Malvern Pa. is attempting to run the problem on their UNIX workstation but I do not yet have any results from them. We believe that we need your review of our results to ascertain whether there is a real problem or simply a problem with our script file that is causing a loss of control on different DOS platforms. I believe the differences in results are too far apart to be attributable to round off error in the individual CPU's. We are enclosing summaries of the different results along with the characteristics of the computers involved. COMPUTERS AT HOLTEC - all pentium's with "good" chip pentium-1 90 MHZ, 64MB ram, 2GB SCSI Hard Drive, no disk compression, running DOS 6.22, or Windows NT 3.5 (workstation) pentium-2 90 MHZ, 32MB ram, 1GB IDE Hard Drive, no disk compression, MSDOS 6.22 pentium-3 90 MHZ, 16MB ram, 540 IDE Hard Drive, Stacker 4.0 disk compression, DOS 6.22 August 4, 1995 1 QA NOTIFICATION TO SASI H-40254 COMPUTERS AT AEC, San Francisco (3rd party independent reviewer) pentium 100 (2 different chips gave same result) 100 MHZ, 64MB ram, ? Hard Drive, other features unknown 486 machine of unknown vintage and components For your information, I am enclosing the autoexec.bat and config.sys files from the Holtec computers along with the summary page created by the script after the postprocessing operations. Please note that the script file uses a tolerance of .005 on the force only but this tolerance has been used consistently on ad platforms. Also note that the solutions on the same computer 5.0A95 and 5.1 are not in reasonable agreement We have rerun the scripts using the default tolerances of .001 on force and moment and notice some differences in the results but the platform differences still remain. I would expect that even if the tolerance setting was not small enough, there should still be a consistent solution across platforms. Since this project is safety related and falls within the purview of 10CFR21, we need your help in resolving this issue. August 4, 1995 2 Winco rack analysis, 6x5x10, transient, use pass, full, MU=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1use.sum. Maximum vertical load in any single pedestal: 55952.0425 lb. Maximum total vertical pedestal load: 111457.107 lb. Maximum shear load in any single pedestal: 44078.667 lb. Maximum fuel-box contact load at any one location: 4981-69403 lb. Maximum rack/rack vertical intertie load at any one location: 54532.7439 lb. Maximum rack/rack horizontal intertie load at any one location: 4176.09255 lb. Maximum corner x-displacement, top plate: 0.470916033 in. Maximum corner x-displacement, baseplate: 0.285784018 in. Maximum corner y-displacement, top plate: 0.629147025 in. Maximum corner y-displacement, baseplate: 0.338490697 in. Pedestal stress factor r1: 9.531621629E-02. Pedestal stress factor r2: 0.11047849. Pedestal stress factor r3: 0.136828741. Pedestal stress factor r4: 0.152901297. Pedestal stress factor r5: 0.223969783. Pedestal stress factor r6: 0.225287353. Corresponding rip = 9.38618944E-02, r2p = 0.11047849, r3p = 3.023009224E-02, r4p - 0.152901297. Pedestal stress factor r7: 9.88653008E-02. Winco rack analysis, 6x5x10, transient, exp pass, full, mu=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1exp.sum. /OUTPUT FILE= rac1exp.lst Box stress factor r1: 0.434341914 (elem = 240, lstep = 149, node = 1). Box stress factor r2: 1.654037667E-02 (elem = 216, 756, node = 1). Box stress factor r3: 0.266888685 (elem = 240, lstep = 604, node = 2). Box stress factor r4: 0.281520977 (elem = 216, lstep = 756, node = 2). Box stress factor r5: 0.632717929 (elem = 240, lstep = 149, node = 1). Box stress factor r6: 0.726763171 (elem = 240, lstep = 757, node = 2). Box stress factor r7: 1.511570275E-02 (elem = 240, lstep = 604, node = 1): Maximum axial stress in cross-bracing: 21116.8328 (elem = 427, lstep = 636, node = 1). *DO LOOP ON PARAMETER= IJK FROM 1.0000 TO 30.000 BY 1.0000 Element 211 (top plate) weld factor: 0.239587619. *ENDDO INDEX= IJK Element 216 (top plate) weld stress factor: 0.274100152. Element 217 (top plate) weld stress factor: 0.227870332. Element 218 (top plate) weld stress factor: 2.449996732E-02. Element 219 (top plate) weld stress factor: 3.036511304E-02. Element 220 (top plate) weld stress factor: 2.799095942E-02. Element 221 (top plate) weld stress factor: 2.316616858E-02. Element 222 (top plate) weld stress factor: 0.2489655. Element 223 (top plate) weld stress factor: 0.185316081. Element 224 (top plate) weld stress factor: 3.124990964E-02. Element 225 (top plate) weld stress factor: 2.975341719E-02. Element 226 (top plate) weld stress factor: 2.394640208E-02. Element 227 (top plate) weld stress factor: 2.099066315E-02. Element 228 (top plate) weld stress factor: 0.218016762. Element 229 (top plate) weld stress factor: 0.219761917. Element 230 (top plate) weld stress factor: 2.571485935E-02. Element 231 (top plate) weld stress factor: 3.005403619E-02. Winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1use-sum. Maximum vertical load in any single pedestal: 53042.1672 lb. Maximum total vertical pedestal load: 111389.239 lb. Maximum shear load in any single pedestal:39960.046 lb. Maximum fuel-box contact load at any one location: 4373.12154 lb. Maximum rack/rack vertical intertie load at any one location: 54532.7439 lb. Maximum rack/rack horizontal intertie load at any one location: 3765.4236 lb. Maximum corner x-displacement, top plate: 0.188436313 in. Maximum corner x-displacement, baseplate: 4.021895771E-02 in. Maximum corner y-displacement, top plate: 0.408960187 in. Maximum corner y-displacement, baseplate: 0.19487087 in. Pedestal stress factor r1: 9.035914428E-02. Pedestal stress factor r2: 9.113358987E-02. Pedestal stress factor r3: 0.126141079. Pedestal stress factor r4: 0.126128119. Pedestal stress factor r5: 0.203041265. Pedestal stress factor r6: 0.204236961. Corresponding r1p = 8.509162978E-02, r2p = 6.845012825E-02, r3p = 0.101375541, r4p = 9.473439946E-02. Pedestal stress factor r7: 9.114295418E-02. rac1use.sum Winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1use.sum. Maximum vertical load in any single pedestal: 51223.4354 lb. Maximum total vertical pedestal load: 111482.836 lb. Maximum shear load in any single pedestal: 39977.5779 lb. Maximum fuel-box contact load at any one location: 11768.5629 lb. Maximum rack/rack vertical intertie load at any one location: 48930.2236 lb. Maximum rack/rack horizontal intertie load at any one location: 3478.63585 lb. Maximum corner x-displacement, top plate: .251264165 in. Maximum corner x-displacement, baseplate: .120556722 in. Maximum corner y-displacement, top plate: .382567581 in. Maximum corner y-displacement, baseplate: .138739342 in. Pedestal stress factor r1: 8.726087255E-02. Pedestal stress factor r2: 9.111082452E-02. Pedestal stress factor r3: .118681005. Pedestal stress factor r4: .126096612. Pedestal stress factor r5: .203130352. Pedestal stress factor r6: .204326567. Corresponding r1p = 8.512896262E-02, r2p = 5.313108243-02, r3p = .117734471, r4p = 7.353296941E-02. Pedestal stress factor r7: 8.575269414E-02. Page (illegible) c:\qemm\loadhi /r:1 c:\mouse\mscmouse /a5 path c:\;c:\dos;c:\qemm;c:\emrcnissa;c:\ansys50a\bin;c:\lantasti; c:\nu;c:\wp51 smartdrv SET SYMANTEC=C:\SYMANTEC SET NU-C:\DOS SET TEMP=C:DOS SET NISA=C:\EMRCNISA c:\qemm\loadhi /r:1 c:\dos\doskey.com c:\nu\image c: device=c:\qemm\dosdata.sys device=c:\qemm\qemm386.sys ram be:n r:3 x=ec00-ecff x=f400-f4ff device=c:\qemm\dos-up.sys @c:\qemm\dos-up.dat dos=high.umb files=100 buffers=50 lastdrive=e shell=c:\qemm\loadhi.com /r:2 c:\command.com c:\ /p C:\WINDOWS\SMARTDRV.EXE @REM THE CHECK LINE BELOW PROVIDES ADDITIONAL SAFETY FOR STACKER DRIVES. @REM PLEASE DO NOT REMOVE IT. @REM C:\STACKER\CHECK /WP PATH C:\WINDOWS;C:\QEMM;C:\;C:\DOS;C:\WP51\;C:\NU;C:\STACKER; c:\dyna2d;c:\lanta: call C:\LANTASTI\STARTNET.BAT C:\CDROMDRV\MSCDEX.EXE /V /E /D:CDROM001 /M:15 REM C:\DOS\SMARTDRV.EXE 2048 PROMPT $p$g set mouse=C:\MOUSE C:\MOUSE\mouse.EXE /Q SET TEMP=C:\DOS c:\dos\doskeyREM DEVICE=C:\DOS\SETVER.EXE DEVICE=C:\QEMM\QEMM386.SYS RAM rem device=c:\dos\himem.sys DOS=HIGH FILES=100 BUFFERS=30 SHELL=C:\DOS\COMMAND.COM C:\DOS\ /E:512 /P DEVICE=C:\CDROMDRV\GSPDX.SYS /D:CDROM001 /P:310 /N:1 LASTDRIVE=g [COMMON] DEVICE=C:\WINDOWS\MOUSE.SYS /Y DEVICE=C:\WINDOWS\SMARTDRV.EXE /DOUBLE_BUFFER STACKS=9,256 c:\windows\vqautil\clmode t6=0 t8=1 t1=2 t2=0 @REM THE CHECK LINE BELOW PROVIDES ADDITIONAL SAFETY FOR STACKER DRIVES. @REM PLEASE DO NOT REMOVE IT. (Illegible):\STACKER\CHECK /WP (Illegible):\qemm\loadhi /r:2 C:\SMARTDRV.EXE /L /X c:\qemm\loadhi /r:3 c:\qemm\files +50 c:\qemm\loadhi /r:1 c:\dos\share.exe /1:65 c:\qemm\loadhi c:\qemm\buffers +15 c:\qemm\loadhi c:qemm\lastdriv +2 PROMPT $P$G PATH C:\WINDOWS;C:\;C:\QEMM;C:\DOS;C:\LANTASTI;C:\LANTASTI.NET;C:\WP51\; path=%path%;c:\f77515\bin;c:\xtgold;c:\norton SET TEMP=C:\TEMP SET TMP=C:\ SET SYMANTEC=C:\SYMANTEC set NU=c:\NORTON rem date rem C:\NORTON\NDD c:/Q rem C:\NORTON\IMAGE c: rem C:\QEMM\LOADHI /R:3 C:\DOS\DOSKEY /INSERT call C:\LANTASTI\STARTNET.BAT c:\qemm\loadhi c:\mouse\mouse date call c:lantasti\laserwin.bat SET LFDIR=C:\F77L5\BIN SET INIT=C:\F77L5\BIN SET LFDIR=C:\F77L5\BIN SET LBCFG=C:\F77L5\BIN\LB.CFG SET LBKEY=C:\F77L5\BIN\LB.KEY SET LBMAC=C:\F77L5\BIN\LB.MAC SET LBPST=C:\F77L5\BIN\LB.PST SET LBLNG=C:\F77L5\BIN\LB.LNG DEVICE=C:\QEMM\QEMM386.SYS RAM DBF:0 DEVICE=C\PLUGPLAY\DRIVERS\DOS\DWCFGMG.SYS /NOLOCK DEVICE=c:qemm\loadhi.sys C:\QUICKCMD\CMD640X.SYS /A33 /L FILES=15 BUFFERS=15,0 STACKS=9,256 SHELL=C:\DOS\COMMAND.COM /P /E:1024 LASTDRIVE=d device=C:\STACKER\DPMS.EXE device=c:\qemm\loadhi.sys C:\STACKER\STACHIGH.SYS dos=high,umb rem DEVICE=C:\SMARTDRV.EXE /DOUBLE_BUFFER rem DEVICE=C:\IFSHLP.SYS Winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1use.sum. Maximum vertical load in any single pedestal: 52251.4672 lb. Maximum total vertical pedestal load: 111523.6 lb. Maximum shear load in any single pedestal: 38805.1286 lb. Maximum fuel-box contact load at any one location: 11614.9373 lb. Maximum rack/rack vertical intertie load at any one location: 52801.9823 lb. Maximum rack/rack horizontal intertie load at any one location: 3843.42255 lb. Maximum corner x-displacement, top plate: .206584866 in. Maximum corner x-displacement, baseplate: 9.776364748E-02 in. Maximum corner x-displacement, top plate: .400667815 in. Maximum corner x-displacement. baseplate: .14144151 in. Pedestal stress factor r1: 8.901215991E-02. Pedestal stress factor r2: 9.732297694E-02. Pedestal stress factor r3: .132169207. Pedestal stress factor r4: .134694178. Pedestal stress factor r5: .19717267. Pedestal stress factor r6: .198334145. Corresponding r1p = 8.263232835E-02, r2p = 9.732297694E-02, r3p = 1.965009542E-02, r4p = .134694178. Pedestal stress factor r7: 9.549856499E-02. Winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2). File \\jrb-66\c-drive\winco\rack\rac1use.sum. Maximum vertical load in any single pedestal: 53194.6549 lb. Maximum total vertical pedestal load: 111510.902 lb. Maximum shear load in any single pedestal: 42319.7815 lb. Maximum fuel-box contact load at any one location: 4262.67681 lb. Maximum rack/rack vertical intertie load at any one location: 56250.5999 lb. Maximum rack/rack horizontal intertie load at any one location: 4125.97222 lb. Maximum corner x-displacement, top plate: 0.237648136 in. Maximum corner x-displacement, baseplate: 8.166648269E-02 in. Maximum corner x-displacement, top plate: 0.512978527 in. Maximum corner x-displacement. baseplate: 0.274416501 in. Pedestal stress factor r1: 9.061891233E-02. Pedestal stress factor r2: 9.542331278E-02. Pedestal stress factor r3: 0.146763177. Pedestal stress factor r4: 0.132065059. Pedestal stress factor r5: 0.215032084. Pedestal stress factor r6: 0.216297638. Corresponding r1p = 9.011649244E-02, r2p = 3.53207394E-02, r3p = 0.146763177, r4p = 4.888360505E-02. Pedestal stress factor r7: 0.106043405. Form "HOLTEC INTERNATIONAL CALCULATION SHEET" omitted. ! file \\jrb-66\c-drive\winco\rack\rac1gen.inp 11/16/94 14:28,mod2/95 ! incorporates modifications discussed at client meeting of 05/26/94 ! incorporates new top and middle plate models 09/30/94 ! incorporates new cross-bracing model 11/10/94 ! incorporates changes 2/27/95 /output,rac1gen,1st /clear /filnam,rac1gen /title,winco rack analysis, 6x5x10, transient, gen pass, full /units,bin /prep7 ! input parameters bxnx=6 ! number of boxes in x-direction (bxnx <= 10) bxny=5 ! number of boxes in y-direction (bxny <= 10) bxtk=0.160 ! box thickness bxht=219.5 ! box height bxid=10.19 ! box interior dimension bxpt=14.13 ! box pitch tptk=0.25 ! top plate thickness mptk=0.25 ! middle plate thickness bptk=0.75 ! baseplate thickness taaa=2.5000 ! top plate apron cross-sectional area taiy=0.540 ! top plate apron iyy box section taiz=4.667 ! top plate apron izz taht=1.000 ! top plate apron effective height maaa=1.9375 ! middle plate apron cross-sectional area maiy=2.989 ! middle plate apron iyy waiz=2.989 ! middle plate apron izz maht=5.815 ! middle plate apron effective height cbwd=3.0 ! cross-bracing width cbtk=0.5 ! cross-bracing thickness pdht=6.625 ! pedestal height (to center of baseplate) ssem=28.14e6 ! stainless steel elastic modulus @ temp ssnu=0.3 ! stainless steel poisson's ratio ssmd=7.51e-4 ! stainless steel mass density ssmd=ssmd*1.0862 ! ssmd * unmodeled mass adjustment factor ssdr=0.0018 ! stainless steel damping ratio fawt=800.0 ! fuel assembly weight ofset=.5*bxid+bxtk+1.08 ! intertie offset abs value ofset1=.5*bxpt ! computed parameters gg=386.4 ! acceleration due to gravity bxht2=bxht/2 ! 1/2 of box height bxht4=bxht/4 ! 1/4 of box height bxht8=bxht/8 ! 1/8 of box height bxbt10=bxht/10 ! 1/10 of box height bxod=bxid+2*bxtk ! box outer dimension bxaa=bxod**2-bxid**2 ! box cross-sectional area bxiy=(bxod**4-bxid**4)/12 ! box iyy bxiz=(bxod**4-bxid**4)/12 ! box izz cbaa=cbwd*cbtk ! cross-bracing cross-sectional area cbiy=cbtk*cbwd**3/12 ! cross-bracing iyy cbiz=cbwd*cbtk**3/12 ! cross-bracing izz fgap=bxpt-bxod ! flux gap ssee=ssem*(bxpt/fgap)**3 ! plate effective elastic modulus geff=ssem*bxpt/fgap*(1/2.6) ! plate effective in plane shear n101x=-(bxnx-1)*bxpt)/2 ! node location references n101y=-(bxny-1)*bxpt)/2 n101z=pdht nn1=100+bxnx ! node number references nn2=101+10*(bxny-1) nn3=nn1+10*(bxny-1) parsav,scalar ! element types et,1,beam4 ! elastic beam et,2,shell63 ! elastic shell et,3,mass21 ! point mass et,10,beam4 ! offset extension ! real constants r,1,bxaa,bxiz,bxiy,bxod,bxod ! boxes r,2,tptk ! top plate r,3,mptk ! middle plate r,4,bptk ! baseplate r,5,0,0,3*fawt/gg,0,0,0 ! vertical fhu mass r,6,taaa,taiz,taiy,taht,taht ! top plate apron r,7,maaa,maiz,maiy,maht,maht ! middle plate apron r,8,cbaa,cbiz,cbiy,cbtk,cbwd ! cross-bracing r,15,bxaa,bxiz,b)dy,bxod,bxod ! off-set elements ! material properties mp,ex,1,ssem ! 240-304L stainless steel mp,nuxy,1,ssnu mp,dens,1,ssmd mp,damp,1,ssdr mp,ex,2,ssee ! plate effective material mp,nuxy,2,ssnu mp,dens,2,ssmd mp,damp,2,ssdr mp,gxy,2,geff mp,ex,4,10*ssem ! offset properties mp,nuxy,4,.3 mp,dens,4,ssmd ! node definitions n,101,n101x,n101y,n101z ngen,bxnx,1,101,101,1,bxpt,0,0 ngen,bxny,10,101,nn1,1,0,bxpt,0 ngen,3,100,101,nn3,1,0,0,bxht2 nsel,s,loc,z,nz(101) ngen,2,2000,all,,,0,0,0 ngen,9,100,2101,2199, 1,0,0,bxht8 nsel,s,node,,2101,2999,1 nsel,r,loc,y,ny(101) nsel,r,loc,x,nx(101)+1,nx(nn1)-1 nsel,r,loc,z,nz(101)+1,nz(301)+1 ndele,all nsel,s,node,,2101,2999,1 nsel,r,loc,y,ny(nn2) nsel,r,loc,x,nx(101)+1,nx(nn1)-1 nsel,r,loc,z,nz(101)+1,nz(301)+1 ndele,all nsel,all ngen,11,100,2102,nn1+1999, 1,0,0,bxht10 ngen,11,100,nn2+2001,nn3+1999, 1,0,0,bxht10 nsel,all ! added offset nodes start with 901 ngen,2,600,301,301,1,-ofset,+ofset1,0 ngen,2,600,nn1+200,nn1+200,1,ofset,ofset1,0 ngen,2,600,nn2+200,nn2+200,1,-ofset,-ofset1,0 ngen,2,600,nn3+200,nn3+200,1,ofset,-ofset1,0 ngen,2,601,301,301,1,ofset1,-ofset0 ngen,2,601,nn1+200,nn1+200,1,-ofset1,-ofset0 ngen,2,601,nn2+200,nn2+200,1,ofset1,ofset0 ngen,2,601,nn3+200,nn3+200,1,-ofset1,ofset,0 ! element definitions type,1 ! boxes real,1 mat,1 e,2101,2201 egen,bxnx,1,all egen,bxny,10,all egen,8,100,all nsel,s,loc,z,nz(2902) esln,s,0 egen,3,100,all nsel,all type,2 ! top plate real,2 mat,2 e,301,302,312,311 esel,s,real,,2 egen,bxnx-1,1,all esel,s,real,,2 egen,bxny-1,10,all type,2 ! middle plate real,3 mat,2 e,201,,202,212,211 esel,s,real,,3 egen,bxnx-1,1,all esel,s,real,,3 egen,bxny-1,10,all type,2 ! baseplate real,4 mat,1 e,101,102,112,111 esel,s,real,,4 egen,bxnx-1,1,all esel,s,real,,4 egen,bxny-1,10,all type,3 ! vertical fhu masses real,5 e,101 *get,mxel,elem,,num,max egen,bxnx,1,mxel egen,bxny,10,mxel,mxel+bxnx-1 type, ! top plate apron real,6 mat,1 e,301,302 esel,s,real,6 egen,bxnx-1,1,all esel,s,real,,6 egen,2,10*(bxny-1),all e,301,311 *get,mxel,elem,num,max egen,bxny-1,10,mxel egen,2,bxnx-1,mxel,mxel+bxny-2 type,1 ! middle plate apron real,7 mat,1 e,201,202 esel,s,real,7 egen,bxnx-1,1,all esel,s,real,7 egen,2,10*(bxny-1),all e-201,211 *get,mxel,elem,num,max egen,bxny-1,10,mxel egen,2,bxnx-1,mxel,mxel+bxny-2 type,1 ! cross-bracing real,8 mat,1 e,2101,2202,2102 *get,mxel,elem,num,max egen,bxnx-2,101,mxel e,nn1-1+2400,nn1+2400,nn1+2300 e,nn1+2000,nn1-1+2100,nn1-1+2000 *get,mxel,elem,num,max egen,bxnx-2,99,mxel e,2502,2501,2401 e,2501,2702,2602 e,2702,2803,2703 get,mxel,elem,num,max egen,bxnx-3,101,mxel e,nn1-1+2900,nn1+2800,nn1+2700 e,nn1+2400,nn1-1+2600,nn1-1+2500 e,nn1-1+2600,nn1-2+2700,nn1-2+2600 *get,mxel,elem,num,max egen,bxnx-3,99,mxel e,3002,2901,2801 esel,s,reg,8 egen,2,(bxny-1)*10,all esel,none e,2101,2211,2111 *get,mxel,elem,,num,max egen,bxny-1,110,mxel e,nn2+2000,nn2-10+2100,nn2-10+2000 *get,mxel,elem,num,max egen,bxny-1,90,mxel egen,2,400,all egen,2,bxnx-1,all esel,all ! added offset elements 2/95 type,10 mad, 15 mat,4 e,301,901 e,311,901 e,321,901 e,301,902 e,302,902 e,303,902 e,nn1+200,nn1+800 e,nn1+210,nn1+800 e,nn1+220,nn1+800 e,nn1+200,nn1+801 e,nn1+199,nn1+801 e,nn1+198,nn1+801 e,nn2+200,nn2+800 e,nn2+190,nn2+800 e,nn2+180,nn2+800 e,nn2+200,nn2+801 e,nn2+201,nn2+801 e,nn2+202,nn2+801 e,nn3+200,nn3+800 e,nn3+190,nn3+800 e,nn3+180,nn3+800 e,nn3+200,nn3+801 e,nn3+199,nn3+801 e,nn3+198,nn3+801 nummrg,node /output,rac1gen,afl aflist etlist,all rlist,all mplist,all nlist,all elist,all /output,rac1gen,1st,,append finish /solution antype,substr seopt,rac1gen,3,0,0 nsel,s,node,,101,301,100 ! define mdofs nsel,s,node,,nn1,nn1+200,100 nsel,s,node,,nn2,nn2+200,100 nsel,s,node,,nn3,nn3+200,100 nsel,s,node,,901,999,1 ! add 12 more offset nodes as dof m,all,ux,,,uy,uz nsel,all m,node(0,0,pdht),uz outpr,basic,last save ! (required for expansion pass) solve finish /show,rac1gen,grp ! produce plots /vup,,z /view,,-1,-1,1 /color,elem,whit esel,s,type,,2 /color,elem,blac esel,all /pnum,node,1 esel,u,real,,5 esel,u,real,,6 esel,u,real,,7 esel,u,real,,8 /title,rac1gen nodes eplot ! plot 1 esel,s,real,1 /pnum,node,0 /pnum,elem,1 /title,rac1gen beam elements - ports eplot ! plot 2 esel,s,real,,2 esel,s,real,,3 esel,s,real,,4 /title,rac1gen plate elements eplot ! plot 3 esel,s,real,,5 /title,rac1gen mass elements eplot ! plot 4 esel,s,real,,6 esel,a,real,,7 /title,rac1gen beam elements - aprons eplot ! plot 5 esel,s,real,,8 /title,rac1gen beam elements - cross-bracing eplot ! plot 6 /show,term /title,winco rack analysis, 6x5x10,transient, gen pass, full esel,all /output,term ! file \\jrb-66\c-drive\winco\rack\rac1use.inp 11/21/94 13-29 ! incorporates modifications discussed at client meeting of 05/26/94 ! incorporates use pass postprocessors 09/12/94 ! incorporates revised rack/rack intertie model 11/21/94 ! incorporates 2/95 modifications /output,rac1use,1st /clear /menu,no /filnam,rac1use /title,winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2) /units,bin /config,nres,2002 /config,nbuf,16 /config,szbio,32768 /config,ncont,100 /config,nvpag,48 /prep7 parres,new,rac1gen,par ! additional input parameters pdod=9.0 ! pedestal outer diameter pdid=5.5 ! pedestal inner diameter pdag=50.85 ! pedestal effective cross-sectional area (see calc) pdiy=626.4 ! pedestal effective iyy (see calc) pdiz=626.4 ! pedestal effective izz (see calc) pds1=2.5e6 ! pedestal normal stiffness (see calc) pds2=2.5e7 ! pedestal sticking stiffness faod=7.00 ! fhu outside dimension plgf=0.8 ! pedestal/liner gap friction coefficient wtmd=9.06e-5 ! water mass density rwgp=5.5 ! rack/wall gaps instv=7.00e6 ! rack/rack intertie stiffness, vertical mod 2/95 insth=4.03e6 ! rack/rack intertie stiffness, horizontal nls=1995 ! number of transient load steps to be run ! additional computed parameters pi=acos(-1) ! pi pdor=pdod/2 ! pedestal outer radius pdir=pdid/2 ! pedestal inner radius pdaa=pi*(pdor**2-pdir**2) ! pedestal cross-sectional area (female portion) igap=(bxid-faod)/2 ! box/fhu initial gap bxtt=bxnx*bxny ! total number of boxes famt=5/6*fawt*bxtt/4/gg ! lumped fhu masses at top and bottom of boxes famm=4/3*fawt*bxtt/4/gg ! lumped fhu masses at middle of boxes ! element types et,3,mass21 ! structural mass et,4,contac52 ! contact keyopt,4,7,1 et,5,fluid38 ! dynamic fluid coupling keyopt,5,3,2 keyopt,5,6,3 et,6,matrix50 ! superelement et,7,combin40 ! combination, x-displacement keyopt,7,3,1 et,8,combin40 ! combination, y-displacement keyopt,8,3,2 et,9,combin40 ! combination, z-displacement keyopt9,3,3 et,11,combin14 ! spring-damper, x-displacement keyopt,11,2,1 et,12,combin14 ! spring-damper, y-displacement keyopt, 12,2,2 et,13,combin14 ! spring-damper, z-displacement keyopt, 13,2,3 ! real constants r,6,famt,famt,0,0,0,0 ! horizontal fhu masses - top and bottom r,7,famm,famm,0,0,0,0 ! horizontal fhu masses - middle r,8,pds1,0,1,pds2 ! pedestal/liner gaps cc=(faod+bxid)/2 ! box/fhu fluid couplings - top and bottom hh=(bxid-faod)/2 mmh=2*wtmd*cc**3*bxht4/3/hh mm1=wtmd*faod**2*bxht4 mm2=wtmd*bxid**2*bxht4 mmhc=mmh*bxtt/4 mm1c=mm1*bxtt/4 r,9,mm2c,mm1c,mmhc,mmhc rr0=0.1 ! box/fhu gap stiffness (see roark, 6 ed, table 24, case 10a) aa=bxid/2 rat1=-rr0/aa ll17t1=(1-ssnu)*(1-rat1**4)/4 ll17t2=rat1**2*(1+(1+ssnu)*log(1/rat1)) ll17=(1-ll17t1-ll17t2)/4 ll11t1=4*rat1**2*(2+rat1**2)log(1/rat1) ll11=(1+4*rat1**2-5*rat1**4-ll11t1)/64 dd=ssem*bxtk**3/12/(1-ssnu**2) kk0=2*pi*dd/aa**2/(ll 17/(1+ssnu)-2*ll11) ! multiply by number of boxes involved rrss=bxnx*bxny/4 kk1=rrss**kk0 r,10,0,0.0045*kk1,0,igap,0,kk1 mmh=2*wtmd*cc**3bxht2/3/hh ! box/fhu fluid couplings - middle mm1=wtmd*faod**2*bxht2 mm2=wtmd*bxid**2*bxht2 mmhc=mmh*bxtt/4 mm1c=mm1*bxtt/4 mm2c=mm2*bxtt/4 r,11,mm2c,mm1c,mmhc,mmhc bb=bxid+2*bxtk ! rack/ground fluid couplings - top and bottom hh=(bxpt-bxod)/2 mm1=bxtt*wtmd*bb**2*bxht4 mm2=mm1*(1+2*hh/bb) mmh=mm1*bb*(1+hh/bb)**2/2/hh mm1c=mm1/4 mm2c=mm2/4 mmhc=mmh/4 r,12,mm2c,mm1c,mmhc,mmhc mm1=bxtt*wtmd*bb**2*bxht2 ! rack/ground fluid couplings - middle mm2=mm1(1+2*hh/bb) mmh=mm1*bb*(1+hh/bb)**2/2/hh mm1c=mm1/4 mm2c=mm2/4 mmhc=mmh/4 r,13,mm2c,mm1c,mmhc,mmhc r,14,0,0-0045*pds1,0,-0.0001,0,1000 ! pedestal/liner dampers r,16,instv ! rack/rack connections, vertical r,17,insth ! rack/rack connections, horizontal r,18,0,0.0045*1e6,0,rwgp,0,1000 ! rack/wall gaps ! material properties mp,mu,3,plgf ! pedestal/liner gaps ! superelement definition type,6 se,rac1gen selist,rac1gen ! node definitions nsel,u,node,,node(0,0,pdht) ngen,2,1000,all,,,0,0,0 n,1,nx(101),ny(101),0 n,nn1-100,nx(nn1),ny(nn1),0 n,nn2-100,nx(nn2),ny(nn2),0 n,nn3-100,nx(nn3),ny(nn3),0 nsel,s,loc,z,nz(301) nsel,r,node,,301,399 ngen,2,100,all,,,0,0,0 nsel,all ! element definitions type,3 ! fhu masses - top and bottom (2-9) real,6 e,1101 e,nn1+1000 e,nn2+1000 e,nn3+1000 esel.s,real,,6 egen,2,200,all type,3 ! fhu masses - middle (10-13) real,7 e,1201 e,nn1+1100 e,nn2+1100 e,nn3+1100 type,4 ! pedestal/liner gaps (14-17) real,8 mat,3 e,1,101 e,nn1-100,nn1 e,nn2-100,nn2 e,nn3-100,nn3 type,9 ! pedestal/liner dampers (18-21) real,14 e,1,101 e,nn1-100,nn1 e,nn2-100,nn2 e,nn3-100,nn3 type,7 ! box/fhu gaps, x-displacement (22-29) real,10 e,101,1101 e,1101,101 e,nn1,nn1+1000 e,nn1+1000,nn1 e,nn2,nn2+1000 e,nn2+1000,nn2 e,nn3,nn3+1000 e,nn3+1000,nn3 type,8 ! box/fhu gaps, y-displacement (30-37) real,10 e,101,1101 e,1101,101 e,nn1,nn1+1000 e,nn1+1000,nn1 e,nn2,nn2+1000 e,nn2+1000,nn2 e,nn3,nn3+1000 e,nn3+1000,nn3 esel,s,type,,7 esel,a,type,,8 egen,3,100,all ! (38-69) type,5 ! box/fhu fluid couplings - top and bottom (70-73) real,9 e,1101,101 e,nn1+1000,nn1 e,nn2+1000,nn2 e,nn3+1000,nn3 esel,s,real,,9 egen,2,200,all ! (74-77) type,5 ! box/fhu fluid couplings - middle (78-81) real, 11 e,1201,201 e,nn1+1100,nn1+100 e,nn2+1100,nn2+100 e,nn3+1100,nn3+100 type,5 ! rack/ground fluid couplings - top and bottom (82-89) real,12 e,101,1 e,nn1,nn1-100 e,nn2,nn2-100 e,nn3,nn3-100 e,301,1 e,nn1+200,nn1-100 e,nn2+200,nn2-100 e,nn3+200,nn3-100 type,5 ! rack/ground fluid couplings - middle (90-93) real,13 e,201,1 e,nn1+100,nn1-100 e,nn2+100,nn2-100 e,nn3+100,nn3-100 type,13 ! rack/rack interties, vertical (94-97) real 16 e,901,nn1+800 e,nn2+800,nn3+800 e,902,nn2+801 e,nn1+801,nn3+801 type, 11 ! rack/rack interties, horizontal (98-101) real 17 e,902,nn2+801 e,nn1+801,nn3+801 type,12 real,17 e,901,nn1+800 e,nn2+800,nn3+800 type,7 ! rack/wall gaps, x-displacement (102-105) real,18 e,401,301 e,nn2+300,nn2+200 e,nn1+200,nn1+300 e,nn3+200,nn3+300 type,8 ! rack/wall gaps, y-displacement (106-109) real,18 e,401,301 e,nn1+300,nn1+200 e,nn2+200,nn2+300 e,nn3+200,nn3+300 esel,all /output,racluse,afl aflist etlist,all rlist,all mplist,all nlist,all elist,all cplist,all /output,racluse,1st,,append finish /solution antype,trans nsel,s,loc,z,0 nsel,a,node,,401,499 d,all,all,0 nsel,all acel,0,0,386.4*0.879 ! gravity factored by buoyancy /output,racluse,afl,,append dlist,all /output,racluse,1st,,append trnopt,full ! full transient analysis nropt,auto ! program-chosen newton-raphson option ncnv,0,50 ! equilibrium iteration convergence controls timint, off ! static load step kbc, 1 ! stepped load application nsubst,2 time,1.0e-6 outpr,basic,last outres,all solve /output,term *msg,info > static equilibrium solution complete /output,racluse,1st,append timint,on ! dynamic load steps deltim,0.0025,0.00001,0.005,on ! specify size of substeps cnvtol,f,,0.005,,1 neqit,100 autots,on *dim,inparr,,nls*2 ! process accel data *dim,xyzarr,,nls,4 *vlen,nls*2 *vread,inparr(1),p2ew3,ans ! read and process x-accel data (1x,f8.4,1x,f12.5) *do,ii,1,nls,1 xyzarr(ii,1)=inparr(ii*2-1) xyzarr(ii,2)=inparr(ii*2) *enddo vread,inparr(1),p2ns3,ans ! read and process y-accel data (1xf8.4,1x,f12.5) *do,ii,1,nls,1 *if,inparr(ii*2-1),ne,xyzarr(ii,1),then /output,term *msg,warn,ii *****DATA ERROR ON TIME NUNBER %i Y-ACC FILE***** /output,racluse,1st,append *endif xyzarr(ii,3)=inparr(ii*2) *enddo *vread,inparr(1),p2vt3,ans ! read and process z-accel data (1x,f8.4,1x,f12.5) *do,ii,1,nls,1 *if,inparr(ii*2-1),ne,xyzarr(ii,1),then /output,term *msg,warn,ii *****DATA ERROR ON TIME NUMBER %i, Z-ACC FILE***** /output,racluse,1st,,append *endif xyzarr(ii,4)=inparr(ii*2) *enddo /output,term msg,info >seismic data preprocessing complete /output,tenn outpr,all,none outres,nsol,last outres,rsol,last outres,misc,last *do,ii,1,nls,1 ! begin solution loop kbc,0 ! ramped load application time,xyzarr(ii,1) acel,xyzarr(ii,2),xyzarr(ii,3),(xyzarr(ii,4)+386.4)*0.879 solve *enddo /output,racluse,xyz *status,xyzarr,1,nls,1,4 /output,racluse,1st,,append finish save /sys,del racluse.esa /sys,del racluse.ema /sys,del racluse.osa /sys,del racluse.tri ! additional input parameters nlp=1995 ! number of load steps to process sy=25000 ! 3041 min yield strength @ 100 deg f su=70000 ! 3041 min tensile strength @ 100 deg f pdkk=2.0 ! pedestal length coefficient pdll=7.25 ! pedestal unsupported length (max) ! additional computed parameters pdry=sqrt(pdiy/pdaa) ! pedestal radius of gyration yy pdrz=sqrt(pdiz/pdaa) ! pedestal radius of gyration zz ft=0.6*sy ! allowable stress in tension fv=0.4*sy ! allowable stress in shear ! cc=sqrt(2*pi**2*ssem/sy) ! pedestal allowable stress in compression klrpdy=pdkk**pdll/pdry klrpdz=pdkk*pdll/pdrz ! fapdn=(1.0-klrpdy**2/2/cc**2)*sy ! fapd1=5/3+3*klrpdy/8/cc ! fapd2=klrpdy**3/8/cc**3 ! fapd=fapdn/(fapd1-fapd2) ! kl/r of pedestal < 120 fapd=sy*(.47-klrpdy/444) fb=0.6*sy ! allowable stress in bending cmy=0.85 ! pedestal allowable combined bending and compression cmz--0.85 ! feypd=12.0*pi**2*ssem/23-0/klrpdy**2 feypd=pi**2**ssem/(2.15*(klrpdy)**2) fezpd=feypd /post1 *dim,lststp,array,nlp *do,lpv,1,nlp,1 set,lpv *get,dummy,active,,set,sbst lststp(lpv)=dummy *enddo /output,term *msg,info >last substep number processing complete /output,racluse,1st,,append finish parsav,all /post26 /xrange,0,20 /color,curve,whit /show,racluse,,grp /output,racluse,sum *msg,info winco rack analysis, 6x5x10, transient, use pass, full, mu=0.8 (th2) *msg,info file \\jrb-66\c-drive\winco\rack\racluse.sum /output,term *msg,info >file header information complete /output,racluse,1st,,append ! solution summary data ! solu,2,dtime,,tss ! solu,3,eqit,,nei ! solu,4,cnvg,,cvi ! solu,5,resfrq,,rfr ! plvar,2 ! plot 1 ! plvar,3 ! plot 2 ! plvar,4 ! plot 3 ! plvar,5 ! plot 4 /output,term *msg,info >solution summary data complete /output,racluse,1st,,append ! maximum vertical load in any single pedestal esol,2,14,,smisc,1,pfx1 esol,3,15,,smisc,1,pfx2 esol,4,16,,smisc,1,pfx3 esol,5,17,smisc,1,pfx4 plvar,2 ! plot 5 plvar,3 ! plot 6 plvar,4 ! plot 7 plvar,5 ! plot 8 *get,pfx1mn,vari,2,extrem,vmin *get,pfx2mn,vari,3,extrem,vmin *get,pfx3mn,vari,4,extrem,vmin *get,pfx4mn,vari,5,extrem,vmin *use, mnof4,pfx1mn,pfx2mn,pfx3mn,pfx4mn minn=abs(minn) /output,racluse,sum,,append *msg,info,minn Maximum vertical load in any single pedestal: %g lb /output,term *msg,info >maximum vertical load in any single pedestal complete /output,racluse,1st,,append ! maximum total vertical pedestal load add,6,2,3,,temp add,2,4,5,6,pfxt plvar,2 ! plot 9 *get,pfxtmn,vari,2,extrem,vmin pfxtmn=abs(pfxtmn) /output,racluse,sum,,append *msg,info,pfxtmn Maximum total vertical pedestal load: %g lb /output,term *msg,info > maximum total vertical pedestal load complete /output,racluse,1st,,append ! maximum shear load in any single pedestal esol,2,14,,smisc,2,pfy1 esol,3,15,,smisc,2,pfy2 esol,4,16,,smisc,2,pfy3 esol,5,17,,smisc,2,pfy4 plvar,2 ! plot 10 plvar,3 ! plot 11 plvar,4 ! plot 12 plvar,5 ! plot 13 *get,pfy1mx,vari,2,extrem,vmax *get,pfy2mx,vari,3,extrem,vmax *get,pfy3mx,vari,4,extrem,vmax *get,pfy4mx,vari,5,extrem,vmax *use,mxof4,pfy1mx,pfy2mx,pfy3mx,pfy4mx /output,racluse,sum,,append *msg,info,maxx Maximum shear load in any single pedestal: %g lb /output,term *msg,info >maximum shear load in any single pedestal complete /output,racluse,1st,append ! maximum fuel-box contact load at any one location, first eight esol,2,22,,smisc,2,f22 esol,3,23,,smisc,2,f23 esol,4,24,,smisc,2,f24 esol,5,25,,smisc,2,f25 esol,6,26,,smisc,2,f26 esol,7,27,,smisc,2,f27 esol,8,28,,smisc,2,f28 esol,9,29,,smisc,2,f29 *get,f22mn,vari,2,extrem,vmin *get,f23mn,vari,3,extrem,vmin *get,f24mn,vari,4,extrem,vmin *get,f25mn,vari,5,extrem,vmin *get,f26mn,vari,6,extrem,vmin *get,f27mn,vari,7,extrem,vmin *get,f28mn,vari,8,extrem,vmin *get,f29mn,vari,9,extrem,vmin *use,mnof8,f22mn,f23mn,f24mn,f25mn,f26mn,f27mn,f28mn,f29mn max1=abs(minn) /output, term *msg,info > maximum fuel-box contact load at any one location, first eight complete /output,racluse,1st,,append ! maximum fuel-box contact load at any one location, second eight esol,2,30,,smisc,2,f30 esol,3,31,,smisc,2,f31 esol,4,32,,smisc,2,f32 esol,5,33,,smisc,2,f33 esol,6,34,,smisc,2,f34 esol,7,35,,smisc,2,f35 esol,8,36,,smisc,2,f36 esol,9,37,,smisc,2,f37 *get,f30mn,vari,2,extrem,vmin *get,f31mn,vari,3,extrem,vmin *get,f32mn,vari,4,extrem,vmin *get,f33mn,vari,5,extrem,vmin *get,f34mn,vari,6,extrem,vmin *get,f35mn,vari,7,extrem,vmin *get,f36mn,vari,8,extrem,vmin *get,f37mn,vari,9,extrem,vmin *use,mnof8,f30mn,f31mn,f32mn,f33mn,f34mn,f35mn,f36mn,f37mn max2=abs(minn) /output,term *msg,info > maximum fuel-box contact load at any one location, second eight complete /output,racluse,1st,,append ! maximum fuel-box contact load at any one location, third eight esol,2,38,,smisc,2,f38 esol,3,39,,smisc,2,f39 esol,4,40,,smisc,2,f40 esol,5,41,,smisc,2,f41 esol,6,42,,smisc,2,f42 esol,7,43,,smisc,2,f43 esol,8,44,,smisc,2,f44 esol,9,45,,smisc,2,f45 plvar,6 ! plot 14 plvar,7 ! plot 15 plvar,8 ! plot 16 plvar,9 ! plot 17 *get,f38mn,vari,2,extrem,vmin *get,f39mn,vari,3,extrem,vmin *get,f40mn,vari,4,extrem,vmin *get,f41mn,vari,5,extrem,vmin *get,f42mn,vari,6,extrem,vmin *get,f43mn,vari,7,extrem,vmin *get,f44mn,vari,8,extrem,vmin *get,f45mn,vari,9,extrem,vmin use,mnof8,f38mn,f39mn,f4Omn,f41mn,f42mn,f43mn,f44mn,f45mn max3=abs(minn) /output,term *msg,info > maximum fuel-box contact load at any one location, third eight complete /output,racluse,1st,,append ! maximum fuel-box contact load at any one location, fourth eight esol,2,46,,smisc,2,f46 esol,3,47,,smisc,2,f47 esol,4,48,,smisc,2,f48 esol,5,49,,smisc,2,f49 esol,6,50,,smisc,2,f50 esol,7,51,,smisc,2,f51 esol,8,52,,smisc,2,f52 esol,9,53,,smisc,2,f53 *get,f46mn,vari,2,extrem,vmin *get,f47mn,vari,3,extrem,vmin *get,f48mn,vari,4,extrem,vmin *get,f49mn,vari,5,extrem,vmin *get,f50mn,vari,6,extrem,vmin *get,f51mn,vari,7,extrem,vmin *get,f52mn,vari,8,extrem,vmin *get,f53mn,vari,9,extrem,vmin *use,mnof8,f46mn,f47mn,f48mn,f49mn,f50mn,f51mn,f52mn,f53mn max4=abs(minn) /output,term *msg,info > maximum fuel-box contact load at any one location, fourth eight complete /output,racluse,1st,,append ! maximum fuel-box contact load at any one location, fifth eight esol,2,54,,smisc,2,f54 esol,3,55,,smisc,2,f55 esol,4,56,,smisc,2,f56 esol,5,57,,smisc,2,f57 esol,6,58,,smisc,2,f58 esol,7,59,,smisc,2,f59 esol,8,60,,smisc,2,f60 esol,9,61,,smisc,2,f61 *get,f54mn,vari,2,extrem,vmin *get,f55mn,vari,3,extrem,vmin *get,f56mn,vari,4,extrem,vmin *get,f57mn,vari,5,extrem,vmin *get,f58mn,vari,6,extrem,vmin *get,f59mn,vari,7,extrem,vmin *get,f60mn,vari,8,extrem,vmin *get,f61mn,vari,9,extrem,vmin *use,mnof8,f54mn,f55mn,f56mn,f57mn,f58mn,f59mn,f60mn,f61mn max5=abs(minn) /output,term *msg,info > maximum fuel-box contact load at any one location, fifth eight complete /output,racluse,1st,append ! maximum fuel-box contact load at any one location, sixth eight esol,2,62,,smisc,2,f62 esol,3,63,,smisc,2,f63 esol,4,64,,smisc,2,f64 esol,5,65,,smisc,2,f65 esol,6,66,,smisc,2,f66 esol,7,67,,smisc,2,f67 esol,8,68,,smisc,2,f68 esol,9,69,,smisc,2,f69 *get,f62mn,vari,2,extrem,vmin *get,f63mn,vari,3,extrem,vmin *get,f64mn,vari,4,extrem,vmin *get,f65mn,vari,5,extrem,vmin *get,f66mn,vari,6,extrem,vmin *get,f67mn,vari,7,extrem,vmin *get,f68mn,vari,8,extrem,vmin *get,f69mn,vari,9,extrem,vmin *use,mnof8,f62mn,f63mn,f64mn,f65mn,f66mn,f67mn,f68mn,f69mn max6=abs(minn) /output,term *msg,info > maximum fuel-box contact load at any one location, sixth eight complete /output,racluse,1st,,append maxx=max1 *if,max2,gt,maxx,then maxx=max2 *endif *if,max3,gt,maxx,then maxx=max3 *endif *if,max4,gt,maxx,then maxx=max4 *endif *if,max5,gt,maxx,then maxx=max5 *endif *if,max6,gt,maxx,then maxx=max6 *endif maxx=maxx*4/bxnx/bxny ! max at any location in one box /output,racluse,sum,,append *msg,info,maxx Maximum fuel-box contact load at any one location: %g lb /output,term *msg,info > maximum fuel/box contact load at any one location complete /output,racluse,1st,,append ! maximum rack/rack vertical intertie load at any one location esol,2,94,,smisc,1,f94 esol,3,95,,smisc,1,f95 esol,4,96,,smisc,1,f96 esol,5,97,,smisc,1,f97 plvar,2 ! plot 22 plvar,3 ! plot 23 plvar,4 ! plot 24 plvar,5 ! plot 25 *get,f94mx,vari,2,extrem,vmax *get,f95mx,vari,3,extrem,vmax *get,f96mx,vari,4,extrem,vmax *get,f97mx,vari,5,extrem,vmax *use,mxof4,f94mx,f95mx,f96mx,f97mx *get,f94mn,vari,2,extrem,vmin *get,f95mn,vari,3,extrem,vmin *get,f96mn,vari,4,extrem,vmin *get,f97mn,vari,5,extrem,vmin *use,mnof4,f94mn,f95mn,f96mn,f97mn minn=abs(minn) *if,minn,gt,maxx,then maxx=minn *endif intmax=maxx /output,racluse,sum,,append *msg,info,maxx Maximum rack/rack vertical intertie load at any one location: %g lb /output,term *msg,info > maximum rack/rack intertie load at any one location complete /output,racluse,1st,,append ! maximum rack/rack horizontal intertie load at any one location esol,2,98,,smisc,1,f98 esol,3,99,,smisc,1,f99 esol,4,100,,smisc,1,f100 esol,5,101,,smisc,1,f101 plvar,2 ! plot 26 plvar,3 ! plot 27 plvar,4 ! plot 28 plvar,5 ! plot 29 *get,f98mx,vari,2,extrem,vmax *get,f99mx,vari,3,extrem,vmax *get,f100mx,vari,4,extrem,vmax *get,f101mx,vari,5,extrem,vmax *use,mxof4,f98mx,f99mx,f100mx,f101mx *get,f98mn,vari,2,extrem,vmin *get,f99mn,vari,3,extrem,vmin *get,f100mn,vari,4,extrem,vmin *get,f101mn,vari,5,extrem,vmin *use,mnof4,f98mn,f99mn,f100mn,f101mn minn=abs(minn) *if,minn,gt,maxx,then maxx--minn *endif intmax=maxx /output,racluse,sum,,append *msg,info,maxx Maximum rack/rack horizontal intertie load at any one location: %g lb /output,term *msg,info > maximum rack/rack intertie load at any one location complete /output,racluse,1st,,append ! maximum corner x-displacements nsol,2,101+200,u,x,n0tx nsol,3,nn1+200,u,x,n1tx nsol,4,nn2+200,u,x,n2tx nsol,5,nn3+200,u,x,n3tx nsol,6,101,u,x,n0bx nsol,7,nn1,u,x,n1bx nsol,8,nn2,u,x,n2bx nsol,9,nn3,u,x,n3bx plvar,2,6 ! plot 30 plvar,3,7 ! plot 31 plvar,4,8 ! plot 32 plvar,5,9 ! plot 33 *get,n0txmx,vari,2,extrem,vmax *get,n1txmx,vari,3,extrem,vmax *get,n2txmx,vari,4,extrem,vmax *get,n3bxmx,vari,5,extrem,vmax *get,n0bxmx,vari,6,extrem,vmax *get,n1bxmx,vari,7,extrem,vmax *get,n2bxmx,vari,8,extrem,vmax *get,n3bxmx,vari,9,extrem,vmax *get,n0txmn,vari,2,extrem,vmin *get,n1txmn,vari,3,extrem,vmin *get,n2txmn,vari,4,extrem,vmin *get,n3txmn,vari,5,extrem,vmin *get,n0bxmn,vari,6,extrem,vmin *get,n1bxmn,vari,7,extrem,vmin *get,n2bxmn,vari,8,extrem,vmin *get,n3bxmn,vari,9,extrem,vmin *use,mxof4,n0txmx,n1txmx,n2txmx,n3txmx *use,mxof4,n0txmx,n1txmx,n2txmx,n3txmx minn=abs(minn) *if,minn,gt,maxx,then maxx=minn *endif /output,racluse,sum,,append *msg,info,maxx Maximum corner x-displacement, top plate: %g in /output,term *msg,info > maximum comer x-displacement, top plate complete /output,racluse,1st,,append *use,mxof4,n0bxmx,n1bxmx,n2bxmx,n3bxmx *use,mnof4,n0bxmn,n1bxmn,n2bxmn,n3bxmn minn=abs(minn) *if,minn,gt,maxx,then maxx=minn *endif /output,racluse,sum,,append *msg,info,maxx Maximum corner x-displacement, baseplate: %g in /output,term *msg,info > maximum comer x-displacement, baseplate complete /output,racluse,1st,,append ! maximum corner y-displacements nsol,2,101+200,u,y,n0ty nsol,3,nn1+200,u,y,n1ty nsol,4,nn2+200,u,y,n2ty nsol,5,nn3+200,u,y,n3ty nsol,6,101,u,y,n0by nsol,7,nn1,u,y,n1by nsol,8,nn2,u,y,n2by nsol,9,nn3,u,y,n3by plvar,2,6 ! plot 34 plvar,3,7 ! plot 35 plvar,4,8 ! plot 36 plvar,5,9 ! plot 37 *get,n0tymx,vari,2,extrem,vmax *get,n1tymx,vari,3,extrem,vmax *get,n2tymx,vari,4,extrem,vmax *get,n3tymx,vari,5,extrem,vmax *get,n0bymx,vari,6,extrem,vmax *get,n1bymx,vari,7,extrem,vmax *get,n2bymx,vari,8,extrem,vmax *get,n3bymx,vari,9,extrem,vmax *get,n0tymn,vari,2,extrem,vmin *get,n1tymn,vari,3,extrem,vmin *get,n2tymn,vari,4,extrem,vmin *get,n3tymn,vari,5,extrem,vmin *get,n0bymn,vari,6,extrem,vmin *get,n1bymn,vari,7,extrem,vmin *get,n2bymn,vari,8,extrem,vmin *get,n3bymn,vari,9,extrem,vmin *use,mxof4,n0tymx,n1tymx,n2tymx,n3tymx *use,mnof4,n0tymn,n1tymn,n2tymn,n3tymn minn=abs(minn *if,minn,gt,maxx,then maxx=minn *endif /output,racluse,sum,,append msg,info,maxx Maximum corner y-displacement, top plate: %g in /output,term *msg,info > maximum corner y-displacement, top plate complete /output,racluse,1st,,append *use,mxof4,n0bymx,n1bymx,n2bymx,n3bymx *use,mnof4,n0bymn,n1bymn,n2bymn,n3bymn minn=abs(minn) *if,minn,gt,maxx,then maxx=minn *endif /output,racluse,sum,,append *msg,info,maxx Maximum corner y-displacement, baseplate: %g in /output,term *msg,info > maximum corner y-displacement, baseplate complete /output,racluse,1st,,append /show,term finish /post1 r1p=0.0 r2p=0.0 r3p=0.0 r4p=0.0 r5p=0.0 r6p=0.0 r7p=0.0 esel,s,elem,,14,17,1 *do,cls,1,nlp,1 set,cls etable,pfn,smisc,1 etable,pfs,smisc,2 etable,pan,nmisc,9 *get,p14fn,elem,14,etab,pfn *get,p15fn,elem,15,etab,pfn *get,p16fn,elem,16,etab,pfn *get,p17fn,elem,17,etab,pfn *get,p14fs,elem,14,etab,pfs *get,p15fs,elem,15,etab,pfs *get,p16fs,elem,16,etab,pfs *get,p17fs,elem,17,etab,pfs *get,p14an,elem,14,etab,pan *get,p15an,elem,15,etab,pan *get,p16an,elem,16,etab,pan *get,p17an,elem,17,etab,pan ! start r1p sd1=abs(p14fn)/pdag sd2=abs(p15fn)/pdag sd3=abs(p16fn)/pdag sd4=abs(p17fn)/pdag rc1=sd1/fapd rc2=sd2/fapd rc3=sd3/fapd rc4=sd4/fapd *use,mxof4,rc1,rc2,rc3,rc4 r1pcv=maxx *if,maxx,gt,r1p,then r1p=maxx *endif ! start r2p svy1=abs(p14fs*cos(p14an))/pdaa svy2=abs(p15fs*cos(p15an))/pdaa svy3=abs(p16fs*cos(p16an))/pdaa svy4=abs(p17fs*cos(p17an))/pdaa rc1=svy1/fv rc2=svy2/fv rc3=svy3/fv rc4=svy4/fv *use,mxof4,rc1,rc2,rc3,rc4 r2pcur=maxx *if,maxx,gt,r2p,then r2p=maxx *endif ! start r3p (bending about y, stress along z) sby1=abs(p14fs*sin(p14an))*pdll*pdor/pdiy sby2=abs(p15fs*sin(p15an))*pdll*pdor/pdiy sby3=abs(p16fs*sin(p16an))pdll*pdor/pdiy sby4=abs(p17fs*sin(p17an))pdll*pdor/pdiy rc 1=sby1/fb, rc2=sby2/fb rc3=sby3/fb rc4=sby4/fb use,mxof4,rc1,rc2,rc3,rc4 r3pcur=maxx *if,maxx,gt,r3p,then r3p=maxx *endif ! start r4p (bending about z, stress along y) sbz1=abs(p14fs*cos(p14an))*pdll*pdor/pdiz sbz2=abs(p15fs*cos(p15an))*pdll*pdor/pdiz sbz3=abs(p16fs*cos(p16an))*pdll*pdor/pdiz sbz4=abs(p17fs*cos(p17an))*pdll*pdor/pdiz rc1=sbz1/fb rc2=sbz2/fb rc3=sbz3/fb rc4=sbz4/fb *use,mxof4,rc1,rc2.rc3,rc4 r4pcur=maxx *if,maxx,gt,r4p,then r4p=maxx *endif !start r5p *if,p14fn,lt,0.0,then dy=1-sd1/feypd dz=1-sd1/fezpd trm1=sd1/fapd trm2=cmy*sby1/dy/fb trm3=cmz*sbzi/dz/fb rc1=trm1+sqrt(trm2**2+trm2) *else rc1=0.0 *endif *if,p15fn,lt,0.0,then dy=1-sd2/feypd dz=1-sd2/fezpd trm1=sd2/fapd trm2=cmy*sby2/dy/fb trm3=cmz*sbz2/dz/fb rc2=trm1+sqrt(trm2**2+trm3**2) *else rc2=0.0 *endif *if,p16fn,lt,0.0,then dy=1-sd3/feypd dz=1-sd3/fezpd trm1=sd3/fapd trm2=cmy*sby3/dy/fb trm3=cmz*sbz3/dz/fb rc3=trm1+sqrt(trm2**2+trm3**2) *else rc3=0.0 *endif *if,p17fn,lt,0.0,then dy=1-sd3/feypd dz=1-sd3/fezpd trm1=sd3/fapd trm2=cmy*sby4/dy/fb trm3=cmz*sbz4/dz/fb rc4=trm1+sqrt(trm2**2+trm3**2) *else rc4=0.0 *endif *use,mxof4,rc1,rc2,rc3,rc4 *if,maxx,gt,r5p,then r5p=maxx *endif ! start r6p rc1=sd1/ft+sqrt((sby1/fb)**2+(sbz1/fb)**2) rc2=sd2/ft+sqrt((sby2/fb)**2+(sbz2/fb)**2) rc3=sd3/ft+sqrt((sby3/fb)**2+(sbz3/fb)**2) rc4=sd4/ft+sqrt((sby4/fb)**2+(sbz4/fb)**2) *use,mxof4,rc1,rc2,rc3,rc4 *if,maxx,gt,r6p,then r6p=maxx r6pr1p=r1pcur r6pr2p=r2pcur r6pr3p=r3pcur r6pr4p=r4pcur *endif ! start r7p svz1=abs(p14fs*sin(p14an))/pdaa svz2=abs(p15fs*sin(p15an))/pdaa svz3=abs(p16fs*sin(p16an))/pdaa svz4=abs(p17fs*sin(p17an))/pdaa rc1=svz1/fv rc2=svz2/fv rc3=svz3/fv rc4=svz4/fv *use,mxof4,rc1,rc2,rc3,rc4 *if,maxx,gt,r7p,then r7p=maxx *endif /output,term *msg,info,cls load step %i processing complete /output,racluse,1st,append *enddo /output,racluse,sum,,append *msg,info,r1p Pedestal stress factor r1: %g *msg,info,r2p Pedestal stress factor r2: %g *msg,info,r3p Pedestal stress factor r3: %g *msg,info,r4p Pedestal stress factor r4: %g *msg,info,r5p Pedestal stress factor r5: %g *msg,info,r6p Pedestal stress factor r6: %g *msg,info,r6pr1p,r6pr2p,r6pr3p,r6pr4p Corresponding r1p=%g, r2p = %g, r3p=%g r4p= %g *msg,info,r7p Pedestal stress factor r7: %g /output,term *msg,info > all processing complete /output,racluse,1st,append esel,all finish /output,term Handwritten, 2 page letter has been omitted. Volume in drive A has no label Volume Serial Number is 1058-16D1 Directory of A:\ RAO8USE INP 26,780 09-14-95 8:04a RAO9USE INP 26,780 09-18-95 5:26p RA10USE INP 26,797 09-21-95 8:17a RA11USE INP 26,789 09-28-95 5:21p ANS INP 83 10-13-95 4:28p RAO2GEN INP 8,077 09-05-95 6:41p P5EW1 ANS 48,000 09-06-94 4:35p P5EW2 ANS 48,000 09-06-94 4:35p P5EW3 ANS 48,000 09-06-94 4:36p P5NS1 ANS 48,000 09-06-94 4:36p P5NS2 ANS 48,000 09-06-94 4:36p P5NS3 ANS 48,000 09-06-94 4:36p PSVT1 ANS 48,000 09-06-94 4:37p PSVT2 ANS 48,000 09-06-94 4:37p PSVT3 ANS 48,000 09-06-94 4:37p MNOF4 868 10-26-94 2:11p MNOF8 1,900 10-26-94 2:13p MXOF4 868 10-26-94 2:11p MXOF8 1,900 10-26-94 2:14p 19 file(S) 552,842 bytes 900,096 bytes free ANSYS ANSYS, Inc. 201 Johnson Road Telephone 412.746.3304 Houston, PA 15342-1300 Facsimile 412.746.9494 October 18, 1995 Alan I. Soler Executive Vice President Holtec International 2060 Fairfax Avenue Cherry Hill, NJ Dear Mr. Soler: Thank you for your letter dated August 4, 1995 concerning an ANSYS analysis which gave you two distinct sets of answers for a problem run on several different platforms. We appreciate input from our customers concerning the performance of ANSYS. Our Technical Support Group was able to reproduce the phenomena you observed, obtaining different results with different computer systems. The reasons for these differences lie with the precision differences experienced with different computer systems coupled with the stability of your problem. These results are not typical for ANSYS analysis. But they are acceptable for this case. The reasons for these differences lies in the nature of the problem and the role which computer inaccuracies (or machine accuracy) play in the solution of it. Factors which contribute to this condition for your problem are: No damping Severe nonlinearities Path dependent solution (e.g. friction) Large number of small time steps Automatic time stepping Unconverged solutions (input included on NCNV command) Of these, your damping condition and the use of unconverged solutions contributes the most to producing the results you observed. The stability of your problem can also be observed from your results which you explained in your August 4th letter. "We have rerun the scripts using the default tolerances of .001 on force and moment and notice some differences in the results This would indicate that the problem was not very stable. Mr. Alan I. Soler Holtec International Page Two October 18, 1995 During our evaluation, we added a small amount of damping to your problem. We feel that it's addition is justifiable. The beta damping we added represents only 1% of critical damping at 5 hertz. With the addition of this small amount of damping and the use of fixed time steps, the results were essentially the same for all computer systems. Adding the LNSRCH command resulted in convergence of all load steps. Thus, by eliminating the zero damping condition and the unconverged solutions, the problem which you experienced is eliminated. One final item which you should be aware of is the possible non- conservative sampling technique which you used to analyze your results. If you are looking for peak loads it is our recommendation that you obtain solutions not only at the absolute peaks of your load steps but also in areas around those peaks. Depending on the problem, we have found that peak loads experienced by structures may lag or precede the absolute peak load. If you have any other questions concerning this issue, please do not hesitate in calling. Sincerely, ANSYS, Inc. William Bryan Quality Assurance Manager cc: MSTI M. Imgrund 2060 Fairfax Avenue, Cherry Hill, NJ 08003-1666 Telephone: (609) 424-0999 H O L T E C Telex: 910-240-6663 INTERNATIONAL Telefax: (609) 424-1710 October 23, 1995 Mr. William Bryan Q.A. Manager ANSYS, Inc. 201 Johnson Road Houston, PA 15342-1300 Dear Mr. Bryan: We are in receipt of your letter of October 18, 1995 concerning our ANSYS problem. I have the following questions which need your consideration: 1. Our analysis did have beta damping of 4% of critical at 7HZ. This damping was applied as follows: a. In the generation pass, all elements having MAT, 1 or MAT, 2 had damping imposed by MP, DAMP, 1 or 2, .0018 (look at first two pages of generation pass). b. In the use pass, the non-linear gap elements had damping associated with them thru the real constant designator constants 10, 14, 18) (see pages 2, 3 of the use pass). My question would be - since this is the maximum permitted by the specification, is the addition of more really needed - also, did we correctly impose the damping? 2. With regard to the NCNV command, we tracked iterations, and it appeared to us that 50 iterations were never exceeded on any step. Did you run a case with NCNV set to kick the system out of the solution mode and indeed find that it was not converging. 3. With regards to automatic time stepping - is it good or bad? I would think that it should not cause a problem if everything else works as it should. 4. With regard to the LNSRCH, I feel it should be the default, since it does seem to give a more stable result. H O L T E C INTERNATIONAL Mr. William Bryan ANSYS, Inc. October 23, 1995 Page 2 5. With regard to "large number of small time steps", are we to conclude that we should never use AUTOTS on a problem of this kind? There are still some unanswered questions in my mind. a. If you turn displacement convergence on and increase the number of iterations in a step, is the problem cured? b. Do we have the damping that we think we have? c. What should the default mode of running this problem be? Sincerely, Alan I. Soler Executive Vice President AIS:mp Document ID: AS124 Table "HOLTEC INTERNATIONAL DOCUMENT TRANSMITTAL FORM" omitted. Table "HOLTEC INTERNATIONAL DOCUMENT TRANSMITTAL FORM" omitted. *** END OF DOCUMENT ***

Page Last Reviewed/Updated Friday, January 31, 2020