学号：s0姓名：刘翔专业：结构工程。
pc构件截面分析：
对一pc构件截面进行全过程分析并得出其m-φ-n曲线。
算例：某预应力钢筋混凝土梁跨径l=6m；梁长分段m=32；截面分段数n=100；梁宽b=0.4m；梁高h=0.
6m； 预应力钢筋面积asp=139.0×10-6m2； 预应力钢筋坐标yp=0.25m；下缘钢筋坐标ys=0.
25m； 上缘钢筋坐标yss=－0.25m；下缘钢筋面积as=628.0×10-6m2； 上缘钢筋面积ass=628.
0×10-6m2；ess=3.25×104n/mm2；混凝土峰值应变staco=0.002；混凝土弹性模量 ec=3.
25×104n/mm2；混凝土极限压应变stacu=0.0033；混凝土峰值点应力fc=32.4×106pa；下缘钢筋屈服强度fs=335.
0×106pa； 下缘钢筋弹性模量es=2.1×105n/mm2；下缘钢筋屈服应变stasso=0.0017；下缘钢筋极限应变stassu=0.
16；上缘钢筋弹性模量ess=2.1×105n/mm2；上缘钢筋屈服强度fss=335.0×106pa；上缘钢筋屈服应变stasso=0.
0017；上缘钢筋极限应变stassu=0.16；预应力钢筋弹性模量esp=2.1×105n/mm2；预应力钢筋条件屈服强度fsp=1581.
0×106pa；预应力钢筋条件屈服应变stasp=0.0017；预应力钢筋极限屈服强度fsup=1860.0×106pa；预应力钢筋极限应变stasup=0.
0275。
解：理论分析：
1） 基本假定：
1 平截面假定，即梁正截面变形后仍保持平面，截面应变为直线分布，不考虑钢筋与混凝土之间的相对位移。
2 不考虑剪切变形的影响，且截面开裂后忽略受拉区混凝土的抗拉作用。
混凝土受压区顶层钎维应变达到0.003时，即认为达到极限状态。
各材料的本构关系分别满足以下条件：
混凝土本构关系采用：
各参数取值见《混凝土结构设计规范》（gb50010-2002）规定）
普通钢筋本构关系采用：
预应力钢筋本构关系采用：
其中——对应于预应力筋应变为时的预应力筋应力；
——预应力筋屈服极限处的应力值；
——对应于预应力筋应力为时的预应力筋应变；
——预应力筋强度极限处的应力值；
——对应于预应力筋强度极限处的拉应变；
2） 全截面各阶段应变分析。
将混凝土的截面分成有限条带，在外荷载的作用下假定每一条带上的应力均匀分布，全截面各阶段应变为：
全截面消压前的截面应变
在预应力的作用下梁截面的应变见图1。假定拉应变为—，压应变为+。
预应力钢筋在有效预应力（扣除全部预应力损失后）作用下的拉应变为：
其中“－”表示为拉应变。
此时各条带混凝土的压应变为：
为任意一条带距截面中心轴的距离，位于受拉区取负值，受压区取正值。
此时预应力钢筋重心水平处混凝土的压应变为：
全截面消压时的截面应变
全截面消压时预应力钢筋的应变为：
此时受拉钢筋的应变为：
其中“－”表示为拉应变。
全截面消压后的截面应变，图2为消压以后矩形截面应变、应力分布图。
应力计算：由图2可知。
式中———截面受压区边缘混凝土应变；——受拉钢筋应变；
———截面有效高度；——截面曲率；
截面上任意一条条带的应变为，为任意一条带距截面中心轴的距离，位于受拉区取负值，受压区取正值。
为中性轴z—z的应变。
根据所采用的本构关系可得：
各条带混凝土的应力
压区非预应力钢筋的应力
拉区非预应力钢筋的应力
预应力钢筋的应力
每一条带上的混凝土、受压钢筋、受拉钢筋及预应力钢筋的力分别为：
得到平衡方程。
3）截面弯矩—曲率计算分析过程：
取曲率；假定梁截面受压区边缘混凝土应变；
求各混凝土条带和钢筋的应变；
按混凝土和钢筋的应力-应变关系求与应变相对应的应力值；
按式判别是否满足平衡条件；
若不满足平衡条件，则需调整应变值，重复步骤③-⑤
满足平衡条件后，按式。
求得所对应的弯矩；
循环步骤①-⑦直到得出整个的关系。
4）荷载-挠度曲线的计算。
在求得构件各截面的曲率以后，将梁分成m个小段，相应结点有m+1个，结点之内的每一个小段内的曲率假定是线形变化的。运用共扼梁法求的任意截面的转角和挠度（i=0，1，…，m）。
图3（b）为简支梁（图3（a））所对应的虚梁，根据虚荷载平衡条件。
跨中位移为。
图3 简支梁与对应的虚梁
采用分级加曲率的方法，求p-的关系曲线。以跨中弯矩最大截面m处的曲率作为控制值，逐级增加曲率。计算框图见图4.
每一级由确定跨中弯矩；
由计算由p计算各截面；
由计算；由求得挠度，重复步骤①~⑥直到跨中位移达到所预定的值。
利用fortran90语言，编制分析程序如下：
cmain programming
real m(300),c(300),p(300),d(300),mom(100),coc(100)
integer i
real fy,es,esh
real fc,fct
real as1,as
real l,h,b
real ao(2)
integer sn,ln,st
cccccccccccc构件参数 ln 是半跨划分的总段数。
real dc,de,ee,em
real sf1,sf2,dsf
real ffc,fs
real mi,mic,mis
real z,e,s,r
real eo,eu,cc
real lp,hh,hn,aas,etop
real esy,da
real dd,dsn,dl
real mm,mo,co,dp,tan
integer j,k,n,ii,jj
integer jmax1,jmax2,jmax3
do 100,i=1,300
m(i)=0.0
c(i)=0.0
d(i)=0.0
p(i)=0.0
100 continue
do 110,i=1,100
mom(i)=0.0
coc(i)=0.0
110 continue
open(1,file=''status='old')
read(1,*)fy,es,esh
read(1,*)fc,fct
read(1,*)as1,as
read(1,*)l,h,b
read(1,*)ao(1),ao(2)
read(1,*)sn,ln,st
close(1)
esy=fy/es
hn=h/sn
eu=-0.004
eo=-0.002
ee=-0.0001
dc=0.2e-6
de=0.5e-4
sf1=0.0
sf2=0.0
ffc=0.0
fs=0.0
mi=0.0
cc=0.0
jmax1=500
jmax2=0
jmax3=0
j=1ii=0
120 j=j+1
em=dec(j)=c(j-1)+dc
130 ii=ii+1
mic=0.0
mis=0.0
ffc=0.0
do 121,i=1,sn
z=h/2.0-(i-1)*hn-hn/2.0
e=ee+z*c(j)
if( then
s=0.0else if ( then
s=fctelse if ( then
s=(2*fct*e)/(e+0.0001)
else if ( then
s=-0.85*fc*(2*e/eo-e/eo*e/eo)
else if ( then
s=-0.85*fc*(1-100.0*(eo-e))
elses=0.0
end if
ffc=ffc+s*b*(h/sn)
mic=mic+s*b*z*(h/sn)
第一章工程结构非线性分析概述
姓名 朱伟华学号 14102020278 作业1 已知薄板l 30m，1.32m，t 0.02m，e 2.06 105mpa，7.8 103kg m3。求自重作用下的变形。要求用ansys求解，并写出命令流。解答1 按薄板进行计算。finish clearprep7 et,1,shell63 r,1...
混凝土强度理论与结构非线性分析 思考题2019
1.造成钢筋混凝土结构呈现非线性变形特性的原因有哪些？2.利用有限元法进行钢筋混凝土结构的非线性分析有什么优点？3.获取截面弯矩 曲率关系曲线的方法有哪些？4.简述给定轴力情况下用数值积分方法计算任意形状截面的弯矩 曲率全过程。曲线的主要步骤。5.已知截面的弯矩m和轴力n求截面曲率时，怎样避免简单的...
材料非线性有限元法
第四章材料非线性有限元法。以上三章分别研究了线性弹性有限元法，材料非线性本构方程和非线性方程组解法，本章就可以研究材料非线性有限元法了。在材料非线性基本方程中，除第二章所述的本构方程外，与线性弹性一样，而非线性有限元法又归结为一系列线性弹性问题。因此，只要在第一章中改用第二章的本构方程，就可建立材料...