UDEC实例详解[1]
![UDEC实例详解[1]](https://img.sodocs.net/img3c/01tt1pog5nyp2lgetcq9-c1.webp)
![UDEC实例详解[1]](https://img.sodocs.net/img3c/01tt1pog5nyp2lgetcq9-e2.webp)
实例1 荷载作用下边坡稳定问题
朔准线黄河大桥，边坡高105m，层理间距取3m，岩层视倾角6°；主节理间距取2m，视倾角85°，计算模型如下：
DEM计算模型
相关计算参数取值如下：桥基长15m，宽12m，荷载180000kN；灰岩弹性模量E=4×104MPa，泊松比μ=0.23，岩体密度γ=2.7×103kg/m3。
new
round 0.05
set delc off
block 0 0 0 50 50 50 62.6 115.2 68.9 123.9 85.5 132.5 &
105.5 140.1 125.4 156.4 180 156.4 180 0
jset 85,0 300,0 0 0 3,0 (50 50)
jset 6 0 300 0 0 0 4 0 (50 50)
gen auto 5
change jmat=1 range angle 84 86
change jmat=2 range angle 5 7
prop mat=1 den=2700 b=1.5e9 s=0.6e9
prop jmat=1 jkn=2e9 jks=2e9 jcoh=1e9 jfr=25 jten=1e5 节理参数不同prop jmat=2 jkn=2e9 jks=2e9 jcoh=1e10 jfr=35
grav 0 -10
hist n=100 xdis 68.9 123.9 ydis 68.9 123.9
hist n=100 xdis 125.4 156.4 ydis 125.4 156.4
hist unbal
bound xvel=0 range -1 0.2 -1 50.1
bound xvel=0 range 224.9 225.1 -1 157 bound yvel=0 range -1 226 -1 0.1 cyc 5000
save nature.sav
UDEC (Version 3.10)
LEGEND
3-Jun-07 23:46 cycle 100000 history plot
-6.93E-03
0.00E+00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
(e+02)
-0.20
0.00
0.20
0.40
0.60
0.80
1.00
1.20
(e-01) JOB TITLE :
Itasca Consulting Group, Inc. Minneapolis, Minnesota USA
变坡点x 位移
UDEC (Version 3.10)
LEGEND
3-Jun-07 23:46 cycle 100000 history plot
-1.25E-01
0.00E+00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
(e+02)
-1.40
-1.20
-1.00
-0.80
-0.60
-0.40
-0.20
0.00
0.20
(e-01) JOB TITLE :
Itasca Consulting Group, Inc. Minneapolis, Minnesota USA
变坡点y 位移
UDEC (Version 3.10)
LEGEND
3-Jun-07 23:46 cycle 100000 history plot
-3.80E-03
0.00E+00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
(e+02)
-0.20
0.00
0.20
0.40
0.60
0.80
1.00
1.20
1.40
(e-01) JOB TITLE :
Itasca Consulting Group, Inc. Minneapolis, Minnesota USA
坡顶点x 位移
UDEC (Version 3.10)
LEGEND
3-Jun-07 23:46 cycle 100000 history plot
-1.99E-01
0.00E+00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
(e+02)
-2.00
-1.80
-1.60
-1.40
-1.20
-1.00
-0.80
-0.60
-0.40
-0.20
0.00
0.20
(e-01) JOB TITLE :
Itasca Consulting Group, Inc. Minneapolis, Minnesota USA
坡顶点y 位移
UDEC (Version 3.10)
LEGEND
3-Jun-07 23:46 cycle 100000 history plot
0.00E+00
0.00E+00
0.00 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
(e+02)
0.00
0.20
0.40
0.60
0.80
1.00
1.20
1.40
1.60
1.80
2.00
(e+06) JOB TITLE :
Itasca Consulting Group, Inc. Minneapolis, Minnesota USA
不平衡力
实例2 边坡渗流稳定问题
问题定义：
; ================================================
; --- fluid flow test run ---
; --- slope : 10 m high ---
; --- 2 joint sets : 20 and 80 deg.; friction = 30 deg.
;; --- r.h.s. water level : 6 m --- no failure ---
; 8 m
; 9 m
; 10 m --- failure ---
; ================================================
;
round 0.05
定义刚体角边圆角半径，缺省值0.5，建议为块体边长的1%，设定后不可更改。set delc off
如果接触块体张开距离大于0.55倍round定义的长度时，接触将被自动删除（缺省值on），off 则表示不自动删除。
block 0,-5 0,0 5,0 11,10 23,10 23,-5
生成一个块体，每次分析该命令只能用一次。
jset 20,0 100,0 0,0 2,0 (5,1)
设置节理性质，20为角度，100为长度，2为间距，（5，1表示起点位置）
jset 80,0 100,0 0,0 3,0 (5,0)
;
del area 0.1
;删除面积小于0.1的块体
gen auto 10
;
prop mat=1 dens=0.0025 k=16667 g=10000
定义材料属性：密度(de nsity)、体积模量(b ulk)，剪切模量(sh ear)（当cons=3、6时，还需要coh esion粘聚力、di lation剪胀角、fri ction内摩擦角、te nsion拉张强度）
prop mat=1 jkn=10000 jks=10000 jfric=45.0
prop mat=1 jperm=1.0e8 azero=0.0005 ares=0.0002
; 定义材料属性：法向刚度系数，切向刚度系数，摩擦角（度），节理渗透系数、法向应力为零时的裂隙开口大小、高应力作用下裂隙开口大小（长度）
Jcons=2 (面面接触，库仑滑动准则)还有一些参数：如：jcohesion，jdilation，jtension Jcons=1（点接触）参数：ccohesion(内聚力)，cdilation(剪胀角)，cfriction(摩擦系数)，cperm(渗透系数)，ctension(拉张强度)，kn(法向刚度)，ks(切向刚度) Jcons=5(面面接触，库仑滑动准则，残余强度) ，与jcons=2相比，多了以下参数：jrescoh(残余内聚力)，jrfric(残余摩擦角)，jrtens(残余拉张强度)
对于jcons=2和5，可以指定非线性应力应变关系，用property nstable n命令，Table n (x,y) (x1,y2)…….
change jmat=1 jcons=2
;设置材料对应关系，如节理采用编号为1的材料，节理性质为2（1，2，3，5，7）
insitu str -0.125 0 -0.25 ygrad 0.0125 0 0.025
;
初始化应力，str表示(0,0)处应力，sxx0,sxy0,syy0, 分别表示x方向应力，剪应力，y方向应力
Ygrad表示变化梯度，sxxy,sxyy,syyy
Sxx=sxx0+sxxx*x+sxxy*y
Sxy=sxy0+sxyx*x+sxyy*y
Syy=syy0+syyx*x+syyy*y
还可以给出Xgrad，同上。
bound xvel=0 range -1,1 -6 1
bound yvel=0 range -1,31 -6 -4
bound xvel=0 range 21,23 -6 11
;
定义边界条件
外部边界条件：1）荷载边界：bound xload fx range x,y……..
bound yload fy range x,y,………
2) 自由边界：bound xfree
bound yfree
3) 应力边界：bound stress sxx0,sxy0,syy0 xgrad sxxx,sxyx,syyx ygrad sxxy,sxyy,xyyy
4) 位移速度边界bound gvel vx0,vy0,vxx,vxy,vyx,vyy
Vx=vx0+vxx*x+vxy*y
Vy=vy0+vyx*x+vyy*y
Bound nvel vn 变形块体法向位移速度
Bound svel vs 变形块体剪切位移速度
Bound xvel vx 变形体x方向位移速度
Bound yvel vy 变形体y方向位移速度
内部边界条件：1）荷载边界：bound interior xload fx range x,y……..
bound interior yload fy range x,y,………
2) 自由边界：bound interior xfree
bound interior yfree
3) 应力边界：bound interior stress sxx0,sxy0,syy0 xgrad sxxx,sxyx,syyx ygrad sxxy,sxyy,xyyy
4) 位移速度边界bound interior gvel vx0,vy0,vxx,vxy,vyx,vyy
Vx=vx0+vxx*x+vxy*y
Vy=vy0+vyx*x+vyy*y
Bound interior nvel vn 变形块体法向位移速度
Bound interior svel vs 变形块体剪切位移速度
Bound interior xvel vx 变形体x方向位移速度
Bound interior yvel vy 变形体y方向位移速度
grav 0 -10
;x、y方向的重力加速度，同set gravity x, y
hist n=100 xdis 11,10 ydis 11,10
保存变量历史，n=100表示第循环100时间步，存储一次，
Xdis 11,10 表示存储距11，10最近的点的x方向位移；
Ydis 11,10表示存储距11，10最近的点的y方向位移；
hist unbal
; 记录最大不平衡力
Solve
求解
save slfl1.sav
;保存文件
; --------------------------------------
;
; flow --- r.h.s. water at y=6m
;
; --- no failure ---
;
; --------------------------------------
fluid dens 0.001
定义流体密度，注意单位要统一
bound imperm range -1,31 -6 -4
设置无渗流的边界范围，x最小最大范围，y最小最大范围
bound pygrad -0.01 range -1,1 -6 1
设置流体压力梯度范围
bound pp=0.06 pygrad -0.01 range 21,23 -6 6
;设置流体压力梯度范围，pp表示压力初值，pygrad表示梯度range 表示范围。
set flow steady
;设置稳定流模式
hist xdis 5.92 1.54 ydis 5.92 1.54
记录位移历史
reset displ jdis
设置块体和节理位移为零
prop mat=1 jfric=25
设置材料1的节理摩擦角为25度。
cyc 500
执行时间500步。
save slfl2.sav
;
; --------------------------------------
;
; flow --- r.h.s. water raised to y=8m
;
; --- some slip, but no block failure ---
;
; --------------------------------------
bound pp=0.08 pygrad -0.01 range 21,23 -6 8 cyc 500
save slfl3.sav
;
; --------------------------------------
;
; flow --- r.h.s. water raised to y=9m
;
; --- no failure ---
;
; --------------------------------------
bound pp=0.09 pygrad -0.01 range 21,23 -6 9 cyc 1000
save slfl4.sav
;
; --------------------------------------
;
; flow --- r.h.s. water raised to y=10m
;
; --- failure ---
;
; --------------------------------------
bound pp=0.10 pygrad -0.01 range 21,23 -6 11 cyc 15000
save slfl5.sav
return
大玄空挨星法详解
大玄空挨星法详解 大玄空挨星只以坐山（24山）的挨星入中飞星，而论八方零正吉凶；同一山向上下元排盘不同。无运星、向星入中。 凡坐山挨星是1、2、3、4，逢上元运则顺飞九宫，下元运则逆飞九宫。（1、2、3、4上顺下逆） 凡坐山挨星是6、7、8、9，逢上元运则逆飞九宫，下元运则顺飞九宫。（6、7、8、9上逆下顺） 阴阳宅以此顺九宫次序（中五→六→七→八→九→一→二→三→四宫）挨排，飞星到九宫，各宫飞到之星在上元或下元九十年中是不动的。只要山向不变，入中之星永远不变。 如：阳宅寅山申向，寅庚丁三山本山星是八白，即用八白星入中。无论是在上元或下元建造，上元用八白星入中逆行排盘，到下元八白入中顺行排盘。 又如壬山丙向，坤壬乙三山本山星是二黑，即以二黑星入中依次。无论是在上元或下元建造，上元用二黑星入中顺行挨排，到下元二黑星入中逆行挨排。 演示举例： 1、阳宅午山子向，用“午酉丑九紫右弼守”诀，午山用九紫，九紫即为本山星，用九紫入中挨排。 2、在下元时，九紫属阳星，则入中顺飞九宫。
九紫居中5，一白到干6，二黑到兑7，三碧到艮8，四绿到离9，五黄到坎1，六白到坤2，七赤到震3，八白到巽4，各宫飞到之星在下元九十年不动。 3、如是上元时，九紫则入中逆飞九宫。九紫居中5，八白到干6，七赤到兑7，六白到艮8，五黄到离9，四绿到坎1，三碧到坤2，二黑到震3，一白到巽4，各宫所置之星在上元九十年亦不动。阴宅各宫所挨到之星在上元九十年中亦不动。 阳宅：午山子向本山星为九紫入中 上元逆行挨星下元顺行挨星 四、大玄空强调合局法则（正神正位装，拨水入零堂）。正神、零神所到之处，必须与地势虚实空满相配合，结合分析。否则犯零正颠倒，仅会排盘亦是枉然。 上元以1、2、3、4为正神，应高满，喜来水、来风、来气；6、7、8、9为零神，应空虚，喜出水、出风、出气。 下元以6、7、8、9为正神，，应高满，喜来水、来风、来气；1、2、3、4为零神，应空虚，喜出水、出风、出气。 山水形势讲究合局，也就是“水里龙神不上山，山上龙神不下水”，不得形局，妄谈挨星。认金龙，一经一纬妙无穷，动不动，直待高人施妙用。 五、各宫九星的吉凶，依合局法则判断。挨星之九星本身并无吉凶之分，其吉凶结果完全是根据合局（正神正位装，拨水入零堂）或反局（上山
大玄空挨星法详解(1)
大玄空挨星法详解 一、大玄空论元运是分上元和下元（上下二元八运），三元分作二元，五运分46各管前后十年。（五黄前10年归4运，后10年归6运。）其实还是三元九运，只不过将五运切半，前半为上元，后半为下元。《青囊三字经》曰：“大玄空，用九星”。即用九星顺逆挨排法来揭示自然天运对人类社会的影响，以破解二宅吉凶之谜。大玄空藉上下元运之不同，将九星分为阴阳，阳顺阴逆。即阳一路，阴一路。 二、大玄空挨星阴阳宅用星不同。有阳宅和阴宅挨星二种，二者不相混也。 二十四山的挨星，指用挨星真诀排出的坐山之星，阳宅入中立极排盘，挨星用本山星；阴宅入中立极排盘，挨星用父母星。父母星由本山星逆经四位（987654321←）而来。并随上下元运之不同而颠颠倒、顺逆行。用挨星真诀而起九宫飞星掌排。 挨星真诀本山星（阳宅） 坐甲癸申山，贪狼一路行；一白为本山星，上元用一白入中顺行，下元一白入中逆行。 坐坤壬乙山，巨门从头出；二黑为本山星，上元用二黑入中顺行，下元二黑入中逆行。 坐子卯未山，三碧禄存到；三碧为本山星，上元用三碧入中顺行，下元三碧入中逆行。 坐戌干巳山，四绿文曲照；四绿为本山星，上元用四绿入中顺行，下元四绿入中逆行。 坐辰巽亥山，六白武曲位，六白为本山星，上元用六白入中逆行。下元六白入中顺行。 坐艮丙辛山，七赤是破军；七赤为本山星，上元用七赤入中逆行。下元七赤入中顺行。 坐寅庚丁山，八白左辅应；八白为本山星，上元用八白入中逆行。下元八白入中顺行。 坐午酉丑山，九紫右弼星。九紫为本山星，上元用九紫入中逆行。下元九紫入中顺行。 ----------------------------------------------------------------------------- 挨星真诀父母星（阴宅） 七赤为一白之父母星（1→7），上元用七赤入中逆行。下元七赤入中顺行。 八白为二黑之父母星（2→8），上元用八白入中逆行。下元八白入中顺行。 九紫为三碧之父母星（3→9），上元用九紫入中逆行。下元九紫入中顺行。 一白为四绿之父母星（4→1），上元用一白入中顺行，下元一白入中逆行。 二黑为五黄之父母星（5→2），上元用二黑入中顺行，下元二黑入中逆行。 三碧为六白之父母星（6→3），上元用三碧入中顺行，下元三碧入中逆行。 四绿为七赤之父母星（7→4），上元用四绿入中顺行，下元四绿入中逆行。 五黄为八白之父母星（8→5），上元用五黄入中顺行，下元五黄入中逆行。 六白为九紫之父母星（9→6），上元用六白入中逆行。下元六白入中顺行。
《中华大玄空风水学》与《三合法》辨伪
《中华大玄空风水学》与《三合法》辨伪.txt同志们：别炒股，风险太大了，还是做豆腐最安全！做硬了是豆腐干，做稀了是豆腐脑，做薄了是豆腐皮，做没了是豆浆，放臭了是臭豆腐！稳赚不亏呀！《中华大玄空风水学》与《三合法》辨伪 杨公盘诀云：“收山出煞用三元，玄空大卦妙如仙。时师昧却杨曾诀，三合充满宇宙间”。目前三合法依然盛行，传之已久，流毒广深，现附载三合水法辨伪。 长生三合水法 甲生亥库未加卯木局乙生午库戌加寅火局 丙生寅库戌加午火局丁生酉库丑加巳金局 庚生巳库丑加酉金局辛生子库辰加申水局 壬生申库辰加子水局癸生卯库未加亥木局 上法即乙丙交而趋戌，辛壬会而趋辰，斗牛纳丁庚之气，金羊收巽之灵，俱从向上起遗却山龙。此即僧一行李虚中测命之术，所称乾巽艮坤四卦具属勉强牵合三方全无理气。 又向上起黄泉水法。 诀曰：庚丁坤上起黄泉，乙庚损伤癸水先，甲癸向上忧见艮，辛壬水路怕当乾。俱称黄泉之
说本稽谬之名，诀无可取之证。 其反变黄泉水诀： 坤向庚丁切莫逢，巽嫌乙丙祸重重，艮向甲癸须当忌，乾怕辛壬最为凶。 既曰杀人，复曰救贫自相刺谬，不辨自明，又八曜水从本卦起浑天甲子过克爻则为官鬼，迹涉卜筮阴阳元窍，毫无相关，且四季不全，入于无著。如坎龙坤免震山猴，巽鸡乾马兑蛇头，艮虎离猪之例。 皆以克山者为煞，张九仪宗之以伪传，愈传愈讹，使各家假书仍讹袭谬，其实于相互至理，动辙相违断不可用。至于净阴净阳。纳甲水法，生成水法皆言朔望盈虚，河图对待，说似近理，然阴阳不交，卦非一家，且日月系天地中阴阳之端，专就太阴而欲概言天地之太阴太阳则已遗本逐末，尚可望合而凝乎，原学者尚思太极之理：万物一极一物，一极一动，一静互为其根于阴阳交媾中打破元关，望毋视为秘本则幸甚，三合凶水断法例。 乾戌水为鼓盆煞，克妻朦胧又嘴哑，乙辰水来主贼盗，男为盗贼女守寡。 太极定纳甲法也，从龙上本卦乾兑离震坎巽坤互对掌诀，自本卦互对宫起贪狼，次及巨禄文廉武破辅等星辅归本，所谓卦变八而归原是也，向上辅星即从向上起辅，乾山则从巽上起辅，次及武破廉贪巨禄文等星，细查惟乾巽二卦龙与向其所翻三吉顺逆相同，余卦向与龙自相矛盾，颠倒吉凶殊难取信至河洛寓生成法。如一得九而合十，三得七而合十，二得八而合十，四得六而合十。又如坎向得壬水即为后天夫妇生成之数，坤向得丙，乙向得艮，亦然余卦纳甲仿此。 净阴净阳水法，以太阴为定其论与平砂玉尺合然，以此论卦义之八变则可而与主太阴以论元窍则窒滞不通，即一行所著之经流传滋蔓学者毋为其所惑，此为谬理也。