level 11
![[呵呵]](/static/emoticons/u5475u5475.png)
level 6
level 1
level 1
我想问下关于flac3d温度场的问题，土体在降温冻结的情况下，引起的土体膨胀位移。您有关于这方面的例题，或是关于温度场的命令流吗？
2017年02月23日 04点02分
9
仿真科技论坛上有一个帖子:广州地铁冻结法施工中冻胀的数值模拟及全程的力学分析。可以参考。
2017年02月23日 04点02分
好的 谢谢！
2017年02月23日 06点02分
level 1
问一下，混凝土和围岩之间的接触面单元的参数怎么设置才好
2017年03月03日 07点03分
11
混凝土和岩石之间一般无需设置接触面，在基坑支护中与土体间可设置，因为混凝土与岩石的力学参数相差不大
2017年03月04日 01点03分
level 1
还有混凝土用什么本构模拟比较好呢，怎么算混凝土损伤
2017年03月03日 07点03分
12
混凝土一般就用弹性，你说的损伤是指什么？
2017年03月04日 01点03分
回复 如果没有你luc ：就是混凝土破坏。我现在在做的水工引水隧洞，隧洞内的高压水会把混凝土压坏，带来混凝土损伤，带来混凝土渗透系数的变化，可是flac并没有混凝土损伤模型，所以就不能模拟混凝土开裂渗流这一过程
2017年03月04日 02点03分
@hehexj2008 那你就换软件呗，找个能模拟你想要的
2017年03月04日 02点03分
同问![[啊]](/static/emoticons/u554a.png)
2024年10月10日 07点10分
level 1
能不能问一下壳单元上sel shell apply pressure 这个命令，后面隧洞通水，水垂直作用在上面，压力正负号怎么用。一般是压为负。这个地方是不是默认为压力
2017年03月06日 09点03分
13
level 11
如果没有你luc
楼主
这是Drucker Prager屈服条件的
new
rest 123.f3sav ;计算所得的最后一步结果
dele zone ran mo null ;删除开挖单元
def safety_factor ;定义安全系数
p_z=zone_head ;zone_head 指向单元体表
i=0
array z_coef(10000) ;定义数组
loop while p_z # null ;循环
i=i+1
o_coh=z_prop(p_z,'cohesion') ;得到计算中的 c 值
o_fri=z_prop(p_z,'friction') ;得到计算中的 fai 值
p_ten=z_prop(p_z,'tension') ;抗拉强度值
o_major=-1*z_sig1(p_z) ;计算中 sig1 值
o_middel=-1*z_sig2(p_z) ;计算中的 sig2 值
o_minor=-1*z_sig3(p_z) ;计算中的 sig3 值
o_ten=z_sig3(p_z)
II=o_major+o_middel+o_minor
JJ=((o_major-o_middel)^2+(o_middel-o_minor)^2+(o_minor-o_major)^2)/6
cosf=cos(o_fri*pi/180)
sinf=sin(o_fri*pi/180)
alf=(2*sinf)/(sqrt(3)*(3-sinf))
K=(6*o_coh*cosf)/(sqrt(3)*(3-sinf))
f_coef=abs((K-alf*II)/(sqrt(JJ))) ;得到安全系数
zone_id=z_id(p_z)
if o_ten>p_ten then
f_coef=0 ;材料拉破坏，用 0 做标识
endif
o_xpos=z_xcen(p_z)
o_ypos=z_ycen(p_z)
o_zpos=z_zcen(p_z)
z_coef(i)=string(o_xpos)+' '+string(o_ypos)+' '+string(o_zpos)+' '+string(f_coef)
p_z=z_next(p_z)
endloop
status=open('safety_factor.txt',1,1)
status=write(z_coef,i)
status=close
end
@safety_factor
2017年03月09日 05点03分
15
new
rest 123.f3sav ;计算所得的最后一步结果
dele zone ran mo null ;删除开挖单元
def safety_factor ;定义安全系数
p_z=zone_head ;zone_head 指向单元体表
i=0
array z_coef(10000) ;定义数组
loop while p_z # null ;循环
i=i+1
o_coh=z_prop(p_z,'cohesion') ;得到计算中的 c 值
o_fri=z_prop(p_z,'friction') ;得到计算中的 fai 值
p_ten=z_prop(p_z,'tension') ;抗拉强度值
o_major=-1*z_sig1(p_z) ;计算中 sig1 值
o_middel=-1*z_sig2(p_z) ;计算中的 sig2 值
o_minor=-1*z_sig3(p_z) ;计算中的 sig3 值
o_ten=z_sig3(p_z)
II=o_major+o_middel+o_minor
JJ=((o_major-o_middel)^2+(o_middel-o_minor)^2+(o_minor-o_major)^2)/6
cosf=cos(o_fri*pi/180)
sinf=sin(o_fri*pi/180)
alf=(2*sinf)/(sqrt(3)*(3-sinf))
K=(6*o_coh*cosf)/(sqrt(3)*(3-sinf))
f_coef=abs((K-alf*II)/(sqrt(JJ))) ;得到安全系数
zone_id=z_id(p_z)
if o_ten>p_ten then
f_coef=0 ;材料拉破坏，用 0 做标识
endif
o_xpos=z_xcen(p_z)
o_ypos=z_ycen(p_z)
o_zpos=z_zcen(p_z)
z_coef(i)=string(o_xpos)+' '+string(o_ypos)+' '+string(o_zpos)+' '+string(f_coef)
p_z=z_next(p_z)
endloop
status=open('safety_factor.txt',1,1)
status=write(z_coef,i)
status=close
end
@safety_factor
level 1
楼主，我想问下：在模拟桩基础上拔的时候，怎么确定基础的极限上拔力
2017年03月24日 00点03分
16
与材料有关
2017年03月24日 03点03分
能不能具体一点，我现在要对桩基础施加上拔力 从100KN开始 慢慢增加上拔力，监测位移，怎么知道力到了极限上拔力 ？另外问一下 apply psource -17 range psource 是施加什么东西? 感谢！
2017年03月24日 09点03分
@We_Pleasure 看应力应变曲线，发生塑性变形的点就是极限了
2017年03月24日 14点03分
@We_Pleasure 这个没见过
2017年03月24日 14点03分
![[泪]](/static/emoticons/u6cea.png)
level 1
楼主你好，我用5.0做采矿采场的长锚索支护模拟，长锚索是全长锚固加预应力。但是模拟出来的效果和无支护情况下一样，请问是怎么回事？是因为锚索和采场没有建立连接吗？
2020年07月20日 06点07分
19
大佬能跟你学习一下吗？？
2020年12月24日 09点12分