FrontISTR  5.9.0
Large-scale structural analysis program with finit element method
dynamic_output.f90
Go to the documentation of this file.
1 !-------------------------------------------------------------------------------
2 ! Copyright (c) 2019 FrontISTR Commons
3 ! This software is released under the MIT License, see LICENSE.txt
4 !-------------------------------------------------------------------------------
6 module m_dynamic_output
7  implicit none
8 contains
9 
11  !----------------------------------------------------------------------*
12  subroutine fstr_dynamic_output( cstep, istep, t_curr, hecMESH, fstrSOLID, fstrDYNAMIC, fstrPARAM, outflag )
13  !----------------------------------------------------------------------*
14  use m_fstr
15  use m_fstr_nodalstress
16  use m_make_result
17  use m_hecmw2fstr_mesh_conv
18  integer, intent(in) :: cstep
19  integer, intent(in) :: istep
20  real(kind=kreal), intent(in) :: t_curr
21  type(hecmwst_local_mesh), intent(in) :: hecmesh
22  type(fstr_solid), intent(inout) :: fstrSOLID
23  type(fstr_dynamic), intent(in) :: fstrDYNAMIC
24  type(fstr_param), intent(in) :: fstrPARAM
25  logical, intent(in) :: outflag
26 
27  type(hecmwst_result_data) :: fstrRESULT
28  integer(kind=kint) :: i, j, ndof, fnum, is, iE, gid
29 
30  ndof = hecmesh%n_dof
31 
32  if( fstrsolid%TEMP_ngrp_tot>0 .or. fstrsolid%TEMP_irres>0 ) then
33  if( ndof==3 ) then
34  allocate( fstrsolid%tnstrain(hecmesh%n_node*6) )
35  allocate( fstrsolid%testrain(hecmesh%n_elem*6) )
36  else if( ndof==2 ) then
37  allocate( fstrsolid%tnstrain(hecmesh%n_node*3) )
38  allocate( fstrsolid%testrain(hecmesh%n_elem*3) )
39  else if( ndof== 4) then
40  write(*,*)'Error: This routine is not implemented'
41  stop
42  endif
43  endif
44 
45  if( ndof==3 .or. ndof == 4 ) then
46  call fstr_nodalstress3d( hecmesh, fstrsolid )
47  else if( ndof==2 ) then
48  call fstr_nodalstress2d( hecmesh, fstrsolid )
49  else if( ndof==6 ) then
50  call fstr_nodalstress6d( hecmesh, fstrsolid )
51  endif
52 
53  if( associated( fstrsolid%contacts ) ) &
54  & call setup_contact_output_variables( hecmesh, fstrsolid, -1 )
55 
56  if( (mod(istep,fstrsolid%output_ctrl(1)%frequency)==0 .or. outflag) ) then
57  fnum = fstrsolid%output_ctrl(1)%filenum
58  call fstr_dynamic_post( fnum, cstep, istep, hecmesh, fstrsolid, fstrdynamic )
59  endif
60 
61  if( fstrsolid%output_ctrl(2)%outinfo%grp_id>0 .and. &
62  (mod(istep,fstrsolid%output_ctrl(2)%frequency)==0 .or. outflag) ) then
63  is = fstrsolid%output_ctrl(2)%outinfo%grp_id
64  fnum = fstrsolid%output_ctrl(2)%filenum
65  do i = hecmesh%node_group%grp_index(is-1)+1, hecmesh%node_group%grp_index(is)
66  ie = hecmesh%node_group%grp_item(i)
67  gid = hecmesh%global_node_ID(ie)
68  write(fnum,'(2i10,1p6e15.7)') istep,gid,(fstrsolid%unode(ndof*(ie-1)+j),j=1,ndof)
69  enddo
70  endif
71 
72  if( iresult==1 .and. &
73  (mod(istep,fstrsolid%output_ctrl(3)%frequency)==0 .or. outflag) ) then
74  if( associated( fstrsolid%contacts ) ) &
75  & call setup_contact_output_variables( hecmesh, fstrsolid, 3, fstrdynamic%t_delta )
76  call fstr_write_result( hecmesh, fstrsolid, fstrparam, istep, t_curr, 0, fstrdynamic )
77  endif
78 
79  if( ivisual==1 .and. &
80  (mod(istep,fstrsolid%output_ctrl(4)%frequency)==0 .or. outflag) ) then
81 
82  if( associated( fstrsolid%contacts ) ) &
83  & call setup_contact_output_variables( hecmesh, fstrsolid, 4, fstrdynamic%t_delta )
84  call fstr_make_result( hecmesh, fstrsolid, fstrresult, istep, t_curr, fstrdynamic )
85  call fstr2hecmw_mesh_conv( hecmesh )
86  call hecmw_visualize_init
87  call hecmw_visualize( hecmesh, fstrresult, istep )
88  call hecmw_visualize_finalize
89  call hecmw2fstr_mesh_conv( hecmesh )
90  call hecmw_result_free( fstrresult )
91  endif
92 
93  end subroutine fstr_dynamic_output
94 
96  !----------------------------------------------------------------------*
97  subroutine fstr_dynamic_post( fnum, cstep, istep, hecMESH, fstrSOLID, fstrDYNAMIC )
98  !----------------------------------------------------------------------*
99  use m_fstr
100  integer, intent(in) :: fnum
101  integer, intent(in) :: cstep
102  integer, intent(in) :: istep
103  type(hecmwst_local_mesh), intent(in) :: hecMESH
104  type(fstr_solid), intent(in) :: fstrSOLID
105  type(fstr_dynamic), intent(in) :: fstrDYNAMIC
106 
107  real(kind=kreal) :: umax(6), umin(6), vmax(6), vmin(6), amax(6), amin(6)
108  real(kind=kreal) :: emax(14), emin(14), smax(14), smin(14)
109  real(kind=kreal) :: mmax(1), mmin(1), emmax(1), emmin(1)
110  real(kind=kreal) :: eemax(14), eemin(14), esmax(14), esmin(14)
111 
112  real(kind=kreal) :: gumax(6), gumin(6), gvmax(6), gvmin(6), gamax(6), gamin(6)
113  real(kind=kreal) :: gemax(14), gemin(14), gsmax(14), gsmin(14)
114  real(kind=kreal) :: gmmax(1), gmmin(1), gemmax(1), gemmin(1)
115  real(kind=kreal) :: geemax(14), geemin(14), gesmax(14), gesmin(14)
116 
117  integer(kind=kint) :: IUmax(6), IUmin(6), IVmax(6), IVmin(6), IAmax(6), IAmin(6)
118  integer(kind=kint) :: IEmax(14), IEmin(14), ISmax(14), ISmin(14)
119  integer(kind=kint) :: IMmax(1), IMmin(1), IEMmax(1), IEMmin(1)
120  integer(kind=kint) :: IEEmax(14), IEEmin(14), IESmax(14), IESmin(14)
121  integer(kind=kint) :: i, j, k, ndof, mdof, ID_area, idx
122  integer(kind=kint) :: label(6)
123 
124  if( istep==0 ) return
125  if( fstrdynamic%idx_eqa==1 .and. istep>0 ) then
126  idx = 2
127  else
128  idx = 1
129  endif
130  ndof = hecmesh%n_dof
131 
132  write( fnum, '(''#### Result step, substep='',I6,I6)') cstep, istep
133  select case (ndof)
134  case (2)
135  mdof = 3
136  label(1)=11;label(2)=22
137  label(3)=12
138  case (3,4,6)
139  mdof = 6
140  label(1)=11;label(2)=22;label(3)=33;
141  label(4)=12;label(5)=23;label(6)=31;
142  end select
143 
144  j = hecmesh%global_node_ID(1)
145  do k = 1, ndof
146  umax(k) = fstrdynamic%DISP(k,idx); umin(k) = fstrdynamic%DISP(k,idx)
147  vmax(k) = fstrdynamic%VEL(k,idx); vmin(k) = fstrdynamic%VEL(k,idx)
148  amax(k) = fstrdynamic%ACC(k,idx); amin(k) = fstrdynamic%ACC(k,idx)
149  iumax(k)= j; iumin(k)= j; ivmax(k)= j; ivmin(k)= j; iamax(k)= j; iamin(k)= j
150  enddo
151  do k = 1, mdof
152  emax(k) = fstrsolid%STRAIN(k); emin(k) = fstrsolid%STRAIN(k)
153  smax(k) = fstrsolid%STRESS(k); smin(k) = fstrsolid%STRESS(k)
154  eemax(k) = fstrsolid%ESTRAIN(k); eemin(k) = fstrsolid%ESTRAIN(k)
155  esmax(k) = fstrsolid%ESTRESS(k); esmin(k) = fstrsolid%ESTRESS(k)
156  iemax(k)= j; iemin(k)= j; ismax(k)= j; ismin(k)= j
157  ieemax(k)= j; ieemin(k)= j; iesmax(k)= j; iesmin(k)= j
158  enddo
159  mmax(1) = fstrsolid%MISES(1); mmin(1) = fstrsolid%MISES(1)
160  emmax(1) = fstrsolid%EMISES(1); emmin(1) = fstrsolid%EMISES(1)
161  immax(1)= j; immin(1)= j; iemmax(1)= j; iemmin(1)= j
162 
163  !C*** Show Displacement / Velosity / Acc
164  do i = 1, hecmesh%nn_internal
165  if(fstrsolid%is_rot(i)==1)cycle
166  j = hecmesh%global_node_ID(i)
167  do k = 1, ndof
168  if ( fstrdynamic%DISP(ndof*(i-1)+k,idx) > umax(k) ) then
169  umax(k) = fstrdynamic%DISP(ndof*(i-1)+k,idx)
170  iumax(k)= j
171  else if( fstrdynamic%DISP(ndof*(i-1)+k,idx) < umin(k) ) then
172  umin(k) = fstrdynamic%DISP(ndof*(i-1)+k,idx)
173  iumin(k)= j
174  endif
175  if ( fstrdynamic%VEL(ndof*(i-1)+k,idx) > vmax(k) ) then
176  vmax(k) = fstrdynamic%VEL(ndof*(i-1)+k,idx)
177  ivmax(k)= j
178  else if( fstrdynamic%VEL(ndof*(i-1)+k,idx) < vmin(k) ) then
179  vmin(k) = fstrdynamic%VEL(ndof*(i-1)+k,idx)
180  ivmin(k)= j
181  endif
182  if ( fstrdynamic%ACC(ndof*(i-1)+k,idx) > amax(k) ) then
183  amax(k) = fstrdynamic%ACC(ndof*(i-1)+k,idx)
184  iamax(k)= j
185  else if( fstrdynamic%ACC(ndof*(i-1)+k,idx) < amin(k) ) then
186  amin(k) = fstrdynamic%ACC(ndof*(i-1)+k,idx)
187  iamin(k)= j
188  endif
189  enddo
190  enddo
191  !C*** Nodal Strain / Stress / MISES
192  !C @node
193  do i = 1, hecmesh%nn_internal
194  if(fstrsolid%is_rot(i)==1)cycle
195  j = hecmesh%global_node_ID(i)
196  do k = 1, mdof
197  if ( fstrsolid%STRAIN(mdof*(i-1)+k) > emax(k) ) then
198  emax(k) = fstrsolid%STRAIN(mdof*(i-1)+k)
199  iemax(k)= j
200  else if( fstrsolid%STRAIN(mdof*(i-1)+k) < emin(k) ) then
201  emin(k) = fstrsolid%STRAIN(mdof*(i-1)+k)
202  iemin(k)= j
203  endif
204  if ( fstrsolid%STRESS(mdof*(i-1)+k) > smax(k) ) then
205  smax(k) = fstrsolid%STRESS(mdof*(i-1)+k)
206  ismax(k)= j
207  else if( fstrsolid%STRESS(mdof*(i-1)+k) < smin(k) ) then
208  smin(k) = fstrsolid%STRESS(mdof*(i-1)+k)
209  ismin(k)= j
210  endif
211  enddo
212  if ( fstrsolid%MISES(i) > mmax(1) ) then
213  mmax(1) = fstrsolid%MISES(i)
214  immax(1)= j
215  else if( fstrsolid%MISES(i) < mmin(1) ) then
216  mmin(1) = fstrsolid%MISES(i)
217  immin(1)= j
218  endif
219  enddo
220  !C*** Elemental Strain / STRESS
221  do i = 1, hecmesh%n_elem
222  id_area = hecmesh%elem_ID(i*2)
223  if( id_area==hecmesh%my_rank ) then
224  j = hecmesh%global_elem_ID(i)
225  do k = 1, mdof
226  if( fstrsolid%ESTRAIN(mdof*(i-1)+k) > eemax(k) ) then
227  eemax(k) = fstrsolid%ESTRAIN(mdof*(i-1)+k)
228  ieemax(k)= j
229  else if( fstrsolid%ESTRAIN(mdof*(i-1)+k) < eemin(k) ) then
230  eemin(k) = fstrsolid%ESTRAIN(mdof*(i-1)+k)
231  ieemin(k)= j
232  endif
233  if( fstrsolid%ESTRESS(mdof*(i-1)+k) > esmax(k) ) then
234  esmax(k) = fstrsolid%ESTRESS(mdof*(i-1)+k)
235  iesmax(k)= j
236  else if( fstrsolid%ESTRESS(mdof*(i-1)+k) < esmin(k) ) then
237  esmin(k) = fstrsolid%ESTRESS(mdof*(i-1)+k)
238  iesmin(k)= j
239  endif
240  enddo
241  if( fstrsolid%EMISES(i) > emmax(1) ) then
242  emmax(1) = fstrsolid%EMISES(i)
243  iemmax(1)= j
244  else if( fstrsolid%EMISES(i) < emmin(1) ) then
245  emmin(1) = fstrsolid%EMISES(i)
246  iemmin(1)= j
247  endif
248  endif
249  enddo
250 
251 
252  write(ilog,*) '##### Local Summary @Node :Max/IdMax/Min/IdMin####'
253  do i = 1, ndof; write(ilog,1029) ' //U',i, ' ',umax(i),iumax(i),umin(i),iumin(i); end do
254  if (ndof == 4) write(ilog,1009) ' //P ' , umax(4),iumax(4),umin(4),iumin(4)
255  do i = 1, ndof; write(ilog,1029) ' //V',i, ' ',vmax(i),ivmax(i),vmin(i),ivmin(i); end do
256  do i = 1, ndof; write(ilog,1029) ' //A',i, ' ',amax(i),iamax(i),amin(i),iamin(i); end do
257  do i = 1, mdof; write(ilog,1029) ' //E',label(i),' ',emax(i),iemax(i),emin(i),iemin(i); end do
258  do i = 1, mdof; write(ilog,1029) ' //S',label(i),' ',smax(i),ismax(i),smin(i),ismin(i); end do
259  write(ilog,1009) '//SMS ' ,mmax(1),immax(1),mmin(1),immin(1)
260  write(ilog,*) '##### Local Summary @Element :Max/IdMax/Min/IdMin####'
261  do i = 1, mdof; write(ilog,1029) ' //E',label(i),' ',eemax(i),ieemax(i),eemin(i),ieemin(i); end do
262  do i = 1, mdof; write(ilog,1029) ' //S',label(i),' ',esmax(i),iesmax(i),esmin(i),iesmin(i); end do
263  write(ilog,1009) '//SMS ' ,emmax(1),iemmax(1),emmin(1),iemmin(1)
264 
265  !C*** Show Summary
266  gumax = umax; gumin = umin; gvmax = vmax; gvmin = vmin; gamax = amax; gamin = amin;
267  gemax = emax; gemin = emin; geemax = eemax; geemin = eemin;
268  gsmax = smax; gsmin = smin; gesmax = esmax; gesmin = esmin;
269  gmmax = mmax; gmmin = mmin; gemmax = emmax; gemmin = emmin;
270 
271  call hecmw_allreduce_r(hecmesh,gumax,ndof,hecmw_max)
272  call hecmw_allreduce_r(hecmesh,gumin,ndof,hecmw_min)
273  call hecmw_allreduce_r(hecmesh,gvmax,ndof,hecmw_max)
274  call hecmw_allreduce_r(hecmesh,gvmin,ndof,hecmw_min)
275  call hecmw_allreduce_r(hecmesh,gamax,ndof,hecmw_max)
276  call hecmw_allreduce_r(hecmesh,gamin,ndof,hecmw_min)
277  call hecmw_allreduce_r(hecmesh,gemax,mdof,hecmw_max)
278  call hecmw_allreduce_r(hecmesh,gemin,mdof,hecmw_min)
279  call hecmw_allreduce_r(hecmesh,gsmax,mdof,hecmw_max)
280  call hecmw_allreduce_r(hecmesh,gsmin,mdof,hecmw_min)
281  call hecmw_allreduce_r(hecmesh,gmmax,1,hecmw_max)
282  call hecmw_allreduce_r(hecmesh,gmmin,1,hecmw_min)
283  call hecmw_allreduce_r(hecmesh,geemax,mdof,hecmw_max)
284  call hecmw_allreduce_r(hecmesh,geemin,mdof,hecmw_min)
285  call hecmw_allreduce_r(hecmesh,gesmax,mdof,hecmw_max)
286  call hecmw_allreduce_r(hecmesh,gesmin,mdof,hecmw_min)
287  call hecmw_allreduce_r(hecmesh,gemmax,1,hecmw_max)
288  call hecmw_allreduce_r(hecmesh,gemmin,1,hecmw_min)
289 
290  do i=1,ndof
291  if(gumax(i) > umax(i)) iumax(i) = 0
292  if(gvmax(i) > vmax(i)) ivmax(i) = 0
293  if(gamax(i) > amax(i)) iamax(i) = 0
294  if(gumin(i) < umin(i)) iumin(i) = 0
295  if(gvmin(i) < vmin(i)) ivmin(i) = 0
296  if(gamin(i) < amin(i)) iamin(i) = 0
297  enddo
298  do i=1,mdof
299  if(gemax(i) > emax(i)) iemax(i) = 0
300  if(gsmax(i) > smax(i)) ismax(i) = 0
301  if(geemax(i) > eemax(i)) ieemax(i) = 0
302  if(gesmax(i) > esmax(i)) iesmax(i) = 0
303  if(gemin(i) < emin(i)) iemin(i) = 0
304  if(gsmin(i) < smin(i)) ismin(i) = 0
305  if(geemin(i) < eemin(i)) ieemin(i) = 0
306  if(gesmin(i) < esmin(i)) iesmin(i) = 0
307  enddo
308  do i=1,1
309  if(gmmax(i) > mmax(i)) immax(i) = 0
310  if(gemmax(i) > emmax(i)) iemmax(i) = 0
311  if(gmmin(i) < mmin(i)) immin(i) = 0
312  if(gemmin(i) < emmin(i)) iemmin(i) = 0
313  enddo
314 
315  call hecmw_allreduce_i(hecmesh,iumax,ndof,hecmw_max)
316  call hecmw_allreduce_i(hecmesh,iumin,ndof,hecmw_max)
317  call hecmw_allreduce_i(hecmesh,ivmax,ndof,hecmw_max)
318  call hecmw_allreduce_i(hecmesh,ivmin,ndof,hecmw_max)
319  call hecmw_allreduce_i(hecmesh,iamax,ndof,hecmw_max)
320  call hecmw_allreduce_i(hecmesh,iamin,ndof,hecmw_max)
321  call hecmw_allreduce_i(hecmesh,iemax,mdof,hecmw_max)
322  call hecmw_allreduce_i(hecmesh,iemin,mdof,hecmw_max)
323  call hecmw_allreduce_i(hecmesh,ismax,mdof,hecmw_max)
324  call hecmw_allreduce_i(hecmesh,ismin,mdof,hecmw_max)
325  call hecmw_allreduce_i(hecmesh,immax,1,hecmw_max)
326  call hecmw_allreduce_i(hecmesh,immin,1,hecmw_max)
327  call hecmw_allreduce_i(hecmesh,ieemax,mdof,hecmw_max)
328  call hecmw_allreduce_i(hecmesh,ieemin,mdof,hecmw_max)
329  call hecmw_allreduce_i(hecmesh,iesmax,mdof,hecmw_max)
330  call hecmw_allreduce_i(hecmesh,iesmin,mdof,hecmw_max)
331  call hecmw_allreduce_i(hecmesh,iemmax,1,hecmw_max)
332  call hecmw_allreduce_i(hecmesh,iemmin,1,hecmw_max)
333 
334  if( hecmesh%my_rank==0 ) then
335  write(ilog,*) '##### Global Summary @Node :Max/IdMax/Min/IdMin####'
336  do i = 1, ndof; write(ilog,1029) ' //U',i, ' ',gumax(i),iumax(i),gumin(i),iumin(i); end do
337  if (ndof == 4) write(ilog,1009) ' //P ' ,gumax(4),iumax(4),gumin(4),iumin(4)
338  do i = 1, ndof; write(ilog,1029) ' //V',i, ' ',gvmax(i),ivmax(i),gvmin(i),ivmin(i); end do
339  do i = 1, ndof; write(ilog,1029) ' //A',i, ' ',gamax(i),iamax(i),gamin(i),iamin(i); end do
340  do i = 1, mdof; write(ilog,1029) ' //E',label(i),' ',gemax(i),iemax(i),gemin(i),iemin(i); end do
341  do i = 1, mdof; write(ilog,1029) ' //S',label(i),' ',gsmax(i),ismax(i),gsmin(i),ismin(i); end do
342  write(ilog,1009) '//SMS ' ,gmmax(1),immax(1),gmmin(1),immin(1)
343  write(ilog,*) '##### Global Summary @Element :Max/IdMax/Min/IdMin####'
344  do i = 1, mdof; write(ilog,1029) ' //E',label(i),' ',geemax(i),ieemax(i),geemin(i),ieemin(i); end do
345  do i = 1, mdof; write(ilog,1029) ' //S',label(i),' ',gesmax(i),iesmax(i),gesmin(i),iesmin(i); end do
346  write(ilog,1009) '//SMS ' ,gemmax(1),iemmax(1),gemmin(1),iemmin(1)
347  endif
348 
349  1009 format(a7,1pe12.4,i10,1pe12.4,i10)
350  1029 format(a,i0,a,1pe12.4,i10,1pe12.4,i10)
351  end subroutine fstr_dynamic_post
352 
353  !C================================================================C
354  !C-- subroutine dynamic_output_monit
355  !C================================================================C
356  subroutine dynamic_output_monit(cstep, istep, t_curr, hecMESH, fstrPARAM, fstrDYNAMIC, fstrEIG, fstrSOLID)
357  use m_fstr
358  integer, intent(in) :: cstep
359  integer, intent(in) :: istep
360  real(kind=kreal), intent(in) :: t_curr
361  type(hecmwst_local_mesh) :: hecmesh
362  type(fstr_param) :: fstrPARAM
363  type(fstr_dynamic) :: fstrDYNAMIC
364  type(fstr_eigen) :: fstrEIG
365  type(fstr_solid) :: fstrSOLID
366 
367  integer(kind=kint) :: idx, ii, jj, ierr, ncmp
368  integer(kind=kint) :: num_monit, ig, is, iE, ik, iunitS, iunit
369  logical :: yes
370 
371  if( mod(istep,fstrdynamic%nout_monit)/=0 ) return
372 
373  if( fstrdynamic%idx_eqa==1 .and. istep>0 ) then
374  idx = 2
375  else
376  idx = 1
377  endif
378 
379  num_monit = 0
380  ig = fstrdynamic%ngrp_monit
381  is = hecmesh%node_group%grp_index(ig-1)+1
382  ie = hecmesh%node_group%grp_index(ig)
383  do ik=is,ie
384  ii = hecmesh%node_group%grp_item(ik)
385  if (ii > hecmesh%nn_internal) cycle
386  num_monit = num_monit+1
387  jj = hecmesh%global_node_id(ii)
388  iunits = 10*(num_monit-1)
389 
390  !C-- displacement
391  if( fstrdynamic%iout_list(1)==1 ) then
392  iunit = iunits + fstrdynamic%dynamic_IW4
393  write( iunit, '(i10, i10,1pe13.4e3,i10,1p6e13.4e3)' ) &
394  cstep, istep, t_curr, jj, &
395  fstrdynamic%DISP( hecmesh%n_dof*(ii-1)+1 : hecmesh%n_dof*ii , idx )
396  end if
397  !C-- velocity
398  if( fstrdynamic%iout_list(2)==1 ) then
399  iunit = iunits + fstrdynamic%dynamic_IW5
400  write( iunit, '(i10, i10,1pe13.4e3,i10,1p6e13.4e3)' ) &
401  cstep, istep, t_curr, jj, &
402  fstrdynamic%VEL( hecmesh%n_dof*(ii-1)+1 : hecmesh%n_dof*ii , idx )
403  end if
404  !C-- acceleration
405  if( fstrdynamic%iout_list(3)==1 ) then
406  iunit = iunits + fstrdynamic%dynamic_IW6
407  write( iunit, '(i10, i10,1pe13.4e3,i10,1p6e13.4e3)' ) &
408  cstep, istep, t_curr, jj, &
409  fstrdynamic%ACC( hecmesh%n_dof*(ii-1)+1 : hecmesh%n_dof*ii , idx )
410  end if
411  !C-- nodal force
412  if( fstrdynamic%iout_list(4)==1 ) then
413  iunit = iunits + fstrdynamic%dynamic_IW7
414  write( iunit, '(i10, i10,1pe13.4e3,i10,1p6e13.4e3)' ) &
415  cstep, istep, t_curr, jj, &
416  fstrsolid%QFORCE( hecmesh%n_dof*(ii-1)+1 : hecmesh%n_dof*ii )
417  end if
418  !C-- strain
419  if( fstrdynamic%iout_list(5) > 0 ) then
420  if (hecmesh%n_dof == 2 .or. hecmesh%n_dof == 3 .or. hecmesh%n_dof == 4) then
421  ncmp = 6
422  else
423  ncmp = 12
424  endif
425  iunit = iunits + fstrdynamic%dynamic_IW8
426  write( iunit, '(i10, i10,1pe13.4e3,i10,1p6e13.4e3)') &
427  cstep, istep, t_curr, jj, &
428  fstrsolid%STRAIN( ncmp*(ii-1)+1 : ncmp*ii )
429  end if
430  !C-- stress
431  if( fstrdynamic%iout_list(6) > 0 ) then
432  if (hecmesh%n_dof == 2 .or. hecmesh%n_dof == 3 .or. hecmesh%n_dof == 4) then
433  ncmp = 6
434  else
435  ncmp = 12
436  endif
437  iunit = iunits + fstrdynamic%dynamic_IW9
438  write( iunit, '(i10, i10,1pe13.4e3,i10,1p7e13.4e3)') &
439  cstep, istep, t_curr, jj, &
440  fstrsolid%STRESS( ncmp*(ii-1)+1 : ncmp*ii )
441  end if
442  enddo
443 
444  if( hecmesh%my_rank==0 ) then
445  if( any(fstrdynamic%iout_list(1:3)==1) ) then
446  inquire( file='dyna_energy.txt', opened=yes )
447  if( .not. yes ) then
448  open( fstrdynamic%dynamic_IW10, file='dyna_energy.txt', status='replace', iostat=ierr )
449  if( ierr/=0 ) then
450  write(*,*) 'stop due to file opening error <dyna_enrgy.txt>'
451  call hecmw_abort( hecmw_comm_get_comm() )
452  endif
453  write( fstrdynamic%dynamic_IW10, * ) &
454  ' time step', ' time ', ' kinetic energy', ' strain energy', ' total energy'
455  endif
456  if(istep==0) then
457  fstrdynamic%kineticEnergy = 0.0d0
458  do ii = 1, hecmesh%n_node*hecmesh%n_dof
459  fstrdynamic%kineticEnergy = fstrdynamic%kineticEnergy &
460  + 0.5d0 * fstreig%mass(ii) * fstrdynamic%VEL(ii,idx) * fstrdynamic%VEL(ii,idx)
461  enddo
462  endif
463  fstrdynamic%totalEnergy = fstrdynamic%kineticEnergy + fstrdynamic%strainEnergy
464  write( fstrdynamic%dynamic_IW10, '(i10, i10,1pe13.4e3,1p3e16.4e3)' ) &
465  cstep, istep, t_curr, &
466  fstrdynamic%kineticEnergy, fstrdynamic%strainEnergy, fstrdynamic%totalEnergy
467  endif
468  if( istep==fstrdynamic%n_step ) close(fstrdynamic%dynamic_IW10)
469  endif
470  end subroutine dynamic_output_monit
471 
472  !C================================================================C
473  !C-- subroutine matvec
474  !C================================================================C
475  subroutine matvec(y,x,hecMAT,ndof,D,AU,AL)
476  use m_fstr
477  type(hecmwst_matrix) :: hecMAT
478 
479  integer(kind=kint) :: ndof, i, is, ie, j, icol, ki, kj, ix, iy, ip, nn
480  real(kind=kreal) :: d(ndof*ndof*hecmat%NP)
481  real(kind=kreal) :: au(ndof*ndof*hecmat%NPU)
482  real(kind=kreal) :: al(ndof*ndof*hecmat%NPL)
483  real(kind=kreal) :: x(ndof*hecmat%NP)
484  real(kind=kreal) :: y(ndof*hecmat%NP)
485 
486  nn=ndof*ndof
487 
488  y=0.0d0
489 
490  do i=1,hecmat%NP
491  is=hecmat%indexU(i-1)+1
492  ie=hecmat%indexU(i)
493  do j=is,ie
494  icol=hecmat%itemU(j)
495  do ki=1,ndof
496  iy=ndof*(i-1)+ki
497  do kj=1,ndof
498  ix=ndof*(icol-1)+kj
499  ip=nn*(j-1)+ndof*(ki-1)+kj
500  y(iy)=y(iy)+au(ip)*x(ix)
501  enddo
502  enddo
503  enddo
504  enddo
505 
506  do i=1,hecmat%NP
507  is=hecmat%indexL(i-1)+1
508  ie=hecmat%indexL(i)
509  do j=is,ie
510  icol=hecmat%itemL(j)
511  do ki=1,ndof
512  iy=ndof*(i-1)+ki
513  do kj=1,ndof
514  ix=ndof*(icol-1)+kj
515  ip=nn*(j-1)+ndof*(ki-1)+kj
516  y(iy)=y(iy)+al(ip)*x(ix)
517  enddo
518  enddo
519  enddo
520  enddo
521 
522  do i=1,hecmat%NP
523  do ki=1,ndof
524  iy=ndof*(i-1)+ki
525  do kj=1,ndof
526  ix=ndof*(i-1)+kj
527  ip=nn*(i-1)+ndof*(ki-1)+kj
528  y(iy)=y(iy)+d(ip)*x(ix)
529  enddo
530  enddo
531  enddo
532  end subroutine matvec
533 
534 end module m_dynamic_output
m_dynamic_output
This module provides functions to output result.
Definition: dynamic_output.f90:6
m_dynamic_output::fstr_dynamic_output
subroutine fstr_dynamic_output(cstep, istep, t_curr, hecMESH, fstrSOLID, fstrDYNAMIC, fstrPARAM, outflag)
Output result.
Definition: dynamic_output.f90:13
m_dynamic_output::matvec
subroutine matvec(y, x, hecMAT, ndof, D, AU, AL)
Definition: dynamic_output.f90:476
m_dynamic_output::fstr_dynamic_post
subroutine fstr_dynamic_post(fnum, cstep, istep, hecMESH, fstrSOLID, fstrDYNAMIC)
Summarizer of output data which prints out max and min output values.
Definition: dynamic_output.f90:98
m_dynamic_output::dynamic_output_monit
subroutine dynamic_output_monit(cstep, istep, t_curr, hecMESH, fstrPARAM, fstrDYNAMIC, fstrEIG, fstrSOLID)
Definition: dynamic_output.f90:357
m_fstr_nodalstress
This module provides functions to calculation nodal stress.
Definition: fstr_NodalStress.f90:6
m_fstr_nodalstress::fstr_nodalstress3d
subroutine fstr_nodalstress3d(hecMESH, fstrSOLID)
Calculate NODAL STRESS of solid elements.
Definition: fstr_NodalStress.f90:16
m_fstr_nodalstress::fstr_nodalstress6d
subroutine fstr_nodalstress6d(hecMESH, fstrSOLID)
Calculate NODAL STRESS of shell elements.
Definition: fstr_NodalStress.f90:1099
m_fstr_nodalstress::fstr_nodalstress2d
subroutine fstr_nodalstress2d(hecMESH, fstrSOLID)
Calculate NODAL STRESS of plane elements.
Definition: fstr_NodalStress.f90:853
m_fstr
This module defines common data and basic structures for analysis.
Definition: m_fstr.F90:15
m_fstr::iresult
integer(kind=kint), pointer iresult
Definition: m_fstr.F90:124
m_fstr::ilog
integer(kind=kint), parameter ilog
FILE HANDLER.
Definition: m_fstr.F90:109
m_fstr::ivisual
integer(kind=kint), pointer ivisual
Definition: m_fstr.F90:125
m_hecmw2fstr_mesh_conv
HECMW to FSTR Mesh Data Converter. Converting Connectivity of Element Type 232, 342 and 352.
Definition: hecmw2fstr_mesh_conv.f90:8
m_hecmw2fstr_mesh_conv::fstr2hecmw_mesh_conv
subroutine fstr2hecmw_mesh_conv(hecMESH)
Definition: hecmw2fstr_mesh_conv.f90:27
m_hecmw2fstr_mesh_conv::hecmw2fstr_mesh_conv
subroutine hecmw2fstr_mesh_conv(hecMESH)
Definition: hecmw2fstr_mesh_conv.f90:15
m_make_result
This module provide a function to prepare output of static analysis.
Definition: make_result.f90:6
m_make_result::fstr_make_result
subroutine, public fstr_make_result(hecMESH, fstrSOLID, fstrRESULT, istep, time, fstrDYNAMIC)
MAKE RESULT for static and dynamic analysis (WITHOUT ELEMENTAL RESULTS) -----------------------------...
Definition: make_result.f90:611
m_make_result::setup_contact_output_variables
subroutine, public setup_contact_output_variables(hecMESH, fstrSOLID, phase, dtime)
Definition: make_result.f90:1637
m_make_result::fstr_write_result
subroutine, public fstr_write_result(hecMESH, fstrSOLID, fstrPARAM, istep, time, flag, fstrDYNAMIC)
OUTPUT result file for static and dynamic analysis.
Definition: make_result.f90:22
m_fstr::fstr_dynamic
Data for DYNAMIC ANSLYSIS (fstrDYNAMIC)
Definition: m_fstr.F90:530
m_fstr::fstr_eigen
Package of data used by Lanczos eigenvalue solver.
Definition: m_fstr.F90:618
m_fstr::fstr_param
FSTR INNER CONTROL PARAMETERS (fstrPARAM)
Definition: m_fstr.F90:156
m_fstr::fstr_solid
Definition: m_fstr.F90:242