论坛首页    职业区    学术与写作    工程技术区    软件区    资料区    商务合作区    社区办公室
 您好! 欢迎 登录 或 注册 最新帖子 邀请注册 活动聚焦 统计排行 社区服务 帮助
 
  • 帖子
  • 日志
  • 用户
  • 版块
  • 群组
帖子
  • 1934阅读
  • 2回复

关于flac3D书上“15-1冰碛土三轴数值模拟试验”代码的疑问。 [复制链接]

上一主题 下一主题
离线intelna
 

发帖
189
土币
1228
威望
129
原创币
0
只看楼主 倒序阅读 使用道具 楼主  发表于: 2013-01-26
关于flac3D书上“15-1冰碛土三轴数值模拟试验”代码的疑问。

前提:按此代码运行确实可以得到预期的结果
但我按书上的代码自己写了一段,却得不到想要的结果
其中有几句关于边界条件的代码我深有疑问,先把书上的代码附上:

n
;导入分组网格    
impgrid 1
group soil range group rock not
;定义材料性质
model mohr
pro den 2200 bulk 4e7 she 1.5e7 co 8e4 &
fric 28 ten 9e3 dil 10 range group soil
pro den 2900 bulk 2.5e10 she 1.2e10 co 4e6 &
fric 45 ten 2e6 dil 15 range group rock
;设置边界条件    
fix x y z range z -0.1 0.1
fix x y z range z 7.9 8.1; =========Q1=========
fix y
;加载    
app nstress -5e5 range x -0.1 0.1
app nstress -5e5 range x 4.9 5.1;=========Q2======
;设置初始条件
ini zvel -2.5e-5 range z 7.9 8.1
;定义轴向应力求解函数
def axi_stress
。。。。。
。。。。。

关于Q1:此处将实验模型用于加载的上端固定住,貌似不能继续加载了啊?为啥要有这一句呢?
关于Q2:这才是我最纳闷儿的。明明所建模型x方向的尺度为0.1m, 他怎么在固定x值为5的地方呢?

谢谢各位达人


离线intelna

发帖
189
土币
1228
威望
129
原创币
0
只看该作者 1楼 发表于: 2013-01-26
关键是为什么貌似有问题的代码能运行出来!!奇了怪了!
离线yuzhuhebut

发帖
21
土币
1005
威望
14
原创币
0
只看该作者 2楼 发表于: 2013-03-17
你好,还上这个论坛吗?我也在做三轴试验,可是那个元胞自动机是怎么回事啊?找不到的……
努力 加油
快速回复
限100 字节
温馨提示:欢迎交流讨论,请勿纯表情、纯引用!
 
上一个 下一个

      https://beian.mps.gov.cn/ 粤公网安备 44010602012919号 广州半山岩土网络科技有限公司 粤ICP备2024274469号

      工业和信息化部备案管理系统网站