UsesCase_MEDstructElement_1.c

Aller à la documentation de ce fichier.

/*  This file is part of MED.
 *
 *  COPYRIGHT (C) 1999 - 2020  EDF R&D, CEA/DEN
 *  MED is free software: you can redistribute it and/or modify
 *  it under the terms of the GNU Lesser General Public License as published by
 *  the Free Software Foundation, either version 3 of the License, or
 *  (at your option) any later version.
 *
 *  MED is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *  GNU Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public License
 *  along with MED.  If not, see <http://www.gnu.org/licenses/>.
 */
 
/*
 * StructElement use case 1 : write struct elements model in a file
 * STEP 1 : suppport mesh creation
 * STEP 2 : struct element model creation
 * STEP 3 : computation mesh creation
 */
 
#include <med.h>
#define MESGERR 1
#include <med_utils.h>
 
#include <string.h>
 
int main (int argc, char **argv) {
  med_idt fid;
  char elementname[MED_NAME_SIZE+1]="";
  const char ballsupportname[MED_NAME_SIZE+1]="BALL_SUPPORT_MESH";
  const char beamsupportname[MED_NAME_SIZE+1]="BEAM_SUPPORT_MESH";
  const char computmeshname[MED_NAME_SIZE+1]="COMPUT_MESH";
  const char beamsectionname[MED_NAME_SIZE+1]="BEAM_SECTION_MESH";
  const med_int elementdim = 3;
  med_int nnode;
  med_geometry_type geotype=MED_NONE;
  med_int ncomp;
  const med_float ballmeshnodescoo[3] = {0.0, 0.0, 0.0 };
  const med_float beammeshnodescoo[3*7] = {0.0,0.0,0.0,
                                           0.0,0.0,2.0,
                                           0.0,0.0,4.0,
                                           0.0,0.0,5.0,
                                           0.0,0.0,7.0,
                                           0.0,0.0,10.0,
                                           0.0,0.0,11.0 };
  const med_float beamsectioncoo[9*3] = {-0.2,-0.2,0.0,
                                         0.0,-0.2,0.0,
                                         0.2,-0.2,0.0,
                                         -0.2, 0.0,0.0,
                                         0.0, 0.0,0.0,
                                         0.2, 0.0,0.0,
                                         -0.2, 0.2,0.0,
                                         0.0, 0.2,0.0,
                                         0.2, 0.2,0.0 };
  const med_int seg2connectivity[2*6] = {1,2, 2,3, 3,4, 4,5, 5,6, 6,7};
  med_int spacedim, meshdim,nseg2;
  /*                                         123456789012345612345678901234561234567890123456 */
  const char axisname[3*MED_SNAME_SIZE+1] = "x               y               z               "