欧拉折线法绘制电场线 | 寒山烟雨
现在的位置: 首页 > 小试身手 > 正文

欧拉折线法绘制电场线

2014年06月29日 小试身手 ⁄ 共 489字 ⁄ 字号 评论 2 条 ⁄ 阅读 3,323 views 次

利用最常规的模拟方法——欧拉折线法绘制四电荷电场线,fortran写的,下一个版本会有python的Vpython动画版。

 Fortran | 
 
 copy code |
?

01
module conts
02
integer :: N=18,Num=4
03
real(4)	:: q(4)
04
data q/1.0,-1.0,1.0,-1.0/
05
complex	:: p(4)
06
data p/(1.0,1.0),(-1.0,1.0),(-1.0,-1.0),(1.0,-1.0)/
07
real(4)	:: pai=3.1415926
08
real(4)	:: ds=0.01
09
real(4)	:: dr=0.3
10
end
11
 
12
program Main
13
use conts
14
complex	:: r,s
15
integer	:: i,j
16
open(10,file="data.dat")
17
do i=1,Num
18
do j=1,N
19
s=exp(j*pai/N*(0,1)*2.)
20
r=p(i)+dr*s
21
call plot(r)
22
enddo
23
enddo
24
end
25
 
26
subroutine plot(r)
27
use conts
28
complex	:: E,r
29
do
30
write(10,*)	real(r),aimag(r)
31
r=r+ds*E(r)/abs(E(r))
32
if(abs(r-p(1))<dr.or.abs(r-p(2))<dr.or.abs(r-p(3))<dr.or.abs(r-p(4))<dr.or.abs(r)>20.0) exit
33
enddo
34
end
35
 
36
complex function E(r)
37
use conts
38
complex	:: r
39
integer	:: i
40
E=(0,0)
41
do i=1,Num
42
E=E+q(i)*(r-p(i))/abs(r-p(i))**3
43
enddo
44
return
45
end

效果如图所示电场线

0

目前有 2 条留言    访客:2 条, 博主:0 条

  1. 玻璃钢电缆支架 2015年04月28日 下午2:39  @回复  Δ-49楼 回复

    利用最常规的模拟方法——欧拉折线法绘制四电荷电场线,fortran写的,下一个版本会有python的Vpython动画版。

  2. 猎户项圈 2014年06月29日 下午4:37  @回复  Δ-48楼 回复

    留言是种美德,写点什么...

给我留言

留言无头像?


×