[esp-r] Re: surface attributes for adjacent zones

[Marco Manzan]

Alle domenica 14 giugno 2009, esp-r-request at lists.strath.ac.uk ha scritto:

>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Sun, 14 Jun 2009 09:51:23 +0530
> From: "NALIN SWAROOP" <nalin_swaroop at iitb.ac.in>
> Subject: [esp-r]  surface attributes for adjacent zones
> To: esp-r at lists.strath.ac.uk
> Message-ID: <a038aada663f19b608867f42f3c16053.squirrel at gpo.iitb.ac.in>
> Content-Type: text/plain;charset=utf-8
>
> Hello all
>
> I have a 8 floor 43 zones building but the the problem is every time I try
> to to give the boundary conditions for zone adjacency(surface in another
> zone option) for one zone and then move on to next zone, after some time I
> check the boundary conditions for the previously defines zones, they are
> changed randomly to some other surfaces in some other zones and I am not
> able to correctly define boundary conditions and the topology vertex
> connection test gives me errors. Is there a prescribed method for defining
> adjacent zones??
>
> TIA
>
> Nalin
>
>
>
> ------------------------------
>
> Message: 2
> Date: Sun, 14 Jun 2009 10:07:53 +0530
> From: "NALIN SWAROOP" <nalin_swaroop at iitb.ac.in>
> Subject: [esp-r]  surface attributes for adjacent zones keep changing
> To: esp-r at lists.strath.ac.uk
> Message-ID: <ff729665bdecf6160185fcf9cfdfaccc.squirrel at gpo.iitb.ac.in>
> Content-Type: text/plain;charset=utf-8
>
> Hello all
>
> I have a 8 floor 43 zones building but the the problem is every time I try
> to to give the boundary conditions for zone adjacency(surface in another
> zone option) for one zone and then move on to next zone, after some time I
> check the boundary conditions for the previously defines zones, they are
> changed randomly to some other surfaces in some other zones and I am not
> able to correctly define boundary conditions and the topology vertex
> connection test gives me errors. Is there a prescribed method for defining
> adjacent zones??
>
> I have a similar problem to this one
> http://lists.strath.ac.uk/archives/esp-r/2006/000104.html
>
> TIA
>
> Nalin
>
>
>
>
> ------------------------------
>
> _______________________________________________
> esp-r mailing list
> esp-r at lists.strath.ac.uk
> http://lists.strath.ac.uk/mailman/listinfo/esp-r
>
> End of esp-r Digest, Vol 35, Issue 7
> ************************************

I had some problems with topology checking too. I finally spotted a possible 
error in source code:
in file esruprj/edonecon.F I changed the following lines in routine LINDIR

       logical close
       PI = 4.0 * ATAN(1.0)
       DIST=SQRT((X2-X1)**2+(Y2-Y1)**2+(Z2-Z1)**2)
      THETA=ASIN((Z2-Z1)/DIST)
in
with 
 
C
C added a check on The value of ASIN argument
C for some geometries the argument becomes greater then one the
C result is THETA = NaN 
C
      ARG = (Z2-Z1)/DIST
      IF (abs(ARG) .LE. 1) THEN
         THETA=ASIN(ARG)
      ELSE
         THETA=SIGN(PI/2,ARG)
      ENDIF

I hope this helps  solving the topology problem.

regards

Marco

-- 
 _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _        
|Marco Manzan               Email: manzan at units.it            |        
|Tel:   (+39) 040 5583506   Fax  : (+39) 040  572033          |        
|Dipartimento di Ingegneria Navale, del Mare e per l'Ambiente |        
|Universita' degli Studi di Trieste                           |        
|piazzale Europa 1, 34127 TRIESTE                             |        
|_ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _ _|