1! 2! CalculiX - A 3-dimensional finite element program 3! Copyright (C) 1998-2021 Guido Dhondt 4! 5! This program is free software; you can redistribute it and/or 6! modify it under the terms of the GNU General Public License as 7! published by the Free Software Foundation(version 2); 8! 9! 10! This program is distributed in the hope that it will be useful, 11! but WITHOUT ANY WARRANTY; without even the implied warranty of 12! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 13! GNU General Public License for more details. 14! 15! You should have received a copy of the GNU General Public License 16! along with this program; if not, write to the Free Software 17! Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. 18! 19 subroutine sigini(sigma,coords,ntens,ncrds,noel,npt,layer, 20 & kspt,lrebar,rebarn) 21! 22! user subroutine sigini 23! 24! INPUT: 25! 26! coords coordinates of the integration point 27! ntens number of stresses to be defined 28! ncrds number of coordinates 29! noel element number 30! npt integration point number 31! layer currently not used 32! kspt currently not used 33! lrebar currently not used (value: 0) 34! rebarn currently not used 35! 36! OUTPUT: 37! 38! sigma(1..ntens) residual stress values in the integration 39! point. If ntens=6 the order of the 40! components is 11,22,33,12,13,23 41! 42 implicit none 43! 44 character*80 rebarn 45 integer ntens,ncrds,noel,npt,layer,kspt,lrebar 46 real*8 sigma(*),coords(*) 47! 48 sigma(1)=-100.d0*coords(2) 49 sigma(2)=-100.d0*coords(2) 50 sigma(3)=-100.d0*coords(2) 51 sigma(4)=0.d0 52 sigma(5)=0.d0 53 sigma(6)=0.d0 54! 55 return 56 end 57 58