nray119 / matlab-stl-files Goto Github PK

surface_all_nodes.xlsx
File 1: surface_all_nodes.csv
In this file each line has the following format:
ID of the vertex belonging to the lung surface in the Finite Element Mesh, followed by X, Y, Z coordinate of the vertex (in Fortran77 format input).

Please note that the format of the files should be identical to these templates so that we can use the script to get the Boundary Conditions.

I had the same problem for extracting the surfaces in ABAQUS. Think of it this way:
1- a face belongs to the lung surface if and only if it belongs to one element.
2- a face belongs to the tumor surface if and only if it belongs to two elements with two distinct/ unequal materials.

These are the conclusions of several iterations to solve this problem. Please let me know if you have any questions with the underlying reasoning.
You have to write a script for COMSOL that exports these meshes. I think COMSOL is compatible with MATLAB but I am not sure.

There are two files that need to be written. The templates for them are attached to this email. File 1 in .csv format. File 2 is stl.
File 1: surface_all_nodes.csv
In this file each line has the following format:
ID of the vertex belonging to the lung surface in the Finite Element Mesh, followed by X, Y, Z coordinate of the vertex (in Fortran77 format input).

Please note that the format of the files should be identical to these templates so that we can use the script to get the Boundary Conditions.

This is the Fortran code that reads the binary file and assigns the material properties.
The code reads the DATA.bin file records- this is the header of the file.

  READ(FID,'(2I5)',REC=1) N_T, N_MAX
  WRITE(FMT,'("(", I4, "F10.4)")') N_T
  READ(FID,FMT,REC=2)(TIMES(J),J=1,N_T)

It then interpolates the alpha1 and alpha2 coefficients over time and calculates the boundary conditions according to the formula.
I recommend you read this code to understand what the fortran code does. The point of the fortran code is just to read and calculate boundary conditions for the model in ABAQUS. You will have to write something similar for COMSOL but as you see, the code is fairly simple and short.
DISP.for.txt