- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
风能利用技术-chapt1
中国并网风力发电行业发展历程 第一阶段:1986-1990探索和示范阶段。项目规模小、单机容量小。共建了4个风电厂,最大单机容量200kw. 第二阶段:1991-1995示范项目取得成效并逐步推广。共建了5个风电厂,最大单机容量500kw. 第三阶段:1996后,扩大建设规模阶段。最大单机容量1300kw. 到2002年底,全国共有32个风电场。 中国风能资源分布 广东省南澳:风力机132台,总装机容量达5.7万千瓦 达坂城风电场:风电机组157台,总容量8.28万千瓦 辉腾锡勒风电场:规划容量120万千瓦 在德国Buettstedt风电场内,共安装了22台 Enercon E-66(1800kw)和4台E-40(500kw),计26台风电机,41.6兆瓦的功率。 第一章 习题 1、画出南半球低气压和高气压梯度风的受力图,画出并解释低气压气旋方向及成因? 2、已知某处风速随时间的变化关系近似为:v=20*abs(sin(2*pi*t/10)),风力机启动风速10m/s,额定风速15m/s,顺桨风速18m/s,空气密度为1.2kg/m3:(a)计算0~10秒钟的平均风能密度?(b)在(w-t)图中画出有效可用风能的分布示意图? 3*、编写Matlab程序求解表1-3中的威布尔参数(用最小二乘法和风速均值方法)。 习题答案 1 南半球的梯度风如下图所示: 在中心区域是低压的情形下,受气压梯度力的作用,气流从圆周外向圆心运动,这种运动又受到地转偏向力的作用,使得气流向左手边发生偏转,因而形成逆时针方向的旋转。 v Fc Fg Fa 低压 v Fc Fa Fg 高压 低压 习题答案 2 习题答案 2 习题答案3 % 此程序用于求解威布尔分布中的两个参数来拟合风频分布。 % vf.dat数据格式为两列,第一列为离散的风速(递增),第二列为相应的风频数据。 clear; data=load(vf.dat);v=data(:,1);f=data(:,2);n=length(f); %-------------------------------------------------------------------- % ------采用最小二乘法A(a1)C(c1) p(1)=f(1); for i=2:n p(i)=f(i)+p(i-1); end for i=1:n x(i)=log(v(i)); y(i)=log(-log(1.-p(i))); end a=(sum(x.^2)*sum(y)-sum(x)*sum(x.*y))/(n*sum(x.^2)-sum(x)^2); b=(-sum(x)*sum(y)+n*sum(x.*y))/(n*sum(x.^2)-sum(x)^2); a1=exp(-a/b);c1=b; % ------采用风速均值和标准差A(a2)C(c2) mu=sum(v.*f); d=sum((v-mu).^2.*f); sigma=sqrt(d); c2=(sigma/mu)^(-1.086);a2=mu/0.889; %------------------------------------------------------- 习题答案3 %--绘图部分---------------------------------------------- clf; xx=[0:0.2:34]; yy=c1/a1*(xx./a1).^(c1-1).*exp(-(xx./a1).^c1); plot(xx,yy,r-,linewidth,3);hold on; yy=c2/a2*(xx./a2).^(c2-1).*exp(-(xx./a2).^c2); plot(xx,yy,b--,linewidth,3);hold on; plot([0,0,0.45,0.45],[0,f(1),f(1),0],k-,linewidth,2);hold on; dv=v(3)-v(2); for i=2:n plot([v(i)-dv/2.,v(i)-dv/2.,v(i)+dv/2.,v(i)+dv/2.],[0,f(i),f(i),0],... k-,linewidth,2);hold on; end hold off; axis([0,24,0,1.2*max(f)]); set(gca,fontsize,16,linewidth,2); set(gca,xtick,[0:2:24],fontsize,16); xlabel(\it v (m/s),fontsize,20); ylabel( 风频 ,fontsize,20); title(拟合的风频分布,font
原创力文档


文档评论(0)