科学网—xiwangxiaoxue1的留言板

IP: 222.178.10.* [2691]张海权 2017-12-22 22:06

任何的成功都是无条件的，失败比成功需要更多的反思。

IP: 222.178.10.* [2690]张海权 2017-12-22 22:04

和上一个问题相适应的是，既然有了高学历，定位就要高，但是身段一定要低。民国大佬杜月笙说得好，不要怕被人利用，被利用说明你还有用。千万不要端着酸腐的臭架子，觉得有个博士的学历就牛逼得很，非得天天让别人把这玩意儿当头衔挂在口头才有意思。我的老板教育我，Ph.D这玩意儿只能印在名片上，但是如果对方的老人家不叫你小X，那说明你还不够谦恭。现在每年国内毕业的博士中，一多半是中科院985，也就是说，好学校的博士市场上压根不稀罕。甘乃迪总统有言，不要问你的国家能为你做什么，问问你为你的国家能做什么。在公司就得会谈生意，在学校就得会写论文，牛逼的上限可以满天吹，但是实用的下限你得满足得了，你对别人一点儿用都没有，别人怎么会帮助你？所以千万不要把独善其身之类的作为自我标榜的理由，泛道德很多时候只是不作为或者干脆无能的借口罢了。

IP: 222.178.10.* [2689]张海权 2017-12-22 11:37

Numerical simulation of progressive failure in cut slope of soft rock using a soil-water coupled finite element analysis

IP: 222.178.10.* [2688]张海权 2017-12-21 20:35

qualitatively similar, they differed quantitatively.

IP: 222.178.10.* [2687]张海权 2017-12-21 20:13

31--- R. J. Finno, W. W. Harris and G. Viggiani, "Pore pressures and strain localization in stiff
clay," in Fracture and Damage in Quasibrittle Structures, Prague, E & FN SPON,1994, p. 441–449.

IP: 222.178.10.* [2686]张海权 2017-12-21 17:20

Predicting the performance of geosynthetic-reinforced soil retaining walls under working stress conditions and at failure威斯康辛大学2008年博士论文

IP: 222.178.10.* [2685]张海权 2017-12-21 15:31

FLAC 程式用於地下水位之模擬，由於土壤與水之耦合效應為一
複雜問題，因此分析時先將土壤之力學效應關閉（ set mechanic=off），
開啟滲流模式（ set flow=on），待水流趨於穩定後即形成一靜止之地
下水位面，如圖 5.1 所示，此時再開啟力學效應，如此則可開始分析
邊坡之穩定性，當訂定地下水位面後，水位面的孔隙壓力會隨著深度
的增加而增加，如圖 5.2。

IP: 222.178.10.* [2684]张海权 2017-12-21 15:28

由于 STABL 程式未考慮土體之變形性，且假設破壞實為完全剪力破
壞，而 FLAC 在考慮不同之彈性模數與張力強度，由所得結果發現彈
性模數與張力強度對土壤邊坡穩定有其影響性，故 STABL 程式分析
所得之安全係數絕多數大於 FLAC 分析所得，則 STABL 應用在邊坡
穩定分析方面有其限制存在。

IP: 222.178.10.* [2683]张海权 2017-12-21 15:27

STABL 雖有自動輸出所分析邊坡之安全係數與臨界滑動面之優
點，但因其須先假設可能之滑動面，而假設可能滑動面上每一點之安
全係數均相同 (但實際上沿滑動面上之安全係數並非完全相同)，而且
其未考慮土體之變形性，故其分析所得之安全係數可能有高估之情
形。

IP: 222.178.10.* [2682]张海权 2017-12-21 14:36

FLAC 程式基本的咚悴襟E如圖 3.1 所示，首先對每個節點而言，由作用力或應力經由邉臃匠淌疆a生節點的不平衡力 (out-of-balance forces)，利用此不平衡力決定節點的速度分量，並經由對速度積分來決定節點位移；另一方面，對每個區間而言，由節點的速度分量可決定應變率(strain rate)，再利用應力-應變組合律決定新的應力增量，藉此應力增量又可更新節點的不平衡力，如此程式將反覆執行此咚悴襟E，直到節點最大的不平衡力達到一可接受的合理範圍為止。

IP: 222.178.10.* [2681]张海权 2017-12-21 14:16

Flac使用K、G而不使用E、v简述

这是因为人们相信体积模量和剪切模量所对应的材料行为方面比杨氏模量和泊松比所对应的更为基础。（K、G）对于所有不违背热力学原理的弹性模型都有意义，而（E、v）对于某些材料来说是没有意义的。一种极端是材料抵抗体积变化而不抵抗剪切变化，另一个极端是材料抵抗剪切变形而不抵抗体积变形。第一种材料相应于有限的K和G=0，第二种相应于K=0，G为有限值。这两种情况下都不能使用（E、v）来表示。

IP: 222.178.10.* [2680]张海权 2017-12-21 14:15

看来您书上第九章关于接触面的问题。有以下几个问题
1.接触面很容易定义，可以在一个节点还没有两个节点编号之前就定义，也可以在有两个编号以后再定义？
2.接触面定义以后，只要在接触的地方，同一个节点有两个编号，是不是就可以计算。
3.利用切割模型法，当土层有很多层的时候，桩的表面露出土体，时是不是也可以这样做？
4.利用切割模型法，或导来导去法时，能不能用 gen merge 0.001 或是attach face tol 0.001，我对书上那个两个楔形体的例子，做过测试，好像都不能使用这两个命令。

IP: 222.178.10.* [2679]张海权 2017-12-21 14:04

自定义本构模型的基本方法
1 头文件(usermodel.h)中进行新的本构模型派生类的声明
修改模型的ID(>100)、名称和版本
修改派生类的私有成员：私有变量和成员函数
2  C++文件(usermodel.cpp)中修改模型结构
(UserModel::UserModel(bool bRegister): ConstitutiveModel)
const char **UserModel:roperties()函数
3  模型的参数名称字符串：用于在模型中显示
const char **UserModel::States()函数
4  计算过程中的状态指示器：单元屈服状态
5  double UserModel::GetProperty()和void UserModel:: SetProperty()函数
6  const char * UserModel::Initialize()函数
参数和状态指示器的初始化，并对派生类声明中定义的私有变量进行赋值
7  const char * UserModel::Run(unsigned nDin,States *ps) 函数
由应变增量计算得到应力增量，从而获得新的应力
8  const char * UserModel::SaveRestore()函数
对计算结果进行保存。
9 程序的调试
在VC++的工程设置中将FLAC3D软件中的EXE文件路径加入到程序的调试范围中，并将FLAC3D自带的DLL文件加入到附加动态链接库(Additional DLLs)中，然后在Initialize()或Run()函数中设置断点，进行调试；
在程 ...

我的回复(2017-12-21 14:12)：f="http://forum.simwe.com/thread-808415-1-91.html" target="_blank">http://forum.simwe.com/thread-808415-1-91.html

我的回复(2017-12-21 14:04)：all c++6.0编译环境中,利用flac3d自带的udm开发包,里面有构建好的框架,把原来udm里面示例的模型删除掉,然后在添加工程菜单中将usercam.cpp 和usercam.h头文件添加进去，编译生成动态链接库。

我的回复(2017-12-21 14:04)：all c++6.0编译环境中，点击新建，选择创建动态链接库，创建一个工程，然后在添加工程菜单中将usercam.cpp 和usercam.h头文件添加进去，编译生成动态链接库。

IP: 222.178.10.* [2678]张海权 2017-12-21 14:02

apply szz -4.5e6 range z 20
apply szz -5.0e6 range z 0
apply sxx -2.5e6 grad 0 0 1.25e4 range x 0
apply sxx -2.5e6 grad 0 0 1.25e4 range x 20
apply syy -2.5e6 grad 0 0 1.25e4 range y 0
apply syy -2.5e6 grad 0 0 1.25e4 range y 20
在上面的边界条件命令中，为什么“x，y 方向有应力渐变grad,而z方向没有呢？”真心的求指点！！！小弟将不胜感激！！！
你这应该是x、y向应力在z向上有渐变，因为你的grad是一个向量0 0 1.25e4，只有z向有值，x、y向是没有值的。同时，你的z向应力在z向上也是有渐变的，因为你在第一句、第二句中已经明确了szz在上、下面上的不同值，只不过对其之间的值未作具体指定；由于没有指定0～20之间的值，所以在你求解初始应力场的时候，程序还是需要通过重力条件作进一步的求解，需要一定的时间；如果你直接指定了z向应力为渐变，可能就就不需要求解，给的值就是它的初始地应力场了。

IP: 222.178.10.* [2677]张海权 2017-12-21 11:43

非均质土对基础的影响是一个很重要的课题，尤其在海洋地基中。土体的不排水强度从海底开始线形增长，通常对于粘土取摩擦角为0，通过粘聚力的线形增加来模拟这种非均质特性。在用FLAC实现的过程中，会有两种方法，一种是精确的确定不排水强度随深度的变化，另外一种是把土体模型深度方向上等分，用等分段中间深度处的不排水强度来代替等分段的不排水强度。对于第一种方法自认为很难实现，2006年岩土力学上有一篇文章“非均质地基承载力及破坏模式的FLAC数值分析”中提到作者是用FISH语言控制实现土的非均质性，但是没有具体讲实现的过程，不知哪位版主或者高手有这样的程序，希望可以提供建议或者提供个思路。关于第二种方法比较容易实现，但是涉及到一个单元的划分的问题。如果单元比较大，等分段是单元高度的几分之一，那么在计算时一个单元有几个不同的强度值，不知会否影响；如果单元划分过小，使得等分段和单元高度一致，则会出现网格过密影响计算时间的问题。

IP: 222.178.10.* [2676]张海权 2017-12-21 11:35

unbalence force的实际意思或计算表达式
可以参看 theory background部分, 当对系统施加外力时,
所施加的力要逐步传给系统的每个节点,在传递过程中, 就存在不平衡力.
随着计算的进行,不平衡力逐渐消散,
gp_xfunbal(p_gp)，是 x方向的不平衡力,其大小应该是逐渐减小的,
gp_xfapp(p_gp)，是节点上的外力,包括体积力,这个力应该是逐渐变大,或者是不变的...

IP: 222.178.10.* [2675]张海权 2017-12-21 11:32

我模型已经算完，想求出每一个网格域里的孔压与竖向总应力的比值，求高手编一个。我自己编的如下，可两个都不行，应该很简单吧，我没学过fish，所以不大会编，求高手帮忙，谢谢了。
你用这个函数看看，用zone指针把孔压和szz数值分别传给变量aaa,bbb。再用z_extra(pnt,1)记录。我觉得可行，有问题还请指正
config zextra 2
def ccc
pnt=zone_head
loop while pnt#null
bbb=z_szz(pnt)
ccc=aaa/bbb ;aaa为孔隙压力，不记得flac用什么变量表示，先查手册……
z_extra(pnt,1)=ccc
pnt=z_next(pnt)
endloop
end
ccc
plot con zextra 1

IP: 222.178.10.* [2674]张海权 2017-12-21 09:58

大家在计算的时候，电脑总是通宵达旦的工作，有时候即便是算完了，还得开到天亮，显然不利于建立节约型社会嘛。
受论坛上若干帖子的启示，今天找到了解决方法，可以实现计算完毕后即时关机，
其实就是一句命令啦，比如说你想在计算完后5秒钟后关机，只要将下面这句话写在你编的命令流的最下方：Sys start shutdown.exe -s -t 5

IP: 222.178.10.* [2673]张海权 2017-12-21 09:52

fish的解释说明
ef find_add
top_head = null
gp_pnt = gp_head
loop while gp_pnt # null
if gp_zpos(gp_pnt) > 0.05 then
new = get_mem(2)
mem(new) = top_head
mem(new+1) = gp_pnt
top_head = new
endif
gp_pnt = gp_next(gp_pnt)
endloop
end
find_add
这中间的top_head 既然是null 了, 为什么还要把new附给top_head ?new在top_head = new 中是表示数列还是什么啊!!究竟top_head = new是在这段语言中是什么意思啊??谢谢啊!!

我的回复(2017-12-21 09:57)：f="http://forum.simwe.com/thread-727880-1-88.html" target="_blank">http://forum.simwe.com/thread-727880-1-88.html

我的回复(2017-12-21 09:54)：次循环时
new = get_mem(2) 开辟了存放两个fish变量的内存，表的首地址付给new
mem(new) = top_head 把top_head的值（实际为空地址）存到所开辟的第一个内存地址中
mem(new+1) = gp_pnt 把gp_pnt 的值存到第二个内存地址中
top_head = new          把第一次所开辟的内存地址付给top_head
第二次循环
new = get_mem(2) 开辟了存放两个fish变量的内存，首地址付给new
mem(new) = top_head 把top_head的值（实际为上一次开辟内存的首地址）存到这次第一个内存地址中
mem(new+1) = gp_pnt 把gp_pnt 的值存到第二个内存地址中
top_head = new          把第二次开辟内存的地址付给top_head
以此类推，可知，
通过循环得到这样一个链表，第一列放上一个所开辟内存的地址，第二列存放满足条件的gp_pnt
2.对于你的问题，为什么还要把new附给top_head
就是要把上一次所开辟的内存地址付给top_head

IP: 222.178.10.* [2672]张海权 2017-12-21 09:50

请教lakewater版主：
在《FLAC/FLAC3D基础与工程实例》一书中讲到自定义本构的调试时有如下一段话
在VC++的工程设置中，讲FLAC3D软件中的EXE文件路径加入到程序的调试范围中，并将FLAC自带的DLL文件加到附加动态链接库(Additional Dlls)中，
有两个不解之处：
关于自定义本够模型的调试：1、我是这样做的，选择启动调试，出现“请指定用于调试会话的可执行文件名称”，我将FLAC3D.EXE程序所在文件位置加入，但随即出现“无调试信息，或者调试信息不匹配”的对话框。但是当直接把生成的自定义本构的dll 文件拷到flac3d文件目录中，一个单元测试，可以完整执行。这里请教：EXE文件是否是FLAC3D.EXE程序呢？我没找到vs2008的工程设置在哪里，具体该怎样设置“程序的调试范围”？
2、怎么将自带的DLL文件加到Additional Dlls中？我查到vc有三种动态链接库，可没有附加动态链接库，附加动态链接库是什么意思呢？
希望能得到您的解答，谢谢！

xiwangxiaoxue1的个人博客分享 http://blog.sciencenet.cn/u/xiwangxiaoxue1 像我这样笨拙的生活，我只赞美那些一面哭泣一面追求的人！

留言板

张海权