华科《有限元程序设计》作业.doc

第二章2、修改FEM1,计算P31例2-2。一数据文件：E1.TXT91,144,182,8,2E10,0.167,1.0,0.0,1 单元节点编号。节点编号。1,3,5,7,9,11,13,182 二输出结果 NN NE ND NFIX E ANU T GM NTYPE 91 144 182 80.2000E+11 0.167 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.375000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 7 0.000000E+00 0.000000E+00 8 0.000000E+00 -0.750000E+01 9 0.000000E+00 0.000000E+00 10 0.000000E+00 0.000000E+00 11 0.000000E+00 0.000000E+00 12 0.000000E+00 0.000000E+00 13 0.000000E+00 0.000000E+00 14 0.000000E+00 0.000000E+00 15 0.000000E+00 -0.750000E+01 16 0.000000E+00 0.000000E+00 17 0.000000E+00 0.000000E+00 18 0.000000E+00 0.000000E+00 19 0.000000E+00 0.000000E+00 20 0.000000E+00 0.000000E+00 21 0.000000E+00 0.000000E+00 22 0.000000E+00 -0.750000E+01 23 0.000000E+00 0.000000E+00 24 0.000000E+00 0.000000E+00 25 0.000000E+00 0.000000E+00 26 0.000000E+00 0.000000E+00 27 0.000000E+00 0.000000E+00 28 0.000000E+00 0.000000E+00 29 0.000000E+00 -0.750000E+01 30 0.000000E+00 0.000000E+00 31 0.000000E+00 0.000000E+00 32 0.000000E+00 0.000000E+00 33 0.000000E+00 0.000000E+00 34 0.000000E+00 0.000000E+00 35 0.000000E+00 0.000000E+00 36 0.000000E+00 -0.750000E+01 37 0.000000E+00 0.000000E+00 38 0.000000E+00 0.000000E+00 39 0.000000E+00 0.000000E+00 40 0.000000E+00 0.000000E+00 41 0.000000E+00 0.000000E+00 42 0.000000E+00 0.000000E+00 43 0.000000E+00 -0.750000E+01 44 0.000000E+00 0.000000E+00 45 0.000000E+00 0.000000E+00 46 0.000000E+00 0.000000E+00 47 0.000000E+00 0.000000E+00 48 0.000000E+00 0.000000E+00 49 0.000000E+00 0.000000E+00 50 0.000000E+00 -0.750000E+01 51 0.000000E+00 0.000000E+00 52 0.000000E+00 0.000000E+00 53 0.000000E+00 0.000000E+00 54 0.000000E+00 0.000000E+00 55 0.000000E+00 0.000000E+00 56 0.000000E+00 0.000000E+00 57 0.000000E+00 -0.750000E+01 58 0.000000E+00 0.000000E+00 59 0.000000E+00 0.000000E+00 60 0.000000E+00 0.000000E+00 61 0.000000E+00 0.000000E+00 62 0.000000E+00 0.000000E+00 63 0.000000E+00 0.000000E+00 64 0.000000E+00 -0.750000E+01 65 0.000000E+00 0.000000E+00 66 0.000000E+00 0.000000E+00 67 0.000000E+00 0.000000E+00 68 0.000000E+00 0.000000E+00 69 0.000000E+00 0.000000E+00 70 0.000000E+00 0.000000E+00 71 0.000000E+00 -0.750000E+01 72 0.000000E+00 0.000000E+00 73 0.000000E+00 0.000000E+00 74 0.000000E+00 0.000000E+00 75 0.000000E+00 0.000000E+00 76 0.000000E+00 0.000000E+00 77 0.000000E+00 0.000000E+00 78 0.000000E+00 -0.750000E+01 79 0.000000E+00 0.000000E+00 80 0.000000E+00 0.000000E+00 81 0.000000E+00 0.000000E+00 82 0.000000E+00 0.000000E+00 83 0.000000E+00 0.000000E+00 84 0.000000E+00 0.000000E+00 85 0.000000E+00 -0.375000E+01 86 0.000000E+00 0.000000E+00 87 0.000000E+00 0.000000E+00 88 0.000000E+00 0.000000E+00 89 0.000000E+00 0.000000E+00 90 0.000000E+00 0.000000E+00 91 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.493087E-23 -0.292701E-06 2 -0.302796E-23 -0.293732E-06 3 -0.153339E-23 -0.294348E-06 4 -0.570359E-25 -0.294638E-06 5 0.141856E-23 -0.294665E-06 6 0.292623E-23 -0.294481E-06 7 0.547819E-23 -0.294144E-06 8 -0.882862E-08 -0.291130E-06 9 -0.590359E-08 -0.291777E-06 10 -0.298383E-08 -0.292165E-06 11 -0.812352E-10 -0.292297E-06 12 0.283649E-08 -0.292166E-06 13 0.583874E-08 -0.291743E-06 14 0.910186E-08 -0.290941E-06 15 -0.176179E-07 -0.284402E-06 16 -0.116677E-07 -0.284969E-06 17 -0.583989E-08 -0.285259E-06 18 -0.634808E-10 -0.285290E-06 19 0.574333E-08 -0.285058E-06 20 0.116770E-07 -0.284542E-06 21 0.178212E-07 -0.283733E-06 22 -0.261524E-07 -0.272982E-06 23 -0.172542E-07 -0.273504E-06 24 -0.858270E-08 -0.273746E-06 25 -0.105726E-10 -0.273733E-06 26 0.858482E-08 -0.273471E-06 27 0.173217E-07 -0.272953E-06 28 0.263089E-07 -0.272163E-06 29 -0.342899E-07 -0.257136E-06 30 -0.225869E-07 -0.257611E-06 31 -0.111988E-07 -0.257820E-06 32 0.413971E-10 -0.257790E-06 33 0.112943E-07 -0.257528E-06 34 0.227154E-07 -0.257030E-06 35 0.344556E-07 -0.256273E-06 36 -0.419151E-07 -0.237119E-06 37 -0.275860E-07 -0.237539E-06 38 -0.136539E-07 -0.237715E-06 39 0.856123E-10 -0.237674E-06 40 0.138326E-07 -0.237424E-06 41 0.277821E-07 -0.236956E-06 42 0.421258E-07 -0.236246E-06 43 -0.489111E-07 -0.213219E-06 44 -0.321693E-07 -0.213572E-06 45 -0.159078E-07 -0.213711E-06 46 0.121169E-09 -0.213662E-06 47 0.161598E-07 -0.213431E-06 48 0.324407E-07 -0.213004E-06 49 0.491889E-07 -0.212355E-06 50 -0.551484E-07 -0.185773E-06 51 -0.362511E-07 -0.186046E-06 52 -0.179224E-07 -0.186143E-06 53 0.140483E-09 -0.186092E-06 54 0.182286E-07 -0.185892E-06 55 0.366100E-07 -0.185521E-06 56 0.555277E-07 -0.184948E-06 57 -0.604758E-07 -0.155183E-06 58 -0.397401E-07 -0.155359E-06 59 -0.196664E-07 -0.155411E-06 60 0.121612E-09 -0.155362E-06 61 0.199782E-07 -0.155207E-06 62 0.402102E-07 -0.154912E-06 63 0.610538E-07 -0.154432E-06 64 -0.647076E-07 -0.121947E-06 65 -0.425367E-07 -0.122006E-06 66 -0.211158E-07 -0.121993E-06 67 0.254625E-10 -0.121930E-06 68 0.213273E-07 -0.121809E-06 69 0.431601E-07 -0.121597E-06 70 0.657402E-07 -0.121227E-06 71 -0.676296E-07 -0.867390E-07 72 -0.445387E-07 -0.866516E-07 73 -0.222398E-07 -0.865094E-07 74 -0.181988E-09 -0.863209E-07 75 0.221824E-07 -0.861103E-07 76 0.453675E-07 -0.858839E-07 77 0.696715E-07 -0.855926E-07 78 -0.691036E-07 -0.505612E-07 79 -0.456703E-07 -0.503107E-07 80 -0.229548E-07 -0.498750E-07 81 -0.439547E-09 -0.492347E-07 82 0.225301E-07 -0.484400E-07 83 0.467512E-07 -0.476218E-07 84 0.730765E-07 -0.469748E-07 85 -0.694272E-07 -0.148263E-07 86 -0.459715E-07 -0.145085E-07 87 -0.231401E-07 -0.137433E-07 88 -0.460978E-09 -0.121744E-07 89 0.227734E-07 -0.953062E-08 90 0.476889E-07 -0.558120E-08 91 0.766665E-07 -0.492018E-23ELEMENT X-STR Y-STR XY-STR1 -0.235102E+03 0.196460E+01 0.179446E+022 -0.157505E+03 -0.457531E+00 -0.277912E+023 -0.157708E+03 -0.167097E+01 0.223379E+024 -0.791844E+02 0.230318E+01 -0.250900E+025 -0.798635E+02 -0.176499E+01 0.249486E+026 -0.132146E+01 0.505946E+01 -0.229986E+027 -0.204379E+01 0.733173E+00 0.267462E+028 0.769114E+02 0.761440E+01 -0.214570E+029 0.765516E+02 0.545906E+01 0.285467E+0210 0.157258E+03 0.933106E+01 -0.201619E+0211 0.157851E+03 0.128788E+02 0.312907E+0212 0.244167E+03 0.868761E+01 -0.193180E+0213 -0.236666E+03 -0.136781E+02 0.267470E+0214 -0.154222E+03 -0.306593E+01 -0.241972E+0215 -0.155452E+03 -0.104350E+02 0.277389E+0216 -0.763528E+02 -0.114424E+01 -0.209798E+0217 -0.774396E+02 -0.765279E+01 0.291527E+0218 0.699715E+00 0.135491E+01 -0.189443E+0219 -0.411326E+00 -0.529871E+01 0.300479E+0220 0.781434E+02 0.375863E+01 -0.183050E+0221 0.768311E+02 -0.410073E+01 0.297593E+0222 0.156609E+03 0.552584E+01 -0.194234E+0223 0.154641E+03 -0.626368E+01 0.263445E+0224 0.233628E+03 0.666478E+01 -0.229558E+0225 -0.230217E+03 -0.157576E+02 0.284958E+0226 -0.149664E+03 -0.413972E+01 -0.215026E+0227 -0.151253E+03 -0.136528E+02 0.311183E+0228 -0.735740E+02 -0.258980E+01 -0.170711E+0229 -0.750273E+02 -0.112923E+02 0.325445E+0230 0.136307E+01 -0.286377E+00 -0.148682E+0231 -0.144906E+00 -0.931585E+01 0.325239E+0232 0.761457E+02 0.223235E+01 -0.149260E+0233 0.744031E+02 -0.820199E+01 0.306905E+0234 0.151287E+03 0.455364E+01 -0.173275E+0235 0.149287E+03 -0.742037E+01 0.271067E+0236 0.227401E+03 0.635939E+01 -0.218262E+0237 -0.219642E+03 -0.158264E+02 0.285530E+0238 -0.143023E+03 -0.488657E+01 -0.189837E+0239 -0.144621E+03 -0.144540E+02 0.329689E+0240 -0.703245E+02 -0.337098E+01 -0.132118E+0241 -0.718512E+02 -0.125127E+02 0.350499E+0242 0.121490E+01 -0.102348E+01 -0.104710E+0243 -0.375403E+00 -0.105466E+02 0.348550E+0244 0.725290E+02 0.165516E+01 -0.106971E+0245 0.707672E+02 -0.889348E+01 0.324216E+0246 0.144533E+03 0.420781E+01 -0.138029E+0247 0.142525E+03 -0.781439E+01 0.279098E+0248 0.218276E+03 0.617208E+01 -0.196537E+0249 -0.205910E+03 -0.153895E+02 0.281361E+0250 -0.134250E+03 -0.564081E+01 -0.162380E+0251 -0.135694E+03 -0.142882E+02 0.341637E+0252 -0.661368E+02 -0.399002E+01 -0.905710E+0153 -0.675595E+02 -0.125089E+02 0.370772E+0254 0.931111E+00 -0.148364E+01 -0.563873E+0155 -0.583554E+00 -0.105537E+02 0.369771E+0256 0.679125E+02 0.133486E+01 -0.589620E+0157 0.662080E+02 -0.887195E+01 0.339661E+0258 0.135772E+03 0.395353E+01 -0.970998E+0159 0.133786E+03 -0.793730E+01 0.281516E+0260 0.205529E+03 0.592302E+01 -0.170046E+0261 -0.189029E+03 -0.147882E+02 0.274933E+0262 -0.123305E+03 -0.648990E+01 -0.130908E+0263 -0.124515E+03 -0.137391E+02 0.350625E+0264 -0.608717E+02 -0.459789E+01 -0.443387E+0165 -0.621098E+02 -0.120119E+02 0.387876E+0266 0.641255E+00 -0.183853E+01 -0.362068E+0067 -0.743733E+00 -0.101316E+02 0.387457E+0268 0.622501E+02 0.114850E+01 -0.745073E+0069 0.606225E+02 -0.859756E+01 0.350605E+0270 0.124860E+03 0.377640E+01 -0.536556E+0171 0.122914E+03 -0.787540E+01 0.278320E+0272 0.189290E+03 0.563025E+01 -0.140676E+0273 -0.168677E+03 -0.140673E+02 0.266583E+0274 -0.110099E+03 -0.748490E+01 -0.936801E+0175 -0.111015E+03 -0.129720E+02 0.358044E+0276 -0.545907E+02 -0.519900E+01 0.847212E+0077 -0.555978E+02 -0.112302E+02 0.402603E+0278 0.177752E+00 -0.201940E+01 0.543351E+0179 -0.105874E+01 -0.942319E+01 0.401243E+0280 0.553727E+02 0.123115E+01 0.464918E+0181 0.538165E+02 -0.808697E+01 0.356225E+0282 0.111826E+03 0.385724E+01 -0.102116E+0183 0.109908E+03 -0.762601E+01 0.269671E+0284 0.169941E+03 0.543703E+01 -0.110806E+0285 -0.144266E+03 -0.131911E+02 0.256407E+0286 -0.944957E+02 -0.872626E+01 -0.476782E+0187 -0.950346E+02 -0.119536E+02 0.364830E+0288 -0.474860E+02 -0.586206E+01 0.710692E+0189 -0.481935E+02 -0.100979E+02 0.415676E+0290 -0.852685E+00 -0.209350E+01 0.119740E+0291 -0.189467E+01 -0.833265E+01 0.411068E+0292 0.469251E+02 0.161314E+01 0.102872E+0293 0.454487E+02 -0.722774E+01 0.355691E+0294 0.967355E+02 0.437142E+01 0.298641E+0195 0.948183E+02 -0.710870E+01 0.255104E+0296 0.148296E+03 0.557954E+01 -0.856926E+0197 -0.114873E+03 -0.121288E+02 0.243627E+0298 -0.763154E+02 -0.104014E+02 0.111397E+0199 -0.763626E+02 -0.106834E+02 0.370550E+02100 -0.398448E+02 -0.714367E+01 0.146968E+02101 -0.400957E+02 -0.864673E+01 0.426809E+02102 -0.307591E+01 -0.306545E+01 0.196600E+02103 -0.370664E+01 -0.684205E+01 0.416778E+02104 0.361817E+02 0.122671E+01 0.165065E+02105 0.349846E+02 -0.594029E+01 0.348420E+02106 0.794619E+02 0.477417E+01 0.646860E+01107 0.776256E+02 -0.622200E+01 0.234206E+02108 0.126019E+03 0.627213E+01 -0.759995E+01109 -0.797547E+02 -0.109759E+02 0.222984E+02110 -0.555194E+02 -0.127711E+02 0.819923E+01111 -0.550023E+02 -0.967474E+01 0.368195E+02112 -0.318093E+02 -0.109996E+02 0.232564E+02113 -0.312706E+02 -0.777381E+01 0.430957E+02114 -0.698591E+01 -0.870597E+01 0.288131E+02115 -0.651802E+01 -0.590428E+01 0.417690E+02116 0.220097E+02 -0.474860E+01 0.245896E+02117 0.219974E+02 -0.482339E+01 0.337007E+02118 0.589974E+02 0.796125E+00 0.106848E+02119 0.580153E+02 -0.508471E+01 0.210524E+02120 0.105840E+03 0.602229E+01 -0.938215E+01121 -0.410355E+02 -0.103523E+02 0.176122E+02122 -0.327636E+02 -0.154937E+02 0.136073E+02123 -0.320189E+02 -0.110343E+02 0.330478E+02124 -0.226070E+02 -0.212029E+02 0.292610E+02125 -0.209082E+02 -0.110309E+02 0.405333E+02126 -0.114652E+02 -0.275266E+02 0.378563E+02127 -0.851257E+01 -0.984603E+01 0.404412E+02128 0.407786E+01 -0.311056E+02 0.367410E+02129 0.798272E+01 -0.772332E+01 0.330494E+02130 0.323347E+02 -0.273295E+02 0.220567E+02131 0.359555E+02 -0.564842E+01 0.206355E+02132 0.889595E+02 -0.110236E+02 -0.994203E+01133 -0.105990E+02 -0.117920E+02 0.668369E+01134 -0.104442E+02 -0.144557E+02 0.706642E+01135 -0.112543E+02 -0.193069E+02 0.197522E+02136 -0.103426E+02 -0.323346E+02 0.215326E+02137 -0.948443E+01 -0.271958E+02 0.269498E+02138 -0.113690E+02 -0.646560E+02 0.347504E+02139 -0.604855E+01 -0.327967E+02 0.297715E+02140 -0.114945E+02 -0.107671E+03 0.463622E+02141 0.104988E+01 -0.325543E+02 0.294520E+02142 -0.141647E+01 -0.158214E+03 0.533252E+02143 0.212767E+02 -0.223266E+02 0.291645E+02144 0.601269E+02 -0.213207E+03 0.400848E+023、以例1为对象，研究单元细分对计算结果的影响。数据文件：15,16,30,6,1.0E0,0.0,1.0,0.0,1 3,1,25,2,43,2,56,3,58,4,75,4,89,5,86,5,910,6,912,7,118,7,1213,8,129,8,1314,9,1310,9,1415,10,140.0,2.00.0,1.50.5,1.50.0,1.00.5,1.01.0,1.00.0,0.50.5,0,5 1.0,0.51.5,0.50.0,0.00.5,0.01.0,0.01.5,0.02.0,0.01,3,7,13,21,22,24,26,28,30 输出结果： NN NE ND NFIX E ANU T GM NTYPE 15 16 30 100.1000E+01 0.000 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 -0.100000E+01 2 0.000000E+00 0.000000E+00 3 0.000000E+00 0.000000E+00 4 0.000000E+00 0.000000E+00 5 0.000000E+00 0.000000E+00 6 0.000000E+00 0.000000E+00 7 0.000000E+00 0.000000E+00 8 0.000000E+00 0.000000E+00 9 0.000000E+00 0.000000E+00 10 0.000000E+00 0.000000E+00 11 0.000000E+00 0.000000E+00 12 0.000000E+00 0.000000E+00 13 0.000000E+00 0.000000E+00 14 0.000000E+00 0.000000E+00 15 0.000000E+00 0.000000E+00 NODE X-DISP Y-DISP 1 -0.105186E-14 -0.443779E+01 2 -0.605738E-16 -0.243779E+01 3 -0.253047E+00 -0.138593E+01 4 0.761151E-16 -0.133863E+01 5 0.490038E-01 -0.935510E+00 6 -0.353913E-01 -0.371081E+00 7 0.111832E-15 -0.601012E+00 8 0.132137E+00 -0.458571E+00 9 0.174399E+00 -0.216827E+00 10 0.165862E+00 -0.240955E-01 11 0.958747E-16 -0.400674E-15 12 0.143812E+00 -0.307660E-15 13 0.228079E+00 -0.151552E-15 14 0.245153E+00 -0.203323E-16 15 0.245153E+00 0.792914E-16 ELEMENT X-STR Y-STR XY-STR 1 -0.506094E+00 -0.400000E+01 0.105186E+01 2 0.980076E-01 -0.219833E+01 0.403116E+00 3 -0.506094E+00 -0.900836E+00 0.749811E+00 4 -0.168790E+00 -0.900836E+00 0.262378E+00 5 0.264275E+00 -0.147523E+01 0.142440E+00 6 0.980076E-01 -0.953877E+00 0.319983E+00 7 0.845234E-01 -0.953877E+00 0.158611E+00 8 -0.168790E+00 -0.308508E+00 0.354638E+00 9 -0.170745E-01 -0.308508E+00 -0.170584E-01 10 0.287624E+00 -0.120202E+01 0.108971E-15 11 0.264275E+00 -0.917143E+00 0.130766E+00 12 0.168533E+00 -0.917143E+00 -0.116748E-01 13 0.845234E-01 -0.433655E+00 0.188064E+00 14 0.341489E-01 -0.433655E+00 -0.536798E-01 15 -0.170745E-01 -0.481910E-01 0.113440E+0016 0.000000E+00 -0.481910E-01 -0.792914E-014、用程序完成习题3 (取半结构) 一、建立数据文件5,3,10,4,2.0E11,0.3,0.1,0.0,11,2,32,4,34,2,50.0,0.01.414,0.00.707,0.7071.414,1.4142.121,0.7071,2,3,4 输出结果： NN NE ND NFIX E ANU T GM NTYPE 5 3 10 40.2000E+12 0.300 0.1000.0000E+00 1 NODE X-LOAD Y-LOAD 1 0.000000E+00 0.000000E+00 2 0.000000E+