第二章
2-1
两类平面问题的异同点。 解:
相同点:平面应力和平面应变都是起源于简化空间问题而设定的概念,把空间问题简化为近似的平面问题,可以减少分析和计算的工作量,并得到满足工程精度要求的结果。 不同点:
(1)平面应力问题考察的是薄板类问题,与平面垂直的方向应力可以忽略;而平面应变问题考察的是柱形体或水坝形体,与平面垂直方向的应变可以忽略。
(2)平面应力问题只有平行于xoy平面的三个应力分量,即σx、σy、τxy=τyx存在。同时,因为板很薄,这三个应力分量,以及分析问题时需考虑的应变分量和位移分量,都可以认为是不沿厚度变化的。即它们只是x和y的函数,不随z而变化。而平面应变问题只考虑平行于xoy平面的三个应变分量:x,y和xy。需要注意的是:由于z方向的伸缩被阻止,所以σz一般并不等于零。 (3)具有不同的物理方程:
平面应力问题的物理方程:
E()xy12Ey()yx 21Exyxy2(1)x 平面应变问题的物理方程:
12x(xy)E112y(yx)
E12(1)xyxyE可以看出,只需将平面应力问题的物理方程中的E换为理方程。
E,换为,就可以得到平面应变问题的物1122-2
要使下列应变分量成为一种可能的应变状态,试确定常数A0,A1,B0,B1,C0,C1,C2之间的关系。
xA0A1(x2y2)(x4y4)2244yB0B1(xy)(xy) 22xC0C1xy(xyC2)解:C14,A1B12C2
2xy22xy22 xyyx2xyC1xyx2C1xyy2C2C1xy xyxy3C1x23C1y2C2C12xA0A1(x2y2)(x4y4)22A12y 122yy2yB0B1(x2y2)(x4y4)2 2B12x1x2x23C1x23C1y2C2C12A112y22B112x2
3C1x212x23C1y212y2C2C12A12B1
→
C142C2A1B1
2-5
下图所示简支梁。上边界受均布荷载q作用。按材料力学方法算得。
q2(lx2)yAyBx2y2J
12qh2y2xy3()xCxBy2xh82x试问其是否满足平衡微分方程及应力边界条件,并求出y。 解:
q2q2q2(lx2)ylyxy2J2J2J
12qh2y212qh212qy2xy3()x3x3xh82h8h2x平衡微分方程
xxyX0xy
yyxY0yxxxyX0xyyyxY0yxq12qxy3xy0 Jhh3J12y12qh212qy2330yh8h2
y12qh212qy233yh8h212qh212qy2y33h8h2yy12qh12qy3q2qyyCyCh38h362hh3233
12qh212qy33q2qy3y3y3Cy3Ch8h62hh3qq yyhCq442qC2q3q2qy3q34y3yy3(y31)
22hh2hh
2-7
试求题下图所示线弹性梁的弯曲应变能。 (1)表示为杆端力偶MAB,MBA的函数;
(2)表示为杆端转角AB,BA的函数。
解:
答案分别为:
(1)纯弯梁的应变能和余应变能
1Ld2w2UEJ(2)dx
20dxB点平衡方程
xM(x)MBA(MBAMAB)
lxM(MM)BABAAB1LldxU20EJ2xx2222MMM(MM)BAABBABAAB1MBAl1Ll1L2dxldx002EJ2EJ2EJ2MABMBAl1(MM)2l1MBABAAB2EJ2EJ6EJ3MBA23MABMBA3MBA2(MM)2BAABl6EJ6EJ6EJ2BAl1M23MBA23MBA2MBA23MABMBA2MABMBAMAB2l6EJMBA2MABMBAMAB2l6EJ(2)
2w(x)M(x) x2EJM(x)x2x3w(x)MBA(MBAMAB)CxD
EJ26lxw(x)x0x2x3MBA(MBAMAB)CxDD0
26lx0l2l2w(x)xlMBA(MBAMAB)Cl026xl
2MBAlMABlC6w(x)MBA2MBAlMABlx2x3(MBAMAB)x 26l6(x)w(x) x(x)x0AB (x)xlBA(x)x02MBAlMABlx2MBAx(MBAMAB)2l6ABx02MBAlMABlAB62MBAlMABlx2(x)xlMBAx(MBAMAB)2l6BAxl
6MBAl3MABl3MBAl2MBAlMABlBA66MBAl2MBAl3MBAl3MABlMABlMBAl2MABlBA662MBAlMABl6ABMBAl2MABl6BA4MBAl2MABl12ABMBAl2MABl6BA4AB2BAl2MABl6BA
l4AB2BAMBAlM4BA2ABABlMBA2MABMBAMAB2U6EJl 4AB2BA24AB2BA4BA2AB4BA2AB26EJ4BA28BA216BA216AB284AB24AB216BAAB20BAAB16BAAB6EJ12BA212AB212BAAB6EJBA2AB2BAAB2EJ2-10
试用最小势能原理计算如图所示刚架,作M图。
解:
2lw2w1112EJ2dx1EJdxFpw1,xl 20022x1x2l22关于变分和求导的计算,有dyd(y),即变分符号和求导符号可以变换位置。 dxdx结合分部积分则有:
2lw2w2w1[EJ2dx1]EJ22dx0x0x2xllww2wwEJ20EJ2dxxx0xxxl2ww2w2wEJ2EJ2w0EJ22wdxxx0xxxx0l2ww2ww2w2w2wEJ2EJ2EJ2wEJ2w0EJ22wdxxxxlxxx0xxxxxxxlx0ll2l2由此可得:
l2w1w12w1w12w12w12w1EJEJEJEJw1EJ2w1dx1w1220x12x1xlx12x1x0x1x12xxxx1111x1lx10112w2w2EJx22x2xFpw1x1l2l2w2w2EJx22x2x22w2EJw22xx220x2l2w2EJw2x2x22x22w2EJw2dx2220xx220l根据边界条件:
w1x1w2x2w2x10w2x2x2l0x20
x200带入可得:
2w1w12w22w1w12w12w1EJEJ)(FpEJEJw1(EJ2)w12222x1x1xlx2x1x1x1x1x1x1x10x1lx1012w2EJw22x2x2l2w12w2EJ22w1dx1EJw2dx22200x1x1x2x2lx2l0
则有:
2w1(1)EJx12M1xl=01x1l2w22w1(2)(EJEJ)0x22x12x01即M1x0=M21x2l力矩连续2w1(3)FpEJ0x1x12xl1即Fp=M1x1x1l2w1(4)EJ0x1x12x01
即F1x00此处剪力为012w2(5)EJx2x22x即F2x2l02l0此处剪力为02w1(6)220x1x12w2(7)20x2x22最后得到:
M1Fp(lx1)M2constM1x0=M21x2lFpl
M2Fpl2-11
试用最小势能原理求下图所示刚架B点的水平位移。
解:
w(x)Ax3Bx2CxD
w(x)Ax3Bx2CxDw(x)x00w(x)0xx0w(x)0xxlw(x)xluw(x)3Al22Bl0xxl3AlB23w(x)xlAx3Alx2u23Al3Al3u22u A3l2u3uw(x)3x33lx2ll
3u22u3lx3x 3llL6u12uUEJ(3l3x)2dxFpu0llL6u12u6u12uUEJ(3l)2(3)2x22(3l3)xdxFpu0llll112u6u12u6uUEJ(3l)2l(3)2l3(3l3)l2Fpu
3llllw(x)111UEJ36u2348u2372u23FpulllEJU12u23FpulU12u2EJFpu3lEJU(24u3Fp)u
lFpl3u24EJ第三章
3-1
写出图示矩形薄板的边界条件,OA变为简支边,并作用有分布的弯矩M,BC边为固支边,AB为自由边。
OMaABaCXY
解:
对于BC边有:
xa0,0 xxa对于OC边有:
y020,20
yy022D(22)x0Mxy对于OA边有:
x00,Mxx02x00,Mx2x0D
对于AB边有:
Myb0,Mxyyb0,Qxyyb0
3-3
有一个矩形薄板,Axy(xa)(yb)是确定它对应怎样的几何,静力边界条件及怎样的载荷(A位移之常数)。 解:
由x00、xa0、y00、yb0
xAy(yb)(a)0、
x0xAy(yb)(a)0、
xayAx(xa)(b)0、
y0yAx(xa)(b)0
yb可知结构为简支边界条件 边界载荷为
Mxx022D222DAy(yb)
yx0xMxxa22D222DAy(yb)
yxax22D222DAx(xa)
xy0y22D222DAx(xa)
xybyMyy0Myyb板上有均布载荷
444PD422248DA
xyyx3-5
四边简支矩形薄板,边长为a和b,受横向分布载荷qq0sin该板的解,并求挠度、弯矩及反力。
证明: 由x0msinxasinyb,试证:挠度函数msinxasinyb是
0yaysin0,xamsinsin0 ababsin0xb0,yamsinsin0 baby0msinMyx0xa22xyxyD22Dm()2sinsinm()2sinsinxyaabbabx0x011xyD()2()22msinsin0babx0aMyxa22xyxyD22Dm()2sinsinm()2sinsinxyaabbabxaxa
11xyD()2()22msinsin0babxaaMxy022xyxyD22Dm()2sinsinm()2sinsinxyaabbaby0y0
11xyD()2()22msinsin0baby0aMxyb22xyxyD22Dm()2sinsinm()2sinsinxyaabbabybyb
11xyD()2()22msinsin0ababyb
故挠度函数满足简支边界条件
444q 又因为4222xxyy4D4m4m4mxy所以qD42224sin sinabbaba4mq04m4m所以D42224q0m 444aabbD42224abbaq0D42224abba444sinxasinyb
11q0()2()21xyxyb1a MyD()2()22msinsinsinsin222ababab24a2b2b4a11q0()2()21xyxyb1a MxD()2()22msinsinsinsin222bababa24a2b2b4a3-8
如图所示的矩形薄板,四边固支,受横向均匀载荷q0作用,现有以下8种试探位移函数:
1、Cmn(1cosmn(2m1)x(2m1)y)(1cos) ab2、C11(1cosxa)(1cosyb)
3、Cmn(1cosmnxa)(1cosyb)
4、Cmnsinmnmxny sinab225、C1(xa)(yb) 2222226、C1(xa)(yb)
7、Cmn(1cosmnxa)(1cosyb)
22228、C1(xa)(yb)xy
试问哪些位移函数可用于里兹法?哪些可用于伽辽金法?哪些函数这两种方法都适用?
只有位移边界条件:
xaya0 0xxayya(1)解:
xaCmn(1cosmn(2m1)x(2m1)y)(1cos)abxa(2m1)(2m1)yCmn(1cos)(1cos)01bmn
yaCmn(1cosmn(2m1)x(2m1)y)(1cos)abya(2m1)x(2m1)Cmn(1cos)(1cos)0a1mn
(2m1)(2m1)x(2n1)yCmn(sin)(1cos)0
xxamnaabxa(2m1)x(2n1)(2n1)yCmn(1cos)(sin)0
yyamnabbya两种方法都适用。 (2)解:
xyxaC11(1cos)(1cos)0abxaxyybC11(1cos)(1cos)0abyb
xyC11sin(1cos)0xxaaabxa
xyC11(1cos)sin0xybabbyb两种方法都适用。 (3)解:
xaCmn(1cosmnxa)(1cosyb)0
两种方法都不适用。 (4)解:
xaCmnsinybmxnysin0abmn
mxnyCmnsinsin0abmnmmxnycossin0aab
nmxnysincos0babCmnxxamnCmnyybmn两种方法都不适用。 (5)解:
xaC1(xa)2(yb)20
两种方法都不适用。 (6)解:
xaC1(x2a2)2(y2b2)20ybC1(x2a2)2(y2b2)20
C12x(x2a2)(y2b2)20xxa
C1(x2a2)22y(y2b2)0xyb两种方法都适用。 (7)解:
xaCmn(1cosmnxybab
xyCmn(1cos)(1cos)0abmn)(1cosy)0xyCmnsin(1cos)0xxamnaab
xyCmn(1cos)sin)0yybmnabb两种方法都适用。 (8)解:
xaC1(x2a2)(y2b2)xy0ybC1(xa)(yb)xy0
2222
C12x(y2b2)xyC1(x2a2)(y2b2)y0xxa
C12x(y2b2)xyC1(x2a2)(y2b2)y0yyb两种方法都不适用。
第四章
4-1
(1)
K=0 (2)
K=4 (3)
K=0 (4)
瞬变
4-2
求下图空间桁架结构各杆内力。
7581623P34
P2零力杆有:
N510 N120
N130 N840 N340 N420
7583P3
去掉零力杆后结构形式如图62P2在节点处268平面内,在68方向上没有外载荷,仅有28杆在此方向上有分量,所以:N820
7583P3(结构传力路线)
进一步去掉零力杆后结构形式如图62P27由节点3处在378平面内受力图:3P3cos8可知:
N38P3cos cos1N37P3cossin
cos7由节点3处在357平面内受力图:3P3sin5可知:
N35P3sin cos2N37P3sinsin
cossincossinP3sin coscos5N37P3由节点2处在256平面内受力图:62可知:
P2P2 cosP2sin cosN25N26第五章
5-1
算下图在Mx作用下正应力薄板不受正应力
Y6.01Y’0.13.21525XX’155.02.540154解: 1求形心
3
6.0153.255.0152.5156.03.22.55.0
90167537.50.38916.73.2402.540x06.03.22.55.0
22813.65316.7y0Y6.01Y’0.13.21525XX’155.02.5401542求惯性矩
如图建立坐标轴,以形心为原点
3
Ix615.38923.25.38925.014.61122.514.61123114.97
Iy6(13.653)25.0(13.653)23.226.34722.526.34726007.19
Ixy6(13.653)15.3893.226.3475.3895.0(13.653)(14.611)2.526.347(14.611)770.05MxMyMx1I2xy3求等效弯矩
IxyIx1.0327Mx
IxIyIxyIy0.2553Mx
MyMxMy14求正应力
2IxyIxIyzMyIyxMxyIxMxMxx1.0327y
6007.193114.974.250105xMx3.315104Mxy0.2553z14.250105(13.563)Mx3.315104Mx(15.389)45.250104Mx
z24.250105(26.437)Mx3.315104Mx(5.389)29.100104Mx z14.250105(26.437)Mx3.315104Mx(14.611)37.200104Mx z44.250105(13.563)Mx3.315104Mx(14.611)54.200104Mx 5-2
计算下图在My下正应力
10.25Y2MyX25O0.125450
3解:
1形心位置
由于机构的对称性,形心在结构几何中心,即坐标系原点 2求惯性矩
0.25503Iy250.12525229114.6
1223求正应力
zMyIyxMxyIxMy9114.6
x1.097104Myx5-3
图示薄壁梁剖面,它的壁板不能承受正应力。四个缘条的面积标注在图上,求剖面的弯心位置。
求中心主轴。
xoy坐标轴不是中心主轴,y0轴是正应力面积的对称轴,所以是一个中心主轴,剖面所有缘条面积的形心在c点,c点到圆心距离为:
所以,另一中心主轴x0轴很容易在图中画出。 剖面对中心主轴的惯性矩为:
各壁板的静矩为:
计算弯心位置,可选o点为力矩中心,根据扭矩平衡可得:
又因为:
计算弯心位置,可选o点为力矩中心,根据扭矩平衡可得:
第六章
6-1
求图示薄壁剖面在通过剖面弯心的剪力Qy作用下的剪流分布,并求剖面的弯心位置。设壁板承受正应力,壁板厚度为t。
解:
Jx2B0B2B111323stds2stds2sttB 0226302SxqB02212212BststBStSt Bstds222022SSxQy Jx212BStSt22BSSQ 3222Qqyy3134BtB3qmax32Qy 4B
弯心在o处 6-2
求图示剖面的弯心位置,设壁板不承受正应力,缘条面积为f。
解:
Jx2102f202f1000f
2S120f xSx2330f
qSxQy Jxq1220f1QyQy 1000f5030f3QyQy
1000f100q12
Qyxdqds
11Qyx2Qy105208Qy 550在x= - 8处 6-3
剖面的壁板不承受正应力,缘条面积为f =2cm,Q x=Q y=10000N。
2
1求惯性矩
Jx2302f302f3600f Jy2302f1800f
2求静矩
122S130fS0 ,xy23Sx2360f,Sy30f
3求开剖面剪流
qSySxQyQx JxJy30f01QyQyQy 3600f1800f12060f30f11QyQxQyQx
3600f1800f6060q12q234利用力矩平衡求q0
Qyxodqdsdq0ds
111Qyxodqdsdq0ds2Qy3030Qy3060Qx3060dq0ds75Qydq0ds10Qy601206065Qydq0ds65Qydq0ds3030q030304q09004q0q065Qy900413Qy1804
q1213Qy1Qy -17.818120180413Qy11 QyQx232.181860601804q235扭角
单位长度上的扭角(也称相对扭角)θ为:
1qpdsGt
∮只沿闭室周线积分,为面积的2倍。
1qdspGt 11100001313113**30*60*60900(4)1804120180460601804a=(10000/900/7.14)*(-13/180/7.14*3.14*30+(1/120-13/180/7.14)*60+(1/60+1/60-13/180/7.14)*60) =0.5187 6弯心
SxGtds1Xsxds
dsJxsGts沿所有壁板积分,∮只沿闭室周线积分,为面积的2倍。
SxGtds1XsxdsdsJxs Gt130f*30*260f*6030f*30*30*260f*60*30900(4)720030120=0 6-4
剖面的壁板承受正应力,壁板厚度t=0.1cm,Q y=10000N。
231451求惯性矩
6
Jx2*909003RsinRtd2R2Bt2*r2tdr0222RR32*Rtsind2RBt2*t03901cos2R3322*Rtd2RBt2*t0223sin2R122*R3t2RBt2*t403223
R32*Rt2RBt2*t432R3t2R2Bt23322求静矩
2S1RsinRtdR2tsindR2tcosdR2tcos1 x00032Sxs1Rt
2SxR2tBRt
25SxR2tBRtrtdrR2tBRtr1R122Rr1 23求开剖面剪流
qSxQy Jxq12xR2tcos1Qy 232Rt2RBt23s1RtQy
232Rt2RBt23R2tBRt32qx122Rr1252qxQy
232Rt2RBt234利用力矩平衡求q0 Qyxodqdsdq0ds2BRtcos1s1Rt 22QyRdQyRds1q0R2000232322Rt2RBtRt2RBt23232Rtcos1s1Rt22QRdsQRdy1y02302322Rt2RBtRt2RBt2323q0R22BQyR2tRtcosR2t2222sdsQRdQRd11yy000222333222Rt2RBtRt2RBtRt2RBt232323R22BQyR2tQyR2tQRty2222sdsRcosdRd11000222333222Rt2RBtRt2RBtRt2RBt232323R2QyR2tBR2220s1ds1R0cosd2232Rt2RBt23R2QyR2tB2R22R6QytB22R22232222Rt2RBt6B2RQ234Rt12Bt23y2R2R2234R12BRB5扭角
单位长度上的扭角(也称相对扭角)θ为:
1
qpdsGt
∮只沿闭室周线积分,为面积的2倍。
1qdspGt2226B2RRtcos122QyQyRd202R3t2R2Bt234R12BR23212222RtBRtRr16B2RR22QyQydr1202234R12BRR3t2R2Bt1232R(4)Gt12222RtBRtRr1RRtcos1222QyRd2Qydr100232322Rt2RBtRt2RBt23231R2(4)6B22R2R2234R12BRQy26B22R2R2234R12BRGtQy
220cosRtQyRd20232Rt2RBt23R3tQy2R2tBRtR122Rr12Qydr1232Rt2RBt23Qy26B22R2R2234R12BR1R2(4)6B22R2R2234R12BR232Rt2RBt23QyGt20Rr12cosR2t2QyRd2Qydr10223232Rt2RBtRt2RBt2323R22RtBRt6B22R2R6B22R2RR3tQy22RQyQy2Qy22234R12BR234R12BR232322Rt2RBtRt2RBt123232R(4)GtR2tQyRR322Qy2323223Rt2RBtRt2RBt2323 R2222223RtBRt6B2RR6B2RRRtQy22RQyQy2Qy22234R12BR234R12BR232322Rt2RBtRt2RBt123232R(4)Gt2R23326RRtBRtQy3RtQy2RQy6RtQyR6212B22R2R2Qy2234R12BR233Rt2R2Bt1232R(4)Gt2226Rt6BtR3Rt6tR612B2R2234R12BR3Rt2Bt1232QyR(4)Gt22222234R12BR6BtR3Rt3Rt2Bt23612B2RQy23Rt2Bt234R212BRGtR2(4)23
6弯心
1XsxdsJxsSxGtds dsGts沿所有壁板积分,∮只沿闭室周线积分,为面积的2倍。
1Xssxds232Rt2RBt23SxGtds dsGtXf(R,B,t,Qy,G,)
有限元
例9-1
平面桁架系统由四个节点和六个杆单元组成。以序号○1至○4作为节点标号,序号(1)至(6)作为单元标号。每个节点具有分别沿x轴和y轴方向上的两个节点位移和对应的节点力,分别用序号1至8作为标识。则节点位移和节点力的量可以分别表示为:
U[U1U2U3U4U5U6U7U8]T S[S1S2S3S4S5S6S7
S8]T
1、结构的基本信息: 1)节点信息
节点编号 坐标x ( cm ) 坐标y(cm) ○1 10 10 ○2 10 0 ○3 0 10 ○4 0 0 2) 单元信息
单元编号 节点连接关系 截面积(cm2) 弹性模量(kg/cm2) (1) 1-3 10 2105 (2) 1-4 22 2105 (3) 1-2 5 2105 (4) 2-3 22 2105 (5) 2-4 10 2105 (6) 3-4 5 2105 3)载荷约束条件 外力载荷:
P[P1P2P3P4P5P6P7P8]T[200400000000]T边界条件:
U[U1U2U3U4U5U6U7UT8][U1U2U3U40002、整体序号与局部序号
0]T
U(e)[U(2i1)U(2i)U(2j1)U(2j)]T S(e)[S(2i1)S(2i)S(2j1)S(2j)]T
3、结构的总体刚度矩阵
(1)S(1)2(1)S(2)11050(1)S(5)2(1)S0(6)020002000U(1)U0(2) 0U(5)U0(6)(2)S(1)4444U(1)(2)4444US(2)1105(2) (2)S(7)4444U(7)(2)US4444(8)(8)(3)S(1)00(3)01S(2)5110(3)S(3)00(3)S(4)0100U(1)U01(2) 00U(3)U01(4)(4)S(3)4444U(3)(4)4444US(4)1105(4) (4)S(5)4444U(5)(4)US4444(6)(6)(5)S(3)2(5)S(4)11050(5)S(7)2(5)S0(8)020002000U(3)U0(4) 0U(7)U0(8)(6)S(5)00(6)01S(6)5110(6)S(7)00(6)S01(8)00U(5)U01(6) 00U(7)U01(8)R6R7R8
TKU20040000R5645对006称0145 K1105204460044454420006440001454、边界条件和方程的求解
UAU(1)U(2)U(3)U(4)
TUBUUUU[0000] (5)(6)(7)(8)1UAKAPA0.00060.00140.00040.0006
TTTPTBKABUA20080400320
以杆(1)为例:
(1)S(1)2020S(1)0000.0006(2)S(1)(5)1105020200.0014120S(1)(6)000000(1)EBU12(kg/cm2)
N(1)(1)F(1)120(kg)
01200T
因篇幅问题不能全部显示,请点此查看更多更全内容