有限元法计算分析程序编写

Size: px

Start display at page:

Download "有限元法计算分析程序编写"

稀虑蔡
5 years ago
Views:

2 m j y o x

3 ) b) c) d)

4 ) b) ) b)

6 β = = P P P P O O = 0 0 P 0 P P 0 O O O O β = P P P O O O O β

7 = = = = P P P P

8 3 5

9 构异运确异

10 ATAB FORTRAN FORTRAN ATAB

11 ATAB ) [ (7 4)] 3 ATAB (+*(7-4))/3^ Eter 3 ATAB s = [ ] + - * / \ ( ) ^ s swer ATAB

12 ) () A=[3;456;789] () Eter (3) ATAB A = 3 3 A = [ ] Eeter [] ATAB A ATAB Cler ATAB ATAB A

13 3) e-3 4.5e33 3 4) ATAB > >

14 5) ATAB ATAB ATAB

15 6) () fucto () ATAB retur (3) globl ATAB edtor/debugger

16 fucto ple % Ths fucto lyze the problem of ple stress use trgle elemets % the vrbles lsts: % ode ode defto % elemet elemet defto % mterl ---- mterl defto % globl stffess mtrx - % f globl ode force vector []{delt}={f} % delt globl ode dsplcemet vector - % bc boudry codto % f ode force % fle_ put fle me % fle_out ---- output fle me % fd_ put fle ID % fd_out output fle ID % globl ode elemet mterl fle_ = put('put fle me: ' 's' ) ; fle_out = put( 'output fle me: ' 's') ;

17 fd_ = fope( fle_ 'r' ) ; % ope put fle ode_umber = fscf( fd_ '%d' ) ; %red ode umber ode = zeros( ode_umber ) ; for =::ode_umber = fscf( fd_ '%d' ) ; ode( : ) = fscf( fd_ '%f' [] ) ; % red ode defto ed elemet_umber = fscf( fd_ '%d' ) ; % red elemet umber elemet = zeros( elemet_umber 4 ) ; for =::elemet_umber e = fscf( fd_ '%d' ) ; elemet( : ) = fscf( fd_ '%d' [4] ) ; % red elemet defto ed

18 mterl_umber = fscf( fd_ '%d' ) ; % red mterl umber mterl = zeros( mterl_umber 3 ) ; for =::mterl_umber m = fscf( fd_ '%d' ) ; mterl( : ) = fscf( fd_ '%f' [3] ) ; % red mterls defto ed bc_umber = fscf( fd_ '%d' ) ; % red boudry codtos umber bc = zeros( bc_umber 3 ) ; for =::bc_umber % red boudry codto defto bc( ) = fscf( fd_ '%d' ) ; bc( ) = fscf( fd_ '%d' ) ; bc( 3 ) = fscf( fd_ '%f' ) ; ed

19 f_umber = fscf( fd_ '%d' ) ; % red ode force umber f = zeros( f_umber 3 ) ; for =::f_umber % red ode force defto f( ) = fscf( fd_ '%d' ) ; f( ) = fscf( fd_ '%d' ) ; f( 3 ) = fscf( fd_ '%f' ) ; ed fclose( fd_ ) ; %close the put fle whch hs bee opeed = sprse( ode_umber * ode_umber * ) ; % defe globl stffess mtrx f = sprse( ode_umber * ) ; % defe globl ode force vector for e=::elemet_umber = Stffesstrx( e ) ; % clculte stffess mtrx of every elemet AssembleStffesstrx( e ) ; % ssemble elemet stffess mtrx ed

20 for f=::f_umber % ccess the ode force = f( f ) ; d = f( f ) ; f( (-)* + d ) = f( (-)* + d ) +f( f 3 ) ; ed for bc=::bc_umber % ccess the boudry codto = bc(bc ) ; % method: multple bg umber d = bc(bc ) ; m = (-)* + d ; f(m) = bc(bc 3)* (mm) * e5 ; (m m) = (m m) * e5 ; ed delt = / f ; % solve the ler equto set

21 fd_out = fope( fle_out 'w' ) ; % ope the output fle for =::ode_umber ed fprtf( fd_out 'ode = %5d u = %.8e v = %.8e\'... delt((-)*+) delt((-)*+) ) ; % output the ode dsplcemet fclose(fd_out) ; % close the output fle retur ;

22 fucto = Stffesstrx( e ) % ths fucto clculte the stffess mtrx of elemet e % the vrbles lst: % elemet stffess mtrx % h the thcess of the elemet % E Youg's modulus % mu Posso's rto % x y xj yj xm ym the coordtes of the ode j m % b c some vrbles defed s followg (jm) % = xj*ym - xm*yj b = yj - ym c = -(xj-xm) ( j m) % delt the re of the trgle elemet % B B mtrx % D D mtrx

23 globl ode elemet mterl = zeros( 6 6 ) ; h = mterl( elemet(e 4) ) ; E = mterl( elemet(e 4) ) ; mu = mterl( elemet(e 4) 3 ) ; x = ode( elemet( e ) ) ; y = ode( elemet( e ) ) ; xj = ode( elemet( e ) ) ; yj = ode( elemet( e ) ) ; xm = ode( elemet( e 3 ) ) ; ym = ode( elemet( e 3 ) ) ;

24 = xj*ym - xm*yj ; j = xm*y - x*ym ; m = x*yj - xj*y ; b = yj - ym ; bj = ym - y ; bm = y - yj ; c = -(xj-xm) ; cj = -(xm-x) ; cm = -(x-xj) ; delt = bs((+j+m)/) ; B = [b 0 bj 0 bm 0 0 c 0 cj 0 cm c b cj bj cm bm] ; B = B//delt ; D = [ mu 0 mu (-mu)/] ; D = D*E/(-mu^) ; = trspose(b)*d*b*h*delt ; retur

25 fucto AssembleStffesstrx( e ) % ths fucto ssemble stffess mtrx of elemet globl elemet for =::3 for j=::3 for p=:: for q=:: m = (-)*+p ; = (j-)*+q ; = (elemet(e)-)*+p ; N = (elemet(ej)-)*+q ; (N) = (N) + (m) ; ed ed ed ed retur

28 0.0.0E

FZ1.s92

FZ1.s92 第一章服装制作工艺基础知识第一节服装术语服装术语是服装技术专用语, 如服装的每一个品种服装上的每一块裁片服装制作过程中每一种操作, 以及所使用的工具服装制作过程中出现的各种弊病和要达到的某些质量要求等, 都有其