m，外半径15mm，长度200mm，矩10E6 Nmm，mass21单元实常数取零。(1)梁建模划网(2)在梁轴线上生成1节点，偏离端面10mm。(3)选定mass21单元属性，通过第二步的节点直接生成单元(E或EN)。(4)指定刚性化区域(CERIG)，先选取mass21对应的节点，再选取梁端面节点。(5)加约束和载荷，梁另一端全约束，在mass21对应节点上加扭矩。(6)求解后验证结果。(7)在mass21上加弯距，求解验证结果。加扭矩，按公式计算得到最大剪应力235.179Mpa。应力分布连续，无应力集中现象。加弯距，按公式计算轴向应力为470.357Mpa，但在应力云图上梁端部明显出现端部效应，取中间截面轴向应力对比，结果基本吻合。所以加弯距时考虑结构加长以减小端部效应的影响。

&lt xmlnamespace prefix ="v" ns ="urn:schemas-microsoft-com:vml" /&gt &lt xmlnamespace prefix ="o" ns ="urn:schemas-microsoft-com:office:office" /&gt

关于实体单元施加弯矩的方法一、施加方法思路1：矩或扭矩说白了就是矩，所谓矩就是力和力臂的乘积。 施加矩可以等效为施加力；思路2：直接施加弯矩或扭矩，此时需要引入一个具有旋转自由度的节点；二、在ANSYS中实现的方法这里说说3个基本方法，当然可以使用这3个方法的组合方法，组合方法就是对3个基本方法的延伸，但原理仍不变。方法1：引入mass21，利用cerig命令Ex1:

/prep7

block,0,1,0,1,0,2

k,9,0.5,0.5,2.5

mp,ex,1,2e10

mp,prxy,1,0.2

mp,prxy,1,0.3

r,2,1e-6

et,1,45

et,2,21

keyopt,2,3,0

lesize,all,0.2

vmesh,all

ksel,s,,,9

type,2

real,2

kmesh,all

allsel

nsel,s,loc,z,2,3

NPLOT

CERIG,node(0.5,0.5,2.5),ALL,ALL, , , ,

allsel

/SOLU

f,node(0.5,0.5,2.5),my,100e3

FINISH

/SOL

nsel,s,loc,z,0

d,all,all

allsel

solve

方法2：利用mpc184单元/prep7

block,0,1,0,1,0,2

mp,ex,1,2e10

mp,prxy,1,0.2

mp,prxy,1,0.3

et,1,45

et,2,184

keyopt,2,1,1

lesize,all,0.2

vmesh,all

n,1000,0.5,0.5,2.5

type,2

mat,2

*do,i,1,36

e,1000,36+i

*enddo

allsel

allsel

/SOLU

f,node(0.5,0.5,2.5),my,100e3

FINISH

/SOL

nsel,s,loc,z,0

d,all,all

allsel

solve

方法3：使用rbe3命令/prep7

block,0,1,0,1,0,2

k,9,0.5,0.5,2.5

mp,ex,1,2e10

mp,prxy,1,0.2

mp,prxy,1,0.3

r,2,1e-6

et,1,45

et,2,21

keyopt,2,3,0

lesize,all,0.2

vmesh,all

ksel,s,,,9

type,2

real,2

kmesh,all

allsel

*di***a,array,36

*do,i,1,36

sla(i)=i+36

*enddo

*di***a2,array,36

*do,i,1,36

sla2(i)=i+36

*enddo

allsel

rbe3,node(0.5,0.5,2.5),all,sla,sla2

allsel

/SOLU

f,node(0.5,0.5,2.5),my,100e3

FINISH

/SOL

nsel,s,loc,z,0

d,all,all

allsel