24小时热门版块排行榜    

查看: 2199  |  回复: 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

金虫 (正式写手)

【答案】应助回帖


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

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的回帖
信息提示
请填处理意见