sp.out

阐述了建立物质本构关系时应满足的基本公理, 较详细地推导了各种类型粘 弹性物质的本构关系, 论述了变率理论与泛函理论之间的关系, 推导了对于粘弹性问题 求解具有实用意义的对应定理, 并以简单的例子

Test PS amplitudes                                          

  Number of layers              3
  Elastic mode                  3

  
    vp        vs        rho       dz        delta    epsilon

   6000.00   3500.00    3.0000   1000.00    0.0000    0.0000
   7000.00   4200.00    4.0000   1000.00    0.0000    0.0000
   6000.00   3500.00    3.0000   1000.00    0.0000    0.0000
  layer il =      1
  i =     1 sins =    0.00000   0.00000   0.00000   0.00000
  i =     1 coss =    1.00000   1.00000   1.00000   1.00000
  i =     2 sins =    0.00857   0.01000   0.00500   0.00600
  i =     2 coss =    0.99996   0.99995   0.99999   0.99998
  i =     3 sins =    0.01714   0.02000   0.01000   0.01200
  i =     3 coss =    0.99985   0.99980   0.99995   0.99993
  i =     4 sins =    0.02571   0.03000   0.01500   0.01800
  i =     4 coss =    0.99967   0.99955   0.99989   0.99984
  i =     5 sins =    0.03429   0.04000   0.02000   0.02400
  i =     5 coss =    0.99941   0.99920   0.99980   0.99971
  i =     6 sins =    0.04286   0.05000   0.02500   0.03000
  i =     6 coss =    0.99908   0.99875   0.99969   0.99955
  i =     7 sins =    0.05143   0.06000   0.03000   0.03600
  i =     7 coss =    0.99868   0.99820   0.99955   0.99935
  i =     8 sins =    0.06000   0.07000   0.03500   0.04200
  i =     8 coss =    0.99820   0.99755   0.99939   0.99912
  i =     9 sins =    0.06857   0.08000   0.04000   0.04800
  i =     9 coss =    0.99765   0.99679   0.99920   0.99885
  i =    10 sins =    0.07714   0.09000   0.04500   0.05400
  i =    10 coss =    0.99702   0.99594   0.99899   0.99854
  i =    11 sins =    0.08571   0.10000   0.05000   0.06000
  i =    11 coss =    0.99632   0.99499   0.99875   0.99820
  i =    12 sins =    0.09429   0.11000   0.05500   0.06600
  i =    12 coss =    0.99555   0.99393   0.99849   0.99782
  i =    13 sins =    0.10286   0.12000   0.06000   0.07200
  i =    13 coss =    0.99470   0.99277   0.99820   0.99740
  i =    14 sins =    0.11143   0.13000   0.06500   0.07800
  i =    14 coss =    0.99377   0.99151   0.99789   0.99695
  i =    15 sins =    0.12000   0.14000   0.07000   0.08400
  i =    15 coss =    0.99277   0.99015   0.99755   0.99647
  i =    16 sins =    0.12857   0.15000   0.07500   0.09000
  i =    16 coss =    0.99170   0.98869   0.99718   0.99594
  i =    17 sins =    0.13714   0.16000   0.08000   0.09600
  i =    17 coss =    0.99055   0.98712   0.99679   0.99538
  i =    18 sins =    0.14571   0.17000   0.08500   0.10200
  i =    18 coss =    0.98933   0.98544   0.99638   0.99478
  i =    19 sins =    0.15429   0.18000   0.09000   0.10800
  i =    19 coss =    0.98803   0.98367   0.99594   0.99415
  i =    20 sins =    0.16286   0.19000   0.09500   0.11400
  i =    20 coss =    0.98665   0.98178   0.99548   0.99348
  amps - 

   0.00000  -0.00255  -0.00510  -0.00765  -0.01020
  -0.01275  -0.01530  -0.01784  -0.02038  -0.02291
  -0.02544  -0.02797  -0.03049  -0.03301  -0.03552
  -0.03802  -0.04052  -0.04301  -0.04549  -0.04796
  layer il =      2
  i =     1 sins =    0.00000   0.00000   0.00000   0.00000
  i =     1 coss =    1.00000   1.00000   1.00000   1.00000
  i =     2 sins =    0.00833   0.00714   0.00500   0.00417
  i =     2 coss =    0.99997   0.99997   0.99999   0.99999
  i =     3 sins =    0.01667   0.01429   0.01000   0.00833
  i =     3 coss =    0.99986   0.99990   0.99995   0.99997
  i =     4 sins =    0.02500   0.02143   0.01500   0.01250
  i =     4 coss =    0.99969   0.99977   0.99989   0.99992
  i =     5 sins =    0.03333   0.02857   0.02000   0.01667
  i =     5 coss =    0.99944   0.99959   0.99980   0.99986
  i =     6 sins =    0.04167   0.03571   0.02500   0.02083
  i =     6 coss =    0.99913   0.99936   0.99969   0.99978
  i =     7 sins =    0.05000   0.04286   0.03000   0.02500
  i =     7 coss =    0.99875   0.99908   0.99955   0.99969
  i =     8 sins =    0.05833   0.05000   0.03500   0.02917
  i =     8 coss =    0.99830   0.99875   0.99939   0.99957
  i =     9 sins =    0.06667   0.05714   0.04000   0.03333
  i =     9 coss =    0.99778   0.99837   0.99920   0.99944
  i =    10 sins =    0.07500   0.06429   0.04500   0.03750
  i =    10 coss =    0.99718   0.99793   0.99899   0.99930
  i =    11 sins =    0.08333   0.07143   0.05000   0.04167
  i =    11 coss =    0.99652   0.99745   0.99875   0.99913
  i =    12 sins =    0.09167   0.07857   0.05500   0.04583
  i =    12 coss =    0.99579   0.99691   0.99849   0.99895
  i =    13 sins =    0.10000   0.08571   0.06000   0.05000
  i =    13 coss =    0.99499   0.99632   0.99820   0.99875
  i =    14 sins =    0.10833   0.09286   0.06500   0.05417
  i =    14 coss =    0.99411   0.99568   0.99789   0.99853
  i =    15 sins =    0.11667   0.10000   0.07000   0.05833
  i =    15 coss =    0.99317   0.99499   0.99755   0.99830
  i =    16 sins =    0.12500   0.10714   0.07500   0.06250
  i =    16 coss =    0.99216   0.99424   0.99718   0.99804
  i =    17 sins =    0.13333   0.11429   0.08000   0.06667
  i =    17 coss =    0.99107   0.99345   0.99679   0.99778
  i =    18 sins =    0.14167   0.12143   0.08500   0.07083
  i =    18 coss =    0.98991   0.99260   0.99638   0.99749
  i =    19 sins =    0.15000   0.12857   0.09000   0.07500
  i =    19 coss =    0.98869   0.99170   0.99594   0.99718
  i =    20 sins =    0.15833   0.13571   0.09500   0.07917
  i =    20 coss =    0.98739   0.99075   0.99548   0.99686
  amps - 

   0.00000   0.00265   0.00531   0.00797   0.01062
   0.01327   0.01592   0.01856   0.02120   0.02384
   0.02648   0.02910   0.03173   0.03435   0.03696
   0.03956   0.04216   0.04475   0.04733   0.04991
  layer il =      3
  i =     1 sins =    0.00000   0.00000   0.00000   0.00000
  i =     1 coss =    1.00000   1.00000   1.00000   1.00000
  i =     2 sins =    0.00714   0.00714   0.00417   0.00417
  i =     2 coss =    0.99997   0.99997   0.99999   0.99999
  i =     3 sins =    0.01429   0.01429   0.00833   0.00833
  i =     3 coss =    0.99990   0.99990   0.99997   0.99997
  i =     4 sins =    0.02143   0.02143   0.01250   0.01250
  i =     4 coss =    0.99977   0.99977   0.99992   0.99992
  i =     5 sins =    0.02857   0.02857   0.01667   0.01667
  i =     5 coss =    0.99959   0.99959   0.99986   0.99986
  i =     6 sins =    0.03571   0.03571   0.02083   0.02083
  i =     6 coss =    0.99936   0.99936   0.99978   0.99978
  i =     7 sins =    0.04286   0.04286   0.02500   0.02500
  i =     7 coss =    0.99908   0.99908   0.99969   0.99969
  i =     8 sins =    0.05000   0.05000   0.02917   0.02917
  i =     8 coss =    0.99875   0.99875   0.99957   0.99957
  i =     9 sins =    0.05714   0.05714   0.03333   0.03333
  i =     9 coss =    0.99837   0.99837   0.99944   0.99944
  i =    10 sins =    0.06429   0.06429   0.03750   0.03750
  i =    10 coss =    0.99793   0.99793   0.99930   0.99930
  i =    11 sins =    0.07143   0.07143   0.04167   0.04167
  i =    11 coss =    0.99745   0.99745   0.99913   0.99913
  i =    12 sins =    0.07857   0.07857   0.04583   0.04583
  i =    12 coss =    0.99691   0.99691   0.99895   0.99895
  i =    13 sins =    0.08571   0.08571   0.05000   0.05000
  i =    13 coss =    0.99632   0.99632   0.99875   0.99875
  i =    14 sins =    0.09286   0.09286   0.05417   0.05417
  i =    14 coss =    0.99568   0.99568   0.99853   0.99853
  i =    15 sins =    0.10000   0.10000   0.05833   0.05833
  i =    15 coss =    0.99499   0.99499   0.99830   0.99830
  i =    16 sins =    0.10714   0.10714   0.06250   0.06250
  i =    16 coss =    0.99424   0.99424   0.99804   0.99804
  i =    17 sins =    0.11429   0.11429   0.06667   0.06667
  i =    17 coss =    0.99345   0.99345   0.99778   0.99778
  i =    18 sins =    0.12143   0.12143   0.07083   0.07083
  i =    18 coss =    0.99260   0.99260   0.99749   0.99749
  i =    19 sins =    0.12857   0.12857   0.07500   0.07500
  i =    19 coss =    0.99170   0.99170   0.99718   0.99718
  i =    20 sins =    0.13571   0.13571   0.07917   0.07917
  i =    20 coss =    0.99075   0.99075   0.99686   0.99686
  amps - 

   0.00000  -0.00128  -0.00255  -0.00383  -0.00510
  -0.00638  -0.00765  -0.00892  -0.01020  -0.01147
  -0.01274  -0.01402  -0.01529  -0.01656  -0.01783
  -0.01910  -0.02037  -0.02163  -0.02290  -0.02416
 
 layer number      1  np =    100

 
 Hyperbolic fit:  vf =     4.71 err =   0.0009

  
     p         x          t       rc        sin(th)

    0.0000      0.00    0.4524  0.000000  0.000000
    0.0014     13.57    0.4524  -.002550  0.005000
    0.0029     27.15    0.4524  -.005103  0.010000
    0.0043     40.72    0.4525  -.007654  0.015000
    0.0057     54.31    0.4525  -.010202  0.020000
    0.0071     67.90    0.4526  -.012750  0.025000
    0.0086     81.51    0.4527  -.015296  0.030000
    0.0100     95.13    0.4529  -.017837  0.035000
    0.0114    108.77    0.4530  -.020376  0.040000
    0.0129    122.42    0.4532  -.022911  0.045000
    0.0143    136.09    0.4534  -.025442  0.050000
    0.0157    149.79    0.4536  -.027969  0.055000
    0.0171    163.51    0.4538  -.030491  0.060000
    0.0186    177.26    0.4540  -.033007  0.065000
    0.0200    191.05    0.4543  -.035517  0.070000
    0.0214    204.86    0.4546  -.038020  0.075000
    0.0229    218.71    0.4549  -.040517  0.080000
    0.0243    232.60    0.4552  -.043006  0.085000
    0.0257    246.52    0.4556  -.045487  0.090000
    0.0271    260.49    0.4559  -.047960  0.095000
    0.0286    274.51    0.4563  -.050424  0.100000
    0.0300    288.57    0.4567  -.052878  0.105000
    0.0314    302.69    0.4572  -.055322  0.110000
    0.0329    316.86    0.4576  -.057756  0.115000
    0.0343    331.08    0.4581  -.060178  0.120000
    0.0357    345.37    0.4586  -.062589  0.125000
    0.0371    359.72    0.4591  -.064987  0.130000