24小时热门版块排行榜    

查看: 1069  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

肉肉的了

新虫 (小有名气)

[求助] 求助三跨连续梁的静态分析,总是出现问题? 已有2人参与

三跨连续梁的有限元模拟总是出现图下的问题

这是我的命令流


FINI
/CLE
/PREP7
W1=0.2
W2=0.2
W3=0.2
T1=0.012
T2=0.012
T3=0.008   !!!!!!!腹板厚

K,1,,W1/2
K,2,T1,W1/2,
K,3,T1,T3/2,
K,4,W3-T2,T3/2,
K,5,W3-T2,W2/2,
K,6,W3,W2/2,
K,7,W3,-W2/2,
K,8,W3-T2,-W2/2,
K,9,W3-T2,-T3/2,
K,10,T1,-T3/2,
K,11,T1,-W1/2,
K,12,,-W1/2

A,1,2,11,12
A,3,4,9,10
A,5,6,7,8
WPCSYS
AGLUE,ALL

WPROTA,,90,
WPOFF,,,T3/2
ASBW,ALL
WPOFF,,,-T3
ASBW,ALL

ET,1,82
LESIZE,ALL,0.005     !网格尺寸
AMESH,ALL
SECWRITE,MYSECTION

!!!!!!!!!!!!!!!!!!!!!!
FINI
/CLE
/CONFIG,NRES,200000
/PREP7
LB=16.5
NE=330
NN=NE+1

DELTL=LB/NE

EM=2.0E11
DENG=7800
GRA=9.8

ET,1,BEAM188
MP,EX,1,EM
MP,NUXY,1,0.2
MP,DENS,1,DENG

SECTYPE,1,BEAM,MESH   
SECREAD,MYSECTION,,,MESH

*DO,I,1,331
N,I,(I-1)*DELTL
*ENDDO
TYPE,1   
MAT,1
ESYS,0
*DO,I,1,330
E,I,I+1
*ENDDO

/ESHAPE,1   
EPLOT  

!施加约束
D,1,UX,,,,,UY
D,111,UY
D,221,UY
D,331,UY
ALLSEL
FINISH
/SOLU
ANTYPE,0

F,51,FY,-P

SOLVE

求助三跨连续梁的静态分析,总是出现问题?


发自小木虫IOS客户端
回复此楼

» 猜你喜欢

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

shgao20

专家顾问 (职业作家)

【答案】应助回帖

这两个约束是用来防止整个梁沿 z 向的刚体运动,和绕梁轴线的刚体转动。
  每个物体都有 6 个刚体自由度,做静力分析时,必须施加足够的约束以防止任何刚体运动。否则就无法进行计算。
  你原来的命令流缺少对这两个自由度的约束,所以添加这两个约束。约束位置不一定使用节点 51,也可以施加在其它节点上。

» 本帖已获得的红花(最新10朵)

5楼2017-06-04 08:28:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

ftc_简单

银虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
约束不足,仔细核对下你的约束!
承接各类优化设计项目,百度乐仿腾讯课堂,海量教程等你来看
2楼2017-06-01 21:43:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shgao20

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
1  约束不足,缺少 z 向位移约束,以及绕量轴线的转动约束。
  2  施加载荷用的 P 没有赋过值。
修改后的命令流如下,可以试试看:
FINI
/CLE
/filn,  test_beam188, 1

/PREP7
W1=0.2
W2=0.2
W3=0.2
T1=0.012
T2=0.012
T3=0.008   !!!!!!!腹板厚

K,1,,W1/2
K,2,T1,W1/2,
K,3,T1,T3/2,
K,4,W3-T2,T3/2,
K,5,W3-T2,W2/2,
K,6,W3,W2/2,
K,7,W3,-W2/2,
K,8,W3-T2,-W2/2,
K,9,W3-T2,-T3/2,
K,10,T1,-T3/2,
K,11,T1,-W1/2,
K,12,,-W1/2

A,1,2,11,12
A,3,4,9,10
A,5,6,7,8
WPCSYS
AGLUE,ALL

WPROTA,,90,
WPOFF,,,T3/2
ASBW,ALL
WPOFF,,,-T3
ASBW,ALL

ET,1,82
LESIZE,ALL,0.005     !网格尺寸
AMESH,ALL
SECWRITE,MYSECTION

!!!!!!!!!!!!!!!!!!!!!!
FINI
/CLE
/filn,  test_beam188, 1

/CONFIG,NRES,200000
/PREP7
LB=16.5
NE=330
NN=NE+1

DELTL=LB/NE

EM=2.0E11
DENG=7800
GRA=9.8

ET,1,BEAM188
MP,EX,1,EM
MP,NUXY,1,0.2
MP,DENS,1,DENG

SECTYPE,1,BEAM,MESH   
SECREAD,MYSECTION,,,MESH

*DO,I,1,331
N,I,(I-1)*DELTL
*ENDDO
TYPE,1   
MAT,1
ESYS,0
*DO,I,1,330
E,I,I+1
*ENDDO

/ESHAPE,1   
EPLOT  
!*   
FINISH  
!*  

/SOLU
ANTYPE,0
!*  
!施加约束
D,1,UX,,,,,UY
D,111,UY
D,221,UY
D,331,UY
ALLSEL

F,51,FY,-P
EPLOT   
!*  
FLST,2,2,1,ORDE,2   
FITEM,2,1   
FITEM,2,331
D,P51X, , , , , ,UZ, , , , ,
!*  
FLST,2,1,1,ORDE,1   
FITEM,2,1   
D,P51X, , , , , ,ROTX, , , , ,  
!*  
save
! /STATUS,SOLU
SOLVE   
FINISH  
save

» 本帖已获得的红花(最新10朵)

3楼2017-06-02 09:12:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

肉肉的了

新虫 (小有名气)

送红花一朵
引用回帖:
3楼: Originally posted by shgao20 at 2017-06-02 09:12:27
1  约束不足,缺少 z 向位移约束,以及绕量轴线的转动约束。
  2  施加载荷用的 P 没有赋过值。
修改后的命令流如下,可以试试看:
FINI
/CLE
/filn,  test_beam188, 1

/PREP7
W1=0.2
W2=0.2
W3=0.2
T1 ...

谢谢您,可以用了,只是不知道为什么要吧51号节点的uz和rotx约束住呢?

发自小木虫IOS客户端
4楼2017-06-02 21:17:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见