simwe论坛-关于flac+ansys
关于surfer绘制等直线，由于需要结合自己实际需要绘制，sufer只能绘制2D等直线，所以首先确定自己绘制等直线的面，然后flac3d导出的单元信息（坐标及应力，应变），通过简单的fortran程序找出所要求的面上的单元的，并记录其坐标位置，即应力应变值。输出格式为surfer认可的格式即可以了。
授人以鱼，不如授人以渔。我编的程序只能针对自己的问题，所以大家不要关心程序，如果，确实需要，按照上述思路实现应该很容易。
建立模型以后直接solve让程序产生应力分布!虽然时间花费上多一点，可是不用那么伤脑筋了。
我得懒人做法是：1）现把模量×10，c，Fi取很大，然后solve for balance；2）ini disp＝0，ini state ＝0，；3）回复原来得E，G，C，Fi等，
subzone和zone哪个是“最小几何单位”？二者有和区别呢？
我们通常说的“单元element”指的是哪一个？
用户接触到的最小几何单位应该是zone
通常说的element也是指zone
subzone和FLAC3d的Mixed Discretization算法有关
比如说一个brick
计算时又自动划分为2套overlay
每套overlay中各有5个tetra
即subzone
据说这样可以更好地模拟材料的塑性变形
这是flac3d为了避免求解中单元出现“沙漏”现象而采取的一种单元离散技术
Q:gen zone brick size 6 6 6
model elas
prop bulk 1e8 shear 7e7
fix x range x -0.1 0.1
table 1 0,0 100,-1e5 /我不知道table在这里什么意思,有什么作用
apply sxx 1.0 hist table 1 range x 5.9,6.1 y 0,6 z 0,2
hist zone sxx 6,0,0
step 100
我对这个table不是很清楚
gen zone cyl p0 0 0 0 p1 1 0 0 p2 0 2 0 p3 0 0 1 size 4 5 4
gen zone reflect norm 1,0,0
gen zone reflect norm 0,0,1
model ss
prop bulk 1.19e10 shear 1.1e10
prop coh 2.72e5 fric 44 ten 2e5
prop ttab 1
table 1 0,2e5 2e-5,0
请问这里的table和后面的数字到底是什么意思啊?
A:table相当于数已经赋值的数组
请问sel cable id=1 begin=(0.0, 0.5, 3.5) end=(8, 0.5, 3.5) nseg=8 中的nseg是什么意思?
把cable分为几段,比如nseg=8 ,就是分为八段.
gen zone reflect dip 90 dd 270 origin 0,0,0
应该是以90 和270度还有原点组成的平面进行对称拷贝吧？？
ini state＝0
初始化单元为弹性。一般来说，mc材料在初始化应力的时候，由于是手动设置的k0，因此局部单元可能已经进入了塑性状态，这会影响后继的计算步骤，因此需要ini state＝0
我这里编写了一个ansys to flac3d的前处理小程序，用起来比较方便。希望能和大家共享和交流。（可以自动形成flac开挖信息）
ansys-flac2002转换程序使用说明:
本程序用于将ansys 生成的数据格式转换为flac计算所需的数据格式。程序需要输入5个文件名。
1。单元信息文件，格式如下：
总单元数
单元信息。。。，例如：
10（总单元数）
1 1 1 1 0 1 8123 8713 8816 881
2 8122 8826 8831 8827
2 1 1 1 0 1 8122 8826 8831 8827 8127 8687 7791 8456
3 1 1 1 0 1 8812 8816 8817 8813 8827 8831 8832 8828
4 1 1 1 0 1 8827 8831 8832 8828 8456 7791 7790 8455
5 1 1 1 0 1 8813 8817 8818 8814 8828 8832 8833 8829
6 1 1 1 0 1 8828 8832 8833 8829 8455 7790 7789 8454
7 1 1 1 0 1 8814 8818 8819 8815 8829 8833 8834 8830
8 1 1 1 0 1 8829 8833 8834 8830 8454 7789 7788 8453
9 1 1 1 0 1 8815 8819 7629 7632 8830 8834 7699 7704
10 1 1 1 0 1 8830 8834 7699 7704 8453 7788 7700 7703
2。节点信息
节点总数
节点信息。。。。例如：
20（节点总数）
1 210.750000000 228.398805249 -327.458594084
2 210.750000000 227.440685854 -299.440511342
3 210.750000000 226.11835406
4 -279.017920414
4 241.820265858 240.000000000 -331.304948353
5 246.293378743 240.000000000 -331.239755597
6 252.148198534 240.000000000 -331.154425346
7 258.555372995 240.000000000 -331.061044881
8 266.830116362 240.000000000 -330.940445794
9 274.869142734 240.000000000 -330.823282132
10 284.683736459 240.000000000 -330.680240712
11 295.556691364 240.000000000 -330.521774355
12 306.689058734 240.000000000 -330.359527228
13 318.078795786 240.000000000 -330.193529103
14 329.460877328 240.000000000 -330.027642551
15 340.835169125 240.000000000 -329.861869531
16 352.899895166 240.000000000 -329.686033873
17 367.026374896 240.000000000 -329.480149475
18 380.592277707 240.000000000 -329.282435127
19 390.569611461 240.000000000 -329.137021875
20 242.519404127 240.000000000 -307.394667740
3。输入开挖信息
开挖步数
第一步开挖单元数
开挖单元信息
第二步开挖单元数
开挖单元信息
。。。。。
例如
2（两步开挖）
5（第一步开挖5个单元）
12527 5 1 1 0 1 17639 17740 8430 8221 17891 17869 8457 8689
12528 5 1 1 0 1 17891 17869 8457 8689 17894 17872 8127 8126
12529 5 1 1 0 1 17638 17739 17740 17639 17890 17868 17869 17891
12530 5 1 1 0 1 17890 17868 17869 17891 17893 17871 17872 17894
12531 5 1 1 0 1 17637 17738 17739 17638 17889 17867 17868 17890
15（第二步开挖15个单元）
12532 5 1 1 0 1 17889 17867 17868 17890 17892 17870 17871 17893
12533 5 1 1 0 1 17641 17742 17738 17637 17887 17873 17867 17889
12534 5 1 1 0 1 17887 17873 17867 17889 17888 17874 17870 17892
12535 5 1 1 0 1 17643 17730 17742 17641 17885 17875 17873 17887
12536 5 1 1 0 1 17885 17875 17873 17887 17886 17876 17874 17888
12537 5 1 1 0 1 17645 17732 17730 17643 17883 17877 17875 17885
12538 5 1 1 0 1 17883 17877 17875 17885 17884 17878 17876 17886
12539 5 1 1 0 1 17647 17734 17732 17645 17881 17879 17877 17883
12540 5 1 1 0 1 17881 17879 17877 17883 17882 17880 17878 17884
12541 5 1 1 0 1 17633 17697 17734 17647 17832 17836 17879 17881
12542 5 1 1 0 1 17832 17836 17879 17881 17831 17835 17880 17882
12543 5 1 1 0 1 8383 8283 17349 17325 8410 8685 17559 17537
12544 5 1 1 0 1 8410 8685 17559 17537 8411 8207 17562 17540
12545 5 1 1 0 1 17325 17349 17348 17324 17537 17559 17558 17536
12546 5 1 1 0 1 17537 17559 17558 17536 17540 17562 17561 17539
4。输入计算结果文件名：该文件存放各种几何信息，包括材料分组信息
5。输入开挖信息文件名：该文件存放开挖分组信息。
所有材料的分组和开挖分组（group）都是自动定义的在ansys里将要分步开挖的单元分别找出来，按使用说明输入即可对不同开挖步信息分别给出它的group名
本code最大单元和结点数量可达1000000个。转换后的文件可能是大了点，不过很好用，转换后的文件直接call就行了。
将开挖信息列出来,是为了给它们一个group name,以便再开挖计算时使用。在ansys里列开挖单元是不难
的。
第一步，在ansys中建立三维模型，并给每个单元赋予相应的材料编号
第二步，将ansys中的单元信息和结点信息导出，格式如上述例子
第三步，在ansys中找出每一步开挖的单元号信息，格式如上
第四步，执行转换程序即可
材料分组和开挖单元分组由转换程序自动完成。
该程序绝对能用，如果出现问题，检查输入文件是否正确
在ＦＬＡＣ现有的模型，包括摩尔库伦和应力应变软化模型，运用到隧道等开挖支护的时候，并不能考虑到由于岩层土体应力变化对其材料参数的影响，很多情况下，没有锚杆支护与使用锚杆支护的效果并不明显，与实际出入较大．这是模型的缺陷．如果中间调整土体岩层材料的参数后，效果要明显许多，但参数的重新选取只能是经验来决定，误差较大．另外，锚杆支护时有个初锚力，然后再逐渐加载到一定值．模拟时，可以按此分步实现，但发现效果也并不明显．如果用ANSYS模拟锚杆支护，只能对杆端起到作用，不能准确模拟整个杆的受力状况．具体可以查阅中国矿业大学出版社＜锚杆支护的数值模拟＞．
一个煤矿巷道的例子：
new
;100m*100m
gr 120 100
;ROCK
m m j=1,30
m m j=32,43
m m j=45,56
m m j=57,60
m m j=62,65
m m j=67,76
m m j=78,79
m m j=81,84
m m j=86,100
;INTERFACE
m n j=31
m n j=44
m n j=61
m n j=66
m n j=77
m n j=80
m n j=85
;40m
gen 0.00, -32.50 0.0,30.00 40.00,17.00 40.00,-32.50 i=1,41 j=1,31
gen 0.00, 30.00 0.0,38.00 40.00,25.00 40.00,17.00 i=1,41 j=32,44
gen 0.00, 38.00 0.0,44.00 40.00,31.00 40.00,25.00 i=1,41 j=45,57
gen 0.00, 44.00 0.0,46.00 40.00,33.00 40.00,31.00 i=1,41 j=57,61
gen 0.00, 46.00 0.0,48.00 40.00,35.00 40.00,33.00 i=1,41 j=62,66
gen 0.00, 48.00 0.0,55.50 40.00,42.50 40.00,35.00 i=1,41 j=67,77
gen 0.00, 55.50 0.0,57.30 40.00,44.30 40.00,42.50 i=1,41 j=78,80
gen 0.00, 57.30 0.0,61.30 40.00,48.30 40.00,44.30 i=1,41 j=81,85
gen 0.00, 61.30 0.0,67.50 40.00,67.50 40.00,48.30 i=1,41 j=86,101
;20m
gen 40.00, -32.50 40.0,17.00 60.00,10.50 60.00,-32.50 i=41,81 j=1,31
gen 40.00, 17.00 40.0,25.00 60.00,18.50 60.00,10.50 i=41,81 j=32,44
gen 40.00, 25.00 40.0,27.50 47.90,24.7675 47.90,22.4325 i=41,57 j=45,50 gen 40.00, 27.50 40.0,31.00 47.90,28.4325 47.90,24.7675 i=41,57 j=50,57
gen 47.90,22.4325 47.90,24.7675 52.1,24.7675 52.1,21.0675 i=57,65 j=45,50 gen 47.90,24.7675 47.90,28.4325 52.1,27.0675 52.1,24.7675 i=57,65 j=50,57
gen 52.1,21.0675 52.1,24.7675 60,21 60,18.5 i=65,81 j=45,50
gen 52.1,24.7675 52.1,27.0675 60,24.5 60,21 i=65,81 j=50,57
gen 40.00, 31.00 40.0,33.00 60.00,26.50 60.00,24.50 i=41,81 j=57,61
gen 40.00, 33.00 40.0,35.00 60.00,28.50 60.00,26.50 i=41,81 j=62,66
gen 40.00, 35.00 40.0,42.50 60.00,36.00 60.00,28.50 i=41,81 j=67,77
gen 40.00, 42.50 40.0,44.30 60.00,37.80 60.00,36.00 i=41,81 j=78,80
gen 40.00, 44.30 40.0,48.30 60.00,41.80 60.00,37.80 i=41,81 j=81,85
gen 40.00, 48.30 40.0,67.50 60.00,67.50 60.00,41.80 i=41,81 j=86,101
;40m
gen 60.00, -32.50 60.0,10.50 100.00,-2.50 100.00,-32.50 i=81,121 j=1,31 gen 60.00, 10.50 60.0,18.50 100.00,5.50 100.00,-2.50 i=81,121 j=32,44 gen 60.00, 18.50 60.0,24.50 100.00,11.50 100.00,5.50 i=81,121 j=45,57 gen 60.00, 24.50 60.0,26.50 100.00,13.50 100.00,11.50 i=81,121 j=57,61 gen 60.00, 26.50 60.0,28.50 100.00,15.50 100.00,13.50 i=81,121 j=62,66 gen 60.00, 28.50 60.0,36.00 100.00,23.00 100.00,15.50 i=81,121 j=67,77 gen 60.00, 36.00 60.0,37.80 100.00,24.80 100.00,23.00 i=81,121 j=78,80 gen 60.00, 37.80 60.0,41.80 100.00,28.80 100.00,24.80 i=81,121 j=81,85 gen 60.00, 41.80 60.0,67.50 100.00,67.50 100.00,28.80 i=81,121 j=86,101
;INTERFACE
interface 1 aside from 1,31 to 121,31 bside from 1,32 to 121,32
interface 1 c= 0.2e+6 f=28 kn=588.00e+6 ks=14.70e+6
interface 2 aside from 1,44 to 121,44 bside from 1,45 to 121,45
interface 2 c= 0.2e+6 f=28 kn=588.00e+6 ks= 14.70e+6
;interface 3 aside from 1,57 to 121,57 bside from 1,58 to 121,58
;interface 3 c= 0.05e+6 f=20 kn=480.00e+6 ks= 9.00e+6
interface 4 aside from 1,61 to 121,61 bside from 1,62 to 121,62
interface 4 c= 0.1e+6 f=25 kn=588.00e+6 ks= 14.70e+6
interface 5 aside from 1,66 to 121,66 bside from 1,67 to 121,67 interface 5 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6 interface 6 aside from 1,77 to 121,77 bside from 1,78 to 121,78 interface 6 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6 interface 7 aside from 1,80 to 121,80 bside from 1,81 to 121,81 interface 7 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6 interface 8 aside from 1,85 to 121,85 bside from 1,86 to 121,86 interface 8 c= 0.1e+6 f=25 kn=588.00e+6 ks=14.70e+6
;BOTTOM
prop D=2500 C= 2.50e+6 F=24 S=3.6e+9 B=9.6e+9 j=32,43
;COAL
prop D=1400 C= 0.80e+6 F=18 S=8.06e+8 B=1.46e+9 j=45,56
;ROOF
prop D=2500 C= 3.00e+6 F=28 S=4.2e+9 B=13.1e+9 j=57,60 prop D=2600 C= 6.8e+6 F=35 S=13.2e+9 B=19.4e+9 j=62,65 prop D=2500 C= 2.50e+6 F=24 S=3.6e+9 B=9.6e+9 j=67,76 prop D=2600 C= 5.00e+6 F=30 S=7.8e+9 B=16.5e+9 j=78,79 prop D=2500 C= 3.50e+6 F=30 S=4.2e+9 B=13.1e+9 j=81,84 prop D=2550 C= 5.00e+6 F=32 S=7.8e+9 B=16.5e+9 j=86,100 set g=9.81
ini syy=-25.00e6 var 0,2.5e6
ini sxx=-30.00e6 var 0,3.0e6
ini szz=-30.00e6 var 0,3.0e6
app syy=-22.50e+6 j=101
fix x i=1
fix x i=121
fix y j=1
set large
his nstep=100
his unbal
pa
solve force 1e3
save pdsyy.sav
cal pdsold.txt
title
old support technology
ini yd=0
ini xd=0
ini yv=0
ini xv=0
ini x=47.9 i=57
ini x=52.1 i=65
his xd i=57 j=53
his xd i=65 j=53
his yd i=61 j=57
his yd i=61 j=50
M N i=57,64 j=50,56
INTERFACE 9 ASIDE FROM 57,50 TO 65,50 BSIDE FROM 65,50 TO 65,57 INTERFACE 10 ASIDE FROM 65,50 TO 65,57 BSIDE FROM 65,57 TO 57,57 INTERFACE 11 ASIDE FROM 65,57 TO 57,57 BSIDE FROM 57,57 TO 57,50 INTERFACE 12 ASIDE FROM 57,57 TO 57,50 BSIDE FROM 57,50 TO 65,50 interface 9 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 10 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 11 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
interface 12 c= 0.1e+6 f=20 kn=620.00e+6 ks=16.00e+6
;maogan
;left
;beam (tuo pan)
STR BEAM BEG 47.90,25.15 END 47.90,25.25 SEG=2 PROP=1
INTERFACE 21 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 1 TO NODE 3 INTERFACE 21 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 47.90,25.9 END 47.90,26.0 SEG=2 PROP=1
INTERFACE 22 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 4 TO NODE 6 INTERFACE 22 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 47.90,26.65 END 47.90,26.75 SEG=2 PROP=1
INTERFACE 23 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 7 TO NODE 9 INTERFACE 23 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 47.90,27.4 END 47.90,27.5 SEG=2 PROP=1
INTERFACE 24 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 10 TO NODE 12 INTERFACE 24 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 47.90,28.15 END 47.90,28.25 SEG=2 PROP=1
INTERFACE 25 ASIDE FROM 57,50 to 57,57 BSIDE FROM NODE 13 TO NODE 15 INTERFACE 25 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
;cable
STR CAB BEG node 2 END 47.8034,25.1741 SEG=1 PROP=2
STR CAB BEG NODE 16 END 46.9341,24.9412 SEG=4 PROP=3
STR CAB BEG NODE 20 END 45.9681,24.6824 SEG=5 PROP=4 ten 10000
STR CAB BEG node 5 END 47.80,25.95 SEG=1 PROP=2
STR CAB BEG NODE 26 END 46.90,25.95 SEG=4 PROP=3
STR CAB BEG NODE 30 END 45.90,25.95 SEG=5 PROP=4 ten 10000
STR CAB BEG node 8 END 47.80,26.7 SEG=1 PROP=2
STR CAB BEG NODE 36 END 46.90,26.7 SEG=4 PROP=3
STR CAB BEG NODE 40 END 45.90,26.7 SEG=5 PROP=4 ten 10000
STR CAB BEG node 11 END 47.80,27.45 SEG=1 PROP=2
STR CAB BEG NODE 46 END 46.90,27.45 SEG=4 PROP=3
STR CAB BEG NODE 50 END 45.90,27.45 SEG=5 PROP=4 ten 10000
STR CAB BEG node 14 END 47.8034,28.2259 SEG=1 PROP=2
STR CAB BEG NODE 56 END 46.9341,28.4588 SEG=4 PROP=3
STR CAB BEG NODE 60 END 45.9681,28.7176 SEG=5 PROP=4 ten 10000
;ROOF
;(gangjinti)
STR BEAM BEG 47.90,28.4325 END 52.10,27.0675 SEG=10 PROP=5 INTERFACE 26 ASIDE FROM 57,57 to 65,57 BSIDE FROM NODE 66 TO NODE 76 INTERFACE 26 c=2e6 f=22 ks=1e7 kn=3e8 glued
;beam (tuo pan)
STR BEAM BEG 48.00,28.4 END 48.10,28.3675 SEG=2 PROP=1
INTERFACE 27 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 77 TO NODE 79 INTERFACE 27 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 48.65,28.1887 END 48.75,28.1563 SEG=2 PROP=1
INTERFACE 28 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 80 TO NODE 82 INTERFACE 28 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 49.30,27.9775 END 49.40,27.9450 SEG=2 PROP=1
INTERFACE 29 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 83 TO NODE 85 INTERFACE 29 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 49.95,27.7662 END 50.05,27.7337 SEG=2 PROP=1
INTERFACE 30 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 86 TO NODE 88 INTERFACE 30 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 50.60,27.5550 END 50.70,27.5225 SEG=2 PROP=1
INTERFACE 31 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 89 TO NODE 91 INTERFACE 31 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 51.25,27.3437 END 51.35,27.3113 SEG=2 PROP=1
INTERFACE 32 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 92 TO NODE 94 INTERFACE 32 c=2e6 f=22 ks=1e7 kn=3e8 glued
STR BEAM BEG 51.90,27.1325 END 52.00,27.10 SEG=2 PROP=1
INTERFACE 33 ASIDE FROM node 66 to node 76 BSIDE FROM NODE 95 TO NODE 97 INTERFACE 33 c=2e6 f=22 ks=1e7 kn=3e8 glued
;cable
STR CAB BEG node 78 END 48.05,28.4889 SEG=1 PROP=6
STR CAB BEG NODE 98 END 48.05,29.0889 SEG=2 PROP=7
STR CAB BEG NODE 100 END 48.05,30.5889 SEG=7 PROP=8 ten 20000
STR CAB BEG node 81 END 48.7309,28.2677 SEG=1 PROP=6
STR CAB BEG NODE 108 END 48.9163,28.8383 SEG=2 PROP=7
STR CAB BEG NoDE 110 end 49.3798,30.2648 SEG=7 PROP=8 ten 20000
STR CAB BEG node 84 END 49.3809,28.0565 SEG=1 PROP=6
STR CAB BEG NODE 118 end 49.5663,28.6271 SEG=2 PROP=7
STR CAB BEG NoDE 120 END 50.0298,30.0536 SEG=7 PROP=8 ten 20000
STR CAB BEG node 87 END 50.0310,27.8453 SEG=1 PROP=6
STR CAB BEG NODE 128 END 50.2164,28.4159 SEG=2 PROP=7
STR CAB BEG NoDE 130 END 50.6798,29.8423 SEG=7 PROP=8 ten 20000
STR CAB BEG node 90 END 50.6810,27.6341 SEG=1 PROP=6
STR CAB BEG NODE 138 END 50.8664,28.2047 SEG=2 PROP=7
STR CAB BEG NoDE 140 end 51.3298,29.6311 SEG=7 PROP=8 ten 20000
STR CAB BEG node 93 END 51.3310,27.4229 SEG=1 PROP=6
STR CAB BEG NODE 148 end 51.5164,27.9935 SEG=2 PROP=7
STR CAB BEG NoDE 150 END 51.9798,29.4198 SEG=7 PROP=8 ten 20000
STR CAB BEG node 96 END 52.0066,27.20333 SEG=1 PROP=6
STR CAB BEG NODE 158 END 52.3449,27.7243 SEG=2 PROP=7
STR CAB BEG NoDE 160 END 53.1482,28.9613 SEG=7 PROP=8 ten 20000
;Right
;beam (tuo pan)
STR BEAM BEG 52.10,25.3175 END 52.10,25.4175 SEG=2 PROP=1
INTERFACE 34 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 168 TO NODE 170 INTERFACE 34 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 52.10,26.0675 END 52.10,26.1675 SEG=2 PROP=1 INTERFACE 35 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 171 TO NODE 173 INTERFACE 35 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
STR BEAM BEG 52.10,26.8175 END 52.10,26.9175 SEG=2 PROP=1 INTERFACE 36 ASIDE FROM 65,50 to 65,57 BSIDE FROM NODE 174 TO NODE 176 INTERFACE 36 c=0.7e6 f=12 ks=1e6 kn=3e7 glued
;cable
STR CAB BEG node 169 END 52.20,25.3407 SEG=1 PROP=2
STR CAB BEG NODE 177 END 53.10,25.0996 SEG=4 PROP=3
STR CAB BEG NODE 181 END 54.0319,24.8499 SEG=5 PROP=4 ten 10000
STR CAB BEG node 172 END 52.20,26.1175 SEG=1 PROP=2
STR CAB BEG NODE 187 END 53.10,26.1175 SEG=4 PROP=3
STR CAB BEG NoDE 191 end 54.10,26.1175 SEG=5 PROP=4 ten 10000
STR CAB BEG node 175 END 52.20,26.8943 SEG=1 PROP=2
STR CAB BEG NODE 197 end 53.10,27.1354 SEG=4 PROP=3
STR CAB BEG NoDE 201 END 54.0319,27.3851 SEG=3 PROP=4 ten 10000
;tuo pan prop=1 (steel size:100*100*10mm)
STR PROP=1 E=2.0E11 I=1.25e-9 area=5.5e-4 d=7800.00 wid=0.125
;cable prop (20mm*2m)
STR PROP=2 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=2 d=7800.0 kb=10e10 ks=9e7 yield=0.60e5 ycom=3e8
STR PROP=3 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=3 d=7800.0 kb=6e2 ks=7e3 yield=0.6e5 ycom=3e8
STR PROP=4 E=2.0E11 i=1.1e-9 area=3.1e-4
str prop=4 d=7800.0 kb=10.0e9 ks=9.0e6 yield=0.6e5 ycom=3e8
;gangjinti prop
STR PROP=5 E=2.0E11 i=5.2e-9 area=2.0e-4 d=7800.00 wid=0.25
STR PROP=6 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=6 d=7800.0 kb=10e10 ks=9e7 yield=0.60e5 ycom=3e8
STR PROP=7 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=7 d=7800.0 kb=6e2 ks=7e3 yield=0.6e5 ycom=3e8
STR PROP=8 E=2.0E11 i=2.49e-9 area=4.9e-4
str prop=8 d=7800.0 kb=10.0e9 ks=9.0e6 yield=0.6e5 ycom=3e8
step 11756
save pdsold.sav
地板回弹的例子：
restore project.sav
;model elastic
model moh
;prop bulk =8.333e9 shear = 0.3e8 fric = 35 coh=1e3 tens=1e3 prop bulk =8.333e9 shear = 1.1111e10 fric = 35 coh=1e6 tens=1e6 set grav 0,-9.81, 0
ini dens = 3000
ini sxx -10.51e6 gradient 0,1.6e4,0 range y -30,120
ini syy -8.76e6 gradient 0,2.6e4,0 range y -30,120
ini szz -12.71e6 gradient 0,2.2e4,0 range y -30,120
fix x range x -100.1 -99.9
fix x range x 199.9 200.1
fix y range y -150.1 -149.9
fix y range y 149.9 150.1
fix z range z -0.1 0.1
fix z range z -149.9 -150.1
hist unbal
;hist gp zdisp 500,400,200
;set mech force = 500
plot hist 1
solve
save s2.sav
;restore s2.sav
ini xdis 0 ydis 0 zdis 0
model null range group 1
set large
solve
save s3.sav
;
;restore s3.sav
model null range group 2
hist unbal
solve
save s4.sav
;
;restore s4.sav
model null range group 3
solve
save s5.sav
;
;restore s5.sav
model null range group 4
solve
save s6.sav
;
;restore s6.sav
model null range group 5
solve
save s7.sav
;
;restore s7.sav
model null range group 6
solve
save s8.sav
;
;restore s6.sav
model null range group 7
solve
save s9.sav
;
;restore s2.sav
model null range group 8
solve
save s10.sav
;
;restore s6.sav
model null range group 9
solve
save s11.sav
;
;restore s6.sav
model null range group 10
solve
save s12.sav
应该说现在岩土工程专业处于一种蓬勃发展却难于有突破的阶段。作为一个岩土工程专业的研究生面临着一种困境，力学和数学知识越来越需要，但是任何一门不管是力学或者是数学都可以耗掉你的一生。尽管很多现代的数学和力学的知识用到了岩土工程中并产生了作用，但都难以有实质性的突破。记得爱因斯坦曾就物理和数学的关系发表过看法，原话记不请了，大意是他曾下决心掌握最先进的数学理论，但发现任何一个分支数学都能用尽他的一生，这使得他很苦恼，后来他发现他的错误在于对于物理的本质认识不足以清楚到能挑选出合适的数学工具来解决问题。
可能我们岩土人面临的苦恼要比爱因斯坦曾经经历的还要痛苦。
力学、数学都属于基础课程，目前很多岩土界的发展都是将这两方面的技术应用到岩土问题上来，而且有一些直接从力学转行过来的专家，比如石根华。
再说数学，这个东西无论与哪种学科交织一起都可能会产生巨大的成果，比如很多的经济学家。
从我们个人来说，数学、力学都是很重要的一门基础，要是这两样工具都能灵活掌握的话，那你岩土方面的论文做出来肯定档次比较高，说的通俗一点，就是别人很难懂。但是要知道，对于岩土，也有其基础学科，在我们需求着高深的数学力学知识的同时，有没有考虑到自己的土力学、岩石力学、试验土力学、理论土力学、计算……这些科目有没有过关呢？
作为一门专业，尤其是这个行业的工作当中，需要的比较丰富的专业知识和经验，所以在学校中我们也应当从这两方面入手。在本科的时候，曾经有过学计算机的热潮，但到最后真正从是这个行业的又有多少呢？试想，我们不会现在正在掀起一场力学、数学的高潮吧，这同样和计算机一样是消耗一生的学问啊？
其实可以这么说，对于数学和力学的掌握程序对于不同的人有不同的要求，如果是作理论研究，那我想什么都不用说了，如果没有扎实的数学和力学基础，我想你会感到力不从心的。如果是从事工程方面的具体工作，那就另当别论了，但是也并不是说数学和力学不重要了，而是说没有必要搞那么深，基本的一些力学基础还得有，只有有了这些基础的东西，专业的东西才能吃透。才可能会在工程界有所创新。
我个人认为，不管你是从事理论研究，还是从事工程的具体工作，数学和力学一定要学好，这个学好的标准不同。这就像不同的工匠需要的工具不同一样，正所谓杀鸡焉用牛刀。譬如说，对于一个从事本构研究的人来说，如果不懂张量理论，那他的研究工作将是一种艰苦的工作。但对于工程人员来说，掌握不掌握晦涩难懂的张量理论来说对他关系不是很大，他可以用别人研究出来的现成的本构的算法编程解决问题，或用现在的商业软件解决实际问题，对于实际问题的解释更多的是靠专业及基本的理论知识而不是张量理论。所以我认为，掌握到何种程度要看你要用这些工具（数学、力学）到何种程度，这并没有一个确定的标准，应该跟你从事的工作相关。
另我一个需要强调的就是，专业知识是我们立足的基础，要不然，任何力学或数学出身的人就可以代替我们了，数学和力学只是我们的一种工具，就像我们天天灌水离不开PC一样，但我们用PC灌水而不必了解PC的详细构造及设计原理。
有一次导师跟我们交流时，说了这么一句话，说好多人认为专业课不太难，我们的专业课非本专业的人好像也可以看懂，但一细究，好像他们也会是什么也不懂，这是因为我们的专业涉及到的许多东西是跟工程紧密结合在一些的，没有感性认识的人理解起来还是有一定困难的，这也是为什么在我们这个行当，越老越吃香的原因吧！这也是从一个侧面说明了专业知识的重要性。
岩土工程是一个独立的专业，但它更象一个边缘学科，比如岩土耦合研究（渗流）。具有扎实的数理基础会具有优势，但是岩土工程研究的终极目的是解决岩土工程中的问题，使岩土工程更可靠，更经济。因此，真正的岩土专业背景才是根本。数学、物理只是工具而已，对于工具，需要专门的人去研究。岩土工程师需要提出问题，涉及到的数理问题请数理学家参与解决似乎更合理一点。否则，正象楼上所说，耗尽一生去研究方法，岂不是本末倒置么！总之一句话，岩土工程要有所突破，需要多学科的学者来合作。（石根华博士是特例，他本人数学出身，而且有丰富的岩土工程经验，所以有丰硕的成果。受人尊敬。）重要的要有扎实的基础理论。特别是在做研究生的，一定不要图省事那个商业软件来做论文。殊不知，在方便背后，你把应该学习掌握的对于基础理论的训练过程省略了，这是目前令人遗憾的事情。在这一点上，我觉得河海大学要求学生编写程序的传统比较好，虽然可能最后编写的程序比较粗糙，比不上商业软件漂亮，可是编程序的过程所做的都是岩土或数值计算的核心内容，这才是重要的。
我在读书阶段为编写程序付出了很多辛苦，但我觉得挺好。也非常感谢导师对我的严格要求，若非如此，我没有现在的收获。
搞岩土的，我认为从力学角度看应该是用数值方法来求解边界值问题，因为，对于岩土这么复杂的介质几乎不存在解析解。
岩土很多都是非连续问题，成熟的有限元又无能为力，所以只能学习DDA（离散元）,EFGM（无单元）,数值流形，界面元等，而这些离散方法作为初学者很难接受，所以还要学习成熟的有限元，在学习这些不连续问题的离散方法。想学会有限元，想学会有限元，又要有足够的力学功底：弹性力学（铁摩辛柯），弹塑性理论，接触力学，断裂力学，弹性稳定理论。
从数学角度来说应该掌握：数值分析，数值逼近，微分方程的数值解法，范函，变分原理，矩阵论。这都是计算力学的基础
就拿一般的岩土软件计算来说，如果你有较为深厚的岩土专业背景，可能你对计算的一些非常具体的内部细节缺乏了解，但你仍能够大致判断出计算结果的可靠性。但是作为研究人员，要深入了解现有计算模型跟实际工程情况之间的差距，那么就必须深入学习了解众多现有的各种复杂的岩土力学模型(这需要很多数学和力学知识如楼上总结)的细节和出发点，同时也要对这些模型在工程实际中应用存在的模型误差，这种误差可能是假设不当，也可能是理论不甚严密。要想把这些搞清楚绝非一年之功。而在此扎实基础上才能够想办法对原有模型提出改进，或对计算方法提出改进。沈珠江院士曾讲过，创新并不是简单地把别的学科里新鲜的概念原理搬过来，而是要把这些新东西跟岩土学科交融起来。
比如石根华是把拓扑学、谢和平把分形科学引入到岩土工程中，都为岩土工程学科发展作出了大的贡献。
所以说，岩土工程的研究生如果想在理论上有所创新，是颇有挑战性的。这也许是这个岩土工程论坛有足够的繁荣空间的理由。
建模的思路：
1.在ansys 中建立方块，
2.画圆圈，画2横一竖
3.通过线分割方块，生成众多平面
4.画轴线，然后沿着线扫掠生成题，
5.然后划分网格，注意设计不同部位体的材料不同
6.导入flac
7.设置边界条件
8.开挖，设置初期支护，并连接初期支护
9.求解，拷贝结果图片
划分网格的时候，用到面连接，里面的用到了，旋转一个面生成划分网格。
隧道命令流解析
;产生一个体，这个体里面包含隧道
gen zone radbrick &
p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &
size 3,5,5,7 &
ratio 1,1,1,1.5 &
dim 1 4 2 fill
;设置模型服从摩尔库仑屈服准则
mod moh
;定义bulk模量，剪切模量，摩擦角，粘聚力，和tens
prop bulk 1e8 shear 3e8 fric 35 coh 1e3 tens 1e3
;施加约束
fix x range x -0.1 0.1
fix z range z -0.1 0.1
fix y range y 9.9 10.1
;施加重力
set grav 0 10 0
;设置密度
ini dens 1000
;设置初始应力，这个有什么用呢？？？？
ini syy 0.0 grad 0 -10000 0
ini sxx 0.0 grad 0 -5000 0
ini szz 0.0 grad 0 -5000 0
apply sxx 0.0 grad 0 -5000 0 range x 9.9 10.1
apply
szz 0.0 grad 0 -5000 0 range z 9.9 10.1
;挖坑的命令，使用fish语言
def digHole
command
range name trench x 0 1 y 0 4 z 0 2
;挖了
model nul range trench
set large
hist gp xdisp 1 0 0
endcommand
end
;表面蓝色显示
plot surf lblue
;休息
pause
;开挖
digHole
;设置大变形
set large
;这个好像是施加位移
hist gp xdisp 1,0,0
; ... and then calculate the response
;设置计算步为1000
step 1000
;我不知道上述解释是否明白，是否对头？？
在手册中的图中有显示shell单元
但我自己做的时候不知用何命令才能在网格图中将的shell，beam等单元显示出来，高手指点！！
plot add sel geom
另外，还可在上面命令的后面加：cid=on node=off .....
再问一下，shell是否只可以加在网格的表面？怎么在网格内部加不上？
可以加上了，但是不知道厚度t如何表现。
plotitems-add-structural element-geometry-shell
就ok了
save 1.dat
solve
restore 1.dat
边界条件可以改
apply施加边界条件，initial 施加初始条件。
什么是边界、初始条件要查有关力学的专业基础书。
很多学生觉得现有各种程序难学，其实是自己专业基础知识没学好。
学好连续介质力学那一套理论，你就会发现，不是程序与英文难懂，而是专业理论没学好!
可以通过编写fish函数来模拟断层的。
THEORY AND BACKGROUND
Section 3 : Interfaces
ftd132.pdf
里面有例子
flac怎么作一个剖面！就是我建了一个模型，怎么看一个切面的视图。
pl set plane ori 0 0 0 nor 1 0 0
pl con disp plane
本在遇到一些渐进性破坏的问题，而缺少这方面的早期的资料。好象关于岩土方面的渐进怀破坏的是由Skempton引进的，所以想寻求Skempton、Bishop等早期牛人关于这方面的论文。望在这方面有兴趣与我交流!!!
请问各位大侠，grout stiffness 是什么意思？
1.gr_k grout stiffness per unit length, kg [F/L2]
2.Ft = kg* ut
3.The grout annulus is assumed to behave as an elastic-perfectly plastic solid. As a result of
relative shear displacement,
ut , between the tendon surface and the borehole surface, the shear force, Ft , mobilized
per length of cable is related to the grout stiffness, kg — i.e.,
4.可见k_g是grout annulus与zone（soil）之间的剪切刚度。（单位剪切位移引起的应力）
5.详细参考cable 一节，有图形表示，很清楚。
我是模仿手册Example Applications部分中Excavation in a Saturated Soi章的例子来建立连续墙模型的，它是和soil仪器设置zone，然后将连续墙部分zone的density，ini为连续墙密度，不知道chichilv 你是怎么做的，有好的方法么？
我现在头疼的是连续墙的支撑，那位有这方面的例子么？
有一个例子，它是按1e-5m/step来设置Y方向速度的，用FISH写的，可以参考。
new
gen zone brick size 1 1 1
title
Oedometer test on Mohr-Coulomb sample
m m
pro bu 200 sh 200 co 1 fric 10 dil 10 ten 5.67
;pro bu 200 sh 200 co 1 fric 10 dil 0 ten 5.67
def d_sigy
a_sy=0.0
pnt =zone_head
c_k=z_prop(pnt,'bulk')
c_g=z_prop(pnt,'shear')
e1=c_k+4.0*c_g/3
e2=c_k-2.0*c_g/3
sf=z_prop(pnt,'friction')*degrad
nf=sin(sf)
nf=(1+nf)/(1-nf)
sp=z_prop(pnt,'dilation')*degrad
np=sin(sp)
np=(1+np)/(1-np)
r1=(e1-e2*nf)/((e1+e2)*nf*np-2*e2*(nf+np)+2*e1)
vyv=-1e-5
dsigy=vyv*(e1+2*r1*(e2*np-e1))
step1=-2*z_prop(pnt,'cohesion')*sqrt(nf)/((e1-e2*nf)*vyv)
end
def esigy
while_stepping
if step a_sy=a_sy+e1*vyv else a_sy=a_sy+dsigy end_if n_sy=z_syy(pnt) end fix x y z d_sigy ini yvel vyv range y 0.9 1.1 his nstep 50 his gp ydisp 0 1 0 his n_sy his a_sy step 1000 plot his -2 line -3 mark vs -1 xmin 0 ymin 0 ymax 4.5 gen zone brick size 10 3 5 mod mohr prop bulk 1e8 shear .3e8 fric 35 prop coh 1e3 tens 1e3 ini dens 1000 set grav 0,0,-10 fix x y z range z -.1 .1 fix y range y -.1 .1 fix y range y 2.9 3.1 fix x range x -.1 .1 fix x range x 9.9 10.1 set large hist unbal solve save cable.sav ini xdis 0 ydis 0 zdis 0 hist gp xdisp 0,1,5 def place cables loop n (1,5) z_d = 5.5 - float z_t = z_d + 0.5 z_b = z_d - 0.5 command free x range x -.1,.1 z z b,z t solve sel cable beg 0.0,0.5,z d end 7.0,0.5,z d nseg 7 sel cable beg 0.0,1.5,z d end 7.0,1.5,z d nseg 7 sel cable beg 0.0,2.5,z d end 7.0,2.5,z d nseg 7 sel cable prop emod 2e10 ytension 1e8 xcarea 1.0 & gr k 2e10 gr coh 1e10 gr per 1.0 endcommand endloop end 是把锚杆分为多个单元吗。如何分段模拟呢。 用结构单元cable ，当然可以分段，这个手册中有很具体的解释，祥看结构单元一章 sel cable id 1 beg 1 6 39.0 end 1 26 37.0 nseg 20 sel cable id 1 prop emod 2.06e11 ytens 3.15e8 xcarea 6.2e-4 & gr_k 4.2e9 gr_coh 2.5e4 gr_fric 10 gr_per 0.471 range y 16 26 我在做－problem sloving 的中3-94页的例子时怎么不能显示例子中所列出的图呀？ pl set color off pl fos pl add con ssr out on pl add velo FLAC3D 中的 solve fos 可以计算岩质边坡的稳定性吗？ 其实，如果设置了薄层（即已经知道滑动面），可将薄层上下的材料参数放大，上诉的第4种方法就变成第2种了，强度折减法的使用有一个要点，不会发生剪切破坏的区域（如基岩）要将其材料参数放大，使材料参数的折减不影响边坡的分析！ 二 抗剪强度折减系数法的理论 2．1抗剪强度折减系数法的概念 抗剪强度折减系数(SSRF ：Shear Strength Reduction Factor)定义为：在外荷载保持不变的情况下，边坡内土体所发挥的最大抗剪强度与外荷载在边坡内所产生的实际剪应力之比。这里定义的抗剪强度折减系数，与极限平衡分析中所定义的土坡稳定安全系数在本质上是一致的。 2．2抗剪强度折减系数法的具体内容 所谓抗剪强度折减技术就是将土体的抗剪强度指标C 和φ，用一个折减系数s F ，如式 (1)和(2) 所示的形式进行折减，然后用折减后的虚拟抗剪强度指标F C 和F φ，取代原来的抗剪强度指标C 和φ，如式（3）所示。 s F F C C /= （式1） )/)((tan tan 1s F F φφ-= （式 2） F F fF C φστtan += （式3） 式中：F C 是折减后土体虚拟的粘聚力；F φ是折减后土体虚拟的内摩擦角；fF τ是折减后的抗剪强度。 折减系数s F 的初始值取得足够小，以保证开始时是一个近乎弹性的问题。然后不断增加s F 的值，折减后的抗剪强度指标逐步减小，直到某一个折减抗剪强度下整个土坡发生失稳，那么在发生整体失稳之前的那个折减系数值，即土体的实际抗剪强度指标与发生虚拟破坏时折减强度指标的比值，就是这个土坡的稳定安全系数。 2．3抗剪强度折减系数法的优点 结合有限差分法的抗剪强度折减系数法较传统的方法具有如下优点： （1）能够对具有复杂地貌、地质的边坡进行计算； （2）考虑了土体的本构关系，以及变形对应力的影响； （3）能够模拟土坡的边坡过程及其滑移面形状（通常由剪应变增量或者位移增量确定滑移面的形状和位置）； （4）能够模拟土体与支护结构（超前支护、土钉、面层等）的共同作用； （5）求解安全系数时，可以不需要假定滑移面的形状，也无需进行条分。 后处理都做了哪些项目，请大家讨论！ 支护\衬砌的内力,观察位移,等值曲线 对于土体的: 还应该记录一些单元的应力变化,用这些应力值结合土工计算原理,做一些曲线，求一些参数,如最大主应力和最小主应力差和应变的关系等,画几个应力圆,,判断塑性硬化和塑性软化.做相关的评价,博士对于这些应该举手之劳,对于我们简直是天文! 画曲线可以借助tecplot 等软件,看看变形曲线是不是符合规范. 对于支护结构的: 结果后处理,还有就是应该把弯矩和轴力等内力提取出来，用我们习惯的形式表达出来.这个是材料力学\结构力学\钢筋混凝土原理的相关知识,也要做安全性的评价,评价支护结构的安全性,看看支护结构做的够不够. 强度折减法，邓颖人教授有一些文章，宋二祥教授也有一些，flac 在这点上比abaqus 好，最起码flac 提供了一个计算的算法！ １．hist 中换线形如何来换啊，比如我想把实线用虚线来代替，在哪里设置啊！ ２．想只显示开挖部分用什么命令啊！ ３．手册中３－５３中的这个图－shell principle stress-1 这个图用什么命令生成的，请指教啊！1.好像不能换线型，用颜色区分 2.比如说 pl bl gr null on ，plot item 对话框中有选项，点上就行 3.哪张图？ ADAMS中的函数_New ADAMS中的函数 ADAMS/View中系统提供的数 学函数大致分类介绍如下。（1）基本数学函数 ABS(x) 数字表达式x的绝对值 DIM(x1，x2) x1>x2时x1与x2之间的差值，x1<x2时返回0 EXP(x) 数字表达式x的指数值 LOG(x) 数字表达式x的自然对数值LOG10(x) 数字表达式x的以10为底的对数值 MAG(x，y，z) 向量[x，y，z]求模 MOD(x1，x2) 数字表达式x1对另一个数字表达式x2取余数 RAND(x) 返回0到1之间的随机数SIGN(x1，x2) 符号函数，当x2>0时返回ABS(x)，当x2<0时返回－ABS(x) SQRT(x) 数字表达式x的平方根值（2）三角函数 SIN(x) 数字表达式x的正弦值SINH(x) 数字表达式x的双曲正弦值COS(x) 数字表达式x的余弦值COSH(x) 数字表达式x的双曲余弦值TAN(x) 数字表达式x的正切值TANH(x) 数字表达式x的双曲正切值ASIN(x) 数字表达式x的反正弦值ACOS(x) 数字表达式x的反余弦值 ATAN(x) 数字表达式x的反正切值ATAN2(x1，x2) 两个数字表达式x1，x2的四象限反正切值 （3）取整函数 INT(x) 数字表达式x取整 AINT(x) 数字表达式x向绝对值小的方向取整 ANINT(x) 数字表达式x向绝对值大的方向取整 CEIL(x) 数字表达式x向正无穷的方向取整 FLOOR(x) 数字表达式x向负无穷的方向取整 NINT(x) 最接近数字表达式x的整数值RTOI(x) 返回数字表达式x的整数部分 位置/方向函数位置/方向函数用于根据不同输入变量计算有关位置或方向的参数。ADAMS/View中系统提供的位置/方向函数分类介绍如下。（1）位置函数 LOC_ALONG_LINE 返回两点连线上与第一点距离为指定值的点 LOC_CYLINDRICAL 将圆柱坐标系下坐标值转化为笛卡儿坐标系下坐标值LOC_FRAME_MIRROR 返回指定点关于指定坐标系下平面的对称点 LOC_GLOBAL 返回参考坐标系下的点在全局坐标系下的坐标值 贴片天线研究 第一部分天线的基本知识 (2) 第二部分贴片天线设计 (11) 第三部分贴片天线的应用 (24) 第四部分贴片天线的性能 以及SAR的分布 (31) 附录 (38) 小组成员：李黎轩冷继男 钟颐华刘同 2004年1月2日 第一部分 天线的基本知识 总括 天线是我们在设计射频系统时所需考虑得最后一部分内容。然而可不能小视天线的重要作用，轻敌将导致设计前功尽弃。天线作为无线传输的一部分，它的作用概括起来说是传送与接受电磁场能量。在第一部分中，我们将介绍天线的最基本知识，以指导接下来贴片天线的设计。 定义 天线是一个具备传输与发送电磁能量的导电元件。天线能够将电磁能量转化为电磁场传播出去，同时又能够通过将空间中的电磁场转化为电磁能量来接收电磁波。如何在同一天线上实现电磁能量的接收（receive ）与传播(transmit)是天线的一个重要属性 . 天线的主要特征参数有： 天线的中心频率（center frequency ）、带宽(bandwidth)、天线的极化(polarization)、天线增益gain 、辐射模型(radiation pattern)、阻抗(impedance)。 传输线的特征参数 λ Lambda Wavelength （单位：米） 在自由空间中传播的电磁场，速度为光速。即8 3.0010/c m s =?. VSWR Voltage Standing Wave Ratio ，电压驻波系数 dB Decibel 分贝的引入为在使用中表示方便 dBm dBm 表示功率，相对于1 mw 为基准定义 dBi 天线增益，以等方向天线为参考 1.step可能是最常用的: step(time,0,0,1,50)+ step(time,4,0,6,-100)+ step(tme,9,0,10,50) 函数原形STEP(A,x1,h1,x2,h2) 解释：由数组A的x值，生成区间（x1，h1）至（x2，h2）之间的阶梯曲线，返回y值的数据。 举个常用的例子。 比如STEP(time,1,0,2,100) time在adams中是个递增的变量，相当于一个数组。那么step的返回值就是随着time变化的值。 这个例子将表示在time从（1，2）的过程中，返回值将从0，100。看看例子，两个小球，一个使用step 函数设置了位移，另外一个是参考。当然，这个变化过程，adams使用了缓和的图形，从其位移图中可以看出来。step既然是个返回值，就可以使用加减法了。如上例，如果设置下面的小球的位移如下：STEP(time,1,0,2,100)+step(time,2,0,3,400)+step(time,3,0,4,-200) 2.以前用过碰撞函数，有单向和双向函数的区分，其中系统的球面等碰撞为其特例！ IMPACT (Displacement Variable, Velocity Variable, Trigger for Displacement Variable, Stiffness Coefficient, Stiffness Force Exponent, Damping Coefficient, Damping Ramp-up Distance) BISTOP (Displacement Variable, Velocity Variable, Low Trigger for Displacement Variable, High Trigger for Displacement Variable, Stiffness Coefficient, Stiffness Force Exponent, Damping Coefficient, Damping Ramp-up Distance) 3.if函数 这个函数最好不要使用，他的使用会带来突变，会使运算的时候不收敛。不过应急的时候还是可以一用。 if(time-1:1,0,if(time-2:0,-1,-1)) IF(Expression1: Expression2, Expression3, Expression4) adams要计算Expression1的值： 如果他的值小于0，则执行Expression2语句，如果Expression1的值等于0，则执行Expression3语句，如果Expression1的值大于0，则执行Expression4语句 我得if语句的意思是：如果时间小于1的时候，加速度为1，如果时间为1，加速度为0，如果时间大于1小于2，则加速度为0，如果时间大于、等于2则，加速度为-1 4. 我得一个想法 就是利用sign函数构造 比较常用的是给机构加上一个与运动方向相反的作用力等等可以先测量施加力对象的运动速度，然后利用速度的变化，插入measure到sign函数里面就可以获得与运动方向相反的作用力 现代远程教育 《药学文献检索（本科）》 课 程 学 习 指 导 书 作者：刘伟 07年 6 月 第一章绪论 （一）本章学习目标 掌握：文献和药学文献的定义；文献的类型；文献检索的定义和原理；检索语言；检索的方法；检索途径。 熟悉：文献检索工具；检索步骤。 了解：文献的特点；重要性。 （二）本章重点、要点 文献的类型；检索语言；检索的方法；检索途径。 （三）本章练习题 1．名词解释：文献；药学文献；ISBN；核心期刊；ISSN；标准；药学文献检索；检索语言；分类检索语言；索引。 2．简答题： （1）依据检索对象的不同，简述文献检索可分为哪几种类型？ （2）简述文献检索的原理？ （3）简述检索语言在信息检索中起着什么重要的作用？ （4）简述检索工具的特点？ （5）简述题录和文摘型检索工具的联系及区别？ （6）简述索取原始文献的方法？ 3．论述题： （1）请论述文献按加工层次如何划分，及其特点和联系。 （2）论述常用的检索方法有哪些，各有什么优缺点？ （3）论述从文献的内容特征进行文献检索的途径有哪些？各有什么优缺点。 （一）本章学习目标 掌握：《中国药学文摘》、《中文科技资料目录·中草药》、《中文科技资料目录·医药卫生》的编排结构、检索途径与方法； 熟悉：《中国药学文摘》、《中文科技资料目录·中草药》、《中文科技资料目录·医药卫生》的概况；《中文科技资料目录·中草药》的使用技巧；《国外科技资料目录·医药卫生》的概况、编排结构和检索途径及方法；《中国医学文摘》的概况、编排结构和检索途径。 了解：《中国药学文摘》、《中文科技资料目录·中草药》的分类目次。 （二）本章重点、要点 《中国药学文摘》、《中文科技资料目录·中草药》、《中文科技资料目录·医药卫生》的检索途径与方法； （三）本章练习题 1．简答题： （1）简述《中国药学文摘》的编排结构？ （2）简述《中国药学文摘》期和年主题索引的区别？ （3）简述《中文科技资料目录·医药卫生》分类途径的检索方法？ 2．论述题： （１）请说明查找“阿司匹林的不良反应”的相关文献时，应首选什么的手工检索工具，怎样来进行查找？ （２）请说明查找“西洋参黄酮药理作用”的相关文献时，应选择什么的检索工具，怎样来进行查找？ （３）请解释下面著录格式每一部分的含义 ①R284.1 ②031 ③200309444 ④用细胞膜色谱法筛选研究红毛七中的有效成分/⑤高琨（西安交通大学药学院）；贺浪冲，杨广德//⑥中国药学杂志-⑦2003，⑧38（1）-14—16 ⑨采用兔血管细胞膜色谱（CMC）模型筛选…（略） （４）请解释下面著录格式每一部分的含义 ①996550 ②天然甾体皂甙化合物的抗肿瘤活性/③李忌（华西医科大学基础医学院），陈俊杰，巨勇…//④天然产物研究与开发.- ⑤1999,11(1).-14-17 电磁辐射的测量基础知识 电磁辐射的测量方法通常与测量点位和辐射源的距离有关，即，所进行的测量是远场测量还是近场测量。由于远场和近场的情况下，电磁场的性质有所不同，因此，要对远场和近场测量有明确的了解。 1、电磁场的远场和近场划分 电磁辐射源产生的交变电磁场可分为性质不同的两个部分，其中一部分电磁场能量在辐射源周围空间及辐射源之间周期性地来回流动，不向外发射，称为感应场；另一部分电磁场能量脱离辐射体，以电磁波的形式向外发射，称为辐射场。 一般情况下，电磁辐射场根据感应场和辐射场的不同而区分为远区场（辐射场）和近区场（感应场）。由于远场和近场的划分相对复杂，要具体根据不同的工作环境和测量目的进行划分，一般而言，以场源为中心，在三个波长范围内的区域，通常称为近区场，也可称为感应场；在以场源为中心，半径为三个波长之外的空间范围称为远区场，也可称为辐射场。近区场通常具有如下特点： 近区场内，电场强度与磁场强度的大小没有确定的比例关系。即：E=377H。一般情况下，对于电压高电流小的场源(如发射天线、馈线等)，电场要比磁场强得多，对于电压低电流大的场源(如某些感应加热设备的模具)，磁场要比电场大得多。 近区场的电磁场强度比远区场大得多。从这个角度上说，电磁防护的重点应该在近区场。 近区场的电磁场强度随距离的变化比较快，在此空间内的不均匀度较大。 远区场的主要特点如下： 在远区场中，所有的电磁能量基本上均以电磁波形式辐射传播，这种场辐射强度的衰减要比感应场慢得多。 在远区场，电场强度与磁场强度有如下关系：在国际单位制中，E=377H，电场与磁场的运行方向互相垂直，并都垂直于电磁波的传播方向。 远区场为弱场，其电磁场强度均较小 近区场与远区场划分的意义： 通常，对于一个固定的可以产生一定强度的电磁辐射源来说，近区场辐射的电磁场强度较大，所以，应该格外注意对电磁辐射近区场的防护。对电磁辐射近区场的防护，首先是对作业人员及处在近区场环境内的人员的防护，其次是对位于近区场内的各种电子、电气设备的防护。而对于远区场，由于电磁场强较小，通常对人的危害较小。 对我们最经常接触的从短波段30MHz到微波段的3000MHz的频段范围，其波长范围从10米到1米。 2、远区场的测量 在远区场（辐射场区），可引入功率密度矢量（波印廷矢量），电场矢量、磁场矢量、波印廷矢量三者方向互相垂直，波印廷矢量的方向为电磁波传播方向。 在数值上，E=377H，S=EH=E2/377。其中电场强度E的单位是(V/m)，磁场强度H的单位是(A/m)，功率密度的单位是（W/m2）,全部是国际单位制(SI)。 由公式可看出，在远场区，电场与磁场不是独立的，可以只测电场强度，磁场强度及功率密度中的一个项目，其他两个项目均可由此换算出来。 一般情况，关于远场和近场的测量问题可以简化为： 国标规定，当电磁辐射体的工作频率低于300MHz时，应对工作场所的电场强度和磁场强度分别测量。当电磁辐射体的工作频率大于300MHz时，可以只测电场强度。 300MHz频率相应的波长为1米，λ/6为16cm，16cm之外辐射场占优势。如按3λ的划分界限，距辐射源3米之外可认为是远场区。 一般电磁环境是指在较大范围内由各种电磁辐射源，通过各种传播途径造成的电磁辐射背景值，因而属于远区场，辐射的频谱非常宽，电磁场强度均较小。 1GHz以下远区辐射场的测量，可用远区场强仪，也可用干扰场强仪。 中国无线论坛中卫出品 MAC绑定IP，DHCP 关闭，MAC过滤，SSID隐藏 解决方案初探 中国无线论坛中卫出品 声明：任何不经别人同意而进入别人网络都是非法和不道 德的行为。 本教程用于学习和交流，如由此产生一切违法行为与本教程 无关。 题记： 本人是中国无线论坛https://www.360docs.net/doc/dc12505908.html,/的ID ―中卫‖。本教程是根据网上现有的资料还有我的理解加上 实际操作实践编辑整理而成。 由于本人也是初学者，缺乏专业的理论知识，因此文中 不免存在理解的偏差甚至错误，希望各位朋友指正。 如果对教程有任何意见和建议，欢迎各位https://www.360docs.net/doc/dc12505908.html, 论坛提问和交流。 谢谢！！ 中国无线论坛中卫出品 本文的书写是基于一个测试环境。由于测试环境限制，并不能真正的体现实 际情况的复杂程度。为了更接近实际情况，我手动设置尽量复杂。希望更多的朋友用此方法进行测试和应用，最后得出真正的解决方法。 很多朋友都遇到这样的问题，自己很辛苦的破解了对方的WEP 甚至WPA 密码，但是对方的AP设置了SSIDSSID隐藏，MAC过滤，关闭了DHCP，甚至MAC绑定IP。让你蹭网的梦想变成打击。 今天我将和大家一起学习和交流有关这几个问题的解决方案。 前提条件，你已经破解了对方的WEP 或WPA密码，网络上必须有合法客户端，并且客户端在进行通信。 我破解的AP的MAC是00:14:6c:3e:f0:ac 客户端MAC是00:16:b6:9d:10:ad 一．MAC地址绑定： 首先MAC地址绑定，如果对方是无客户端的，那你根本破解不了。因为无 客户端的破解一定要注入攻击，要注入攻击就必须建立虚拟连接。你如果不知道对方的合法客户端的MAC地址，你是不能建立虚拟连接的，不能建立虚拟连接就不能注入。所以MAC 地址绑定的WEP 破解是需要有客户端的，在监听的时候能获得合法客户端的MAC，然后把自己的MAC修改为对方的MAC就能实现正常连接。 二．DHCP关闭，MAC绑定IP，子网掩码 DHCP 关闭的AP，你在连接的时候提示受限，不能正常获得有效的IP。网 上已经有朋友出了用科莱网络分析系统软件来抓包的解决方法，这个方法经过我的测试是可以的。我把方法和各位分享。 首先你正常连接AP，会出现对话框让你输入密码，你输入正确的密码，最 后出现受限，这时候你手动随便设定一个IP 地址，我手动设定一个比较复杂的IP地址。27.122.1.100，子网掩码设为255.0.0.0，网关和DNS空着，如下图所示再连接你已破解的WEP。这时候下面的无线网络连接会显示正常连接。 中国无线论坛中卫出品 附录 COMSOL Multiphysics 的MATLAB 矢量计算基础 W. B. J. ZIMMERMAN 1，J. M. REES 2 1 Department of Chemical and Process Engineering, University of Sheffield, Newcastle Street, Sheffield S1 3JD United Kingdom 2 Department of Applied Mathematics, University of Sheffield, Hicks Building, Sheffield 矢量计算支撑了偏微分方程和它们的数值近似求解。为了很好的使用有限元方法，建模人员应该掌握矢量计算基础知识。本科毕业的工程师可能学过矢量计算的数学课程，但是由于没有碰到过矢量计算的实际应用，这时在工程建模中使用矢量计算就受到限制。本附录介绍了所有COMSOL MULTIPHYSICS WITH MATLAB 中用到的矢量计算基础知识。所以也可以将该附录当作是COMSOL MULTIPHYSICS WITH MATLAB 多变量微分计算的入门读本。当我们写该附录时曾经争论过是否将这部分内容直接加入到第一章（数值分析基础）中，因为导数的数值近似是偏微分方程求解的基础，而偏微分方程是COMSOL MULTIPHYSICS 的基本运算单元。确实，在学习波谱法求解偏微分方程时，基本理论就是“导数理论”——如何使用波变换方法来近似导数。所以通过对比发现，有限元方法的基础就是数值微分。所以争论就不存在了，第一章主要是关于COMSOL MULTIPHYSICS 直接计算的基本问题的。但是不管多有用，近似导数仍然只是建模的一个中间步骤，不是目标本身。 我们这里只考虑用于矢量计算的MATLAB 基础，本附录的重点在于特征值分析和逻辑表达式。这些在整本书中都有体现。应当注意到我们这里介绍的每个功能都可以在COMSOL Script 中实现。本书中唯一不能在COMSOL Script 中实现的Matlab 命令就是fminsearch 。 1．矢量回顾 1.1 矢量表达 FEMLAB 可以处理标量、矢量和矩阵数据，这里简单介绍一下矢量的表达（作为MATLAB 矩阵数据类型的一个特例）。标量可以作为一个单独的数，但是矢量是具有大小和方向的。在如图1所示的右手坐标系系统中，向量a 用以下形式表达： 123123(,,) a a a a a a =++=a i j k a (1) 这里i ，j 和k 是坐标方向的单位矢量，1a ，2a ，3a 是向量a 在各轴方向上的分量。它们是a 对各单位矢量i ，j 和k 的投影。对于坐标系中的P 点（x ，y ，z ），矢量P 对于初始坐标系统O 的位置为： (,,) x y z x y z =++=r i j k (2) MATLAB 用分量的形式描述列矢量或行矢量： >> a = [1; 2; 3]; % column vector 看很多人咨询修改固件问题决定做一个简易教程！使用软件：1.winhex 下载链接：winhex10.rar(493.48 KB, 下载次数: 5) 1.启动界面 启动界面后选择打开你需要的固件文件就可以了 2.固件文件结构简单讲解 编程器固件中包含uboot、fw和art。 uboot就像电脑的bios，是底层的管理系统； fw就像电脑的操作系统，实现路由器的各种功能； art就像电脑的无线驱动程序，是无线校验码； uboot的长度为128KB（0x20000）；art为64KB（0x10000）；fw有4M和8M的区别，4M的为3840K （0x3c0000）,8M的为7936KB（0x7c0000）。在刷机之前要对要刷入flash的uboot、fw、art的文件长度ultraedit或winhex进行校验。尤其是uboot，如果大小不对，千万不要尝试刷入，那是一定会变砖的。 uboot、fw和art在flash中的位置如下： 4M的FLASH：flash地址从0x000000~0x3FFFFF ttl访问flash的地址从0x9F000000~0x9F3FFFFF flash起始地址 TTL起始地址 flash终止地址 TTL终止地址 uboot 0X000000 0X9F000000 0X01FFFF 0X9F01FFFF fw 0X020000 0X9F020000 0X3DFFFF 0X9F3DFFFF art 0X3F0000 0X9F3F0000 0X3FFFFF 0X9F3FFFFF 8M的FLASH：flash地址从0x000000~0x7FFFFF ttl访问flash的地址从0x9F000000~0x9F3FFFFF flash起始地址 TTL起始地址 flash终止地址 TTL终止地址 uboot 0X000000 0X9F000000 0X01FFFF 0X9F01FFFF fw 0X020000 0X9F020000 0X7DFFFF 0X9F7DFFFF art 0X7F0000 0X9F7F0000 0X7FFFFF 0X9F7FFFFF ADAMS/View中系统提供的数学函数大致分类介绍如下。 （1）基本数学函数 ABS(x) 数字表达式x的绝对值 DIM(x1，x2) x1>x2时x1与x2之间的差值，x1<x2时返回0 EXP(x) 数字表达式x的指数值 LOG(x) 数字表达式x的自然对数值 LOG10(x) 数字表达式x的以10为底的对数值 MAG(x，y，z) 向量[x，y，z]求模 MOD(x1，x2) 数字表达式x1对另一个数字表达式x2取余数 RAND(x) 返回0到1之间的随机数 SIGN(x1，x2) 符号函数，当x2>0时返回ABS(x)，当x2<0时返回－ABS(x) SQRT(x) 数字表达式x的平方根值 （2）三角函数 SIN(x) 数字表达式x的正弦值 SINH(x) 数字表达式x的双曲正弦值 COS(x) 数字表达式x的余弦值 COSH(x) 数字表达式x的双曲余弦值 TAN(x) 数字表达式x的正切值 TANH(x) 数字表达式x的双曲正切值 ASIN(x) 数字表达式x的反正弦值 ACOS(x) 数字表达式x的反余弦值 ATAN(x) 数字表达式x的反正切值 ATAN2(x1，x2) 两个数字表达式x1，x2的四象限反正切值 （3）取整函数 INT(x) 数字表达式x取整 AINT(x) 数字表达式x向绝对值小的方向取整 ANINT(x) 数字表达式x向绝对值大的方向取整 CEIL(x) 数字表达式x向正无穷的方向取整 FLOOR(x) 数字表达式x向负无穷的方向取整 NINT(x) 最接近数字表达式x的整数值 RTOI(x) 返回数字表达式x的整数部分 位置/方向函数位置/方向函数用于根据不同输入变量计算有关位置或方向的参数。ADAMS/View中系统提供的位置/方向函数分类介绍如下。 （1）位置函数 LOC_ALONG_LINE 返回两点连线上与第一点距离为指定值的点 LOC_CYLINDRICAL 将圆柱坐标系下坐标值转化为笛卡儿坐标系下坐标值 LOC_FRAME_MIRROR 返回指定点关于指定坐标系下平面的对称点 LOC_GLOBAL 返回参考坐标系下的点在全局坐标系下的坐标值 LOC_INLINE 将一个参考坐标系下的坐标值转化为另一参考坐标系下的坐标值并归一化 LOC_LOC 将一个参考坐标系下的坐标值转化为另一参考坐标系下的坐标值 药学全文电子期刊导航 A Advanced Drug Delivery Reviews《先进药物输送评论》荷兰 ISSN:0169-409X，1987年创刊，全年12期，Elsevier Science 出版，SCI、EI收录期刊，SCI 2003年影响因子6.588，2003年EI收录103篇。综论先进的药物输送系统的设计、研究和进展及其在实验和临床治疗中的应用，以及对兽医学和农业的影响等。 Archiv der Pharmazie《制药文献》德国 ISSN:0365-6233，1822年创刊，全年12期，John Wiley，SCI收录期刊，SCI 2003年影响因子0.624。刊载制药学、药剂学、毒物学、医用化学与药用植物等方面的研究论述。 B Biomedicine & Pharmacotherapy《生物医学与药物疗法》法国 ISSN:0753-3322，1947年创刊，全年10期，Elsevier Science 出版，SCI收录期刊，SCI 2002年影响因子1.035，2003年EI收录552篇。刊载生物医学和药物疗法，包括与之相关的生化、科、药学、外科、生物物理、免疫、血液、神经与分泌等方面的研究论文、评论和札记。 Bioorganic & Medicinal Chemistry《生物有机化学与医药化学》英国 ISSN:0968-0896，1993年创刊，全年24期，Elsevier Science 出版，SCI收录期刊，SCI 2002年影响因子2.043。刊载生物化学包括生物有机化学和医学及其相关学科的研究论文与评论。涉及受体、通路、酶、核苷酸、脂质、糖类等重要生物靶中的分子相互作用。 Bioorganic & Medicinal Chemistry Letters《生物有机化学与医药化学通讯》英国 ISSN:0960-894X，1991年创刊，全年24期，Elsevier Science 出版，SCI收录期刊，SCI 2002年影响因子2.051。刊载生物有机化学、医药化学领域的实验与理论研究成果简报，报道药物设计和发展的最新进展。 BT3破解无线WEP/WPA教程 声明：任何不经别人同意而进入别人网络都是非法和不道德的行为。本教程用于学习和交流，如要实验请拿自已的AP开刀！！ 题记： 本人是中国无线论坛https://www.360docs.net/doc/dc12505908.html,/的ID“中卫”， 无线安全版块是本论坛一个特殊而重要的版块，我们一直非常努力的想把这个板块做好。作为板块现阶段的的一个重点就是无线WEP 和WPA的破解内容。我根据各位坛友的教程和自己的理解整理编辑成这篇《BT3破解无线WEP/WPA教程》。 由于本人也是初学者，缺乏专业的理论知识，因此文中不免存在理解的偏差甚至错误，希望各位朋友指正。 最后希望更多的朋友参与到教程的整理和编辑中，不断把教程修正和完善。 如果对教程有任何意见和建议，欢迎各位到https://www.360docs.net/doc/dc12505908.html,论坛提问和交流。 谢谢！！ 中卫 08年7月13日 开放式WEP破解 1．装备：IMBX60笔记本（内置Intel3945无线网卡）、BT3的U盘系统（需用syslinux命令来指定启动BT3的盘符） 2．用户名：root密码：toor，进入图形界面：startx。启动BT3后,(启动黑屏：xconf再输入startx) 3．加载3945网卡的驱动。打开一个shell 输入modprobe –r iwl3945卸载原来的网卡驱动 输入modprobe ipwraw 加载可监听的网卡驱动 注：不同的网卡有不同的加载方式 LINUX驱动是通过模块进行加载的，可以用lsmod来查看机器已加载的模块 然后可以通过modinfo ipwraw（模块名）来查看所加载驱动模块的版本信息 最新的ipwraw的版本是ipwraw-ng-2.3.4-04022008.tar.bz2的。 ADAMS常用函数的说明 一、几个常用函数的说明 1、 STEP函数 格式：STEP (x, x0, h0, x1, h1) 参数说明： x ―自变量，可以是时间或时间的任一函数 x0 ―自变量的STEP函数开始值，可以是常数或函数表达式或设计变量； x1 ―自变量的STEP函数结束值，可以是常数、函数表达式或设计变量； h0 ― STEP函数的初始值，可以是常数、设计变量或其它函数表达式； h1 ― STEP函数的最终值，可以是常数、设计变量或其它函数表达式。 2、 IF函数 格式：IF(表达式1: 表达式2, 表达式3, 表达式4) 参数说明： 表达式1－ADAMS的评估表达式； 表达式2－如果的Expression1值小于0，IF函数返回的Expression2值； 表达式3－如果表达式1的值等于0，IF函数返回表达式3的值； 表达式4－如果表达式1的值大于0，IF函数返回表达式4的值； 例如：函数IF(time-2.5:0,0.5,1) 结果：0.0 if time < 2.5 0.5 if time = 2.5 1.0 if time > 2.5 3、AKISPL函数 格式：AKISPL (First Independent Variable, Second Independent Variable，Spline Name, Derivati ve Order) 参数说明： First Independent Variable ——spline中的第一个自变量 Second Independent Variable(可选) ——spline中的第二自变量 Spline Name ——数据单元spline的名称 Derivative Order(可选) ——插值点的微分阶数，一般用0就可以了 例如： function = AKISPL(DX(marker_1, marker_2), 0, spline_1) spline_1用下表中的离散数据定义： 综合性医药卫生类核心期刊 1 中华医学杂志19 中国医科大学学报 2 第四军医大学学报20 中山大学学报.医学版 3 北京大学学报。医学版 21 西安交通大学学报.医学版 4 第二军医大学学报22 江苏医药 5 第三军医大学学报23 新医学 6 解放军医学杂志24 天津医药 7 中国医学科学院学报25 山东大学学报.医学版 8 复旦学报.医学版26 南京医科大学学报.自然科学版 9 华中科技大学学报.医学版27 军医进修学院学报 10 吉林大学学报.医学版28 陕西医学杂志 11 中南大学学报.医学版29 中国现代医学杂志 12 四川大学学报.医学版30 广西医科大学学报 13 第一军医大学学报31 山东医药 14 苏州大学学报.医学版32 郑州大学学报.医学版 15 广东医学33 哈尔滨医科大学学报 16 上海医学34 安徽医科大学学报 17 军事医学科学院院刊35 北京医学 18 上海第二医科大学学报 36 医学与哲学 预防医学、卫生学类核心期刊 1 中华流行病学杂志14 辐射防护 2 中华预防医学杂志15 环境与健康杂志 3 中华劳动卫生职业病杂志16 中国媒介生物学及控制杂志 4 卫生研究17 中国卫生统计 5 营养学报18 环境与职业医学 6 工业卫生与职业病19 中国老年医学杂志 7 中国职业医学20 中国学校卫生 8 中国公共卫生21 中国医院管理 9 中国工业医学杂志22 现代预防医学 10 中国计划生育学杂志23 中国医院管理杂志 11 中国人兽共患病杂志24 中国消毒学杂志 12 中国辐射卫生25 中国妇幼保健 13 卫生毒理学杂志 中国医学类核心期刊 1 中草药 10 北京中医药大学学报 2 中国中药杂志11 中国医药学报 3 中国中西医结合杂志12 上海中医药杂志 4 中成药 13 辽宁中医杂志 5 中药材 14 新中医 6 中医杂志15 时针国医国药 7 中国针灸16 陕西中医 8 中国中医基础医学杂志 17 江苏中医药 9 中药药理与临床18 四川中医 能量释放率的计算是基于最小势能原理推导来的，而势能=应变能－外载荷做功－ 摩擦力做功（可能转化为热能）－其他能（如声能等）。一般情况下，如果没有转化为其他能量的话，那么此时的裂纹扩展的断裂力学参量就是应变能释放率；如果不是，则称为能量释放率。 断裂能是材料固有的特性，和断裂韧性是一致的。能量释放率是实际的裂纹扩展参数，是一个动态演化的变量，它表明了裂纹推进一定长度需要的能量。 [断裂与失效] abaqus断裂韧性扫盲贴[复制链接] 最近论坛上很多人都开始搞xfem，里面有几个参数比较令人痛苦，比如断裂韧性KIC（fracture thoughness）和裂纹表面能G（fracture energy）。 首先要知道的是，这两个参数都不是通过abaqus仿真能得到的，而是材料本身的特性，跟密度一样，是先天决定的。所以你要想得到这两个数值，只有两个途径:查文献，或者做实验。 第二，解释一下这两个属性的定义。 早在1920年，格里菲斯就从能量平衡的观点研究了玻璃的脆断。他证明了，对于弹性体中预先存在的一条裂纹，当总位能的减小等于或超过两个新的裂纹表面的表面能时，裂纹就会发生扩展。之后他又研究了裂纹尖端附近的应力场，发现当裂纹前端的应力场强度达到材料的某一临界值时，裂纹就将发生扩展。前者的临界能量就是裂纹表面能，后者的临界应力场就是断裂韧度。 众所周知，裂纹大体分为三类，I,II,III，文献3中有详述。这里只对第一类裂纹举例。 如上图所示的裂纹中，应力强度因子的表达式为 由此可知，应力强度因子与裂纹尖端附近区域内点的坐标无关，它与应力场有关，与裂纹的形状和裂纹的尺寸及方向有关，与载荷的大小和作用方向有关，与材料的某些常数有关，所以应力强度因子可以有效地反映裂纹尖端应力场强度。用abaqus可以算出裂纹尖端的应力强度因子如下图所示，但是临界应力强度因子既断裂韧性却是一个固定的常数，是材料本身的属性，需要做试验确定。 裂纹表面能与应力强度因子的关系为 MAC绑定IP，DHCP关闭，MAC过滤，SSID隐藏 解决方案初探 声明：任何不经别人同意而进入别人网络都是非法和不道德的行为。 本教程用于学习和交流，如由此产生一切违法行为与本教程无关。 题记： 本人是中国无线论坛https://www.360docs.net/doc/dc12505908.html,/的ID “中卫”。本教程是根据网上现有的资料还有我的理解加上实际操作实践编辑整理而成。 由于本人也是初学者，缺乏专业的理论知识，因此文中不免存在理解的偏差甚至错误，希望各位朋友指正。 如果对教程有任何意见和建议，欢迎各位https://www.360docs.net/doc/dc12505908.html, 论坛提问和交流。 谢谢！！ 本文的书写是基于一个测试环境。由于测试环境限制，并不能真正的体现实际情况的复杂程度。为了更接近实际情况，我手动设置尽量复杂。希望更多的朋友用此方法进行测试和应用，最后得出真正的解决方法。 很多朋友都遇到这样的问题，自己很辛苦的破解了对方的WEP甚至WPA 密码，但是对方的AP设置了SSIDSSID隐藏，MAC过滤，关闭了DHCP，甚至MAC绑定IP。让你蹭网的梦想变成打击。 今天我将和大家一起学习和交流有关这几个问题的解决方案。 前提条件，你已经破解了对方的WEP或WPA密码，网络上必须有合法客户端，并且客户端在进行通信。 我破解的AP的MAC是00:14:6c:3e:f0:ac 客户端MAC是00:16:b6:9d:10:ad 一．MAC地址绑定： 首先MAC地址绑定，如果对方是无客户端的，那你根本破解不了。因为无客户端的破解一定要注入攻击，要注入攻击就必须建立虚拟连接。你如果不知道对方的合法客户端的MAC地址，你是不能建立虚拟连接的，不能建立虚拟连接就不能注入。所以MAC地址绑定的WEP破解是需要有客户端的，在监听的时候能获得合法客户端的MAC，然后把自己的MAC修改为对方的MAC就能实现正常连接。 二．DHCP关闭，MAC绑定IP，子网掩码 DHCP关闭的AP，你在连接的时候提示受限，不能正常获得有效的IP。网上已经有朋友出了用科莱网络分析系统软件来抓包的解决方法，这个方法经过我的测试是可以的。我把方法和各位分享。 首先你正常连接AP，会出现对话框让你输入密码，你输入正确的密码，最后出现受限，这时候你手动随便设定一个IP地址，我手动设定一个比较复杂的IP地址。27.122.1.100，子网掩码设为255.0.0.0，网关和DNS空着，如下图所示 再连接你已破解的WEP。这时候下面的无线网络连接会显示正常连接。 ScienceDirect 是世界著名的学术期刊出版商Elsevier公司开发的互联网上最全面的一个全文文献数据库，内容涵盖数学、物理、生命科学、化学、计算机、临床医学、环境科学、材料科学、航空航天、工程与能源技术、地球科学、天文学、及经济、商业管理、社会科学等几乎所有学科领域，提供Elsevier公司出版的1,800多种学术期刊的检索和全文，以及其它著名组织和STM出版商的期刊。 CSA剑桥科学文摘数据库(免费) 括60多个数据库，覆盖的学科范围包括：生命科学、水科学与海洋学、环境科学、计算机科学、材料科学以及社会科学。检索结果为文献的题录文摘信息。 Springer 是世界上著名的科技出版集团, 通过Springer Link系统发行电子图书并提供学术期刊的检索服务。Springer Link系统分11种学科，并由此构成了Springer Link的11个全文电子图书馆。其出版形式包括印刷版和光盘版及网上版数据库。 Chemical Sciences (化学): 37种 UnCover web (免费)IngentaConnect 从17000种期刊中免费检索文章(目前该服务被称为ingenta)(免费) 资源简介【URL】https://www.360docs.net/doc/dc12505908.html,/ 美国化学工程师学会 美国化学工程师学会 (AIChE) 成立于 1908 年，是一家拥有 40,000 多名会员的专业组织，领导着化学工程职业的发展。其会员都是富有创意的问题解决专家，他们利用自己的科学技术知识开发流程、设计和运营工厂，从而在合理的成本内生产有用的产品。 国家科学数字图书馆CSDL 中国科学院国家科学数字图书馆(Chinese National Science Digital Library 简称 BIOS: https://www.360docs.net/doc/dc12505908.html,/main.asp---『bios大本营』 https://www.360docs.net/doc/dc12505908.html,/→ BIOS维修网站←（有点难连的） https://www.360docs.net/doc/dc12505908.html,/---BIOS专页 IT--book: https://www.360docs.net/doc/dc12505908.html,/szulib ... .htm---------------网上电脑教程 IT报纸杂志 https://www.360docs.net/doc/dc12505908.html,/microcomputer/index.pcshow---欢迎访问《微型计算机》网站 https://www.360docs.net/doc/dc12505908.html,/---电脑报 -- 全国发行量第一的计算机报 https://www.360docs.net/doc/dc12505908.html,/index.htm----电脑爱好者 https://www.360docs.net/doc/dc12505908.html,/https://www.360docs.net/doc/dc12505908.html, 编程 https://www.360docs.net/doc/dc12505908.html,/soft/program/book/index.html---[海阔天空下载站]--编程书库 https://www.360docs.net/doc/dc12505908.html,/-----问专家 - VB编程, VC编程, C#编程, Delphi编程, C++编程.. https://www.360docs.net/doc/dc12505908.html,/vvbksd/index.asp-------计算机图书大全 https://www.360docs.net/doc/dc12505908.html,/wenzhang06/defaul t.asp---C语言之家 综合 https://www.360docs.net/doc/dc12505908.html,/---中华电脑书库 https://www.360docs.net/doc/dc12505908.html,/article/sort/c_all_1.htm----所有文章 https://www.360docs.net/doc/dc12505908.html,/jiaocheng/---绿岛----教程天地！！！ https://www.360docs.net/doc/dc12505908.html,/72339069014638592/index.shtml---天极网---软件频道 https://www.360docs.net/doc/dc12505908.html,/softxue/--------☆新时代软件教程：操作系统主页制作服务器设计软件网络技术编程语言文字编辑 https://www.360docs.net/doc/dc12505908.html,/main.htm-------步步高电脑教程站首页：：安阳信息港专业电脑学习栏目 https://www.360docs.net/doc/dc12505908.html,/te/te.htm--------计算机网络教程 https://www.360docs.net/doc/dc12505908.html,/index.php----------软件技巧与提示－－－在这儿你会成为高手 https://www.360docs.net/doc/dc12505908.html,/download.asp--------电子图书免费电子图书电子图书下载 ebook（有的是英文书） https://www.360docs.net/doc/dc12505908.html,/netbook/netdnsj.htm-------电脑书籍（内容有点老，但很多是基础内容） https://www.360docs.net/doc/dc12505908.html,/---------IT电脑顾问专家 -- https://www.360docs.net/doc/dc12505908.html,（有点像是问题解答） https://www.360docs.net/doc/dc12505908.html,/--------IT人的社区 https://www.360docs.net/doc/dc12505908.html, https://www.360docs.net/doc/dc12505908.html,/--------系统分析之窗 DOS: https://www.360docs.net/doc/dc12505908.html,/dosmain.htm---新DOS时代 https://www.360docs.net/doc/dc12505908.html,/dosbbs/----中国DOS联盟之联合DOS论坛--论坛首页 IT考试及认证: https://www.360docs.net/doc/dc12505908.html,/bbs/----------流星技术论坛 https://www.360docs.net/doc/dc12505908.html,/www/index.asp---中国IT认证实验室---国内权威IT认证资讯、技术信息库！ https://www.360docs.net/doc/dc12505908.html,/china/traincert/default.asp---Microsoft培训与认证 SOFT开发编程: C语言 四大要务 第一要务：学习如何发帖 研读《ABAQUS版发帖必读》和《ABAQUS版版面小提示【新手必看】》，学会如何发帖。“入乡随俗”是世间和谐的铁律。所以你初来aba版第一件事应该是看这两贴。有人用“紧急求助”“救命”为题，让人很是生气：别人光看题目也不知道你是什么问题，就算解决了也不能为后人搜索到。所以，题目一定要显示问题之所在。另外，论坛规则，不得讨论盗版问题，也请大家遵守。 第二要务：下载顶置的《ABAQUS常见问题汇总-2.0版》 这是aba-aba斑竹根据常见问题而精心总结的，也是论坛里面各位高手集体智慧的结晶。就我的经验，新手所提的大部分问题都能在其中找到解决方案，这样的话，节约了版面，也使问题解决的速度加快，利版利民！我以前也发过一些2.0里面有板书的帖子，后来仔细读了2.0，对它佩服的五体投地，真实我们新手的指明灯啊！有时，有些简单问题高手不高兴搭理，就是因为这类问题被问过太多次，比如如何给实体单元加弯矩啊（实体单元没有转动自由度得先建reference point，然后耦合）；如何施加随时间变化的载荷（amplitude）2.0中有板书。 第三要务：版内搜索 和第二步一样，很多问题早已经得到解决，所以遇到问题（称之为你的“盲点”）可以先版内搜索，一样是加快问题解决速度。而且，这样你还会发现一个“副产品”：由于同样是这个问题，别人提起的方式和角度不一样，各位高手解答的方式和角度也是各异（虽然异曲同工），你都浏览一次，能使你对这个盲点有全方位的认识！盲点快速变成两点。比如你的帮助文件不能搜索了该怎么办？（我总结过一次）。出现“应变速度大于波速”怎么办？----版上已经有Robert大侠的总结了，新手多半是第一条（单位问题）。 第四要务：学会看帮助文件Documentation 其实帮助文件时最好的辅导教材，可惜是English版，很多人望而生畏（包括我），在论坛下载了好多中文的书籍和例子学习，后来经论坛里面的各位版主提醒就硬着头皮去看Documentation。我才发现Documentation没那么难看懂，而且发现Documentation写得很有条理，很有章法。你一旦掌握了这个章法，无论你看到哪个地方，都很容易上手。所以，看Documentation就两点经验：1，对英文自信一点；2，看看Documentation编写结构，下次一看就明白了。对于你不熟悉比较模糊的专业词汇，我给大家推荐一个专业词汇的网站中国期刊网字典，查看各行各业的专业词汇（中国期刊网字典，可以查各行各业的专业词汇，我就很喜欢，建议大家去试试） （二）学习步骤 我具有有限元初步知识和3D-CAD学习经历，没有其他有限元软件经历，有一定的数学和力学功底（材料力学、弹性力学）。因此我的经验对这类网友比较有借鉴价值，其他网友可参考大米斑竹写得《ABAQUS版版面小提示【新手必看】》（顶置） （1）书籍 我觉得我学习aba过程中，书记方面对我帮助最大的是石亦平博士的《ABAQUS有限元分析实例详解》，《abaqus在机械工程中的应用》也不错。庄老的书也很好，不过我觉得不好懂，并不适合新手。可以通过这些书首先从软件操作的角度掌握分析一个模型的CAE流程，进ADAMS中的函数_New
微波仿真论坛_贴片天线研究
adams常用函数
药学文献检索资料1
微波仿真论坛_电磁场的远场和近场划分
无线论坛档
微波仿真论坛附录COMSOLMultiphysics的MATLAB矢量计算基础18页
路由器固件修改教程(个人理解)-无线路由区-中国无线论坛 -
ADAMS中的函数
药学全文电子期刊导航
BT3研究学习无线WEP和WPA增加版教程
ADAMS部分常用函数的说明
药学类核心期刊
能量释放率&断裂能(摘自simwe论坛)
MAC绑定IP，DHCP关闭，MAC过滤解决方案初探
药学数据库说明
IT论坛
给simwe论坛新手的一些建议和经验