tri3-q.fem

来自「InsightToolkit-1.4.0(有大量的优化算法程序)」· FEM 代码 · 共 69 行

% NOTE: in this element, you have to list the "corner" nodes first and then
%       list the edge nodes.  Note that nodes 0,1,2 below form the triangle
%       and nodes 3,4,5 lie along the edges formed by 0,1,2.  If you list
%       the nodes continuously, Itpack will complain that you have negative
%       values along the diagonal of K.

<Node>
	0	% Global node number
	2 0 0	% Nodal coordinates
<Node>
	1
	2 2 0
<Node>
	2	% Global node number
	2 2 2	% Nodal coordinates
<Node>
	3
	2 1 0
<Node>
	4	% Global node number
	2 2 1	% Nodal coordinates
<Node>
	5
	2 1 1
<END>	% End of nodes

<MaterialLinearElasticity>
	0	% Global material number
	E  : 30000000	% Young modulus
	A  : 0	% Crossection area
	I  : 0	% Moment of inertia
	nu : 0.3	% Poisson's ratio
	h  : 1
	RhoC : 1
	END:	% End of material definition

<END>	% End of materials

<Element2DC0QuadraticTriangularStrain>
	0	% Global element number
	0	% Node 1 ID
	1	% Node 2 ID
	2	% Node 3 ID
	3	% Node 4 ID
	4 	% Node 5 ID
	5	% Node 6 ID
	0	% MaterialLinearElasticity ID
<END>	% End of elements

<LoadBC>
	0
	  0
	  0
	1 0
<LoadBC>
	1
	  0
	  1
	1 0
<LoadNode>
	2
	0
	2
	2 10 0

<END>	% End of loads