fish语言编写预应力锚索群锚赋值方法

whutzll

最近一直被一个问题困扰：
计算完成后使用“plot sel cable foce”显示结果不尽人意
应该在锚索轴向大都是tension，而不应是compression；
使用命令“plot contour displacemnet”显示结果看只是计算了一根锚杆，另一根呢？
以下是命令流：
---------------------------------------------------------------------------------------------------------------------
new
gen zone brick p0 0 0 0 p1 2 0 0 p2 0 0 -0.4 p3 0 1.3 0 size 20 4 13
;
mo mohr
pro bulk 2.68e7 shear 2.92e7 fric 20 coh 2.3e4 ten 3e3
ini dens 2450

;set grav 0 -9.81 0
;set large
;
;ini syy -312448.5 grad 0 24034.5 0
;ini szz -34716.5  grad 0 2670.5 0
;ini sxx -34716.5  grad 0 2670.5 0
;
fix x range x=-0.001 0.001
fix x range x=1.999 2.001
fix z range z=-0.001 0.001
fix z range z=-0.4001 -0.399
fix x y z range y=-0.001 0.001
;
;set mech ratio 1.0e-5
;solve
;sav ini_0.sav
;
;pause
;位移和速度归零
;ini xdisp 0
;ini ydisp 0
;ini zdisp 0
;ini xvel 0
;ini yvel 0
;ini zvel 0
;
;---------------------------------------------------------------------------------------------------------------
;下面设置自由式锚索(Free_Anchored_Cable)
;;---------------------------------------------------------------------------------------------------------------
def constant
    pc=2   ;（自由式锚杆锚固段剖分数）
  a1x=1.01191-0.0001
  a1y=1.3-0.08
  a3x=1.01191-0.0001-0.7*cos(17*degrad)
  a3y=1.3-0.08-0.7*sin(17*degrad)
  a2x=1.01191-0.0001-0.54*cos(17*degrad)
  a2y=1.3-0.08-0.54*sin(17*degrad)
;  a2z=-0.08-0.08*(m-1)

end
constant
;
;生成第1根锚杆
  sel cable id=1 begin a3x,a3y,-0.08 end a2x,a2y,-0.08 nseg=pc
  sel cable id=1 begin a2x,a2y,-0.08 end a1x,a1y,-0.08 nseg=1
  sel cable id=1 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=08.0e8,gr_k=4.1307e5,gr_per=9.42e-3,yt=2.1e8 range cid 1,2
  sel cable id=1 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=1,gr_k=1,gr_per=9.42e-3,yt=2.1e8 range cid 3
  sel delete link id=4
     ;生成自由段端点与岩体的固结点（模拟垫墩效应）
  sel link id=15 4 target zone
  sel link attach xdir=rigid xrdir=rigid ydir=rigid yrdir=rigid &
     zdir=rigid zrdir=rigid
  sel link constit nydeform 6 area=1.0 k=2.84e12 &
     ycomp=4.17e3 ytens=4.17e3 range id=15
  sel cable id=1 pretension=6.67 range cid=3
;
;生成第2根锚杆
  sel cable id=2 begin a3x,a3y,-0.16 end a2x,a2y,-0.16 nseg=pc
  sel cable id=2 begin a2x,a2y,-0.16 end a1x,a1y,-0.16 nseg=1
  sel cable id=2 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=08.0e8,gr_k=4.1307e5,gr_per=9.42e-3,yt=2.1e8 range cid 4,5
  sel cable id=2 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=1,gr_k=1,gr_per=9.42e-3,yt=2.1e8 range cid 6
  sel delete link id=8
     ;生成自由段端点与岩体的固结点（模拟垫墩效应）
  sel link id=16 8 target zone
  sel link attach xdir=rigid xrdir=rigid ydir=rigid yrdir=rigid &
     zdir=rigid zrdir=rigid
  sel link constit nydeform 6 area=1.0 k=2.84e12 &
     ycomp=4.17e3 ytens=4.17e3 range id=16
  sel cable id=2 pretension=6.67 range cid=6

;
step 100
save 0.sav

[ 本帖最后由 whutzll 于 2008-11-21 17:05 编辑 ]

hql1225

第一根锚索受压力啊老大？？

cherrywmh

pretention 好小呀，6.67N？？

whutzll

原帖由 hql1225 于 2008-11-22 11:04 发表
第一根锚索受压力啊老大？？

多谢hql1225 的提醒，自己试着针对锚索设置改动了一下，其中还加了重力，这次是拉力了。
但仍然和上次的问题一样，得不到想要的结果，再把修改后命令流贴出来吧，希望大伙帮忙解决：
new
gen zone brick p0 0 0 0 p1 2 0 0 p2 0 0 -0.4 p3 0 1.3 0 size 20 4 13
;
mo mohr
pro bulk 2.68e7 shear 2.92e7 fric 20 coh 2.3e4 ten 3e3
ini dens 2450
set grav 0 10 0
;
fix x range x=-0.001 0.001
fix x range x=1.999 2.001
fix z range z=-0.001 0.001
fix z range z=-0.4001 -0.399
fix x y z range y=-0.001 0.001
;--------------------------------------------------------------------------------------
;下面设置自由式锚索(Free_Anchored_Cable)
;;-------------------------------------------------------------------------------------
def constant
    pc=2   ;（自由式锚杆锚固段剖分数）
;  a1x=1.0119
;  a1y=1.22
;  a3x=0.3119
;  a3y=0.52
;  a2x=0.4719
;  a2y=0.68
;生成第1根锚杆
COMMAND
sel cable id=1 begin 0.3119,0.52,-0.04 end 0.4719,0.68,-0.04 nseg=pc
sel node cid=4 1.0119,1.22,-0.04
sel cablesel cid=3 id=2 nodes=2,4
  sel cable id=1 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=08.0e8,gr_k=4.1307e5,gr_per=9.42e-3,yt=2.1e8 range cid 1,2
  sel cable id=2 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=1,gr_k=1,gr_per=9.42e-3,yt=2.1e8 range cid 3
  sel delete link id=4
     ;生成自由段端点与岩体的固结点（模拟垫墩效应）
  sel link id=15 4 target zone
  sel link attach xdir=rigid xrdir=rigid ydir=rigid yrdir=rigid &
     zdir=rigid zrdir=rigid
  sel link constit nydeform 6 area=1.0 k=2.84e12 &
     ycomp=4.17e3 ytens=4.17e3 range id=15
  sel cable id=2 pretension=6.67 range cid=3
END_COMMAND
;生成第2根锚杆
COMMAND
sel cable id=3 begin 0.3119,0.52,-0.16 end 0.4719,0.68,-0.16 nseg=pc
sel node cid=8 1.0119,1.22,-0.16
sel cablesel cid=6 id=4 nodes=6,8
  sel cable id=3 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=08.0e8,gr_k=4.1307e5,gr_per=9.42e-3,yt=2.1e8 range cid 4,5
  sel cable id=4 prop den=8960,emod=1.40e11,xcarea=4.91e-6 &
     gr_coh=1,gr_k=1,gr_per=9.42e-3,yt=2.1e8 range cid 6
  sel delete link id=8
     ;生成自由段端点与岩体的固结点（模拟垫墩效应）
  sel link id=16 8 target zone
  sel link attach xdir=rigid xrdir=rigid ydir=rigid yrdir=rigid &
     zdir=rigid zrdir=rigid
  sel link constit nydeform 6 area=1.0 k=2.84e12 &
     ycomp=4.17e3 ytens=4.17e3 range id=16
  sel cable id=4 pretension=6.67 range cid=6
END_COMMAND
end
constant
;
pl sel cable force
step 100
save 0.sav

caikewei1984

太强了人啊，学习中

caikewei1984

zenm hai mei

zhjweix

顶一下:D:D:D:D

wwt900

这个我倒是经历过类似的，最后检查多次，主要还是校对cid对应的情况，锚固端，自由端，尤其是预应力施加之处对应的cid，要完全正确

xiaohaizzz

:):):):):):):):):):)

zzsqzxj

太强了，楼主，学习了