msp10.m
%1:entrefer, 2:culasse stator, 3: fer rotor, %4?7: conducteurs stator, 28 et 29 conducteurs rotor entref=1e-3; %entrefer
rs=39.385e-3; %rayon int閞ieur stator rc=71.75e-3; %rayon culasse (ext閞ieur stator) rr=rs-entref; %rayon rotor nb=140; %nbre spires par phase long=0.125; %longueur utile machine p=1; f=50;
beta=35*pi/180; beta=beta/p; %culasse rotor beta1=50*pi/180; beta1=beta1/p; % 閜anouissement polaire dessinp1 a=0;
gr=groto(gr,(aa)); g=[gsf gr]; limites [p,e,t]=initmesh(g);
%[p,e,t]=refinemesh(g,p,e,t,[1 30]); %pdegplot(g); axis equal %hold on ; pdemesh(p,e,t) %pause % d閒inition des constantes kexc=0;
sigmaexc=kexc*5e6;
nuo=1/(4e-7*pi); murs=500; murr=500; ks=1; kr=0;
jex1=1e6*ks; jex11=-jex1; jex2=-jex1/2; jex21=-jex2; jex3=-jex1/2; jex31=-jex3; jexr1=1e6*kr; jexr2=-jexr1; garnissagep1 % assemblage des matrices u=assempde(bond,p,e,t,nu,0,j); limites.m ng=size(g,2);
ng4=ng-5; bb10=[1 1 1 1 1 1 48 48 49 48]; bb00=[0 1 1 1 1 1 48 48 49 48]; bb1=bb10; bb0=bb00; for ib=1:3
bb1=[bb10 bb1];
end
for ib=1:ng4; bb0=[bb00 bb0]; end
bond=[bb1 bb0]; groto.m
function [gb]=groto(ga,teta) % cette fonction renvoie la matrice de g閛m閠rie gb qui repr閟ente la %
rotation dangle teta et de centre o(0,0) de celle pass閑 en param鑤re csr=cos(teta); snr=sin(teta); grot=sparse(... [ 1 0 0 0 0 0 0 0 0 0 ;...
0 csr 0 -snr 0 0 0 0 0 0 ;... 0 0 csr 0 -snr 0 0 0 0 0 ;... 0 snr 0 csr 0 0 0 0 0 0 ;... 0 0 snr 0 csr 0 0 0 0 0 ;... 0 0 0 0 0 1 0 0 0 0 ;... 0 0 0 0 0 0 1 0 0 0 ;...
0 0 0 0 0 0 0 csr -snr 0 ;... 0 0 0 0 0 0 0 snr csr 0 ;... 0 0 0 0 0 0 0 0 0
1 ]); gb=grot*ga; couple.m tic; clear;clc; close all; aa=22.5*pi/180; a=aa+pi/2; pas=0.2
nom=3*pi/pas nom=round(nom)+1 nnn=1:nom
ener=zeros(1,nom) for s=1:nom msp10 aa=aa+pas energie1
ener(1,s)=energ_total end
coup(1,1)=0 for s=2:nom
coup(1,s)=ener(1,s)-ener(1,s-1) coup(1,s)=coup(1,s)/pas
end subplot(2,1,1),plot(nnn,ener,r),title(energie),ylabel(energie magnetique
stockee),xlabel(k),hold on
subplot(2,1,2),plot(nnn,coup,b),title(couple),ylabel(couple),xlabel(k) toc garnissagep1.m npt=size(p,2); ntrg=size(t,2); k=sparse(npt,npt); b=sparse(npt,1); ke=zeros(3);me=zeros(3); % d閒inition de j j=sparse(1,ntrg);
indice1=find((t(4,:)>=4) & (t(4,:)<=7)); indice11=find((t(4,:)>=16) & (t(4,:)<=19)); indice2=find((t(4,:)>=12) & (t(4,:)<=15)); indice21=find((t(4,:)>=24) & (t(4,:)<=27)); indice3=find((t(4,:)>=20) & (t(4,:)<=23)); indice31=find((t(4,:)>=8) & (t(4,:)<=11)); indicer1=find(t(4,:)==28); indicer2=find(t(4,:)==29); ij=1; j(indice1)=jex1*ij*ones(1,length(indice1)); j(indice11)=jex11*ij*ones(1,length(indice11)); j(indice2)=jex2*ij*ones(1,length(indice2)); j(indice21)=jex21*ij*ones(1,length(indice21)); j(indice3)=jex3*ij*ones(1,length(indice3)); j(indice31)=jex31*ij*ones(1,length(indice31)); j(indicer1)=jexr1*ij*ones(1,length(indicer1)); j(indicer2)=-jexr1*ij*ones(1,length(indicer2)); % d閒inition de nu nu=nuo*ones(1,ntrg); %rotor mur=1000;
indice=find(t(4,:)==3);
nu(indice)=nuo/murr*ones(1,length(indice)); %stator mur =1000; indice=find(t(4,:)==2);
nu(indice)=nuo/murs*ones(1,length(indice)); sigma=sparse(1,ntrg);
% indice=find((t(4,:)==28)|(t(4,:)==29)); % sigma(indice)=sigmaexc*ones(size(indice)); energie1.m êlcul de lenergie [ux,uy]=pdegrad(p,t,u); bx=uy; by=-ux;
for i=1:length(bx)
b(i)=(bx(i)^2+by(i)^2)^.5; end
nu_e=1/(4e-7*pi);%inverse de la permeabilite de lair nu_c=nu_e; nu_s=nu_e/500*ones(1,ntrg);%inverse de la permeabilite du fer
nu_r=nu_e/500*ones(1,ntrg); h=sparse(1,ntrg);
ind_e=find(t(4,:)==2); ind_s=find(t(4,:)==1); ind_r=find(t(4,:)==3);
ind_c=find((t(4,:)>=4)&(t(4,:)<=27)); bstr=b(ind_s)%matrice dinduction stator transpos閑(ligne) brtr=b(ind_r)%matrice dinduction rotor transpos閑(ligne) betr=b(ind_e)%matrice dinduction entrefer transpos閑(ligne) bctr=b(ind_c)%matrice dinduction conducteur
transposee(ligne) h(ind_e)=nu_e*b(ind_e); h(ind_s)=nu_s(ind_s).*bstr; h(ind_r)=nu_r(ind_r).*brtr; h(ind_c)=nu_c*b(ind_c); aire=pdetrg(p,t); vol_trg=long*aire;
energ_elem_e=0.5*h(ind_e).*betr.*vol_trg(ind_e) energ_tot_e=sum(energ_elem_e); energ_elem_s=0.5*h(ind_s).*bstr.*vol_trg(ind_s) energ_tot_s=sum(energ_elem_s); energ_elem_r=0.5*h(ind_r).*brtr.*vol_trg(ind_r) energ_tot_r=sum(energ_elem_r); energ_elem_c=0.5*h(ind_c).*bctr.*vol_trg(ind_c)篇二：关于srtp 关于srtp
一、什么是srtp 1、大学生科研训练计划（student research training program,简称srtp）是专为在校本科生设计的一种科研项目资助计划。，是针对在校本科生开展的科学研究训练项目，是在本科教育阶段实施实践教学改革的一项措施。 在国外，有美国麻省理工学院实行的 urop 计划。 urop ，意为 undergraduate research opportunities program ，即给本科生一个参与科学研究训练的“机会”， mit 已经有 70% 的本科生在大学四年期间参加过这个活动。 在国内， 清华大学首先借鉴 urop 计划，从 1996 年开始创建并实施的“ srt ”计划。 目前已经有众多国内高校开展了“ srt ”计划, srtp 采取项目化的运作模式，通过设立创新基金和本科生自主申报的方式确定立项并给予以资金支持，鼓励学生在导师指导下独立完成项目研究. 2、“ srt ”计划 项目的形式是在教师指导下，以学生为主体开展课外科学研究活动。参加对象主要为本科生。“ srt ”计划实行导师和同学双向选择，学生可以根据自己的情况选择项目。与课堂教学相比，“ srt ”计划项目中涉及的知识领域更广泛。在这个过程中能充分发挥学生的独立工作能力和能动性，培养学生独立思考和敢于怀疑的批判精神。学生能做到“以我为主”，进行调查研究、查阅文献、分析论证、制定方案、设计或实验、分析总结
等方面的独立能力训练，导师则发挥其主导作用。完成“ srt ”计划的学生可以获得相应的
学分和成绩，其中达到一定水平的还可以取代其相关的课程设计乃至毕业设计（综 合论文训练）。 3、“ srt ”计划 的开展已经取得了很好的成效。以清华大学为例，从 1996 年到 2002 年，已设立“ srt ”项目近 1800 项，参加学生人数 4000 多人。尤其是 2002 年，“ srt ”计划的项目立项数及参加学生人数（共设立“ srt ”项目 632 项，参加学生人数达到 1200 多人）均已增至 1996 年的 6 倍多。 2000 年以来参加“ srt ”计划的同学目前有 1000 多人已经获得课程学分，并涌现了一批优秀指导教师和优秀学生项目，如完成软硬件研制并实际应用获得好评、参加学术会议、发表核心刊物和 sci 等文章、获得专利、参加“挑战杯”科展及系列赛事活动等。 在清华大学的本科生中“ srt ”计划的参加比例已经达到了 30% 左右 。 4、srtp 的核心是支持本科生开展科研训练，学生参与srtp的过程本质上是在进行研究性学习。所以，它注重学生参与研究的学习过程，而并非期望本科生创造出多少原创性成果。srtp 为学有余力的大学生提供直接参与科学研究的机会，引导学生进入科学前沿，了解社 会发展动态。同学们通过发现问题、激发创新思维、独立完成课题等过程，积极主动的探索新
的知识领域，从而体验到一种全新的研究性学习的乐趣。 二、立项来源 1、srtp立项来自：①教师的教学、科研、生产、管理不同领域中的基础性、应用性和开发性的研究课题，根据学生的实际情况，细化后再转化为srtp项目；②学生自定科研项目
和研究课题。 2、srtp是高等学校培养的人才的一种计划，不仅要学习和掌握本专业的基本知识与技
能，而且学习具备创造性地解决所学专业领域内理和实践问题的基本能力。 三、实施“ srt ”计划的目的 第一，优化本科培养模式。开展 “ srt ”计划，是为了进行大学生实践教学改革，培养学生创新意识和创新能力 , 锻炼才干。进一步优化面向社会需要，培养复合型、适应型人才的本科教育培养模式。 第二，增进师生之间的科研交流合作。通过“ srt ”计划，促使教师的教学与科研有机融合 , 形成师生之间的指导、合作关系。 第三，促进产、学、研合作。“ srt ”计划要向社会团体、企业开放，体现产、学、研合作。逐步实现计划的项目来
源、经费资助多样化。 四、哪些研究课题可以作为“ srt ”计划的项目 “ srt ”计划项目的成果要突出应用性、创新性和实用性，特别要适合本科生的研究能力，可以在一年时间内将课题完成。
“ srt ”计划项目的成果形式为新技术、新产品开发和新设计形成的专利、成果，计算机软件，实验研究成果，有一定理论价值与实际应用价值的调查报告，统计分析报告，咨询
报告，学术论文等；有特色的网站建设；相关专业主要课程的 cai 课件或网络课件。 “ srt ”计划的课题按照理工、经管、人文社科等分成三个项目类型，其中经管、人文社科的项目主要是那些需要进行社会调查、统计分析的课题，要求学生能够走向社会，调查
研究。 三、如何申报 “ srt ”计划项目 1、申请时间
2005 年首先在 2003 级本科生中开展 “ srt ”计划项目。 2005 年的申请时间是 9 月份，具体申请时间由教务处实习科下发通知到各学院办公室。 注意：凡按大类招生进行了专业选择的 03 级学生，在申报项目时需按照新的专业班级填写。
SRTP实施方案
