COMSOL Multiphysics™

elgpspec
elgpspec can be added in the fem.elem field unless meshextend is called with
property 'standard' set to 'off'.
Note that the sum of weights in an explicit pattern specification are supposed to be
equal to the element’s volume in the element local coordinate system, that is, to 1/
2 for triangles and 1/6 for tetrahedra.
Examples
Given a fem structure with an xmesh field, extract elements and add an additional
pattern which uses fourth order integration on curved simplices and explicit zeroth
order integration on other simplices. Finally, the xmesh is updated.
elstr = fem.xmesh.getElements;
clear elem;
for i=1:length(elstr)
elem{i} = eval(elstr(i));
if strcmp(elem{i}.elem,'elgpspec')
igpspec = i;
end
end
elstr = fem.xmesh.getInitElements;
clear eleminit;
for i=1:length(elstr)
eleminit{i} = eval(elstr(i));
end
newptrn = length(elem{igpspec}.geom{1}.ep)+1;
elem{igpspec}.geom{1}.ep{newptrn} = {'s(2)', '4', ...
'ls(2)', {{'0.333333333','0.333333333','0.5'}}};
fem.elem = elem;
fem.eleminit = eleminit;
Here, additional elements using the integration point pattern with index newptrn
can be added to the fem.elem field.
fem.xmesh = meshextend(fem,'standard','off');
See Also
elempty, elepspec, elshape, eleqc, eleqw
81