24小时热门版块排行榜    

查看: 2330  |  回复: 5

实体建模初学

新虫 (初入文坛)

[求助] 如何位移加载

我做的模型用力加载的话 荷载位移曲线后面段不正确 想改用位移加载

我就用王新敏那本书上p485 那个简支梁的例子试试看

现在有个问题是 力加载的时候是加在一个面上
原来的命令流是Asel,s,loc,z,0,50
Asel,r,loc,y,300
sfa,all,1,pres,q0
那我位移加载的时候是不是就该成
Asel,s,loc,z,0,50
Asel,r,loc,y,300
da,all,uy,-6

求节点反力的时候 支座上的节点就211 213 216 三个点
不知道我下面这样命令流哪里错了 我得到的荷载位移曲线还是错了 节点反力求出来的荷载比实际小了一倍  可是我明明乘过2了
/POST26
NSOL,2,205,U,Y
Prod,3,2,,,,,,        -1

RFORCE,4,211,F,Y
RFORCE,5,213,F,Y
ADD,6,4,5

RFORCE,7,216,F,Y
ADD,8,6,7
Prod,9,8,,,,,,        2/1000

PRVAR,3,9
Xvar,3                  !定义变量3为x轴
Plvar,9


下面是整条命令流

finish
/clear
/config,nres,2000
/prep7
et,1,solid65,,,,,,,1
et,2,link8
mp,ex,1,13585
mp,prxy,1,0.2
fc=14.3
ft=1.43
tb,concr,1
tbdata,,0.5,0.95,ft,-1  !关闭压碎
tb,miso,1,,11
tbpt,,0.0002,fc*0.19
tbpt,,0.0004,fc*0.36
tbpt,,0.0006,fc*0.51
tbpt,,0.0008,fc*0.64
tbpt,,0.0010,fc*0.75
tbpt,,0.0012,fc*0.84
tbpt,,0.0014,fc*0.91
tbpt,,0.0016,fc*0.96
tbpt,,0.0018,fc*0.99
tbpt,,0.002,fc
tbpt,,0.0033,fc*0.85
mp,ex,2,2.0e5
mp,prxy,2,0.3
tb,biso,2
tbdata,,300,0
pi=acos(-1)
r,1,0.25*pi*22*22
r,2,0.25*pi*22*22/2
r,3,0.25*pi*10*10
r,4,0.25*pi*10*10/2
r,5
!  创建几何模型
blc4,,,150/2,300,2000/2
*do,i,1,9
wpoff,,,100
vsbw,all
*enddo             !切出箍筋位置
wpcsys,-1
wpoff,,,50
vsbw,all            !切出拟加载面
wpcsys,-1
wprota,,-90
wpoff,,,30
vsbw,all           !切出拉区钢筋竖向位置
wpoff,,,240
vsbw,all           !切出压区钢筋竖向位置
wpcsys,-1
wpoff,30
wprota,,,90
vsbw,all           !切出钢筋水平位置
wpcscy,-1
!划分钢筋网格
elemsiz=50
lsel,s,loc,x,30
lsel,r,loc,y,30
cm,zj,line
latt,2,1,2
lesize,all,elemsiz    ! 拉区外侧钢筋
lsel,s,loc,x,75
lsel,r,loc,y,30
cm,zjb,line
latt,2,2,2
lesize,all,elemsiz    !拉区中间钢筋
lsel,s,loc,x,30
lsel,r,loc,y,270
cm,jlj,line
latt,2,3,2
lesize,all,elemsiz    !压区钢筋
lsel,s,tan1,z
lsel,r,loc,y,30,270
lsel,r,loc,x,30,70
lsel,u,loc,z,50   
lsel,u,loc,z,0      
cm,gj,line
latt,2,3,2
lesize,all,elemsiz    !对称位置外的箍筋
lsel,s,loc,z,0
lsel,r,loc,y,30,270
lsel,r,loc,x,30,70
cm,gjb,line
latt,2,4,2
lesize,all,elemsiz    !对称位置处的箍筋
lsel,all
cmsel,s,zj
cmsel,a,zjb
cmsel,a,jlj
cmsel,a,gj
cmsel,a,gjb
cm,gj,line
lmesh,all

lsel,all
/eshape,1
eplot
!划分混凝土网格
vatt,1,5,1
mshkey,1
esize,elemsiz
vmesh,all
allsel,all

lsel,s,loc,y,0
lsel,r,loc,z,900
dl,all,,uy           
asel,s,loc,z,0
da,all,symm
asel,s,loc,x,75
da,all,symm

Asel,s,loc,z,0,50
Asel,r,loc,y,300
da,all,uy,-6
Allsel,all

/solu
antype,0                 !静态分析
nsubst,100                !子步数
outres,all,all              !输出所有求解选项
autos,on                !使用自动时间步
neqit,50                !设置每个荷载子步的迭代次数为50

solve

[ Last edited by 实体建模初学 on 2012-12-24 at 22:46 ]
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldhltw

金虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
csgt0: 金币+2, 欢迎交流 2012-12-25 11:49:37
如果位移和力士对应的话,可以考虑那么来处理
或是可以检查一下你的命令流总是都有问题,
2楼2012-12-25 10:36:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldhltw

金虫 (正式写手)

【答案】应助回帖

检查对应单位以及数量级
应该不会有那么大的差距
3楼2012-12-25 10:38:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

实体建模初学

新虫 (初入文坛)

引用回帖:
3楼: Originally posted by ldhltw at 2012-12-25 10:38:25
检查对应单位以及数量级
应该不会有那么大的差距

说错了  荷载是小了一半。。。。
4楼2012-12-25 11:12:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ldhltw

金虫 (正式写手)

【答案】应助回帖


csgt0: 金币+1, 谢谢 2012-12-26 13:32:45
看你进行的是混凝土的计算,但是你说只有三个点被约束
可能计算有一定问题吧
那样对应的约束点的力很大
5楼2012-12-25 11:18:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

asdf1

新虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
csgt0: 金币+1, 谢谢 2012-12-26 13:32:52
先加力,后加位移这种,
你应该是在加力结束时,已有位移的基础上,继续增大位移来加载,
如果在加力结束时,直接用da,all,uy,-6,程序会认为你是从0加到这么多。
此外,这种位移加载,最好在求解过程中进行控制。就相当于是多步分析。一步一步的增大位移。
我的任务就是每天来这里胡说八道一番
6楼2012-12-26 12:28:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 实体建模初学 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 G6APbkg8SA6w 2026-08-29 3/150 2026-08-30 02:29 by ZPa0EcMwuECS
[论文投稿] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 gy1nBQXYQJqL 2026-08-29 4/200 2026-08-30 02:09 by ZPa0EcMwuECS
[论文投稿] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 gy1nBQXYQJqL 2026-08-29 6/300 2026-08-30 01:58 by ZPa0EcMwuECS
[公派出国] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 ASdOkHsho7FD 2026-08-28 7/350 2026-08-30 01:12 by ZPa0EcMwuECS
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 ASdOkHsho7FD 2026-08-28 5/250 2026-08-30 00:52 by ZPa0EcMwuECS
[基金申请] 有没有仍没收到信息的 +5 德尚中行 2026-08-27 6/300 2026-08-29 18:10 by manplx
[基金申请] 为什么到现在没收到通知? +4 tannykie 2026-08-29 4/200 2026-08-29 17:14 by lmz0216
[基金申请] 麻烦专家们看看评委们的意见(F口面上) +5 gdd2018 2026-08-28 10/500 2026-08-29 14:10 by Jacob678
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 ASdOkHsho7FD 2026-08-28 7/350 2026-08-29 09:26 by G6APbkg8SA6w
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 申请删除本帖 +6 lyz123lyz 2026-08-27 7/350 2026-08-27 17:31 by 宁静致远sy
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜? 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +7 process2012 2026-08-23 10/500 2026-08-26 19:23 by hmhminy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 今天务委会开完了,明天出结果吗 +19 angus9576 2026-08-25 23/1150 2026-08-26 10:03 by zp519
[基金申请] 国合现在查不到了吗? +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 如果此刻你正在为国基感到焦虑,不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
信息提示
请填处理意见