有限元设计实例

1、设计条件某地下岩石洞室剖面如图所示,洞室序号24,岩石地弹性参数 E = 1.0 105MPa,i = 0.25,a=5.6,b=12.32(b/a=2.2), qy=6.4, qx=1.28(q%/ qy =0.2).设计要求：按平面应变问题取四分之一洞室计算岩石地位移、应变和应力，并对计算结果进行必要地讨论_J%图.计算洞室剖面示意图单元划分网格54484236输入数据文件54,80,12,10,16,1100000,0.25,0.0,15.60000,12.37274,20.50002,30.25276,41.95605,56.00000,5.60000,12.37274,20.5000

2、2,30.25276,41.95605,56.00000,5.60000,12.37274,20.50002,30.25276,41.95605,56.00000,5.60000,12.37274,20.50002,30.25276,41.95605,56.00000,5.60000,12.37274,20.50002,30.25276,41.95605,56.00000,4.20000,9.27955,15.37502,22.68957,31.46704,42.00000,2.80000,6.18637,10.25001,15.12638,20.97803,28.00000,1.40000,

3、3.09318,5.12501,7.56319,36.42492010.48901,47.86814.0000014.00000,61.60000.00000.00000,12.32000.00000.00000,18.94223.00000.00000,26.88891.00000.00000,36.42492.00000.00000,47.868143.08000.00000,61.600004.73556-9.86,116.72223-19.72,239.10623-19.72,3511.96704-19.72,4715.40000-9.86,596.16000-44.8,609.471

4、12-89.6,7213.44445-89.6,8418.21246-89.6,9623.93407-44.8,10830.800001, 2,79.240002,8,714.206672,3,820.166683,9,827.318693,4,935.901114,10,946.200004,5,1012.320005,11,1018.942235,6,1126.888916,12,1136.424927,8,1347.868148,14,1361.600008,9,1412.320009,15,1418.942239,10,1526.8889110,16,1536.4249210,11,1

5、647.8681411,17,1661.6000011,12,1712.3200012,18,1718.9422313,14,1926.8889114,20,1936.4249214,15,2047.8681415,21,2061.6000015,16,2112.3200016,22,2118.9422316,17,2226.8889117,23,2247,48,5447,54,532,4,6,8,10,1297,99,101,103,105,10717,18,2318,24,2319,20,2520,26,2520,21,2621,27,2621,22,2722,28,2722,23,282

6、3,29,2823,24,2924,30,2925,26,3225,32,3126,27,3326,33,3227,28,3427,34,3328,29,3528,35,3429,30,3629,36,3531,32,3831,38,3732,33,3932,39,3833,34,4033,40,3934,35,4134,41,4035,36,4235,42,4137,38,4437,44,4338,39,4538,45,4439,40,4639,46,4540,41,4740,47,4641,42,4841,48,4743,44,5043,50,4944,45,5144,51,5045,46

7、,5245,52,5146,47,5346,53,52输出数据文件PLANE PROBLEM/LXM=1 NJ= 54 NE= 80NZ=12NPJ=10LDD=16 E0=1.0e5 EMU= .25 OU= .00TE= 1.00NODAL CO-ORDINATESNO= 1X= 5.600 'X= 12.373Y= .0001Y= .000NO= 2NO= 3X= 20.500Y= .000NO= 4X= 30.253Y= .000NO= 5X= 41.956Y= .000NO= 6X= 56.000Y= .000NO= 7X= 5.600 Y= 3.080NO= 8X= 12.

8、373Y= 4.736NO= 9X= 20.500Y= 6.722NO= 10X= 30.253Y= 9.106NO= 11X= 41.956Y= 11.967NO= 12X= 56.000Y= 15.400NO= 13X= 5.600Y= 6.160NO= 14X= 12.373Y= 9.471NO= 15X= 20.500Y= 13.444NO= 16X= 30.253Y= 18.212NO= 17X= 41.956Y= 23.934NO= 18X= 56.000Y= 30.800NO= 19X= 5.600Y= 9.240NO= 20X= 12.373Y= 14.207NO= 21X=

9、20.500Y= 20.167NO= 22X= 30.253Y= 27.319NO= 23X= 41.956Y= 35.901NO= 24X= 56.000Y= 46.200NO= 25X= 5.600Y= 12.320NO= 26X= 12.373Y= 18.942NO= 27X= 20.500Y= 26.889NO= 28X= 30.253Y= 36.425NO= 29X= 41.956Y= 47.868NO= 30X= 56.000Y= 61.600NO= 31X= 4.200Y= 12.320NO= 32X= 9.280Y= 18.942NO= 33X= 15.375Y= 26.889

10、NO= 34X= 22.690Y= 36.425NO= 35X= 31.467Y= 47.868NO= 36X= 42.000Y= 61.600NO= 37X=2.800Y= 12.320NO= 38X=6.186Y= 18.942NO= 39X= 10.250Y= 26.889NO= 40X= 15.126Y= 36.425NO= 41X= 20.978Y= 47.868NO= 42X= 28.000Y= 61.600NO= 43X=1.400Y= 12.320NO= 44X=3.093Y= 18.942NO= 45X=5.125Y= 26.889NO= 46X=7.563Y= 36.425

11、NO= 47X= 10.489Y= 47.868NO= 48X= 14.000Y= 61.600NO= 49X=.000Y= 12.320NO= 50X=.000Y= 18.942NO= 51X=.000Y= 26.889NO= 52X=.000Y= 36.425NO= 53X=.000Y= 47.868NO= 54X= .000 Y= 61.600NODAL NUMBERS OF ELEMENTLE=11 27LE= 2287LE=32 38LE= 4398LE=53 49LE= 64109LE=74 510LE= 851110LE=95 611LE= 10(61211LE= 117 8I

12、13LE= 1281413LE= 138 9)14LE= 1491514LE= 159 10 15LE= 16101615LE= 1710 11 16LE= 18111716LE= 1911 12 17LE= 20121817LE= 2113 14 19LE= 22142019LE= 2314 15 20LE= 24152120LE=25151621LE=26162221LE=27161722LE=28172322LE=29171823LE=30182423LE=31192025LE=32202625LE=33202126LE=34212726LE=35212227LE=36222827LE=

13、37222328LE=38232928LE=39232429LE=40243029LE=41252632LE=42253231LE=43262733LE=44263332LE=45272834LE=46273433LE=47282935LE=48283534LE=49293036LE=50293635LE=51313238LE=52313837LE=53323339LE=54323938LE=55333440LE=56334039LE=57343541LE=58344140LE=59353642LE=60354241LE=61373844LE=62374443LE=63383945LE=643

14、84544LE=65394046LE=66394645LE=67404147LE=68404746LE=69414248LE=70414847LE=71434450LE=72435049LE=73444551LE=74445150LE=75454652LE=76455251LE=77464753LE=78465352LE=79474854LE=804754 53NODAL DISPLACEMENTSNO= 1 U= -.0000968V= .0000000NO= 2 U=.0000600V= .0000000NO= 3 U=.0001800V= .0000000NO= 4 U=.0002808

15、V= .0000000NO= 5 U=.0003837V= .0000000NO= 6 U= .0005031V=.0000000NO= 7U=-.0000770.0002575V=NO= 8U=.0000729.0003422V=NO= 9U=.0001810.0004303V=NO= 10U=.0002729V=.0005451NO= 11U=.0003705V=.0006920NO= 12U=.0004803.0008671V=NO= 13U=-.0000495.0005212V=NO= 14U=.0000886.0006733V=NO= 15U=.0001714.0008487V=NO

16、= 16U=.0002501.0010850V=NO= 17U=.0003419V=.0013827NO= 18U=.0004429V=.0017395NO= 19U=-.0000162V=.0007939NO= 20U=.0000882V=.0009851NO= 21U=.0001451.0012556V=NO= 22U=.0002267.0016214V=NO= 23U=.0003123.0020699V=NO= 24U=.0004018.0026101V=NO= 25U=.0000117.0010766V=NO= 26U=.0000562.0012738V=NO= 27U=.000131

17、6.0016552V=NO= 28U=.0002102.0021491V=NO= 29U=.0002837.0027495V=NO= 30U=.0003545.0034743V=NO= 31U=-.0000050.0011890V=NO= 32U=.0000217.0013378V=NO= 33U=.0000812.0017019V=NO= 34U=.0001446.0021916V=NO= 35U=.0002006.0027936V=NO= 36U=.0002427.0035205V=NO= 37U=-.0000107.0012861V=NO= 38U=-.0000030.0014211V=

18、NO= 39U=.0000373.0017605V=NO= 40U=.0000834.0022418V=NO= 41U=.0001227.0028462V=NO= 42U=.0001386.0035778V=NO= 43U=-.0000071.0013506V=NO= 44U=-.0000088.0014943V=NO= 45U=.0000100.0018186V=NO=46 U=NO=47 U=NO=48 U=NO=49 U=NO= 50 U=NO= 51 U=NO= 52 U=NO= 53 U=NO= 54 U=.0000342 V=.0022890.0000556 V=.0028949.

19、0000569 V=.0036354.0000000 V=.0013726.0000000 V=.0015254.0000000 V=.0018448.0000000 V=.0023092.0000000 V=.0029130.0000000 V=.0036641RFACTIVE FORCES OF,NODAL CONSTRAINTR( 2)=30.4445100R( 4)=59.6905900R( 6)=64.9455200R( 8)=73.3995600R( 10)=85.0538300R( 12)=44.8661200R( 97)=4.6294510R( 99)=12.2499500R(

20、101)=14.3128300R(103)=16.2350700R(105)=19.4074400R(107)=12.0452200STRESS OF ELEMENTLE= 1XX= -.566 YY= -9.107 XY= .258X1=3.995 X2=-4.561 PX=-.283CETA= .000LE= 2XX=-.314YY=-7.813XY=.316X1 =3.605X2=-3.920PX=-.157CETA= .000LE= 3XX=-1.118YY=-8.080XY=.109X1 =2.924X2=-4.042PX=-.559CETA= .000LE= 4XX=-.970YY

21、=-7.152XY=.198X1 =2.613X2=-3.582PX=-.485CETA= .000X1 =3.764X2=-4.794 PX=CETA= .000-.515LE= 12XX=-.543YY= -7.639 XY=.601X1 =3.327X2=-3.870 PX=-.272CETA= .000LE= 13XX=-1.299YY= -7.890 XY=.382X1 =2.669X2=-3.967 PX=-.649CETA= .000LE= 14XX=-1.183YY= -7.032 XY=.297X1 =2.348X2= -3.531 PX=-.592CETA= .000LE=

22、 5XX=-1.321YY= -7.269XY=.006X1 =2.314X2=-3.634PX=-.660CETA= .000LE= 15XX=-1.317YY= -7.077XY=.081X1 =2.223X2= -3.539PX=-.658CETA= .000LE= 6XX=-1.238YY=-6.798XY=.080X1 =2.162X2=-3.400PX=-.619CETA= .000LE= 7XX=-1.340YY=-6.832XY=-.035X1 =2.077X2=-3.416PX=-.670CETA= .000LE= 16XX=-1.257YY= -6.743XY=.090X1

23、 =2.116X2= -3.373PX=-.628CETA= .000LE= 17XX=-1.298YY=-6.757XY=-.022X1 =2.081X2=-3.378PX=-.649CETA= .000LE= 8XX=-1.280YY=-6.594XY=.019X1 =2.017X2=-3.297PX=-.640CETA= .000LE= 18XX=-1.228YY=-6.566XY=.016X1 =2.055X2=-3.283PX=-.614CETA= .000LE= 9XX=-1.292YY=-6.598XY=-.044X1 =2.007X2=-3.300PX=-.646CETA= .

24、000LE= 10XX=-1.270YY=-6.429XY=-.007X1 =1.944X2=-3.214PX=-.635CETA= .000LE= 19XX=-1.300YY=-6.591XY=-.030X1 =1.995X2=-3.295PX=-.650CETA= .000LE= 20XX=-1.261YY=-6.463XY=-.006X1 =1.971X2=-3.232PX=-.630CETA= .000LE= 21LE= 11XX= -1.029 YY=-9.474 XY=.694XX=-1.729 YY= -10.021 XY=1.266X1 =3.471X2=CET-5.200A=

25、.0(PX=)0-.864LE= 22XX=-.776YY=-7.282XY=.799X1 =2.962X2=-3.738PX=-.388CETA= .000LE= 23XX=-1.406YY=-7.492XY=.421X1 =2.369X2=-3.775PX=-.703CETA= .000X1 =3.151X2=-5.768 PX=CETA= .000-1.309LE= 32XX=-.856YY= -6.789 XY=.949X1 =2.686X2=-3.543 PX=-.428CETA= .000LE= 33XX=-1.004YY= -6.838 XY=.187X1 =2.421X2=-3

26、.425 PX=-.502CETA= .000XX=LE= 24-1.237 YY= -6.869 XY=2.213 X2= -3.449 PX=CETA= .000.288-.618LE= 34XX=-1.030YY=-6.685XY=.368X1 =2.337X2=-3.366PX=-.515CETA= .000LE= 25XX=-1.223YY=-6.864XY=.058X1 =2.210X2=-3.433PX=-.612CETA= .000LE= 35XX= -1.198 YY= -6.741 XY= .163X1=2.177 X2= -3.375 PX= -.599CETA= .00

27、0LE= 26XX= -1.126 YY= -6.658 XY= .124X1=2.206 X2=-3.332 PX= -.563LE= 36XX= -1.139 YY= -6.560 XY= .168X1=2.147 X2=-3.285 PX= -.569CETA= .000CETA= .000LE= 37LE= 27XX=-1.264 YY= -6.704 XY=.031XX=-1.280 YY= -6.608 XY=.094X1 =2.088 X2=-3.352 PX=-.632X1 =2.025 X2= -3.305 PX=-.640CETA= .000CETA= .000LE= 28

28、XX=-1.201YY=-6.525XY=.053X1 =2.062X2=-3.263PX=-.601CETA= .000LE= 38XX=-1.237YY=-6.469XY=.073X1 =1.999X2=-3.236PX=-.618CETA= .000LE= 29XX=-1.289YY=-6.555XY=.008X1 =1.988X2=-3.277PX=-.644CETA= .000LE= 39XX=-1.296YY=-6.489XY=.031X1 =1.948X2=-3.245PX=-.648CETA= .000LE= 30XX=-1.262YY=-6.451XY=.013X1 =1.9

29、63X2=-3.225PX=-.631CETA= .000LE= 40XX=-1.279YY=-6.412XY=.007X1 =1.927X2=-3.206PX=-.640CETA= .000LE= 41XX=-2.617 YY= -10.662 XY=1.925XX=-.698 YY= -5.665 XY=.639LE= 312.216 X2= -2.913 PX=-.349X1 =1.889 X2=-2.658 PX= -.384CETA= .000CETA= .000LE= 42XX=-1.938YY=-9.612XY=3.009X1 =3.907X2=-5.845PX=-.969CET

30、A= .000LE= 43XX=-1.113YY=-6.484XY=.342X1 =2.151X2=-3.264PX=-.556CETA= .000LE= 52XX=-1.745YY=-6.539XY=2.737X1 =2.765X2=-4.511PX=-.873CETA= .000LE= 53XX=-1.154YY=-6.207XY=.493X1 =1.998X2=-3.151PX=-.577CETA= .000LE= 44XX=-1.128YY= -6.955 XY=.784X1 =2.453X2=-3.581 PX=-.564CETA= .000LE= 45XX=-1.261YY= -6

31、.557 XY=.200X1 =2.025X2= -3.286 PX=-.630CETA= .000LE= 54XX= -1.302 YY= -6.459 XY= 1.117X1=2.159 X2=-3.461 PX= -.651CETA= .000LE= 55XX= -1.287 YY= -6.449 XY= .283X1=1.953 X2=-3.240 PX= -.643CETA= .000LE= 46XX= -1.154 YY= -6.607 XY= .329X1=2.169 X2=-3.323 PX= -.577CETA= .000LE= 56XX=-1.223YY=-6.415XY=

32、.475X1 =2.028X2=-3.250PX=-.611CETA= .000LE= 47XX=-1.321YY=-6.496XY=.101X1 =1.929X2=-3.250PX=-.660CETA= .000LE= 57XX= -1.367 YY= -6.478 XY= .169X1=1.877 X2=-3.245 PX= -.684CETA= .000LE= 48XX=-1.236YY=-6.484XY=.155X1 =2.010X2=-3.246PX=-.618CETA= .000LE= 58XX=-1.278YY=-6.422XY=.238X1 =1.944X2=-3.222PX=

33、-.639CETA= .000LE= 49XX=-1.288YY=-6.420XY=.012X1 =1.922X2=-3.210PX=-.644CETA= .000LE= 59XX=-1.351YY=-6.432XY=.058X1 =1.865X2=-3.217PX=-.676CETA= .000LE= 50XX=-1.296YY=-6.423XY=.048X1 =1.916X2=-3.212PX=-.648CETA= .000LE= 60XX=-1.342YY=-6.404XY=.095X1 =1.861X2=-3.204PX=-.671CETA= .000LE= 61LE= 51XX=-.

34、768 YY= -4.857 XY=.994XX=-1.076 YY= -3.823 XY=.954X1 =1.134X2=-2.210 PX=CETA= .000-.538LE= 62XX=-1.645YY= -4.119 XY=1.859X1 =1.410X2=-3.056 PX=-.823CETA= .000LE= 63XX=-1.302YY= -5.607 XY=.547X1 =1.570X2=-2.872 PX=-.651CETA= .000LE= 64XX=-1.651YY=-5.548XY=1.021X1 =1.374X2=-3.025PX=-.825CETA= .000LE=

35、65XX=-1.366YY=-6.179XY=.310X1 =1.744X2=-3.110PX=-.683CETA= .000X1=.300 X2= -1.611 PX= -.655CETA= .000LE= 72XX= -1.531 YY= -2.972 XY= .628X1=.190 X2=-1.721 PX= -.765CETA= .000LE= 73XX=X1 =-1.450 YY= -4.976 XY=1.060 X2= -2.510 PX=CETA= .000.279-.725LE= 74XX=-1.947 YY= -4.937 XY=.402X1 =.574 X2=-2.522

36、PX=-.974CETA= .000LE= 75XX=-1.457 YY= -5.821 XY=.162X1 =1.459 X2= -2.916 PX=-.729CETA= .000LE= 66XX=-1.451YY=-6.055XY=.501X1 =1.630X2=-3.081PX=-.726CETA= .000LE= 67XX=-1.440YY=-6.367XY=.192X1 =1.751X2=-3.191PX=-.720CETA= .000LE= 76XX=X1 =-1.713 YY= -5.766 XY=.205-.8571.180X2= -2.893 PX=CETA= .000LE=

37、 77XX=-1.500YY= -6.195XY=.090X1 =1.600X2= -3.099PX=-.750CETA= .000LE= 68XX=-1.401YY=-6.285XY=.258X1 =1.755X2=-3.156PX=-.701CETA= .000LE= 69XX=-1.515YY=-6.412XY=.092X1 =1.693X2=-3.208PX=-.757LE= 78XX=-1.567YY=-6.151XY=.107X1 =1.510X2=-3.078PX=-.784CETA= .000LE= 79XX=-1.690YY=-6.372XY=.044X1 =1.496X2=-3.186PX=-.845CETA= .000LE= 80CETA= .000LE= 70XX=-1.437YY=-6.358XY=.124X1 =1.746X2=-3.182PX=-.718CETA= .000XX=-1.552YY=-6.352XY=.069X1 =1.625X2=-3.177PX=-.776CETA= .000LE= 71THE ENDXX= -1.310 YY