FEAP User Forum

FEAP => Input File Issues => Topic started by: yolandacruz12 on February 17, 2019, 07:01:48 AM

Title: Problem I calculate solid mesh
Post by: yolandacruz12 on February 17, 2019, 07:01:48 AM

Hello,

I am new here I explain my problem. I have this practice whose objective is to enter the generation of solid meshes. The solid to be studied is two-dimensional, and is the one defined in the following figure attachment
(https://drive.google.com/open?id=1bac0UA3-aNMQszYZopPfzeMN7-IuiJ0v)
Since the geometry and the loads are symmetric, we are going to model only a quarter of the problem, as shown in the following figure


To solve this problem we went to the input file Ip4 for
feap that contains the definition of the problem.
Analyze the contents of the file, which is the following:
feap * * MEF - Practice 4 - BLOCK - curved sides * *
0,0,0,2,2,4! Nod, elt, mat, dim, gdl, npe
! - end of line
PARAmeter! Parameters
a = 100! width (mm)
b = 100! long (mm)
12
G. Rus and A. Martínez
chap. 1. Management of FEAP; meshing
r = 30! radius (mm)
n = 15! division 1
m = 15! division 2
p = atan (1.0)! pi rooms
s = r * sin (p)! coordinates of a bisecting point
c = r * cos (0.5 * p)! coordinates for quadratic side
d = r * sin (0.5 * p)! coordinates for quadratic side
! - end of line
BLOCk 1! Block of nodes and elements
CARTesian, n, m, 1,1,1,0,0! Type, jump1, jump2, nod1, elt
1, mat1, skip, type
1, r, 0! Corner 1
2, a, 0! Corner 2
3, a, b! Corner 3
4, s, s! Corner 4
8, c, d! To do quadratic side
! - end of line
BLOCk 2! Block of nodes and elements
CARTesian, n, m, m * (n + 1) +1,0,1,0,0! Type, jump1, jump2,
nod1, elt1, mat1, skip, type
1, s, s! Corner 1
2, a, b! Corner 2
3.0, b! Corner 3
4.0, r! Corner 4
8, d, c! To do quadratic side
! - end of line
MATErial, 1! Material 1
SOLID! Type of problem
ELAStic, ISOTropic, 210.0.3! Material properties
! - end of line
END! End of the preprocess
! - end of line
BATCh! Starts post-processing
PLOT, MESH! Draw the mesh
END! Finalizes postprocess
! - end of line
INTEractive! Gives the user control
! - end of line
STOP! End of file

I put everything but the program does not return results, can you help me?

Title: Re: Problem I calculate solid mesh
Post by: Prof. S. Govindjee on February 17, 2019, 08:53:11 AM

Have you looked in you output file?  What has been read correctly?  What error is shown?

Title: Re: Problem I calculate solid mesh
Post by: yolandacruz12 on February 17, 2019, 09:20:07 AM

Hi Prof,

The truth is that I am a little novice with this program I come from matlab and I am with interfeace matfeap. The error I got is [1] PETSC ERROR: Null Object: Parameter # 1!

I'm looking at the input parameters and settings

Thanks

Title: Re: Problem I calculate solid mesh
Post by: Prof. S. Govindjee on February 17, 2019, 05:17:02 PM

I suggest trying to run FEAP by itself before adding the complexity of the matfeap interface.

Title: Re: Problem I calculate solid mesh
Post by: Prof. R.L. Taylor on February 18, 2019, 01:35:01 PM

You need to check your input file carefully.  You are missing many ',' commas -- you have '.' period instead. 
You will need loads and boundary conditions before you solve the problem.

Title: Re: Problem I calculate solid mesh
Post by: Prof. R.L. Taylor on February 18, 2019, 07:21:30 PM

One other error.
In expressions do not leave blanks.   Instead of  m * (m - 1) write m*(m-1)  the blanks confuse the  parser on different parameters sometimes.

Always  check the output file (the one with O....)  for any differences in what you expect an input to create.