49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
75
76
77
78#include "implicit_f.inc"
79
80
81
82#include "radioss_maptable.inc"
83#include "com04_c.inc"
84#include "param_c.inc"
85#include "scr17_c.inc"
86
87
88
89 INTEGER, INTENT(IN) :: IPART(LIPART1,*),SH4TREE(*),
90 . SH3TREE(*),IXS(NIXS,*),IXC(NIXC,*),IXTG(NIXTG,*),IXS10(6,*),
91 . IXQ(NIXQ,*),IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),IPARTS(*),
92 . IPARTQ(*),IPARTC(*),IPARTG(*),IPARTT(*),IPARTP(*),IPARTR(*),
93 . KNOD2ELS(*),KNOD2ELC(*),KNOD2ELTG(*),NOD2ELS(*),NOD2ELC(*),
94 . NOD2ELTG(*),NOD2ELQ(*),KNOD2ELQ(*)
95 INTEGER ELTYP,JCLAUSE,OPT_A,OPT_O,OPT_E
96
97 INTEGER, INTENT(IN), DIMENSION(NUMNOD,2) :: ITABM1
98 INTEGER, INTENT(IN), DIMENSION(NRBODY,2) :: IRBODYM
99 INTEGER, INTENT(IN), DIMENSION(2,NRBODY) :: RBY_MSN
100
101
102
103 LOGICAL :: IS_AVAILABLE
105 . x(3,*),skew(lskew,*)
106 CHARACTER(LEN=NCHARFIELD) :: KEYSET
107 CHARACTER(LEN=NCHARTITLE) :: SET_TITLE
108 TYPE (SET_SCRATCH) ::
109
110 TYPE (SET_) :: CLAUSE
111 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(*)
112 TYPE (BOX_) , DIMENSION(NBBOX) :: IBOX
113
114
115
116 INTEGER BOXTYPE,SURF_LIST,LINE_LIST_1D,ELTYP_ALL
117
118 INTEGER ADMBID
119 DATA admbid/0/
120
121
122
123
124 IF (keyset == 'box') BOXTYPE = 1
125 IF (KEYSET == 'box2') BOXTYPE = 2
126!
127!
128 ! ---------------------
129 ! Fill Boxes
130 ! ---------------------
131!
132!
133 ! NODE box
134 CALL CREATE_NODE_BOX(
135 . CLAUSE ,ITABM1 ,JCLAUSE ,IS_AVAILABLE ,LSUBMODEL,
136 . IBOX ,X ,SKEW ,SET_TITLE ,KEYSET )
137
138
139 ! Solid's in box
140 IF (numels > 0)
142 * clause ,iparts ,jclause ,is_available ,lsubmodel,
143 * ibox ,x ,skew ,set_title ,keyset ,
144 * boxtype ,numels ,nixs ,ixs ,8 ,
145 * ipart ,admbid ,admbid ,admbid ,elt_solid)
146
147
148 IF (numelq > 0)
150 * clause ,ipartq ,jclause ,is_available ,lsubmodel,
151 * ibox ,x ,skew ,set_title ,keyset ,
152 * boxtype ,numelq ,nixq ,ixq ,4 ,
153 * ipart ,admbid ,admbid ,admbid ,elt_quad )
154
155
156 IF (numelc > 0)
158 * clause ,ipartc ,jclause ,is_available ,lsubmodel,
159 * ibox ,x ,skew ,set_title ,keyset ,
160 * boxtype ,numelc ,nixc ,ixc ,4 ,
161 * ipart ,sh4tree ,3 ,ksh4tree ,elt_sh4n )
162
163
164 IF (numeltg > 0 .AND. numeltria == 0)
166 * clause ,ipartg ,jclause ,is_available ,lsubmodel,
167 * ibox ,x ,skew ,set_title ,keyset ,
168 * boxtype ,numeltg ,nixtg ,ixtg ,3 ,
169 * ipart ,sh3tree ,3 ,ksh3tree ,elt_sh3n )
170
171
172 IF (numeltria > 0)
174 * clause ,ipartg ,jclause ,is_available ,lsubmodel,
175 * ibox ,x ,skew ,set_title ,keyset ,
176 * boxtype ,numeltria,nixtg ,ixtg ,3 ,
177 * ipart ,admbid ,admbid ,admbid ,elt_tria )
178
179
180 IF (numelt > 0)
182 * clause ,ipartt ,jclause ,is_available ,lsubmodel,
183 * ibox ,x ,skew ,set_title ,keyset ,
184 * boxtype ,numelt ,nixt ,ixt ,2 ,
185 * ipart ,admbid ,admbid ,admbid ,elt_truss)
186
187
188 IF (numelp > 0)
190 * clause ,ipartp ,jclause ,is_available ,lsubmodel,
191 * ibox ,x ,skew ,set_title ,keyset ,
192 * boxtype ,numelp ,nixp ,ixp ,2 ,
193 * ipart ,admbid ,admbid ,admbid ,elt_beam )
194
195
196 IF (numelr > 0)
198 * clause ,ipartr ,jclause ,is_available ,lsubmodel,
199 * ibox ,x ,skew ,set_title ,keyset ,
200 * boxtype ,numelr ,nixr ,ixr ,2 ,
201 * ipart ,admbid ,admbid ,admbid ,elt_spring)
202
203
204
205 line_list_1d = clause%NB_TRUSS + clause%NB_BEAM + clause%NB_SPRING
206 IF (line_list_1d > 0)
207
208
210 * .false. )
211
212
213
214 surf_list = clause%NB_SOLID + clause%NB_QUAD +
215 + clause%NB_SH4N + clause%NB_SH3N
216 IF (surf_list > 0) THEN
217
218
220 * ixs ,ixs10 ,sh4tree ,sh3tree ,ixc ,
221 * ixtg ,knod2els ,nod2els ,knod2elc ,nod2elc ,
222 * knod2eltg ,nod2eltg ,ipartc ,ipartg ,iparts ,
223 * ipart ,clause ,opt_a ,opt_o ,ixq ,
224 * knod2elq ,nod2elq ,x ,keyset ,delbuf ,
225 * .false. )
226
227
229 . .false.)
230 ENDIF
231
232
233
234 IF (clause%NB_RBODY > 0) THEN
236 . ibox ,x ,skew ,set_title ,keyset
237 . rby_msn )
238 ENDIF
239
subroutine create_elt_box(clause, iparte, jclause, is_available, lsubmodel, ibox, x, skew, set_title, keyset, boxtype, numel, nix, ix, nix1, ipart, eltree, klevtree, keltree, eltyp)
subroutine create_line_from_element(ixt, ixp, ixr, clause, delbuf, go_in_array)
subroutine create_line_from_surface(clause, keyset, opt_a, opt_e, delbuf, go_in_array)
subroutine create_rbody_box(clause, irbodym, jclause, is_available, lsubmodel, ibox, x, skew, set_title, keyset, rby_msn)
subroutine create_surface_from_element(ixs, ixs10, sh4tree, sh3tree, ixc, ixtg, knod2els, nod2els, knod2elc, nod2elc, knod2eltg, nod2eltg, ipartc, ipartg, iparts, ipart, clause, opt_a, opt_o, ixq, knod2elq, nod2elq, x, keyset, delbuf, go_in_array)
integer, parameter nchartitle
integer, parameter ncharfield