价带结构

这是pscad运行的电弧炉模型，是学习pscad很好的例子(This is pscad running electric arc furnace model, is learning pscad good example)

用于abaqus二次开发vumat子程序，适用于三维梁单元的线弹性模型(Secondary development of abaqus vumat subroutine, be applicable to the linear elastic model of three-dimensional beam element )

用于abaqus模拟三维马氏体相变的UEL子程序。(3D martensitic transformation UEL subroutine.)

说明：  sw差分格式求解一维黎曼问题，深入理解黎曼问题，学习计算流体力学(Sw difference scheme to solve one-dimensional Riemann problem, in-depth understanding of Riemann problem, learning computational fluid dynamics)

二阶迎风差分格式。代码简单用于计算二阶迎风差分格式的fortran程序。(Fortran program about 2nd-upwind.)

平面波方法计算声子晶体能带结构的源程序（复式格子圆柱体按正方形排列,水、四氯化碳为散射体,水银为基体）(Plane wave method phononic crystal band structure of the source (double-cylinder square array grid, water, carbon tetrachloride, the scattering body, mercury as a matrix))

PSCAD仿真程序，Buck降压电路，有PI控制，很有参考价值(PSACD program,Buck Chopper)

ABAQUS 二维fortran粘弹性边界单元；可用于动力计算(ABAQUS 2D fortran viscoelastic boundary element for dynamic calculation)

用于计算不同阶数的第一类贝塞尔函数Fortran代码!第一类n阶BESSEL函数bessjfunction bessj(n,x)integer n,IACCreal(8) bessj,x,BIGNO,BIGNIPARAMETER (IACC=40,BIGNO=1.e10,BIGNI=1.e-10)INTEGER j,jsum,mreal(8) ax,bj,bjm,bjp,sum,tox,bessj0,bessj1if(nfloat(n))then tox=2./ax bjm=bessj0(ax) bj=bessj1(ax) do j=1,n-1    bjp=j*tox*bj-bjmbjm=bjbj=bjp end do bessj=bjelse tox=2./ax m=2*((n+int(sqrt(float(IACC*n))))/2) bessj=0. jsum=0 sum=0. bjp=0. bj=1. do j=m,1,-1    bjm=j*tox*bj-bjpbjp=bjbj=bjmif(abs(bj)>BIGNO)then  bj=bj*BIGNI  bjp=bjp*BIGNI  bessj=bessj*BIGNI  sum=sum*BIGNIendifif(jsum/=0) s

利用快速解耦法进行电力系统潮流计算，对配电网进行潮流计算(Use of PQ decomposition method for the electric power flow calculation system)