模型說明：

使用ANSYS軟件，模擬一個截面為0.2*0.3m的矩形梁，長為4m，然后再約1/3處切開為兩個體然后再進行體的粘結(jié)，對小部分的底面進行全約束，在梁端施加一個2kN的集中荷載，然后求解并求內(nèi)力，然后繪制內(nèi)力圖。

結(jié)果說明：

在約束處有一個力的突變，但是力的大小卻遠大于施加的力，不符合理論分析結(jié)果，相差甚遠。

分析命令流如下：

finish

/clear

/prep7

*CREAT,ZSuNQM,mac!

SUCR,suna,cplane,3!

SUSEL,s,suna!

!!!!!映射截面正應(yīng)力和剪應(yīng)力

SUMAP,susx,s,x!

SUMAP,susy,s,y!

SUMAP,susz,s,z!

SUMAP,susxy,s,xy!

SUMAP,susyz,s,yz!

SUMAP,susxz,s,xz!

!!!!!對面積重心取矩

SUEVAL,Aarea,da,sum!

SUEVAL,sscx,gcx,intg!

Sucx=sscx/Aarea!

SUEVAL,sscy,gcy,intg!

Sucy=sscy/Aarea!

SUEVAL,sscz,gcz,intg!

Sucz=sscz/Aarea!

!!!!!積分求截面軸力、剪力、彎矩和扭矩

SUEVAL,suNz,susz,intg!

SUEVAL,suQy,susyz,intg!

SUEVAL,suQx,susxz,intg!

SUCALC,sxgcx,susz,mult,gcx!

SUEVAL,MZX,sxgcx,intg!

SUCALC,sxgcy,susz,mult,gcy!

SUEVAL,MZY,sxgcy,intg!

SUCALC,sxgcx,susxz,mult,gcy!

SUEVAL,TZX,sxgcx,intg!

SUCALC,sxgcy,susyz,mult,gcx!

SUEVAL,TZY,sxgcy,intg!

!!!!!將彎矩和扭矩簡化到面積中心上

MZX=MZX-SuNz*sucx!

MZY=MZY-SuNz*sucy!

TZ=-TZX+TZY-suQy*sucx+suQx*sucy!

*END

et,1,solid95

mp,ex,1,2e11

mp,prxy,1,0.3

blc4,2,3,0.2,0.3,4

wpoffs,,,1.3

vsbw,all

vglue,all

da,10,all

esize,0.05

vmesh,all

n1=node(2.1,3.15,4)

f,n1,fy,-2e3

finish

/solu

solve

finish

/post1

*dim,ny1,,6,41

*dim,y1,,41,6

*do,i,1,41

csys,0

wpcsys,-1

wpoffs,,,(i-1)*0.1

ZSuNQM

ny1(1,i)=suNz,suQx,suQy,MZX,MZY,TZ

*enddo

*do,i,1,6

*do,j,1,41

y1(j,i)=ny1(i,j)

*enddo

*enddo

*create,sc,mac

*cfopen,4_5,txt

*vwrite,sequ,y1(1,1),y1(1,2),y1(1,3),y1(1,4),y1(1,5),y1(1,6)

(F3.0,F10.0,F10.0,F10.0,F10.0,F10.0,F10.0)

*cfclos

*end

sc

模型圖如下：

圖一：梁模型圖

圖二：約束情況（對圖示面施加一個全約束）

圖三：計算應(yīng)力云圖結(jié)果

得到的內(nèi)力圖如下：

圖四：軸力圖（未施加軸力，但在約束處軸力卻為16kN，遠大于梁端施加的荷載2kN）

圖五：剪力圖（在自由端剪力為施加的外力值，與理論分析相符，但是在約束處卻反向突變?yōu)?2kN）

圖六：彎矩圖（與理論分析比較接近）