预应力锚杆建立的简便方法，欢迎各位大神交流指导。。。
本帖最后由 涛声依旧2015 于 2015-9-24 17:13 编辑刚开始发的格式有问题，再发一次：
想要问的问题在最后面希望大家能一起探讨，交流，能在百忙之中给予指导意见。。
三种方式模拟预应力锚索：
方式1、通过删除－建立link模拟。
删除掉锚索端头原来CABLE和ZONE之间自动建立的连接，然后建立钢性连接。示例命令如下：
sel cable id=1 beg0, 0, 0 end 0 ,29, 0 nseg 10
sel cable id=1 beg0,29,0 end 0,35,0 nseg 6
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e5 gr_k 2e7 range cid 11,16
sel delete link range id 1
sel link id=100 1 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel cable id=1 pretension 60e3 range cid 1,10
方式2、通过赋三段属性模拟。即将锚索的端头、自由段、锚固段赋不同的属性来模拟预应力锚索。示例命令如下：
sel cable id=1 beg0, 0, 0 end 0 ,29, 0nseg 10
sel cable id=1 beg0,29,0end 0,35,0 nseg 6
sel cableprop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10
sel cableprop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e5 gr_k 2e7 range cid 11,16
sel cableprop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e8 gr_k 2e10 range cid 1,16
sel cable id=1 pretension 60e3 range cid 1,10
方式3、借助别的结构单元（如liner单元）来模拟托盘。删除掉锚索端头的link，然后建立新的link，新link的target为liner上的node。示例命令如下：
sel cable id=1 beg0, 0, 0 end 0 ,29, 0nseg 10
sel cable id=1 beg0,29,0end 0,35,0 nseg 6
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10
sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e5 gr_k 2e7 range cid 11,16
sel liner range y=-.1, .1 x=-1,1 z=-1,1
sel liner PROP iso=( 25e9, 0.15) thick=0.1
sel liner PROP cs_nk=8e8 cs_sk=8e8 &
cs_ncut=0.0 cs_scoh=0.0 cs_scohres=0.0 cs_sfric=0.0
sel delete link range id 15
sel link id=100 1 target node tgt_num 188
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel cable id=1 pretension 60e3 range cid 1,10
我重点想向各位大神请教的是还有其它建立预应力锚杆的方法么？这几种方法大体相似，缺点是当你需要建立很多预应力锚杆时，这样建立很费劲，写的程序不简洁，，哪位大神有建立预应力锚杆更好的方法求交流指导（比如像彭文斌实用教程最后一章讲的最后一个基坑开挖实例，实例过程中建立的预应力锚杆方法就感觉方法很好，很简洁方便，但是看的不大懂，有哪位高手懂这个的程序同样希望你能在百忙之中给予指导），不胜感激
都是参数化建模。
例如我之前自己算的一个基坑例子：
1. 首先定义参数，数据结构怎么个构造看自己了。
def para_sel
; 略
; 略
; 略...
; 14 锚索 cable
e_14 = 2.1e11
xcarea_14 = 2*7*pi*0.0025^2
ytension_14 = 1.35e6
gr_fric_14 = 0
gr_k_14_1 = 1 ;2e10 ; 锚头段 ; 单位长度上水泥浆的刚度
gr_coh_14_1 = 1 ;1e10 ; 单位长度上水泥浆粘结力
gr_k_14_2 = 1 ; 中间段
gr_coh_14_2 = 1
gr_k_14_3 = 2e7 ; 锚固段
gr_coh_14_3 = 2e5 ;
gr_per_14 = pi*0.15
; 锚索预应力 pretension
arr_pre = get_array(6,2)
arr_pre(1,1)=270e3*0.75
arr_pre(2,1)=280e3*0.75
arr_pre(3,1)=290e3*0.75
arr_pre(4,1)=270e3*0.75
; 略
; 略
; 略...
; 锚索位置参数 position
arr_cable_w = get_array(6, 4) ; 高程，长度, 锚固长度, 角度
arr_cable_w(1,1) = 41.0
arr_cable_w(1,2) = 19.5
arr_cable_w(1,3) = 6.0
arr_cable_w(1,4) = 10.0
; 略
; 略
; 略...
end
@para_sel
2. 然后安装，参数、数据结构根据自己情况定。
def install_cable_w(flag)
delta_L1 = 1.5
delta_X1 = -0.7817
delta_Y1 = -37.4919
alpha_1 = atan(delta_X1/delta_Y1)
x1 = 0.00
y1 = 0.00
z1 = arr_cable_w(flag, 1)
z2 = z1-arr_cable_w(flag,2)*sin(arr_cable_w(flag,4)*pi/180.0)
z3 = z1-(arr_cable_w(flag,2)-arr_cable_w(flag,3))*sin(arr_cable_w(flag,4)*pi/180.0)
cable_id = 1100 + flag
loop n(1,24)
x1 = x1 - delta_L1 * sin(alpha_1)
y1 = y1 - delta_L1 * cos(alpha_1)
p_gp1 = gp_near(x1,y1,z1)
x11 = gp_xpos(p_gp1)
y11 = gp_ypos(p_gp1)
z11 = gp_zpos(p_gp1)
x2 = x11 - arr_cable_w(flag,2)*cos(arr_cable_w(flag,4)*pi/180.0)
y2 = y11
x3 = x11-(arr_cable_w(flag,2)-arr_cable_w(flag,3))*cos(arr_cable_w(flag,4)*pi/180.0)
y3 = y11
command
sel cable id @cable_id begin @x11 @y11 @z11 end @x3 @y3 @z3 nseg 5
sel cable id @cable_id begin @x3@y3@z3 end @x2 @y2 @z2 nseg 10
endcommand
endloop
pre_w = arr_pre(flag,1)
command
sel cable id @cable_id prop emod @e_14 xcarea @xcarea_14 ytension @ytension_14 &
gr_k @gr_k_14_3 gr_coh @gr_coh_14_3 gr_per @gr_per_14 range z @z2 @z3
sel cable id @cable_id prop emod @e_14 xcarea @xcarea_14 ytension @ytension_14 &
gr_k @gr_k_14_2 gr_coh @gr_coh_14_2 gr_per @gr_per_14 range z @z3 @z11
sel cable id @cable_id pretension @pre_w range z @z3 @z11
endcommand
end
3. 然后用起来就是这样
; 开挖西部第一层 z1 = 40.5 z2 = 42.5
@excav(40.5,42.5)
@install_cable_w(1)
@relink(30000)
solve
save quarter_excav_w1
彭文斌实用教程最后一章讲的最后一个基坑开挖实例，实例过程中建立的预应力锚杆方法跟平时用的传统方法都不同，有哪位高手指导是怎么施加的么？求指点迷津。 三种方法建立的锚杆运行结果一样么？会不会在运行的时候出现锚杆折断这种现象啊？ 569051814 发表于 2015-9-24 16:32
三种方法建立的锚杆运行结果一样么？会不会在运行的时候出现锚杆折断这种现象啊？ ...
三种方式建立运行的结果基本都是一样的~"锚杆折断"???不太理解你的意思！！在有些情况下地质条件可能太差，地质参数可能过小，只用锚杆支护可能仍然会失稳~~~还有施加的预应力也不可能太大，最起码不能超过规范上锚杆使用时的最大抗拉强度~~~ 涛声依旧2015 发表于 2015-9-24 17:18
三种方式建立运行的结果基本都是一样的~"锚杆折断"???不太理解你的意思！！在有些情况下地质条件可能太差 ...
假设锚杆水平建立的，运行后看到锚杆形状这样了，一一一∧一一一？这是怎么回事啊？是锚杆参数太小了？ 569051814 发表于 2015-9-26 15:51
假设锚杆水平建立的，运行后看到锚杆形状这样了，一一一∧一一一？这是怎么回事啊？是锚杆参数太小了？ ...
我以前做的没有出现过你这种情况，觉得不应该出现这种问题，可能是你的锚杆设置的方向有问题~~ 涛声依旧2015 发表于 2015-9-26 21:46
我以前做的没有出现过你这种情况，觉得不应该出现这种问题，可能是你的锚杆设置的方向有问题~~ ...
运行的时候会提示锚杆大变形之类的，不知道你有没有这种情况。 569051814 发表于 2015-9-28 10:19
运行的时候会提示锚杆大变形之类的，不知道你有没有这种情况。
一般情况下都没有，可能是你的地址参数设置的太小，或者锚杆设置的有问题，好好检查检查~~~ 同求锚杆问题 echowasd 发表于 2015-10-1 01:01
都是参数化建模。
例如我之前自己算的一个基坑例子：
这些命令流绝对是好东西呀，要好好学学一下，谢谢分享呀，必须点个赞！！！ :lol:lol:lol:lol:lol:lol:lol:lol:lol:lol:lol 感谢分享，新手下载学习 echowasd 发表于 2015-10-1 01:01
都是参数化建模。
例如我之前自己算的一个基坑例子：
您好厉害，学习了，谢谢分享 为什么我模拟预应力锚杆老是报错，就是用的上面三个方法中的第一个。
页:
[1]